CN106683063A - 一种图像去噪的方法与装置 - Google Patents

一种图像去噪的方法与装置 Download PDF

Info

Publication number
CN106683063A
CN106683063A CN201710046553.4A CN201710046553A CN106683063A CN 106683063 A CN106683063 A CN 106683063A CN 201710046553 A CN201710046553 A CN 201710046553A CN 106683063 A CN106683063 A CN 106683063A
Authority
CN
China
Prior art keywords
image
region
pixel
noise point
chromatic noise
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
Application number
CN201710046553.4A
Other languages
English (en)
Other versions
CN106683063B (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.)
Xiamen Meitu Technology Co Ltd
Original Assignee
Xiamen Meitu Technology 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 Xiamen Meitu Technology Co Ltd filed Critical Xiamen Meitu Technology Co Ltd
Priority to CN201710046553.4A priority Critical patent/CN106683063B/zh
Publication of CN106683063A publication Critical patent/CN106683063A/zh
Application granted granted Critical
Publication of CN106683063B publication Critical patent/CN106683063B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00Image enhancement or restoration
    • G06T5/70Denoising; Smoothing
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/20Special algorithmic details
    • G06T2207/20004Adaptive image processing
    • G06T2207/20012Locally adaptive
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/20Special algorithmic details
    • G06T2207/20021Dividing image into blocks, subimages or windows

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Image Processing (AREA)

Abstract

本发明公开一种图像去噪的方法与装置,涉及图像处理领域,为解决去噪后图像模糊的程度高的问题。该图像去噪方法,包括:将图像分割为多个图像区域;将多个图像区域中像素值分布差异表征因子满足预定条件的图像区域确定为具有彩色噪声点的图像区域;对具有彩色噪声点的图像区域进行滤波。

Description

