资源描述:
《数字图像处理实验报告图像增强实验》由会员上传分享,免费在线阅读,更多相关内容在工程资料-天天文库。
1、课程名称专业班级:字像处理导论名号电气与信息学陵和谐勤奋求是创新实验题目图像增强实验-空域增强实验室DSP室&信号室实验时间实验类别设计同组人数2成绩指导教师签字:一.实验目的1.掌握图像增强的基本定义及0的。2.理解空间域增强的基本原理及方法。3.掌握进行图像的空域增强的方法。二.实验内容1.2.读出blurry_moon.tif这幅图像,采用3X3的拉普拉斯算子w=[1,1,1;1-81;1,1,1]对其进行滤波。编写函数w=genlaplacian(n),自动产生任一奇数尺寸n的拉普拉斯算子,如5X5的拉普拉斯算
2、子[111111111111-24111111111111]3.分别采用5X5,9X9,15X15和25X25大小的拉普拉斯算子对blurry_moon.tif进行锐化滤波,并利用式=完成图像的锐化增强,观察其有何不同,要求在同一窗口中显示。4.采用不同的梯度算子对blurry_moon.tif进行锐化滤波,并比较其效果。5.自己设计锐化空间滤波器,并将其对噪声图像进行处理,示处理后的图像;三.实验具体实现1.读出blurry_moon.tif这幅图像,采用3X3的拉普拉斯算子w=[1,1,1;1-81;1,1,1]对
3、其进行滤波。I=imread(’trees,tif’);T=double(T);subplot(1,2,1),imshow(T,[]);title(’OriginalImage’);w=[1,1,1;1,-8,1;1,1,1];K=conv2(T,w,’same’);subplot(1,2,2)imshow(K);titleCLaplacianTransformation’);2.编写函数wgenlaplacian(n)w=[1111自动产生任一奇数尺寸n的拉普拉斯算子111111111-24111111如5X5的拉普
4、拉斯算子11111]I=imread(’cameraman,tif’);T=double(I);subplot(121),imshow(T,[]);title(,OriginalTmage’);1;1,w=[1,1,1,1,1;1,1,1,1,1;1,1,-24,1,1;1,1,1,1,1,1,1,1;];K=conv2(T,w,’same');subplot(122),imshow(K);titleClaplaciantransform’);3.分别采用5X5,9X9,15X15和25X25大小的拉普拉斯算子对blu
5、rry_moon.tif进行锐化滤波,并利完成隴的锐化增强,观察贿何不同,要求在同-窗口屮显示。functionw=genlaplacian(n)w=ones(n);x=ceil(n/2);w(x,x)=-1*(n*n-1);end;wl=genlaplacian(5);w2=genlaplacian(9);w3=genlaplacian(15);w4二genlaplacian(25);I=imread(,blurry^moon.tif’);T=double(I);subplot(231),imshow(T,[]);t
6、itleCOriginalTmage’);Kl=conv2(T,wl,’same’);subplot(232),imshow(Kl);titleClaplacianoperator5氺5’);K2=conv2(T,w2,’same’);subplot(233),imshow(K2);titleClaplacianoperator9氺9’);K3=conv2(T,w3,’same’);subplot(234),imshow(K3);title(’laplacianoperator15本15’);K4=conv2(T,w4
7、,’same’);subplot(235),imshow(Kd);title(,laplacianoperator25*25’);OriginalImagelaplacianoperator5W5laplacianoperator9*9laplacianoperator15*15laplacianoperator25*251.采用不同的梯度算子对bliiiTy_moon.tif进行锐化滤波,并比较其效果。[l,map]=imread(’trees.tif);l=double(l);subplot(2,3,1)imsho
8、w(l’map);title(’OriginalImage’);[Gx,Gy]=gradient(l);%gradientcalculationG=sqrt(Gx.*Gx+Gy.*Gy);%matrixJ1=G;%gradientlsubplot(2,3,2)imshow(J1,map);title(*OperatorlImage