图像对比度增强实验

电子科技大学通信学院 《图像对比度增强实验指导书》

2011年6月 图像对比度增强实验 一、 实验目的 本实验是针对拟开设的《图像处理》课程而开发的综合性实验。该实验与概率统计、高等数学、信号处理等课程都有一定程度的联系。通过本实验，促使学生利用上述课程相关知识，更好地掌握图像增强的基本方法和原理。要求学生达到：掌握经典的数字图像空域增强算法的实现方法；了解实验所用的图像处理算法的运行效果；学会分析与比较各种方法的差异。 二、 实验内容与要求

1. 2学时：掌握图像对比度增强原理； 2. 2学时：掌握图像变换原理； 3. 2学时：给定输入图片，在MATLAB上实现线性变换增强和直方图增强仿真；

4. 8学时：构建图像采集系统，在已有的采集软件平台上将仿真算法用C语言实现；

5. 2学时：测试实际输入图片增强效果。 三、 实验开设方式 本实验开设方式为两人一组实验，实验授课时间2学时，上机14学时。 四、 实验器材 硬件环境：

计算机一台； USB摄像头一个； 三脚架一个；

软件环境： Winxp操作系统； VC6.0开发环境；Open CV1.0，gsl-1.8 五、 实验原理 图像增强的目的是针对应用或人们主观需求，对输入图像进行某种处理，使得处理后的图像在特定结构或对比度等方面有明显的改善，其源头可以追溯到1969年。随着应用需求的不断发展，世界各国对此问题展开了广泛的研究。 特定结构或对比度不够清晰，可能来自质量不佳的成像设备、恶劣的大气条件、不恰当图像压缩算法等，其特点是图像中像素灰度集中在一个相对较小的范围，导致后端处理设备的观察者不能有效对图像进行判读。 造成图像对比度不够清晰的原因是多方面，包括电子系统的热噪声、光照过强、光照过弱、目标反射率过低、大雾天气、逆光拍摄、压缩等等。由于其产生原因的多样性，导致图像对比度不够清晰在图像中体现出来的现象也不同，可以分为以下3类：整体偏暗、整体偏亮、分布在亮和暗的两端，并且有可能在空间分布上存在多个区域。 现有的对比度增强技术根据其处理方法理论依据不同可分为：直方图均衡化、基于Retinex理论的图像增强、基于梯度场重建的图像增强；根据其处理范围，又可以分为全局处理与局部处理两大类。本试验将主要论述直方图均衡化和基于梯度场重建的图像增强两种方法，其中直方图均衡化仅描述基本原理，其实现由学生独立完成，而基于梯度场重建的图像增强方法，本实验将重点阐述，并给出参考代码，要求学生在此基础上进行进一步的完善。 直方图均衡化的基本原理就是对图像进行灰度变换。灰度变换有逆反处理、阈值变换、灰度拉伸、灰度切分、灰度级修正、动态范围调整等方法。虽然它们对图像的处理效果不同，但处理过程中都运用了点运算，通常可分为线性变换、分段线性变换、非线性变换。其缺点是需要用户根据不同的图像调整不同的变换函数。灰度变换是最简单的对比度增强技术，它可增大图像动态范围，扩展对比度，使图像清晰、特征明显，是图像增强的重要手段之一。它主要利用点运算来修正像素灰度，由输入像素点的灰度值确定相应输出点的灰度值，是一种基于图像变换的操作。灰度变换不改变图像内的空间关系，灰度级的改变是根据某种特定的灰度变换函数进行。 实验所需基本结构如图1所示。 图1 实验所需基本结构图 本实验把成像设备（即摄像头）采集的一幅图像，传入计算机，由图像增强技术算法实现的编程软件处理后，对图像进行增强，然后实时显示增强后的图像。 图像对比度增强技术，不仅要保持图像整体的一致性，还需要对图像的局部区域进行增强处理，使其具有最佳的表现力。因此，本实验采用梯度场方法。技术方案特点有： a) 梯度域增强避免了亮度不同对增强算子的影响 b) 重建图像是基于最小二乘法,与原始的图像在亮度方面不同 c) 重建图像在梯度域与原始图像具有强烈的相似性 d) 重建图像具有亮度平均值的相对保持性 e) 重建图像的边界条件周期延拓(采用DST变换要求) f) 所有算子都是直接计算 g) 较为复杂的DST变换有快速算法

