图像增强与MATLAB实现

1、数字图像处理数字图像处理课程设计课程设计课设题目课设题目：图像增强与 MATLAB 实现学校学院学校学院：华东交通大学理学院学生班级学生班级：13 级信息计算（2）班学生姓名学生姓名：魏超学生学号学生学号：20130810010216指导老师：指导老师：范自柱图像增强与 MATLAB 实现摘摘 要要数字图像处理是指将图像信号转换成数字格式并利用计算机对其进行处理的过程。图像增强是数字图像处理的过程中经常采用的一种方法，它对提高图像质量起着重要的作用。本文先对图像增强的原理进行概述，然后对图像增强的方法分类并给出直方图增强、对比度增强、平滑和锐化等几种常用的增强方法的理论基础，通过 Matlab

2、 实验得出的实际处理效果来对比各种算法的优缺点，讨论不同的增强算法的技术要点，并对其图像增强方法进行性能评价。关键字：关键字：图像；图像增强；算法目 录一、一、MATLABMATLAB 的简介的简介 .1 11.1 MATLAB 主要功能.1二、二、MATLABMATLAB 的主要功能的主要功能 .1 12.1 数字增强技术概述 .12.2 数字图像的表示 .1三、直方图的均衡化三、直方图的均衡化.2 23.1 图像的灰度 .23.2 灰度直方图 .23.3 直方图均衡化 .3四、图像二值化四、图像二值化.5 54.1 图像二值化 .5五、对比度增强五、对比度增强.6 65.1 对比度增强 .

3、65.2 灰度调整 .65.3 对数变换 .7六、滤波六、滤波.8 86.1 平滑滤波 .86.2 线性平滑滤波程序： .86.3 非线性滤波 .9七、锐化七、锐化.1313八、参考文献八、参考文献.1414九、自我评价九、自我评价.1414一、一、MatlabMatlab 的简介的简介1.1 MATLAB 主要功能MATLAB 是建立在向量、数组和矩阵基础上的一种分析和仿真工具软件包，包含各种能够进行常规运算的“工具箱”，如常用的矩阵代数运算、数组运算、方程求根、优化计算及函数求导积分符号运算等；同时还提供了编程计算的编程特性，通过编程可以解决一些复杂的工程问题；也可绘制二维、三维图形，输出

4、结果可视化。目前，已成为工程领域中较常用的软件工具包之一。二、二、MATLAB 的主要功能的主要功能2.1 数字增强技术概述图像增强是按特定的需要突出一幅图像中的某些信息，同时，消弱或去除某些信息使得图像更加实用。图像增强技术主要包含直方图修改处理、图像平滑处理、图像尖锐化处理等。图像增强技术主要包括：直方图修改处理，图像平滑处理，图像尖锐化处理，彩色图像处理。从纯技术上讲主要有两类：频域处理法和空域处理法。频域处理法主要是卷积定理，采用修改图像傅立叶变换的方法实现对图像的增强处理技术；空域处理法：是直接对图像中的像素进行处理，基本上是以灰度映射变换为基础的。2.2 数字图像的表示 图像并不能

5、直接用计算机来处理，处理前必须先转化成数字图像。由于从外界得到的图像多是二维（2-D）的，一幅图像可以用一个 2-D 数组表示。),(yxf这里 x 和 y 表示二维空间 XY 中一个坐标点的位置，而 f 则代表图像在点的某种性质数值。为了能够用计算机对图像进行处理，需要坐标空间和性),(yx质空间都离散化。三、直方图的均衡化三、直方图的均衡化3.1 图像的灰度常用的图像一般是灰度图，这时 f 表示灰度值，反映了图像上对应点的亮度。亮度是观察者对所看到的物体表面反射光强的量度。作为图像灰度的量度函数应大于零。人们日常看到的图像一般是从目标上反射出来的光组成),(yxf的，所以可看成由两部分构成

6、：入射到可见场景上光的量；场景中目标),(yxf对反射光反射的比率。确切地说它们分别称为照度成分和反射成分),(yxi。与和都成正比，可表示成),(yxr),(yxf),(yxi),(yxr。),(yxf),(yxi),(yxr3.2 灰度直方图灰度变换是图像增强的一种重要手段，使图像对比度扩展，图像更加清晰，特征更加明显。灰度级的直方图给出了一幅图像概貌的描述，通过修改灰度直方图来得到图像增强。程序如下：%灰度直方图I=imread(D:imagerice.bmp);%读取图像subplot(2,1,1);Imshow(I);%显示图像title(a)原图);subplot(2,1,2);i

