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

数字图像处理作业（第三次）直方图处理
第一题第二题(例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);
end
end
(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-1
for j=2:n-1
X2=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);
end
end
end
figure,imshow(X1);
放大了约130倍的钨丝的SME图像：全局直方图均衡的结果：使用局部直方图统计增强后的图像：
1. 问题：
用自己编写编程实现中值滤波，用96页例3.14的图测试。

2．程序：
(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;Y=mid_filter(X,n);
figure,imshow(Y);
被椒盐噪声污染的电路板X射线图像：用3×3中值滤波器降噪后的图像：。