内容发布更新时间 : 2024/12/29 1:04:01星期一 下面是文章的全部内容请认真阅读。

《五邑大学学生实验报告》

五 邑 大 学

《数字图像处理与分析》

实验报告

实验时间： 实验班级：

实验报告总份（片）数： 份（片） 实验教师： 王 鹏

数学学院(系) 实验室

数学学院学院（系）信息与计算科学专业 班____组、学号 3116004982 姓名_______闫勇志_______协作者______________

实验项目 空间域图像增强（四）

教师评语：

一、 实验目的与要求 实验目的：

1、熟悉中值滤波器、Gauss滤波器、Laplace算子、Sobel算子等常

见空间域图像算法原理；

2、熟悉一般滤波器的工作原理，能够自己设计实现某一特定功能的

滤波器。

3、用上述滤波器的方法实现图像的增强。

实验要求：

1、对图像库中的Fig3.35(a).jpg，分别用matlab库函数实现中值滤波器、Gauss

滤波器、Laplace滤波器、Sobel滤波器对噪声图像进行滤波。比较滤波的结果，说

明各种滤波器的优缺点及适用场合；

2、根据中值滤波算法的思想，对图像库中的Fig3.35(a).jpg，自己编写一个实

现中值滤波的程序，同屏显示自己编写程序的结果和matlab库函数的结果。

3、对图像库中的Fig3.40(a).jpg，用空间滤波器的方法进行增强处理；

二、 实验方案 1.

I = double(imread('E:\\ images_chapter_03\\Fig3.35(a).jpg'));

[m,n] = size(I); X1 = zeros(m,n); X2 = zeros(m,n); X3 = zeros(m,n); X4 = zeros(m,n);

H2 = fspecial('gaussian'); %高斯 H3 = fspecial('laplacian');%拉普拉斯 H4 = fspecial('sobel'); %Sobel算子

X1 = medfilt2(I); %中值滤波 X2 = imfilter(I,H2); %高斯 X3 = imfilter(I,H3); %拉普拉斯 X4 = imfilter(I,H4); %Sobel算子

subplot(2,3,1),imshow(uint8(I));title('原图像'); subplot(2,3,2),imshow(uint8(X1));title('中值滤波'); subplot(2,3,3),imshow(uint8(X2));title('高斯滤波'); subplot(2,3,4),imshow(uint8(X3));title('拉普拉斯滤波'); subplot(2,3,6),imshow(uint8(X4));title('Sobel滤波');

2.

%%自写滤波

I = double(imread('E:\\ images_chapter_03\\Fig3.35(a).jpg')); [m,n] = size(I); X = zeros(m,n); for i = 3:m-3

for j = 3:n-3 a =

[I(i-1,j-1),I(i-1,j),I(i-1,j+1),I(i,j-1),I(i,j),I(i,j+1),I(i+1,j-1),I(i+1,j),I(i+1,j+1)];

X(i,j) = median(reshape(a,numel(a),1)); end

end

subplot(1,2,1),imshow(uint8(I));title('原图像');

subplot(1,2,2),imshow(uint8(X));title('3X3中值滤波器');

3.

I = double(imread('E:\\ images_chapter_03\\Fig3.35(a).jpg')); [m,n] = size(I); X = zeros(m,n);

A = fspecial('sobel'); X = imfilter(I,A); A1 = X+I;

subplot(1,2,1),imshow(uint8(I));title('原图像');

subplot(1,2,2),imshow(uint8(A1));title('sobel滤波器增强图像');