7、mhist(I);%绘制图像的灰度直方图title(b)原图的灰度直方图);图 4.13.3 直方图均衡化直方图均衡化过程如下：（1）计算原图像的灰度直方图；)(KrrP（2）计算原图像的灰度累积分布函数，进一步求出灰度变换表；ks（3）根据灰度变换表，将原图像各灰度级映射为新的灰度级。灰度直方图均衡化程序：%灰度直方图均衡化I=imread(D:imagebubbles.bmp);%读取图像subplot(2,2,1);Imshow(I);%显示图像title(原图);subplot(2,2,2);imhist(I);%绘制图像的灰度直方图title(原图的灰度直方图);subplot(2,

8、2,3);J=histeq(I,64);%对图像进行均衡化处理，返回有 64 级灰度的图像 JImshow(J);%显示图像title(原图直方图均衡化);subplot(2,2,4);imhist(J);%绘制图像的灰度直方图title(均衡后的灰度直方图)以下展示了直方图均衡化的效果：图4.2优势：能够使得处理后图像的概率密度函数近似服从均匀分布，其结果扩张了像素值的动态范围，是一种常用的图像增强算法。不足：不能抑制噪声。四、图像二值化四、图像二值化4.1 图像二值化图像的二值化处理就是将图像上的像素点的灰度值设置为 0 或 255，也就是讲整个图像呈现出明显的黑白效果。首先，图像的二值化

9、有利于图像的进一步处理，使图像变得简单，而且数据量减小，能凸显出感兴趣的目标的轮廓。其次，要进行二值图像的处理与分析，首先要把灰度图像二值化，得到二值化图像。图像二值化程序：%图像二值化 （选取一个域值，(5)将图像变为黑白图像）I=imread(beauty.tif);bw=im2bw(I,0.5);%选取阈值为 0.5subplot(1,3,1);imshow(I);title(原图);subplot(1,3,2);imshow(bw);title(显示二值图像);J=find(I=150);I(J)=255;subplot(1,3,3);imshow(I);title( 图像二值化 (

10、域值为 150 ） );下图为图像二值化的效果：图 4.3五、对比度增强五、对比度增强5.1 对比度增强对比度增强是按一定的规则修改输入图像每一个像素的灰度，从而改变图像灰度的动态范围。例如，观察图 4.4 可以发现，该图的对比度不高其灰度直方图没有低于 35 或高于 210 的值，如果将图像数据映射到整个灰度范围内，则图像的对比度将大大增大。图 4.45.2 灰度调整灰度调整程序：%imadjust 函数I=imread(D:imagerice.bmp);%读取图像subplot(2,2,1);Imshow(I);%显示图像title(原图);subplot(2,2,2);imhist(I)

11、;%绘制图像的灰度直方图title(原图的灰度直方图);subplot(2,2,3);J=imadjust(I,0.3 0.7,);%对图像进行灰度变换Imshow(J);%显示图像title(原图直方图均衡化);subplot(2,2,4);imhist(J);%绘制图像的灰度直方图title(均衡后的灰度直方图)以下展示了常用对比度扩展法的结果： 图 4.5从图 4.5【原图】可以看出原始图像动态范围较小，整体较暗，反映在直方图上像素主要集中在低灰度的一侧。经过对比度调整，图像变亮。5.3 对数变换对数变化常用来扩展低值灰度，压缩高值灰度，这样可以使低值灰度的图像细节更容易看清。对数变换程

12、序：对数变换程序：%对数变换I=imread(beauty.tif);%读取图像I=mat2gray(I);%对数变换不支持 uint8 类型数据，将一个矩阵转化为灰度图像的数据格式(double)J=log(I+1);subplot(1,2,1);Imshow(I);%显示图像title(原图);subplot(1,2,2);Imshow(J);title(对数变换后的图像)图 4.6对数变换确实能够扩展低值灰度，而压缩高值灰度，使低值灰度的图像细节更容易看清。六、滤波六、滤波6.1 平滑滤波平滑技术用于平滑图像中的噪声。平滑噪声可以在空间域中进行，基本方法是求像素灰度的平均或中值。为了既平

13、滑噪声又保护图像信号。6.2 线性滤波输出图像的值等于输入图像滤波后值的局部平均，各个项具有相同的权。下面是平滑窗口分别为矩形和圆形的情况。 1111111111111111111111111251,kjhrect0111011111111111111101110211,kjhcirc对一些图像进行线性滤波可以去除图像中某些类型的噪声，如采用邻域平均法的均值滤波器就非常适用于去除通过扫描得到的图像中的颗粒噪声。6.2 线性平滑滤波程序：%线性平滑滤波I=imread(beauty.tif);J=imnoise(I,salt & pepper,0.02);subplot(221),ims

14、how(I)title(原图像)subplot(222),imshow(J)title(添加椒盐噪声图像)K1=filter2(fspecial(average,3),J)/255;%应用 3*3 邻域窗口法subplot(223),imshow(K1)title(3x3 窗的邻域平均滤波图像)K2=filter2(fspecial(average,7),J)/255;%应用 7*7 邻域窗口法subplot(224),imshow(K2)title(7x7 窗的邻域平均滤波图像)如图 4.8(a)为线性平滑滤波的例子：6.3 非线性滤波中值滤波是一种最常用的图像增强技术，是非线性滤波。对椒盐

15、噪声有很好的去噪效果。下图是加高斯噪声后，中值滤波和平均滤波的滤波效果程序：%中值滤波和平均滤波I=imread(girl.bmp);I=rgb2gray(I);J=imnoise(I,gaussian,0,0.01);subplot(2,2,1);imshow(I);title(原图);subplot(2,2,2);imshow(J);title(noise);K=fspecial(average,5);K1=filter2(K,J)/255;subplot(2,2,3);imshow(K1);title(平均滤波);L=medfilt2(J,3 5);subplot(2,2,4);imsh

16、ow(L);title(中值滤波);对比：图 4.9(a)但对于椒盐噪声，中值滤波效果比均值滤波效果好。中值滤波程序：%二维中值滤波hood=3;P,map=imread(leno.bmp);I=rgb2gray(P);imshow(I,map);noisy=imnoise(I,salt & pepper,0.05); %加入强度为 0.05 的椒盐噪声subplot(221);imshow(noisy,map);title(加入椒盐噪声);filtered1=medfilt2(noisy,hood hood);%33 窗口二维中值滤波subplot(222);imshow(filte

17、red1,map);title(33 窗口);hood=5;filtered2=medfilt2(noisy,hood hood); %55 窗口二维中值滤波subplot(223);imshow(filtered2,map);title(55 窗口);hood=7;filtered3=medfilt2(noisy,hood hood); %77 窗口二维中值滤波subplot(224);imshow(filtered3,map); title(77 窗口);I=imread(leno.bmp);figure;imshow(I);title(原图);实现中值滤波的效果图如下：图 4.9(b)比

18、较发现，图像的噪声点被去除；图像的边缘稍微的变得模糊。高通滤波边缘增强程序：%高通滤波边缘增强I=imread(girl.bmp);I=rgb2gray(I);subplot(2,2,1);imshow(I);title(original pic);J=fspecial(average,3);J1=conv2(I,J)/255;subplot(2,2,2);imshow(J1);title(3*3lowpass);K=fspecial(prewitt);K1=filter2(K,J1)*5;subplot(2,2,3);imshow(K1);title(prewitt);L=fspecial(

19、sobel);L1=filter2(L,J1)*5;subplot(2,2,4);imshow(L1);title(sibel);下图为高通滤波边缘增强的例子：图 4.9(c)优势：去噪效果明显，并且能够较好的保持图像边缘位置和细节。不足：非线性滤波算法的实现相对线性滤波比较困难。七、锐化七、锐化图像锐化处理的作用是使灰度反差增强，从而使模糊图像变得更加清晰。图像模糊的实质就是图像受到平均运算或积分运算，因此可以对图像进行逆运算，如微分运算以突出图像细节使图像变得更为清晰。锐化程序：%锐化a=imread(dowels.tif);subplot(131);imshow(a);title(原图)

20、;b=double(a);%将图像矩阵转化为 double 类型s=size(b);c=zeros(s(1,1),s(1,2);for x=2:s(1,1)-1 for y=2:s(1,2)-1 c(x,y)=(-b(x+1,y)-b(x-1,y)-b(x,y+1)-b(x,y-1)+4*b(x,y); endend%用拉氏算子对图像进行滤波，这个过程相当于运用了一个 33 的掩膜0，1，0；1，4，1；0，1，0subplot(132);imshow(c);title(Laplace 锐化滤波图像);d=b+c;%当拉普拉斯掩膜中心系数为正时，增强图像为原图像于拉氏算子滤波图像之和d=uint8(d);%将图像矩阵变回 uint8 格式subplot(133);imshow(d);title(Laplace 锐化滤波结果);图 4.10 展示了在 MAT