数字图像处理实验报告——图像复原实验

.'.实验报告课程名称数字图像处理导论专业班级_______________姓名_______________学号_______________电气与信息学院和谐勤奋求是创新.'.实验题目图像复原实验-空域滤波复原实验室DSP室&信号室实验时间2015年10月13日实验类别设计同组人数2成绩指导教师签字：一．实验目的1.掌握图像滤波的基本定义及目的。2.理解空间域滤波的基本原理及方法。3.掌握进行图像的空域滤波的方法。二．实验内容1.读出eight.tif这幅图像，给这幅图像分别加入椒盐噪声和高斯噪声后并与前一张图显示在同一图像窗口中。2.对加入噪声图像选用不同的平滑（低通）模板做运算，对比不同模板所形成的效果，要求在同一窗口中显示。3.使用函数imfilter时，分别采用不同的填充方法（或边界选项，如零填充、’replicate’、’symmetric’、’circular’）进行低通滤波，显示处理后的图像。4.运用for循环，将加有椒盐噪声的图像进行10次，20次均值滤波，查看其特点,显示均值处理后的图像（提示:利用fspecial函数的’average’类型生成均值滤波器）。5.对加入椒盐噪声的图像分别采用均值滤波法，和中值滤波法对有噪声的图像做处理，要求在同一窗口中显示结果。6.自己设计平滑空间滤波器，并将其对噪声图像进行处理，显示处理后的图像。三．实验具体实现1.读出（自己选定.tif）这幅图像，给这幅图像分别加入椒盐噪声和高斯噪声后并与前一张图显示在同一图像窗口中。I=imread('trees.tif');subplot(1,3,1)imshow(I);title('OriginalImage');J=imnoise(I,'salt&pepper',0.05);%noisedensity=0.05subplot(1,3,2)imshow(J);title('salt&pepper');K=imnoise(I,'gaussian',0.01,0.01);subplot(1,3,3)imshow(K);title('gaussian').'.2.对加入噪声图像选用不同的平滑（低通）模板做运算，对比不同模板所形成的效果，要求在同一窗口中显示。I=imread('moon.tif');H=fspecial('sobel');subplot(2,2,1)imshow(I);title('QriginalImage');Sobel=imfilter(I,H,'replicate');subplot(2,2,2)imshow(Sobel);title('SobelImage')H=fspecial('laplacian',0.4);lap=imfilter(I,H,'replicate');subplot(2,2,3)imshow(lap);title('LaplacianImage').'.H=fspecial('gaussian',[33],0.5);gaussian=imfilter(I,H,'replicate');subplot(2,2,4)imshow(gaussian);title('GaussianImage')3.使用函数imfilter时，分别采用不同的填充方法（或边界选项，如零填充、’replicate’、’symmetric’、’circular’）进行低通滤波，显示处理后的图像。originalRGB=imread('trees.tif');subplot(3,2,1)imshow(originalRGB);title('QriginalImage');h=fspecial('motion',50,45);%motionblurredfilteredRGB=imfilter(originalRGB,h);subplot(3,2,2)imshow(filteredRGB);title('MotionBlurredImage');boundaryReplicateRGB=imfilter(originalRGB,h,'replicate');subplot(3,2,3).'.imshow(boundaryReplicateRGB);title('0-Padding');boundary0RGB=imfilter(originalRGB,h,0);subplot(3,2,4)imshow(boundary0RGB);title('Replicate');boundarysymmetricRGB=imfilter(originalRGB,h,'symmetric');subplot(3,2,5)imshow(boundarysymmetricRGB);title('Symmetric');boundarycircularRGB=imfilter(originalRGB,h,'circular');subplot(3,2,6)imshow(boundarycircularRGB);title('Circular');.'.4.运用for循环，将加有椒盐噪声的图像进行10次，20次均值滤波，查看其特点,显示均值处理后的图像（提示:利用fspecial函数的’average’类型生成均值滤波器）。I=imread('kids.tif');J=imnoise(I,'salt&pepper',0.05);subplot(1,3,1)imshow(J);title('salt&pepperNoise');h=fspecial('average');%AveragingFilteringJ1=imfilter(J,h);fori=1:10J1=imfilter(J,h);subplot(1,3,2)imshow(J1);title('10AveragingFiltering');endJ2=imfilter(J,h);fori=1:20J2=imfilter(J,h);subplot(1,3,3)imshow(J2);title('20AveragingFiltering');end.'.5.对加入椒盐噪声的图像分别采用均值滤波法，和中值滤波法对有噪声的图像做处理，要求在同一窗口中显示结果。I=imread('trees.tif');J=imnoise(I,'salt&pepper',0.05);subplot(1,3,1)imshow(J);title('OriginalImage');h=fspecial('average');%AveragingFilteringJ1=imfilter(J,h);subplot(1,3,2)imshow(J1);.'.title('AveragingFiltering');J2=medfilt2(J);%MedianFilteringsubplot(1,3,3)imshow(J2);title('MedianFiltering');.'.6.自己设计平滑空间滤波器，并将其对噪声图像进行处理，显示处理后的图像。domain=[00800;00800;88888;00800;00800];I=imread('trees.tif');J=imnoise(I,'salt&pepper',0.05);subplot(1,2,1)imshow(J);title('OriginalImage');K1=ordfilt2(J,5,domain);subplot(1,2,2)imshow(K1);title('5*5SmoothingFiteredImage');.'..'.附录：可能用到的函数和参考结果**************报告里不能用参考结果中的图像1)读出eight.tif这幅图像，给这幅图像分别加入椒盐噪声和高斯噪声后并与前一张图显示在同一图像窗口中。I=imread('cameraman.tif');subplot(1,3,1)imshow(I);.'.title('QriginalImage');J=imnoise(I,'salt&pepper',0.05);%noisedensity=0.05subplot(1,3,2)imshow(J);title('salt&pepper');K=imnoise(I,'gaussian',0.01,0.01);subplot(1,3,3)imshow(K);title('gaussian');图2.1初始图像及椒盐噪声图像、高斯噪声污染图2)对加入噪声图像选用不同的平滑（低通）模板做运算，对比不同模板所形成的效果，要求在同一窗口中显示。I=imread('trees.tif');H=fspecial('sobel');subplot(2,2,1)imshow(I);title('QriginalImage');Sobel=imfilter(I,H,'replicate');subplot(2,2,2)imshow(Sobel);title('SobelImage')H=fspecial('laplacian',0.4);lap=imfilter(I,H,'replicate');subplot(2,2,3)imshow(lap);title('LaplacianImage')H=fspecial('gaussian',[33],0.5);gaussian=imfilter(I,H,'replicate');subplot(2,2,4)imshow(gaussian);title('GaussianImage').'.图2.2原图像及各类低通滤波处理图像3)使用函数imfilter时，分别采用不同的填充方法（或边界选项，如零填充、’replicate’、’symmetric’、’circular’）进行低通滤波，显示处理后的图像。originalRGB=imread('sedemo_onion.png');subplot(3,2,1)imshow(originalRGB);title('OriginalImage');h=fspecial('motion',50,45);%motionblurredfilteredRGB=imfilter(originalRGB,h);subplot(3,2,2)imshow(filteredRGB);title('MotionBlurredImage');boundaryReplicateRGB=imfilter(originalRGB,h,'replicate');subplot(3,2,3)imshow(boundaryReplicateRGB);title('0-Padding');boundary0RGB=imfilter(originalRGB,h,0);subplot(3,2,4)imshow(boundary0RGB);title('Replicate');boundarysymmetricRGB=imfilter(originalRGB,h,'symmetric');subplot(3,2,5)imshow(boundarysymmetricRGB);title('Symmetric');boundarycircularRGB=imfilter(originalRGB,h,'circular');subplot(3,2,6)imshow(boundarycircularRGB);title('Circular');.'.图2.3原图像及运动模糊图像图2.4函数imfilter各填充方式处理图像4)运用for循环，将加有椒盐噪声的图像进行10次，20次均值滤波，查看其特点,显示均值处理后的图像。I=imread('kids.tif');J=imnoise(I,'salt&pepper',0.05);subplot(1,3,1)imshow(J);title('salt&pepperNoise');h=fspecial('average');%AveragingFilteringJ