2.845

2023影响因子

(CJCR)

  • 中文核心
  • EI
  • 中国科技核心
  • Scopus
  • CSCD
  • 英国科学文摘

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

基于最大加权投影求解的彩色图像灰度化对比度保留算法

卢红阳 刘且根 熊娇娇 王玉皞 邓晓华

陈兴国, 俞扬. 强化学习及其在电脑围棋中的应用. 自动化学报, 2016, 42(5): 685-695. doi: 10.16383/j.aas.2016.y000003
引用本文: 卢红阳, 刘且根, 熊娇娇, 王玉皞, 邓晓华. 基于最大加权投影求解的彩色图像灰度化对比度保留算法. 自动化学报, 2017, 43(5): 843-854. doi: 10.16383/j.aas.2017.c160077
CHEN Xing-Guo, YU Yang. Reinforcement Learning and Its Application to the Game of Go. ACTA AUTOMATICA SINICA, 2016, 42(5): 685-695. doi: 10.16383/j.aas.2016.y000003
Citation: LU Hong-Yang, LIU Qie-Gen, XIONG Jiao-Jiao, WANG Yu-Hao, DENG Xiao-Hua. Maximum Weighted Projection Solver for Contrast Preserving Decolorization. ACTA AUTOMATICA SINICA, 2017, 43(5): 843-854. doi: 10.16383/j.aas.2017.c160077

基于最大加权投影求解的彩色图像灰度化对比度保留算法

doi: 10.16383/j.aas.2017.c160077
基金项目: 

国家自然科学基金 61362001

详细信息
    作者简介:

    卢红阳  南昌大学信息工程学院硕士研究生.主要研究方向为稀疏表示在图像处理的应用, 计算机视觉, 遥感数据融合和重建.E-mail:luhongyang6890@163.com

    熊娇娇  南昌大学信息工程学院硕士研究生.主要研究方向为稀疏表示在图像处理的应用, 磁共振图像重建.E-mail:xiongjiaojiao0126@163.com

    王玉皞  博士, 南昌大学信息工程学院教授.主要研究方向为信道建模和无线电测量, 软件无线电, 复杂场景感知监控, 非线性信号处理, 多媒体, 传感网络以及无线通信.E-mail:wangyuhao@ncu.edu.cn

    邓晓华  博士, 南昌大学教授.主要研究方向为空间物理卫星观测和计算机数值研究.E-mail:dengxhua@gmail.com

    通讯作者:

    刘且根  博士, 南昌大学信息工程学院副教授.主要研究方向为稀疏表示在图像处理的应用, 计算机视觉, 磁共振图像重建.E-mail:liuqiegen@ncu.edu.cn

Maximum Weighted Projection Solver for Contrast Preserving Decolorization

Funds: 

National Natural Science Foundation of China 61362001