一种图像去噪的方法与装置
技术领域
本发明涉及图像处理领域,尤其涉及一种图像去噪的方法与装置。
背景技术
随着电子设备的发展,越来越多的电子设备都具备有拍照摄影功能,但拍摄过程中会引入彩色噪声,使拍摄得到的图像中出现分散的颜色单一的彩色噪声点,影响图像的美观程度。
为了保证图像的美观,现阶段采用对整个图像进行滤波,将彩色噪声去除的方法。但是,图像包括大量的像素点,由于图像中像素点的数量众多,因此图像的信息较为复杂,在对整个图像进行滤波后,图像中不含彩色噪声点的部分会被过度模糊,导致整个图像出现过度模糊的情况,降低了图像的美观程度。
发明内容
本发明提供了一种图像去噪的方法与装置,能够降低图像模糊的程度,从而保证图像的美观程度。
一方面,本发明提供了一种图像去噪方法,包括:将图像分割为多个图像区域;将多个图像区域中像素值分布差异表征因子满足预定条件的图像区域确定为具有彩色噪声点的图像区域;对具有彩色噪声点的图像区域进行滤波。
另一方面,本发明提供了一种图像去噪装置,包括:图像分割单元,被配置为将图像分割为多个图像区域;彩噪确定单元,被配置为将多个图像区域中像素值分布差异表征因子满足预定条件的图像区域确定为具有彩色噪声点的图像区域;滤波单元,被配置为对具有彩色噪声点的图像区域进行滤波。
本发明提供的图像去噪方法与装置,将图像分割为多个图像区域,将这多个图像区域中像素值分布差异表征因子满足预定条件的图像区域确定为具有彩色噪声的图像区域,对具有彩色噪声点的图像区域进行滤波,从而滤去图像中的彩色噪声。像素分布差异表征因子能够表征像素点在图像区域中的像素值分布差异,若图像区域的像素分布差异因子满足预定条件,表明在该图像区域中具有像素值与图像区域中的大部分像素点有较大差异的彩色噪声点,对具有彩色噪声点的图像区域进行滤波,去除具有彩色噪声点的图像区域中的彩色噪声,并不对不具有彩色噪声的图像区域进行滤波,避免了对整个图像进行滤波,图像中不含彩色噪声点的部分被过度模糊的情况,从而避免了整个图像出现过度模糊的情况,提高了图像的美观程度。
附图说明
从下面结合附图对本发明的具体实施方式的描述中可以更好地理解本发明。其中,相同或相似的附图标记表示相同或相似的特征。
图1为本发明实施例提供的图像去噪的方法的流程图之一;
图2为本发明实施例提供的图像去噪的方法的流程图之二;
图3为本发明实施例提供的图像去噪的方法的流程图之三;
图4为本发明实施例提供的图像去噪的方法的流程图之四;
图5为本发明实施例提供的图像去噪的装置的结构示意图之一;
图6为本发明实施例提供的图像去噪的装置的结构示意图之二。
具体实施方式
下面将详细描述本发明的各个方面的特征和示例性实施例。在下面的详细描述中,提出了许多具体细节,以便提供对本发明的全面理解。但是,对于本领域技术人员来说很明显的是,本发明可以在不需要这些具体细节中的一些细节的情况下实施。下面对实施例的描述仅仅是为了通过示出本发明的示例来提供对本发明的更好的理解。本发明决不限于下面所提出的任何具体配置和算法,而是在不脱离本发明的精神的前提下覆盖了元素、部件和算法的任何修改、替换和改进。在附图和下面的描述中,没有示出公知的结构和技术,以便避免对本发明造成不必要的模糊。
本发明提供一种图像去噪方法和装置,应用于图像处理中,在图像的拍摄过程中,受到拍摄环境的影响或相机本身的影响,拍摄得到的图像中会存在彩色噪声点,影响图像的美观程度。本发明提供的图像去噪方法和装置在图像处理过程中,将图像划分为多个图像区域,根据彩色噪声点的特性,在多个图像区域中确定出具有彩色噪声点的区域,对具有彩色噪声点的图像区域进行滤波处理,将具有彩色噪声点的图像区域中的彩色噪声滤除,并不对不具有彩色噪声点的图像区域进行滤波处理,保证图像中不含彩色噪声点的部分的纹理信息清晰,降低整个图像模糊的程度。
图1为本发明实施例提供的图像去噪的方法的流程图之一,如图1所示,图像去噪的方法包括步骤101-步骤103。
在步骤101中,将图像分割为多个图像区域。
其中,将图像分割为多个图像区域,可以将图像均匀分割为多个图像区域,也可以将图像分为多个超像素,每个超像素作为一个图像区域,在这里并不限定图像分割的具体方法。比如,一张分辨率为1280×960的图片,每个图像区域约有80个像素点,则该图片可以分割为15360个图像区域。具体的,可以将图像的RGB通道(三原色色彩模式通道)转化为Lab通道(明度亮度模式通道),再进行超像素分割,将图像分割为多个超像素。
在步骤102中,将多个图像区域中像素值分布差异表征因子满足预定条件的图像区域确定为具有彩色噪声点的图像区域。
其中,图像区域的像素值分布差异表征因子表征图像区域中像素点的像素值的分布差异,也就是说,根据像素值分布差异表征因子的大小,能够得知图像区域中的像素点的颜色分布状况,像素值分布差异表征因子的值越大,表示该图像区域中的各个像素点之间的颜色差异越大;像素值分布差异表征因子的值越小,则表示该图像区域中的各个像素点之间的颜色差异越小。预定条件可以根据用户对图像处理的效果要求设定,若像素值分布差异表征因子满足预定条件,表示该图像区域中具有和周边像素点的颜色差异很大的像素点,与周边像素点的颜色差异很大的像素点往往就是彩色噪声点;若像素值分布差异表征因子不满足预定条件,表示该图像区域中的所有像素点的颜色相近或过渡自然,并不存在彩色噪声点。
在步骤103中,对具有彩色噪声点的图像区域进行滤波。
其中,为了滤除具有彩色噪声点的图像区域中的彩色噪声,一般采用保边滤波(即边缘保持的滤波),比如双边滤波等滤波方法。为了滤波过程能够更加便捷,会将图像区域的颜色模式转换便于处理的模式,比如YUV模式(Y为明亮度,U和V均为色度),在YUV模式中,彩色噪声体现在U通道和V通道上,并不体现在Y通道上,因此只需要对具有彩色噪声点的图像区域的U通道和V通道进行保边滤波,就可以去除彩色噪声。需要说明的是,在进行滤波时,也可以将图像的颜色模式转换为其他模式,在此并不限定本方案中利用的颜色模式。
本发明提供的图像去噪方法,将图像分割为多个图像区域,将这多个图像区域中像素值分布差异表征因子满足预定条件的图像区域确定为具有彩色噪声的图像区域,对具有彩色噪声点的图像区域进行滤波,从而滤去图像中的彩色噪声。像素分布差异表征因子能够表征像素点在图像区域中的像素值分布差异,若图像区域的像素分布差异因子满足预定条件,表明在该图像区域中具有像素值与图像区域中的大部分像素点有较大差异的彩色噪声点,对具有彩色噪声点的图像区域进行滤波,去除具有彩色噪声点的图像区域中的彩色噪声,并不对不具有彩色噪声的图像区域进行滤波,避免了对整个图像进行滤波,图像中不含彩色噪声点的部分被过度模糊的情况,从而避免了整个图像出现过度模糊的情况,提高了图像的美观程度。
图2为本发明实施例提供的图像去噪的方法的流程图之二,为了能够准确的得到图像区域的像素分布差异表征因子,如图2所示,上述步骤102可以具体细化为步骤1021-步骤1023。
在步骤1021中,对于多个图像区域中的任意一个图像区域,构建图像区域的像素值的协方差矩阵。
需要说明的是,需要说明的是,为了便于计算,在步骤1021之前,可以将图像转换为RGB图像模式,也就是说,像素点的像素值可以利用RGB模式(R为红色通道,G为绿色通道,B为蓝色通道)表示,下面以图像区域的颜色模式为RGB模式为例进行说明,但除RGB模式外的适用于本发明的其他颜色模式表示的像素值也在本发明的保护范围内。
具体的,可以根据公式(1),计算构建图像区域的像素值的协方差矩阵:
其中,Φ为多个图像区域中的任意一个图像区域,cov(Φ)为图像区域Φ的协方差矩阵,N为图像区域Φ中的像素点个数,pj为图像区域Φ中的第j个像素点,像素点pj的像素值为(rj,gj,bj),图像区域Φ的像素平均值为矩阵为矩阵的转置矩阵,其中,j和N均为正整数,且1≤j≤N。值得一提的是,现阶段中,像素点的像素值一般利用三维点坐标表示,但利用其他维数的点坐标表示也适用于本发明的方案中,也属于本发明的保护范围内,比如,像素值用四维点坐标表示,像素点pj的像素值为(Aj,Bj,Cj,Dj)。
在步骤1022中,利用图像区域的像素值的协方差矩阵的特征根,计算图像区域的像素值分布差异表征因子。
其中,根据步骤1021构建的图像区域的像素值的协方差矩阵,能够得到图像区域的像素值的协方差矩阵的特征根,根据图像区域的像素值的协方差矩阵的特征根,能够计算得到该图像区域的像素值分布差异表征因子。本发明实施例以像素值表示为三维点坐标的模式(如RGB模式,YUV模式等)下的图像为例进行说明,像素值表示为三维点坐标的模式下图像区域的像素值的协方差矩阵为3×3矩阵,特征根有3个,图像区域的像素值的协方差矩阵的3个特征根分别为第一特征根λ0、第二特征根λ1和第三特征根λ2,λ0≤λ1≤λ2。具体的,分布差异表征因子包括α和λ0,其中,
在步骤1023中,在图像区域的像素值分布差异表征因子满足预定条件时,确定图像区域为具有彩色噪声点的图像区域。
其中,预定条件包括和λ0>a2,当且λ0>a2时,确定该图像区域为具有彩色噪声点的图像区域,a1为第一判断阈值,a2为第二判断阈值,第一判断阈值和第二判断阈值可以根据对图像进行滤波的预期目标设定。比如,第一判断阈值a1为1/6,第二判断阈值a2为255/6,也就是说,当且λ0>255/6时,可以确定该图像区域为具有彩色噪声点的图像区域。当该图像区域为具有彩色噪声点的图像区域时,彩色噪声点与周边像素点的颜色分布差异较大,图像区域中的像素点颜色分布差异越大,则得到的三个特征根的值越大,且越接近。因此,可以通过的范围来确定图像区域中是否含有彩色噪声点。
需要说明的是,由于图像边缘区域的像素点的颜色分布差异也较大,为了避免误判,在步骤101中采用超像素分割的方法,由于超像素分割本身会沿着图像边缘分割,因此能够避免误判。
图3为本发明实施例提供的图像去噪的方法的流程图之三,如图3所示,上述步骤103可以具体细化为步骤1031-步骤1033。
在步骤1031中,确定具有彩色噪声点的图像区域中的像素点的滤波范围。
其中,滤波范围为以具有彩色噪声点的图像区域的像素点为中心,以预设距离L构造的(2L+1)×(2L+1)的范围。滤波范围需要包含彩色噪声点,因此也可以采用其他方式圈定滤波范围,在此并不限定。需要说明的是,预设距离L可以在具有彩色噪声点的图像区域进行椭圆拟合,椭圆拟合对应的短轴的长度为预设距离L,在同一图像区域中的像素点的预设距离L均一致。
在步骤1032中,根据滤波范围内的像素点的像素值、具有彩色噪声点的图像区域的像素平均值、以及具有彩色噪声点的图像区域的像素值的协方差矩阵的特征向量,计算滤波范围内的像素点的调整系数。
其中,现有技术中对图像进行滤波时,依赖滤波半径、空间距离系数和像素差系数,上述三个系数若选取不当,可能会造成滤波的过度模糊,或者滤波不干净。本发明实施例提供的调整系数能够削弱彩色噪声点的影响权重,从而达到较优的滤波效果。
具体的,可以根据公式获取具有彩色噪声点的图像区域中的像素点沿特征向量的投影。其中,dImi为滤波范围内的第m个像素点沿具有彩色噪声点的图像区域的协方差矩阵的第i个特征向量的投影,滤波范围内的第m个像素点的像素值为(rqm,gqm,bqm),具有彩色噪声点的图像区域的像素平均值为vi为具有彩色噪声点的图像区域的协方差矩阵的第i个特征向量,m与i均为正整数。还可以根据公式获取滤波范围内的第m个像素点的调整系数。其中,β(qm)为滤波范围内的第m个像素点的调整系数,|dIm1|、|dIm2|和|dIm3|分别为滤波范围内的第m个像素点沿具有彩色噪声点的图像区域的协方差矩阵的第一、第二、和第三特征向量的投影的绝对值。由于含彩色噪声点的图像区域的第一特征根、第二特征根和第三特征根比较接近,因此,|dIm1|、|dIm2|和|dIm3|均较大。|dIm1|+|dIm2|+|dIm3|的值越大,滤波范围内的第m个像素点属于彩色噪声点的可能性越大。调整系数的大小与|dIm1|+|dIm2|+|dIm3|的大小成反比,当滤波范围内的第m个像素点属于彩色噪声点的可能性较大时,对应的利用滤波范围内的第m个像素点的调整系数,削弱滤波范围内的第m个像素点对滤波的影响力。
在步骤1033中,根据调整系数,对具有彩色噪声点的图像区域进行滤波。
其中,为了便于进行滤波过程,在步骤1033之前,可以将图像转换为YUV图像模式,也就是说,进行滤波的时候,可以将RGB通道转换为YUV通道,由于彩色噪声点均体现在U通道和V通道,因此只需在U通道和V通道中进行滤波即可。当完成滤波后,可以将滤波后的U通道和V通道结合原始的Y通道,再转换为RGB通道。
具体的,在具有彩色噪声点的图像区域中获取接近像素点和最远像素点,接近像素点的像素值与具有彩色噪声点的图像区域的像素平均值的欧式距离最小,最远像素点的像素值与具有彩色噪声点的图像区域的像素平均值的欧式距离最大。根据公式σs=dmax/3和得到空间距离系数σs和像素差系数σr,dmax为接近像素点的位置和最远像素点的位置的欧式距离,与图像区域的像素平均值的欧式距离最大的像素值为(rmax,gmax,bmax),具有彩色噪声点的图像区域的像素平均值为根据公式得到具有彩色噪声点的图像区域中第k个像素点pk滤波后的像素值C1(pk),其中, C0(pk)为具有彩色噪声点的图像区域中第k个像素点pk滤波前的像素值,为具有彩色噪声点的图像区域中第k个像素点pk的滤波范围,β(q)为滤波范围内的像素点q的调整系数,||pk-q||为像素点pk的位置和像素点q的位置的欧氏距离,||I(pk)-I(q)||为像素点pk的像素值和像素点q的像素值的欧式距离。
值得一提的是,上述C0(pk)和C1(pk)的计算是在YUV通道下进行的,而||I(pk)-I(q)||的计算实在RGB通道下进行的。
需要说明的是,在滤波过程中引入了调整系数,利用调整系数来调整滤波过程中的彩色噪声点的权重的影响,削弱滤波过程中彩色噪声点的权重,增加了非彩色噪声点的权重,从而尽可能地将彩色噪声点过滤干净,并避免出现过度模糊的情况。而且,在现有技术中对整个图像进行滤波时采用的空间距离系数σs和像素差系数σr均为固定值,在本发明实施例中,σs是根据接近像素点的位置和最远像素点的位置的欧式距离计算得到的,σr是根据与图像区域的像素平均值的欧式距离最大的像素值为(rmax,gmax,bmax)以及具有彩色噪声点的图像区域的像素平均值计算得到的。也就是说,本发明实施例中的空间距离系数σs和像素差系数σr是根据具有彩色噪声点的图像区域自适应变化的,在滤波过程中,能够避免空间距离系数σs和像素差系数σr固定引起的滤波过度模糊或者滤波不干净的情况。
图4为本发明实施例提供的图像去噪的方法的流程图之四,如图4所示,上述实施例还可以包括步骤104。
在步骤104中,对具有彩色噪声点的图像区域相邻的图像区域进行滤波。
其中,由于具有彩色噪声点的图像区域相邻的图像区域也可能部分含有彩色噪声点,为了尽量将图像中的彩色噪声点过滤干净,因此需要对具有彩色噪声点的图像区域相邻的图像区域也进行滤波。
在本发明实施例中,能够通过对确定的具有彩色噪声点的图像区域相邻的图像区域进行滤波,而进一步将整个图像中的彩色噪声点过滤的更干净。
图5为本发明实施例提供的图像去噪的装置的结构示意图之一,如图5所示,图像去噪的装置200包括图像分割单元201、彩噪确定单元202和滤波单元203。
其中,图像分割单元201,被配置为将图像分割为多个图像区域。
彩噪确定单元202,被配置为将多个图像区域中像素值分布差异表征因子满足预定条件的图像区域确定为具有彩色噪声点的图像区域。
滤波单元205,被配置为对具有彩色噪声点的图像区域进行滤波。
本发明提供的图像去噪的装置200,图像分割单元201将图像分割为多个图像区域,彩噪确定单元202将这多个图像区域中像素值分布差异表征因子满足预定条件的图像区域确定为具有彩色噪声的图像区域,滤波单元203对具有彩色噪声点的图像区域进行滤波,从而滤去图像中的彩色噪声。像素分布差异表征因子能够表征像素点在图像区域中的像素值分布差异,若图像区域的像素分布差异因子满足预定条件,表明在该图像区域中具有像素值与图像区域中的大部分像素点有较大差异的彩色噪声点,对具有彩色噪声点的图像区域进行滤波,去除具有彩色噪声点的图像区域中的彩色噪声,并不对不具有彩色噪声的图像区域进行滤波,避免了对整个图像进行滤波,图像中不含彩色噪声点的部分被过度模糊的情况,从而避免了整个图像出现过度模糊的情况,提高了图像的美观程度。
具体的,上述实施例中的彩噪确定单元202具体被配置为对于多个图像区域中的任意一个图像区域,构建图像区域的像素值的协方差矩阵;利用图像区域的像素值的协方差矩阵的特征根,计算图像区域的像素值分布差异表征因子;以及在图像区域的像素值分布差异表征因子满足预定条件时,确定图像区域为具有彩色噪声点的图像区域。
其中,分布差异表征因子包括α和λ0λ0、λ1和λ2分别为图像区域的像素值的协方差矩阵的第一、第二、和第三特征根,λ0≤λ1≤λ2
更进一步地,彩噪确定单元202具体被配置为在且λ0>a2时,确定图像区域为具有彩色噪声点的图像区域,a1为第一判断阈值,a2为第二判断阈值。
其中,第一判断阈值a1可以为1/6,第二判断阈值a2可以为255/6。
具体的,上述实施例中的滤波单元203具体被配置为确定具有彩色噪声点的图像区域中的像素点的滤波范围;根据滤波范围内的像素点的像素值、具有彩色噪声点的图像区域的像素平均值、以及具有彩色噪声点的图像区域的像素值的协方差矩阵的特征向量,计算滤波范围内的像素点的调整系数;根据调整系数,对具有彩色噪声点的图像区域进行滤波。
其中,滤波范围为以具有彩色噪声点的图像区域的像素点为中心,以预设距离L构造的(2L+1)×(2L+1)的范围。
更进一步地,滤波单元203具体被配置为根据公式获取具有彩色噪声点的图像区域中的像素点沿特征向量的投影,dImi为滤波范围内的第m个像素点沿具有彩色噪声点的图像区域的协方差矩阵的第i个特征向量的投影,滤波范围内的第m个像素点的像素值为(rqm,gqm,bqm),具有彩色噪声点的图像区域的像素平均值为vi为具有彩色噪声点的图像区域的协方差矩阵的第i个特征向量,m与i均为正整数;根据公式获取滤波范围内的第m个像素点的调整系数,β(qm)为滤波范围内的第m个像素点的调整系数,|dIm1|、|dIm2|和|dIm3|分别为滤波范围内的第m个像素点沿具有彩色噪声点的图像区域的协方差矩阵的第一、第二、和第三特征向量的投影的绝对值。
上述滤波单元203具体被配置为在具有彩色噪声点的图像区域中获取接近像素点和最远像素点,接近像素点的像素值与具有彩色噪声点的图像区域的像素平均值的欧式距离最小,最远像素点的像素值与具有彩色噪声点的图像区域的像素平均值的欧式距离最大;根据公式σs=dmax/3和得到空间距离系数σs和像素差系数σr,dmax为接近像素点的位置和最远像素点的位置的欧式距离,与图像区域的像素平均值的欧式距离最大的像素值为(rmax,gmax,bmax),具有彩色噪声点的图像区域的像素平均值为根据公式得到具有彩色噪声点的图像区域中第k个像素点pk滤波后的像素值C1(pk),其中, C0(pk)为具有彩色噪声点的图像区域中第k个像素点pk滤波前的像素值,为具有彩色噪声点的图像区域中第k个像素点pk的滤波范围,β(q)为滤波范围内的像素点q的调整系数,||pk-q||为像素点pk的位置和像素点q的位置的欧氏距离,||I(pk)-I(q)||为像素点pk的像素值和像素点q的像素值的欧式距离。
在本发明实施例中,滤波单元203在滤波过程中引入了调整系数,利用调整系数来调整滤波过程中的彩色噪声点的权重的影响,削弱滤波过程中彩色噪声点的权重,增加了非彩色噪声点的权重,从而尽可能地将彩色噪声点过滤干净,并避免出现过度模糊的情况。而且,在现有技术中对整个图像进行滤波时采用的空间距离系数σs和像素差系数σr均为固定值,在本发明实施例中,σs是根据接近像素点的位置和最远像素点的位置的欧式距离计算得到的,σr是根据与图像区域的像素平均值的欧式距离最大的像素值为(rmax,gmax,bmax)以及具有彩色噪声点的图像区域的像素平均值计算得到的。也就是说,本发明实施例中的空间距离系数σs和像素差系数σr是根据具有彩色噪声点的图像区域自适应变化的,在滤波过程中,能够避免空间距离系数σs和像素差系数σr固定引起的滤波过度模糊或者滤波不干净的情况。
图6为本发明实施例提供的图像去噪的装置的结构示意图之二,如图6所示,在图5所示的图像去噪的装置200的基础上,图像去噪的装置200还包括周边滤波单元204,被配置为对具有彩色噪声点的图像区域相邻的图像区域进行滤波。
在本发明实施例中,周边滤波单元204能够通过对确定的具有彩色噪声点的图像区域相邻的图像区域进行滤波,而进一步将整个图像中的彩色噪声点过滤的更干净。
以上所述的结构示意图中所示的功能单元可以实现为硬件、软件、固件或者它们的组合。当以硬件方式实现时,其可以例如是电子电路、专用集成电路(ASIC)、适当的固件、插件、功能卡等等。当以软件方式实现时,本发明的元素是被用于执行所需任务的程序或者代码段。程序或者代码段可以存储在机器可读介质中,或者通过载波中携带的数据信号在传输介质或者通信链路上传送。“机器可读介质”可以包括能够存储或传输信息的任何介质。机器可读介质的例子包括电子电路、半导体存储器设备、ROM、闪存、可擦除ROM(EROM)、软盘、CD-ROM、光盘、硬盘、光纤介质、射频(RF)链路,等等。代码段可以经由诸如因特网、内联网等的计算机网络被下载。
本发明并不局限于上文所描述并在图中示出的特定步骤和结构。并且,为了简明起见,这里省略对已知方法技术的详细描述。例如,在上述实施例中,描述和示出了若干具体的步骤作为示例。但是,本发明的方法过程并不限于所描述和示出的具体步骤,本领域的技术人员可以在领会本发明的精神之后,做出各种改变、修改和添加,或者改变步骤之间的顺序。

