untitled5.m

来自「图象处理中的直方图程序,有些文件名需自己改动一下」· M 代码 · 共 19 行

g=imread ('meinv.jpg');
I=rgb2gray(g);
J=I;
L=256; %灰度级 
nk=zeros(L,1);%出现次数
[row,col]=size(I);
n=row*col; %总像素个数
for i = 1:row
for j = 1:col
num = double(I(i,j))+1; %获取像素点灰度级
nk(num) = nk(num)+1; %统计nk 
end
end
figure;
subplot(3,1,1);imshow(J),title('原图');
subplot(3,1,2),plot(nk),title('直方图（nk）');