图像处理的MATLAB实现实验一 空域图像增强

图像处理的MATLAB实现实验一 空域图像增强

一、实验目的

(1)掌握基本的空域图像增强方法,观察图像增强的效果,加深理解;

(2)了解空域平滑模板的特性及其对不同噪声的影响;

(3)了解空域锐化模板的特性及其对边缘的影响。

二、实验内容

(1)直方图处理:直方图均衡

(2)空域平滑:均值滤波、中值滤波;

实验要求

(1)用matlab语言进行仿真实验;

(2)递交实验报告,要求给出实验原理、源程序、实验结果及分析。

具体实验内容及要求

4.1 实验内容

4.1.1 直方图均衡

(1)读入原图像pollen.png并显示原图像以及直方图

(2)对原图像进行直方图均衡处理

(3)显示均衡后图像以及直方图。

4..1.2 图像空域平滑

(1)读入原图像lena.bmp并显示;

(2)对原图像分别添加高斯噪声和椒盐噪声,并显示加噪图像;

(3)采用均值滤波进行去噪处理,并显示去噪图像;

(4)采用中值滤波进行去噪处理,并显示去噪图像。

4.1.3 空域锐化

(1)读入原图像bridge.gif并显示;

(2)采用sobel算子对图像进行处理,并显示结果;

(3)尝试采用其他锐化模板进行处理。

4.2 实验原理

4.2.1 直方图均衡实验原理

对图像像素个数多的灰度级进行展宽,而对图像中像素个数少的灰度级进行压缩。而且,输入灰度级r与输出灰度级s的概率密度函数和有如下关系

积分形式如下 4.2.2 图像滤波

、椒盐噪声的中值滤波

由于椒盐噪声的出现使该点的像素比周围的亮或暗许多,如果在某个模板中,对像素由小到大重新排列,那么最暗或最亮的点一定被排在两侧,取模板中间位置的灰度值像素代替待处理图像像素的灰度值,从而达到滤除噪声的目的。

、高斯噪声的均值滤波

均值滤波是一种空域线性的滤波方法,用像素邻域内各像素的灰度平均值代替该像素原来的灰度值;均值滤波采用的是模板操作,将模板在图像中从左到右,从上到下的顺序移动将模板中心与每个像素重合;将模板中个系数与其对应的像素一并相乘,然后再经所有的结果一并相加;将上面相加的结果重新付给模板中心对应的像素点,那么该灰度值,就是经均值滤波后平滑后的灰度值。

4.2.3 图像锐化

对于图像锐化空域中有两种方式,第一个是使用sobel算子,即梯度锐化;而另一种是采用laplace算子,对于的八邻域内,在水平方向上和竖直方向上的差分方程为:

那么对于该像素点的梯度为;而对于采用laplace算子的差分方程为。采用这两种方式就可以进行图像的锐化。

4.3 实验源程序,结果及分析

4.3.1 直方图均衡

(1) 实验源程序及图像:

I=imread('C:\Documents and Settings\Administrator\桌面\实验\图像\pollen.png'); //由于是彩色图形,所以再求直方图时要先求RGB

imshow(I,[])
X=1:366
Y=1:366
I(X,Y,1:3)
>>R=ans(:,:,1);
>> G=ans(:,:,2);
>> B=ans(:,:,3);
imhist(R)
Imhist(G)
Imhist(B)
J=histeq(R,64)
imshow(J,[])
J=histeq(G,64)
Imshow(J)
J=histeq(B,64)
Imshow(J)
Imhist(J)

(2) 实验结果分析

变换以后,灰度级减少了,原直方图中频率较小的灰度级被合并到一个或几个灰度级中,即频率小的灰度级被压缩,频率大的灰度级被增强,这也是均衡化的实质,即减少图像的灰度级换取对比度的加大。本次实验选取的是一张彩色的图像,而彩色图像有R、G、B三个分量,因此先得将一个三维图像变成一幅二维的灰度图像,通过比较均衡以前和均衡以后的灰度图像,明显均衡以后的图像虽然对比度大大增加,但是图像的主要信息没有丢失,且更加符合人眼的视觉特性。

4.3.2 图像的滤波

(1)实验源程序

I=imread('C:\Documents and Settings\Administrator\桌面\实验\图像\lena.bmp');
subplot(331)
imshow(I);
title('原图像')
G1=imnoise(I,'gaussian',0,0.02);
G2=imnoise(I,'salt&pepper',0.08);
subplot(332)
imshow(G1);
title('被高斯噪声污染的图像')
subplot(333)
imshow(G2);
title('被椒盐噪声污染的图像'