低照度彩色图像的自适应亮度增强

低照度彩色图像的自适应亮度增强 Adaptive brightness correction of dim-lightening color images 1 first-author 陈　刚（1995-），男，江苏南京人，硕士研究生，研究方向是光电成像与图像处理。E-mail： 1219023223@njupt.edu.cn http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19370079&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19370083&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19370082&type= 22.01332855 32.00399780 陈　刚 （1995-），男，江苏南京人，硕士研究生，研究方向是光电成像与图像处理。E-mail： 1219023223@njupt.edu.cn 1219023223@njupt.edu.cn 1 1 2 1 corresp E-mail： yucy@njupt.edu.cn E-mail： yucy@njupt.edu.cn 喻春雨（1976-），女，辽宁沈阳人，副教授，2006年于南京理工大学获光学工程博士学位，2006~2008年在北京大学信息科学技术学院从事博士后工作，2014~2015年在宾夕法尼亚大学医学院访学。E-mail： yucy@njupt.edu.cn http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19370088&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19370092&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19370089&type= 22.01332855 32.00400543 喻春雨 （1976-），女，辽宁沈阳人，副教授，2006年于南京理工大学获光学工程博士学位，2006~2008年在北京大学信息科学技术学院从事博士后工作，2014~2015年在宾夕法尼亚大学医学院访学。E-mail： yucy@njupt.edu.cn yucy@njupt.edu.cn 南京邮电大学 电子与光学工程学院、微电子学院，江苏 南京 210023 南京邮电大学 电子与光学工程学院、微电子学院，江苏 南京 210023 Nanjing Univ Posts & Telecommun，College of Electronic and Optical Engineering & College of Microelectronics， Nanjing 210023， China Nanjing Univ Posts & Telecommun，College of Electronic and Optical Engineering & College of Microelectronics， Nanjing 210023， China 南京邮电大学 自动化院，江苏 南京 210023 南京邮电大学 自动化院，江苏 南京 210023 Nanjing Univ Posts & Telecommun， School of Automation， Nanjing 210023， China Nanjing Univ Posts & Telecommun， School of Automation， Nanjing 210023， China 为恢复低照度场景图像的原有色彩信息和细节信息，本研究提出自适应亮度调节的低照度彩色图像增强方法。该方法首先对低照度场景连续拍摄多帧图像，并对其进行自适应伽马亮度校正；然后将多帧亮度校正后图像转换到YUV色彩空间并行两种处理，一种是提取Y通道分量分组进行基于权值调整二阶盲辨识的盲源分离降噪，一种是进行帧平均后提取Y通道分量与多个盲源分离降噪的Y通道分量依次进行结构匹配，并选出匹配最佳Y通道分量；再将最佳Y通道分量进行基于皮尔生长曲线的亮度调整后与经帧平均处理的U、V通道分量重新组合；最后将重组图像转换回RGB空间，即可得到视觉效果显著提升的彩色图像。本文提出的图像增强方法满足实时处理要求，可将原彩色图像的极低亮度提高54.4倍、中等亮度提高3.5倍；并将图像信息熵提高1.3~2.9倍。与典型的图像增强算法相比，本文提出的方法对低照度彩色图像在降低噪声、均衡光照和恢复细节方面有较大改善。 To improve the imaging performance in dim-lighting scenes， we propose adaptive brightness correction of dim-lighting color images. First， frame-averaging and Gamma correction are performed on a sequence of images. Second， the processed images are grouped before the following parallel processing steps are performed. In one of these steps， the images are converted into the YUV color space， and then， based on blind source separation （BSS）， noise is removed from Y-channel components group by group； in another step， the best BSS-denoised Y-channel component is selected by matching with the frame-averaged Y-channel. Third， the best Y-channel component is adjusted based on the peel growth curve and then recombined with the averaged U-channel and V-channel components. Finally， the recombined image is converted back to the RGB color space. Experiment results show that the proposed image enhancement algorithm could increase the brightness of the low-luminosity image by 3.5-54.4 times and increase the information entropy by 1.3-2.9 times. The proposed algorithm outperforms the classical image enhancement algorithms in terms of noise reduction， brightness balancing， and color information restoration. 低照度 伽马校正 色彩空间 盲源分离 皮尔生长曲线 dim-lightening gamma correction color space blind source separation peel growth curve 1 引　言 1 由于低照度场景中所拍摄的图像具有亮度低、噪声强、色彩差、细节信息缺失的特点，所以还原与增强低照度场景图像的研究吸引了越来越多的关注。目前，该方面的研究工作很多依赖于专用设备，如通过微光像增强器 ［ 1 1 ］ 、红外热成像仪等对低照度场景进行成像并辅以图像处理工作。但这些专用设备不仅价格高，而且其工作范围主要在可见光和近红外波段，对更低光照度场景拍摄图像效果并不理。因此，在提高专用成像设备性能并降低其成本的同时，研究不依赖于专用设备的低照度场景图像增强算法是非常有研究意义的。 常用的图像增强算法主要涉及直方图均衡、Retinex算法、去雾算法和伽马校正。一些经典的图像增强算法均具有其自身优势及局限所在，如：Zuiderveld ［ 2 2 ］ 提出的限制对比度自适应直方图均衡（Contrast Limited Adaptive Histogram Equalization， CLAHE）能够很好解决图像对比度低的问题，但存在图像对比度过度增强、颜色失真的问题；Jobson等人 ［ 3 3 ］ 提出了一种带有色彩恢复的多尺度Retinex算法（Multi-Scale Retinex with Color Restoration， MSRCR），抑制了图像明暗对比强烈的边缘交界处易产生“光晕伪影”的情形，但在低对比度区域会丢失细节信息；Dong等提出了一种基于去雾模型的低照度图像增强方法 ［ 4 4 ］ ，主要利用暗通道反转类似于有雾图像的特性，反转暗通道图像去雾后再转回实现图像增强，但图像轮廓尤其高亮度区域的轮廓出现严重伪影；Guo等提出了一种基于光照图估计的低照度图像增强算法LIME（Low-Light Image Enhancement via Illumination Map Estimat-ion） ［ 5 5 ］ ，先提取原低照度图像RGB三个通道的最大值得到初始光照图，再通过先验结构和伽马校正修正初始光照图，方法的缺点是处理后图像存在颜色过饱和；还有利用自适应伽马校正改善光照不均匀对图像的影响，其中增强效果较好的是刘志成等人提出的一种基于二维伽马函数的光照不均匀图像自适应校正算法 ［ 6 6 ］ ，利用多尺度高斯函数提取场景的光照分量来构造二维伽马函数，并利用光照分量的分布特性调整二维伽马函数的参数，有效降低了光照过强区域的图像亮度，提高光照过暗区域的图像亮度，但在亮度归一化后值为1的区域无法得到校正，且对场景亮度极低的图像增强效果不理想；近年来，Zohair提出了一种结合简单对数、指数运算、双曲正割累积分布函数和归一化的图像增强算法 ［ 7 7 ］ ，能够有效防止图像对比度过度增强，具有良好的视觉效果，但该算法不能自适应增强，对不同亮度动态范围的需要人工干预。 以上所述为典型的图像亮度增强方法，均是针对拍摄的单张图像进行处理，因此增强前进行的图像降噪也受单帧图像处理方法的局限，而低照度图像往往含有较强程度的噪声，进行高效降噪有助于提升图像亮度增强效果。基于以上分析，本文提出一种基于自适应亮度调节的低照度彩色图像增强方法，为克服对单帧图像进行增强在算法和效果上的局限，该方法采用的盲源分离方法（Blind Source Separation，BSS）和帧平均（Frame Averaging， FA）两种基于多帧的处理方法对低照度彩色图像在增强过程中进行降噪，所选用的BSS方法是对时序信号处理敏感、运算快的权值调整二阶盲辨识（Weight-Adjusted Second-order Blind Identification，WASOBI） ［ 8 8 ］ ，并对降噪后图像进行自适应参数的亮度伽马校正及基于皮尔生长曲线 ［ 9 9 ］ 的亮度矫正以避免图像局部亮度过增强，实验分析表明：该方法可以有效避免对比度过大和局部区域亮度饱和问题，使整幅图像亮度均匀、能够有效恢复原场景图像中丰富的色彩信息和细节信息。 2 方法原理与实验步骤 1 2.1　自适应伽马校正 2 本文提出的图像增强算法是基于多帧图像进行处理的，为避免单一的伽马校正对不同帧图像会引入光波动影响，因此采用自适应伽马校正方法进行亮度校正，其数学表达如下： <math display="block" id="M1"><msub><mrow><mi>I</mi></mrow><mrow><mi mathvariant="normal">o</mi><mi mathvariant="normal">u</mi><mi mathvariant="normal">t</mi></mrow></msub><mo>=</mo><mi>c</mi><mo>×</mo><mn mathvariant="normal">255</mn><mo>×</mo><msup><mrow><mfenced separators="|"><mrow><mi>F</mi></mrow></mfenced></mrow><mrow><mi>γ</mi></mrow></msup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19370098&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19370095&type= 32.25799942 4.74133301 ， 其中： <math display="block" id="M2"><mi>F</mi><mo>=</mo><mfrac><mrow><msub><mrow><mi>I</mi></mrow><mrow><mi mathvariant="normal">i</mi><mi mathvariant="normal">n</mi></mrow></msub></mrow><mrow><mn mathvariant="normal">255</mn></mrow></mfrac><mo>,</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19370103&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19370101&type= 16.59466553 7.53533268 <math display="block" id="M3"><mi>γ</mi><mo>=</mo><mi mathvariant="normal">l</mi><mi mathvariant="normal">o</mi><msub><mrow><mi mathvariant="normal">g</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub><mfenced separators="|"><mrow><mi>m</mi><mo>-</mo><msub><mrow><mi>F</mi></mrow><mrow><mi mathvariant="normal">m</mi><mi mathvariant="normal">e</mi><mi mathvariant="normal">d</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">a</mi><mi mathvariant="normal">n</mi></mrow></msub></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19370109&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19370106&type= 32.08866501 4.82600021 ， 式（1） 式（1） 中， <math id="M4"><msub><mrow><mi>I</mi></mrow><mrow><mi mathvariant="normal">i</mi><mi mathvariant="normal">n</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19370164&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19370162&type= 2.87866688 3.72533321 为输入图像亮度； <math id="M5"><msub><mrow><mi>I</mi></mrow><mrow><mi mathvariant="normal">o</mi><mi mathvariant="normal">u</mi><mi mathvariant="normal">t</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19370123&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19370120&type= 3.97933316 3.72533321 表示输出图像亮度； <math id="M6"><mi>c</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19370127&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19370125&type= 1.52400005 2.79399991 为比例系数，通常取1， <math id="M7"><mi>F</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19370131&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19370129&type= 2.53999996 2.79399991 表示 <math id="M8"><msub><mrow><mi>I</mi></mrow><mrow><mi mathvariant="normal">i</mi><mi mathvariant="normal">n</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19370144&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19370141&type= 2.87866688 3.72533321 的归一化亮度； <math id="M9"><msub><mrow><mi>F</mi></mrow><mrow><mi mathvariant="normal">m</mi><mi mathvariant="normal">e</mi><mi mathvariant="normal">d</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">a</mi><mi mathvariant="normal">n</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19370139&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19370137&type= 8.46666718 3.72533321 为取 <math id="M10"><msub><mrow><mi>I</mi></mrow><mrow><mi mathvariant="normal">i</mi><mi mathvariant="normal">n</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19370144&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19370141&type= 2.87866688 3.72533321 中值； <math id="M11"><mi>γ</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19370149&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19370146&type= 1.94733346 3.64066648 为伽马校正指数，其值由目标亮度均值 m 确定； m 取值如下： <math display="block" id="M12"><mi>m</mi><mo>=</mo><mfenced open="{" close="" separators="|"><mrow><mtable columnalign="left"><mtr><mtd><mn mathvariant="normal">0.4</mn><mo>,</mo><mtext>  </mtext><mi mathvariant="normal">i</mi><mi mathvariant="normal">f</mi><mtext> </mtext><mover><mrow><mi>F</mi></mrow><mo>¯</mo></mover><mo>≤</mo><mn mathvariant="normal">0.01</mn></mtd></mtr><mtr><mtd><mfrac><mrow><mi>r</mi><mi>o</mi><mi>u</mi><mi>n</mi><mi>d</mi><mfenced separators="|"><mrow><mover><mrow><mi>F</mi></mrow><mo>¯</mo></mover><mi mathvariant="normal">*</mi><mn mathvariant="normal">1000</mn></mrow></mfenced></mrow><mrow><mn mathvariant="normal">100</mn></mrow></mfrac><mo>+</mo><mn mathvariant="normal">0.3</mn></mtd></mtr><mtr><mtd><mn mathvariant="normal">0.5</mn><mo>,</mo><mtext> </mtext><mi mathvariant="normal">i</mi><mi mathvariant="normal">f</mi><mtext> </mtext><mover><mrow><mi>F</mi></mrow><mo>¯</mo></mover><mo>≥</mo><mn mathvariant="normal">0.02</mn><mtext> </mtext></mtd></mtr></mtable></mrow></mfenced><mo>,</mo><mtext> </mtext><mi mathvariant="normal">i</mi><mi mathvariant="normal">f</mi><mtext> </mtext><mn mathvariant="normal">0.01</mn><mo><</mo><mover><mrow><mi>F</mi></mrow><mo>¯</mo></mover><mo><</mo><mn mathvariant="normal">0.02</mn></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19370154&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19370152&type= 74.33734131 20.91266632 ， 式（2） 式（2） 中， <math id="M13"><mover><mrow><mi>F</mi></mrow><mo>¯</mo></mover></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19370159&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19370156&type= 2.53999996 3.38666677 为 <math id="M14"><msub><mrow><mi>I</mi></mrow><mrow><mi mathvariant="normal">i</mi><mi mathvariant="normal">n</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19370164&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19370162&type= 2.87866688 3.72533321 的均值。 图1 图1 为当输入图像帧数为1~200时，分别经过固定参数伽马校正和自适应参数伽马校正的图像亮度均值对比，由对比可见：经自适应参数伽马校正后，图像亮度波动较小，有效降低由光照条件、曝光时间等不稳定的成像因素造成的图像亮度波动。 2.2　色彩空间变换 2 由于基于RGB色彩空间的图像处理需要分别在R，G，B三个通道进行，占用较大计算资源，这里将由1.1获得的初步亮度调节图像在YUV色彩空间进行图像降噪 ［ 10 10 ］ ，RGB色彩空间到YUV色彩空间的转换公式为： <math display="block" id="M15"><mfenced open="{" close="" separators="|"><mrow><mtable columnalign="left"><mtr><mtd><mi>Y</mi><mo>=</mo><mn mathvariant="normal">0.299</mn><mi>R</mi><mo>+</mo><mn mathvariant="normal">0.587</mn><mi>G</mi><mo>+</mo><mn mathvariant="normal">0.114</mn><mi>B</mi></mtd></mtr><mtr><mtd><mi>U</mi><mo>=</mo><mo>-</mo><mn mathvariant="normal">0.147</mn><mi>R</mi><mo>-</mo><mn mathvariant="normal">0.289</mn><mi>G</mi><mo>+</mo><mn mathvariant="normal">0.437</mn><mi>B</mi></mtd></mtr><mtr><mtd><mi>V</mi><mo>=</mo><mn mathvariant="normal">0.615</mn><mo>-</mo><mn mathvariant="normal">0.515</mn><mi>G</mi><mo>-</mo><mn mathvariant="normal">0.100</mn><mi>B</mi></mtd></mtr></mtable></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19370174&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19370172&type= 56.13399887 13.12333298 ， 式（3） 式（3） 中， Y 为亮度； U ， V 为色度。 2.3　基于YUV色彩空间的图像降噪 2 将由2.2获得的YUV色彩空间图像并行进行两种处理：一种是帧平均降噪；另一种是提取代表亮度信息的Y通道分量后，对其分组进行基于盲源分离BSS的分离降噪。而UV是表示图像色彩的恒常量，且值的动态范围较小，不适合进行基于盲源分离BSS的分离降噪。 2 . 3 . 1　盲源分离降噪 3 对多个Y通道分量进行盲源分离降噪如下： <math display="block" id="M16"><msub><mrow><mi>Y</mi></mrow><mrow><mi>x</mi></mrow></msub><mfenced separators="|"><mrow><mi>t</mi></mrow></mfenced><mo>=</mo><mi>A</mi><msub><mrow><mi>Y</mi></mrow><mrow><mi>s</mi></mrow></msub><mfenced separators="|"><mrow><mi>t</mi></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19370176&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19370175&type= 24.89200020 4.06400013 ， 式（4） 式（4） 中， <math id="M17"><msub><mrow><mi>Y</mi></mrow><mrow><mi>x</mi></mrow></msub><mfenced separators="|"><mrow><mi>t</mi></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19370178&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19370177&type= 8.55133343 4.06400013 是作为混合观测信号的多个Y通道分量； <math id="M18"><msub><mrow><mi>Y</mi></mrow><mrow><mi>s</mi></mrow></msub><mfenced separators="|"><mrow><mi>t</mi></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19370180&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19370179&type= 8.12800026 4.06400013 是源信号的多个Y通道分量； <math id="M19"><mi>A</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19370183&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19370181&type= 2.87866688 2.79399991 是未知混合系数矩阵 ［ 11 11 ］ ，定义观测信号的时延协方差矩阵： <math display="block" id="M20"><msub><mrow><mi>R</mi></mrow><mrow><msub><mrow><mi>Y</mi></mrow><mrow><mi>x</mi></mrow></msub></mrow></msub><mfenced separators="|"><mrow><mi>τ</mi></mrow></mfenced><mo>=</mo><mi>A</mi><msub><mrow><mi>R</mi></mrow><mrow><msub><mrow><mi>Y</mi></mrow><mrow><mi>s</mi></mrow></msub></mrow></msub><mfenced separators="|"><mrow><mi>τ</mi></mrow></mfenced><msup><mrow><mi>A</mi></mrow><mrow><mi mathvariant="normal">T</mi></mrow></msup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19370185&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19370184&type= 33.27399826 4.48733330 ， 其中： <math display="block" id="M21"><msub><mrow><mi>R</mi></mrow><mrow><msub><mrow><mi>Y</mi></mrow><mrow><mi>s</mi></mrow></msub></mrow></msub><mfenced separators="|"><mrow><mi>τ</mi></mrow></mfenced><mo>=</mo><mi>E</mi><mo stretchy="false">[</mo><msub><mrow><mi>Y</mi></mrow><mrow><mi>s</mi></mrow></msub><mfenced separators="|"><mrow><mi>t</mi></mrow></mfenced><mtext> </mtext><msup><mrow><msub><mrow><mi>Y</mi></mrow><mrow><mi>s</mi></mrow></msub></mrow><mrow><mi mathvariant="normal">T</mi></mrow></msup><mfenced separators="|"><mrow><mi>t</mi><mo>+</mo><mi>τ</mi></mrow></mfenced><mo stretchy="false">]</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19370187&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19370186&type= 48.93733215 4.48733330 ， 式（5） 式（5） 中， <math id="M22"><msub><mrow><mi>R</mi></mrow><mrow><msub><mrow><mi>Y</mi></mrow><mrow><mi>s</mi></mrow></msub></mrow></msub><mfenced separators="|"><mrow><mi>τ</mi></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19370190&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19370188&type= 9.82133293 4.48733330 为Y通道分量的时延协方差矩阵， <math id="M23"><msub><mrow><mi>R</mi></mrow><mrow><msub><mrow><mi>Y</mi></mrow><mrow><mi>x</mi></mrow></msub></mrow></msub><mfenced separators="|"><mrow><mi>τ</mi></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19370195&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19370193&type= 10.07533360 4.48733330 为时延为 <math id="M24"><mi>τ</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19370198&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19370196&type= 1.60866666 2.79399991 的Y通道混合观测矩阵的时延协方差矩阵。根据Yeredor提出的将联合对角化转换为非线性权值最小化二乘的方法求解权值矩阵 <math id="M25"><mi>W</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19370202&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19370200&type= 3.80999994 2.87866688 （ <math id="M26"><mi>W</mi><mo>=</mo><msup><mrow><mi>A</mi></mrow><mrow><mo>-</mo><mn mathvariant="normal">1</mn></mrow></msup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19370207&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19370204&type= 14.73200035 3.04800010 ） ［ 12 12 ］ 。由此，初步估计降噪图序列为： <math display="block" id="M27"><mi>Y</mi><mfenced separators="|"><mrow><mi>t</mi></mrow></mfenced><mo>=</mo><mi>W</mi><msub><mrow><mover accent="true"><mi>Y</mi><mo>̂</mo></mover></mrow><mrow><mi>x</mi></mrow></msub><mfenced separators="|"><mrow><mi>t</mi></mrow></mfenced><mo>+</mo><mi>W</mi><msub><mrow><mover accent="true"><mi>Y</mi><mo>¯</mo></mover></mrow><mrow><mi>x</mi></mrow></msub><mfenced separators="|"><mrow><mi>t</mi></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19370212&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19370210&type= 42.84133148 4.57200003 ， 其中： <math display="block" id="M28"><msub><mrow><mover accent="true"><mi>Y</mi><mo>̂</mo></mover></mrow><mrow><mi>x</mi></mrow></msub><mfenced separators="|"><mrow><mi>t</mi></mrow></mfenced><mo>=</mo><mi>A</mi><msub><mrow><mi>Y</mi></mrow><mrow><mi>x</mi></mrow></msub><mfenced separators="|"><mrow><mi>t</mi></mrow></mfenced><mo>-</mo><mi>A</mi><msub><mrow><mover accent="true"><mi>Y</mi><mo>¯</mo></mover></mrow><mrow><mi>x</mi></mrow></msub><mfenced separators="|"><mrow><mi>t</mi></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19370215&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19370214&type= 41.90999985 4.57200003 ， 式（6） 式（6） 中， <math id="M29"><mi>Y</mi><mfenced separators="|"><mrow><mi>t</mi></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19370219&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19370217&type= 7.53533268 3.97933316 为初步估计的降噪图序列； <math id="M30"><msub><mrow><mover accent="true"><mi>Y</mi><mo>̂</mo></mover></mrow><mrow><mi>x</mi></mrow></msub><mfenced separators="|"><mrow><mi>t</mi></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19370223&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19370221&type= 8.55133343 4.57200003 表示Y通道中心化后的混合观测信号； <math id="M31"><msub><mrow><mover accent="true"><mi>Y</mi><mo>¯</mo></mover></mrow><mrow><mi>x</mi></mrow></msub><mfenced separators="|"><mrow><mi>t</mi></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19370227&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19370225&type= 8.55133343 4.23333359 表示Y通道混合观测信号的均值。 2 . 3 . 2　图像质量评价排序 3 将由2.3获得的YUV色彩空间图像帧平均降噪后，提取Y通道分量作为参考分量，和由2.3.1得到的多个盲源分离降噪的Y通道分量同时进行系数为 <math id="M32"><mn mathvariant="normal">0.1</mn><mo>≤</mo><mi>r</mi><mo>≤</mo><mn mathvariant="normal">1</mn></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19370231&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19370229&type= 17.78000069 3.30200005 的尺寸缩放，并依次计算参考分量和每个Y通道分量的结构相似性（Structural Similarity Index Metric，SSIM） ［ 13 13 ］ ，选出结构最匹配的Y通道分量。其中，系数r根据原图像尺寸确定，文中图像尺寸为480×640，缩放系数r取值0.1。 2.4　结合皮尔生长曲线的亮度调节 2 对由2.3.2选出的Y通道分量采用对比度受限的自适应直方图均衡（Contrast-Limited Adaptive Histogram Equalization，CLAHE） ［ 14 14 ］ 进行动态范围扩展，并将其处理后和处理前进行融合系数为0.5的加权融合 ［ 15 15 ］ 。为避免融合图像在高亮度区域的过度增强，对该区域进行基于皮尔生长曲线的亮度校正，其数学模型如下： <math display="block" id="M33"><msub><mrow><mi>I</mi></mrow><mrow><mi mathvariant="normal">o</mi><mi mathvariant="normal">u</mi><mi mathvariant="normal">t</mi></mrow></msub><mfenced separators="|"><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></mfenced><mo>=</mo><mn mathvariant="normal">255</mn><mo>×</mo><mfenced separators="|"><mrow><mfrac><mrow><msub><mrow><mi>I</mi></mrow><mrow><mi mathvariant="normal">i</mi><mi mathvariant="normal">n</mi></mrow></msub><mfenced separators="|"><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></mfenced></mrow><mrow><mn mathvariant="normal">255</mn></mrow></mfrac><mo>-</mo><mfrac><mrow><mi>K</mi></mrow><mrow><mn mathvariant="normal">1</mn><mo>+</mo><mi>α</mi><msup><mrow><mi>e</mi></mrow><mrow><mo>-</mo><mi>β</mi><mi mathvariant="normal">Δ</mi></mrow></msup></mrow></mfrac></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19370235&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19370233&type= 61.21400070 11.00666618 ， 式（7） 式（7） 中， I in （ i ， j ）为输入图像即加权融合处理后的图像亮度； I out （ i ， j ）为输出图像亮度； <math id="M34"><mfrac><mrow><mi>K</mi></mrow><mrow><mn mathvariant="normal">1</mn><mo>+</mo><mi>α</mi><msup><mrow><mi>e</mi></mrow><mrow><mo>-</mo><mi>β</mi><mi mathvariant="normal">Δ</mi></mrow></msup></mrow></mfrac></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19370250&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19370248&type= 15.15533447 7.87400007 为皮尔生长曲线模型； α ， β ， K 为大于0的参数，Δ=1 <math id="M35"><mo>-</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19370243&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19370241&type= 3.72533321 2.79399991 I in （ i ， j ）。 以 <math id="M36"><mi mathvariant="normal">Δ</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19370246&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19370245&type= 2.53999996 2.79399991 为输入值， <math id="M37"><mfrac><mrow><mi>K</mi></mrow><mrow><mn mathvariant="normal">1</mn><mo>+</mo><mi>α</mi><msup><mrow><mi>e</mi></mrow><mrow><mo>-</mo><mi>β</mi><mi mathvariant="normal">Δ</mi></mrow></msup></mrow></mfrac></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19370250&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19370248&type= 15.15533447 7.87400007 为输出值，则： <math display="block" id="M38"><mfrac><mrow><mi>K</mi></mrow><mrow><mn mathvariant="normal">1</mn><mo>+</mo><mi>α</mi><msup><mrow><mi>e</mi></mrow><mrow><mo>-</mo><mi>β</mi><mi mathvariant="normal">Δ</mi></mrow></msup></mrow></mfrac><mo>=</mo><mfenced open="[" close="" separators="|"><mrow><mtable><mtr><mtd><mtable><mtr><mtd><mtable><mtr><mtd><mn mathvariant="normal">0.15</mn></mtd><mtd><mn mathvariant="normal">0.145</mn></mtd><mtd><mn mathvariant="normal">0.136</mn></mtd></mtr></mtable></mtd><mtd><mn mathvariant="normal">0.12</mn></mtd><mtd><mn mathvariant="normal">0.085</mn></mtd><mtd><mn mathvariant="normal">0.15</mn></mtd></mtr></mtable></mtd><mtd><mn mathvariant="normal">0.085</mn></mtd><mtd><mn mathvariant="normal">0.01</mn></mtd><mtd><mfenced open="" close="]" separators="|"><mrow><mtable><mtr><mtd><mn mathvariant="normal">0</mn></mtd><mtd><mn mathvariant="normal">0</mn></mtd></mtr></mtable></mrow></mfenced></mtd></mtr></mtable></mrow></mfenced><mo>,</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19370254&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19370253&type= 123.52866364 7.87400007 <math display="block" id="M39"><mi mathvariant="normal">Δ</mi><mo>=</mo><mfenced open="[" close="" separators="|"><mrow><mtable><mtr><mtd><mtable><mtr><mtd><mtable><mtr><mtd><mn mathvariant="normal">0</mn></mtd><mtd><mn mathvariant="normal">0.02</mn></mtd><mtd><mn mathvariant="normal">0.04</mn></mtd></mtr></mtable></mtd><mtd><mn mathvariant="normal">0.06</mn></mtd><mtd><mn mathvariant="normal">0.1</mn></mtd><mtd><mn mathvariant="normal">0.15</mn></mtd></mtr></mtable></mtd><mtd><mn mathvariant="normal">0.2</mn></mtd><mtd><mn mathvariant="normal">0.3</mn></mtd><mtd><mfenced open="" close="]" separators="|"><mrow><mtable><mtr><mtd><mn mathvariant="normal">0.4</mn></mtd><mtd><mn mathvariant="normal">0.5</mn></mtd></mtr></mtable></mrow></mfenced></mtd></mtr></mtable></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19370259&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=19370257&type= 96.09667206 4.74133301 图像经过基于皮尔生长曲线的亮度校正后的归一化亮度值如 图2 图2 所示。由图可见，图像的高亮度值在变换后被有效降低，防止图像局部亮度过饱和。 经以上处理后，再将输出图像转换回RGB色彩空间，即完成对低照度彩色图像的增强和去噪。 2.5　算法流程 2 相应于以上的算法原理和方法，本文提出的算法步骤如下： Step 1：在固定视场内连续获取 N 帧低照度彩色图像，采用帧平均法，获得 n 张初次降噪图； Step 2：对 N 张降噪图像依次采用自适应伽马校正； Step 3：将由Step 2获得的图像序列依次从RGB色彩空间转到YUV色彩空间； Step 4：对由Step 3获得的YUV色彩空间图像并行进行两种处理：一个是对彩色图像进行帧平均处理；另一个是对Y通道分量分 n 组进行基于权值调整二阶盲辨识的分离降噪； Step 5：以由Step4 获得的帧平均处理图像的Y通道分量作为参考，采用缩小图结构匹配法对Step 4获得的盲源分离降噪后 n 个Y通道分量质量排序，并提取出最优质的 Y 通道分量； Step 6：对由Step 5获得的最优质Y通道分量进行结合皮尔生长曲线的亮度调节； Step 7：将由Step 6获得的处理结果与由Step 4获得的帧平均处理的U，V通道分量重组，并转回RGB色彩空间； Step 8：对Step 7获得的RGB图像进行最大值受限的归一化处理，即完成低照度彩色图像增强。 相应于以上的算法步骤，本文提出算法的流程图如下： 3 实验与性能分析 1 3.1　算法性能分析 2 选择如图 6（a） 6（a） ~ 6（f） 6（f） 六种不同的低照度彩色场景进行算法性能分析，不同场景增强前具有不同程度的低亮度，如 图6 图6 第1行所示。兼顾运算时间和增强效果的因素，对每种场景连续拍摄200帧并进行均匀分组，分组数 n 取值为1~100。 所拍摄图像的尺寸为480×640，运算环境为Intel（R） Core（TM） i7-4720HQ @2.60 GHz 2.60 GHz四核处理器、8GB内存的笔记本电脑。对本文提出的图像增强算法采用均值和运行时间进行性能评价，前者反映图像的亮度水平，后者反映算法的运行速度。 由 图4 图4 的增强后图像亮度均值随分组数 n 的变化曲线可知：在分组数 n 取值为（1，100）时，图像亮度均值呈上升趋势；在分组数 n 取值为（20，100）时，图像亮度均值上升缓慢。 图5 图5 为分组数 n 取不同值时，本文算法对不同低照度彩色场景进行图像增强的运行时间随分组数 n 的变化曲线； 图5 图5 （a）选图像缩放系数为1； 图5 图5 （b）选图像缩放系数为0.1。由 图5 图5 可知：运行时间随分组数 n 增加而增加，且采用缩小图像进行质量评价排序可大大降低运行时间。 3.2　实验结果对比 2 综合以上图像增强效果和运行时间的分析，本文提出算法在连续拍摄200帧时选取分组数 n 为25。 图6 图6 给出不同低照度场景经5种典型算法和本文算法增强处理后的图像性能对比，5种典型增强算法分别为MSRCR算法 ［ 3 3 ］ 、Dong算法 ［ 4 4 ］ 、LIME ［ 5 5 ］ 算法、Liu算法 ［ 6 6 ］ 和Zohair ［ 7 7 ］ 算法。 图6 图6 展示的是室内、外6种不同弱光照场景下的低照度彩色图像经5种经典算法和本文提出方法的增强效果对比。 图6 图6 中第1行是低照度场景原图，从左向右场景亮度逐渐增强； 图6 图6 中第2行-第7行依次为经典方法MSRCR ［ 3 3 ］ ，Dong ［ 4 4 ］ ，LIME ［ 5 5 ］ ，Liu ［ 6 6 ］ ，Zohair ［ 7 7 ］ 和本文方法的增强处理效果。由对比可见，图像亮度在经不同方法增强后均得到提升，场景中景物的清晰度均得到改善。但对比5种典型图像增强算法和本文方法，可以看出：Dong ［ 4 4 ］ ，LIME ［ 5 5 ］ 和本文方法对图像亮度的提升较大，MSRCR ［ 3 3 ］ ，Liu ［ 6 6 ］ 和Zohair ［ 7 7 ］ 对图像亮度的提升不足，且经MSRCR ［ 3 3 ］ 增强后图像有光晕、细节扭曲现象；而对比Dong ［ 4 4 ］ ，LIME ［ 5 5 ］ 和本文方法，经Dong ［ 4 4 ］ 增强后在高亮度区的物体边缘处有明显的轮廓，且Dong ［ 4 4 ］ 和LIME ［ 5 5 ］ 增强后图像存在局部亮度过饱和、图像亮度不柔和现象；经LIME ［ 5 5 ］ 增强后图像，视觉效果较好，但颜色失真、偏绿；经Zohair ［ 7 7 ］ 增强的图像细节保留较好，但颜色失真、有亮度过增强区域。 由 图6 图6 中图像增强处理前、后所对应的直方图可见，当图像中存在如场景d，e和f中的高亮度区域时，该区域的灰度值经5种典型算法增强后均存在“过增强”，其中LIME ［ 5 5 ］ 和Zohair ［ 7 7 ］ 算法对场景f的增强凸显这一问题；而本文算法可以有效避免高亮度值的过度增强，使处理后的整幅图像亮度分布均匀，符合人眼视觉特性。 相应于 图6 图6 的视觉效果对比， 表1 表1 给出本文方法和5种典型算法对低照度彩色图像增强后的评价参数对比，其中：均值反映图像亮度的增强能力；信息熵反映图像包含的细节信息量；运行时间反映算法的实效性。此外，还采用了基于自然场景统计的通用失真无参考图像质量评价BRISQUE ［ 16 16 ］ 和衡量结构相似程度的有参考图像质量评价SSIM ［ 14 14 ］ 。BRISQUE是从图像中提取MSCN（Mean Subtracted Contrast Normalized）系数，然后将MSCN系数拟合成非对称性广义高斯分布，并提取拟合的高斯分布特征，输入到支持向量机SVM中进行回归，从而得到图像质量的评分，其值越低说明对图像信息保留越好，具体算法描述见文献 ［ 16 16 ］ ；SSIM是一种衡量两张图像相似度的指标，它从图像组成的角度反映场景中物体结构的属性，并将失真建模为亮度、对比度和结构三个不同因素的组合，SSIM值越大说明图像失真越少，具体算法描述见文献 ［ 13 13 ］ 。 Scene Method Dim-lightening MSRCR ［3］ Dong ［4］ LIME ［5］ Liu ［6］ Zohair ［ 7 7 ］ The proposed 5 Scene a Mean 1.347 2 6.777 0 8.755 5 34.753 0 12.945 2 14.699 8 73.312 2 Entropy 2.221 9 2.624 7 4.469 1 6.108 8 4.433 0 3.087 9 6.441 2 BRISQUE 40.865 3 49.075 4 36.832 6 25.646 3 47.520 8 30.969 3 14.654 9 SSIM 0.013 1 0.074 3 0.128 3 0.417 7 0.218 1 0.236 0 0.691 7 Runtime/s — 2.162 3 2.080 4 2.069 6 2.442 7 1.998 6 4.355 6 5 Scene b Mean 2.696 8 8.991 7 22.455 1 46.414 3 19.448 3 29.353 0 88.453 4 Entropy 2.929 6 3.554 4 5.576 9 6.482 2 4.958 5 3.926 1 6.178 9 BRISQUE 46.892 6 48.037 1 39.595 9 38.535 3 47.528 3 14.004 1 18.921 2 SSIM 0.031 4 0.125 5 0.374 5 0.521 8 0.347 1 0.468 8 0.694 0 Runtime/s — 2.032 4 2.110 7 1.948 0 2.278 9 2.003 8 4.376 6 5 Scene c Mean 4.229 9 18.695 3 71.173 1 78.042 6 36.884 7 75.731 8 111.203 8 Entropy 3.122 1 5.001 2 6.622 8 6.891 3 5.379 6 4.766 3 6.589 7 BRISQUE 49.664 8 52.356 2 38.473 1 38.354 9 44.984 26.947 1 37.985 SSIM 0.068 4 0.316 8 0.742 1 0.663 8 0.619 6 0.758 0 0.563 8 Runtime/s — 2.087 8 2.002 9 2.005 5 2.178 3 2.006 8 4.488 4 5 Scene d Mean 4.571 3 39.873 3 25.771 3 50.371 6 22.656 1 31.925 6 108.593 7 Entropy 2.960 1 5.021 4 5.829 3 6.584 1 5.091 8 4.038 0 6.719 6 BRISQUE 50.418 3 51.823 3 46.730 4 48.715 2 48.048 48.786 5 40.048 5 SSIM 0.042 2 0.546 1 0.362 0 0.520 0 0.363 7 0.438 6 0.855 5 Runtime/s — 2.445 4 2.376 2 2.413 6 2.473 2 2.300 0 4.562 0 5 Scene e Mean 13.391 4 45.618 1 68.827 0 68.827 0 39.758 9 87.356 8 123.404 5 Entropy 4.934 7 6.079 7 7.286 8 7.286 8 6.250 6 5.659 4 7.376 2 BRISQUE 51.998 4 42.998 6 43.482 5 43.482 5 51.349 5 38.339 1 38.645 3 SSIM 0.156 7 0.606 8 0.609 3 0.619 7 0.572 5 0.729 7 0.839 7 Runtime/s — 2.418 8 2.390 7 2.390 7 2.511 0 2.375 7 4.669 3 5 Scene f Mean 23.007 2 55.725 8 69.325 1 83.660 7 48.767 2 102.222 7 80.935 9 Entropy 5.521 7 6.253 8 7.275 0 7.105 1 6.332 9 5.721 5 7.399 2 BRISQUE 21.750 7 8 25.849 2 25.190 4 23.891 6 43.125 4 16.349 9 17.392 3 SSIM 0.220 6 0.678 3 0.574 9 0.656 2 0.646 0 0.704 8 0.744 0 Runtime/s — 2.145 3 2.065 8 2.015 8 2.356 2 2.075 0 4.283 1 由 表1 表1 中成像性能评价参数对比得出：相比较5种典型算法，对低照度场景拍摄200张序列图像，并结合亮度提升效果和运行时间选择分组数为25时，本文算法对低照度场景图像在亮度增强上占绝对优势，可将极低亮度的原彩色图像的亮度均值提高54.4倍、中等亮度的原彩色图像的亮度均值提高3.5倍；并可以将图像信息熵提高1.3-2.9倍，和典型算法中最优的LIME ［ 5 5 ］ 方法相近；由基于自然场景统计的通用失真评价参数BRISQUE比较可知，经本文算法处理后，BRISQUE降低20%~64%，使图像视觉效果更自然；在运行时间上满足实时处理要求。 4 结　论 1 本文提出的低照度彩色图像增强算法是一种基于多帧的处理方法，针对图像降噪采用盲源分离结合帧平均进行分组处理；针对图像亮度增强采用自适应参数的伽马亮度校正和基于皮尔生长曲线高亮度调节。经本文提出方法处理后。原场景中极低亮度可被提高54.42倍，是典型算法增强效果的2倍以上；原图像细节信息可被恢复2.9倍，和增强效果最好LIME ［ 5 5 ］ 方法相近；且BRISQUE降低20%~64%，SSIM提高2.37~51.8倍；算法运行可满足实时处理需要，因此本文提出了一种高效的、适用于日常低照度场景的图像增强方法。 参考文献 汪子君 ， 罗渊贻 ， 蒋尚志 ， 等 . 基于引导滤波的自适应红外图像增强改进算法 ［J］. 光谱学与光谱分析 ， 2020 ， 40 （ 11 ）： 3463 - 3467 . 基于引导滤波的自适应红外图像增强改进算法 WANG Z J ， LUO Y Y ， JIANG SH Z ， et al . Adaptive improved infrared image enhancement algorithm based on guided filtering ［J］. Spectroscopy and spectral analysis ， 2020 ， 40 （ 11 ）： 3463 - 3467 . （in Chinese） Adaptive improved infrared image enhancement algorithm based on guided filtering ZUIDERVELD K . Contrast Limited Adaptive Histogram Equalization ［J］. Graphics Gems ， 1994 ： 474 - 485 . Contrast Limited Adaptive Histogram Equalization LI M D ， LIU J Y ， YANG W H ， et al . Structure-revealing low-light image enhancement via robust retinex model ［J］. IEEE Transactions on Image processing ， 2018 ， 27 （ 6 ） ： 2828 - 2841 .‏ Structure-revealing low-light image enhancement via robust retinex model DONG X ， WANG G ， PANG Y ， et al . Fast efficient algorithm for enhancement of low lighting video ［C］. 2011 IEEE International Conference on Multimedia and Expo. IEEE ， 2011 ： 1 - 6 . Fast efficient algorithm for enhancement of low lighting video GUO X ， LI Y ， LING H . LIME： Low-light image enhancement via illumination map estimation ［J］. IEEE Transactions on image processing ， 2016 ， 26 （ 2 ）： 982 - 993 . LIME： Low-light image enhancement via illumination map estimation 刘志成 ， 王殿伟 ， 刘颖 ， 等 . 基于二维伽马函数的光照不均匀图像自适应校正算法 ［J］. 北京理工大学学报 ， 2016 ， 36 （ 2 ）： 191 - 196+214 . 基于二维伽马函数的光照不均匀图像自适应校正算法 LIU ZH CH ， WANG D W ， LIU Y ， et al . Adaptive correction algorithm for uneven illumination image based on two-dimensional gamma function ［J］. Journal of Beijing University of technology ， 2016 ， （ 2 ）： 191 - 196+214 . （in Chinese） Adaptive correction algorithm for uneven illumination image based on two-dimensional gamma function ZOHAIR A . Nighttime image enhancement using a new illumination boost algorithm ［J］. IET Image Processing ， 2019 ， 13 （ 8 ） ： 1314 - 1320 . Nighttime image enhancement using a new illumination boost algorithm YEREDOR A . Blind separation of Gaussian sources via second-order statistics with asymptotically optimal weighting ［J］. IEEE Signal Processing Letters ， 2000 ， 7 （ 7 ）： 197 - 200 . Blind separation of Gaussian sources via second-order statistics with asymptotically optimal weighting 李华夏 . 趋势预测及其应用 ［J］. 职大学报 ， 2001 ，（ 2 ）： 9 - 10+78 . 趋势预测及其应用 LI H Y . Trend prediction and its Application ［J］. Journal of Vocational University ， 2001 ，（ 2 ）： 9 - 10+78 . （in Chinese） Trend prediction and its Application 田会娟 ， 蔡敏鹏 ， 关涛 ， 等 . 基于YCbCr颜色空间的Retinex低照度图像增强方法研究 ［J］. 光子学报 ， 2020 ， 49 （ 02 ）： 173 - 184 . 基于YCbCr颜色空间的Retinex低照度图像增强方法研究 TIAN H J ， CAI M P ， GUANN TAO ， et al . Research on Retinex Low-illumination Image Enhancement Based on YCbCr color space ［J］. Acta Photonica Sinica ， 2020 ， 49 （ 2 ）： 173 - 184 . （in Chinese） Research on Retinex Low-illumination Image Enhancement Based on YCbCr color space 沈帆 ， 李翰林 ， 孙斌 ， 等 . 基于Anscombe变换的X射线图像序列盲源分离降噪 ［J］. 光学 精密工程 ， 2020 ， 28 （ 1 ）： 244 - 250 . 基于Anscombe变换的X射线图像序列盲源分离降噪 SHEN F ， LI H L ， SUN B ， et al . Blind source separation and noise reduction of X-ray image sequences based on Anscombe transform ［J］. Opt. Precision Eng. ， 2020 ， 28 （ 1 ）： 244 - 250 . （in Chinese） Blind source separation and noise reduction of X-ray image sequences based on Anscombe transform 刘邦 ， 肖涵 ， 易灿灿 . 基于核函数的二阶盲辨识的单通道信号盲分离方法研究 ［J］. 机械强度 ， 2018 ， 40 （ 5 ）： 1043 - 1049 . 基于核函数的二阶盲辨识的单通道信号盲分离方法研究 LIU B ， XIAO H ， YI C C . Blind separation method of single channel signal based on second-order blind identification based on kernel function ［J］. Mechanical strength ， 2018 ， 40 （ 5 ）： 1043 - 1049 . （in Chinese） Blind separation method of single channel signal based on second-order blind identification based on kernel function WANG J J ， CHEN P ， ZHENNG N N ， et al . Associations between MSE and SSIM as cost functns in linear decomposition with application to bit allocation for sparse coding ［J］. Neurocomputing ， 2021 ， 422 ： 139 - 149 . Associations between MSE and SSIM as cost functns in linear decomposition with application to bit allocation for sparse coding JIN-SAN K ， KI-DOO K . Dynamic Range Enhancement Using New MSF Image and CLAHE . 2020 ， 45 （ 2 ）： 264 - 273 . Dynamic Range Enhancement Using New MSF Image and CLAHE 李郁峰 ， 蒋勇 ， 陈念年 ， 等 . RGB颜色空间的夜视图像彩色融合 ［J］. 光学 精密工程 ， 2018 ， 26 （ 11 ）： 2795 - 2804 . RGB颜色空间的夜视图像彩色融合 LI Y F ， JIANG Y ， CHEN N N ， et al .. Color Fusion of Night Vision Images under RGB Color Space ［J］. Opt. Precision Eng. ， 2018 ， 26 （ 11 ）： 2795 - 2804 . （in Chinese） Color Fusion of Night Vision Images under RGB Color Space MITTAL A ， MOORTHY A K ， BOVIK A C . No-reference image quality assessment in the spatial domain ［J］. IEEE Transactions on image processing ， 2012 ， 21 （ 12 ）： 4695 - 4708 . No-reference image quality assessment in the spatial domain 10.37188/OPE.20212908.1999 TP391.9 A 1004-924X（2021）08-1999-09 陈刚,刘言,杨贺超等.低照度彩色图像的自适应亮度增强[J].光学精密工程,2021,29(08):1999-2007. 陈刚,刘言,杨贺超等.低照度彩色图像的自适应亮度增强[J].光学精密工程,2021,29(08):1999-2007. CHEN Gang,LIU Yan,YANG He-chao,et al.Adaptive brightness correction of dim-lightening color images[J].Optics and Precision Engineering,2021,29(08):1999-2007. CHEN Gang,LIU Yan,YANG He-chao,et al.Adaptive brightness correction of dim-lightening color images[J].Optics and Precision Engineering,2021,29(08):1999-2007. 国家自然科学基金资助项目（No.61801239）；中央高校基本科研业务费专项资金资助项目（No.30918014106）；南京邮电大学校企合作项目（No.2018外002， No.2019外157） 2021-01-06 2021-04-08 2021-08-15 2021-09-01 光学精密工程 8 29