Claims (14)

1.一种图像去噪方法,包括:
将图像分割为多个图像区域;
将所述多个图像区域中像素值分布差异表征因子满足预定条件的图像区域确定为具有彩色噪声点的图像区域;
对所述具有彩色噪声点的图像区域进行滤波。
2.根据权利要求1所述的图像去噪方法,其中,将所述多个图像区域中像素值分布差异表征因子满足预定条件的图像区域确定为具有彩色噪声点的图像区域的处理包括:
对于所述多个图像区域中的任意一个图像区域,构建所述图像区域的像素值的协方差矩阵;
利用所述图像区域的像素值的协方差矩阵的特征根,计算所述图像区域的像素值分布差异表征因子;以及
在所述图像区域的像素值分布差异表征因子满足所述预定条件时,确定所述图像区域为所述具有彩色噪声点的图像区域。
3.根据权利要求2所述的图像去噪方法,所述分布差异表征因子包括α和λ0,其中,其中,λ0、λ1和λ2分别为所述图像区域的像素值的协方差矩阵的第一、第二、和第三特征根,λ0≤λ1≤λ2
在所述图像区域的像素值分布差异表征因子满足所述预定条件时,确定所述图像区域为所述具有彩色噪声点的图像区域的处理包括:
且λ0>a2时,确定所述图像区域为所述具有彩色噪声点的图像区域,a1为所述第一判断阈值,a2为所述第二判断阈值。
4.根据权利要求1所述的图像去噪方法,其中,对所述具有彩色噪声点的图像区域进行滤波的处理包括:
确定所述具有彩色噪声点的图像区域中的像素点的滤波范围;
根据所述滤波范围内的像素点的像素值、所述具有彩色噪声点的图像区域的像素平均值、以及所述具有彩色噪声点的图像区域的像素值的协方差矩阵的特征向量,计算所述滤波范围内的像素点的调整系数;
根据所述调整系数,对所述具有彩色噪声点的图像区域进行滤波。
5.根据权利要求4所述的图像去噪方法,其中,根据所述滤波范围内的像素点的像素值、所述具有彩色噪声点的图像区域的像素平均值、以及所述具有彩色噪声点的图像区域的像素值的协方差矩阵的特征向量,计算所述滤波范围内的像素点的调整系数的处理包括:
根据公式获取所述具有彩色噪声点的图像区域中的像素点沿所述特征向量的投影,dImi为所述滤波范围内的第m个像素点沿所述具有彩色噪声点的图像区域的协方差矩阵的第i个特征向量的投影,所述滤波范围内的第m个像素点的像素值为(rqm,gqm,bqm),所述具有彩色噪声点的图像区域的像素平均值为vi为所述具有彩色噪声点的图像区域的协方差矩阵的第i个特征向量,m与i均为正整数;
根据公式获取所述滤波范围内的第m个像素点的调整系数,β(qm)为所述滤波范围内的第m个像素点的调整系数,|dIm1|、|dIm2|和|dIm3|分别为所述滤波范围内的第m个像素点沿所述具有彩色噪声点的图像区域的协方差矩阵的第一、第二、和第三特征向量的投影的绝对值。
6.根据权利要求4所述的图像去噪方法,其中,根据所述滤波范围内的像素点的调整系数,对所述具有彩色噪声点的图像区域进行滤波的处理包括:
在所述具有彩色噪声点的所述图像区域中获取接近像素点和最远像素点,所述接近像素点的像素值与具有彩色噪声点的所述图像区域的像素平均值的欧式距离最小,所述最远像素点的像素值与具有彩色噪声点的所述图像区域的像素平均值的欧式距离最大;
根据公式σs=dmax/3和得到空间距离系数σs和像素差系数σr,dmax为所述接近像素点的位置和最远像素点的位置的欧式距离,与所述图像区域的像素平均值的欧式距离最大的像素值为(rmax,gmax,bmax),具有彩色噪声点的所述图像区域的像素平均值为
根据公式得到具有彩色噪声点的所述图像区域中第k个像素点pk滤波后的像素值C1(pk),其中,C0(pk)为具有彩色噪声点的所述图像区域中第k个像素点pk滤波前的像素值,为具有彩色噪声点的所述图像区域中第k个像素点pk的滤波范围,β(q)为所述滤波范围内的像素点q的调整系数,||pk-q||为像素点pk的位置和像素点q的位置的欧氏距离,||I(pk)-I(q)||为像素点pk的像素值和像素点q的像素值的欧式距离。
7.根据权利要求1所述的图像去噪的方法,还包括:
对所述具有彩色噪声点的图像区域相邻的图像区域进行滤波。
8.一种图像去噪装置,包括:
图像分割单元,被配置为将图像分割为多个图像区域;
彩噪确定单元,被配置为将所述多个图像区域中像素值分布差异表征因子满足预定条件的图像区域确定为具有彩色噪声点的图像区域;
滤波单元,被配置为对所述具有彩色噪声点的图像区域进行滤波。
9.根据权利要求8所述的图像去噪装置,其中,所述彩噪确定单元具体被配置为:
对于所述多个图像区域中的任意一个图像区域,构建所述图像区域的像素值的协方差矩阵;
利用所述图像区域的像素值的协方差矩阵的特征根,计算所述图像区域的像素值分布差异表征因子;以及
在所述图像区域的像素值分布差异表征因子满足所述预定条件时,确定所述图像区域为所述具有彩色噪声点的图像区域。
10.根据权利要求9所述的图像去噪装置,所述分布差异表征因子包括α和λ0,其中,其中,λ0、λ1和λ2分别为所述图像区域的像素值的协方差矩阵的第一、第二、和第三特征根,λ0≤λ1≤λ2
所述彩噪确定单元具体被配置为:
且λ0>a2时,确定所述图像区域为所述具有彩色噪声点的图像区域,a1为所述第一判断阈值,a2为所述第二判断阈值。
11.根据权利要求8所述的图像去噪装置,其中,所述滤波单元具体被配置为:
确定所述具有彩色噪声点的图像区域中的像素点的滤波范围;
根据所述滤波范围内的像素点的像素值、所述具有彩色噪声点的图像区域的像素平均值、以及所述具有彩色噪声点的图像区域的像素值的协方差矩阵的特征向量,计算所述滤波范围内的像素点的调整系数;
根据所述调整系数,对所述具有彩色噪声点的图像区域进行滤波。
12.根据权利要求11所述的图像去噪装置,其中,所述滤波单元具体被配置为:
根据公式获取所述具有彩色噪声点的图像区域中的像素点沿所述特征向量的投影,dImi为所述滤波范围内的第m个像素点沿所述具有彩色噪声点的图像区域的协方差矩阵的第i个特征向量的投影,所述滤波范围内的第m个像素点的像素值为(rqm,gqm,bqm),所述具有彩色噪声点的图像区域的像素平均值为vi为所述具有彩色噪声点的图像区域的协方差矩阵的第i个特征向量,m与i均为正整数;
根据公式获取所述滤波范围内的第m个像素点的调整系数,β(qm)为所述滤波范围内的第m个像素点的调整系数,|dIm1|、、|dIm2|和|dIm3|分别为所述滤波范围内的第m个像素点沿所述具有彩色噪声点的图像区域的协方差矩阵的第一、第二、和第三特征向量的投影的绝对值。
13.根据权利要求11所述的图像去噪装置,其中,所述滤波单元具体被配置为:
在所述具有彩色噪声点的所述图像区域中获取接近像素点和最远像素点,所述接近像素点的像素值与具有彩色噪声点的所述图像区域的像素平均值的欧式距离最小,所述最远像素点的像素值与具有彩色噪声点的所述图像区域的像素平均值的欧式距离最大;
根据公式σs=dmax/3和得到空间距离系数σs和像素差系数σr,dmax为所述接近像素点的位置和最远像素点的位置的欧式距离,与所述图像区域的像素平均值的欧式距离最大的像素值为(rmax,gmax,bmax),具有彩色噪声点的所述图像区域的像素平均值为
根据公式得到具有彩色噪声点的所述图像区域中第k个像素点pk滤波后的像素值C1(pk),其中,C0(Pk)为具有彩色噪声点的所述图像区域中第k个像素点pk滤波前的像素值,为具有彩色噪声点的所述图像区域中第k个像素点pk的滤波范围,β(q)为所述滤波范围内的像素点q的调整系数,||pk-q||为像素点pk的位置和像素点q的位置的欧氏距离,||I(pk)-I(q)||为像素点pk的像素值和像素点q的像素值的欧式距离。
14.根据权利要求8所述的图像去噪的装置,还包括:
周边滤波单元,被配置为对所述具有彩色噪声点的图像区域相邻的图像区域进行滤波。
CN201710046553.4A 2017-01-22 2017-01-22 一种图像去噪的方法与装置 Active CN106683063B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710046553.4A CN106683063B (zh) 2017-01-22 2017-01-22 一种图像去噪的方法与装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710046553.4A CN106683063B (zh) 2017-01-22 2017-01-22 一种图像去噪的方法与装置

