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

  • 格式:doc
  • 大小:920.50 KB
  • 文档页数:17

下载文档原格式

  / 17
  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

实验报告

课程名称数字图像处理导论

专业班级_______________

姓名_______________

学号_______________

电气与信息学院

和谐勤奋求是创新

附录:可能用到的函数和参考结果**************报告里不能用参考结果中的图像

1)读出eight.tif这幅图像,给这幅图像分别加入椒盐噪声和高斯噪声后并与前一张

图显示在同一图像窗口中。

I=imread('cameraman.tif');

subplot(1,3,1)

imshow(I);

title(' Qriginal Image ');

J = imnoise(I,'salt & pepper',0.05); %noise density=0.05

subplot(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(' Qriginal Image ');

Sobel = imfilter(I,H,'replicate');

subplot(2,2,2)

imshow(Sobel);

title(' Sobel Image ')

H = fspecial('laplacian',0.4);

lap = imfilter(I,H,'replicate');

subplot(2,2,3)

imshow(lap);

title(' Laplacian Image ')

H = fspecial('gaussian',[3 3],0.5);

gaussian = imfilter(I,H,'replicate');

subplot(2,2,4)

imshow(gaussian);

title(' Gaussian Image ')

图2.2 原图像及各类低通滤波处理图像

3)使用函数imfilter时,分别采用不同的填充方法(或边界选项,如零填

充、’replicate’、’symmetric’、’circular’)进行低通滤波,显示处理后的图像。originalRGB = imread('sedemo_onion.png');

subplot(3,2,1)

imshow(originalRGB);

title(' Original Image ');

h = fspecial('motion', 50, 45); %motion blurred

filteredRGB = imfilter(originalRGB, h);

subplot(3,2,2)

imshow(filteredRGB);

title(' Motion Blurred Image ');

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 & pepper Noise');

h=fspecial('average'); %Averaging Filtering

J1=imfilter(J,h);

for i=1:10