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

:..---------------------------------作者:_____________-----------------------------日期::_____________数字图像处理作业(第三次)直方图处理数字图像处理作业(第三次)直方图处理第一题第二题():1)编写函数Y=MyHisteq(X),对灰度图象X进行直方图均衡化,得到均衡化后的图像Y,假设X的数据类型为uint8。2)利用自己编写的直方图均衡化函数,,。:(1)函数程序:functionY=MyHisteq(X)%对灰度图象X进行直方图均衡化,得到均衡化后的图像Y[m,n]=size(X);A=zeros(1,256);fork=0:255A(k+1)=length(find(X==k))/(m*n);endS1=zeros(1,256);fori=1:256forj=1:iS1(i)=A(j)+S1(i);%算SkendendS2=round(S1*256);fori=1:256B(i)=sum(A(find(S2==i)));endfigure,bar(0:255,B,'b');title('均衡化后的直方图');xlabel('灰度'),ylabel('出现概率')Y=X;fori=0:255Y(find(X==i))=S2(i+1);endend(2)主程序:X=imread('Fig0316(1)(top_left).tif');figure,imshow(X);Y=MyHisteq(X);figure,imshow(Y);:第三题():编写程序,,。:X=imread('Fig0327(a)(tungsten_original).tif');figure,imshow(X);Y=histeq(X);figure,imshow(Y);E=;k0=;k1=;k2=;[m,n]=size(X);X1=X;Mg=sum(X1(:))/(m*n);Vg=sum(sum((X1-Mg).^2))/(m*n);fori=2:m-1forj=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&&Vs<=k2*Vg)X1(i,j)=E*X1(i,j);endendendfigure,imshow(X1);:放大了约130倍的钨丝的SME图像:全局直方图均衡的结果:使用局部直方图统计增强后的图像::用自己编写编程实现中值滤波,。:(1)函数程序:functionY=mid_