1、多尺度形态学的局部对比度增强优化算法 刘莉 刘珊 钱雪飞 北华航天工业学院电子与控制工程学院 北华航天工业学院机电工程学院 中国石油管道局工程有限公司国际事业部 摘 要： 为了解决灰度图像对比度低而难以增强目标区域的问题, 对传统的局部对比度增强算法进行优化, 提出了多尺度形态学的局部对比度增强优化算法。该算法利用多尺度形态学算子, 对原有局部对比度拉伸算法进行改进, 通过建立形态学塔形和差别塔型, 最后重建图像, 以实现局部目标对比度增强。实验结果表明, 提出的算法可以很好地增强灰度图像目标区域的对比度, 应用前景广泛。关键词： 形态学; 多尺度; 局部对比度; 图像增强; 作者简介：刘莉
2、(1991-) , 女, 河北唐山人, 硕士研究生, 助教, 研究方向:动态测试与智能仪器, 电子邮箱:。作者简介：刘珊, 电子邮箱:。基金：河北省科技创新引导计划项目 17247673DOptimization algorithm for local contrast enhancement based on multi-scale morphologyLIU Li LIU Shan QIAN Xuefei Electronic and Control Engineering School, North China Institute of Aerospace Engineering; Me
3、chanical and Electrical Engineering School, North China Institute of Aerospace Engineering; Petroleum Pipeline Engineering Co.LTD International; Abstract： In order to solve the problem that it is difficult to enhance the target region of the gray image with low contrast, this paper makes some improv
4、ements on the traditional local contrast enhancement algorithm, and puts forward an Optimization algorithm for local contrast enhancement based on multi-scale morphology. To realize the local contrast enhancement of the images, The algorithm uses the multi-scale morphological operators to improve th
5、e original local contrast stretching algorithm, establishes morphological and difference pyramid tower, and finally reconstructs the image. Experimental results indicate that the proposed algorithm can enhance the contrast of the target area of gray image very well and has a wide application prospec
6、t.Keyword： morphology; multi-scale; local contrast; image enhancement; 0 引言图像增强是数字图像处理的一个重要分支, 它能改善人的视觉效果, 比如突出图像中目标物体的某些特点、提取目标的特征参数等, 这些都有利于对图像中目标的识别、跟踪和理解1。图像在形成过程中最常见的一种图像问题是对比度降低, 因此, 需要对图像进行增强处理, 即突出图像中感兴趣的部分, 减弱或去除不需要的信息, 增加图像的对比度, 从而加强有用信息, 得到一种更加实用的图像或者便于进行分析处理的图像。为了提高图像对比度, 学者们提出了很多基于图像对比度
7、的增强算法, 如全局对比度拉伸算法、直方图均衡化等2-3, 这样的算法可以使图像灰度级的全部动态范围得到充分地利用, 但是并不能有效增强某一目标区域的对比度。为此, 又有许多专家提出局部对比度增强算法, 例如小波理论4、局部对比度拉伸5、数学形态学6等。小波算法可以提取不同尺度上有用的图像特征, 但一些图像细节并不能得到有效的增强效果, 而且结果图像并不是很清晰;局部对比度拉伸能够改变某个区域的像素灰度值, 但这种变换仅仅能够改变局部灰度值的统计数据, 而不能区分一块区域上灰度值的连续变化和在一块区域中由于某个特征 (明亮或暗淡) 存在引起的灰度变化两者之间的区别, 使得图像在两类情况下都均衡
8、地拉伸对比度, 导致不需要的对比度增强结果;数学形态学, 是实现对比度增强的一类重要算子。它利用结构元素收集图像的信息, 考察图像各个部分之间的相互关系, 从而了解图像的结构特征, 进而进行某一区域的对比度增强。仅运用一个结构元素的形态学对比度算子将仅仅能增强对应于所用结构元素尺寸的尺度上的图像特征, 并不能使各个尺度大小的区域得到有效增强。为了使目标图像局部区域的对比度增强效果更好, 更容易进行目标识别, 图像更清晰, 本文提出了一种基于多尺度形态学的局部对比度增强优化算法。1 传统局部对比度增强算法传统局部对比度增强算法一般采用局部对比度拉伸算法7, 利用局部统计数据, 通过在一个预先确定
9、的包围某一像素的邻域上扩大灰度值和灰度均值间的不同, 进而改变图像的灰度值, 如式 (1) 所示:式中, k (k1) 是一个全局扩大因子, (x, y) 为某一像素, f (x, y) 为图像灰度值函数, mean l (x, y) 为全局灰度均值函数, f En (x, y) 为结果图像灰度值函数。后来, 有学者提出扩大因子也可以是一个基于相同邻域上与像素 (x, y) 相关的函数。假设在这个邻域的上灰度值方差是 l (x, y) , 那么因子 k (x, y) 可定义如式 (2) 所示:式中, 是一个用户自定义变量。上面的变换仅利用局部灰度值统计数据对图像中某一区域进行均衡的对比度拉伸,
10、 不能区分区域上灰度值的连续变化和区域中由于明亮或暗淡特征存在引起的灰度变化两者之间的区别, 使得在期望的平滑区域内产生不需要的对比度强化, 以致需要更进一步的抑制图像对比度增强的处理。数学形态学不仅可以增强图像的局部对比度, 而且能够处理空间特征存在的局部对比度问题, 因此, 可采用数学形态学算子来改进等式 (1) 。2 多尺度形态学的局部对比度增强算法2.1 算法的提出2.1.1 多尺度数学形态学算子本文研究的基础是灰度数学形态学8-10, 其基本算子主要包括四种, 腐蚀、膨胀、开运算和闭运算, 基本算子如式 (3) 至 (6) 所示:设 f 是输入图像, B 为结构元素, D f和 DB
11、分别是 f 和 B 的定义域, 则:(1) 灰度膨胀(2) 灰度腐蚀(3) 灰度开运算(4) 灰度闭运算在上述基本算子的基础上, Meyer 提出了高帽变换算子, 包括白高帽变换 (WTH) 和黑高帽变换 (BTH) , WTH 用来提取明亮图像区域, BTH 用来提取暗淡图像区域。定义如式 (7) 和 (8) 所示:(5) 白高帽变换(6) 黑高帽变换令 B0, B1, , Bn-1。代表 n 个尺度上具有相同形状和增长大小的结构元素。运用结构元素 Bi (u, v) 的多尺度的膨胀、腐蚀、开运算、闭运算以及高帽变换可表示如式 (9) 至式 (12) :2.1.2 利用多尺度形态学算子的局部
12、对比度增强设一幅灰度图像为 f (x, y) , 用灰值形态学算子可表示为式中, (fB) (x, y) 是灰度图像 f (x, y) 被灰度值结构元素 B (x, y) 的开运算, 结合式 (1) 和等式 (15) 可以得到下面的局部对比度拉伸变换:这里 k (k1) 依然是一个全局扩大因子, 该变换使明亮特征变的更亮, 进而提高了局部对比度。现在设 k=2。等式 (16) 变成令 FB (x, y) =f (x, y) - (foB) (x, y) , 即通过白高帽变换, 得到了大小比 B 小的特征。相似地, 把在尺度 i 下的明亮图像特征表示为 , 因为它包括了所有 f (x, y) 的
13、小于 i B 的特征。因此, 在点 (x, y) 处明亮图像特征的值 由于显示了在尺度 i 下明亮特征的存在, 给出了一种原始图像局部对比度的增强方法。定义:显然, i (x, y) 包含图像 f (x, y) 中比尺度 (i-1) 大比尺度 i 小的明亮特征。因此, 运用多尺度方法, 可得到明亮特征的局部对比度拉伸, 如式 (19) 所示。理论上, m 对应于结构元素, 可以大到全部图像。然而, 为了实际的应用目的, m 比较小。其次, 随着 i 的增加, i (x, y) 变为空图像的可能性更大。所以我们忽略当 im 时的 i (x, y) 。令 ki-1=ki+1 对于所有的 i 都成立
14、, 并且令k m=1, 最后得到明亮特征的局部对比度拉伸公式, 如式 (20) 所示:同理, 在多尺度黑高帽变换的基础上以类似的方式处理可以得到暗淡特征的局部对比度拉伸, 如式 (21) 所示:将两种特征的局部对比度拉伸公式进行结合, 以获得改进图像, 我们将式 (20) 和式 (21) 相加除 2 得到局部对比度增强公式, 如式 (22) 所示:常系数 0.5 被用作将明亮和暗淡两种特征分配均等或许以公正的权重。为了推广到更一般的情况, 将等式 (22) 改写为:这里所有特征, 不论暗淡还是明亮, 比尺度 n 小的都假定为图像中的噪声。2.2 算法的实现(1) 形态学塔形的建立在等式 (23
15、) 的执行中, 描述了基于局部对比度增强的特征涉及许多形态学塔形的建立, 待增强图像运用高阶相似的圆盘型结构元素进行一连串的灰值形态学开运算操作。图像的结果序列被存入一个塔形中, 塔形称作开运算塔形。同理, 闭运算塔形是由输入图像的灰度闭运算的一系列结果构成。因此在开运算 (闭运算) 塔形的第 i 个条目包含图像被结构元素 i B 开运算 (闭运算) 的结果。(2) 差别塔形的建立开运算后图像的像素灰度值小于或者等于原始图像的灰度值, 将其分别保存在形态学塔形的相应条目里, 将原始图像灰度值与开运算塔形的每个条目分别做差, 得出在特定结构元素尺度上的明亮图像特征, 并保存在另一塔形的相应条目里
16、, 这个塔形叫做差别塔形。同理, 另一个差别塔型由闭运算塔形及其相应的暗淡图像特征构成。因此, 在对应于开运算 (闭运算) 的差别塔形的第 i 个条目包含了小于或者等于 i B 的明亮 (暗淡) 特征。(3) 重建图像为了重建最终图像, 所需步骤如下:(1) 将对应于开运算操作的差别塔形的所有条目相加, 得到多尺度下的明亮特征组成的结果图像。(2) 同理, 将对应于闭运算操作的差别塔形的所有条目相加, 得到多尺度下的暗淡特征组成的结果图像。(3) 将原始图像、多尺度下的明亮特征组成的结果图像和暗淡特征组成的结果图像相结合, 得到最终的局部增强图像。“+”和“-”操作被用于三幅不同图像的对应像素
17、之间。但在一些情况下, 用等式 (26) 计算的 fEn (x, y) 值可能会超过允许的像素灰度级范围。在这种情况下, 权值 0.5 可以替换成 , 从而使少于总像素数 1%的点的灰度值被截断, 其中, 的值利用二值搜索技术可以在 0 到 0.5 之间选择。具体算法执行框图如图 1 所示。图 1 算法执行框图 Fig.1 The block diagram of executable algorithm 下载原图3 实验结果及分析3.1 实验结果为了显示提出算法的效力, 选用 Walkbridge 图像, 对其进行增强处理, 并用直方图均衡化、基于局部数据的对比度拉伸以及常系数对比度拉伸等不
18、同的方法与之作比较。本实验中, 局部窗口被固定成 1313, 选用的尺度数为 n=6, 结构元素 B 尺寸大小为 33, 常系数局部直方图拉伸算法中扩大因子被固定为 5, 利用局部统计数据的局部对比度增强中, 扩大因子为 0.8。试验在 Matlab7.1上运行11, 实验结果如图 2 所示。图 2 Walkbridge 图像增强结果 Fig.2 The enhancement results of Walkbridge 下载原图图 2 中, 基于局部数据的对比度拉伸算法一定程度上增强了对比度, 但是也产生了一定数量的噪声;而局部直方图均衡化算法的结果有严重的噪声, 完全影响了图像的应用;常系
19、数局部对比度拉伸算法增强了图像, 噪声较少, 但是与提出算法相比, 噪声和失真仍然存在, 尤其在明亮图像区域上, 存在过度增强的情况;而提出的算法可以有效增强图像, 生成最少噪声, 与原始图像的对比度的相比, 效果更好。3.2 实验结果分析本文选用图像峰值信噪比对实验中的算法作一个定量比较。峰值信噪比是一种评价图像的客观标准, 广泛用来检测图像质量, 其表示为 PSNR, 定义如式 (27) 所示:式中, (x, y) 是图像像素的坐标, f 和 f En分别表示原始图像和增强后的图像, PSNR 值越大, 表明图像增强算法的性能越好, 产生的噪声越少。通过计算几种局部算法和提出算法的 PSN
20、R 值, 将每种算法的 PSNR 平均值作比较, 如表 1 所示。从表 1 中可以看出, 提出算法的 PSNR 值最大, 这表明提出算法的图像对比度增强的性能最好, 可将其运用在不同的应用中。表 1 各个局部算法的 PSNR 值的定量比较 下载原表 4 结论本文提出了多尺度形态学的局部对比度增强优化算法, 该方法利用多尺度形态学理论, 对传统的局部对比度增强算法进行了改进, 通过建立形态学塔形和差别塔型, 最后对局部图像进行重建, 进而实现图像对比度增强。实验结果表明, 同其他局部增强算法相比, 提出的算法对图像的局部对比度增强效果最好, 产生的噪声最少。通过计算几种算法的峰值信噪比, 对几种
21、局部增强算法的效果进行定量比较得知, 提出算法对图像局部增强的性能最好。因此, 可根据需要, 将其应用在不同图像的局部小目标增强和识别中, 应用前景广泛。参考文献1李艳梅.图像增强的相关技术及应用研究D.电子科技大学, 2013. 2完颜丹丹.红外图像对比度增强算法研究D.河南科技大学, 2012. 3郑婷.基于人眼视觉感知图像对比度增强算法的研究D.电子科技大学, 2016. 4尹士畅, 喻松林.基于小波变换和直方图均衡的红外图像增强J.激光与红外, 2013, 43 (02) :225-228. 5刘德全, 崔涛, 杨雅宁.局部对比度自适应直方图均衡化图像增强的算法研究J.信息与电脑 (理
22、论版) , 2016, 12 (07) :79-80. 6安静.基于数学形态学的图像增强算法及其应用D.西北师范大学, 2016. 7周启双.图像局部增强算法研究D.长江大学, 2016. 8闫海霞.基于数学形态学的图像边缘检测和增强算法的研究D.吉林大学, 2009. 9J.Serra.Image Mathematical Morphology and its Applications to Image ProcessingM.Kluwer Academic Publishers, 1994:28-41. 10J.Serra.Image Analysis Using Mathematical MorphologyM.London:Academic Press, 1982:61-80. 11Stephen J.Chapman.MATLAB 编程M.北京:科学出版社出版, 2007:58-76.
