改进MSRCR的电力设备图像增强算法研究

改进MSRCR的电力设备图像增强算法研究作者：束江崔昊杨刘晨斐

来源：《物联网技术》2017年第10期

摘要：针对传统多尺度颜色恢复算法处理电力设备图像时存在的光晕现象、颜色失真等问题，文中提出了一种改进的retinex图像增强方法。通过傅里叶变换将原图像转换至频域后进行卷积计算获取原图像反射分量，计算原图像的亮度信息并与反射分量进行线性加权，采用直方图均衡对合成后的图像进行处理以增加图像对比度。与传统MSRCR算法相比，该算法处理后的图片不存在光晕现象与颜色失真问题，颜色恢复能力得到加强。图像评价函数指标表明，本算法在均值，标准差，信息熵，平均梯度上都表现较佳，运算速率更快，更利于电力设备图像自动化处理与分析。

关键词：图像增强；MSRCR算法；光晕现象；图像融合；线性加权

中图分类号：TP391.41 文献标识码：A 文章编号：2095-1302（2017）10-00-04

0 引言

图像增强技术是电力设备远程视频监控系统中的一种常见技术[1，2]，用以解决受光照环境不均匀和设备成像质量差产生的图像色彩暗淡、对比度小、饱和度低等问题。retinex算法是基于人类视觉系统提出的一种图像增强方法[3]，该理论提出人类视觉图像由入射分量L（x，y）和反射分量R（x，y）合成得到，通过估计入射分量L（x，y）计算出反射分量R（x，y）。在传统retinex算法的基础上，逐渐发展出单尺度retinex（SSR）算法[3]、多尺度retinex （MSR）算法[4]、带色彩恢复的多尺度retinex（MSRCR）算法[5]等方法。这类方法在实现色感一致性与动态压缩等方面具备一定的优势，但在处理饱和度低、图像灰化与明暗对比强的图片时，易出现噪声增强，光晕和迭代求解效率低等问题[6]，从而给目标特征的检测和提取造成了一定困难。

颜色失真、噪声及光晕现象是retinex图像增强处理过程中普遍存在的问题。为此，常戬等人[7]提出利用双边滤波对反射分量进行处理来解决光晕现象，占必超等人[8]提出了平稳小波变换结合retinex算法、吴一全等人[9]提出了Contourlet变换与retinex结合的算法，来解决光晕现象与颜色失真问题，李锦等人[10]提出采用直方图均衡与retinex算法加权融和的方法，Petro等人[11]提出MSRCP算法（即对图像的光度数据进行retinex处理，然后再把数据根据原始RGB的比例映射到每个通道）在保留原始颜色分布的基础上增强图像。上述方法可在一定程度上消除光晕现象，但算法的具体实现过程较复杂。由于电力设备的种类和数量众多，特别是设备群多处于野外复杂、恶劣环境中，拍摄的电力设备图像的光晕现象与颜色失真问题较为突出，这就对图像增强技术的分析能力、增强效果、批量处理以及处理效率提出了更高的要求。

本文运用频域变换与光照补偿的方式对传统MSRCR算法进行了改进，消除了MSRCR算法中的光晕现象，提高了算法效率。该方法将电力设备原图像转换到频域内进行卷积计算，达到抑制新噪声产生和消除光晕现象的目的。将处理后得到的反射分量与原图像的亮度进行线性加权实现了原图像色彩基调的还原。选用各通道的均值和方差来自适应拉伸范围，对合成后的图像进行直方图均衡[12]以增加图像对比度。实验结果表明，经本文算法增强后的电力设备图像细节清晰，较好地解决了颜色失真及光晕现象，相对于传统MSRCR算法，图像处理的运算速率得到较大提升。

1 retinex算法的改进

1.1 retinex算法基础

retinex理论认为视觉图像信号I（x，y）由环境光学分量L（x，y）和目标物体反射分量R（x，y）合成，通过估计入射分量L（x，y）计算出反射分量R（x，y） [3， 4]：

其中，F（x，y）为环绕函数，δ为高斯函数的标准值，它决定着retinex算法的增强效果。为进一步提升MSR算法的颜色恢复效果，对传统MSR算法进行改进得到MSRCR算法[5]：

