基于 matlab 的直方均衡化代码 2007-04-15 20:15 clear all %一，像的预处理，读入彩像将其灰度化 PS=imread(1.jpg); %读入 JPG 彩像文件 imshow(PS) %显示出来 title(输入的彩 JPG 像) imwrite(rgb2gray(PS),PicSampleGray.bmp); %将彩片灰度化并保存 PS=rgb2gray(PS); %灰度化后的数据存入数组%二，绘制直方 m,n=size(PS); %测量像尺寸参数 GP=zeros(1,256); %预创建存放灰度出现概率的 向量 for k=0:255GP(k+1)=length(find(PS=k)/(m*n); %计算每级灰度出现的概率， 将其存入 GP 中相应位置 end figure,bar(0:255,GP,g) %绘制直方 title(原像直方) xlabel(灰度值) ylabel(出现概率)%三，直方均衡化 S1=zeros(1,256); for i=1:256for j=1:iS1(i)=GP(j)+S1(i); %计算 Skend end S2=round(S1*256)+0.5); %将 Sk 归到相近级的灰 度 for i=1:256GPeq(i)=sum(GP(find(S2=i); %计算现有每个灰度级出现 的概率 end figure,bar(0:255,GPeq,b) %显示均衡化后的直方 title(均衡化后的直方) xlabel(灰度值) ylabel(出现概率)%四，像均衡化 PA=PS;for i=0:255PA(find(PS=i)=S2(i+1); %将各个像素归一化后的灰 度值赋给这个像素 end figure,imshow(PA) %显示均衡化后的像 title(均衡化后像) imwrite(PA,PicEqual.bmp);