空间域数字图像的平滑与锐化.doc

实验三空间域数字图像的平滑与锐化
实验目的和要求
掌握空间域数字图像的平滑与锐化。
实验内容和原理
利用加权平均掩模实现数字图像的平滑；
利用拉普拉斯算子实现数字图像的锐化
实验环境
硬件：一般PC机
操作系统：WindowsXP
编程平台：MATLAB或高级语言
算法描述及实验步骤
Code：
X=imread('moon.tif');
subplot(2,2,1);imshow(X);title原图
b=size(X);
X=double(X);
%f=[0-10;-14-1;0-10;];%用四领域
f=[-1-1-1;-18-1;-1-1-1;];%用八领域
g=[121;242;121;];%模糊用的算子
Y=zeros(b);
for(i=2:b(1)-1)
for(j=2:b(2)-1)Y(i,j)=X(i,j)*g(2,2)+X(i+1,j)*g(3,2)+X(i,j+1)*g(2,3)+X(i+1,j+1)*g(3,3)+X(i+1,j-1)*g(3,1)+X(i-1,j+1)*g(1,3)+X(i-1,j-1)*g(1,1)+X(i-1,j)*g(1,2)+X(i,j-1)*g(2,1);
end;
end;
Y=mat2gray(Y/16);
subplot(2,2,2);imshow(Y);title模糊后
Z=zeros(b);
for(i=2:b(1)-1)
for(j=2:b(2)-1)Z(i,j)=Y(i,j)*f(2,2)+Y(i+1,j)*f(3,2)+Y(i,j+1)*f(2,3)+Y(i+1,j+1)*f(3,3)+Y(i+1,j-1)*f(3,1)+Y(i-1,j+1)*f(1,3)+Y(i-1,j-1)*f(1,1)+Y(i-1,j)*f(1,2)+Y(i,j-1)*f(2,1);
end;
end;
Z=mat2gray(Z);
subplot(2,2,3);imshow(Z);title锐化后
M=zeros(b);
for(i=2:b(1)-1)
for(j=2:b(2)-1)
M(i,j)=X(i,j)+Y(i,j);
end;
end;
M=mat2gray(M);
subplot(2,2,4);imshow(M);title增强后

调试过程

实验结果

总结