More Information
    Author Bio:

     Master student at the School of Information Engineering, Nanchang University. Her research interest covers sparse representations theory and its applications in image processing, computer vision, remote sensing data fusion and reconstruction

     Master student at the School of Information Engineering, Nanchang University. Her research interest covers sparse representations theory and its applications in image processing, and magnetic resonance imaging (MRI) reconstruction

     Ph.D., professor at the School of Information Engineering, Nanchang University. His research interest covers channel modeling and radio measurement, software-defined radio, complex environment sensing and monitoring, nonlinear signal processing, multimedia, sensor networks, and wireless communication

     Ph.D., professor at Nanchang University. His research interest covers physical space satellite observations and computer numerical study

    Corresponding author: DENG Xiao-Hua  Ph.D., professor at Nanchang University. His research interest covers physical space satellite observations and computer numerical study
  • 摘要: 针对目前彩色图像灰度化难以充分保留原彩色图像对比度的问题,本文提出了基于最大加权投影求解的彩色图像灰度化模型及算法.首先,在好的彩色图像灰度化算法应使灰度化图像具有最大对比度的假设下,本模型提出最大加权投影的目标优化函数,并且将原始彩色图像梯度权重引入到最大化函数中,使得原彩色图像中对比度较小的区域也能够在灰度化后的图像中得到保持.每个彩色通道梯度的高斯加权系数反映灰度图像的对比度和原彩色图像的颜色顺序.其次,对所提模型使用参数离散搜索策略求解,通过对线性离散参数模型产生的候选图像进行搜索,由于只有几个算术运算,计算速度较快.最后,为评价所提出算法在复杂场景下图像灰度化对比度保持性能,本文对Cadik、CSDD和COLOR250数据集分别进行灰度化实验.定性和定量实验结果表明,所提算法相比于其他算法能较好地保留原彩色图像颜色对比度,同时具有对噪声鲁棒和运算速度快的优势.
  • 在图像处理和现代科技应用中, 伴随着数字图像处理技术的逐渐进步, 彩色图像灰度化技术越来越受到重视.首先在一些图像处理和模式识别等技术系统的预处理阶段, 例如人脸识别、运动目标跟踪等, 为了减少输入图像的信息量或者是减少后续的运算量, 都需要将彩色图像灰度化, 这样不仅能够增加后续算法的处理速度, 还能够大大提高系统综合应用实效.另外, 因为打印成本的原因, 很多课本报刊及论文等资料仍然将彩色图像转换成灰度图像打印出来.最后灰度图像在艺术美学方面也有广泛应用, 如中国水墨画渲染等.因此, 研究图像灰度化技术具有重要意义.

    上述应用需求推动着各种彩色图像灰度化方法的发展.彩色图像的灰度化过程可以看成是一个将彩色图像从三维颜色矢量空间变换到一维灰度标量空间的降维过程.在这种空间转换过程中, 不可避免的会有许多彩色图像的颜色信息因空间转换而丢失.各种算法试图利用有限范围的灰度值呈现输入彩色图像的对比度和颜色顺序等细节.因此, 为得到比较好的彩色图像灰度化效果, 应该考虑在映射彩色图像亮度信息的基础上, 尽量多地向灰度图像映射原彩色图像的颜色对比信息.

    图 1  彩色图像灰度化的应用示意图
    Fig. 1  Application diagram of color-to-gray conversion

    若原输入图像是RGB格式, 常用的彩色图像灰度化的方法是以固定权重线性地对RGB通道求和(即Matlab中的rgb2gray函数).然而, 对于一些有不平衡颜色值区域的彩色图像, 仅使用高亮度通道图像不能如实地表示它们的结构和对比度.近几年, 有关认知驱动灰度化的理论和算法取得显著性进展, 从而产生了大量新方法.这些灰度化算法整体上主要分为两大类:局部映射算法和全局映射算法.在局部算法映射中, 多维度色彩到单维度灰度映射的像素值取决于颜色的局部空间分布.局部映射类算法运算量小、计算简单、算法处理速度很快, 基于这些优点, 局部映射类的算法对于大部分的彩色图像都具有比较好的处理效果. Bala等[1]提出了一种将高频色度分量加入亮度分量的空间域方法, 在色度分量进行高通滤波, 并将滤波器输出量与亮度相关的项加权, 将结果加入到亮度分量. Neumann等[2]通过将彩色空间中的色彩和亮度对比度作为梯度对比度从而得到灰度图. Smith等[3]将图像分解成几段频率成分, 并用色度分量调整组合权重.这些局部映射算法具有有效地保护局部色差的特点, 但往往只是将一部分彩色图像信息通过空间映射的方法映射到灰度图像, 并没有比较好地平衡原彩色图像中的亮度值信息和色度值信息, 导致灰度图像只能体现原彩色图像的部分对比信息.

    在全局的映射方法中, 输入图像的所有像素采用相同的灰度值化映射, 其主要分为基于变换的降维类和像素点颜色对比的方法.对于变换的降维类的方法, 以主成分分析(Principle component analysis, PCA)变换降维为主要代表.对于像素点颜色对比的方法, 这一类算法的思想是在构造从彩色图像到灰度图像的映射函数时, 综合利用彩色图像像素点的亮度值信息和颜色对比度信息, 并尽量多地向灰度图像映射原彩色图像相邻区域的不同颜色对比信息, 从而增加灰度图像的对比度.对此Gooch等[4]提出一种基于彩色图像相邻像素点间颜色差值的灰度化算法, 首先, 通过求解彩色图像相邻像素点之间的颜色差值和亮度差值; 然后, 通过得到的颜色信息和亮度信息构造灰度级目标函数; 最后, 通过构造和求解优化方程得到灰度图像. Rasche等[5]直接对不同的颜色对集进行约束, 并且通过建立二次多项式求解出灰度图像. Kim等[6]分别从参数化和扩充梯度点集的角度对Gooch的工作进行了扩展.在国内, 上海交通大学的Du等基于显著性引导的图像灰度化算法, 他们使用区域显著性保持作为优化目标并将灰度图像参数化为依赖局部特征和不依赖局部特征两步部分的权重和进行优化求解来实现[7].浙江大学的宋明黎等主要侧重于在图像灰度化过程中对空域中各种视觉效果特征的保持进行建模, 在建立优化模型后采用参数优化的方式进行求解[8-10].全局映射算法缺点是将彩色图像信息映射到灰度图像后, 可能会出现原彩色图像中不同区域相同的颜色信息映射为相同的灰度值的情况, 还可能会出现彩色图像中不同区域不同颜色信息映射为同一灰度值的情况[11-12].

    现有的方法普遍存在计算成本高和对原图像对比度小的区域不能很好保持的缺点.为了克服这些缺陷, 研究人员重新考虑简单和传统的rgb2gray模型.具体地讲, 它假定输出的灰度图像是彩色图像中RGB通道的线性组合, 即分别代表RGB颜色通道分量.例如在经典Matlab的rgb2gray函数, 所有图像的权重都是固定的.最近, 一些研究人员通过引入一些度量对离散的候选灰度图像集合进行比较, 自适应地选择通道权重. Lu等提出基于梯度误差最小化的CP (Contrast preserving)[13]及RTCP (Real-time contrast preserving)模型[14]算法, 其中RTCP算法[14]使用66个候选值离散化线性参数模型的解空间, 然后基于能量值选取其中最小梯度误差的候选值作为最优解.用双峰分布函数自动选择合适的符号来保持图像的部分局部对比度信息, 但是这种方法依然会损失很多明显的局部对比度信息[15]. Song等[16]探讨现存方法的鲁棒性, 研究利用联合双边滤波来保持离散候选图像集合中灰度图像的多尺度对比度.为了得到比较好的彩色图像灰度化效果, 须考虑在映射彩色图像亮度信息的基础上, 尽量多地向灰度图像映射原彩色图像的颜色对比信息. Liu等[17]通过深入分析目前最为普遍的以梯度差值拟合作为准则所衍生的算法, 发现彩色图像灰度化效果依赖于图像的梯度模值, 因此提出梯度相关性相似测度模型, 并从线性参数模型中得到离散候选图像, 将求得最大梯度相关值的那个灰度图像作为最后结果.上述方法都向着更少的时间和更好的对比度保留发展, 特别是原彩色图像局部梯度的保留.本文针对目前算法计算成本高和难以保留对比度小的区域信息的特点, 提出可以使得对原图像对比度小的区域也能很好保持的最大加权函数模型.另外, 本文使用参数离散搜索策略求解, 算法计算速度较快.

    局部保留投影(Locality preserving projections, LPP)近年来在人脸识别领域应用广泛[18-21]. LPP原理是将原高维的样本特征空间投影到低维特征空间, 并且降维后的特征向量能保持原高维数据空间的局部特性[22-24].本文借鉴这一思想, 将其用在彩色图像灰度化的降维过程, 使灰度化后的图像能最大程度的保留原彩色图像的局部特征及颜色对比度.本文提出了基于最大加权投影求解的彩色图像灰度化模型及算法.在好的彩色图像灰度化算法能最大地保留图像的梯度变化的假设下, 将原始图像梯度权重引入到最大化函数中, 使得原彩色图像中对比度较小的区域也能够在灰度化后图像中得到保持.并且对输入彩色图像的每个彩色通道梯度用高斯加权系数来表示, 使得变换后灰度图像可保留原彩色图像的对比度和颜色顺序.对于计算机视觉领域, 保留灰度化后的图像的对比度比颜色顺序对人眼感知效果上更重要; 另一方面保留颜色顺序在对感知信息细节要求比较高的地方也能增加图像对比度.在图像灰度化模型中, 如何快速地进行求解通道权重参数是一个重要问题.目前一些灰度化算法利用现有的优化技术, 使算法尽可能快地收敛, 并且尽量少地依赖于参数.因此, 本文模型利用参数化离散搜索策略求解, 对线性参数模型产生的候选图像进行离散搜索.该非迭代算法由于只有几步简单的算术运算, 没有参数优化, 使实时运算成为可能.

    本节回顾目前常采用梯度信息进行彩色图像灰度化的方法.为了保留彩色图像灰度化转换的特征显著性, 常见策略是最小化输入的彩色图像和灰度图像的像素差.假设输入的彩色图像是RGB格式, ${{\delta }_{x, y}}$ ( $\left| {{\delta}_{x, y}}\right|=\sqrt{\sum\limits_{c=\{r, g, b\}}{{{({{I}_{c, x}}-{{I}_{c, y}})}^{2}}}}$ )为具有符号值的一对彩色像素的色彩对比度, 用 $g$ 表示灰度化后的图像, ${{g}_{x}}-{{g}_{y}}$ 表示像素 ${{g}_{x}}$ 和 ${{g}_{y}}$ 的灰度差值, 比较经典的 $L_2$ 范数能量函数定义如下:

    $ \begin{align} \underset{g}{\mathop{\min }}\sum\limits_{{(x, y)\in P}}{{{({{g}_{x}}-{{g}_{y}}-{{\delta }_{x, y}})}^{2}}} \end{align} $

    (1)

    灰度化后的图像 $g$ 可以是参数形式或者无参数形式, $P$ 表示包含局部和非局部的像素对集.

    Lu等[13]提出了由两个高斯函数组成的双模式对比度保留模型, 其目标函数如下:

    $ \begin{align} \underset{g}{\mathop{\min }}&-\sum\limits_{{(x, y)\in P}} \ln \{{{\alpha }_{x, y}}{{N}_{\sigma }}({{g}_{x}}-{{g}_{y}}+{{\delta }_{x, y}})+ \notag\\ & (1-{{\alpha }_{x, y}}){{N}_{\sigma }}({{g}_{x}}-{{g}_{y}}-{{\delta }_{x, y}})\} \end{align} $

    (2)

    其中, ${{N}_{\sigma}}(\cdot)$ 表示高斯分布, 表示弱彩色排序参数.上式通过双高斯模型构建像素对之间的误差能量函数来保持原始图像的对比度信息. Lu等[15]提出由有限多变量多项式函数构成的映射函数 $g$ , 即令 $g=\sum\limits_{{{{m}_{l}}\in Z}}{{{w}_{l}}{{m}_{l}}}$ , 并使用二阶多变量多项式函数集.因得到如下的参数优化模型:

    $ \begin{align} & \underset{{{w}_{l}}}{\mathop{\min }}\,-\sum\limits_{\left( x,y \right)\in P}{\ln }\{{{\alpha }_{x,y}}{{N}_{\sigma }}\left( \sum\limits_{{{m}_{l}}\in Z}{{{w}_{l}}\left( {{m}_{lx}}-{{m}_{ly}} \right)+{{\delta }_{x,y}}} \right)+ \\ & \left( 1-{{\alpha }_{x,y}} \right){{N}_{\sigma }}\left( \sum\limits_{{{m}_{l}}\in Z}{{{w}_{l}}\left( {{m}_{lx}}-{{m}_{ly}} \right)-{{\delta }_{x,y}}} \right)\} \\ \end{align} $

    (3)

    由于映射函数是由权重 $\left\{ {{w}_{l}} \right\}$ 唯一确定, 因此寻找灰度图像函数 $g$ 变成了最小化目标函数(3) 关于未知权重 $\left\{ {{w}_{l}} \right\}$ 的优化问题.式(2) 和(3) 本质上是考虑相邻像素对之间的差异.这种方法虽然能够整体上较好地保持原始图像的全局对比度信息的局部, 但彩色图像灰度化效果依赖于图像的梯度模值, 会导致局部色彩对比度信息损失.

    鉴于实时性的考虑, Lu等[4]进一步提出了对权重 $\left\{ {{w}_{l}} \right\}$ 进行离散搜索, 寻找最优解的策略.即

    $ \begin{align} \underset{g}{\mathop{\min }} -&\sum\limits_{{(x, y)\in P}}\ln \{ {{N}_{\sigma }}({{g}_{x}}-{{g}_{y}}+{{\delta }_{x, y}})+\notag\\ &{{N}_{\sigma }}({{g}_{x}}-{{g}_{y}}-{{\delta }_{x, y}})\} \end{align} $

    (4)

    式(4) 是式(2) 的松弛形式, 从式(4) 中可以看出, 算法最后选取梯度误差测量值最小时所对应的灰度图像.由于认识到权重微小的变化不会使灰度图像改变太多, 他们使用离散搜索策略近似求式(4) 的最小解.通过缩减求解空间, 并在由66个候选集组成的离散范围内搜索可能的最优解.该策略很大程度上减少候选值集合并且节约算法的运行时间.

    He等[25]提出了局部保持投影算法, 算法的主要思想为:通过降维使原始空间中相互离得近的点在低维空间中仍然保持较近的距离, 以此保留原始数据的局部结构. LPP算法是一种能将高维数据投影到低维空间上的非线性降维方法, 能使降维后的数据较好地保持数据原有的流型结构.该算法假设每一个数据点与它的邻近数据点位于流型的一个线性或者近似线性的区域内, 再将高维数据的全局非线性转换为局部线性, 而相互重叠的局部邻域便能提供全局结构信息.其具体算法如下:

    假设高维空间 ${{ R}^{n}}$ 中存在数据点集, 在低维空间中寻找一组对应点, 使得 ${{y}_{i}}$ 在某种意义上可以很好地表示. LPP目的在于通过一个转换矩阵, 把高维空间样本 ${{ R}^{n}}$ 中的样本数据集 $X$ 映射到低维空间中, 以实现对高维彩色图像三个通道数据的降维, 得到一维的灰度图像数据空间.

    $ \begin{align} Y={{W}^{\rm T}}X\in {{ R}^{d}} \end{align} $

    (5)

    高维数据 ${{x}_{i}}$ 经过投影后得到低维空间中对应的点 ${{y}_{i}}$ . LPP中, 相似度权重矩阵 ${{S}_{i, j}}$ 表示两个样本点 ${{x}_{i}}$ 和 ${{x}_{j}}$ 的局部结构相似度, 采用 $K$ 邻近法来定义:

    $ \begin{align} {{S}_{ij}}=\left\{\!\!\begin{array}{ll} \exp \left(\dfrac{ {-\left\| x_i-x_j \right\|}^{2}}{t} \right), &\mbox{若}\ x_j\in N_k\left( x_i \right) \\ 0, &\mbox{其他} \\ \end{array} \right. \end{align} $

    (6)

    其中, ${{x}_{i}}$ 在 ${{x}_{j}}$ 的 $K$ 邻近邻域中, 或者 ${{x}_{j}}$ 在 ${{x}_{i}}$ 的 $K$ 邻近邻域中, 式中 $t$ 是一个大于0的常量, 作为一个经验值.因此, LPP的目标函数为

    $ \begin{align} & {{F}_{\text{LPP}}}=\sum\limits_{i=1}^{N}{\sum\limits_{j=1}^{N}{{{\left\| {{y}_{i}}-{{y}_{j}} \right\|}^{2}}{{s}_{ij}}}}=\sum\limits_{i=1}^{N}{\sum\limits_{{{x}_{j}}\in {{N}_{K}}\left( {{x}_{i}} \right)}{}}{{\left\| {{W}^{\text{T}}}{{x}_{i}}-{{W}^{\text{T}}}{{x}_{j}} \right\|}^{2}}{{s}_{ij}}= \\ & {{\sum\limits_{i=1}^{N}{\sum\limits_{{{x}_{j}}\in {{N}_{K}}\left( {{x}_{i}} \right)}{\left\| {{W}^{\text{T}}}{{x}_{i}}-{{W}^{\text{T}}}{{x}_{j}} \right\|}}}^{2}}\exp \left( -\frac{{{\left\| {{x}_{i}}-{{x}_{j}} \right\|}^{2}}}{t} \right) \\ \end{align} $

    (7)

    其中, ${{x}_{i}}$ 和 ${{x}_{j}}$ 是接近的, 最小化 ${{F}_{LPP}}$ 是为了保证 ${{y}_{i}}$ 和 ${{y}_{j}}$ 也是接近的, 相似度系数矩阵 ${{S}_{i, j}}$ 值的范围是.如果 ${{S}_{i, j}}$ 趋近于1, 则表示 ${{x}_{i}}$ 和 ${{x}_{j}}$ 是接近的; 如果 ${{S}_{i, j}}$ 趋近于0, 则表示 ${{x}_{i}}$ 和 ${{x}_{j}}$ 是相对远离的.

    本文借鉴局部保留投影算法思想进行彩色图像灰度化, 像素点 $x$ 和 $y$ 的灰度值表示为 ${{g}_{x}}$ 和 ${{g}_{y}}$ , 其由相应的彩色像素点 ${{I}_{x}}, {{I}_{y}}$ 经过转换矩阵 $W$ 后得到转换后灰度图像的像素点 ${{g}_{x}}={{W}^{T}}{{I}_{x}}$ , ${{g}_{y}}={{W}^{T}}{{I}_{y}}$ .灰度图像梯度为.对于彩色图像(RGB格式), 每个通道的梯度图像表示为 $\mathop{\delta_{xy}}_{c=\{r, g, b\}}={I}_{c, x}-{I}_{c, y}$ .受局部保留投影算法思想的启发, 本文计算任意像素点的相似性, 最后彩色图像三个通道总的相似性为

    $ \begin{align} {{S}_{x, y}}=\sum\limits_{{c=\{r, g, b\}}}{\exp \left( -\frac{{{\left( \left\| {{I}_{c, x}}-{{I}_{c, y}} \right\| \right)}^{2}}}{t} \right)} \end{align} $

    (8)

    本文分别对彩色图像三个通道内任意像素求相似系数值, 三个通道的相似系数值之和为.最大加权投影模型计算每个通道的相似度系数而不是一次性计算整个彩色图像的相似度系数, 则最大加权投影模型目标函数如下:

    $ \begin{align} & \max \sum\limits_{\left( {{I}_{x}},{{I}_{y}}\in P \right)}{{{\left\| {{g}_{x}}-{{g}_{y}} \right\|}^{2}}{{S}_{x,y}}}=\ \ \\ & \max \sum\limits_{\left( {{I}_{x}},{{I}_{y}}\in P \right)}{{{\left\| {{g}_{x}}-{{g}_{y}} \right\|}^{2}}\sum\limits_{c=\{r,g,b\}}{\exp \left( -\frac{{{\left( \left\| {{I}_{c,x}}-{{I}_{c,y}} \right\| \right)}^{2}}}{t} \right)}} \\ & \ \text{s}\text{.t}\text{.g=}\sum\limits_{\text{c=}\{\text{r},\text{g},\text{b}\}}{{{\text{w}}_{\text{c}}}{{\text{I}}_{\text{c}}}};\sum\limits_{\text{c=}\{\text{r},\text{g},\text{b}\}}{{{\text{w}}_{\text{c}}}}\text{=1} \\ \end{align} $

    (9)

    其中, ${{g}_{x}}-{{g}_{y}}$ 表示灰度图像中像素值为 ${{g}_{x}}$ 和 ${{g}_{y}}$ 的梯度, $P$ 表示包含局部和非局部候选集的像素对集合. 表示每个彩色通道图像乘以权重系数 ${{w}_{c}}$ 后求和, 得到最后的灰度化图像(此处约束 $\sum\limits_{c=\{r, g, b\}}{{{w}_{c}}}=1$ ).

    局部保留投影算法思想和最大投影保持算法思想主要有两点区别.首先, LPP中样本点 ${{x}_{i}}$ 在 ${{x}_{j}}$ 的 $K$ 最近邻域内则有相似度系数矩阵 ${{S}_{ij}}$ , 其为, 而当 ${{x}_{i}}$ 不在 ${{x}_{j}}$ 的 $K$ 最近邻域中时, 相似权重矩阵 ${{S}_{ij}}$ 则为0.而最大投影保持算法中, 在整个通道所有像素对的相似度系数表示为, ${{x}_{i}}$ 则不一定限制在 ${{x}_{j}}$ 的 $K$ 最近邻域范围内.其次, LPP是最小化目标函数来保证降维后的 ${{y}_{i}}$ 与 ${{y}_{j}}$ 相似度最小.然而, 最大投影保持算法则是最大化目标函数, 其目的是使原图像对比度小的区域也能很好保持.本质思想是用每个通道的梯度作为先验引导, 计算原彩色图像的相似度系数矩阵 ${{S}_{ij}}$ , 同时相似度系数矩阵 ${{S}_{ij}}$ 和 $t$ 的取值也有很大相关性.如图 2所示指数变换图, 若 $t$ 固定时, 相似度矩阵 ${{S}_{x, y}}$ 取决于每个通道彩色图像任意两个像素点之间梯度.假如彩色图像中的像素点 ${{I}_{c, x}}$ 和 ${{I}_{c, y}}$ 比较接近, ${{\delta }_{xy}}$ 值小, 则相似度系数 ${{S}_{x, y}}$ 值大, 即当 $t$ 固定时其相似度系数 $ {{S}_{x, y}}$ 趋近于1.由于求解的是目标函数的最大化, 相似度系数 ${{S}_{x, y}}$ 大, 则迫使灰度图像梯度的 $\nabla g={{g}_{x}}-{{g}_{y}}$ 的变化的动力越大, 也就是说灰度图像像素点 ${{g}_{x}}$ 和 ${{g}_{y}}$ 的梯度值 $\nabla g$ 变化程度相对更大.因此, 所提算法引入的加权函数可以使得对原彩色图像对比度小的区域在灰度化后的图像中对比度大或者得到保持.此外, 假如彩色图像中的像素点和 ${{I}_{c, y}}$ 相距比较远时, 则 ${{\delta }_{xy}}$ 值大, 当 $t$ 固定时, 则两者之间的相似度系数 ${{S}_{x, y}}$ 趋近于0, 由于整个目标函数最大化, 迫使 $\nabla g={{g}_{x}}-{{g}_{y}}$ 的值变化的动力相对较小, 结果使得灰度图像像素点 ${{g}_{x}}$ 和 ${{g}_{y}}$ 的梯度值 $\nabla g$ 的变化程度相对较小.可以发现, 灰度化后图像梯度值 $\nabla g$ 变化趋势和原彩色图像梯度值的变化趋势相反, 从而使原来对比度小的像素对区域在整个目标函数中优化的比重更大, 期望达到对其更好体现的目的.

    图 2  指数变换图
    Fig. 2  Exponential transformation diagram

    为了使算法具有更好的实时性, 和文献[14]类似, 本文将连续的参数优化问题转化为将用离散的参数空间搜索的方式估计灰度图像.首先, 因为系数的变化对结果的影响不是特别大, 我们将搜索空间缩小, 三个彩色通道的权重系数为 ${{w}_{r}}, {{w}_{g}}, {{w}_{b}}$ , 分别在的空间范围内以0.1的间距分成10等份并且约束三个颜色通道各权重和为1, 将整个空间划分成个搜索状态, 并在66候选集组成的离散范围内搜索可能的最优解.其次, 像素对集合 $P$ 也是影响算法有效性的关键.本质上, $P$ 越大算法实现的性能越好, 但相应的计算时间也增加.考虑到自然图像固有的颜色冗余, Lu等下采样输入的高分辨率图像, 使图像变为 $64\times 64$ 大小.在缩小后的图像中随机的采样 ${{64}^{2}}$ 的像素对, 形成像素对集合 $P$ (本文删除小于0.07像素值).本文采取下采样方案提高算法运算速度.因此本文用有限的候选集和在目标函数最大化下选取权重参数估计模型(9).整个算法求解如下:

    最大加权投影算法(MWPDe)

    1) 穷举搜索离散66个权重候选集;

    2) 当 $j=1$ 到66时, 进行以下循环;

    3) 搜索权重 $w$ 计算式(9);

    4) 结束.

    为了定量衡量灰度化算法的对比度保留效果, 使用彩色对比度增强保持率(Color contrast preserving ratio, CCPR)作为评价指标[15]. CCPR定义如下:

    $ \begin{align} {\rm CCPR}=\frac{\#\{(x, y)|(x, y)\in \Omega, \left| {{g}_{x}}-{{g}_{y}} \right|\ge \tau \}}{\left\| \Omega \right\|}\; \end{align} $

    (10)

    其中, $\Omega $ 是其原始彩色图像梯度 ${{\delta }_{x, y}}\ge \tau $ 的集合. 是集合 $\Omega $ 中像素对的个数. 表示集合中经过灰度化后仍然保留差异的像素对个数.从式(10) 可以看出, CCPR值在0到1之间, 其值越接近1, 说明灰度化效果越好.

    目前彩色图像灰度化主要有三个测试图像集.第一个为经典的Cadik图像集, 该图像集含有各种不同的图像特征和图案, 包含24幅测试图像.第二个是Du等整理收集复杂场景的CSDD图像集, 由22幅颜色分布较复杂的图像组成.最后一个是Lu等提出的COLOR250测试图像集, 包括从Berkeley图像分割BSDS测试库和显著性SD测试库中选出的200幅自然图像以及在文档打印中经常用到的50幅与曲线图表等有关的图像.这三个图像集涵盖了简单仿真、自然场景等各种复杂图像.我们采用人员参与的主观评价和定量分析的客观评价对本文算法的灰度化效果进行分析.

    本节在各种合成的自然图像集中, 对所提出来MWPDe算法与基本的Gooch算法[4]; rgb2gray算法; Lu等提出的二阶对比度保留算法(CP)[13]及实时的对比度保留算法(RTCP)[14]对比. MWPDe中参数 $t$ 设为经验值0.01.

    Cadik数据集图像相对简单并且包含的颜色信息量少.本节对Cadik数据集用Gooch、rgb2gray、CP、RTCP和MWPD五种算法进行测试, $\tau $ 取值范围均为从1到15.从表 1可以看出MWPDe算法得到的灰度图像CCPR值比rgb2gray、CP、RTCP算法的CCPR值大, 并且随着 $\tau $ 的增大, MWPDe算法得到的灰度图像CCPR值比比较算法的优势越来越明显.从图 3观察可知本文算法不仅能够很好地保留细节特征, 也能使彩色图像的特征在灰度化后的图像中仍能辨别, 仍有很好的颜色排序保留效果. CP和RTCP算法灰度化后的图像颜色对比度相当, 而如图 3(b)(c)列Gooch算法和rgb2gray算法灰度化后的图像颜色细节特征保留少, 对比度保留差. RTCP算法尽管也是使用离散搜索求解, 但也丢失部分对比度信息.可以对比看到本文MWPDe算法在颜色特征显著的地方, 对比度保留效果较好, 如图 3中第2、7、8行:第2行中能明显地看出天蓝色和绿色的区别; 第7行MWPDe算法相比Gooch、rgb2gray、CP算法能较好地保留青提上同一颜色之间的细节对比度, 而无法区分灰度化后绿色的青提和黄色的梨; 第8行中能保留绿色叶子和红色叶子的对比度.这些示例说明MWPDe算法能较好地保留原彩色图像中对比度小区域的对比度.

    表 1  Cadik数据集的CCPR值
    Table 1  CCPR value of Cadik dataset
    τ Gooch rgb2gray CP RTCP MWPDe
    1 0.96 0.93 0.97 0.96 0.96
    2 0.93 0.88 0.94 0.94 0.94
    3 0.91 0.84 0.92 0.92 0.92
    4 0.89 0.8 0.91 0.9 0.91
    5 0.87 0.77 0.89 0.89 0.9
    6 0.85 0.75 0.87 0.87 0.88
    7 0.83 0.73 0.86 0.86 0.87
    8 0.81 0.71 0.84 0.85 0.86
    9 0.79 0.69 0.83 0.83 0.85
    10 0.77 0.68 0.82 0.82 0.84
    11 0.75 0.66 0.8 0.8 0.84
    12 0.72 0.65 0.79 0.79 0.83
    13 0.7 0.63 0.77 0.77 0.82
    14 0.69 0.62 0.76 0.76 0.81
    15 0.67 0.61 0.75 0.75 0.81
    下载: 导出CSV 
    | 显示表格
    图 3  Cadik数据集彩色图像灰度化转化结果
    Fig. 3  Color-to-gray conversion results of Cadik dataset
    ((a)输入彩色图像; (b) Gooch算法; (c) rgb2gray算法; (d) CP算法; (e) RTCP算法; (f) MWPDe算法)

    相比于Cadik数据集图像相对简单且包含的颜色信息量少, CSDD数据集图像包含更为丰富的颜色和图案信息.本节, 在CSDD数据集上对Gooch、rgb2gray、CP、RTCP和MWPDe五种算法进行测试, $\tau $ 取值范围均为从1到15.从表 2可以看出 $\tau $ 取值大于5时, MWPDe算法得到的灰度图像的CCPR值比Gooch、rgb2gray、CP、RTCP算法的CCPR值大, 但MWPDe算法相比于RTCP算法等方法的优势没有像在表 1中对于Cadik数据那么大.对大部分测试图, 如第二行红色箭头所指部分MWPDe算法效果的对比度比较明显, 另外四种算法则不能区分该区域的颜色顺序和局部的对比度. MWPDe算法对于第3、4、5行的灰度化效果相比于Gooch、rgb2gray、CP、RTCP算法效果较好, 能将颜色信息映射到对比度更明显灰度值上, 有较好的视觉效果.但也有少部分图像如第一行中红色箭头所指部分, MWPDe算法丢失了局部对比度, 而其他三种算法的视觉效果都相对较好.从图 4可以看出对于CSDD数据集场景下的复杂图像, 图像灰度化的难度更大, 可以在以后的研究中进一步分析改进MWPDe算法.

    表 2  CSDD数据集的CCPR值
    Table 2  CCPR value of CSDD dataset
    τ Gooch rgb2gray CP RTCP MWPDe
    1 0.97 0.96 0.96 0.96 0.95
    2 0.93 0.94 0.93 0.93 0.92
    3 0.91 0.91 0.91 0.91 0.91
    4 0.89 0.89 0.89 0.9 0.89
    5 0.87 0.87 0.87 0.88 0.88
    6 0.85 0.85 0.85 0.87 0.86
    7 0.83 0.83 0.83 0.85 0.85
    8 0.81 0.81 0.81 0.84 0.84
    9 0.79 0.8 0.79 0.82 0.83
    10 0.77 0.78 0.77 0.81 0.82
    11 0.76 0.76 0.75 0.8 0.81
    12 0.74 0.75 0.74 0.78 0.8
    13 0.73 0.73 0.72 0.77 0.79
    14 0.71 0.71 0.7 0.76 0.78
    15 0.7 0.68 0.7 0.75 0.77
    下载: 导出CSV 
    | 显示表格
    图 4  CSDD数据集彩色图像灰度化转化结果
    Fig. 4  Color-to-gray conversion results of CSDD dataset
    ((a)输入彩色图像; (b) Gooch算法; (c) rgb2gray算法; (d) CP算法; (e) RTCP算法; (f) MWPDe算法)

    本节使用COLOR250数据集进行实验.这个更大的数据集图像包括突出的检测数据集以及其他场景下的图像.这些图像包含很多的重要结构特征如大规模的边缘、平滑区域和细小的纹理等.在没有噪声的环境下, Gooch、rgb2gray、CP、RTCP和MWPDe五种算法对COLOR250数据集的灰度化结果的CCPR值如表 3所示, 可以看出当 $\tau$ 值大于6时MWPDe算法的CCPR值比其他四种算法高.图 5可以看出rgb2gray算法灰度化效果比较稳定, 缺乏局部和全局的对比度, 如第一行其将同一亮度但不同颜色的红色和绿色映射为接近的灰度值, 使得灰度化后很难区分两种颜色.本文MWPDe算法能较好地保留颜色对比度, 比如红色的花朵和绿色的草地, RTCP和本文MWPDe算法灰度化后的图像对比度比较明显.

    表 3  COLOR250数据集的CCPR值
    Table 3  CCPR value of COLOR250
    τ Gooch rgb2gray CP RTCP MWPDe
    1 0.95 0.96 0.96 0.96 0.95
    2 0.92 0.93 0.93 0.93 0.92
    3 0.89 0.9 0.9 0.9 0.9
    4 0.86 0.87 0.88 0.88 0.88
    5 0.83 0.85 0.86 0.86 0.86
    6 0.81 0.83 0.84 0.84 0.85
    7 0.79 0.81 0.82 0.83 0.83
    8 0.76 0.79 0.8 0.81 0.82
    9 0.74 0.77 0.79 0.8 0.81
    10 0.72 0.75 0.77 0.78 0.8
    11 0.7 0.74 0.76 0.77 0.79
    12 0.68 0.72 0.74 0.76 0.78
    13 0.67 0.7 0.73 0.75 0.77
    14 0.65 0.69 0.72 0.74 0.76
    15 0.63 0.67 0.7 0.72 0.75
    下载: 导出CSV 
    | 显示表格
    图 5  COLOR250数据集彩色图像灰度化转化结果
    Fig. 5  Color-to-gray conversion results of COLOR250 dataset
    ((a)输入彩色图像; (b) Gooch算法; (c) rgb2gray算法; (d) CP算法; (e) RTCP算法; (f) MWPDe算法)

    为了验证本文提出的MWPDe算法灰度化效果, 本节进行用户主观对比实验.借鉴文献[26]中的实验模式, 使用Cadik图像集的24幅图像, 进行用户选择性实验和准确性实验, 并且和rgb2gray、CP、RTCP算法对比.我们邀请50位年龄在20到40岁之间的参与者, 他们没有视觉缺陷并且从事不同行业.

    用户主观对比实验由两部分组成:选择性实验和准确性实验.首先, 选择性实验按如下流程进行:在没有彩色图像作为参考的情况下, 每位参与者每次从2幅灰度图像中选择他们认为相对较好的(一幅是本文算法得到的灰度图像, 另一幅是随机地从另外三种方法得到的结果图).其次, 准确性实验按如下流程进行:在有彩色图像作为参考的情况下, 观察者每次在rgb2gray、CP、RTCP和MWPDe四种算法得到的灰度图像中选择最能表示原始彩色图像的一幅.一般地, 准确性实验(带参考图像)相对于选择性实验(无参考图像)耗时间比较长.用户主观对比实验统计结果如图 6图 7所示, 从选择性实验和准确性实验来看, 提出的MWPDe算法效果相对于另外三种方法要好.从选择性实验看, 参与者对本文算法的灰度化效果有较优先的选择, 可以看到本文的算法得到的灰度化效果相对其他三种算法有较好的视觉效果, 反映出所提的算法更接近人眼视觉感知.

    图 6  选择性实验结果
    Fig. 6  Preference experiment result
    图 7  准确性实验结果
    Fig. 7  Accuracy experiment result

    本节对提出算法在有噪声情形下的灰度化效果进行分析.所有的算法中输入带均值为0, 方差为0.03的高斯噪声的图像.采用Cadik图像集进行抗噪性实验分析, 取值从1到15时各算法所得图像的CCPR值如表 4所示.可以看出RTCP和MWPDe算法的灰度化效果更好, 说明这两种方法具有良好的抗噪性能.当 $\tau $ 值从5到15, RTCP和MWPDe算法的抗噪性能相当, 原因可能是二者都采用了离散搜索策略求解.图 8是部分彩色图像经过灰度化后的结果, 可以发现本文算法得到的灰度图像视觉效果更好. CP算法的抗噪性最差, CP算法灰度化后的图像都残留有严重噪声.

    表 4  噪声情况下Cadik数据集的CCPR值
    Table 4  CCPR value of Cadik dataset with noise added
    τ rgb2gray CP RTCP MWPDe
    1 0.92 0.95 0.94 0.96
    2 0.86 0.91 0.91 0.92
    3 0.82 0.87 0.88 0.89
    4 0.78 0.84 0.85 0.86
    5 0.74 0.82 0.84 0.84
    6 0.72 0.79 0.82 0.82
    7 0.7 0.77 0.81 0.81
    8 0.68 0.75 0.8 0.79
    9 0.66 0.73 0.78 0.78
    10 0.65 0.71 0.77 0.77
    11 0.63 0.7 0.76 0.76
    12 0.62 0.68 0.75 0.75
    13 0.61 0.66 0.74 0.74
    14 0.6 0.65 0.73 0.74
    15 0.58 0.64 0.72 0.73
    下载: 导出CSV 
    | 显示表格
    图 8  噪声方差为0.03时对Cadik数据集进行抗噪性能分析
    Fig. 8  Antinoise performance analysis of Cadik dataset with 0.03 noise variance
    ((a)输入彩色图像; (b) rgb2gray算法; (c) CP算法; (d) RTCP算法; (e) MWPDe算法)

    除分析噪声对算法的影响, 我们也分析WMPDe算法对低分辨率图像的灰度化效果.我们用Cadik图像集进行低分辨率实验分析, 将Cadik图像集中的测试图像下采样2个因子, 其他的参数设置不变, 对算法进行实验.图 9是部分彩色图像经过降采样变成低分辨率图像后灰度化的结果, 可以发现相比于另外三种对比算法, 本文WMPDe算法也能较好地保持原彩色图像的对比度, 得到的灰度图像视觉效果更好.

    图 9  算法对低分辨率图像实验效果
    Fig. 9  Experiment results on the low resolution images
    ((a)输入低分辨率彩色图像(下采样2个因子); (b) rgb2gray算法; (c) CP算法; (d) RTCP算法; (e) MWPDe算法)

    由式(9) 可知, WMPDe算法只有一个参数 $t$ . $t$ 值变化引起权重函数 ${{S}_{x, y}}(t)$ 的变化, 本节测试 $t$ 值变化对图像灰度化效果的影响.如图 10(a)为 $t$ 值分别为0.01和0.00001时的相似度权重系数 ${{S}_{x, y}}(t)$ 随 $t$ 的变化, 当 $t$ 值小时, ${{S}_{x, y}}(t)$ 的变化相对比较陡峭. (b)列为输入的彩色图像, (c)列为 $t=0.01$ 时输入彩色图像的灰度化效果; (d)列为 $t=0.00001$ 时输入彩色图像的灰度化效果.对比(c)和(d)列可以得到第一幅蝴蝶在 $t=0.01$ 时边缘处对比度比较明显, 线条处的纹理信息基本丢失, 而 $t=0.00001$ 效果却相反, 边缘处对比度信息丢失, 线条处纹理信息明显, 说明 $t$ 值小时对原彩色图像对比度小的区域能有较好地保留对比度.第二幅轮子在 $t=0.01$ 时灰度化效果基本没有区分度, $t=0.00001$ 时保留了原彩色图像的亮度对比度, 能相对缓解(c)列中的不足.第三幅在 $t=0.01$ 时灰度化效果能较好地将叶子和花瓣颜色区分, 而无法区分花瓣中的颜色对比度, 在 $t=0.00001$ 保留了花瓣的颜色对比度, 但无法区分红色的花瓣和绿色的叶子.由这些现象可以看出当 $t$ 值小时且原彩色图像对比度不明显时, 即颜色有过渡处的相似系数 ${{S}_{x, y}}(t)$ 变化较大, 由于本文求解的目标函数最大化, ${{S}_{x, y}}(t)$ 则会引导灰度化后的图像在相应的区域保留相应的细节对比度, 实现局部对比度保留.

    图 10  参数敏感度分析
    Fig. 10  Sensitivity study of parameters
    ((a)相似系数 ${{S}_{x, y}}(t)$ $t$ 的变化趋势; (b)输入彩色图像; (c) $t=0.01$ 时WMPDe算法灰度化效果; (d) $t=0.00001$ 时WMPDe算法灰度化效果)

    好的彩色图像灰度化算法不仅要求灰度化后图像仍能保持原图像的对比度, 并且要求计算速度快.目前RTCP算法的计算时间最为实时, 其计算过程包括矩阵的点乘、指数和对数的计算. WMPDe算法求解过程借鉴其思想, 只需要5个矩阵点乘的计算过程, 使得其计算过程比其他的算法甚至RTCP都简洁和实时, 并且灰度化后的图像能保持原彩色图像的颜色对比度.表 5为输入 $390\times 293$ 彩色图像时不同算法的运行时间, 可以看到WMPDe算法计算时间比CP算法, RTCP算法的时间短, Gooch算法的运算时间最长, 原因是该算法需要将彩色图像从RGB颜色空间转换到Lab颜色空间, 这个空间转换的过程会增加算法的处理时间. rgb2gray算法的计算时间最短, 这也是该算法至今相对较流行的原因.总的来看, WMPDe算法的速度接近当前RTCP算法的两倍.

    表 5  输入390×293彩色图像时不同算法的运行时间
    Table 5  Runtime of different algorithms with the input 390×293 color image
    Methods Gooch rgb2gray CP RTCP MWPDe
    Runtime (s) 3.00E+04 0.0048 0.5233 0.0728 0.0343
    下载: 导出CSV 
    | 显示表格

    本文提出基于加权投影最大化的彩色图像灰度化模型及算法.针对目前其他灰度化算法不能较好地保留图像纹理色差及局部对比度的缺点, 借鉴目前模式识别领域中局部保留投影思想, 建立最大化的加权局部保留投影模型.该模型提出最大加权投影的目标优化函数, 并且将原始图像梯度的权重引入到最大化函数中, 使得原彩色图像中对比度比较小的区域在灰度化后也能得到较好地保留.并且对输入彩色图像的每个彩色通道梯度用高斯加权系数表示, 使得变换后灰度图像可保留原彩色图像的对比度和颜色顺序.进一步地, 对所提模型使用参数离散搜索方法求解, 通过对线性参数模型产生的候选图像进行离散搜索, 计算过程只需要几个算术运算, 计算速度较快.从对Cadik、CSDD和COLOR250数据集的定性和定量实验可以看出, WMPDe算法相比于其他算法能较好地保留颜色对比度, 并且具有对噪声鲁棒和计算实时的优势.

  • 图  1  彩色图像灰度化的应用示意图

    Fig.  1  Application diagram of color-to-gray conversion

    图  2  指数变换图

    Fig.  2  Exponential transformation diagram

    图  3  Cadik数据集彩色图像灰度化转化结果

    ((a)输入彩色图像; (b) Gooch算法; (c) rgb2gray算法; (d) CP算法; (e) RTCP算法; (f) MWPDe算法)

    Fig.  3  Color-to-gray conversion results of Cadik dataset

    ((a) Input images; (b) Gooch; (c) rgb2gray; (d) CP; (e) RTCP; (f) MWPDe)

    图  4  CSDD数据集彩色图像灰度化转化结果

    ((a)输入彩色图像; (b) Gooch算法; (c) rgb2gray算法; (d) CP算法; (e) RTCP算法; (f) MWPDe算法)

    Fig.  4  Color-to-gray conversion results of CSDD dataset

    ((a) Input images; (b) Gooch; (c) rgb2gray; (d) CP; \\(e) RTCP; (f) MWPDe)

    图  5  COLOR250数据集彩色图像灰度化转化结果

    ((a)输入彩色图像; (b) Gooch算法; (c) rgb2gray算法; (d) CP算法; (e) RTCP算法; (f) MWPDe算法)

    Fig.  5  Color-to-gray conversion results of COLOR250 dataset

    ((a) Input images; (b) Gooch; (c) rgb2gray; \\(d) CP; (e) RTCP; (f) MWPDe)

    图  6  选择性实验结果

    Fig.  6  Preference experiment result

    图  7  准确性实验结果

    Fig.  7  Accuracy experiment result

    图  8  噪声方差为0.03时对Cadik数据集进行抗噪性能分析

    ((a)输入彩色图像; (b) rgb2gray算法; (c) CP算法; (d) RTCP算法; (e) MWPDe算法)

    Fig.  8  Antinoise performance analysis of Cadik dataset with 0.03 noise variance

    ((a) Input images; (b) rgb2gray; (c) CP; (d) RTCP; (e) MWPDe)

    图  9  算法对低分辨率图像实验效果

    ((a)输入低分辨率彩色图像(下采样2个因子); (b) rgb2gray算法; (c) CP算法; (d) RTCP算法; (e) MWPDe算法)

    Fig.  9  Experiment results on the low resolution images

    ((a) Input low resolution color images (2 downsampling factor); (b) rgb2gray; (c) CP; (d) RTCP; (e) MWPDe)

    图  10  参数敏感度分析

    ((a)相似系数 ${{S}_{x, y}}(t)$ 随 $t$ 的变化趋势; (b)输入彩色图像; (c) $t=0.01$ 时WMPDe算法灰度化效果; (d) $t=0.00001$ 时WMPDe算法灰度化效果)

    Fig.  10  Sensitivity study of parameters

    ((a) Similar coefficients variation tendency ${{S}_{x, y}}(t)$ with $t$ ; (b) Input color images; (c) $t=0.01$ ; (d) $t=0.00001$ )

    表  1  Cadik数据集的CCPR值

    Table  1  CCPR value of Cadik dataset

    τ Gooch rgb2gray CP RTCP MWPDe
    1 0.96 0.93 0.97 0.96 0.96
    2 0.93 0.88 0.94 0.94 0.94
    3 0.91 0.84 0.92 0.92 0.92
    4 0.89 0.8 0.91 0.9 0.91
    5 0.87 0.77 0.89 0.89 0.9
    6 0.85 0.75 0.87 0.87 0.88
    7 0.83 0.73 0.86 0.86 0.87
    8 0.81 0.71 0.84 0.85 0.86
    9 0.79 0.69 0.83 0.83 0.85
    10 0.77 0.68 0.82 0.82 0.84
    11 0.75 0.66 0.8 0.8 0.84
    12 0.72 0.65 0.79 0.79 0.83
    13 0.7 0.63 0.77 0.77 0.82
    14 0.69 0.62 0.76 0.76 0.81
    15 0.67 0.61 0.75 0.75 0.81
    下载: 导出CSV

    表  2  CSDD数据集的CCPR值

    Table  2  CCPR value of CSDD dataset

    τ Gooch rgb2gray CP RTCP MWPDe
    1 0.97 0.96 0.96 0.96 0.95
    2 0.93 0.94 0.93 0.93 0.92
    3 0.91 0.91 0.91 0.91 0.91
    4 0.89 0.89 0.89 0.9 0.89
    5 0.87 0.87 0.87 0.88 0.88
    6 0.85 0.85 0.85 0.87 0.86
    7 0.83 0.83 0.83 0.85 0.85
    8 0.81 0.81 0.81 0.84 0.84
    9 0.79 0.8 0.79 0.82 0.83
    10 0.77 0.78 0.77 0.81 0.82
    11 0.76 0.76 0.75 0.8 0.81
    12 0.74 0.75 0.74 0.78 0.8
    13 0.73 0.73 0.72 0.77 0.79
    14 0.71 0.71 0.7 0.76 0.78
    15 0.7 0.68 0.7 0.75 0.77
    下载: 导出CSV

    表  3  COLOR250数据集的CCPR值

    Table  3  CCPR value of COLOR250

    τ Gooch rgb2gray CP RTCP MWPDe
    1 0.95 0.96 0.96 0.96 0.95
    2 0.92 0.93 0.93 0.93 0.92
    3 0.89 0.9 0.9 0.9 0.9
    4 0.86 0.87 0.88 0.88 0.88
    5 0.83 0.85 0.86 0.86 0.86
    6 0.81 0.83 0.84 0.84 0.85
    7 0.79 0.81 0.82 0.83 0.83
    8 0.76 0.79 0.8 0.81 0.82
    9 0.74 0.77 0.79 0.8 0.81
    10 0.72 0.75 0.77 0.78 0.8
    11 0.7 0.74 0.76 0.77 0.79
    12 0.68 0.72 0.74 0.76 0.78
    13 0.67 0.7 0.73 0.75 0.77
    14 0.65 0.69 0.72 0.74 0.76
    15 0.63 0.67 0.7 0.72 0.75
    下载: 导出CSV

    表  4  噪声情况下Cadik数据集的CCPR值

    Table  4  CCPR value of Cadik dataset with noise added

    τ rgb2gray CP RTCP MWPDe
    1 0.92 0.95 0.94 0.96
    2 0.86 0.91 0.91 0.92
    3 0.82 0.87 0.88 0.89
    4 0.78 0.84 0.85 0.86
    5 0.74 0.82 0.84 0.84
    6 0.72 0.79 0.82 0.82
    7 0.7 0.77 0.81 0.81
    8 0.68 0.75 0.8 0.79
    9 0.66 0.73 0.78 0.78
    10 0.65 0.71 0.77 0.77
    11 0.63 0.7 0.76 0.76
    12 0.62 0.68 0.75 0.75
    13 0.61 0.66 0.74 0.74
    14 0.6 0.65 0.73 0.74
    15 0.58 0.64 0.72 0.73
    下载: 导出CSV

    表  5  输入390×293彩色图像时不同算法的运行时间

    Table  5  Runtime of different algorithms with the input 390×293 color image

    Methods Gooch rgb2gray CP RTCP MWPDe
    Runtime (s) 3.00E+04 0.0048 0.5233 0.0728 0.0343
    下载: 导出CSV
  • [1] Bala R, Eschbach R. Spatial color-to-grayscale transform preserving chrominance edge information. In: Proceedings of the 12th Color Imaging Conference: Color Science and Engineering Systems, Technologies, and Applications. Scottsdale, USA: The Society for Imaging Science and Technology, 2004. 82-86
    [2] Neumann L, Čadík M, Nemcsics A. An efficient perception-based adaptive color to gray transformation. In: Proceedings of the 3rd Eurographics conference on Computational Aesthetics in Graphics, Visualization and Imaging. Switzerland: Eurographics Association, 2007. 73-80
    [3] Smith K, Landes P-E, Thollot J, Myszkowski K. Apparent greyscale: a simple and fast conversion to perceptually accurate images and video. Computer Graphics Forum, 2008, 27(2): 193-200 doi: 10.1111/j.1467-8659.2008.01116.x
    [4] Gooch A A, Olsen S C, Tumblin J, Gooch B. Color2gray: salience-preserving color removal. ACM Transactions on Graphics, 2005, 24(3): 634-639 doi: 10.1145/1073204
    [5] Rasche K, Geist R, Westall J. Detail preserving reproduction of color images for monochromats and dichromats. IEEE Computer Graphics and Applications, 2005, 25(3): 22-30 doi: 10.1109/MCG.2005.54
    [6] Kim Y, Jang C, Demouth J, Lee S. Robust color-to-gray via nonlinear global mapping. ACM Transactions on Graphics, 2009, 28(5): Article No.161
    [7] Du H, He S F, Sheng B, Ma L Z, Lau R W H. Saliency-guided color-to-gray conversion using region-based optimization. IEEE Transactions on Image Processing, 2015, 24(1): 434-443 doi: 10.1109/TIP.2014.2380172
    [8] Song M L, Tao D P, Chen C, Li X L, Chen C W. Color to gray: visual cue preservation. IEEE Transactions on Pattern Analysis and Machine Intelligence, 2010, 32(9): 1537-1552 doi: 10.1109/TPAMI.2009.74
    [9] Song M L, Tao D P, Chen C, Bu J J, Yang Y Z. Color-to-gray based on chance of happening preservation. Neurocomputing, 2013, 119: 222-231 doi: 10.1016/j.neucom.2013.03.037
    [10] 宋明黎, 王慧琼, 孙陈纯.基于高斯混合模型的色彩转换.计算机辅助设计与图形学学报, 2008, 20(11): 1471-1476 http://www.cnki.com.cn/Article/CJFDTOTAL-JSJF200811014.htm

    Song Ming-Li, Wang Hui-Qiong, Sun Chen-Chun. Gaussian mixture model based approach on color transfer. Journal of Computer-Aided Design and Computer Graphics, 2008, 20(11): 1471-1476 http://www.cnki.com.cn/Article/CJFDTOTAL-JSJF200811014.htm
    [11] 王辉, 陈小雕, 王毅刚.改进的非线性全局映射灰度化方法.计算机辅助设计与图形学学报, 2013, 25(10): 1476-1479 http://www.cnki.com.cn/Article/CJFDTOTAL-JSJF201310007.htm

    Wang Hui, Chen Xiao-Diao, Wang Yi-Gang. The improved color-to-gray via nonlinear global mapping. Journal of Computer-Aided Design and Computer Graphics, 2013, 25(10): 1476-1479 http://www.cnki.com.cn/Article/CJFDTOTAL-JSJF201310007.htm
    [12] 朱薇, 刘利刚. IHLS颜色空间的灰度化全局非线性三向映射.计算机辅助设计与图形学学报, 2013, 25(2): 154-159 http://www.cnki.com.cn/Article/CJFDTOTAL-JSJF201302003.htm

    Zhu Wei, Liu Li-Gang. Decolorize nonlinear tri-mapping in IHLS color space. Journal of Computer-Aided Design and Computer Graphics, 2013, 25(2): 154-159 http://www.cnki.com.cn/Article/CJFDTOTAL-JSJF201302003.htm
    [13] Lu C W, Xu L, Jia J Y. Contrast preserving decolorization. In: Proceedings of the 2012 IEEE International Conference on Computational Photography. Seattle, WA, USA: IEEE, 2012. 1-7
    [14] Lu C W, Xu L, Jia J Y. Real-time contrast preserving decolorization. In: Proceedings of the SIGGRAPH Asia 2012 Technical Briefs. Singapore, Singapore: ACM, 2012. Article No.34
    [15] Lu C W, Xu L, Jia J Y. Contrast preserving decolorization with perception-based quality metrics. International Journal of Computer Vision, 2014, 110(2): 222-239 doi: 10.1007/s11263-014-0732-6
    [16] Song Y B, Bao L C, Yang Q X. Real-time video decolorization using bilateral filtering. In: Proceedings of the 2014 IEEE Winter Conference on Applications of Computer Vision. Steamboat Springs, CO, USA: IEEE, 2014. 159-166
    [17] Liu Q G, Liu X P, Xie W S, Wang Y H, Liang D. GcsDecolor: gradient correlation similarity for efficient contrast preserving decolorization. IEEE Transactions on Image Processing, 2015, 24(9): 2889-2904 doi: 10.1109/TIP.2015.2423615
    [18] Lu G F, Lin Z, Jin Z. Face recognition using regularised generalised discriminant locality preserving projections. IET Computer Vision, 2011, 5(2): 107-116 doi: 10.1049/iet-cvi.2009.0138
    [19] Vasuhi S, Vaidehi V. Identification of human faces using orthogonal locality preserving projections. In: Proceedings of the 2009 International Conference on Signal Processing Systems. Singapore, Singapore: IEEE, 2009. 718-722
    [20] Gu X H, Yang L P, Peng J. Face recognition with discriminant locality preserving projections in complete kernel space. In: Proceedings of the 2011 International Conference on Computer Science and Network Technology. Harbin, China: IEEE, 2011. 1166-1169
    [21] Yi J, Ruan Q Q. Gabor-based improved locality preserving projections for face recognition. In: Proceedings of the 2007 IEEE International Conference on Image Processing. San Antonio, TX, USA: IEEE, 2007. I-153-I-156
    [22] Lin C, Sun X R, Zhou K J, Xia F. Supervised gabor-based kernel locality preserving projections for face recognition. In: Proceedings of the 2011 International Conference on and 4th International Conference on Cyber, Physical and Social Computing Internet of Things (iThings/CPSCom). Dalian, China: IEEE, 2011. 252-257
    [23] Wang H X. A new feature extraction method for image recognition using structural two-dimensional locality preserving projections. In: Proceedings of the 16th IEEE International Conference on Image Processing. Cairo, Egypt: IEEE, 2009. 2037-2040
    [24] Yi J, Ruan Q Q. An image matrix compression based supervised locality preserving projections for face recognition. In: Proceedings of the 2007 International Symposium on Intelligent Signal Processing and Communications Systems. Xiamen, China: IEEE, 2007. 738-741
    [25] He X F, Yan S C, Hu Y X, Niyogi P, Zhang H J. Face recognition using Laplacianfaces. IEEE Transactions on Pattern Analysis and Machine Intelligence, 2005, 27(3): 328-340 doi: 10.1109/TPAMI.2005.55
    [26] Čadík M. Perceptual evaluation of color-to-grayscale image conversions. Computer Graphics Forum, 2008, 27(7): 1745-1754 doi: 10.1111/cgf.2008.27.issue-7
  • 期刊类型引用(20)

    1. 叶荣,何云,高泉,章广传,邵郭奇,李彤. 基于OpenCV和模糊数学的茶叶病害分析方法研究. 北方园艺. 2024(04): 145-153 . 百度学术
    2. 雷永伟,雷蕾,宋冬利,林含潇. 复杂背景下连续帧的车轮廓形光条提取方法. 机车电传动. 2024(05): 124-131 . 百度学术
    3. 张俊升,王洪磊,李佳城. 基于双目结构光视觉的煤流量测量研究. 工矿自动化. 2023(07): 19-26 . 百度学术
    4. 廖一帆,李子豪,伍春花,汪国有,刘且根. 基于辅助变量增强的可逆彩色图像灰度化. 电子与信息学报. 2023(12): 4448-4457 . 百度学术
    5. 李秋菊,徐海平,苟春. 基于CIELab和HIS空间的声呐图像伪彩色增强方法及质量评价. 水下无人系统学报. 2022(02): 204-208+230 . 百度学术
    6. 刘亚,张团善,王恩芝. 基于高效边缘检测的彩色图像灰度化算法. 轻工机械. 2022(06): 52-58 . 百度学术
    7. 顾梅花,王苗苗,李立瑶,冯婧. 彩色图像多尺度融合灰度化算法. 计算机工程与应用. 2021(04): 209-215 . 百度学术
    8. 程鑫,同军超,施恢胜,张帆,吴臻冕,张卫超. 弱光条件下非规则曲面上浅刻蚀字符的识别研究. 数字制造科学. 2021(02): 81-86 . 百度学术
    9. 赵宇晓,李剑飞,朱志亮,段润涛,黄晓生. 基于四元数色度距离的对比度保持灰度化方法. 南京理工大学学报. 2021(04): 472-480 . 百度学术
    10. 张俊康,金正猛,冯灿. 基于原始对偶方法的图像去色算法. 南京邮电大学学报(自然科学版). 2021(04): 51-58 . 百度学术
    11. 陈扬,李旦,张建秋. 互补色小波颜色恒常性/白平衡方法. 自动化学报. 2020(07): 1378-1389 . 本站查看
    12. 迟臣鑫,陈伟强,朱鹏程,查剑锋,殷和健,马晨阳. 采煤塌陷积水区面积无人机采集方法. 金属矿山. 2020(08): 136-141 . 百度学术
    13. 王焱,吕猛,孟祥福,李宇浩. 具有光照鲁棒的图像匹配方法. 计算机应用. 2019(01): 262-266 . 百度学术
    14. 顾梅花,苏彬彬,王苗苗,王志磊. 彩色图像灰度化算法综述. 计算机应用研究. 2019(05): 1286-1292 . 百度学术
    15. 张超,李小平. 一种灰度化混合法在集装箱箱号识别中的运用. 计算机与现代化. 2019(05): 41-45 . 百度学术
    16. 张斌,曾碧,林伟. 光照不均环境下的手机外壳缺陷检测研究. 机电工程技术. 2019(09): 40-42+77 . 百度学术
    17. 曾绍华,唐文密,詹林庆,黄秀芬. 基于自适应密度峰值聚类的野外紫色土彩色图像分割. 农业工程学报. 2019(19): 200-208 . 百度学术
    18. 高直,吴宝福,王媛媛. 基于数学形态学的车牌自动识别系统. 计算机时代. 2018(07): 56-59+63 . 百度学术
    19. 吴建斌,周紫薇. 一种低成本平台下的水表码字识别方法研究. 软件导刊. 2018(07): 75-78 . 百度学术
    20. 丁李,曾水玲,熊涛. 基于BP神经网络的湘西民间方块苗文特征提取方法. 吉首大学学报(自然科学版). 2018(05): 14-18 . 百度学术

    其他类型引用(49)

  • 加载中
图(10) / 表(5)
计量
  • 文章访问数:  1996
  • HTML全文浏览量:  456
  • PDF下载量:  499
  • 被引次数: 69
出版历程
  • 收稿日期:  2016-01-22
  • 录用日期:  2016-04-28
  • 刊出日期:  2017-05-01

目录

/

返回文章
返回