CN108632583B - 数字图像处理方法及装置 - Google Patents

数字图像处理方法及装置 Download PDF

Info

Publication number
CN108632583B
CN108632583B CN201710177711.XA CN201710177711A CN108632583B CN 108632583 B CN108632583 B CN 108632583B CN 201710177711 A CN201710177711 A CN 201710177711A CN 108632583 B CN108632583 B CN 108632583B
Authority
CN
China
Prior art keywords
pixel
value
data
filtering
transient response
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201710177711.XA
Other languages
English (en)
Other versions
CN108632583A (zh
Inventor
胡毅
蔡进
王浩
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Spreadtrum Communications Shanghai Co Ltd
Original Assignee
Spreadtrum Communications Shanghai Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Spreadtrum Communications Shanghai Co Ltd filed Critical Spreadtrum Communications Shanghai Co Ltd
Priority to CN201710177711.XA priority Critical patent/CN108632583B/zh
Publication of CN108632583A publication Critical patent/CN108632583A/zh
Application granted granted Critical
Publication of CN108632583B publication Critical patent/CN108632583B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04N—PICTORIAL COMMUNICATION, e.g. TELEVISION
    • H04N23/00—Cameras or camera modules comprising electronic image sensors; Control thereof
    • H04N23/10—Cameras or camera modules comprising electronic image sensors; Control thereof for generating image signals from different wavelengths
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04N—PICTORIAL COMMUNICATION, e.g. TELEVISION
    • H04N9/00—Details of colour television systems
    • H04N9/64—Circuits for processing colour signals

Landscapes

  • Engineering & Computer Science (AREA)
  • Multimedia (AREA)
  • Signal Processing (AREA)
  • Image Processing (AREA)
  • Color Image Communication Systems (AREA)

Abstract

一种数字图像处理方法及装置。所述方法包括：对待处理的RGB数字图像进行瞬态响应滤波，得到经瞬态响应滤波后的RGB图像数据；对所述经瞬态响应滤波后的RGB图像数据进行伪彩色滤波，得到经伪彩色滤波后的RGB图像数据；对所述经瞬态响应滤波后的RGB图像数据及所述经伪彩色滤波后的RGB图像数据进行融合，输出所述待处理的RGB数字图像对应的去除伪彩色的RGB数字图像。应用上述方案，可以在不损失图像色度域的信息的情况下，有效去除数字图像中的伪彩色。

Description

数字图像处理方法及装置
技术领域
本发明实施例涉及图像处理技术领域，具体涉及一种数字图像处理方法及装置。
背景技术
为了将图像数字化，出现了由光学系统及图像传感器组成的数字图像处理系统。
在数字图像处理系统中，由于光学系统对于不同波长的光折射率不同以及图像传感器在高光下容易溢出等原因，经数字图像处理系统处理后的数字图像经常出现伪彩色，严重影响了数字图像的品质。
为了去除数字图像中的伪彩色，目前通常采用的做法是：将数字图像中所有接近于灰色的色彩抑制为灰色。由于人眼对于接近于灰色的伪彩色会较敏感，由此可以使得数字图像中用户肉眼可见的伪彩色被显著地减少。
然而，在上述伪彩色去除方法中，伪彩色的去除是以牺牲图像色度域的信息为代价的。若要不损失图像色度域的信息，数字图像中的伪彩色就会有所残留，数字图像的品质就会下降。
发明内容
本发明实施例要解决的问题是如何在不损失图像色度域的信息的情况下，有效去除数字图像中的伪彩色。
为解决上述问题，本发明实施例提供了一种数字图像处理方法，所述方法包括：对待处理的RGB数字图像进行瞬态响应滤波，得到经瞬态响应滤波后的RGB图像数据；对所述经瞬态响应滤波后的RGB图像数据进行伪彩色滤波，得到经伪彩色滤波后的RGB图像数据；对所述经瞬态响应滤波后的RGB图像数据及所述经伪彩色滤波后的RGB图像数据进行融合，输出所述待处理的RGB数字图像对应的去除伪彩色的RGB数字图像。
可选地，所述对待处理的RGB数字图像进行瞬态响应滤波，包括：采用如下步骤对所述待处理的RGB数字图像中第n个像素点进行瞬态响应滤波：以所述第n个像素点为中心的预设第一大小的图像块作为第一图像处理窗口，分别计算所述第一图像处理窗口中所述第n个像素点水平方向上以及垂直方向上R通道数据及B通道数据的最大值和最小值；分别计算所述第n个像素点水平方向上以及垂直方向上R通道数据及B通道数据瞬态响应滤波的最终结果，其中，任一方向上的任一通道数据均采用如下方式计算：基于所述第n个像素点相应方向上相应通道数据的最大值和最小值，对所述第n个像素点进行相应方向相应通道数据的瞬态响应滤波；计算所述第n个像素点相应方向相应通道数据瞬态响应滤波的上限值及下限值，并基于计算结果对第n个像素点相应方向相应通道数据瞬态响应滤波的结果进行限制，得到所述第n个像素点相应方向相应通道数据瞬态响应滤波的最终结果；分别对所述第n个像素点水平方向及垂直方向中相应通道数据瞬态响应滤波的最终结果进行融合，得到所述第n个像素点相应通道数据的瞬态响应滤波结果；其中，n为正整数。
可选地，所述计算所述第一图像处理窗口中所述第n个像素点水平方向上以及垂直方向上R通道数据及B通道数据的最大值和最小值，包括：计算相应方向上，位于所述第n个像素点两侧的像素点中相应通道数据的最大值及最小值；将所述第n个像素点第一侧相应通道数据的最大值与第二侧相应通道数据的最小值，和第二侧相应通道数据的最大值与第一侧相应通道数据的最小值中差异大的一组所对应的最大值及最小值，分别作为所述第n个像素点相应方向上相应通道数据的最大值及最小值。
可选地，所述基于所述第n个像素点相应方向上相应通道数据的最大值和最小值，对所述第n个像素点进行相应方向相应通道数据的瞬态响应滤波，包括：使用预设的瞬态响应滤波系数，对所述相应方向上相应通道数据的最大值和最小值及所述第n个像素点进行卷积运算。
可选地，所述使用预设的瞬态响应滤波系数，对所述相应方向上相应通道数据的最大值和最小值及所述第n个像素点进行卷积运算，包括：当XIN(i，j)>GIN(i，j)时，采用如下方法获得第n个像素点进行相应方向相应通道数据的瞬态响应滤波的结果：XPF(i，j)＝[p0，p1，p2]*[XMAX，XIN(i，j)，XMIN]T；当XIN(i，j)≤GIN(i，j)时，采用如下方法获得第n个像素点进行相应方向相应通道数据的瞬态响应滤波的结果：XPF(i，j)＝[p0，p1，p2]*[XMIN，XIN(i，j)，XMAX]T；其中，[p0，p1，p2]为所述瞬态响应滤波系数，X为R通道或B通道，XMAX为相应方向上R通道数据的最大值或者为相应方向上B通道数据的最大值，(i，j)为第n个像素点的坐标，XIN(i，j)为第n个像素点R通道或B通道数据，GIN(i，j)为第n个像素点的G通道数据，XPF(i，j)为第n个像素点R通道或B通道数据的瞬态响应滤波结果。
可选地，[p0，p1，p2]＝[-0.25，1.375，-0.125]。
可选地，所述计算所述第n个像素点相应方向相应通道数据瞬态响应滤波的上限值及下限值，包括：根据所述第n个像素点相应方向相应通道数据与所述第n个像素点G通道数据的比较结果，确定所述第n个像素点相应方向相应通道数据瞬态响应滤波的上限值及下限值。
可选地，所述对所述第n个像素点同一方向中各通道瞬态响应滤波的最终结果进行融合，包括：对所述第n个像素点水平方向上相应通道数据瞬态响应滤波的最终结果，与垂直方向上相应通道数据瞬态响应滤波的最终结果进行比较，将二者中较小者作为所述第n个像素点相应通道数据的瞬态响应滤波结果。
可选地，所述对所述经瞬态响应滤波后的RGB图像数据进行伪彩色滤波，包括：采用如下步骤对所述经瞬态响应滤波后的RGB图像数据中第k个像素点进行水平方向及垂直方向伪彩色滤波，k为正整数：以所述第k个像素点为中心的预设第二大小的图像块作为第二图像处理窗口，对所述经瞬态响应滤波后的RGB图像数据进行色彩空间转换，得到所述第k个像素点相应方向上各个像素点R通道数据对应的色差信号KR的值，以及B通道数据对应的色差信号KB的值；基于所述第k个像素点相应方向上各个像素点的相应色差信号值，对所述第k个像素点相应方向上相应色差信号值进行伪彩色滤波；对所述第k个像素点相应方向上相应色差信号值的伪彩色滤波结果进行融合，得到所述第k个像素点相应色差信号值的伪彩色滤波结果。
可选地，所述基于所述第k个像素点相应方向上各个像素点的相应色差信号值，对所述第k个像素点相应方向上相应色差信号值进行伪彩色滤波，包括：基于所述第k个像素点相应方向上各个像素点相应色差信号值的符号，对所述第k个像素点相应方向上各个像素点相应色差信号值进行滤波处理；计算所述第k个像素点相应方向上各个像素点相应色差信号值的权重系数；基于所述第k个像素点相应方向上各个像素点相应色差信号值的滤波处理结果，以及所述各个像素点的权重系数，对所述第k个像素点相应方向上相应色差信号值进行相应方向伪彩色滤波。
可选地，所述计算所述第k个像素点相应方向上各个像素点的权重系数，包括：计算所述第k个像素点相应方向上各个像素点的R、B及G通道数据对应的梯度；基于第k个像素点相应方向上各个像素点的R、B及G通道数据对应的梯度，计算所述第k个像素点相应方向上各个像素点的权重系数。
可选地，所述计算所述第k个像素点相应方向上各个像素点的权重系数，包括：计算所述第k个像素点相应方向上各个像素点的R、B、G及Y通道数据对应的梯度；基于第k个像素点相应方向上各个像素点的R、B、G及Y通道数据对应的梯度，计算所述第k个像素点相应方向上各个像素点的权重系数。
可选地，所述基于所述第k个像素点相应方向上各个像素点相应色差信号值的滤波处理结果，以及所述各个像素点的权重系数，对所述各个像素点进行相应方向伪彩色滤波，包括：采用如下方法获得所述第k个像素点相应方向相应色差信号值的伪彩色滤波结果：其中，(p,q)为所述第k个像素点的坐标，KX(p,q)为第k个像素点相应方向色差信号KX值的伪彩色滤波结果，r∈[-L，L]且r为整数，L为所述第二图像处理窗口相应方向的半径，wk(p,q+r)为像素点(p,q+r)的权重系数，CLIP_KXTI(p,q+r)为像素点(p,q+r)X通道数据的滤波处理结果。
可选地，所述对所述经瞬态响应滤波后的RGB图像数据中第k个像素点进行水平方向及垂直方向伪彩色滤波，还包括：对所述第k个像素点相应方向上各个像素点相应色差信号值进行中值滤波，得到所述经瞬态响应滤波后的RGB图像数据中各个像素点相应色差信号值的中值滤波结果；将第k个像素点相应色差信号值的伪彩色滤波结果，与第k个像素点相应方向相应色差信号值的中值滤波结果进行融合，融合结果作为第k个像素点相应色差信号值伪彩色滤波的最终结果。
可选地，所述对所述经瞬态响应滤波后的RGB图像数据及所述经伪彩色滤波后的RGB图像数据进行融合，包括：计算所述待处理的RGB数字图像中各个像素点相应通道数据的局部对比度值；基于所述待处理的RGB数字图像中各个像素点相应通道数据的局部对比度值，计算相应像素点相应通道数据的融合参数；基于所述融合参数，对相应像素点相应通道数据对应的瞬态响应滤波结果及伪彩色滤波结果进行融合，作为所述相应像素点去除伪彩色后相应通道的数据值。
可选地，所述基于所述待处理的RGB数字图像中各个像素点的局部对比度值，计算相应像素点的融合参数，包括：采用如下方法计算像素点(h,k)的X通道数据的融合参数Xalpha(i，j)：
其中，Xcontrast(i,j)为像素点(i，j)的局部对比度值，Xmax为像素点(i，j)X通道数据的最大值，Xmin为像素点(i，j)X通道数据的最小值，sigma1及sigma2为预设的调节参数，X为R通道或B通道。
可选地，所述基于所述融合参数，对相应像素点相应通道数据对应的瞬态响应滤波结果及伪彩色滤波结果进行融合，包括：采用如下方法得到像素点(h,k)去除伪彩色后X通道数据的数据值Xfcsoutput(h,k)：
Xfcsoutput(i,j)＝(1-Xalpha(i,j))*XTI(i,j)+Xalpha(i,j)*Xfcs(i,j)
其中，XTI(i,j)为像素点(i，j)X通道的瞬态响应滤波结果，Xfcs(i,j)为像素点(i，j)X通道数据的伪彩色滤波结果，Xfcs(i,j)＝KX(i，j)+GIN(i，j)，GIN(i，j)为像素点(i，j)的G通道数据。
本发明实施例还提供了一种数字图像处理装置，所述装置包括：瞬态响应滤波器，适于对待处理的RGB数字图像进行瞬态响应滤波，得到经瞬态响应滤波后的RGB图像数据；伪彩色滤波器，适于对所述经瞬态响应滤波后的RGB图像数据进行伪彩色滤波，得到经伪彩色滤波后的RGB图像数据；融合处理器，适于对所述经瞬态响应滤波后的RGB图像数据及所述经伪彩色滤波后的RGB图像数据进行融合，输出所述待处理的RGB数字图像对应的去除伪彩色的RGB数字图像。
可选地，所述瞬态响应滤波器，适于对所述待处理的RGB数字图像中第n个像素点进行瞬态响应滤波，n为正整数，包括：第一计算单元，适于以所述第n个像素点为中心的预设第一大小的图像块作为第一图像处理窗口，分别计算所述第一图像处理窗口中所述第n个像素点水平方向上以及垂直方向上R通道数据及B通道数据的最大值和最小值；第二计算单元，适于分别计算所述第n个像素点水平方向上以及垂直方向上R通道数据及B通道数据瞬态响应滤波的最终结果；包括：瞬态响应滤波子单元、第一计算子单元及处理子单元，其中：所述瞬态响应滤波子单元，适于基于所述第n个像素点相应方向上相应通道数据的最大值和最小值，对所述第n个像素点进行相应方向相应通道数据的瞬态响应滤波；所述第一计算子单元，适于计算所述第n个像素点相应方向相应通道数据瞬态响应滤波的上限值及下限值；所述处理子单元，适于基于计算结果对第n个像素点相应方向相应通道数据瞬态响应滤波的结果进行限制，得到所述第n个像素点相应方向相应通道数据瞬态响应滤波的最终结果；第一融合单元，适于分别对所述第n个像素点水平方向及垂直方向中相应通道数据瞬态响应滤波的最终结果进行融合，得到所述第n个像素点相应通道数据的瞬态响应滤波结果。
可选地，所述第一计算单元，适于计算相应方向上，位于所述第n个像素点两侧的像素点中相应通道数据的最大值及最小值；将所述第n个像素点第一侧相应通道数据的最大值与第二侧相应通道数据的最小值，和第二侧相应通道数据的最大值与第一侧相应通道数据的最小值中差异大的一组所对应的最大值及最小值，分别作为所述第n个像素点相应方向上相应通道数据的最大值及最小值。
可选地，所述瞬态响应滤波子单元，适于使用预设的瞬态响应滤波系数，对所述相应方向上相应通道数据的最大值和最小值及所述第n个像素点进行卷积运算。
可选地，所述瞬态响应滤波子单元，适于当XIN(i，j)>GIN(i，j)时，采用如下方法获得第n个像素点进行相应方向相应通道数据的瞬态响应滤波的结果：XPF(i，j)＝[p0，p1，p2]*[XMAX，XIN(i，j)，XMIN]T；当XIN(i，j)≤GIN(i，j)时，采用如下方法获得第n个像素点进行相应方向相应通道数据的瞬态响应滤波的结果：XPF(i，j)＝[p0，p1，p2]*[XMIN，XIN(i，j)，XMAX]T；其中，[p0，p1，p2]为所述瞬态响应滤波系数，X为R通道或B通道，XMAX为相应方向上R通道数据的最大值或者为相应方向上B通道数据的最大值，(i，j)为第n个像素点的坐标，XIN(i，j)为第n个像素点R通道或B通道数据，GIN(i，j)为第n个像素点的G通道数据，XPF(i，j)为第n个像素点R通道或B通道数据的瞬态响应滤波结果。
可选地，[p0，p1，p2]＝[-0.25，1.375，-0.125]。
可选地，所述第一计算子单元，适于根据所述第n个像素点相应方向相应通道数据与所述第n个像素点G通道数据的比较结果，确定所述第n个像素点相应方向相应通道数据瞬态响应滤波的上限值及下限值。
可选地，所述第一融合单元，适于对所述第n个像素点水平方向上相应通道数据瞬态响应滤波的最终结果，与垂直方向上相应通道数据瞬态响应滤波的最终结果进行比较，将二者中较小者作为所述第n个像素点相应通道数据的瞬态响应滤波结果。
可选地，所述伪彩色滤波器，适于对所述经瞬态响应滤波后的RGB图像数据中第k个像素点进行水平方向及垂直方向伪彩色滤波，k为正整数，包括：色彩空间转换单元，适于以所述第k个像素点为中心的预设第二大小的图像块作为第二图像处理窗口，对所述经瞬态响应滤波后的RGB图像数据进行色彩空间转换，得到所述第k个像素点相应方向上各个像素点R通道数据对应的色差信号KR的值，以及B通道数据对应的色差信号KB的值；伪彩色滤波单元，适于基于所述第k个像素点相应方向上各个像素点的相应色差信号值，对所述第k个像素点相应方向上相应色差信号值进行伪彩色滤波；第二融合单元，适于对所述第k个像素点相应方向上相应色差信号值的伪彩色滤波结果进行融合，得到所述第k个像素点相应色差信号值的伪彩色滤波结果。
可选地，所述伪彩色滤波单元包括：第一滤波子单元，适于基于所述第k个像素点相应方向上各个像素点相应色差信号值的符号，对所述第k个像素点相应方向上各个像素点相应色差信号值进行滤波处理；第二计算子单元，适于计算所述第k个像素点相应方向上各个像素点相应色差信号值的权重系数；第二滤波子单元，适于基于所述第k个像素点相应方向上各个像素点相应色差信号值的滤波处理结果，以及所述各个像素点的权重系数，对所述第k个像素点相应方向上相应色差信号值进行相应方向伪彩色滤波。
可选地，所述第二计算子单元，适于计算所述第k个像素点相应方向上各个像素点的R、B及G通道数据对应的梯度；基于第k个像素点相应方向上各个像素点的R、B及G通道数据对应的梯度，计算所述第k个像素点相应方向上各个像素点的权重系数。
可选地，所述第二计算子单元，适于计算所述第k个像素点相应方向上各个像素点的R、B、G及Y通道数据对应的梯度；基于第k个像素点相应方向上各个像素点的R、B、G及Y通道数据对应的梯度，计算所述第k个像素点相应方向上各个像素点的权重系数。
可选地，所述第二滤波子单元，适于采用如下方法获得所述第k个像素点相应方向相应色差信号值的伪彩色滤波结果：
其中，(p,q)为所述第k个像素点的坐标，KX(p,q)为第k个像素点相应方向色差信号KX值的伪彩色滤波结果，r∈[-L，L]且r为整数，L为所述第二图像处理窗口相应方向的半径，wk(p,q+r)为像素点(p,q+r)的权重系数，CLIP_KXTI(p,q+r)为像素点(p,q+r)X通道数据的滤波处理结果。
可选地，所述伪彩色滤波器，还包括：第三滤波单元，对所述第k个像素点相应方向上各个像素点相应色差信号值进行中值滤波，得到所述经瞬态响应滤波后的RGB图像数据中各个像素点相应色差信号值的中值滤波结果；第三融合单元，适于将第k个像素点相应色差信号值的伪彩色滤波结果，与第k个像素点相应方向相应色差信号值的中值滤波结果进行融合，融合结果作为第k个像素点相应色差信号值伪彩色滤波的最终结果。
可选地，所述融合处理器，包括：第三计算单元，适于计算所述待处理的RGB数字图像中各个像素点相应通道数据的局部对比度值；第四计算单元，适于基于所述待处理的RGB数字图像中各个像素点相应通道数据的局部对比度值，计算相应像素点相应通道数据的融合参数；第四融合单元，适于基于所述融合参数，对相应像素点相应通道数据对应的瞬态响应滤波结果及伪彩色滤波结果进行融合，作为所述相应像素点去除伪彩色后相应通道的数据值。
可选地，所述第四计算单元，适于采用如下方法计算像素点(h,k)的X通道数据的融合参数Xalpha(i，j)：
其中，Xcontrast(i,j)为像素点(i，j)的局部对比度值，Xmax为像素点(i，j)X通道数据的最大值，Xmin为像素点(i，j)X通道数据的最小值，sigma1及sigma2为预设的调节参数，X为R通道或B通道。
可选地，所述第四融合单元，适于采用如下方法得到像素点(h,k)去除伪彩色后X通道数据的数据值Xfcsoutput(h,k)：
Xfcsoutput(i,j)＝(1-Xalpha(i,j))*XTI(i,j)+Xalpha(i,j)*Xfcs(i,j)
其中，XTI(i,j)为像素点(i，j)X通道的瞬态响应滤波结果，Xfcs(i,j)为像素点(i，j)X通道数据的伪彩色滤波结果，Xfcs(i,j)＝KX(i，j)+GIN(i，j)，GIN(i，j)为像素点(i，j)的G通道数据。
相对于现有技术，本发明实施例的优点在于：
采用上述方案，对待处理的RGB数字图像整体去除伪彩色，相对于仅对数字图像中接近于灰色的色彩区域去除伪彩色，可以避免数字图像中低饱和度的色彩信息的损失，也就可以避免造成图像色度域的信息的损失，故可以提高数字图像的品质。并且，先对待处理的RGB数字图像进行瞬态响应滤波后，再进行伪彩色滤波，可以压缩伪彩色在局部图像上的范围，在需要去除同样宽度的伪彩色的前提下，有效减小计算量。
进一步地，使用像素点不同通道数据的梯度信息，计算该像素点的相应方向上各个像素点的权重系数，进而基于计算得到的权重系数进行伪彩色滤波，可以在去除数字图像中伪彩色的同时，更好地保持数字图像中带颜色纹理区域的色彩。
进一步地，基于待处理的RGB数字图像中各个像素点的局部对比度值，计算像素点的融合参数，可以使得局部对比度强的区域(也就是伪彩色更容易产生的区域)更多使用伪彩色滤波的结果，而其他区域则倾向于使用瞬态响应滤波的结果，从而可以在更好地去除伪彩色的同时，保护数字图像的色彩平衡。
附图说明
图1是本发明实施例中一种数字图像的处理方法的流程图；
图2是本发明实施例中一种对待处理的RGB数字图像进行瞬态响应滤波的方法的流程图；
图3是本发明实施例中一种第一图像处理窗口的示意图；
图4是本发明实施例中一种伪彩色滤波的方法的流程图；
图5是本发明实施例中一种第二图像处理窗口的示意图；
图6是本发明实施例中一种对第k个像素点相应方向上相应色差信号值进行伪彩色滤波的方法的流程图；
图7是本发明实施例中一种图像融合的方法的流程图；
图8是本发明实施例中另一种伪彩色滤波的方法的流程图；
图9是本发明实施例中另一种数字图像的处理方法的流程图；
图10是本发明实施例中一种数字图像处理装置的结构示意图；
图11是本发明实施例中一种瞬态响应滤波器的结构示意图；
图12是本发明实施例中一种伪彩色滤波器的结构示意图；
图13是本发明实施例中一种融合处理器的结构示意图。
具体实施方式
现有的伪彩色去除方法中，仅将数字图像中接近于灰色的色彩抑制为灰色。若要不损失图像色度域的信息，数字图像中的伪彩色就会有所残留，数字图像的品质就会下降。
针对上述问题，本发明的实施例提供了一种数字图像处理方法，应用所述方法对待处理的RGB数字图像整体去除伪彩色，相对于仅对数字图像中接近于灰色的色彩区域去除伪彩色，可以避免数字图像中低饱和度的色彩信息的损失，也就可以避免造成图像色度域的信息的损失，故可以提高数字图像的品质。并且，先对待处理的RGB数字图像进行瞬态响应滤波后，再进行伪彩色滤波，可以压缩伪彩色在局部图像上的范围，在需要去除同样宽度的伪彩色的前提下，有效减小计算量。
为使本发明的上述目的、特征和优点能够更为明显易懂，下面结合附图对本发明的具体实施例作详细地说明。
参照图1，本发明实施例提供了一种数字图像的处理方法，所述方法可以包括如下步骤：
步骤11，对待处理的RGB数字图像进行瞬态响应滤波，得到经瞬态响应滤波后的RGB图像数据。
对待处理的RGB数字图像整体进行瞬态响应滤波，可以压缩伪彩色在局部图像上的范围，在需要去除同样宽度的伪彩色的前提下，有效减小计算量。
步骤12，对所述经瞬态响应滤波后的RGB图像数据进行伪彩色滤波，得到经伪彩色滤波后的RGB图像数据。
对所述经瞬态响应滤波后的RGB图像数据进行伪彩色滤波，可以在保留图像中低饱和度色彩信息的前提下，有效去除数字图像中的伪彩色，提高数字图像品质。
步骤13，对所述经瞬态响应滤波后的RGB图像数据及所述经伪彩色滤波后的RGB图像数据进行融合，输出所述待处理的RGB数字图像对应的去除伪彩色的RGB数字图像。
对经瞬态响应滤波后的RGB图像数据及经伪彩色滤波后的RGB图像数据进行融合，可以在去除伪彩色的同时，保护数字图像中的色彩均衡。
参照图2，本发明实施例还提供了一种对待处理的RGB数字图像进行瞬态响应滤波的方法，利用所述方法可以对所述待处理的RGB数字图像中第n个像素点进行R通道数据及B通道数据的瞬态响应滤波。
具体地，所述方法可以包括如下步骤：
步骤21，以所述第n个像素点为中心的预设第一大小的图像块作为第一图像处理窗口，分别计算所述第一图像处理窗口中所述第n个像素点水平方向上以及垂直方向上R通道数据及B通道数据的最大值和最小值，n为正整数。
在具体实施中，所述第一图像处理窗口的大小可以根据RGB数字图像的大小进行选择。例如，参照图3，将像素点CUV作为RGB数字图像的第n个像素点，所述第一图像处理窗口的大小可以为以像素点CUV为中心的7*13的图像块。
在具体实施中，当RGB数字图像中无法构建以某一像素点为中心的第一图像处理窗口时，可以对所述RGB数字图像进行图像边界扩展或图像边界镜像处理，进而获得以该像素点为中心的第一图像处理窗口。比如，当该像素点为RGB数字图像的第一个像素点或者最后一个像素点时，可以对所述RGB数字图像进行图像边界扩展或图像边界镜像处理，进而获得第一个像素点或者最后一个像素点为中心的第一图像处理窗口。
下面以图3中示出的第一图像处理窗口为例，对如何计算像素点CUV水平方向上以及垂直方向上R通道数据及B通道数据的最大值和最小值进行说明：
参照图3，像素点CUV水平方向上的像素点共12个，分别为像素点L5、L4、L3、L2、L1、L0、R0、R1、R2、R3、R4及R5。像素点CUV垂直方向上的像素点共6个，分别为像素点T2、T1、T0、B0、B1及B2。将像素点CUV及其水平方向的12个像素点的像素值，作为对像素点CUV进行水平方向瞬态响应滤波的输入数据。将像素点CUV及其垂直方向上的6个像素点的像素值，作为对像素点CUV进行垂直方向瞬态响应滤波的输入数据。
在本发明的一实施例中，可以先计算相应方向上，位于所述第n个像素点两侧的像素点中相应通道数据的最大值及最小值，进而将所述第n个像素点第一侧相应通道数据的最大值与第二侧相应通道数据的最小值，和第二侧相应通道数据的最大值与第一侧相应通道数据的最小值中差异大的一组所对应的最大值及最小值，分别作为所述第n个像素点相应方向上相应通道数据的最大值及最小值。
以图3中的像素点CUV为所述第n个像素点为例，计算像素点CUV水平方向上R通道数据的最大值R1MAX及最小值R1MIN时，可以先计算像素点CUV左侧的各个像素点中R通道数据的最大值REMAX及最小值REMIN，再计算像素点CUV右侧的各个像素点中R通道数据的最大值RWMAX及最小值RWMIN。
若(RWMAX-REMIN)＞(REMAX-RWMIN)，则R1MAX＝RWMAX，R1MIN＝REMIN。
若(REMAX-RWMIN)＞(RWMAX-REMIN)，则R1MAX＝REMAX，R1MIN＝RWMIN。
相似地，可以计算得到像素点CUV水平方向上B通道数据的最大值B1MAX及最小值B1MIN。
计算像素点CUV垂直方向上R通道数据的最大值R2MAX及最小值R2MIN时，可以先计算像素点CUV上方的各个像素点中R通道数据的最大值RUMAX及最小值RUMIN，再计算像素点CUV下方的各个像素点中R通道数据的最大值RDMAX及最小值RDMIN。
若(RUMAX-RDMIN)＞(RDMAX-RUMIN)，则R2MAX＝RUMAX，R2MIN＝RDMIN。
若(RDMAX-RUMIN)＞(RUMAX-RDMIN)，则R2MAX＝RDMAX，R2MIN＝RUMIN。
相似地，可以计算得到像素点CUV垂直方向上B通道数据的最大值B2MAX及最小值B2MIN。
步骤22，分别计算所述第n个像素点水平方向上以及垂直方向上R通道数据及B通道数据瞬态响应滤波的最终结果。
在具体实施中，任一方向上的任一通道数据均采用如下方式计算：基于所述第n个像素点相应方向上相应通道数据的最大值和最小值，对所述第n个像素点进行相应方向相应通道数据的瞬态响应滤波；计算所述第n个像素点相应方向相应通道数据瞬态响应滤波的上限值及下限值，并基于计算结果对第n个像素点相应方向相应通道数据瞬态响应滤波的结果进行限制，得到所述第n个像素点相应方向相应通道数据瞬态响应滤波的最终结果。
在本发明的一实施例中，可以使用预设的瞬态响应滤波系数，对所述相应方向上相应通道数据的最大值和最小值及所述第n个像素点进行卷积运算，得到所述第n个像素点相应方向相应通道数据的瞬态响应滤波结果。
在具体实施中，使用预设的瞬态响应滤波系数，可以采用多种方式，对第n个像素点相应方向相应通道数据的瞬态响应滤波，得到第n个像素点水平方向相应通道数据的瞬态响应滤波的结果X1PF(i，j)，以及垂直方向相应通道数据的瞬态响应滤波的结果X2PF(i，j)。
本发明的实施例中，为了便于描述，X表示水平或垂直方向相应通道，X1表示水平方向R通道或B通道，X2为垂直方向R通道或B通道。X可以为X1，也可以为X2。
在本发明的一实施例中，当XIN(i，j)>GIN(i，j)时，可以采用如下方法获得第n个像素点进行相应方向相应通道数据的瞬态响应滤波的结果：
XPF(i，j)＝[p0，p1，p2]*[XMAX，XIN(i，j)，XMIN]T (1)
当XIN(i，j)≤GIN(i，j)时，可以采用如下方法获得第n个像素点进行相应方向相应通道数据的瞬态响应滤波的结果：
XPF(i，j)＝[p0，p1，p2]*[XMIN，XIN(i，j)，XMAX]T (2)
其中，[p0，p1，p2]为所述瞬态响应滤波系数，X为R通道或B通道，XMAX为相应方向上R通道数据的最大值或者为相应方向上B通道数据的最大值，(i，j)为第n个像素点的坐标，XIN(i，j)为第n个像素点R通道或B通道数据，GIN(i，j)为第n个像素点的G通道数据，XPF(i，j)为第n个像素点R通道或B通道数据的瞬态响应滤波结果。
以图3中的像素点CUV为所述第n个像素点为例，像素点CUV的坐标为(2，6)。像素点CUV的G通道数据为GIN(2，6)，R通道数据为RIN(2，6)，B通道数据为BIN(2，6)。像素点CUV水平方向上R通道数据的最大值R1MAX及最小值R1MIN。像素点CUV垂直方向上R通道数据的最大值R2MAX及最小值R2MIN。
假设RIN(2，6)＞GIN(2，6)，则利用公式(1)可以获得像素点CUV水平及垂直方向上R通道数据的瞬态响应滤波的结果R1PF(2，6)及R2PF(2，6)，以及获得像素点CUV水平及垂直方向上B通道数据的瞬态响应滤波的结果B1PF(2，6)及B2PF(2，6)，其中：
R1PF(2，6)＝[p0，p1，p2]*[R1MAX，RIN(2，6)，R1MIN]T；
R2PF(2，6)＝[p0，p1，p2]*[R2MAX，RIN(2，6)，R2MIN]T；
B1PF(2，6)＝[p0，p1，p2]*[B1MAX，RIN(2，6)，B2MIN]T；
B2PF(2，6)＝[p0，p1，p2]*[B2MAX，RIN(2，6)，B2MIN]T。
在具体实施中，所述瞬态响应滤波系数[p0，p1，p2]可以根据实际情况进行设置，在本发明的一实施例中，为了获得较好的瞬态响应滤波效果，可以设置[p0，p1，p2]＝[-0.25，1.375，-0.125]。
在具体实施中，可以采用多种方法计算所述第n个像素点相应方向相应通道数据瞬态响应滤波的上限值XTIMAX(i，j)及下限值XTIMIN(i，j)。
在本发明的一实施例中，可以根据所述第n个像素点相应方向相应通道数据与所述第n个像素点G通道数据的比较结果，确定所述第n个像素点相应方向相应通道数据瞬态响应滤波结果的上限值XTIMAX(i，j)及下限值XTIMIN(i，j)。其中，第n个像素点水平方向相应通道数据瞬态响应滤波结果的上限值为X1TIMAX(i，j)，下限值为X1TIMIN(i，j)。第n个像素点垂直方向相应通道数据瞬态响应滤波结果的上限值为X2TIMAX(i，j)，下限值为X2TIMIN(i，j)。
在本发明的一实施例中，当XIN(i，j)>GIN(i，j)时，XTIMAX(i，j)＝XIN(i，j)，XTIMIN(i，j)为XMIN及GIN(i，j)中的较小值，即XTIMIN(i，j)＝min(XMAX，GIN(i，j))。
当XIN(i，j)≤GIN(i，j)时，XTIMAX(i，j)为XMIN及GIN(i，j)中的较大值，即XTIMAX(i，j)＝max(XMIN，GIN(i，j))；XTIMAX(i，j)＝XMIN。
利用相应的上限值及下限值将第n个像素点水平方向相应通道数据瞬态响应滤波结果进行限制，得到第n个像素点水平方向相应通道数据瞬态响应滤波最终结果为X1TI(i，j):
当XPF(i，j)>XTIMAX(i，j)时，X1TI(i，j)＝XTIMAX(i，j)；
当XPF(i，j)≤XTIMIN(i，j)时，X1TI(i，j)＝XTIMIN(i，j)；
当XTIMAX≥XPF(i，j)>XTIMIN(i，j)时，X1TI(i，j)＝XPF(i，j)。
相似地，利用相应的上限值及下限值将第n个像素点垂直方向相应通道数据瞬态响应滤波结果进行限制，可以得到第n个像素点垂直方向相应通道数据瞬态响应滤波最终结果为X2TI(i，j)。
以图3中的像素点CUV(2，6)为所述第n个像素点为例，像素点CUV水平方向R通道数据瞬态响应滤波结果的上限值为R1TIMAX(2，6)，下限值为R1TIMIN(2，6)。像素点CUV水平方向B通道数据瞬态响应滤波结果的上限值为B1TIMAX(2，6)，下限值为B1TIMIN(2，6)。像素点CUV垂直方向R通道数据瞬态响应滤波结果的上限值为R2TIMAX(2，6)，下限值为R2TIMIN(2，6)。像素点CUV垂直方向B通道数据瞬态响应滤波结果的上限值为B2TIMAX(2，6)，下限值为B2TIMIN(2，6)。
假设RIN(2，6)＞GIN(2，6)，则可以得到：
R1TIMAX(2，6)＝R1IN(2，6)；
R1MIN(2，6)＝min(R1MAX，GIN(2，6))；
B1TIMAX(2，6)＝B1IN(2，6)；
B1TIMIN(2，6)min(B1MAX，GIN(2，6))；
R2TIMAX(2，6)＝R2IN(2，6)；
R2MIN(2，6)＝min(R2MAX，GIN(2，6))；
B2TIMAX(2，6)＝B2IN(2，6)；
B2TIMIN(2，6)min(B2MAX，GIN(2，6))。
假设R1PF(2，6)>R1TIMAX，像素点CUV水平方向R通道数据瞬态响应滤波R1TI(2，6)＝R1TIMAX(2，6)＝R1IN(2，6)；
假设R2PF(2，6)≤R1TIMAX(2，6)，像素点CUV垂直方向R通道数据瞬态响应滤波R2TI(2，6)＝R2MIN(2，6)＝min(R2MAX，GIN(2，6))；
假设B1PF(2，6)>R1TIMAX，像素点CUV水平方向B通道数据瞬态响应滤波B1TI(2，6)＝B1TIMAX(2，6)＝B1IN(2，6)；
假设B2TIMAX≥B2PF(2，6)>B2TIMIN(2，6)，像素点CUV垂直方向B通道数据瞬态响应滤波B2TI(2，6)＝B2PF(2，6)。
步骤23，分别对所述第n个像素点水平方向及垂直方向中相应通道数据瞬态响应滤波的最终结果进行融合，得到所述第n个像素点相应通道数据的瞬态响应滤波结果。
在具体实施中，可以采用多种方法将第n个像素点水平方向及垂直方向中相应通道数据瞬态响应滤波的最终结果进行融合，具体不作限制。
在本发明的一实施例中，可以先对所述第n个像素点水平方向上相应通道数据瞬态响应滤波的最终结果，与垂直方向上相应通道数据瞬态响应滤波的最终结果进行比较，进而将二者中较小者作为所述第n个像素点相应通道数据的瞬态响应滤波结果XTI(i，j)＝min(X1TI(i，j)，X2TI(i，j))。
具体地，令第n个像素点水平方向R通道数据瞬态响应滤波最终结果为R1TI(i，j)，垂直方向R通道数据瞬态响应滤波最终结果为R2TI(i，j)，则第n个像素点R通道数据瞬态响应滤波最终结果RTI(i，j)＝min(R1TI(i，j)，R2TI(i，j))。
令第n个像素点水平方向B通道数据瞬态响应滤波最终结果为B1TI(i，j)，垂直方向B通道数据瞬态响应滤波最终结果为B2TI(i，j)，则第n个像素点B通道数据瞬态响应滤波最终结果BTI(i，j)＝min(B1TI(i，j)，B2TI(i，j))。
参照图4，本发明实施例还提供了一种伪彩色滤波的方法，利用所述方法可以对经瞬态响应滤波后的RGB图像数据中第k个像素点进行水平方向及垂直方向的伪彩色滤波，进而得到经瞬态响应滤波后的RGB图像数据中各个像素点的伪彩色滤波结果，其中k为正整数。
具体地，所述方法可以包括如下步骤：
步骤41，以所述第k个像素点为中心的预设第二大小的图像块作为第二图像处理窗口，对所述经瞬态响应滤波后的RGB图像数据进行色彩空间转换，得到所述第k个像素点相应方向上各个像素点R通道数据对应的色差信号KR的值，以及B通道数据对应的色差信号KB的值。
在具体实施中，所述第二图像处理窗口的大小可以根据经瞬态响应滤波后RGB数字图像的大小进行选择。例如，参照图4，将像素点NEW作为经瞬态响应滤波后RGB数字图像的第k个像素点，所述第一图像处理窗口的大小可以为以像素点NEW为中心的5*13的图像块。
在具体实施中，当经瞬态响应滤波后RGB数字图像中无法构建以某一像素点为中心的第二图像处理窗口时，可以对所述经瞬态响应滤波后RGB数字图像进行图像边界扩展或图像边界镜像处理，进而获得以该像素点为中心的第二图像处理窗口。比如，当该像素点为经瞬态响应滤波后RGB数字图像的第一个像素点或者最后一个像素点时，可以对所述RGB数字图像进行图像边界扩展或图像边界镜像处理，进而获得第一个像素点或者最后一个像素点为中心的第二图像处理窗口。
参照图5，像素点NEW水平方向上的像素点共12个，分别为像素点F5、F4、F3、F2、F1、F0、H0、H 1、H 2、H 3、H 4及H 5。像素点NEW垂直方向上的像素点共4个，分别为像素点S2、S1、A0及A1。将像素点NEW及其水平方向的12个像素点的像素值，作为对像素点NEW进行水平方向伪彩色滤波的输入数据。将像素点NEW及其垂直方向上的4个像素点的像素值，作为对像素点NEW进行垂直方向伪彩色滤波的输入数据。
在具体实施中，瞬态响应滤波后RGB数字图像所在的色彩空间为RGB空间，对其进行色彩空间转换，将其转换至KR+KB空间。以所述第k个像素点的坐标为(p，q)为例，所述第k个像素点R通道数据瞬态响应滤波最终结果为RTI(p，q)，所述第k个像素点B通道数据瞬态响应滤波最终结果为BTI(p，q)，所述第k个像素点原始输入的G通道数据GIN(p，q)，采用公式(3)可以得到所述第k个像素点R通道数据对应的色差信号KR的值KRTI(p，q)，采用公式(4)可以得到所述第k个像素点B通道数据对应的色差信号KB的值KBTI(p，q)：
KRTI(p，q)＝RTI(p，q)-GIN(p，q) (3)
KBTI(p，q)＝BTI(p，q)-GIN(p，q) (4)
参照上述获得第k个像素点色差信号KB及KR值的描述，可以获得第二图像处理窗口中第k个像素点水平及垂直方向各个像素点对应的色差信号KB及KR值。
本发明的实施例中，为例便于描述，KX表示水平或垂直方向相应色差信号，KX1表示水平方向色差信号KR或KB，KX2为垂直方向色差信号KR或KB。KX可以为KX1，也可以为KX2。
步骤42，基于所述第k个像素点相应方向上各个像素点的相应色差信号值，对所述第k个像素点相应方向上相应色差信号值进行伪彩色滤波。
在具体实施中，可以采用多种方法对所述第k个像素点相应方向上相应色差信号值进行伪彩色滤波，具体不作限制。对所述第k个像素点相应方向上相应色差信号值进行伪彩色滤波，可以得到所述第k个像素点水平方向上色差信号KR及KB值的伪彩色滤波结果KX1(p，q)，以及垂直方向上色差信号KR及KB值的伪彩色滤波结果KX2(p，q)。
步骤43，对所述第k个像素点相应方向上相应色差信号值的伪彩色滤波结果进行融合，得到所述第k个像素点相应色差信号值的伪彩色滤波结果。
在具体实施中，可以采用多种方法对所述第k个像素点相应方向上相应色差信号值的伪彩色滤波结果进行融合，具体不作限制。
在本发明的一实施例中，可以先对所述第k个像素点水平方向上相应色彩信号值的伪彩色滤波结果，与垂直方向上相应色差信号值的伪彩色滤波结果进行比较，进而将二者中较小者作为所述第k个像素点相应色差信号值的伪彩色滤波结果KX(p，q)＝min(KX1(p，q)，KX2(p，q))。
具体地，令第k个像素点水平方向色差信号KR值的伪彩色滤波结果为KR1(i，j)，垂直方向色差信号KR值的伪彩色滤波结果为KR2(p，q)，则第k个像素点色差信号KR值的伪彩色滤波结果KR(p，q)＝min(KR1(p，q)，KR2(p，q))。
令第k个像素点水平方向色差信号KB值的伪彩色滤波结果为KB1(p，q)，垂直方向色差信号KB值的伪彩色滤波结果为KB2(p，q)，则第k个像素点色差信号KB值的伪彩色滤波结果KB(p，q)＝min(KB1(p，q)，KB2(p，q))。
参照图6，本发明实施例还提供了一种对第k个像素点相应方向上相应色差信号值进行伪彩色滤波的方法，利用所述方法可以得到第k个像素点水平及垂直方向上色差信号KR及KB值的伪彩色滤波结果。
具体地，所述方法可以包括如下步骤：
步骤61，基于所述第k个像素点相应方向上各个像素点相应色差信号值的符号，对所述第k个像素点相应方向上各个像素点相应色差信号值进行滤波处理。
在具体实施中，可以采用如下方法对所述第k个像素点水平方向上各个像素点相应色差信号值进行滤波处理：
当KX1TI(p，q)>0时，CLIP_KX1TI(p，q+r1)＝min(KX1TI(p，q+r1)，KX1TI(p，q))；
当KX1TI(p，q)＜0时，CLIP_KX1TI(p，q+r1)＝max(KX1TI(p，q+r1)，KX1TI(p，q))；
当KX1TI(p，q)＝0时，CLIP_KX1TI(p，q+r1)＝KX1TI(p，q)。
其中，KX1TI(p，q)为第k个像素点水平方向上相应色差信号值，KX1TI(p，q+r1)为坐标为(p，q+r1)的像素点在水平方向上的相应色差信号值，CLIP_KX1TI(p，q+r1)为坐标为(p，q+r1)的像素点在水平方向上相应色差信号值的滤波处理结果，-L1≤r1≤L1且r1为整数，L1为第二图像处理窗口在水平方向的半径。
在具体实施中，可以采用如下方法对所述第k个像素点垂直方向上各个像素点相应色差信号值进行滤波处理：
当KX2TI(p，q)>0时，CLIP_KX2TI(p，q+r2)＝min(KX2TI(p，q+r2)，KX2TI(p，q))；
当KX2TI(p，q)＜0时，CLIP_KX2TI(p，q+r2)＝max(KX2TI(p，q+r2)，KX2TI(p，q))；
当KX2TI(p，q)＝0时，CLIP_KX2TI(p，q+r2)＝KX2TI(p，q)。
其中，KX2TI(p，q)为第k个像素点垂直方向上相应色差信号值，KX2TI(p，q+r1)为坐标为(p，q+r1)的像素点在垂直方向上的相应色差信号值，CLIP_KX2TI(p，q+r1)为坐标为(p，q+r1)的像素点在垂直方向上相应色差信号值的滤波处理结果，-L2≤r2≤L2且r1为整数，L2为第二图像处理窗口在垂直方向的半径。
以图5中像素点NEW作为所述第k个像素点为例，其中，L1＝6，L2＝2。
假设第k个像素点垂直方向上色差信号KR的值KR2TI(2，6)＞0，则像素点H1(2，8)对应的色差信号KR的值KR2TI(2，8)的滤波结果CLIP_KR2TI(2，8)＝min(KR2TI(2，8)，KR2TI(2，6))。
步骤62，计算所述第k个像素点相应方向上各个像素点相应色差信号值的权重系数。
在具体实施中，可以采用多种方法计算所述第k个像素点相应方向上各个像素点的权重系数，具体不作限制。
在本发明的一实施例中，可以先计算所述第k个像素点相应方向上各个像素点的R、B及G通道数据对应的梯度，再基于第k个像素点相应方向上各个像素点的R、B及G通道数据对应的梯度，计算所述第k个像素点相应方向上各个像素点的权重系数。
在本发明的另一实施例中，除计算第k个像素点相应方向上各个像素点的R、B及G通道数据对应的梯度外，还可以计算第k个像素点相应方向上各个像素点的Y通道数据对应的梯度，进而基于第k个像素点相应方向上各个像素点的R、B、G及Y通道数据对应的梯度，计算所述第k个像素点相应方向上各个像素点的权重系数。
具体地，以DX表示相应通道数据对应的梯度，其中，DX可以为DR、DB、DY或DG，采用如下公式可以计算第k个像素点水平方向上的像素点(p，q+r1)的梯度DX(p，q+r1)：
DX(p，q+r1)＝abs(X(p，q+r1)–X(p，q)) (5)
其中，X(p，q+r1)为像素点(p，q+r1)相应通道数据，X(p，q)为第k个像素点的相应通道数据，通过对X(p，q+r1)及X(p，q)的差值取绝对值，得到DX(p，q+r1)。
相似地，可以获得第k个像素点垂直方向上像素点(p，q+r2)的梯度DX(p，q+r2)。
获得第k个像素点相应方向上各个像素点的R、B及G通道数据的梯度后，可以采用公式(6)获得像素点(p，q+r1)色差信号KB值的权重系数wk1(p，q+r1)，采用公式(7)获得像素点(p，q+r1)色差信号KR值的权重系数wk2(p，q+r1)，采用公式(8)获得像素点(p，q+r2)色差信号KB值的权重系数wk1(p，q+r2)，采用公式(9)获得像素点(p，q+r2)色差信号KR值的权重系数wk1(p，q+r2)：
wk1(p，q+r1)＝1.0\(DG(p，q+r1)+DB(p，q+r1)) (6)
wk2(p，q+r1)＝1.0\(DG(p，q+r1)+DR(p，q+r1)) (7)
wk1(p，q+r2)＝1.0\(DG(p，q+r2)+DR(p，q+r2)) (8)
wk2(p，q+r2)＝1.0\(DG(p，q+r2)+DB(p，q+r2)) (9)
获得第k个像素点相应方向上各个像素点的R、B、G及Y通道数据的梯度后，可以采用公式(10)获得像素点(p，q+r1)色差信号KB值的权重系数wk1(p，q+r1)，采用公式(11)获得像素点(p，q+r1)色差信号KR值的权重系数wk2(p，q+r1)，采用公式(12)获得像素点(p，q+r2)色差信号KB值的权重系数wk1(p，q+r2)，采用公式(13)获得像素点(p，q+r2)色差信号KR值的权重系数wk1(p，q+r2)：
wk1(p，q+r1)＝1.0\(DG(p，q+r1)+DB(p，q+r1)+DY(p，q+r1)) (10)
wk2(p，q+r1)＝1.0\(DG(p，q+r1)+DB(p，q+r1)+DY(p，q+r1)) (11)
wk1(p，q+r2)＝1.0\(DG(p，q+r2)+DB(p，q+r2)+DY(p，q+r2)) (12)
wk2(p，q+r2)＝1.0\(DG(p，q+r2)+DB(p，q+r2)+DY(p，q+r2)) (13)
步骤63，基于所述第k个像素点相应方向上各个像素点相应色差信号值的滤波处理结果，以及所述各个像素点的权重系数，对所述第k个像素点相应方向相应色差信号值进行伪彩色滤波。
在本发明的一实施例中，可以采用如下公式获得所述第k个像素点相应方向相应色差信号值的伪彩色滤波结果：
其中，(p,q)为所述第k个像素点的坐标，KX(p,q)为第k个像素点相应方向色差信号KX值的伪彩色滤波结果，r∈[-L，L]且r为整数，L为所述第二图像处理窗口相应方向的半径，wk(p,q+r)为像素点(p,q+r)的权重系数，CLIP_KXTI(p,q+r)为像素点(p,q+r)X通道数据的滤波处理结果。
具体地，应用公式(14)可以得到：
第k个像素点水平方向色差信号KR值的伪彩色滤波结果KR1(p,q)：
第k个像素点水平方向色差信号KB值的伪彩色滤波结果KB1(p,q)：
第k个像素点垂直方向色差信号KB值的伪彩色滤波结果KR2(p,q)：
第k个像素点垂直方向色差信号KR值的伪彩色滤波结果KB2(p,q)：
参照图7，本发明实施例还提供了一种图像融合的方法，利用所述方法可以对经瞬态响应滤波后的RGB图像数据及经伪彩色滤波后的RGB图像数据进行融合。
具体地，所述方法可以包括如下步骤：
步骤71，计算所述待处理的RGB数字图像中各个像素点相应通道数据的局部对比度值。
在具体实施中，某一像素点的局部对比度值通常由该像素点的水平方向对比度值及垂直方向对比度值合并而成，其中具体合并方式可以包括多种。
在本发明的一实施例中，可以将该像素点水平方向相应通道数据局部对比度值及垂直方向相应通道数据局部对比度值中较大者，作为该像素点相应方向相应通道数据的局部对比度值。
具体地，以第一图像处理窗口为例，第n个像素点水平方向X通道的局部对比度值X1contrast(i，j)＝X1MAX-X1MIN，垂直方向X通道的局部对比度值X2contrast(i，j)＝X2MAX-X2MIN，则第n个像素点X通道的局部对比度值Xcontrast(i，j)＝max(X1contrast(i，j)，X2contrast(i，j))。其中，Xcontrast(i，j)可以为R通道数据的局部对比度值Rcontrast(i，j)，也可以为B通道数据的局部对比度值Bcontrast(i，j)。
步骤72，基于所述待处理的RGB数字图像中各个像素点相应通道数据的局部对比度值，计算相应像素点相应通道数据的融合参数。
在本发明的一实施例中，可以采用如下方法计算像素点(i，j)的R通道数据或B通道数据的融合参数Xalpha(i，j)：
其中，Xcontrast(i,j)为像素点(i，j)的局部对比度值，Xmax为像素点(i，j)X通道数据的最大值，Xmin为像素点(i，j)X通道数据的最小值，sigma1及sigma2为预设的调节参数，X为R通道或B通道。
在具体实施中，Xmax＝max(X1MAX，X2MAX)，其中，X1MAX为像素点(i，j)水平方向X通道数据的最大值，X2MAX为像素点(i，j)垂直方向X通道数据的最大值。
在具体实施中，Xmin＝max(X1MIN，X2MIN)，其中，X1MIN为像素点(i，j)水平方向X通道数据的最小值，X2MIN为像素点(i，j)垂直方向X通道数据的最小值。
应用公式(15)，可以得到像素点(i，j)R通道数据的融合参数Ralpha(i，j)：
以及得到像素点(i，j)B通道数据的融合参数Balpha(i，j)：
步骤73，基于所述融合参数，对相应像素点相应通道数据对应的瞬态响应滤波结果及伪彩色滤波结果进行融合，作为所述相应像素点去除伪彩色后相应通道的数据值。
在本发明的一实施例中，可以先将像素点(i，j)色差信号KX值的伪彩色滤波结果KX(i，j)转换至RGB色差空间，得到像素点(i，j)X通道数据的伪彩色滤波结果Xfcs(i,j)，再采用如下方法得到像素点(i，j)去除伪彩色后R通道数据或B通道数据的数据值Xfcsoutput(i，j)：
Xfcsoutput(i,j)＝(1-Xalpha(i,j))*XTI(i,j)+Xalpha(i,j)*Xfcs(i,j) (16)
其中，XTI(i,j)为像素点(i，j)X通道的瞬态响应滤波结果，Xfcs(i,j)为像素点(i，j)X通道数据的伪彩色滤波结果，Xfcs(i,j)＝KX(i，j)+GIN(i，j)，GIN(i，j)为像素点(i，j)的G通道数据。
具体地，应用公式(16)可以得到像素点(i，j)去除伪彩色后R通道数据的数据值Rfcsoutput(i,j)：
Rfcsoutput(i,j)＝(1-Ralpha(i,j))*RTI(i,j)+Ralpha(i,j)*KR(i,j)。
应用公式(17)可以得到像素点(i，j)去除伪彩色后B通道数据的数据值Bfcsoutput(i,j)：
Bfcsoutput(i,j)＝(1-Balpha(i,j))*BTI(i,j)+Balpha(i,j)*KB(i,j)。
基于待处理的RGB数字图像中各个像素点的局部对比度值，来计算像素点的融合参数，可以使得局部对比度强的区域(也就是伪彩色更容易产生的区域)更多使用伪彩色滤波的结果，而其他区域则倾向于使用瞬态响应滤波的结果，从而可以在更好地去除伪彩色的同时，保护数字图像的色彩平衡。
参照图8，本发明实施例还提供了另一种伪彩色滤波的方法，利用所述方法可以对经瞬态响应滤波后的RGB图像数据中第k个像素点进行水平方向及垂直方向的伪彩色滤波，进而得到经瞬态响应滤波后的RGB图像数据中各个像素点的伪彩色滤波结果。
具体地，所述方法可以包括如下步骤：
步骤81，以所述第k个像素点为中心的预设第二大小的图像块作为第二图像处理窗口，对所述经瞬态响应滤波后的RGB图像数据进行色彩空间转换，得到所述第k个像素点相应方向上各个像素点R通道数据对应的色差信号KR的值，以及B通道数据对应的色差信号KB的值。
关于步骤81，具体可以参照上述关于步骤41的描述，此处不再赘述。
步骤82，基于所述第k个像素点相应方向上各个像素点相应色差信号值的符号，对所述第k个像素点相应方向上各个像素点相应色差信号值进行滤波处理。
步骤83，计算所述第k个像素点相应方向上各个像素点相应色差信号值的权重系数。
步骤84，基于所述第k个像素点相应方向上各个像素点相应色差信号值的滤波处理结果，以及所述各个像素点的权重系数，对所述第k个像素点相应方向相应色差信号值进行伪彩色滤波。
关于步骤82～84，具体可以分别参照上述关于步骤61～63的描述，此处不作赘述。
步骤85，对所述第k个像素点相应方向上相应色差信号值的伪彩色滤波结果进行融合，得到所述第k个像素点相应色差信号值的伪彩色滤波结果。
关于步骤85，具体可以参照上述关于步骤43的描述，此处不再赘述。
步骤86，对所述经瞬态响应滤波后的RGB图像数据中各个像素点相应色差信号值进行中值滤波，得到所述经瞬态响应滤波后的RGB图像数据中各个像素点相应色差信号值的中值滤波结果。
在具体实施中，可以选取预设第三大小的图像块作为中值滤波的处理窗口，比如，可以选取3*3的图像块作为中值滤波的处理窗口，也可以选取5*5的图像块作为中值滤波的处理窗口，当然也可以选取其它大小的图像块作为中值滤波的处理窗口，具体处理窗口的大小不受限制。
采用中值滤波器对经瞬态响应滤波后的RGB图像数据中各个像素点相应色差信号值进行中值滤波，可以得到每个像素点水平及垂直方向色差信号KR及KB值的中值滤波结果。
以第k个像素点为例，对第k个像素点进行中值滤波后，可以获得第k个像素点水平方向色差信号KR值的中值滤波结果KR1-med(i，j)，垂直方向色差信号KR值的中值滤波结果KR2-med(i，j)，水平方向色差信号KB值的中值滤波结果KB1-med(i，j)，垂直方向色差信号KB值的中值滤波结果KB2-med(i，j)。
步骤87，将第k个像素点相应色差信号值的伪彩色滤波结果，与第k个像素点相应方向相应色差信号值的中值滤波结果进行融合，融合结果作为第k个像素点相应色差信号值伪彩色滤波的最终结果。
在具体实施中，可以采用多种方法对第k个像素点相应色差信号值的伪彩色滤波结果，与第k个像素点相应方向相应色差信号值的中值滤波结果进行融合，具体不作限制。
在本发明的一实施例中，可以采用中值滤波的方法对第k个像素点相应色差信号值的伪彩色滤波结果，与第k个像素点相应方向相应色差信号值的中值滤波结果进行融合。
比如，对第k个像素点水平方向色差信号KR值的中值滤波结果KR1-med(i，j)、垂直方向色差信号KR值的中值滤波结果KR2-med(i，j)、以及第k个像素点色差信号KR值的伪彩色滤波结果KR(i，j)进行中值滤波，即分别去掉三者中的最大值及最小值，将三者中的中间值作为第k个像素点相应色差信号值伪彩色滤波的最终结果KR-med(i，j)。
将KR-med(i，j)作为像素点(i，j)色差信号KX值的伪彩色滤波结果KX(i，j)，转换至RGB色差空间后，代入公式(16)，可以得到像素点(i，j)去除伪彩色后R通道的数据值。
通过对瞬态响应滤波后的RGB图像数据中各个像素点相应色差信号值进行中值滤波，进而基于中值滤波结果获得各个像素点相应色差信号值伪彩色滤波的最终结果，可以保持去除伪彩色后的数字图像边缘完整，提高图像质量。
参照图9，本发明实施例还提供了另一种数字图像处理方法，所述方法可以包括如下步骤：
步骤91，对待处理的RGB数字图像进行瞬态响应滤波；
步骤92，对所述经瞬态响应滤波后的RGB图像数据进行伪彩色滤波；
步骤93，对所述经瞬态响应滤波后的RGB图像数据及所述经伪彩色滤波后的RGB图像数据进行融合；
步骤94，输出所述待处理的RGB数字图像对应的去除伪彩色的RGB数字图像。
在具体实施中，可以采用如下步骤对待处理的RGB数字图像中第n个像素点进行水平方向瞬态响应滤波：
步骤911a，计算第一图像处理窗口中第n个像素点水平方向上R通道数据及B通道数据的最大值和最小值；
步骤912a，对第n个像素点水平方向上R通道数据及B通道数据进行瞬态响应滤波；
步骤913a，根据第n个像素点水平方向相应通道数据瞬态响应滤波的上限值及下限值，对第n个像素点水平方向上R通道数据及B通道数据进行瞬态响应滤波结果进行限制。
为了方便描述，将步骤911c简称为CLIP过程。
在具体实施中，可以采用如下步骤对待处理的RGB数字图像中第n个像素点进行垂直方向瞬态响应滤波：
步骤911b，计算第一图像处理窗口中第n个像素点垂直方向上R通道数据及B通道数据的最大值和最小值；
步骤912b，对第n个像素点垂直方向上R通道数据及B通道数据进行瞬态响应滤波；
步骤913b，根据第n个像素点垂直方向相应通道数据瞬态响应滤波的上限值及下限值，对第n个像素点垂直方向上R通道数据及B通道数据进行瞬态响应滤波结果进行限制。
获得第n个像素点垂直方向相应通道数据瞬态响应滤波结果，以及水平方向相应通道数据瞬态响应滤波结果后，执行步骤914：
步骤914，对所述第n个像素点水平方向及垂直方向中相应通道数据瞬态响应滤波的最终结果进行融合，得到所述第n个像素点相应通道数据的瞬态响应滤波结果。
在具体实施中，可以采用如下步骤对经瞬态响应滤波后的RGB图像数据中第k个像素点进行水平方向及垂直方向伪彩色滤波：
步骤921，将所述经瞬态响应滤波后的RGB图像数据转换至KR/KB色彩空间。
在KR/KB色彩空间中，执行中值滤波操作及伪彩色滤波操作。
其中，对KR/KB色彩空间中第k个像素点(p，q)执行中值滤波操作，可以包括如下步骤：
步骤922a，对第k个像素点(p，q)水平方向色差信号KR及KB值进行中值滤波，得到k个像素点(p，q)水平方向色差信号KR值的中值滤波结果KR1-med(i，j)，以及水平方向色差信号KR值的中值滤波结果KB1-med(i，j)。
步骤922b，对第k个像素点(p，q)垂直方向色差信号KR及KB值进行中值滤波，得到k个像素点(p，q)垂直方向色差信号KR值的中值滤波结果KR2-med(i，j)，以及水平方向色差信号KR值的中值滤波结果KB2-med(i，j)。
对KR/KB色彩空间中第k个像素点(p，q)执行伪彩色滤波操作，可以包括如下步骤：
步骤923a，对第k个像素点(p，q)水平方向色差信号KR及KB值伪彩色滤波，得到伪彩色滤波结果KR1(p，q)及KB1(p，q)。
步骤923b，对第k个像素点(p，q)垂直方向色差信号KR及KB值伪彩色滤波，得到伪彩色滤波结果KR2(p，q)及KB2(p，q)。
步骤923c，融合第k个像素点(p，q)水平及垂直方向伪彩色滤波结果。
具体地，对第k个像素点(p，q)水平及垂直方向色差信号KR值的伪彩色滤波结果进行融合，得到第k个像素点(p，q)色差信号KR值的伪彩色滤波结果KR(p，q)；以及对第k个像素点(p，q)水平及垂直方向色差信号KB值的伪彩色滤波结果进行融合，得到第k个像素点(p，q)色差信号KB值的伪彩色滤波结果KB(p，q)。
获得第k个像素点(p，q)相应色差信号值的伪彩色滤波结果及中值滤波结果后，执行步骤924。
步骤924，融合三种计算结果。
具体地，将第k个像素点(p，q)色差信号KR值的伪彩色滤波结果KR(p，q)、水平方向色差信号KR值的中值滤波结果KR1-med(i，j)，与k个像素点(p，q)垂直方向色差信号KR值的中值滤波结果KR2-med(i，j)进行融合，得到第k个像素点(p，q)色差信号KR值的最终伪彩色滤波结果KR-med(i，j)。
将第k个像素点(p，q)色差信号KB值的伪彩色滤波结果KB(p，q)、水平方向色差信号KB值的中值滤波结果KB1-med(i，j)，与k个像素点(p，q)垂直方向色差信号KB值的中值滤波结果KB2-med(i，j)进行融合，得到第k个像素点(p，q)色差信号KB值的最终伪彩色滤波结果KB-med(i，j)。
步骤925，将第k个像素点(p，q)色差信号KB值的最终伪彩色滤波结果KB-med(i，j)及色差信号KR值的最终伪彩色滤波结果KR-med(i，j)转换至RGB色彩空间。
由上述内容可知，本发明实施例中的数字图像处理方法，对待处理的RGB数字图像整体去除伪彩色，可以避免造成图像色度域的信息的损失。并且，先对待处理的RGB数字图像进行瞬态响应滤波后，再进行伪彩色滤波，可以压缩伪彩色在局部图像上的范围，在需要去除同样宽度的伪彩色的前提下，有效减小计算量。
为了使本领域技术人员更好地理解和实现本发明，以下对上述数据业务的处理方法对应的移动终端进行详细描述。
参照图10，本发明实施例还提供了一种数字图像处理装置100，所述数字图像处理装置100可以包括：瞬态响应滤波器101，伪彩色滤波器102及融合处理器103。其中：
所述瞬态响应滤波器101，适于对待处理的RGB数字图像进行瞬态响应滤波，得到经瞬态响应滤波后的RGB图像数据；
所述伪彩色滤波器102，适于对所述经瞬态响应滤波后的RGB图像数据进行伪彩色滤波，得到经伪彩色滤波后的RGB图像数据；
所述融合处理器103，适于对所述经瞬态响应滤波后的RGB图像数据及所述经伪彩色滤波后的RGB图像数据进行融合，输出所述待处理的RGB数字图像对应的去除伪彩色的RGB数字图像。
在本发明的一实施例中，参照图11，所述瞬态响应滤波器101，适于对所述待处理的RGB数字图像中第n个像素点进行瞬态响应滤波，n为正整数。
具体地，所述瞬态响应滤波器101可以包括：第一计算单元111，第二计算单元112以及第一融合单元113。其中：
所述第一计算单元111，适于以所述第n个像素点为中心的预设第一大小的图像块作为第一图像处理窗口，分别计算所述第一图像处理窗口中所述第n个像素点水平方向上以及垂直方向上R通道数据及B通道数据的最大值和最小值；
所述第二计算单元112，适于分别计算所述第n个像素点水平方向上以及垂直方向上R通道数据及B通道数据瞬态响应滤波的最终结果。
所述第一融合单元113，适于分别对所述第n个像素点水平方向及垂直方向中相应通道数据瞬态响应滤波的最终结果进行融合，得到所述第n个像素点相应通道数据的瞬态响应滤波结果。
在具体实施中，所述第二计算单元112可以包括：瞬态响应滤波子单元1121、第一计算子单元1122及处理子单元1123，其中：
所述瞬态响应滤波子单元1121，适于基于所述第n个像素点相应方向上相应通道数据的最大值和最小值，对所述第n个像素点进行相应方向相应通道数据的瞬态响应滤波；
所述第一计算子单元1122，适于计算所述第n个像素点相应方向相应通道数据瞬态响应滤波的上限值及下限值；
所述处理子单元1123，适于基于计算结果对第n个像素点相应方向相应通道数据瞬态响应滤波的结果进行限制，得到所述第n个像素点相应方向相应通道数据瞬态响应滤波的最终结果。
在具体实施中，所述第一计算单元111，适于计算相应方向上，位于所述第n个像素点两侧的像素点中相应通道数据的最大值及最小值；将所述第n个像素点第一侧相应通道数据的最大值与第二侧相应通道数据的最小值，和第二侧相应通道数据的最大值与第一侧相应通道数据的最小值中差异大的一组所对应的最大值及最小值，分别作为所述第n个像素点相应方向上相应通道数据的最大值及最小值。
在具体实施中，所述瞬态响应滤波子单元1121，适于使用预设的瞬态响应滤波系数，对所述相应方向上相应通道数据的最大值和最小值及所述第n个像素点进行卷积运算。
在具体实施中，所述瞬态响应滤波子单元1121，适于当XIN(i，j)>GIN(i，j)时，采用如下方法获得第n个像素点进行相应方向相应通道数据的瞬态响应滤波的结果：XPF(i，j)＝[p0，p1，p2]*[XMAX，XIN(i，j)，XMIN]T；当XIN(i，j)≤GIN(i，j)时，采用如下方法获得第n个像素点进行相应方向相应通道数据的瞬态响应滤波的结果：XPF(i，j)＝[p0，p1，p2]*[XMIN，XIN(i，j)，XMAX]T；其中，[p0，p1，p2]为所述瞬态响应滤波系数，X为R通道或B通道，XMAX为相应方向上R通道数据的最大值或者为相应方向上B通道数据的最大值，(i，j)为第n个像素点的坐标，XIN(i，j)为第n个像素点R通道或B通道数据，GIN(i，j)为第n个像素点的G通道数据，XPF(i，j)为第n个像素点R通道或B通道数据的瞬态响应滤波结果。
在发明的一实施例中，[p0，p1，p2]＝[-0.25，1.375，-0.125]。
在具体实施中，所述第一计算子单元1122，适于根据所述第n个像素点相应方向相应通道数据与所述第n个像素点G通道数据的比较结果，确定所述第n个像素点相应方向相应通道数据瞬态响应滤波的上限值及下限值。
在具体实施中，所述第一融合单元113，适于对所述第n个像素点水平方向上相应通道数据瞬态响应滤波的最终结果，与垂直方向上相应通道数据瞬态响应滤波的最终结果进行比较，将二者中较小者作为所述第n个像素点相应通道数据的瞬态响应滤波结果。
在本发明的一实施例中，参照图12，所述伪彩色滤波器102，适于对所述经瞬态响应滤波后的RGB图像数据中第k个像素点进行水平方向及垂直方向伪彩色滤波，k为正整数。
具体地，所述伪彩色滤波器102可以包括：色彩空间转换单元121，伪彩色滤波单元122及第二融合单元123。其中：
所述色彩空间转换单元121，适于以所述第k个像素点为中心的预设第二大小的图像块作为第二图像处理窗口，对所述经瞬态响应滤波后的RGB图像数据进行色彩空间转换，得到所述第k个像素点相应方向上各个像素点R通道数据对应的色差信号KR的值，以及B通道数据对应的色差信号KB的值；
所述伪彩色滤波单元122，适于基于所述第k个像素点相应方向上各个像素点的相应色差信号值，对所述第k个像素点相应方向上相应色差信号值进行伪彩色滤波；
所述第二融合单元123，适于对所述第k个像素点相应方向上相应色差信号值的伪彩色滤波结果进行融合，得到所述第k个像素点相应色差信号值的伪彩色滤波结果。
在具体实施中，所述伪彩色滤波单元122可以包括：第一滤波子单元1221，第二计算子单元1222及第二滤波子单元1223。其中：
所述第一滤波子单元1221，适于基于所述第k个像素点相应方向上各个像素点相应色差信号值的符号，对所述第k个像素点相应方向上各个像素点相应色差信号值进行滤波处理；
所述第二计算子单元1222，适于计算所述第k个像素点相应方向上各个像素点相应色差信号值的权重系数；
所述第二滤波子单元1223，适于基于所述第k个像素点相应方向上各个像素点相应色差信号值的滤波处理结果，以及所述各个像素点的权重系数，对所述第k个像素点相应方向上相应色差信号值进行相应方向伪彩色滤波。
在本发明的一实施例中，所述第二计算子单元1222，适于计算所述第k个像素点相应方向上各个像素点的R、B及G通道数据对应的梯度；基于第k个像素点相应方向上各个像素点的R、B及G通道数据对应的梯度，计算所述第k个像素点相应方向上各个像素点的权重系数。
在本发明的另一实施例中，所述第二计算子单元，适于计算所述第k个像素点相应方向上各个像素点的R、B、G及Y通道数据对应的梯度；基于第k个像素点相应方向上各个像素点的R、B、G及Y通道数据对应的梯度，计算所述第k个像素点相应方向上各个像素点的权重系数。
在具体实施中，所述第二滤波子单元1223，适于采用如下方法获得所述第k个像素点相应方向相应色差信号值的伪彩色滤波结果：
其中，(p,q)为所述第k个像素点的坐标，KX(p,q)为第k个像素点相应方向色差信号KX值的伪彩色滤波结果，r∈[-L，L]且r为整数，L为所述第二图像处理窗口相应方向的半径，wk(p,q+r)为像素点(p,q+r)的权重系数，CLIP_KXTI(p,q+r)为像素点(p,q+r)X通道数据的滤波处理结果。
在本发明的一实施例中，参照图12，所述伪彩色滤波器，还可以包括：第三滤波单元124及第三融合单元125，其中：
所述第三滤波单元124，对所述第k个像素点相应方向上各个像素点相应色差信号值进行中值滤波，得到所述经瞬态响应滤波后的RGB图像数据中各个像素点相应色差信号值的中值滤波结果；
所述第三融合单元125，适于将第k个像素点相应色差信号值的伪彩色滤波结果，与第k个像素点相应方向相应色差信号值的中值滤波结果进行融合，融合结果作为第k个像素点相应色差信号值伪彩色滤波的最终结果。
在本发明的一实施例中，参照图13，所述融合处理器103可以包括：第三计算单元131，第四计算单元132及第四融合单元133。其中：
所述第三计算单元131，适于计算所述待处理的RGB数字图像中各个像素点相应通道数据的局部对比度值；
所述第四计算单元132，适于基于所述待处理的RGB数字图像中各个像素点相应通道数据的局部对比度值，计算相应像素点相应通道数据的融合参数；
所述第四融合单元133，适于基于所述融合参数，对相应像素点相应通道数据对应的瞬态响应滤波结果及伪彩色滤波结果进行融合，作为所述相应像素点去除伪彩色后相应通道的数据值。
在本发明的一实施例中，所述第四计算单元132，适于采用如下方法计算像素点(h,k)的X通道数据的融合参数Xalpha(i，j)：
其中，Xcontrast(i,j)为像素点(i，j)的局部对比度值，Xmax为像素点(i，j)X通道数据的最大值，Xmin为像素点(i，j)X通道数据的最小值，sigma1及sigma2为预设的调节参数，X为R通道或B通道。
在本发明的一实施例中，所述第四融合单元133，适于采用如下方法得到像素点(h,k)去除伪彩色后X通道数据的数据值Xfcsoutput(h,k)：
Xfcsoutput(i,j)＝(1-Xalpha(i,j))*XTI(i,j)+Xalpha(i,j)*Xfcs(i,j)
其中，XTI(i,j)为像素点(i，j)X通道的瞬态响应滤波结果，Xfcs(i,j)为像素点(i，j)X通道数据的伪彩色滤波结果，Xfcs(i,j)＝KX(i，j)+GIN(i，j)，GIN(i，j)为像素点(i，j)的G通道数据。
本发明实施例中的数字图像处理装置100，不仅可以有效去除数字图像中的伪彩色，而且可以保留数字图像低饱和度色彩信息，并且实现代价上也是经济可行的。
本领域普通技术人员可以理解上述实施例的各种方法中的全部或部分步骤是可以通过程序来指令相关的硬件来完成，该程序可以存储于一计算机可读存储介质中，存储介质可以包括：ROM、RAM、磁盘或光盘等。
虽然本发明披露如上，但本发明并非限定于此。任何本领域技术人员，在不脱离本发明的精神和范围内，均可作各种更动与修改，因此本发明的保护范围应当以权利要求所限定的范围为准。

Claims (34)

1.一种数字图像处理方法，其特征在于，包括：
对待处理的RGB数字图像进行瞬态响应滤波，得到经瞬态响应滤波后的RGB图像数据；
对所述经瞬态响应滤波后的RGB图像数据进行伪彩色滤波，得到经伪彩色滤波后的RGB图像数据；
对所述经瞬态响应滤波后的RGB图像数据及所述经伪彩色滤波后的RGB图像数据进行融合，输出融合后的RGB数字图像。
2.如权利要求1所述的数字图像处理方法，其特征在于，所述对待处理的RGB数字图像进行瞬态响应滤波，包括：
采用如下步骤对所述待处理的RGB数字图像中第n个像素点进行瞬态响应滤波：
以所述第n个像素点为中心的预设第一大小的图像块作为第一图像处理窗口，分别计算所述第一图像处理窗口中所述第n个像素点水平方向上以及垂直方向上R通道数据及B通道数据的最大值和最小值；
分别计算所述第n个像素点水平方向上以及垂直方向上R通道数据及B通道数据瞬态响应滤波的最终结果，其中，任一方向上的任一通道数据均采用如下方式计算：基于所述第n个像素点相应方向上相应通道数据的最大值和最小值，对所述第n个像素点进行相应方向相应通道数据的瞬态响应滤波；计算所述第n个像素点相应方向相应通道数据瞬态响应滤波的上限值及下限值，并基于计算结果对第n个像素点相应方向相应通道数据瞬态响应滤波的结果进行限制，得到所述第n个像素点相应方向相应通道数据瞬态响应滤波的最终结果；
分别对所述第n个像素点水平方向及垂直方向中相应通道数据瞬态响应滤波的最终结果进行融合，得到所述第n个像素点相应通道数据的瞬态响应滤波结果；
其中，n为正整数。
3.如权利要求2所述的数字图像处理方法，其特征在于，所述计算所述第一图像处理窗口中所述第n个像素点水平方向上以及垂直方向上R通道数据及B通道数据的最大值和最小值，包括：
计算相应方向上，位于所述第n个像素点两侧的像素点中相应通道数据的最大值及最小值；
将所述第n个像素点第一侧相应通道数据的最大值与第二侧相应通道数据的最小值，和第二侧相应通道数据的最大值与第一侧相应通道数据的最小值中差异大的一组所对应的最大值及最小值，分别作为所述第n个像素点相应方向上相应通道数据的最大值及最小值。
4.如权利要求2所述的数字图像处理方法，其特征在于，所述基于所述第n个像素点相应方向上相应通道数据的最大值和最小值，对所述第n个像素点进行相应方向相应通道数据的瞬态响应滤波，包括：
使用预设的瞬态响应滤波系数，对所述相应方向上相应通道数据的最大值和最小值及所述第n个像素点进行卷积运算。
5.如权利要求4所述的数字图像处理方法，其特征在于，所述使用预设的瞬态响应滤波系数，对所述相应方向上相应通道数据的最大值和最小值及所述第n个像素点进行卷积运算，包括：
当XIN(i，j)>GIN(i，j)时，采用如下方法获得第n个像素点进行相应方向相应通道数据的瞬态响应滤波的结果：XPF(i，j)＝[p0，p1，p2]*[XMAX，XIN(i，j)，XMIN]T；
当XIN(i，j)≤GIN(i，j)时，采用如下方法获得第n个像素点进行相应方向相应通道数据的瞬态响应滤波的结果：XPF(i，j)＝[p0，p1，p2]*[XMIN，XIN(i，j)，XMAX]T；
其中，[p0，p1，p2]为所述瞬态响应滤波系数，X为R通道或B通道，XMAX为相应方向上R通道数据的最大值或者为相应方向上B通道数据的最大值，(i，j)为第n个像素点的坐标，XIN(i，j)为第n个像素点R通道或B通道数据，GIN(i，j)为第n个像素点的G通道数据，XPF(i，j)为第n个像素点R通道或B通道数据的瞬态响应滤波结果；XMIN为相应方向上R通道数据的最小值或者为相应方向上B通道数据的最小值。
6.如权利要求5所述的数字图像处理方法，其特征在于，[p0，p1，p2]＝[-0.25，1.375，-0.125]。
7.如权利要求2所述的数字图像处理方法，其特征在于，所述计算所述第n个像素点相应方向相应通道数据瞬态响应滤波的上限值及下限值，包括：
根据所述第n个像素点相应方向相应通道数据与所述第n个像素点G通道数据的比较结果，确定所述第n个像素点相应方向相应通道数据瞬态响应滤波的上限值及下限值。
8.如权利要求2所述的数字图像处理方法，其特征在于，所述对所述第n个像素点水平方向及垂直方向中相应通道数据瞬态响应滤波的最终结果进行融合，包括：
对所述第n个像素点水平方向上相应通道数据瞬态响应滤波的最终结果，与垂直方向上相应通道数据瞬态响应滤波的最终结果进行比较，将二者中较小者作为所述第n个像素点相应通道数据的瞬态响应滤波结果。
9.如权利要求1所述的数字图像处理方法，其特征在于，所述对所述经瞬态响应滤波后的RGB图像数据进行伪彩色滤波，包括：
采用如下步骤对所述经瞬态响应滤波后的RGB图像数据中第k个像素点进行水平方向及垂直方向伪彩色滤波，k为正整数：
以所述第k个像素点为中心的预设第二大小的图像块作为第二图像处理窗口，对所述经瞬态响应滤波后的RGB图像数据进行色彩空间转换，得到所述第k个像素点相应方向上各个像素点R通道数据对应的色差信号KR的值，以及B通道数据对应的色差信号KB的值；
基于所述第k个像素点相应方向上各个像素点的相应色差信号值，对所述第k个像素点相应方向上相应色差信号值进行伪彩色滤波；
对所述第k个像素点相应方向上相应色差信号值的伪彩色滤波结果进行融合，得到所述第k个像素点相应色差信号值的伪彩色滤波结果。
10.如权利要求9所述的数字图像处理方法，其特征在于，所述基于所述第k个像素点相应方向上各个像素点的相应色差信号值，对所述第k个像素点相应方向上相应色差信号值进行伪彩色滤波，包括：
基于所述第k个像素点相应方向上各个像素点相应色差信号值的符号，对所述第k个像素点相应方向上各个像素点相应色差信号值进行滤波处理；
计算所述第k个像素点相应方向上各个像素点相应色差信号值的权重系数；
基于所述第k个像素点相应方向上各个像素点相应色差信号值的滤波处理结果，以及所述各个像素点的权重系数，对所述第k个像素点相应方向上相应色差信号值进行相应方向伪彩色滤波。
11.如权利要求10所述的数字图像处理方法，其特征在于，所述计算所述第k个像素点相应方向上各个像素点相应色差信号值的权重系数，包括：
计算所述第k个像素点相应方向上各个像素点的R、B及G通道数据对应的梯度；
基于第k个像素点相应方向上各个像素点的R、B及G通道数据对应的梯度，计算所述第k个像素点相应方向上各个像素点相应色差信号值的权重系数。
12.如权利要求10所述的数字图像处理方法，其特征在于，所述计算所述第k个像素点相应方向上各个像素点相应色差信号值的权重系数，包括：
计算所述第k个像素点相应方向上各个像素点的R、B、G及Y通道数据对应的梯度；
基于第k个像素点相应方向上各个像素点的R、B、G及Y通道数据对应的梯度，计算所述第k个像素点相应方向上各个像素点相应色差信号值的权重系数。
13.如权利要求10所述的数字图像处理方法，其特征在于，所述基于所述第k个像素点相应方向上各个像素点相应色差信号值的滤波处理结果，以及所述各个像素点的权重系数，对所述第k个像素点相应方向上相应色差信号值进行相应方向伪彩色滤波，包括：
采用如下方法获得所述第k个像素点相应方向相应色差信号值的伪彩色滤波结果：
其中，(p,q)为所述第k个像素点的坐标，KX(p,q)为第k个像素点相应方向色差信号KX值的伪彩色滤波结果，r∈[-L，L]且r为整数，L为所述第二图像处理窗口相应方向的半径，wk(p,q+r)为像素点(p,q+r)的权重系数，CLIP_KXTI(p,q+r)为像素点(p,q+r)X通道数据的滤波处理结果。
14.如权利要求9所述的数字图像处理方法，其特征在于，所述对所述经瞬态响应滤波后的RGB图像数据中第k个像素点进行水平方向及垂直方向伪彩色滤波，还包括：
对所述第k个像素点相应方向上各个像素点相应色差信号值进行中值滤波，得到所述经瞬态响应滤波后的RGB图像数据中各个像素点相应色差信号值的中值滤波结果；
将第k个像素点相应色差信号值的伪彩色滤波结果，与第k个像素点相应方向相应色差信号值的中值滤波结果进行融合，融合结果作为第k个像素点相应色差信号值伪彩色滤波的最终结果。
15.如权利要求1所述的数字图像处理方法，其特征在于，所述对所述经瞬态响应滤波后的RGB图像数据及所述经伪彩色滤波后的RGB图像数据进行融合，包括：
计算所述待处理的RGB数字图像中各个像素点相应通道数据的局部对比度值；
基于所述待处理的RGB数字图像中各个像素点相应通道数据的局部对比度值，计算相应像素点相应通道数据的融合参数；
基于所述融合参数，对相应像素点相应通道数据对应的瞬态响应滤波结果及伪彩色滤波结果进行融合，作为所述相应像素点去除伪彩色后相应通道的数据值。
16.如权利要求15所述的数字图像处理方法，其特征在于，所述基于所述待处理的RGB数字图像中各个像素点相应通道数据的局部对比度值，计算相应像素点相应通道数据的融合参数，包括：
采用如下方法计算像素点(h,k)的X通道数据的融合参数Xalpha(i，j)：
其中，Xcontrast(i,j)为像素点(i，j)的局部对比度值，Xmax为像素点(i，j)X通道数据的最大值，Xmin为像素点(i，j)X通道数据的最小值，sigma1及sigma2为预设的调节参数，X为R通道或B通道。
17.如权利要求15所述的数字图像处理方法，其特征在于，所述基于所述融合参数，对相应像素点相应通道数据对应的瞬态响应滤波结果及伪彩色滤波结果进行融合，包括：
采用如下方法得到像素点(h,k)去除伪彩色后X通道数据的数据值Xfcsoutput(h,k)：
Xfcsoutput(i,j)＝(1-Xalpha(i,j))*XTI(i,j)+Xalpha(i,j)*Xfcs(i,j)
其中，XTI(i,j)为像素点(i，j)X通道的瞬态响应滤波结果，Xfcs(i,j)为像素点(i，j)X通道数据的伪彩色滤波结果，Xfcs(i,j)＝KX(i，j)+GIN(i，j)，GIN(i，j)为像素点(i，j)的G通道数据；Xalpha(i，j)为像素点(i，j)相应通道数据的融合参数Xalpha(i，j)；KX(i，j)为像素点(i，j)色差信号KX值的伪彩色滤波结果。
18.一种数字图像处理装置，其特征在于，包括：
瞬态响应滤波器，适于对待处理的RGB数字图像进行瞬态响应滤波，得到经瞬态响应滤波后的RGB图像数据；
伪彩色滤波器，适于对所述经瞬态响应滤波后的RGB图像数据进行伪彩色滤波，得到经伪彩色滤波后的RGB图像数据；
融合处理器，适于对所述经瞬态响应滤波后的RGB图像数据及所述经伪彩色滤波后的RGB图像数据进行融合，输出融合后的RGB数字图像。
19.如权利要求18所述的数字图像处理装置，其特征在于，所述瞬态响应滤波器，适于对所述待处理的RGB数字图像中第n个像素点进行瞬态响应滤波，n为正整数，包括：
第一计算单元，适于以所述第n个像素点为中心的预设第一大小的图像块作为第一图像处理窗口，分别计算所述第一图像处理窗口中所述第n个像素点水平方向上以及垂直方向上R通道数据及B通道数据的最大值和最小值；
第二计算单元，适于分别计算所述第n个像素点水平方向上以及垂直方向上R通道数据及B通道数据瞬态响应滤波的最终结果；包括：瞬态响应滤波子单元、第一计算子单元及处理子单元，其中：
所述瞬态响应滤波子单元，适于基于所述第n个像素点相应方向上相应通道数据的最大值和最小值，对所述第n个像素点进行相应方向相应通道数据的瞬态响应滤波；
所述第一计算子单元，适于计算所述第n个像素点相应方向相应通道数据瞬态响应滤波的上限值及下限值；
所述处理子单元，适于基于计算结果对第n个像素点相应方向相应通道数据瞬态响应滤波的结果进行限制，得到所述第n个像素点相应方向相应通道数据瞬态响应滤波的最终结果；
第一融合单元，适于分别对所述第n个像素点水平方向及垂直方向中相应通道数据瞬态响应滤波的最终结果进行融合，得到所述第n个像素点相应通道数据的瞬态响应滤波结果。
20.如权利要求19所述的数字图像处理装置，其特征在于，所述第一计算单元，适于计算相应方向上，位于所述第n个像素点两侧的像素点中相应通道数据的最大值及最小值；将所述第n个像素点第一侧相应通道数据的最大值与第二侧相应通道数据的最小值，和第二侧相应通道数据的最大值与第一侧相应通道数据的最小值中差异大的一组所对应的最大值及最小值，分别作为所述第n个像素点相应方向上相应通道数据的最大值及最小值。
21.如权利要求19所述的数字图像处理装置，其特征在于，所述瞬态响应滤波子单元，适于使用预设的瞬态响应滤波系数，对所述相应方向上相应通道数据的最大值和最小值及所述第n个像素点进行卷积运算。
22.如权利要求21所述的数字图像处理装置，其特征在于，所述瞬态响应滤波子单元，适于当XIN(i，j)>GIN(i，j)时，采用如下方法获得第n个像素点进行相应方向相应通道数据的瞬态响应滤波的结果：XPF(i，j)＝[p0，p1，p2]*[XMAX，XIN(i，j)，XMIN]T；当XIN(i，j)≤GIN(i，j)时，采用如下方法获得第n个像素点进行相应方向相应通道数据的瞬态响应滤波的结果：XPF(i，j)＝[p0，p1，p2]*[XMIN，XIN(i，j)，XMAX]T；其中，[p0，p1，p2]为所述瞬态响应滤波系数，X为R通道或B通道，XMAX为相应方向上R通道数据的最大值或者为相应方向上B通道数据的最大值，(i，j)为第n个像素点的坐标，XIN(i，j)为第n个像素点R通道或B通道数据，GIN(i，j)为第n个像素点的G通道数据，XPF(i，j)为第n个像素点R通道或B通道数据的瞬态响应滤波结果；XMIN为相应方向上R通道数据的最小值或者为相应方向上B通道数据的最小值。
23.如权利要求22所述的数字图像处理装置，其特征在于，[p0，p1，p2]＝[-0.25，1.375，-0.125]。
24.如权利要求19所述的数字图像处理装置，其特征在于，所述第一计算子单元，适于根据所述第n个像素点相应方向相应通道数据与所述第n个像素点G通道数据的比较结果，确定所述第n个像素点相应方向相应通道数据瞬态响应滤波的上限值及下限值。
25.如权利要求19所述的数字图像处理装置，其特征在于，所述第一融合单元，适于对所述第n个像素点水平方向上相应通道数据瞬态响应滤波的最终结果，与垂直方向上相应通道数据瞬态响应滤波的最终结果进行比较，将二者中较小者作为所述第n个像素点相应通道数据的瞬态响应滤波结果。
26.如权利要求18所述的数字图像处理装置，其特征在于，所述伪彩色滤波器，适于对所述经瞬态响应滤波后的RGB图像数据中第k个像素点进行水平方向及垂直方向伪彩色滤波，k为正整数，包括：
色彩空间转换单元，适于以所述第k个像素点为中心的预设第二大小的图像块作为第二图像处理窗口，对所述经瞬态响应滤波后的RGB图像数据进行色彩空间转换，得到所述第k个像素点相应方向上各个像素点R通道数据对应的色差信号KR的值，以及B通道数据对应的色差信号KB的值；
伪彩色滤波单元，适于基于所述第k个像素点相应方向上各个像素点的相应色差信号值，对所述第k个像素点相应方向上相应色差信号值进行伪彩色滤波；
第二融合单元，适于对所述第k个像素点相应方向上相应色差信号值的伪彩色滤波结果进行融合，得到所述第k个像素点相应色差信号值的伪彩色滤波结果。
27.如权利要求26所述的数字图像处理装置，其特征在于，所述伪彩色滤波单元包括：
第一滤波子单元，适于基于所述第k个像素点相应方向上各个像素点相应色差信号值的符号，对所述第k个像素点相应方向上各个像素点相应色差信号值进行滤波处理；
第二计算子单元，适于计算所述第k个像素点相应方向上各个像素点相应色差信号值的权重系数；
第二滤波子单元，适于基于所述第k个像素点相应方向上各个像素点相应色差信号值的滤波处理结果，以及所述各个像素点的权重系数，对所述第k个像素点相应方向上相应色差信号值进行相应方向伪彩色滤波。
28.如权利要求27所述的数字图像处理装置，其特征在于，所述第二计算子单元，适于计算所述第k个像素点相应方向上各个像素点的R、B及G通道数据对应的梯度；基于第k个像素点相应方向上各个像素点的R、B及G通道数据对应的梯度，计算所述第k个像素点相应方向上各个像素点相应色差信号值的权重系数。
29.如权利要求27所述的数字图像处理装置，其特征在于，所述第二计算子单元，适于计算所述第k个像素点相应方向上各个像素点的R、B、G及Y通道数据对应的梯度；基于第k个像素点相应方向上各个像素点的R、B、G及Y通道数据对应的梯度，计算所述第k个像素点相应方向上各个像素点相应色差信号值的权重系数。
30.如权利要求27所述的数字图像处理装置，其特征在于，所述第二滤波子单元，适于采用如下方法获得所述第k个像素点相应方向相应色差信号值的伪彩色滤波结果：
其中，(p,q)为所述第k个像素点的坐标，KX(p,q)为第k个像素点相应方向色差信号KX值的伪彩色滤波结果，r∈[-L，L]且r为整数，L为所述第二图像处理窗口相应方向的半径，wk(p,q+r)为像素点(p,q+r)的权重系数，CLIP_KXTI(p,q+r)为像素点(p,q+r)X通道数据的滤波处理结果。
31.如权利要求26所述的数字图像处理装置，其特征在于，所述伪彩色滤波器，还包括：
第三滤波单元，对所述第k个像素点相应方向上各个像素点相应色差信号值进行中值滤波，得到所述经瞬态响应滤波后的RGB图像数据中各个像素点相应色差信号值的中值滤波结果；
第三融合单元，适于将第k个像素点相应色差信号值的伪彩色滤波结果，与第k个像素点相应方向相应色差信号值的中值滤波结果进行融合，融合结果作为第k个像素点相应色差信号值伪彩色滤波的最终结果。
32.如权利要求18所述的数字图像处理装置，其特征在于，所述融合处理器，包括：
第三计算单元，适于计算所述待处理的RGB数字图像中各个像素点相应通道数据的局部对比度值；
第四计算单元，适于基于所述待处理的RGB数字图像中各个像素点相应通道数据的局部对比度值，计算相应像素点相应通道数据的融合参数；
第四融合单元，适于基于所述融合参数，对相应像素点相应通道数据对应的瞬态响应滤波结果及伪彩色滤波结果进行融合，作为所述相应像素点去除伪彩色后相应通道的数据值。
33.如权利要求32所述的数字图像处理装置，其特征在于，所述第四计算单元，适于采用如下方法计算像素点(h,k)的X通道数据的融合参数Xalpha(i，j)：
其中，Xcontrast(i,j)为像素点(i，j)的局部对比度值，Xmax为像素点(i，j)X通道数据的最大值，Xmin为像素点(i，j)X通道数据的最小值，sigma1及sigma2为预设的调节参数，X为R通道或B通道。
34.如权利要求32所述的数字图像处理装置，其特征在于，所述第四融合单元，适于采用如下方法得到像素点(h,k)去除伪彩色后X通道数据的数据值Xfcsoutput(h,k)：
Xfcsoutput(i,j)＝(1-Xalpha(i,j))*XTI(i,j)+Xalpha(i,j)*Xfcs(i,j)
其中，XTI(i,j)为像素点(i，j)X通道的瞬态响应滤波结果，Xfcs(i,j)为像素点(i，j)X通道数据的伪彩色滤波结果，Xfcs(i,j)＝KX(i，j)+GIN(i，j)，GIN(i，j)为像素点(i，j)的G通道数据；Xalpha(i，j)为像素点(i，j)相应通道数据的融合参数Xalpha(i，j)；KX(i，j)为像素点(i，j)色差信号KX值的伪彩色滤波结果。
CN201710177711.XA 2017-03-23 2017-03-23 数字图像处理方法及装置 Active CN108632583B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710177711.XA CN108632583B (zh) 2017-03-23 2017-03-23 数字图像处理方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710177711.XA CN108632583B (zh) 2017-03-23 2017-03-23 数字图像处理方法及装置

Publications (2)

Publication Number Publication Date
CN108632583A CN108632583A (zh) 2018-10-09
CN108632583B true CN108632583B (zh) 2019-10-22

Family

ID=63707414

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710177711.XA Active CN108632583B (zh) 2017-03-23 2017-03-23 数字图像处理方法及装置

Country Status (1)

Country Link
CN (1) CN108632583B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111325700B (zh) * 2020-02-26 2023-07-21 无锡久仁健康云科技有限公司 基于彩色图像的多维度融合方法和系统

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1937781A (zh) * 2005-09-20 2007-03-28 恩益禧电子股份有限公司 图像处理设备和图像处理方法
CN101483784A (zh) * 2008-01-08 2009-07-15 华晶科技股份有限公司 数字影像的伪色抑制方法
CN104093011A (zh) * 2014-07-22 2014-10-08 浙江宇视科技有限公司 一种伪彩色抑制方法和装置
CN105809644A (zh) * 2016-03-15 2016-07-27 深圳英飞拓科技股份有限公司 图像边缘伪色抑制方法及装置

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP3862613B2 (ja) * 2002-06-05 2006-12-27 キヤノン株式会社 画像処理装置及び画像処理方法並びにコンピュータプログラム

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1937781A (zh) * 2005-09-20 2007-03-28 恩益禧电子股份有限公司 图像处理设备和图像处理方法
CN101483784A (zh) * 2008-01-08 2009-07-15 华晶科技股份有限公司 数字影像的伪色抑制方法
CN104093011A (zh) * 2014-07-22 2014-10-08 浙江宇视科技有限公司 一种伪彩色抑制方法和装置
CN105809644A (zh) * 2016-03-15 2016-07-27 深圳英飞拓科技股份有限公司 图像边缘伪色抑制方法及装置

Also Published As

Publication number Publication date
CN108632583A (zh) 2018-10-09

Similar Documents

Publication Publication Date Title
CN106375740B (zh) 生成rgb图像的方法、装置和系统
CN106780311B (zh) 一种结合皮肤粗糙度的快速人脸图像美化方法
US9495582B2 (en) Digital makeup
CN101072289B (zh) 一种图像特效的自动生成方法及装置
CN104899905B (zh) 人脸图像处理方法及装置
CN106127709A (zh) 一种低照度彩色眼底图像判断方法与增强方法
WO2016045242A1 (zh) 一种图像放大方法、图像放大装置及显示设备
CN105913373B (zh) 图像处理方法和装置
CN109040720B (zh) 一种生成rgb图像的方法及装置
JP2008113222A (ja) 画像処理装置、撮影装置、および、これらにおける画像処理方法ならびに当該方法をコンピュータに実行させるプログラム
TW201208349A (en) Method and apparatus for converting 2D image into 3D image
CN104318535B (zh) 图像去雾的方法、装置及移动终端
CN103714520A (zh) 基于fpga实现数字视频图像增强系统及其方法
CN103957396A (zh) 利用智能设备进行舌诊时的图像处理方法、装置及设备
CN110246087A (zh) 参考多通道的多分辨率去除图像色度噪声的系统和方法
CN104392425B (zh) 一种基于人脸的自动调整对比度的图像增强方法
CN108632583B (zh) 数字图像处理方法及装置
CN105763768A (zh) 一种图像处理方法、装置和系统
Peng et al. Single image restoration using scene ambient light differential
CN106097261A (zh) 图像处理方法和装置
CN105389835B (zh) 一种图像处理方法、装置及终端
CN106446859B (zh) 利用手机前置摄像头自动识别人眼中黑点和血丝的方法
CN105654435A (zh) 一种人脸皮肤柔化美白方法
CN108156397A (zh) 一种处理监控画面的方法与设备
KR101513931B1 (ko) 구도의 자동보정 방법 및 이러한 구도의 자동보정 기능이 탑재된 영상 장치

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant