CN109993706A - 数字图像处理方法及装置、计算机可读存储介质 - Google Patents
数字图像处理方法及装置、计算机可读存储介质 Download PDFInfo
- Publication number
- CN109993706A CN109993706A CN201810001923.7A CN201810001923A CN109993706A CN 109993706 A CN109993706 A CN 109993706A CN 201810001923 A CN201810001923 A CN 201810001923A CN 109993706 A CN109993706 A CN 109993706A
- Authority
- CN
- China
- Prior art keywords
- color difference
- difference channel
- pixel point
- nth pixel
- image processing
- 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.)
- Granted
Links
- 238000003672 processing method Methods 0.000 title claims abstract description 30
- 238000003860 storage Methods 0.000 title claims abstract description 14
- 238000001914 filtration Methods 0.000 claims abstract description 170
- 230000001052 transient effect Effects 0.000 claims abstract description 168
- 230000004044 response Effects 0.000 claims abstract description 167
- 239000003086 colorant Substances 0.000 claims abstract description 54
- 230000004927 fusion Effects 0.000 claims description 32
- 238000000034 method Methods 0.000 claims description 14
- 235000013399 edible fruits Nutrition 0.000 claims description 12
- 239000000284 extract Substances 0.000 claims description 9
- 102100037651 AP-2 complex subunit sigma Human genes 0.000 claims description 7
- 101000806914 Homo sapiens AP-2 complex subunit sigma Proteins 0.000 claims description 7
- 241000208340 Araliaceae Species 0.000 claims description 3
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 claims description 3
- 235000003140 Panax quinquefolius Nutrition 0.000 claims description 3
- 235000008434 ginseng Nutrition 0.000 claims description 3
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 3
- 230000002708 enhancing effect Effects 0.000 claims 2
- 238000004364 calculation method Methods 0.000 description 6
- 230000007423 decrease Effects 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 230000003287 optical effect Effects 0.000 description 3
- 230000014759 maintenance of location Effects 0.000 description 2
- 238000012876 topography Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 210000003296 saliva Anatomy 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/77—Retouching; Inpainting; Scratch removal
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10024—Color image
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20024—Filtering details
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20212—Image combination
- G06T2207/20221—Image fusion; Image merging
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Image Processing (AREA)
- Color Image Communication Systems (AREA)
Abstract
一种数字图像处理方法及装置、计算机可读存储介质，所述方法包括：对待处理的YUV数字图像进行瞬态响应滤波处理；对经过瞬态响应滤波处理后的YUV数字图像在kr kb色彩空间进行伪彩色滤波；对所述经过瞬态响应滤波处理后的YUV数字图像及经过伪彩色滤波后的YUV数字图像进行融合，输出所述待处理的YUV数字图像对应的去除伪彩色的YUV数字图像。上述方案可以在不损失图像色度域的信息的情况下，有效去除数字图像中的伪彩色。
Description
技术领域
本发明涉及图像处理技术领域，尤其涉及一种数字图像处理方法及装置、计算机可读存储介质。
背景技术
为了将图像数字化，出现了由光学系统及图像传感器组成的数字图像处理系统。
在数字图像处理系统中，由于光学系统对于不同波长的光折射率不同以及图像传感器在高光下容易溢出等原因，经数字图像处理系统处理后的数字图像经常出现伪彩色，严重影响了数字图像的品质。
为了去除数字图像中的伪彩色，目前通常采用的做法是：将数字图像中所有接近于灰色的色彩抑制为灰色。由于人眼对于接近于灰色的伪彩色会较敏感，由此可以使得数字图像中用户肉眼可见的伪彩色被显著地减少。
然而，在上述伪彩色去除方法中，伪彩色的去除是以牺牲图像色度域的信息为代价的。若要不损失图像色度域的信息，数字图像中的伪彩色就会有所残留，数字图像的品质就会下降。
发明内容
本发明实施例解决的是如何在不损失图像色度域的信息的情况下，有效去除数字图像中的伪彩色。
为解决上述技术问题，本发明实施例提供一种数字图像处理方法，包括：对待处理的YUV数字图像进行瞬态响应滤波处理；对经过瞬态响应滤波处理后的YUV数字图像在krkb色彩空间进行伪彩色滤波；对所述经过瞬态响应滤波处理后的YUV数字图像及经过伪彩色滤波后的YUV数字图像进行融合，输出所述待处理的YUV数字图像对应的去除伪彩色的YUV数字图像。
可选的，所述对待处理的YUV数字图像进行瞬态响应滤波，包括：采用如下步骤对所述待处理的YUV数字图像中的第n个像素点进行瞬态响应滤波处理：以所述第n个像素点为中心的预设第一大小的图像块作为第一图像处理窗口，并根据所述第一图像处理窗口，分别计算所述第一图像处理窗口相应方向上的色差通道X的最大值和最小值；所述色差通道X包括色差通道U和色差通道V；根据所述第一图像处理窗口相应方向上的色差通道X的最大值和最小值，计算在相应方向上，所述第n个像素点的色差通道X对应的瞬态响应滤波结果；对所述第n个像素点的色差通道X在水平方向上的瞬态响应滤波结果与垂直方向上的瞬态响应滤波结果进行融合，得到所述第n个像素点的色差通道X对应的最终瞬态响应滤波结果；所述相应方向包括水平方向和垂直方向。
可选的，所述计算所述第一图像处理窗口相应方向上色差通道X的最大值和最小值，包括：计算相应方向上，位于所述第n个像素点两侧的像素点的色差通道X对应的最大值和最小值；将所述第n个像素点第一侧的色差通道X的最大值与第二侧的色差通道X的最小值，与第二侧的色差通道X的最大值与第一侧的色差通道X的最小值中差值最大的一组所对应的最大值及最小值，作为所述第一图像处理窗口相应方向上色差通道X的最大值和最小值；所述相应方向包括水平方向和垂直方向。
可选的，所述计算在相应方向上，所述第n个像素点的色差通道X对应的瞬态响应滤波结果，包括：计算在相应方向上，所述第n个像素点的色差通道X对应的二阶梯度；根据所述二阶梯度，计算增强因子；根据所述增强因子，计算相应方向上的色差通道X对应的瞬态响应滤波结果；所述相应方向包括水平方向和垂直方向。
可选的，所述计算在相应方向上，所述第n个像素点的色差通道X对应的二阶梯度，包括：在相应方向上，根据所述第n个像素点的色差通道X，以及在相应方向上与所述第n个像素点相邻两侧的像素点的色差通道X，计算所述第n个像素点的色差通道X对应的二阶梯度。
可选的，采用如下公式计算所述增强因子：其中，g(X_dev)为所述增强因子，c为常数，abs(X_dev)为取X_dev的绝对值，X_dev为第n个像素点的色差通道X对应的二阶梯度。
可选的，所述计算在相应方向上，所述第n个像素点的色差通道X对应的瞬态响应滤波结果，包括：当X_dev＞0时，所述第n个像素点的色差通道X对应的瞬态响应滤波结果为：X-g(X_dev)*(X-XMIN)；当X_dev＜0时，所述第n个像素点的色差通道X对应的瞬态响应滤波结果为：X+g(X_dev)*(XMAX-X)；当X_dev＝0时，所述第n个像素点的色差通道X对应的瞬态响应滤波结果为X；其中，X为第n个像素点的色差通道X对应的值，XMIN为第n个像素点的色差通道X在相应方向上的最小值，XMAX为第n个像素点的色差通道X在相应方向上的最大值。
可选的，所述对所述第n个像素点的色差通道X在水平方向上的瞬态响应滤波结果与垂直方向上的瞬态响应滤波结果进行融合，包括：选取所述色差通道X在水平方向上的瞬态响应滤波结果与垂直方向上的瞬态响应滤波结果中的最小值，作为所述色差通道X对应的最终瞬态响应滤波结果。
可选的，所述对经过瞬态响应滤波处理后的YUV数字图像在kr kb色彩空间进行伪彩色滤波，包括：对所述第n个像素点对应的色差通道X进行裁剪处理，包括：从所述第一图像处理窗口中的相应方向上所有像素点的色差通道X中，选择绝对值最小的色差通道X，作为相应方向上裁剪处理后的第n个像素点对应的色差通道X；将所述待处理的YUV数字图像转换成对应的RGB图像，并分别计算第一图像处理窗口中相应方向上的所有像素点的R通道、G通道、D通道以及Y通道对应的梯度；根据所述R通道、G通道、D通道以及Y通道对应的梯度，计算第一图像处理窗口中所有像素点的色差通道X对应的权重系数；根据所述色差通道X对应的权重系数以及所述相应方向上裁剪处理后的第n个像素点对应的色差通道X，计算相应方向上的伪彩色滤波结果；根据所述相应方向上的伪彩色滤波结果，确定所述第n个像素点对应的伪彩色滤波结果。
可选的，所述根据所述相应方向上的伪彩色滤波结果，确定所述第n个像素点对应的伪彩色滤波结果，包括：将水平方向上的伪彩色滤波结果与垂直方向上的伪彩色滤波结果进行比较，从中选择最小值作为所述第n个像素点对应的伪彩色滤波结果。
可选的，所述确定第n个像素点对应的伪彩色滤波结果，还包括：以所述第n个像素点为中心点，提取出大小为3*3的第二图像处理窗口；以所述第n个像素点为中心点，提取出大小为5*5的第三图像处理窗口；分别计算所述第二图像处理窗口的色差通道X对应的第一中值滤波结果，以及所述第三图像处理窗口的色差通道X对应的第二中值滤波结果；将所述水平方向上的伪彩色滤波结果、所述垂直方向上的伪彩色滤波结果、所述第一中值滤波结果、所述第二中值滤波结果进行比较，从中选择中间值作为所述第n个像素点对应的伪彩色滤波结果。
可选的，所述对所述经过瞬态响应滤波处理后的YUV数字图像及经过伪彩色滤波后的YUV数字图像进行融合，包括：计算所述待处理的YUV数字图像中第n个像素点的色差通道X的局部对比度值；基于所述待处理的YUV数字图像中第n个像素点的色差通道X的局部对比度值，计算融合参数；基于所述融合参数，对第n个像素点的色差通道X对应的瞬态响应滤波结果及伪彩色滤波结果进行融合，作为第n个像素点去除伪彩色后的色差通道X。
可选的，所述计算所述待处理的YUV数字图像中第n个像素点的色差通道X的局部对比度值，包括：计算水平方向上第n个像素点的色差通道X对应的局部对比度值，以及垂直方向上第n个像素点的色差通道X对应的局部对比度值；对所述水平方向上第n个像素点的色差通道X对应的局部对比度值，与所述垂直方向上第n个像素点的色差通道X对应的局部对比度值进行比较，从中选择最大值作为所述第n个像素点的色差通道X的局部对比度值。
可选的，所述计算水平方向上第n个像素点的色差通道X对应的局部对比度值，包括：根据所述第一图像处理窗口水平方向上的色差通道X的最大值和最小值，计算所述水平方向上第n个像素点的色差通道X对应的局部对比度值；所述计算垂直方向上第n个像素点的色差通道X对应的局部对比度值，包括：根据所述第一图像处理窗口垂直方向上的色差通道X的最大值和最小值，计算所述垂直方向上第n个像素点的色差通道X对应的局部对比度值。
可选的，所述基于所述待处理的YUV数字图像中第n个像素点的色差通道X的局部对比度值，计算融合参数，包括：采用如下公式计算第n个像素点的色差通道X的融合参数Xalpha(i，j)：
其中，Xcontrast(i,j)为第n个像素点的局部对比度值，Xmax为第n个像素点的X通道数据的最大值，Xmin为第n个像素点的色差通道X的最小值，sigma1及sigma2为预设的调节参数。
可选的，所述基于所述融合参数，对第n个像素点的色差通道X对应的瞬态响应滤波结果及伪彩色滤波结果进行融合，包括：采用如下方法对第n个像素点的色差通道X对应的瞬态响应滤波结果及伪彩色滤波结果进行融合：Xfcsoutput(i,j)＝(1-Xalpha(i,j))*XTI(i,j)+Xalpha(i,j)*Xfcs(i,j)；其中，Xfcsoutput(i,j)为第n个像素点对应的去除伪彩色后相应通道的数据值，XTI(i,j)为第n个像素点的色差通道X的瞬态响应滤波结果，Xfcs(i,j)为第n个像素点的色差通道X的伪彩色滤波结果。
可选的，在对经过瞬态响应滤波处理后的YUV数字图像在kr kb色彩空间进行伪彩色滤波之前，还包括：确定经过滤波处理后的YUV数字图像处于预设的待处理色彩空间内。
本发明实施例还提供了一种数字图像处理装置，包括：瞬态响应滤波处理单元，用于对待处理的YUV数字图像进行瞬态响应滤波处理；伪彩色滤波单元，用于对经过瞬态响应滤波处理后的YUV数字图像在kr kb色彩空间进行伪彩色滤波；融合单元，用于对所述经过瞬态响应滤波处理后的YUV数字图像及经过伪彩色滤波后的YUV数字图像进行融合；输出单元，用于输出所述待处理的YUV数字图像对应的去除伪彩色的YUV数字图像。
可选的，所述瞬态响应滤波处理单元，包括：第一计算子单元，用于以第n个像素点为中心的预设第一大小的图像块作为第一图像处理窗口，并根据所述第一图像处理窗口，分别计算所述第一图像处理窗口相应方向上的色差通道X的最大值和最小值；所述色差通道X包括色差通道U和色差通道V；第二计算子单元，用于根据所述第一图像处理窗口相应方向上的色差通道X的最大值和最小值，计算在相应方向上，所述第n个像素点的色差通道X对应的瞬态响应滤波结果；瞬态响应滤波融合子单元，用于对所述第n个像素点的色差通道X在水平方向上的瞬态响应滤波结果与垂直方向上的瞬态响应滤波结果进行融合，得到所述第n个像素点的色差通道X对应的最终瞬态响应滤波结果；所述相应方向包括水平方向和垂直方向。
可选的，所述第一计算子单元，用于计算相应方向上，位于所述第n个像素点两侧的像素点的色差通道X对应的最大值和最小值；将所述第n个像素点第一侧的色差通道X的最大值与第二侧的色差通道X的最小值，与第二侧的色差通道X的最大值与第一侧的色差通道X的最小值中差值最大的一组所对应的最大值及最小值，作为所述第一图像处理窗口相应方向上色差通道X的最大值和最小值；所述相应方向包括水平方向和垂直方向。
可选的，所述第二计算子单元，用于计算在相应方向上，所述第n个像素点的色差通道X对应的二阶梯度；根据所述二阶梯度，计算增强因子；根据所述增强因子，计算相应方向上的色差通道X对应的瞬态响应滤波结果；所述相应方向包括水平方向和垂直方向。
可选的，所述第二计算子单元，用于在相应方向上，根据所述第n个像素点的色差通道X，以及在相应方向上与所述第n个像素点相邻两侧的像素点的色差通道X，计算所述第n个像素点的色差通道X对应的二阶梯度。
可选的，采用如下公式计算所述增强因子：其中，g(X_dev)为所述增强因子，c为常数，abs(X_dev)为取X_dev的绝对值，X_dev为第n个像素点的色差通道X对应的二阶梯度。
可选的，所述第二计算子单元，用于当X_dev＞0时，确定所述第n个像素点的色差通道X对应的瞬态响应滤波结果为：X-g(X_dev)*(X-XMIN)；当X_dev＜0时，确定所述第n个像素点的色差通道X对应的瞬态响应滤波结果为：X+g(X_dev)*(XMAX-X)；当X_dev＝0时，确定所述第n个像素点的色差通道X对应的瞬态响应滤波结果为X；其中，X为第n个像素点的色差通道X对应的值，XMIN为第n个像素点的色差通道X在相应方向上的最小值，XMAX为第n个像素点的色差通道X在相应方向上的最大值。
可选的，所述瞬态响应滤波融合子单元，用于选取所述色差通道X在水平方向上的瞬态响应滤波结果与垂直方向上的瞬态响应滤波结果中的最小值，作为所述色差通道X对应的最终瞬态响应滤波结果。
可选的，所述伪彩色滤波单元，包括：裁剪子单元，用于对所述第n个像素点对应的色差通道X进行裁剪处理，包括：从所述第一图像处理窗口中的相应方向上所有像素点的色差通道X中，选择绝对值最小的色差通道X，作为相应方向上裁剪处理后的第n个像素点对应的色差通道X；梯度计算子单元，用于将所述待处理的YUV数字图像转换成对应的RGB图像，并分别计算第一图像处理窗口中相应方向上的所有像素点的R通道、G通道、D通道以及Y通道对应的梯度；权重系数计算子单元，用于根据所述R通道、G通道、D通道以及Y通道对应的梯度，计算第一图像处理窗口中所有像素点的色差通道X对应的权重系数；伪彩色滤波结果子单元，用于根据所述色差通道X对应的权重系数以及所述相应方向上裁剪处理后的第n个像素点对应的色差通道X，计算相应方向上的伪彩色滤波结果；确定子单元，用于根据所述相应方向上的伪彩色滤波结果，确定所述第n个像素点对应的伪彩色滤波结果。
可选的，所述确定子单元，用于将水平方向上的伪彩色滤波结果与垂直方向上的伪彩色滤波结果进行比较，从中选择最小值作为所述第n个像素点对应的伪彩色滤波结果。
可选的，所述确定子单元，用于以所述第n个像素点为中心点，提取出大小为3*3的第二图像处理窗口；以所述第n个像素点为中心点，提取出大小为5*5的第三图像处理窗口；分别计算所述第二图像处理窗口的色差通道X对应的第一中值滤波结果，以及所述第三图像处理窗口的色差通道X对应的第二中值滤波结果；将所述水平方向上的伪彩色滤波结果、所述垂直方向上的伪彩色滤波结果、所述第一中值滤波结果、所述第二中值滤波结果进行比较，从中选择中间值作为所述第n个像素点对应的伪彩色滤波结果。
可选的，所述融合单元，包括：局部对比度值计算子单元，用于计算所述待处理的YUV数字图像中第n个像素点的色差通道X的局部对比度值；融合参数计算子单元，用于基于所述待处理的YUV数字图像中第n个像素点的色差通道X的局部对比度值，计算融合参数；伪彩色滤波结果融合子单元，用于基于所述融合参数，对第n个像素点的色差通道X对应的瞬态响应滤波结果及伪彩色滤波结果进行融合，作为第n个像素点去除伪彩色后的色差通道X。
可选的，所述局部对比度值计算子单元，用于计算水平方向上第n个像素点的色差通道X对应的局部对比度值，以及垂直方向上第n个像素点的色差通道X对应的局部对比度值；对所述水平方向上第n个像素点的色差通道X对应的局部对比度值，与所述垂直方向上第n个像素点的色差通道X对应的局部对比度值进行比较，从中选择最大值作为所述第n个像素点的色差通道X的局部对比度值。
可选的，所述局部对比度值计算子单元，用于根据所述第一图像处理窗口水平方向上的色差通道X的最大值和最小值，计算所述水平方向上第n个像素点的色差通道X对应的局部对比度值；根据所述第一图像处理窗口垂直方向上的色差通道X的最大值和最小值，计算所述垂直方向上第n个像素点的色差通道X对应的局部对比度值。
可选的，所述融合参数计算子单元，用于采用如下公式计算第n个像素点的色差通道X的融合参数Xalpha(i，j)：其中，Xcontrast(i,j)为第n个像素点的局部对比度值，Xmax为第n个像素点的X通道数据的最大值，Xmin为第n个像素点的色差通道X的最小值，sigma1及sigma2为预设的调节参数。
可选的，所述伪彩色滤波结果融合子单元，用于采用如下方法对第n个像素点的色差通道X对应的瞬态响应滤波结果及伪彩色滤波结果进行融合：Xfcsoutput(i,j)＝(1-Xalpha(i,j))*XTI(i,j)+Xalpha(i,j)*Xfcs(i,j)；其中，Xfcsoutput(i,j)为第n个像素点对应的去除伪彩色后相应通道的数据值，XTI(i,j)为第n个像素点的色差通道X的瞬态响应滤波结果，Xfcs(i,j)为第n个像素点的色差通道X的伪彩色滤波结果。
可选的，所述数字图像处理装置还包括：确定单元，用于确定经过滤波处理后的YUV数字图像是否处于预设的待处理色彩空间内。
本发明实施例还提供了一种计算机可读存储介质，其上存储有计算机指令，所述计算机指令运行至执行上述任一种所述的数字图像处理方法的步骤。
本发明实施例还提供了一种数字图像处理装置，包括存储器和处理器，所述存储器上存储有计算机指令，所述计算机指令运行至执行上述任一种所述的数字图像处理方法的步骤。
与现有技术相比，本发明实施例的技术方案具有以下有益效果：
对待处理的YUV数字图像整体去除伪彩色，相对于仅对数字图像中接近于灰色的色彩区域去除伪彩色，可以避免数字图像中低饱和度的色彩信息的损失，也就可以避免造成图像色度域的信息的损失，故可以提高数字图像的品质。并且，先对待处理的YUV数字图像进行瞬态响应滤波后，再进行伪彩色滤波，可以压缩伪彩色在局部图像上的范围，在需要去除同样宽度的伪彩色的前提下，有效减小计算量。
进一步，在对经过瞬态响应滤波处理后的YUV数字图像在kr kb色彩空间进行伪彩色滤波之前，先确定经过滤波处理后的YUV数字图像处于预设的待处理色彩空间内，可以减少进行伪彩色滤波的图像的面积，故可以进一步降低计算量。
附图说明
图1是本发明实施例中的一种数字图像处理方法的流程图；
图2是本发明实施例中的一种第一图像处理窗口的示意图；
图3是本发明实施例中的一种数字图像处理装置的结构示意图。
具体实施方式
为了去除数字图像中的伪彩色，目前通常采用的做法是：将数字图像中所有接近于灰色的色彩抑制为灰色。由于人眼对于接近于灰色的伪彩色会较敏感，由此可以使得数字图像中用户肉眼可见的伪彩色被显著地减少。然而，伪彩色的去除是以牺牲图像色度域的信息为代价的。若要不损失图像色度域的信息，数字图像中的伪彩色就会有所残留，数字图像的品质就会下降。
在本发明实施例中，对待处理的YUV数字图像整体去除伪彩色，相对于仅对数字图像中接近于灰色的色彩区域去除伪彩色，可以避免数字图像中低饱和度的色彩信息的损失，也就可以避免造成图像色度域的信息的损失，故可以提高数字图像的品质。并且，先对待处理的YUV数字图像进行瞬态响应滤波后，再进行伪彩色滤波，可以压缩伪彩色在局部图像上的范围，在需要去除同样宽度的伪彩色的前提下，有效减小计算量。
为使本发明的上述目的、特征和有益效果能够更为明显易懂，下面结合附图对本发明的具体实施例做详细的说明。
本发明实施例提供了一种数字图像处理方法，参照图1，以下通过具体步骤进行详细说明。
步骤S101，对待处理的YUV数字图像进行瞬态响应滤波处理。
对待处理的YUV数字图像进行瞬态响应滤波处理，可以压缩伪彩色在局部图像上的范围，在需要去除同样宽度的伪彩色前提下，有效降低YUV数字图像处理的计算复杂度。
步骤S102，对经过瞬态响应滤波处理后的YUV数字图像在kr kb色彩空间进行伪彩色滤波。
在具体实施中，可以在kr kb色彩空间对经过瞬态响应滤波后的YUV图像进行伪彩色滤波，在保留图像中低饱和度色彩信息的前提下，有效去除数字图像中的伪彩色，提高YUV数字图像的品质。
步骤S103，对所述经过瞬态响应滤波处理后的YUV数字图像集经过伪彩色滤波后的YUV数字图像进行融合，输出所述待处理的YUV数字图像对应的去除伪彩色的YUV数字图像。
在具体实施中，对经过瞬态响应滤波处理后的YUV数字图像与经过伪彩色滤波后的YUV数字图像进行融合，可以在去除YUV数字图像中的伪彩色的同时，保护YUV数字图像中的色彩均衡。
下面对步骤S101进行详细说明。
在具体实施中，在对待处理的YUV数字图像进行瞬态滤波响应处理时，由于待处理的YUV数字图像的数据量较大，若一次性完成对待处理的YUV数字图像进行瞬态响应滤波处理，则运算量可能较大。
在本发明实施例中，在对待处理的YUV数字图像进行瞬态响应滤波处理时，可以从待处理的YUV数字图像中选取一块图像区域进行瞬态响应滤波处理。在完成所选取的图像区域的瞬态响应滤波处理之后，从待处理的YUV数字图像中选取另一块图像区域进行瞬态响应滤波处理。以此类推，直至完成对待处理的YUV数字图像中的所有图像区域完成瞬态响应滤波处理，从而完成对待处理的YUV数字图像进行瞬态响应滤波处理。
在具体实施中，为便于表述，可以将从待处理YUV数字图像中选取的图像区域称之为第一图像处理窗口。第一图像处理窗口的大小可以根据YUV数字图像的大小进行选择。
参照图2，给出了本发明实施例中的一种第一图像处理窗口的示意图，第一图像处理窗口的大小为：13*7，也即：第一图像处理窗口水平方向上的长度为13个像素点，垂直方向上的长度为7个像素点。图2中，在水平方向上，像素点CUR左侧像素点个数和右侧像素点个数均为6个，且像素点CUR左侧像素点分别为：L5、L4、L3、L2、L1及L0，像素点CUR右侧像素点分别为：R0、R1、R2、R3、R4及R5。在垂直方向上，像素点CUR上方像素点个数和下方像素点个数均为3个，且像素点CUR上方像素点分别为：T2、T1及T0，像素点CUR下方像素点分别为：B0、B1及B2。
可以理解的是，在实际应用中，第一图像处理窗口的大小还可以为其他值，并不仅限于本发明上述实施例中所提供的举例。例如，第一图像处理窗口的大小为：水平长度为9个像素点，垂直方向为5个像素点。在具体实施中，可以根据实际的应用场景来选择第一图像处理窗口的大小，此处不做赘述。
在本发明实施例中，将像素点CUR作为待处理的YUV数字图像的第n个像素点。
在具体实施中，可以根据第一图像处理窗口，分别计算相应方向上色差通道X的最大值和最小值，相应方向包括水平方向和垂直方向。具体而言，计算第一图像处理窗口相应方向上色差通道X的最大值和最小值，包括：计算第一图像处理窗口水平方向上所有像素点的色差通道X的最大值和最小值，以及计算第一图像处理窗口垂直方向上所有像素点的色差通道X的最大值和最小值。
在本发明实施例中，色差通道X包括色差通道U和色差通道V。因此，计算第一图像处理窗口的色差通道X在水平方向和垂直方向上的最大值，包括：计算第一图像处理窗口的水平方向上所有像素点的色差通道U的最大值、在垂直方向上所有像素点的色差通道U的最大值，以及第一图像处理窗口水平方向上所有像素点的色差通道V的最大值、在垂直方向上所有像素点的色差通道V的最大值。
相应地，计算第一图像处理窗口水平方向和垂直方向上的色差通道X最小值，包括：计算第一图像处理窗口水平方向上所有像素点的色差通道U最小值、在垂直方向上所有像素点的色差通道U的最小值，以及第一图像处理窗口水平方向上所有像素点的色差通道V的最小值、在垂直方向上所有像素点的色差通道V的最小值。
在计算第一图像处理窗口水平方向和垂直方向上的色差通道X的最大值和最小值时，可以采用如下步骤进行：计算在相应方向上，位于第n个像素点两侧的像素点的色差通道X对应的最大值和最小值；将第n个像素点第一侧像素点对应的色差通道X的最大值与第二侧像素点对应的色差通道X的最小值，与第二侧像素点对应的色差通道X的最大值与第一侧像素点对应的色差通道X的最小值中，差值最大的一组所对应的最大值和最小值，作为第一图像处理窗口相应方向上的色差通道X的最大值和最小值；相应方向包括水平方向和垂直方向。
具体而言，计算第一图像处理窗口水平方向上的色差通道U的最大值和最小值，包括：计算在水平方向上，位于第n个像素点左侧的像素点的色差通道U对应的最大值XEMAX和最小值XEMIN，以及位于第n个像素点右侧的像素点的色差通道U对应的最大值XWMAX和最小值XWMIN。
当(XEMAX-XWMIN)＞(XWMAX-XEMIN)时，第一图像处理窗口水平方向上的色差通道U的最大值为XEMAX，最小值为XWMIN；当(XEMAX-XWMIN)≤(XWMAX-XEMIN)时，第一图像处理窗口水平方向上的色差通道U的最大值为XWMAX，最小值为XEMIN。
相应地，计算第一图像处理窗口垂直方向上的色差通道U的最大值和最小值，包括：计算在垂直方向上，位于第n个像素点上方的像素点的色差通道U对应的最大值XUMAX和最小值XUMIN，以及位于第n个像素点下方的像素点的色差通道U对应的最大值XDMAX和最小值XDMIN。
当(XUMAX-XDMIN)＞(XDMAX-XUMIN)时，第一图像处理窗口垂直方向上的色差通道U的最大值为XUMAX，最小值为XDMIN；当(XUMAX-XDMIN)≤(XDMAX-XUMIN)时，第一图像处理窗口垂直方向上的色差通道U的最大值为XDMAX，最小值为XUMIN。
在本发明实施例中，在计算第一图像处理窗口水平方向上的色差通道V对应的最大值和最小值时，可以参照上述计算第一图像处理窗口水平方向上的色差通道U的最大值和最小值的计算过程，只需要将其中的色差通道U替换为色差通道V即可。
相应地，在计算第一图像处理窗口垂直方向上的色差通道V对应的最大值和最小值时，可以参照上述计算第一图像处理窗口垂直方向上的色差通道U的最大值和最小值的计算过程，只需要将其中的色差通道U替换为色差通道V即可。
以图2为例，对第一图像处理窗口水平方向上色差通道U的最大值和最小值、垂直方向上色差通道U的最大值和最小值的计算方法进行说明。
设定像素点L5对应的色差通道U为UL5，像素点L4对应的色差通道U为UL4，像素点R0对应的色差通道U为UR0，以此类推，像素点R5对应的色差通道U为UR5。
在第n个像素点右侧的像素点的色差通道U对应的最大值XEMAX为：
XEMAX＝max(UR0，UR1，UR2，UR3，UR4，UR5)；
在第n个像素点右侧的像素点的色差通道U对应的最小值XEMIN为：
XEMIN＝min(UR0，UR1，UR2，UR3，UR4，UR5)；
第n个像素点左侧的像素点的色差通道U对应的最大值XWMAX为：
XWMAX＝max(UL0，UL1，UL2，UL3，UL4，UL5)；
第n个像素点左侧的像素点的色差通道U对应的最小值XWMIN为：
XWMIN＝min(UL0，UL1，UL2，UL3，UL4，UL5)。
将(XEMAX-XWMIN)与(XWMAX-XEMIN)进行比较，当(XUMAX-XDMIN)＞(XDMAX-XUMIN)时，第一图像处理窗口水平方向上色差通道U的最大值为XUMAX，最小值为XDMIN；当(XUMAX-XDMIN)≤(XDMAX-XUMIN)时，第一图像处理窗口水平方向上色差通道U的最大值为XDMAX，最小值为XUMIN。
设定像素点T0对应的色差通道U为UT0，像素点B0对应的色差通道U为UB0，以此类推。
在第n个像素点上方的像素点的色差通道U对应的最大值XUMAX为：
XUMAX＝max(UT0，UT1，UT2)；
在第n个像素点上方的像素点的色差通道U对应的最小值XUMIN为：
XUMIN＝min(UT0，UT1，UT2)；
第n个像素点下方的像素点的色差通道U对应的最大值XDMAX为：
XDMAX＝max(UB0，UB1，UB2)；
第n个像素点下方的像素点的色差通道U对应的最小值XDMIN为：
XDMIN＝min(UB0，UB1，UB2)。
将(XUMAX-XDMIN)与(XDMAX-XUMIN)进行比较，当(XUMAX-XDMIN)＞(XDMAX-XUMIN)时，第一图像处理窗口垂直方向上色差通道U的最大值为XUMAX，最小值为XDMIN；当(XUMAX-XDMIN)≤(XDMAX-XUMIN)时，第一图像处理窗口垂直方向上色差通道U的最大值为XDMAX，最小值为XUMIN。
在具体实施中，在得到第一图像处理窗口水平方向及垂直方向上色差通道X的最大值和最小值之后，可以根据第一图像处理窗口水平方向上的色差通道X的最大值和最小值，计算水平方向上第n个像素点的色差通道X对应的瞬态响应滤波结果；相应地，可以根据第一图像处理窗口垂直方向上的色差通道X的最大值和最小值，计算垂直方向上第n个像素点的色差通道X对应的瞬态响应滤波结果。
在具体实施中，在计算水平方向和垂直方向上，第n个像素点的色差通道X对应的瞬态响应滤波结果时，可以采用如下步骤：计算在相应方向上，第n个像素点的色差通道X对应的二阶梯度；根据色差通道X对应的二阶梯度，计算增强因子；根据增强因子，计算相应方向上的色差通道X对应的瞬态响应滤波结果。
在本发明实施例中，计算得到的瞬态滤波响应结果包括：水平方向上第n个像素点的色差通道U对应的瞬态响应滤波结果、垂直方向上第n个像素点的色差通道U对应的瞬态响应滤波结果、水平方向上第n个像素点的色差通道V对应的瞬态响应滤波结果、垂直方向上第n个像素点的色差通道V对应的瞬态响应滤波结果。
具体而言，计算在相应方向上，第n个像素点的色差通道X对应的二阶梯度，得到的二阶梯度包括：水平方向上第n个像素点的色差通道U对应的二阶梯度、垂直方向上第n个像素点的色差通道U对应的二阶梯度、水平方向上第n个像素点的色差通道V对应的二阶梯度以及垂直方向上第n个像素点的色差通道V对应的二阶梯度。
针对上述的四组二阶梯度，可以分别计算得到一一对应的增强因子。根据一一对应的增强因子，计算一一对应的瞬态响应滤波结果。
在具体实施中，在计算相应方向上，第n个像素点的色差通道X对应的二阶梯度时，可以根据第n个像素点的色差通道X以及在相应方向上与第n个像素点相邻两侧的像素点的色差通道X，计算第n个像素点的色差通道X对应的二阶梯度。
具体而言，在计算水平方向上第n个像素点的色差通道X对应的二阶梯度时，可以根据第n个像素点的色差通道X以及在水平方向上与第n个像素点相邻两侧的像素点的色差通道X计算得到。相应地，在计算垂直方向上第n个像素点的色差通道X对应的二阶梯度时，可以根据第n个像素点的色差通道X以及在垂直方向上与第n个像素点相邻两侧的像素点的色差通道X计算得到。
在具体实施中，可以根据实际的应用需求来选取第n个像素点相邻两侧的像素点的个数。在水平方向上选取的相邻两侧像素点的个数可以多于垂直方向上选取的相邻两侧像素点的个数。
下面结合图2，对第n个像素点的色差通道X对应的二阶梯度的计算过程进行说明，色差通道X为色差通道U或色差通道V。
在水平方向上，选取的第n个像素点相邻两侧的像素点为L1、L0、R0以及R1。
在水平方向上第n个像素点的色差通道X对应的二阶梯度为：
X_dev＝(X_dev_1-X_dev_0)+2*(X_dev_2-X_dev_1)+(X_dev_3-X_dev_2)； (1)
其中，X_dev_0＝XL0-XL1，X_dev_1＝XCUR-XL0，X_dev_2＝XR0-XCUR，X_dev_2＝XR1-XR0，XL0为像素点L0对应的色差通道X，以此类推，XL1为像素点L1对应的色差通道X，XCUR为第n个像素点对应的色差通道X，XR0为像素点R0对应的色差通道X，XR1为像素点R1对应的色差通道X。
将上式(1)中的X替代为U，即可得到水平方向上第n个像素点的色差通道U对应的二阶梯度；相应地，将上式(1)中的X替代为V，即可得到水平方向上第n个像素点的色差通道V对应的二阶梯度。
在垂直方向上，选取的第n个像素点相邻两侧的像素点为T0以及B0。
在垂直方向上选取的第n个像素点的色差通道X对应的二阶梯度为：
X_dev＝(X_dev_1-X_dev_0)； (2)
其中，X_dev_0＝XCUR-XT0，X_dev_1＝XB0-XCUR，XT0为像素点T0对应的色差通道X，XB0为像素点B0对应的色差通道X。
将上式(2)中的X替代为U，即可得到垂直方向上第n个像素点的色差通道U对应的二阶梯度；相应地，将上式(2)中的X替代为V，即可得到垂直方向上第n个像素点的色差通道V对应的二阶梯度。
在具体实施中，在计算得到第n个像素点的色差通道X对应的二阶梯度之后，可以采用下式计算对应的增强因子：
其中，g(X_dev)为增强因子，c为常数，abs(X_dev)为取X_dev的绝对值，X_dev为第n个像素点的色差通道X对应的二阶梯度。
在具体应用中，c可以为可调常数，可以根据实际的应用场景进行设置。在本发明一实施例中，设定c＝0.3。
在具体实施中，在计算得到第一图像处理窗口相应方向上的色差通道X的最大值和最小值，以及相应方向上第n个像素点对应的增强因子之后，可以计算在相应方向上第n个像素点的色差通道X对应的瞬态响应滤波结果。
具体而言，在计算得到第一图像处理窗口水平方向上的色差通道X的最大值和最小值，以及水平方向上第n个像素点对应的增强因子之后，可以计算在水平方向上第n个像素点的色差通道X对应的瞬态响应滤波结果。
相应地，在计算得到第一图像处理窗口垂直方向上的色差通道X的最大值和最小值，以及垂直方向上第n个像素点对应的增强因子之后，可以计算在垂直方向上第n个像素点的色差通道X对应的瞬态响应滤波结果。
在具体实施中，可以根据相应方向上第n个像素点的色差通道X的二阶梯度，来计算相应方向上的第n个像素点的色差通道X对应的瞬态响应滤波结果。
在具体实施中，当X_dev＞0时，第n个像素点的色差通道X对应的瞬态响应滤波结果为：X-g(X_dev)*(X-XMIN)；当X_dev＜0时，第n个像素点的色差通道X对应的瞬态响应滤波结果为：X+g(X_dev)*(XMAX-X)；当X_dev＝0时，第n个像素点的色差通道X对应的瞬态响应滤波结果为X；其中，X为第n个像素点的色差通道X对应的值，XMIN为第n个像素点的色差通道X在相应方向上的最小值，XMAX为第n个像素点的色差通道X在相应方向上的最大值。
换而言之，在水平方向上，XMIN为第n个像素点的色差通道X在水平方向上的最小值，XMAX为第n个像素点的色差通道X在水平方向上的最大值。在垂直方向上，XMIN为第n个像素点的色差通道X在垂直方向上的最小值，XMAX为第n个像素点的色差通道X在垂直方向上的最大值。XMIN以及XMAX的具体计算过程可以参照本发明上述实施例中提供的计算过程，此处不做赘述。
在具体实施中，在得到水平方向上色差通道X的瞬态响应滤波结果以及垂直方向上色差通道X的瞬态响应滤波结果之后，可以对两个方向上色差通道X的瞬态响应滤波结果进行融合。
在本发明实施例中，对两个方向上色差通道X的瞬态响应滤波结果进行融合，可以为：对水平方向上色差通道X的瞬态响应滤波结果与垂直方向上色差通道X的瞬态响应滤波结果进行比较，选择其中的最小值作为色差通道X对应的最终瞬态响应滤波结果。
例如，水平方向上第n个像素点的色差通道U对应的瞬态响应滤波结果为U_HOR，垂直方向上第n个像素点的色差通道U对应的瞬态相应滤波结果为U_VER，且U_HOR＜U_VER，则选择U_HOR作为第n个像素点的色差通道U对应的最终瞬态响应滤波结果。
下面对步骤S102进行详细说明。
在具体实施中，在获取到瞬态响应滤波处理后的YUV数字图像之后，可以对经过瞬态响应滤波处理后的YUV数字图像在kr kb色彩空间进行伪彩色滤波。
在具体实施中，可以先对第n个像素点的色差通道X进行裁剪(CLIP)处理，得到裁剪处理后的第n个像素点对应的色差通道X。
在本发明实施例中，可以获取第一图像处理窗口的相应方向上所有像素点的色差通道X，并从中选取绝对值最小的色差通道X，作为相应方向上裁剪处理后的第n个像素点对应的色差通道X。色差通道X对应的绝对值可以为：色差通道X对应的值与128相减并取绝对值。
以图2为例，在水平方向上，分别计算像素点L5、L4、L3、L2、L1、L0、CUR、R0、R1、R2、R3、R4以及R5对应的色差通道U对应的绝对值，并取其中的最小值作为水平方向上第n个像素点对应的裁剪后的色差通道U。在垂直方向上，分别计算像素点T1、T0、CUR、B0、B1对应的色差通道U对应的绝对值，并取其中的最小值作为垂直方向上第n个像素点对应的裁剪后的色差通道U。
将待处理的YUV数字图像信号转换成对应的RGB数字图像信号。根据转换得到的RGB数字图像信号，分别计算第一图像处理窗口中所有像素点的R通道对应的梯度DR、G通道对应的梯度DG、B通道对应的梯度DB以及Y通道对应的梯度DY。根据梯度DR、DG、DB以及DY，计算第一图像处理窗口中每一个像素点的色差通道X对应的权重系数。
在具体实施中，第一图像处理窗口中的某一个像素点对应的梯度DR可以为：该像素点对应的R通道的通道数据与第n个像素点对应的R通道的通道数据之差的绝对值；第一图像处理窗口中的某一个像素点对应的梯度DG可以为：该像素点对应的G通道的通道数据与第n个像素点对应的G通道的通道数据之差的绝对值；第一图像处理窗口中的某一个像素点对应的梯度DB可以为：该像素点对应的B通道的通道数据与第n个像素点对应的B通道的通道数据之差的绝对值；第一图像处理窗口中的某一个像素点对应的梯度DY可以为：该像素点对应的Y通道的通道数据与第n个像素点对应的Y通道的通道数据之差的绝对值。
在本发明实施例中，可以采用如下公式来计算第l个像素点的R通道对应的梯度DR：
DR(i,j+l)＝abs(R(i,j+l)-R(i,j))；
其中，DR(i,j+l)为第l个像素点的R通道对应的梯度，R(i,j+l)为第l个像素点对应的R通道，R(i,j)为第n个像素点的R通道；在水平方向上，l的取值范围与第一图像处理窗口水平方向上的像素点个数相关；在垂直方向上，l的取值范围与第一图像处理窗口垂直方向上的像素点的个数相关。
结合图2，l的取值范围为-5≤l≤5，-2≤l≤2。当l为负值时，第l个像素点位于第n个像素点的左侧；当l为正值时，第l个像素点位于第n个像素点的右侧。
例如，l＝-1时，第l个像素点为像素点L1。l＝1时，第l个像素点为像素点R1。
可以采用如下公式来计算第l个像素点的G通道对应的梯度DG：
DG(i,j+l)＝abs(G(i,j+l)-G(i,j))；
其中，DG(i,j+l)为第l个像素点的G通道对应的梯度，G(i,j+l)为第l个像素点对应的G通道，G(i,j)为第n个像素点的G通道。
可以采用如下公式来计算第l个像素点的G通道对应的梯度DB：
DB(i,j+l)＝abs(B(i,j+l)-B(i,j))；
其中，DB(i,j+l)为第l个像素点的B通道对应的梯度，B(i,j+l)为第l个像素点对应的B通道，B(i,j)为第n个像素点的B通道。
可以采用如下公式来计算第l个像素点的G通道对应的梯度DY：
DY(i,j+l)＝abs(Y(i,j+l)-Y(i,j))；
其中，DY(i,j+l)为第l个像素点的Y通道对应的梯度，Y(i,j+l)为第l个像素点对应的Y通道，Y(i,j)为第n个像素点的Y通道。
在具体实施中，在计算色差通道X对应的权重系数时，色差通道U对应的权重系数的计算方法与色差通道V对应的权重系数的计算方法不完全一致。
在本发明实施例中，相应方向上第l个像素点的色差通道U对应的权重系数为：
wk(I,j+l)＝1/(DG(i,j+l)+DB(i,j+l)+DY(i,j+l))；
相应方向上第l个像素点的色差通道V对应的权重系数为：
wk(I,j+l)＝1/(DG(i,j+l)+DR(i,j+l)+DY(i,j+l))；
当相应方向为水平方向时，l的取值范围与第一图像处理窗口水平方向上的像素点个数相关；当相应方向为垂直方向时，l的取值范围与第一图像处理窗口垂直方向上的像素点的个数相关。
结合图2，l的取值范围为-5≤l≤5，-2≤l≤2。当l为负值时，第l个像素点位于第n个像素点的左侧；当l为正值时，第l个像素点位于第n个像素点的右侧。
例如，l＝-1时，第l个像素点为像素点L1。l＝1时，第l个像素点为像素点R1。
在得到权重系数之后，可以根据权重系数以及相应方向上裁剪处理后的第n个像素点对应的色差通道X，计算相应方向上的伪彩色滤波结果。
在具体实施中，可以采用下式来计算得到相应方向上第n个像素点对应的伪彩色滤波结果：
其中，KX(i,j)为第n个像素点对应的伪彩色滤波结果，CLIP_XTI(i,j+l)为裁剪处理后的第n个像素点对应的色差通道X。
KX(i,j)包括水平方向上第n个像素点对应的伪彩色滤波结果以及垂直方向上第n个像素点对应的伪彩色滤波结果。在计算水平方向上第n个像素点对应的伪彩色滤波结果时，根据水平方向上各像素点的权重系数计算；在计算垂直方向上第n个像素点对应的伪彩色滤波结果时，根据垂直方向上各像素点的权重系数计算。
在具体实施中，在计算得到水平方向上第n个像素点对应的伪彩色滤波结果，以及垂直方向上第n个像素点对应的伪彩色滤波结果之后，可以将水平方向上的伪彩色滤波结果与垂直方向上的伪彩色滤波结果进行比较，从中选择最小值作为第n个像素点对应的伪彩色滤波结果。
在具体实施中，在执行步骤S102之前，还可以先判断是否需要对经过瞬态响应滤波处理后的YUV数字图像在kr kb色彩空间进行伪彩色滤波处理。
在具体实施中，可以确定第n个像素点是否处于预设的待处理的色彩空间之内。当第n个像素点处于预设的待处理色彩空间之内时，则执行步骤S102；反之，若第n个像素点处于预设的待处理色彩空间之外时，则无需执行步骤S102，也即无需进行伪彩色滤波。
在实际应用中，可以预先定义待处理色彩空间。例如，预先定义待处理色彩空间为存在伪彩色问题的色彩空间。换而言之，可以预先将存在伪彩色问题的色彩空间定义为待处理色彩空间。针对其他的不存在伪彩色问题的色彩空间，则无需对其进行伪彩色滤波，因此可以提高图像处理的速度，降低计算复杂度。
在具体实施中，在确定第n个像素点对应的伪彩色滤波结果时，还可以结合中值滤波结果。
在本发明实施例中，可以第n个像素点为中心点，提取出大小为3*3的第二图像处理窗口以及大小为5*5的第三图像处理窗口。通过中值滤波器，对第二图像处理窗口进行中值滤波处理，计算得到第二图像处理窗口的色差通道X对应的第一中值滤波结果。相应地，通过中值滤波器，对第三图像处理窗口进行中值滤波处理，计算得到第三图像处理窗口的色差通道X对应的第二中值滤波结果。
将水平方向上第n个像素点的伪彩色滤波结果、垂直方向上第n个像素点的伪彩色滤波结果、第一中值滤波结果以及第二中值滤波结果进行比较，从中选择中间值作为第n个像素点对应的伪彩色滤波结果。
下面对步骤S103进行详细说明。
在对经过瞬态响应滤波处理后的YUV数字图像及经过伪彩色滤波后的YUV数字图像进行融合时，可以先计算待处理的YUV数字图像中第n个像素点的色差通道X的局部对比度值；基于待处理的YUV数字图像中第n个像素点的色差通道X的局部对比度值，计算第n个像素点色差通道X的融合参数；基于所述融合参数，对第n个像素点的色差通道X对应的瞬态响应滤波结果及伪彩色滤波结果进行融合，作为第n个像素点去除伪彩色后的色差通道X。
在具体实施中，在计算水平方向上待处理的YUV数字图像中第n个像素点的色差通道X对应的局部对比度值时，可以采用下式进行计算：
Xcontrasthor(i,j)＝max((XEMAX-XWMIN),(XWMAX-XEMIN))；
在计算垂直方向上待处理的YUV数字图像中第n个像素点的色差通道X对应的局部对比度值时，可以采用下式进行计算：
Xcontrastver(i,j)＝max((XUMAX-XDMIN),(XDMAX-XUMIN))；
其中，Xcontrasthor(i,j)为水平方向上待处理的YUV数字图像中第n个像素点的色差通道X对应的局部对比度值，Xcontrastver(i,j)为垂直方向上待处理的YUV数字图像中第n个像素点的色差通道X对应的局部对比度值。
在具体实施中，在计算得到水平方向和垂直方向上待处理的YUV数字图像中第n个像素点的色差通道X对应的局部对比度值之后，可以从中选取最大值作为第n个像素点的局部对比度值，也即：
Xcontrast(i,j)＝max(Xcontrasthor(i,j)，Xcontrastver(i,j))。
根据第n个像素点的局部对比度值计算融合参数，得到第n个像素点对应的融合参数为：
其中，Xcontrast(i,j)为第n个像素点的局部对比度值，Xmax为第n个像素点的X通道数据的最大值，Xmin为第n个像素点的色差通道X的最小值，sigma1及sigma2为预设的调节参数，sigma1及sigma2均可调。
在计算得到融合参数后，即可对第n个像素点的色差通道X对应的瞬态响应滤波结果及伪彩色滤波结果进行融合。
在具体实施中，可以采用下式对第n个像素点的色差通道X对应的瞬态响应滤波结果及伪彩色滤波结果进行融合：
Xfcsoutput(i,j)＝(1-Xalpha(i,j))*XTI(i,j)+Xalpha(i,j)*Xfcs(i,j)；
其中，Xfcsoutput(i,j)为第n个像素点对应的去除伪彩色后相应通道的数据值，XTI(i,j)为第n个像素点的色差通道X的瞬态响应滤波结果，Xfcs(i,j)为第n个像素点的色差通道X的伪彩色滤波结果。
参照图3，给出了本发明实施例中的一种数字图像处理装置30，包括：瞬态响应滤波处理单元301、伪彩色滤波单元302、融合单元303以及输出单元304，其中：
瞬态响应滤波处理单元301，用于对待处理的YUV数字图像进行瞬态响应滤波处理；
伪彩色滤波单元302，用于对经过瞬态响应滤波处理后的YUV数字图像在kr kb色彩空间进行伪彩色滤波；
融合单元303，用于对所述经过瞬态响应滤波处理后的YUV数字图像及经过伪彩色滤波后的YUV数字图像进行融合；
输出单元304，用于输出所述待处理的YUV数字图像对应的去除伪彩色的YUV数字图像。
可选的，所述瞬态响应滤波处理单元301，可以包括：
第一计算子单元，可以用于以第n个像素点为中心的预设第一大小的图像块作为第一图像处理窗口，并根据所述第一图像处理窗口，分别计算所述第一图像处理窗口相应方向上的色差通道X的最大值和最小值；所述色差通道X包括色差通道U和色差通道V；
第二计算子单元，可以用于根据所述第一图像处理窗口相应方向上的色差通道X的最大值和最小值，计算在相应方向上，所述第n个像素点的色差通道X对应的瞬态响应滤波结果；
瞬态响应滤波融合子单元，可以用于对所述第n个像素点的色差通道X在水平方向上的瞬态响应滤波结果与垂直方向上的瞬态响应滤波结果进行融合，得到所述第n个像素点的色差通道X对应的最终瞬态响应滤波结果；所述相应方向包括水平方向和垂直方向。
在具体实施中，所述第一计算子单元，可以用于计算相应方向上，位于所述第n个像素点两侧的像素点的色差通道X对应的最大值和最小值；将所述第n个像素点第一侧的色差通道X的最大值与第二侧的色差通道X的最小值，与第二侧的色差通道X的最大值与第一侧的色差通道X的最小值中差值最大的一组所对应的最大值及最小值，作为所述第一图像处理窗口相应方向上色差通道X的最大值和最小值；所述相应方向包括水平方向和垂直方向。
在具体实施中，所述第二计算子单元，可以用于计算在相应方向上，所述第n个像素点的色差通道X对应的二阶梯度；根据所述二阶梯度，计算增强因子；根据所述增强因子，计算相应方向上的色差通道X对应的瞬态响应滤波结果；所述相应方向包括水平方向和垂直方向。
在具体实施中，所述第二计算子单元，可以用于在相应方向上，根据所述第n个像素点的色差通道X，以及在相应方向上与所述第n个像素点相邻两侧的像素点的色差通道X，计算所述第n个像素点的色差通道X对应的二阶梯度。
在具体实施中，可以采用如下公式计算所述增强因子：
其中，g(X_dev)为所述增强因子，c为常数，abs(X_dev)为取X_dev的绝对值，X_dev为第n个像素点的色差通道X对应的二阶梯度。
在具体实施中，所述第二计算子单元，可以用于当X_dev＞0时，确定所述第n个像素点的色差通道X对应的瞬态响应滤波结果为：X-g(X_dev)*(X-XMIN)；当X_dev＜0时，确定所述第n个像素点的色差通道X对应的瞬态响应滤波结果为：X+g(X_dev)*(XMAX-X)；当X_dev＝0时，确定所述第n个像素点的色差通道X对应的瞬态响应滤波结果为X；其中，X为第n个像素点的色差通道X对应的值，XMIN为第n个像素点的色差通道X在相应方向上的最小值，XMAX为第n个像素点的色差通道X在相应方向上的最大值。
在具体实施中，所述瞬态响应滤波融合子单元，可以用于选取所述色差通道X在水平方向上的瞬态响应滤波结果与垂直方向上的瞬态响应滤波结果中的最小值，作为所述色差通道X对应的最终瞬态响应滤波结果。
在具体实施中，所述伪彩色滤波单元302，包括：
裁剪子单元，可以用于对所述第n个像素点对应的色差通道X进行裁剪处理，包括：从所述第一图像处理窗口中的相应方向上所有像素点的色差通道X中，选择绝对值最小的色差通道X，作为相应方向上裁剪处理后的第n个像素点对应的色差通道X；
梯度计算子单元，可以用于将所述待处理的YUV数字图像转换成对应的RGB图像，并分别计算第一图像处理窗口中相应方向上的所有像素点的R通道、G通道、D通道以及Y通道对应的梯度；
权重系数计算子单元，可以用于根据所述R通道、G通道、D通道以及Y通道对应的梯度，计算第一图像处理窗口中所有像素点的色差通道X对应的权重系数；
伪彩色滤波结果子单元，可以用于根据所述色差通道X对应的权重系数以及所述相应方向上裁剪处理后的第n个像素点对应的色差通道X，计算相应方向上的伪彩色滤波结果；确定子单元，用于根据所述相应方向上的伪彩色滤波结果，确定所述第n个像素点对应的伪彩色滤波结果。
在具体实施中，所述确定子单元，可以用于将水平方向上的伪彩色滤波结果与垂直方向上的伪彩色滤波结果进行比较，从中选择最小值作为所述第n个像素点对应的伪彩色滤波结果。
在具体实施中，所述确定子单元，可以用于以所述第n个像素点为中心点，提取出大小为3*3的第二图像处理窗口；以所述第n个像素点为中心点，提取出大小为5*5的第三图像处理窗口；分别计算所述第二图像处理窗口的色差通道X对应的第一中值滤波结果，以及所述第三图像处理窗口的色差通道X对应的第二中值滤波结果；将所述水平方向上的伪彩色滤波结果、所述垂直方向上的伪彩色滤波结果、所述第一中值滤波结果、所述第二中值滤波结果进行比较，从中选择中间值作为所述第n个像素点对应的伪彩色滤波结果。
在具体实施中，所述融合单元303，包括：
局部对比度值计算子单元，可以用于计算所述待处理的YUV数字图像中第n个像素点的色差通道X的局部对比度值；
融合参数计算子单元，可以用于基于所述待处理的YUV数字图像中第n个像素点的色差通道X的局部对比度值，计算融合参数；
伪彩色滤波结果融合子单元，可以用于基于所述融合参数，对第n个像素点的色差通道X对应的瞬态响应滤波结果及伪彩色滤波结果进行融合，作为第n个像素点去除伪彩色后的色差通道X。
在具体实施中，所述局部对比度值计算子单元，可以用于计算水平方向上第n个像素点的色差通道X对应的局部对比度值，以及垂直方向上第n个像素点的色差通道X对应的局部对比度值；对所述水平方向上第n个像素点的色差通道X对应的局部对比度值，与所述垂直方向上第n个像素点的色差通道X对应的局部对比度值进行比较，从中选择最大值作为所述第n个像素点的色差通道X的局部对比度值。
在具体实施中，所述局部对比度值计算子单元，可以用于根据所述第一图像处理窗口水平方向上的色差通道X的最大值和最小值，计算所述水平方向上第n个像素点的色差通道X对应的局部对比度值；根据所述第一图像处理窗口垂直方向上的色差通道X的最大值和最小值，计算所述垂直方向上第n个像素点的色差通道X对应的局部对比度值。
在具体实施中，所述融合参数计算子单元，可以用于采用如下公式计算第n个像素点的色差通道X的融合参数Xalpha(i，j)：
其中，Xcontrast(i,j)为第n个像素点的局部对比度值，Xmax为第n个像素点的X通道数据的最大值，Xmin为第n个像素点的色差通道X的最小值，sigma1及sigma2为预设的调节参数。
在具体实施中，所述伪彩色滤波结果融合子单元，可以用于采用如下方法对第n个像素点的色差通道X对应的瞬态响应滤波结果及伪彩色滤波结果进行融合：Xfcsoutput(i,j)＝(1-Xalpha(i,j))*XTI(i,j)+Xalpha(i,j)*Xfcs(i,j)；其中，Xfcsoutput(i,j)为第n个像素点对应的去除伪彩色后相应通道的数据值，XTI(i,j)为第n个像素点的色差通道X的瞬态响应滤波结果，Xfcs(i,j)为第n个像素点的色差通道X的伪彩色滤波结果。
在具体实施中，所述数字图像处理装置还可以包括：确定单元(图3中未示出)，用于确定经过滤波处理后的YUV数字图像是否处于预设的待处理色彩空间内。
本发明实施例还提供了一种计算机可读存储介质，计算机可读存储介质为非易失性存储介质或非瞬态存储介质，其上存储有计算机指令，所述计算机指令运行时执行本发明上述任一实施例提供的数字图像处理方法的步骤，此处不做赘述。
本发明实施例还提供了一种数字图像处理装置，包括存储器和处理器，所述存储器上存储有计算机指令，所述计算机指令运行时执行本发明上述任一实施例提供的数字图像处理方法的步骤，此处不做赘述。
本领域普通技术人员可以理解上述实施例的各种方法中的全部或部分步骤是可以通过程序来指示相关的硬件来完成，该程序可以存储于一计算机可读存储介质中，存储介质可以包括：ROM、RAM、磁盘或光盘等。
虽然本发明披露如上，但本发明并非限定于此。任何本领域技术人员，在不脱离本发明的精神和范围内，均可作各种更动与修改，因此本发明的保护范围应当以权利要求所限定的范围为准。
Claims (36)
1.一种数字图像处理方法，其特征在于，包括：
对待处理的YUV数字图像进行瞬态响应滤波处理；
对经过瞬态响应滤波处理后的YUV数字图像在kr kb色彩空间进行伪彩色滤波；
对所述经过瞬态响应滤波处理后的YUV数字图像及经过伪彩色滤波后的YUV数字图像进行融合，输出所述待处理的YUV数字图像对应的去除伪彩色的YUV数字图像。
2.如权利要求1所述的数字图像处理方法，其特征在于，所述对待处理的YUV数字图像进行瞬态响应滤波，包括：
采用如下步骤对所述待处理的YUV数字图像中的第n个像素点进行瞬态响应滤波处理：
以所述第n个像素点为中心的预设第一大小的图像块作为第一图像处理窗口，并根据所述第一图像处理窗口，分别计算所述第一图像处理窗口相应方向上的色差通道X的最大值和最小值；所述色差通道X包括色差通道U和色差通道V；
根据所述第一图像处理窗口相应方向上的色差通道X的最大值和最小值，计算在相应方向上，所述第n个像素点的色差通道X对应的瞬态响应滤波结果；对所述第n个像素点的色差通道X在水平方向上的瞬态响应滤波结果与垂直方向上的瞬态响应滤波结果进行融合，得到所述第n个像素点的色差通道X对应的最终瞬态响应滤波结果；所述相应方向包括水平方向和垂直方向。
3.如权利要求2所述的数字图像处理方法，其特征在于，所述计算所述第一图像处理窗口相应方向上色差通道X的最大值和最小值，包括：
计算相应方向上，位于所述第n个像素点两侧的像素点的色差通道X对应的最大值和最小值；
将所述第n个像素点第一侧的色差通道X的最大值与第二侧的色差通道X的最小值，与第二侧的色差通道X的最大值与第一侧的色差通道X的最小值中差值最大的一组所对应的最大值及最小值，作为所述第一图像处理窗口相应方向上色差通道X的最大值和最小值；所述相应方向包括水平方向和垂直方向。
4.如权利要求2所述的数字图像处理方法，其特征在于，所述计算在相应方向上，所述第n个像素点的色差通道X对应的瞬态响应滤波结果，包括：
计算在相应方向上，所述第n个像素点的色差通道X对应的二阶梯度；
根据所述二阶梯度，计算增强因子；
根据所述增强因子，计算相应方向上的色差通道X对应的瞬态响应滤波结果；所述相应方向包括水平方向和垂直方向。
5.如权利要求4所述的数字图像处理方法，其特征在于，所述计算在相应方向上，所述第n个像素点的色差通道X对应的二阶梯度，包括：
在相应方向上，根据所述第n个像素点的色差通道X，以及在相应方向上与所述第n个像素点相邻两侧的像素点的色差通道X，计算所述第n个像素点的色差通道X对应的二阶梯度。
6.如权利要求5所述的数字图像处理方法，其特征在于，采用如下公式计算所述增强因子：
其中，g(X_dev)为所述增强因子，c为常数，abs(X_dev)为取X_dev的绝对值，X_dev为第n个像素点的色差通道X对应的二阶梯度。
7.如权利要求6所述的数字图像处理方法，其特征在于，所述计算在相应方向上，所述第n个像素点的色差通道X对应的瞬态响应滤波结果，包括：
当X_dev＞0时，所述第n个像素点的色差通道X对应的瞬态响应滤波结果为：X-g(X_dev)*(X-XMIN)；
当X_dev＜0时，所述第n个像素点的色差通道X对应的瞬态响应滤波结果为：X+g(X_dev)*(XMAX-X)；
当X_dev＝0时，所述第n个像素点的色差通道X对应的瞬态响应滤波结果为X；
其中，X为第n个像素点的色差通道X对应的值，XMIN为第n个像素点的色差通道X在相应方向上的最小值，XMAX为第n个像素点的色差通道X在相应方向上的最大值。
8.如权利要求2所述的数字图像处理方法，其特征在于，所述对所述第n个像素点的色差通道X在水平方向上的瞬态响应滤波结果与垂直方向上的瞬态响应滤波结果进行融合，包括：
选取所述色差通道X在水平方向上的瞬态响应滤波结果与垂直方向上的瞬态响应滤波结果中的最小值，作为所述色差通道X对应的最终瞬态响应滤波结果。
9.如权利要求2所述的数字图像处理方法，其特征在于，所述对经过瞬态响应滤波处理后的YUV数字图像在kr kb色彩空间进行伪彩色滤波，包括：
对所述第n个像素点对应的色差通道X进行裁剪处理，包括：从所述第一图像处理窗口中的相应方向上所有像素点的色差通道X中，选择绝对值最小的色差通道X，作为相应方向上裁剪处理后的第n个像素点对应的色差通道X；将所述待处理的YUV数字图像转换成对应的RGB图像，并分别计算第一图像处理窗口中相应方向上的所有像素点的R通道、G通道、D通道以及Y通道对应的梯度；
根据所述R通道、G通道、D通道以及Y通道对应的梯度，计算第一图像处理窗口中所有像素点的色差通道X对应的权重系数；
根据所述色差通道X对应的权重系数以及所述相应方向上裁剪处理后的第n个像素点对应的色差通道X，计算相应方向上的伪彩色滤波结果；
根据所述相应方向上的伪彩色滤波结果，确定所述第n个像素点对应的伪彩色滤波结果。
10.如权利要求9所述的数字图像处理方法，其特征在于，所述根据所述相应方向上的伪彩色滤波结果，确定所述第n个像素点对应的伪彩色滤波结果，包括：
将水平方向上的伪彩色滤波结果与垂直方向上的伪彩色滤波结果进行比较，从中选择最小值作为所述第n个像素点对应的伪彩色滤波结果。
11.如权利要求10所述的数字图像处理方法，其特征在于，所述确定第n个像素点对应的伪彩色滤波结果，还包括：
以所述第n个像素点为中心点，提取出大小为3*3的第二图像处理窗口；
以所述第n个像素点为中心点，提取出大小为5*5的第三图像处理窗口；
分别计算所述第二图像处理窗口的色差通道X对应的第一中值滤波结果，以及所述第三图像处理窗口的色差通道X对应的第二中值滤波结果；
将所述水平方向上的伪彩色滤波结果、所述垂直方向上的伪彩色滤波结果、所述第一中值滤波结果、所述第二中值滤波结果进行比较，从中选择中间值作为所述第n个像素点对应的伪彩色滤波结果。
12.如权利要求2所述的数字图像处理方法，其特征在于，所述对所述经过瞬态响应滤波处理后的YUV数字图像及经过伪彩色滤波后的YUV数字图像进行融合，包括：
计算所述待处理的YUV数字图像中第n个像素点的色差通道X的局部对比度值；
基于所述待处理的YUV数字图像中第n个像素点的色差通道X的局部对比度值，计算融合参数；
基于所述融合参数，对第n个像素点的色差通道X对应的瞬态响应滤波结果及伪彩色滤波结果进行融合，作为第n个像素点去除伪彩色后的色差通道X。
13.如权利要求12所述的数字图像处理方法，其特征在于，所述计算所述待处理的YUV数字图像中第n个像素点的色差通道X的局部对比度值，包括：计算水平方向上第n个像素点的色差通道X对应的局部对比度值，以及垂直方向上第n个像素点的色差通道X对应的局部对比度值；
对所述水平方向上第n个像素点的色差通道X对应的局部对比度值，与所述垂直方向上第n个像素点的色差通道X对应的局部对比度值进行比较，从中选择最大值作为所述第n个像素点的色差通道X的局部对比度值。
14.如权利要求13所述的数字图像处理方法，其特征在于，所述计算水平方向上第n个像素点的色差通道X对应的局部对比度值，包括：
根据所述第一图像处理窗口水平方向上的色差通道X的最大值和最小值，计算所述水平方向上第n个像素点的色差通道X对应的局部对比度值；
所述计算垂直方向上第n个像素点的色差通道X对应的局部对比度值，包括：根据所述第一图像处理窗口垂直方向上的色差通道X的最大值和最小值，计算所述垂直方向上第n个像素点的色差通道X对应的局部对比度值。
15.如权利要求12所述的数字图像处理方法，其特征在于，所述基于所述待处理的YUV数字图像中第n个像素点的色差通道X的局部对比度值，计算融合参数，包括：
采用如下公式计算第n个像素点的色差通道X的融合参数Xalpha(i，j)：
其中，Xcontrast(i,j)为第n个像素点的局部对比度值，Xmax为第n个像素点的X通道数据的最大值，Xmin为第n个像素点的色差通道X的最小值，sigma1及sigma2为预设的调节参数。
16.如权利要求15所述的数字图像处理方法，其特征在于，所述基于所述融合参数，对第n个像素点的色差通道X对应的瞬态响应滤波结果及伪彩色滤波结果进行融合，包括：
采用如下方法对第n个像素点的色差通道X对应的瞬态响应滤波结果及伪彩色滤波结果进行融合：
Xfcsoutput(i,j)＝(1-Xalpha(i,j))*XTI(i,j)+Xalpha(i,j)*Xfcs(i,j)；
其中，Xfcsoutput(i,j)为第n个像素点对应的去除伪彩色后相应通道的数据值，XTI(i,j)为第n个像素点的色差通道X的瞬态响应滤波结果，Xfcs(i,j)为第n个像素点的色差通道X的伪彩色滤波结果。
17.如权利要求1所述的数字图像处理方法，其特征在于，在对经过瞬态响应滤波处理后的YUV数字图像在kr kb色彩空间进行伪彩色滤波之前，还包括：确定经过滤波处理后的YUV数字图像处于预设的待处理色彩空间内。
18.一种数字图像处理装置，其特征在于，包括：
瞬态响应滤波处理单元，用于对待处理的YUV数字图像进行瞬态响应滤波处理；
伪彩色滤波单元，用于对经过瞬态响应滤波处理后的YUV数字图像在kr kb色彩空间进行伪彩色滤波；
融合单元，用于对所述经过瞬态响应滤波处理后的YUV数字图像及经过伪彩色滤波后的YUV数字图像进行融合；
输出单元，用于输出所述待处理的YUV数字图像对应的去除伪彩色的YUV数字图像。
19.如权利要求18所述的数字图像处理装置，其特征在于，所述瞬态响应滤波处理单元，包括：第一计算子单元，用于以第n个像素点为中心的预设第一大小的图像块作为第一图像处理窗口，并根据所述第一图像处理窗口，分别计算所述第一图像处理窗口相应方向上的色差通道X的最大值和最小值；所述色差通道X包括色差通道U和色差通道V；第二计算子单元，用于根据所述第一图像处理窗口相应方向上的色差通道X的最大值和最小值，计算在相应方向上，所述第n个像素点的色差通道X对应的瞬态响应滤波结果；瞬态响应滤波融合子单元，用于对所述第n个像素点的色差通道X在水平方向上的瞬态响应滤波结果与垂直方向上的瞬态响应滤波结果进行融合，得到所述第n个像素点的色差通道X对应的最终瞬态响应滤波结果；所述相应方向包括水平方向和垂直方向。
20.如权利要求19所述的数字图像处理装置，其特征在于，所述第一计算子单元，用于计算相应方向上，位于所述第n个像素点两侧的像素点的色差通道X对应的最大值和最小值；将所述第n个像素点第一侧的色差通道X的最大值与第二侧的色差通道X的最小值，与第二侧的色差通道X的最大值与第一侧的色差通道X的最小值中差值最大的一组所对应的最大值及最小值，作为所述第一图像处理窗口相应方向上色差通道X的最大值和最小值；所述相应方向包括水平方向和垂直方向。
21.如权利要求19所述的数字图像处理装置，其特征在于，所述第二计算子单元，用于计算在相应方向上，所述第n个像素点的色差通道X对应的二阶梯度；根据所述二阶梯度，计算增强因子；根据所述增强因子，计算相应方向上的色差通道X对应的瞬态响应滤波结果；所述相应方向包括水平方向和垂直方向。
22.如权利要求21所述的数字图像处理装置，其特征在于，所述第二计算子单元，用于在相应方向上，根据所述第n个像素点的色差通道X，以及在相应方向上与所述第n个像素点相邻两侧的像素点的色差通道X，计算所述第n个像素点的色差通道X对应的二阶梯度。
23.如权利要求22所述的数字图像处理装置，其特征在于，采用如下公式计算所述增强因子：
其中，g(X_dev)为所述增强因子，c为常数，abs(X_dev)为取X_dev的绝对值，X_dev为第n个像素点的色差通道X对应的二阶梯度。
24.如权利要求23所述的数字图像处理装置，其特征在于，所述第二计算子单元，用于当X_dev＞0时，确定所述第n个像素点的色差通道X对应的瞬态响应滤波结果为：X-g(X_dev)*(X-XMIN)；当X_dev＜0时，确定所述第n个像素点的色差通道X对应的瞬态响应滤波结果为：X+g(X_dev)*(XMAX-X)；当X_dev＝0时，确定所述第n个像素点的色差通道X对应的瞬态响应滤波结果为X；其中，X为第n个像素点的色差通道X对应的值，XMIN为第n个像素点的色差通道X在相应方向上的最小值，XMAX为第n个像素点的色差通道X在相应方向上的最大值。
25.如权利要求19所述的数字图像处理装置，其特征在于，所述瞬态响应滤波融合子单元，用于选取所述色差通道X在水平方向上的瞬态响应滤波结果与垂直方向上的瞬态响应滤波结果中的最小值，作为所述色差通道X对应的最终瞬态响应滤波结果。
26.如权利要求19所述的数字图像处理装置，其特征在于，所述伪彩色滤波单元，包括：裁剪子单元，用于对所述第n个像素点对应的色差通道X进行裁剪处理，包括：从所述第一图像处理窗口中的相应方向上所有像素点的色差通道X中，选择绝对值最小的色差通道X，作为相应方向上裁剪处理后的第n个像素点对应的色差通道X；梯度计算子单元，用于将所述待处理的YUV数字图像转换成对应的RGB图像，并分别计算第一图像处理窗口中相应方向上的所有像素点的R通道、G通道、D通道以及Y通道对应的梯度；权重系数计算子单元，用于根据所述R通道、G通道、D通道以及Y通道对应的梯度，计算第一图像处理窗口中所有像素点的色差通道X对应的权重系数；伪彩色滤波结果子单元，用于根据所述色差通道X对应的权重系数以及所述相应方向上裁剪处理后的第n个像素点对应的色差通道X，计算相应方向上的伪彩色滤波结果；确定子单元，用于根据所述相应方向上的伪彩色滤波结果，确定所述第n个像素点对应的伪彩色滤波结果。
27.如权利要求26所述的数字图像处理装置，其特征在于，所述确定子单元，用于将水平方向上的伪彩色滤波结果与垂直方向上的伪彩色滤波结果进行比较，从中选择最小值作为所述第n个像素点对应的伪彩色滤波结果。
28.如权利要求27所述的数字图像处理装置，其特征在于，所述确定子单元，用于以所述第n个像素点为中心点，提取出大小为3*3的第二图像处理窗口；以所述第n个像素点为中心点，提取出大小为5*5的第三图像处理窗口；分别计算所述第二图像处理窗口的色差通道X对应的第一中值滤波结果，以及所述第三图像处理窗口的色差通道X对应的第二中值滤波结果；将所述水平方向上的伪彩色滤波结果、所述垂直方向上的伪彩色滤波结果、所述第一中值滤波结果、所述第二中值滤波结果进行比较，从中选择中间值作为所述第n个像素点对应的伪彩色滤波结果。
29.如权利要求19所述的数字图像处理装置，其特征在于，所述融合单元，包括：局部对比度值计算子单元，用于计算所述待处理的YUV数字图像中第n个像素点的色差通道X的局部对比度值；融合参数计算子单元，用于基于所述待处理的YUV数字图像中第n个像素点的色差通道X的局部对比度值，计算融合参数；伪彩色滤波结果融合子单元，用于基于所述融合参数，对第n个像素点的色差通道X对应的瞬态响应滤波结果及伪彩色滤波结果进行融合，作为第n个像素点去除伪彩色后的色差通道X。
30.如权利要求29所述的数字图像处理装置，其特征在于，所述局部对比度值计算子单元，用于计算水平方向上第n个像素点的色差通道X对应的局部对比度值，以及垂直方向上第n个像素点的色差通道X对应的局部对比度值；对所述水平方向上第n个像素点的色差通道X对应的局部对比度值，与所述垂直方向上第n个像素点的色差通道X对应的局部对比度值进行比较，从中选择最大值作为所述第n个像素点的色差通道X的局部对比度值。
31.如权利要求30所述的数字图像处理装置，其特征在于，所述局部对比度值计算子单元，用于根据所述第一图像处理窗口水平方向上的色差通道X的最大值和最小值，计算所述水平方向上第n个像素点的色差通道X对应的局部对比度值；根据所述第一图像处理窗口垂直方向上的色差通道X的最大值和最小值，计算所述垂直方向上第n个像素点的色差通道X对应的局部对比度值。
32.如权利要求29所述的数字图像处理装置，其特征在于，所述融合参数计算子单元，用于采用如下公式计算第n个像素点的色差通道X的融合参数Xalpha(i，j)：
其中，Xcontrast(i,j)为第n个像素点的局部对比度值，Xmax为第n个像素点的X通道数据的最大值，Xmin为第n个像素点的色差通道X的最小值，sigma1及sigma2为预设的调节参数。
33.如权利要求32所述的数字图像处理装置，其特征在于，所述伪彩色滤波结果融合子单元，用于采用如下方法对第n个像素点的色差通道X对应的瞬态响应滤波结果及伪彩色滤波结果进行融合：
Xfcsoutput(i,j)＝(1-Xalpha(i,j))*XTI(i,j)+Xalpha(i,j)*Xfcs(i,j)；
其中，Xfcsoutput(i,j)为第n个像素点对应的去除伪彩色后相应通道的数据值，XTI(i,j)为第n个像素点的色差通道X的瞬态响应滤波结果，Xfcs(i,j)为第n个像素点的色差通道X的伪彩色滤波结果。
34.如权利要求18所述的数字图像处理装置，其特征在于，还包括：确定单元，用于确定经过滤波处理后的YUV数字图像是否处于预设的待处理色彩空间内。
35.一种计算机可读存储介质，其上存储有计算机指令，其特征在于，所述计算机指令运行时执行权利要求1～17任一项所述的数字图像处理方法的步骤。
36.一种数字图像处理装置，包括存储器和处理器，所述存储器上存储有计算机指令，其特征在于，所述计算机指令运行时执行权利要求1～17任一项所述的数字图像处理方法的步骤。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810001923.7A CN109993706B (zh) | 2018-01-02 | 2018-01-02 | 数字图像处理方法及装置、计算机可读存储介质 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810001923.7A CN109993706B (zh) | 2018-01-02 | 2018-01-02 | 数字图像处理方法及装置、计算机可读存储介质 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109993706A true CN109993706A (zh) | 2019-07-09 |
| CN109993706B CN109993706B (zh) | 2021-02-26 |
Family
ID=67128420
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810001923.7A Active CN109993706B (zh) | 2018-01-02 | 2018-01-02 | 数字图像处理方法及装置、计算机可读存储介质 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109993706B (zh) |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1933557A (zh) * | 2005-09-16 | 2007-03-21 | 索尼株式会社 | 成像方法和成像装置 |
| CN103020919A (zh) * | 2013-01-09 | 2013-04-03 | 西安电子科技大学 | 基于非局部Lee的极化SAR相干斑噪声抑制方法 |
| US20130112879A1 (en) * | 2011-11-09 | 2013-05-09 | B.E. Meyers & Co., Inc. | Laser Viewing System |
| CN106507066A (zh) * | 2016-11-22 | 2017-03-15 | 深圳艾科创新微电子有限公司 | 抑制图像伪彩色的方法、装置 |
| CN106920216A (zh) * | 2015-12-25 | 2017-07-04 | 展讯通信（上海）有限公司 | 一种消除图像噪声的方法及装置 |
-
2018
- 2018-01-02 CN CN201810001923.7A patent/CN109993706B/zh active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1933557A (zh) * | 2005-09-16 | 2007-03-21 | 索尼株式会社 | 成像方法和成像装置 |
| US20130112879A1 (en) * | 2011-11-09 | 2013-05-09 | B.E. Meyers & Co., Inc. | Laser Viewing System |
| CN103020919A (zh) * | 2013-01-09 | 2013-04-03 | 西安电子科技大学 | 基于非局部Lee的极化SAR相干斑噪声抑制方法 |
| CN106920216A (zh) * | 2015-12-25 | 2017-07-04 | 展讯通信（上海）有限公司 | 一种消除图像噪声的方法及装置 |
| CN106507066A (zh) * | 2016-11-22 | 2017-03-15 | 深圳艾科创新微电子有限公司 | 抑制图像伪彩色的方法、装置 |
Non-Patent Citations (3)
| Title |
|---|
| MOHAMMAD REZA KHOSRAVI ET AL.: "A novel fake color scheme based on depth protection for MR Passive/Optical Sensors", 《2015 2ND INTERNATIONAL CONFERENCE ON KNOWLEDGE-BASED ENGINEERING AND INNOVATION(KBEI)》 * |
| 丁丽媛: "基于稀疏表示的Bayer彩色滤波阵列插值", 《中国优秀硕士学位论文全文数据库 信息科技辑》 * |
| 朱波 等: "改进的Bayer插值算法及其硬件实现", 《光电子· 激光》 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109993706B (zh) | 2021-02-26 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108205804B (zh) | 图像处理方法、装置及电子设备 | |
| CN104252700B (zh) | 一种红外图像的直方图均衡化方法 | |
| CN110276767A (zh) | 图像处理方法和装置、电子设备、计算机可读存储介质 | |
| CN110473185A (zh) | 图像处理方法和装置、电子设备、计算机可读存储介质 | |
| CN104134191B (zh) | 图像去噪方法及其装置 | |
| CN106127709A (zh) | 一种低照度彩色眼底图像判断方法与增强方法 | |
| CN105809644B (zh) | 图像边缘伪色抑制方法及装置 | |
| CN106375740A (zh) | 生成rgb图像的方法、装置和系统 | |
| CN110263641A (zh) | 疲劳检测方法、装置及可读存储介质 | |
| CN104318535B (zh) | 图像去雾的方法、装置及移动终端 | |
| WO2008145072A1 (fr) | Procédé et dispositif de suppression d'un bruit d'image basé sur le filtre médian | |
| CN101867832A (zh) | 阈值确定装置和阈值确定方法 | |
| CN111340732B (zh) | 一种低照度视频图像增强方法及装置 | |
| CN109389560A (zh) | 一种自适应加权滤波图像降噪方法、装置及图像处理设备 | |
| CN110136144A (zh) | 一种图像分割方法、装置及终端设备 | |
| CN110288670A (zh) | 一种ui描边特效的高性能渲染方法 | |
| CN105809633B (zh) | 去除颜色噪声的方法及装置 | |
| CN106027854B (zh) | 一种应用于相机中适于fpga实现的联合滤波降噪方法 | |
| CN105763768A (zh) | 一种图像处理方法、装置和系统 | |
| CN107767356A (zh) | 一种图像处理方法和装置 | |
| CN106657816A (zh) | 一种基于orb算法的图像配准和图像融合并行的多路视频快速拼接算法 | |
| CN110290370A (zh) | 图像处理方法及装置 | |
| CN106709862B (zh) | 一种图像处理方法及装置 | |
| CN101087364A (zh) | 视频增强系统和方法 | |
| CN105427265A (zh) | 一种彩色图像对比度增强方法及系统 |
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 | ||
| CB02 | Change of applicant information |
Address after: 100089 18 / F, block B, Zhizhen building, No.7, Zhichun Road, Haidian District, Beijing Applicant after: Beijing Ziguang zhanrui Communication Technology Co.,Ltd. Address before: 100084, Room 516, building A, Tsinghua Science Park, Beijing, Haidian District Applicant before: BEIJING SPREADTRUM HI-TECH COMMUNICATIONS TECHNOLOGY Co.,Ltd. |
|
| CB02 | Change of applicant information | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |