基本小函数(8)

图像处理的有用资料

clear all

%一，图像的预处理，读入彩色图像将其灰度化

PS=imread('PicSample.jpg'); %读入JPG彩色图像文件 imshow(PS) %显示出来 figure NO 1 title('输入的彩色JPG图像')

imwrite(rgb2gray(PS),'PicSampleGray.bmp'); %将彩色图片灰度化并保存 PS=rgb2gray(PS); %灰度化后的数据存入数组 figure,imshow(PS) %显示灰度化后的图像，也是均衡化前的样品 figure NO 2 title('灰度化后的图像')

%二，绘制直方图

[m,n]=size(PS); GP=zeros(1,256); 向量

for k=0:255

GP(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: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('灰度变化曲线')

%测量图像尺寸参数

%预创建存放灰度出现概率的%计算每级灰度出现的概率，将%绘制直方图 figure NO 3 %计算Sk %将Sk归到相近级的灰度 %计算现有每个灰度级出现的%显示均衡化后的直方%显示灰度变化曲