数字图像处理作业(第三次)直方图处理.docx

计143 吴彦博 143063数字图像处理作业（第三次）直方图处理第一题第二题(例3.6)1. 问题：1）编写函数Y=MyHisteq(X), 对灰度图象X进行直方图均衡化，得到均衡化后的图像Y，假设X的数据类型为uint8。 2）利用自己编写的直方图均衡化函数，重现例3.6中的实验，即对图3.16中的图像进行直方图均衡化。2程序：(1) 函数程序：function Y = MyHisteq(X)%对灰度图象X进行直方图均衡化，得到均衡化后的图像Y m,n=size(X); A=zeros(1,256); for k=0:255 A(k+1)=length(find(X=k)/(m*n); end S1=zeros(1,256); for i=1:256 for j=1:i S1(i)=A(j)+S1(i); %算Sk end end S2=round(S1*256); for i=1:256 B(i)=sum(A(find(S2=i); end figure,bar(0:255,B,b); title(均衡化后的直方图); xlabel(灰度),ylabel(出现概率) Y=X; for i=0:255 Y(find(X=i)=S2(i+1); endend(2) 主程序：X=imread(Fig0316(1)(top_left).tif);figure,imshow(X);Y=MyHisteq(X);figure,imshow(Y);3. 实验结果：第三题(例3.12)1. 问题：编写程序，重现例3.12中的实验，产生图3.27。2程序：X=imread(Fig0327(a)(tungsten_original).tif);figure,imshow(X);Y=histeq(X);figure,imshow(Y);E=4.0;k0=0.4;k1=0.02;k2=0.4;m,n=size(X);X1=X;Mg=sum(X1(:)/(m*n);Vg=sum(sum(X1-Mg).2)/(m*n);for i=2:m-1for j=2:n-1X2=X1(i-1:i+1,j-1:j+1);Ms=sum(X2(:)/9;Vs=sum(sum(X2-Ms).2)/9;if(Ms=k0*Mg&k1*Vg=Vs&Vsn,qn X1=double(X); X2=X1; for i=1:height-n+1 for j=1:height-n+1 C=X1(i:i+(n-1),j:j+(n-1); %取出X1中从(i,j)开始的n行n列元素,即模板(nn的) E=C(1,:); %是C矩阵的第一行 for k=2:n E=E,C(k,:); %将C矩阵变为一个行矩阵 end mm=median(E); %mm是中值 X2(i+(n-1)/2,j+(n-1)/2)=mm;%将模板各元素的中值赋给模板中心位置的元素 end end Y=uint8(X2);end(2) 主程序：X=imread(Fig0335(a)(ckt_board_saltpep_prob_pt05).tif);figure,imshow(X);n=3;Y=mid_