（一） 算法原理 计算机处理的图像是离散的网格数据，必须对前面的连续函数描述形式进行离散化处理。 令原始图像为(,)Ixy，增强后图像为(,)Ixy，尺寸为（行X列）MN，1,xN，1,yM，则增强后的梯度图像为：

),(),(),(yxyxIyxG (1)

),(yxI是图像(,)Ixy的梯度，yx，定义图像的一阶前向差分： ),()1,(),(),1(),(yxIyxIyxIyxIyxI 令(,)xyIxyII (,)xy是梯度图像提升因子，









CyxIyx2),(

1),(

(2) 其中是梯度场增益因子，取值为0.3~0.45，2),(maxyxIC，是梯度场归一化因子，是梯度场归一化偏移因子，保证(,)xy为有界值，取值为，10C。

增强后的图像应满足：),(2),(minarg),(~yxIdxyGyxIyxI，根据变分原理可得，重建的图像与增强后的梯度图像满足如下泊松方程： GdivyxI),(

~

2 (3)

其中 ),(~4)1,(~)1,(~),1(~),1(~),(~2yxIyxIyxIyxIyxIyxI (4)

增强梯度场G梯度的散度用后向一阶差分逼近： )1,(),(),1(),(yxGyxGyxGyxGGdiv

yyxx (5)

求解泊松方程的边界条件为：将输入的原始图像(,)Ixy扩展为（行X列）22MN，0,1xN，0,1yM，且满足：











NxyIMxINxyIxIMyyNIyNIMyyIyI1),1()1,(1),1()0,(1),(),1(1),1(),0(

(6) 重建图像(,)Ixy的取值可能超出8bit范围，需要对其进行归一化处理 ),(~maxmI~yxIax，),(~minmI~yxIin (7)

归一化后的图像为： inaxinyxIImI~mI~mI~),(~*255

(8)

下面将阐述如何从泊松方程中求取重建图像

令2T，),(yxIU，FGdiv)(，向量U是重构图像),(~yxI的一个逼近，矩阵T是一个‘1 1 -4 1 1’的对角矩阵，F是包含边界条件的散度矩阵。

2cos22cos2





NxMyBAxyxy

(,)xy分别是矩阵的行列索引值。 离散正弦变换： 给定向量011(,,...,)Nxxxx，其正弦变换011(,,...,)NXXXX为：

正变换：101)1)(1(sinNnnkNknxX

逆变换：101)1)(1(sin12NkknNknXNx （二） OpenCV 1.0安装手册 1) 安装OpenCV 从http://www.opencv.org.cn下载OpenCV安装程序。假如要将OpenCV安装到C:\Program Files\OpenCV。（图2，图3为OpenCV 1.0rc1的安装界面，OpenCV 1.0安装界面与此基本一致。）在安装时选择"将\OpenCV\bin加入系统变量"（Add\OpenCV\bin to the systerm PATH）。

Div(G) 2D IDST 计算特征值 2D DST

G(x,y) Bij Aij

重建图像 图2 软件安装目录选择

图3 软件安装路径设置 2) 配置Windows环境变量 检查C:\Program Files\OpenCV\bin是否已经被加入到环境变量PATH，如果没有，请加入。加入后需要注销当前Windows用户（或重启）后重新登陆才生效。(可以在任务管理器里重启explorer.exe) 图4 软件安装环境变量设置 3) 配置Visual C++ 6.0 全局设置：菜单Tools->Options->Directories：先设置lib路径，选择Library files，在下方填入路径：

C:\Program Files\OpenCV\lib 然后选择include files，在下方填入路径： C:\Program Files\OpenCV\cxcore\include C:\Program Files\OpenCV\cv\include C:\Program Files\OpenCV\cvaux\include C:\Program Files\OpenCV\ml\include C:\Program Files\OpenCV\otherlibs\highgui C:\Program Files\OpenCV\otherlibs\cvcam\include

图5 VC6.0开发软件环境路径设置