Publications (2)

Publication Number Publication Date
CN106683063A true CN106683063A (zh) 2017-05-17
CN106683063B CN106683063B (zh) 2019-07-05

Family

ID=58859613

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710046553.4A Active CN106683063B (zh) 2017-01-22 2017-01-22 一种图像去噪的方法与装置

Country Status (1)

Country Link
CN (1) CN106683063B (zh)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109410132A (zh) * 2018-09-29 2019-03-01 深圳市牧月科技有限公司 一种基于边缘保护滤波器统一框架的图像处理方法
CN109754376A (zh) * 2018-12-28 2019-05-14 深圳美图创新科技有限公司 图像去噪方法及装置
CN110472466A (zh) * 2019-02-27 2019-11-19 鑫喆喆 基于大数据分析的模式切换系统及方法
CN112419161A (zh) * 2019-08-20 2021-02-26 RealMe重庆移动通信有限公司 图像处理方法及装置、存储介质及电子设备
CN113052767A (zh) * 2019-12-26 2021-06-29 深圳市万普拉斯科技有限公司 滤波矩阵生成方法、装置、终端设备和可读存储介质

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20050018889A1 (en) * 2003-07-09 2005-01-27 Jianying Li Systems and methods for filtering images
US20070201064A1 (en) * 2006-02-28 2007-08-30 Jae-Sung Heo Method and apparatus for removing color noise in image signal
CN101375588B (zh) * 2006-01-27 2011-07-06 微软公司 用于基于区域的图像去噪的方法和数码相机
CN102349090A (zh) * 2009-03-16 2012-02-08 株式会社理光 噪声减少设备、噪声减少方法、噪声减少程序和记录介质
CN104103082A (zh) * 2014-06-06 2014-10-15 华南理工大学 一种基于区域描述和先验知识的图像显著性检测方法
CN105354599A (zh) * 2015-11-25 2016-02-24 安徽创世科技有限公司 一种基于改进的slic超像素分割算法的颜色识别方法
CN105931197A (zh) * 2016-04-13 2016-09-07 石家庄铁道大学 一种基于模糊度理论的图像去噪方法

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20050018889A1 (en) * 2003-07-09 2005-01-27 Jianying Li Systems and methods for filtering images
CN101375588B (zh) * 2006-01-27 2011-07-06 微软公司 用于基于区域的图像去噪的方法和数码相机
US20070201064A1 (en) * 2006-02-28 2007-08-30 Jae-Sung Heo Method and apparatus for removing color noise in image signal
CN102349090A (zh) * 2009-03-16 2012-02-08 株式会社理光 噪声减少设备、噪声减少方法、噪声减少程序和记录介质
CN104103082A (zh) * 2014-06-06 2014-10-15 华南理工大学 一种基于区域描述和先验知识的图像显著性检测方法
CN105354599A (zh) * 2015-11-25 2016-02-24 安徽创世科技有限公司 一种基于改进的slic超像素分割算法的颜色识别方法
CN105931197A (zh) * 2016-04-13 2016-09-07 石家庄铁道大学 一种基于模糊度理论的图像去噪方法

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109410132A (zh) * 2018-09-29 2019-03-01 深圳市牧月科技有限公司 一种基于边缘保护滤波器统一框架的图像处理方法
CN109754376A (zh) * 2018-12-28 2019-05-14 深圳美图创新科技有限公司 图像去噪方法及装置
CN110472466A (zh) * 2019-02-27 2019-11-19 鑫喆喆 基于大数据分析的模式切换系统及方法
CN110472466B (zh) * 2019-02-27 2020-05-29 南京奇趣数字科技有限公司 基于大数据分析的模式切换系统及方法
CN112419161A (zh) * 2019-08-20 2021-02-26 RealMe重庆移动通信有限公司 图像处理方法及装置、存储介质及电子设备
CN112419161B (zh) * 2019-08-20 2022-07-05 RealMe重庆移动通信有限公司 图像处理方法及装置、存储介质及电子设备
CN113052767A (zh) * 2019-12-26 2021-06-29 深圳市万普拉斯科技有限公司 滤波矩阵生成方法、装置、终端设备和可读存储介质

