实验五 图像的空间域增强

一、实验目的

1. 1. 熟悉图像空间域增强方法，掌握增强模板使用方法
   2. 掌握均值滤波器、中值滤波器的理论基础和实现方法
   3. 掌握图像锐化的基本理论和实现方法
   4. 验证图像滤波处理结果

二、实验原理与方法

图像增强是数字图像处理的基本内容之一，其目的是根据应用需要突出图像中的某些“有用”信息，削弱或去除不需要的信息，以改善图像的视觉效果,或突出图像的特征，便于计算机处理。图像增强可以在空间域进行， 也可以在频率域中进行。

空间域滤波主要利用空间模板进行，如 3´3，5´5 模板等。

 

一般来说，使用大小为 m×n 的滤波器对大小为 M×N 的图像 f 进行空 间滤波，可表示成：

其中，m=2a+1, n=2b+1, �(�,�)是滤波器系数，�(�,�)是图像值

均值滤波器是一种空间平滑滤波器，它是对包含噪声的图像上的每个像素点， 用它邻域内像素的平均值替代原来的像素值。

例如，采用一个 3×3 的模板，待处理的像素为 f(i,j)，则处理后图像对应的像素值为 g(i-1,j+1)

g(i,j)=1/9*(f(i-1,j-1)+f(i-1,j)+f(i-1,j+1)+f(i,j-1)+f(i,j)+ f(i,j+1)+f(i+1,j-1)+f(i+1,j)+f(i+1,j+1));

中值滤波器也是一种空间平滑滤波器，它是对以图像像素点为中心的一个滑动窗口内的诸像素灰度值排序，用中值代替窗口中心像素的原来灰度值，因此它是一种非   线性的图像平滑法。

采用 Laplacian 锐化算子进行图像边缘的锐化，是采用二阶差分运算获得像素间的差异值，由此，获得对图像景物边界的锐化。Laplacian 也可以算子也可以写成是模板作用的方式，如下：

设待处理的像素为 f(i,j)，则处理后图像对应的像素值为 g(i-1,j+1)，则g(i,j)=4*f(i,j) -(f(i-1,j-1)+f(i,j-1)+f(i,j+1)+f(i+1,j));

常用的锐化算子还有 Roberts、Prewitt 和 Sobel 算子等

三、实验内容与步骤

1. 读入一幅 256 级灰度的数字图像
2. 利用邻域平均法，分别采用3´3，5´5, 7´7, 9´9模板对加噪声图像进行平滑处理，显示原图像、加噪图像和处理后的图像。
3. 利用中值滤波法，分别采用3´3，5´5, 7´7, 9´9模板对加噪声图像进行去噪处理，显示原图像、加噪图像和处理后的图像。
4. 利用Laplacian 锐化算子（α=-1）对256 级灰度的数字图像进行锐化处理，显示处理前、后图像。

5. 分别利用Roberts、Prewitt 和Sobel 边缘检测算子，对数字图像进行边缘检测，显示处理前、后图像。

i = imread('20140423043802258.jpg');
i = rgb2gray(i);
imshow(i),title('原始图像');

%添加噪声
g1 = imnoise(i,'salt & pepper',0.02);%椒盐噪声占比图像的'2%'
g2 = imnoise(i,'gaussian',0,0.01);%高斯,均值为 0，方差为 0.01
figure(2);
subplot(121),imshow(g1),title('椒盐噪声污染图像');
subplot(122),imshow(g2),title('高斯噪声污染图像');

%用fspecial构建滤波器,用imfilter使用滤波器
h3 = fspecial('average',[3,3]); %average指定均值滤波器
h5 = fspecial('average',[5,5]);
h7 = fspecial('average',[7,7]);
fil3 = imfilter(g1,h3);
fil5 = imfilter(g1,h5);
fil7 = imfilter(g1,h7);
figure(3);
subplot(131),imshow(fil3),title('average3*3'),
subplot(132),imshow(fil5),title('average5*5'),
subplot(133),imshow(fil7),title('average7*7');

%使用中值滤波器，因为fspecial不能定义中值滤波器，中值滤波器有专门的函数medfilt2
fi3m = medfilt2(g1,[3,3]);
fi5m = medfilt2(g1,[5,5]);
fi7m = medfilt2(g1,[7,7]);
figure(4);
subplot(131),imshow(fi3m),title('median3*3')
subplot(132),imshow(fi5m),title('median5*5')
subplot(133),imshow(fi7m),title('median7*7');

%拉普拉斯算子
lap = fspecial('laplacian',0);
lapg1 = imfilter(i,lap);
lapg2 = i - (-1)*lapg1;
figure(5),imshow(lapg2),title('laplacian');

%边缘检测算子，用edge函数
b1 = edge(i,'sobel');
b2 = edge(i,'prewitt');
b3 = edge(i,'roberts');
b4 = edge(i,'log');%log就是拉普拉斯算子
figure(6);
subplot(141),imshow(b1),title('sobel');
subplot(142),imshow(b2),title('prewitt');
subplot(143),imshow(b3),title('roberts');
subplot(144),imshow(b3),title('laplacian');

 

不想写了哈哈哈，这几篇写得有点敷衍，记录生活😊

希望看到的人能获得想要的爱情😘🤣😎