基于matlab的直方图均衡化代码

《基于matlab的直方图均衡化代码》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、基于matlab的直方图均衡化代码2007-04-1520:15clearall%一，图像的预处理，读入彩色图像将其灰度化PS=imread('1.jpg');                %读入JPG彩色图像文件imshow(PS)                                 %显示出来title('输入的彩色JPG图像')imwrite(rgb2gray(PS),'PicSampleGray.bmp');%将彩色图片灰度化并保存PS=rgb2gray(PS);                           %灰度化后的数据存入数组%二，绘制直方图[m,n

2、]=size(PS);                            %测量图像尺寸参数GP=zeros(1,256);                           %预创建存放灰度出现概率的向量fork=0:255    GP(k+1)=length(find(PS==k))/(m*n);     %计算每级灰度出现的概率，将其存入GP中相应位置endfigure,bar(0:255,GP,'g')                   %绘制直方图title('原图像直方图')xlabel('灰度值')ylabel('出现概率')%三，直方图均衡化S1=zeros(1

3、,256);fori=1:256    forj=1:i         S1(i)=GP(j)+S1(i);                %计算Sk    endendS2=round((S1*256)+0.5);                         %将Sk归到相近级的灰度fori=1:256    GPeq(i)=sum(GP(find(S2==i)));          %计算现有每个灰度级出现的概率endfigure,bar(0:255,GPeq,'b')                 %显示均衡化后的直方图title('均衡化后的直方图')xlabel('

4、灰度值')ylabel('出现概率')%四，图像均衡化PA=PS;fori=0:255    PA(find(PS==i))=S2(i+1);               %将各个像素归一化后的灰度值赋给这个像素endfigure,imshow(PA)                          %显示均衡化后的图像title('均衡化后图像')imwrite(PA,'PicEqual.bmp');