Also Published As

Publication number Publication date
CN106683063B (zh) 2019-07-05

Similar Documents

Publication Publication Date Title
CN106683063B (zh) 一种图像去噪的方法与装置
CN107403421B (zh) 一种图像去雾方法、存储介质及终端设备
CN106375740B (zh) 生成rgb图像的方法、装置和系统
Lee et al. A space-variant luminance map based color image enhancement
CN111429370B (zh) 一种煤矿井下的图像增强方法、系统及计算机存储介质
CN101867832B (zh) 阈值确定装置和阈值确定方法
CN109389560B (zh) 一种自适应加权滤波图像降噪方法、装置及图像处理设备
CN110730336B (zh) 一种去马赛克方法及装置
CN102509294B (zh) 一种基于单幅图像的全局深度估计方法
Chaudhry et al. A framework for outdoor RGB image enhancement and dehazing
CN109191406A (zh) 图像处理方法、装置及设备
CN104680482A (zh) 一种图像处理的方法和装置
CN107169973A (zh) 一种图像的背景去除和合成方法及装置
Hanumantharaju et al. Color image enhancement using multiscale retinex with modified color restoration technique
Liu et al. Low-light video image enhancement based on multiscale retinex-like algorithm
CN107424124B (zh) 一种图像增强方法及装置
CN107993189B (zh) 一种基于局部分块的图像色调动态调节方法和装置
CN105427265B (zh) 一种彩色图像对比度增强方法及系统
CN110378893A (zh) 图像质量评价方法、装置和电子设备
CN110570384A (zh) 一种对场景图像进行光照均衡处理的方法、装置、计算机设备以及计算机存储介质
CN105574826B (zh) 遥感影像的薄云去除方法
CN108377372B (zh) 一种白平衡处理方法、装置、终端设备和存储介质
US7978910B2 (en) Method and apparatus for adaptively filtering input image in color domains
CN112435182B (zh) 图像降噪方法及装置
Kim et al. On recent results in demosaicing of Samsung 108MP CMOS sensor using deep learning

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