基本小函数(2)

图像处理的有用资料

GP(k+1)=length(find(PS==k))/(m*n); %计算每级灰度出现的概率，将其存入GP中相应位置 end

figure,bar(0:255,GP,'g') %绘制直方图 figure NO 3 title('原图像直方图') xlabel('灰度值') ylabel('出现概率') %三，直方图均衡化 S1=zeros(1,256); for i=1:256 for j=1:i

S1(i)=GP(j)+S1(i); end end

S2=round(S1*256); for i=1:256

GPeq(i)=sum(GP(find(S2==i))); 率 end

figure,bar(0:255,GPeq,'b') 图 figure NO 4 title('均衡化后的直方图') xlabel('灰度值') ylabel('出现概率')

figure,plot(0:255,S2,'r') 线 figure NO 5 legend('灰度变化曲线') xlabel('原图像灰度级') ylabel('均衡化后灰度级') %四，图像均衡化 PA=PS; for i=0:255

PA(find(PS==i))=S2(i+1); 这个像素 end

figure,imshow(PA) %计算Sk %将Sk归到相近级的灰度 %计算现有每个灰度级出现的概%显示均衡化后的直方%显示灰度变化曲%将各个像素归一化后的灰度值赋给%显示均衡化后的图