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

数字图像处理作业(第三次)直方图处理
计143 吴彦博 143063
2
计143 吴彦博 143063
3
计143 吴彦博 143063
4
Y=MyHisteq(X
数字图像处理作业(第三次)直方图处理
计143 吴彦博 143063
2
计143 吴彦博 143063
3
计143 吴彦博 143063
4
Y=MyHisteq(X);
figure,imshow(Y);
3. 实验结果：
计143 吴彦博 143063
5
第三题()
1. 问题：
编写程序，，。
2．程序：
X=imread('Fig0327(a)(tungsten_original).tif');
figure,imshow(X);
计143 吴彦博 143063
6
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);
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);
计143 吴彦博 143063
7
3. 实验结果：
放大了约130倍的钨丝的SME图像：
全局直方图均衡的结果：
使用局部直方图统计增强后的图像：
1. 问题：
计143 吴彦博 143063
8
用自己编写编程实现中值滤波，。
程序：
(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矩阵变为一个行矩阵
计143 吴彦博 143063
9
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);
3. 实验结果：
被椒盐噪声污染的电路板X射线图像：
用3×3中值滤波器降噪后的图像：
计143 吴彦博 143063
10
计143 吴彦博 143063
11