新编数字图像处理实验

w =

1 1 1 1 1 1 1 1 1 1 1 1 -24 1 1 1 1 1 1 1 1 1 1 1 1

3) 分别采用5×5，9×9，15×15和25×25大小的拉普拉斯算子对图像进行滤波，并利用公式

2g) ?f ( x , y ) 完成对图像的锐化增强，观察其有何不同，要求在同一窗口显示 ( x , y ) ? f ( x , y ?

拉普拉斯算子滤波及图像增强 n=5;

w=ones(n); x=ceil(n/2);

w(x,x)=-1*(n*n-1);

I=imread('C:\\Program Files\\MATLAB\\R2009a\\work\\lena','bmp'); T=double(I);

K=conv2(T,w,'same'); J=T-K;

subplot(2,3,1);imshow(I); subplot(2,3,2);imshow(K); n=9;

w=ones(n); x=ceil(n/2);

w(x,x)=-1*(n*n-1);

I=imread('C:\\Program Files\\MATLAB\\R2009a\\work\\lena','bmp'); T=double(I);

K=conv2(T,w,'same'); J=T-K;

subplot(2,3,3);imshow(K); >> n=15; w=ones(n); x=ceil(n/2);

w(x,x)=-1*(n*n-1);

I=imread('C:\\Program Files\\MATLAB\\R2009a\\work\\lena','bmp'); T=double(I);

K=conv2(T,w,'same'); J=T-K;

subplot(2,3,4);imshow(K); >> n=25; w=ones(n); x=ceil(n/2);

w(x,x)=-1*(n*n-1);

I=imread('C:\\Program Files\\MATLAB\\R2009a\\work\\lena','bmp'); T=double(I);

K=conv2(T,w,'same');

J=T-K;

subplot(2,3,5);imshow(K);

初始图像与不同拉普拉斯算子锐化图像

close all n=5;

w=ones(n); x=ceil(n/2);

w(x,x)=-1*(n*n-1);

I=imread('C:\\Program Files\\MATLAB\\R2009a\\work\\lena','bmp'); T=double(I);

K=conv2(T,w,'same'); J=T-K;

subplot(1,3,1);imshow(I); subplot(1,3,2);imshow(K); subplot(1,3,3);imshow(J)

图像锐化的实质是将原图像与梯度信息叠加，相当于对目标物的边缘进行了增强。

拉普拉斯算子锐化与锐化增强图像

4) 采用不同的梯度算子对图像进行锐化滤波，并比较其效果 close all

>> [I,map]=imread('C:\\Program Files\\MATLAB\\R2009a\\work\\lena','bmp'); >> I=double(I);

subplot(2,3,1);imshow(I,map); >> [Gx,Gy]=gradient(I); >> G=sqrt(Gx.*Gx+Gy.*Gy); >> J1=G;

>> subplot(2,3,2);imshow(J1,map); >> J2=I;

>> K=find(G>=7); >> J2(K)=G(K);

>> subplot(2,3,3);imshow(J2,map); >> J3=I;

K=find(G>=7); J3(K)=255;

subplot(2,3,4);imshow(J3,map); >> J4=I;

K=find(G<=7); J4(K)=255;

subplot(2,3,5);imshow(J4,map); J5=I;

K=find(G<=7); J5(K)=0;

Q=find(G>=7);

>> subplot(2,3,6); imshow(J5,map);

5)自己设计锐化滤波器，并将其对噪声图像进行处理，显示处理后的图像 close all

I=imread('C:\\Program Files\\MATLAB\\R2009a\\work\\lena','bmp'); J=imnoise(I,'salt & pepper',0.5);

>> domain=[8 8 0 8 8;8 8 0 8 8;0 0 0 0 0;8 8 0 8 8;8 8 0 8 8]; >> K1=ordfilt2(J,5,domain); subplot(1,3,1);imshow(I) subplot(1,3,2);imshow(J) subplot(1,3,3);imshow(K1)

五、思考题

1. 均值滤波和中值滤波有那些优缺点？主要区别？对所给图像，哪种方法效果更佳？ 2. 锐化滤波算子各有什么特点？