数字图像处理MATLAB方法..docVIP

图像处理 一．计算图像统计参数：最大值：扫描法：设第一个象素灰度级为最大值，然后与后一个象素灰度级比较，如果后面的大，则把后一象素值赋给最大值，依次逐一比较到图象最后的一个象素，找出最大值。程序：[A,map]=imread(E:\lena.bmp);??????? %显示原图imshow(A,map),???????????????????????????????????? title(原图);image=double(A);%求负片k=image(1,1);for i=1:256,?? for j=1:256,????? if kimage(i,j);????? k=image(i,j);end,?? end,end,k%显示负片figure,imshow(rimage,map),title(负片); 最小值：最大值求法类似，把第一个象素灰度级设为最小值，然后逐一比较。程序：[A,map]=imread(E:\lena.bmp);??????? %显示原图imshow(A,map),???????????????????????????????????? title(原图);image=double(A);%求负片k=image(1,1);for i=1:256,?? for j=1:256,????? if kimage(i,j);????? k=image(i,j);end,?? end,end,k%显示负片figure,imshow(rimage,map),title(负片); 均值：把图中的所有的象素值累加起来再除以图象的象素的个数。程序：[A,map]=imread(E:\lena.bmp);??????? %显示原图imshow(A,map),???????????????????????????????????? title(原图);image=double(A);%求负片k=image(1,1);for i=1:256,?? for j=1:256,????? if kimage(i,j);????? k=image(i,j);end,?? end,end,k%显示负片figure,imshow(rimage,map),title(负片); 直方图：把图中的不同灰度等级的象素分别累加起来，再分别除以图象象素的总数。即：p(sk)=nk/n。可以得到原图的灰度值分布图。参考程序：[A,map]=imread(e:\lena.bmp);??????? imshow(A,map),???????????????????????????????????? title(原图);image=double(A);for i=1:256,??? a(i)=0;end,for i=1:256,?? for j=1:256,?????? b=image(i,j)+1;???? ?????? a(b)=a(b)+1;????????? ?? end,end,for k=1:256,??? p(k)=a(k)/(256^2);end,figure,stem(p);title(直方图); 二．对比度增强：直方图均衡：原始图的直方图变换为均匀分布的形式，这样就增加了象素灰度值的动态范围从而达到增强对比度的效果。程序：[A,map]=imread(e:\lena.bmp);??????? imshow(A,map),???????????????????????????????????? title(原图);image=double(A);for i=1:256,??? n(i)=0;end,for i=1:256,?? for j=1:256,?????? s=image(i,j)+1;???? ?????? n(s)=n(s)+1;????????? ?? end,end,for k=1:256,??? p(k)=n(k)/(256^2);end,figure,bar(p,r);title(直方图);for k=1:256,??? q(k)=0;end,for k=1:256,???? ??? for j=1:k,??????? q(k)=q(k)+p(j);??? end,end,figure,bar(q,y