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

1、精选优质文档-倾情为你奉上数字图像处理作业（第三次）直方图处理第一题第二题(例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

2、 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('Fi

3、g0316(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

4、(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&&Vs<=k2*Vg)X1(i,j)=E*X1(i,j);endendendfigure,imshow(X1);3. 实验结果：放大了约130倍的钨丝的SME图像：全局直方图均衡的结果：使用局部直方图统计增强后的图像：1. 问题：用自己编写编程实现中值滤波，用96页例3.14的图测试。2 程序

5、：(1) 函数程序：function Y=mid_filter(X,n)%自编的中值滤波函数。X是需要滤波的图像,n是模板大小(即n×n) height, width=size(X); %输入图像是p×q的,且p>n,q>n 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列元素,即模板(n×n的) 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;