其中，Ci（x，y）是色彩恢复函数，β为增益常数，α是受控制的非线性强度。利用MSRCR算法处理后的图像像素值一般会出现负值，进行修正后表示为[5]：

其中，G为增益常数，b为增益补偿，一般都为常数。

传统MSRCR的卷积计算在空域进行，在retinex理论框架内，光照是一种缓变因素而非突变因素，即照射光的变化呈现平滑性，在明暗对比度强烈的边缘区域容易产生光晕现象。

1.2 retinex算法的改进

1.2.1 光晕消除及效率提升

针对传统retinex算法存在的边缘光晕现象，本文利用快速傅里叶变换（FFT）将原图像频谱移频到圆心（FFTshift）分辨出图像频率分布，分离出包含光晕信息的低频部分并保留体现图像本质的高频信息，以此可以较好地解决光晕问题。

由公式（2）～（4）可以看出，在MSR算法中，L（x，y）的估算是利用高斯函数与原图像反复迭代计算得到，因而效率低且会引入不可调节的新噪声。采用FFT在频域内完成计算，缩小数据范围，可以使得计算过程中引入的额外噪声减小，原本的卷积计算简化为乘积计算且时间复杂度将降低logn/n倍，从而缩短运算时间。在算法的具体实现过程中，公式（4）利用图1方法进行简化：

图2以电力设备的绝缘子为例，对比了传统MSRCR与带FFT变换MSRCR（FFT-MSRCR）算法图像增强后的效果。从图2（a）可以看出，原图中因照度低产生的阴影导致设备视觉细节较为模糊。图2（b）所示为经过传统MSRCR处理后的设备图像增强效果图，尽管设备的部分细节信息较原图有所提升，但也出现了颜色失真以及绝缘子瓷瓶处的光晕现象。图2（c）所示为FFT-MSRCR算法处理后的图像增强效果，相对于原图及传统MSRCR处理后的图像而言，该图中绝缘子纹理则更为清晰，细节色彩分明且处理速度大幅提升。

但该方法处理后图像仍存在部分构件边缘处的颜色失真问题。这是由于在RGB色彩空间的图像增强中，FFT虽然能提升处理速度，消除光晕影响，但估算光照分量被剔除的本质未变，从而在一定程度上导致了颜色失真。

1.2.2 色彩保真性改进

为保留原图的颜色基调，使图像色彩恢复得到提升，本文引入HSI颜色空间的亮度来补偿被剔除的光照分量[13]。HSI色彩空间的亮度信息与RGB色彩空间的转换关系如下[14]：

S（x，y）=[r（x，y）+g（x，y）+b（x，y）]/3 （8）

其中，S（x，y）为亮度信息，r（x，y），g（x，y），b（x，y）代表原图像的三色分量。将S（x，y）与三通道计算得到的反射分量RR（x，y），RG（x，y），RB（x，y）进行权重值为r的加权融合，得到改进后的反射分量RR'（x，y），RG'（x，y），RB'（x，y）。

图像的信息熵在很大程度上受r值选取的影响，图3所示为10幅不同图像的信息熵对比图，每幅图片取4个r值。对于图像编号为3，5，8的黑白图像，r值取0.7时信息熵最大，图像包含的细节较多。其它编号为彩色图像，r值取0.4时信息熵最大，但实际处理后的图像色调偏暗。仿真试验表明，对本文所处理的电力设备图像，r取0.6较为合理。

分别计算出反射分量RR'（x，y），R'G（x，y），RB'（x，y）各数据的均值Mean和均方差Var。通过公式Min=Mean-K×Var，Max=Mean+K×Var，K的取值一般为2～3，计算各通道的Min和Max值，利用如下公式对Log[R（x，y）]的每一个值Value进行线性映射：

具体求解中同时需要增加溢出判断条件进行限定，即Value>255，R（x，y）=255；Value

为进一步增强图像对比度及视觉效果，本文在上述算法的基础上增加了直方图均衡。通过直方图均衡截取反射分量两端各1%的像素再进行压缩，增加了图像对比度。在对原电力设备图像的颜色修复工作中，先通过亮度与三色通道的反射分量进行加权叠加，然后对三色通道合成后的图像进行直方图均衡。通过此方式可以较好地保留原图的颜色基调，并有较强的对比度。

1.3 本文算法具体步骤