CN104680483B - 图像的噪声估计方法、视频图像去噪方法及装置 - Google Patents
图像的噪声估计方法、视频图像去噪方法及装置 Download PDFInfo
- Publication number
- CN104680483B CN104680483B CN201310603614.4A CN201310603614A CN104680483B CN 104680483 B CN104680483 B CN 104680483B CN 201310603614 A CN201310603614 A CN 201310603614A CN 104680483 B CN104680483 B CN 104680483B
- Authority
- CN
- China
- Prior art keywords
- image
- image block
- determining
- current frame
- value
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 61
- 238000009499 grossing Methods 0.000 claims abstract description 26
- 238000001514 detection method Methods 0.000 claims description 37
- 230000003068 static effect Effects 0.000 claims description 24
- 230000000903 blocking effect Effects 0.000 claims description 23
- 238000005192 partition Methods 0.000 abstract 1
- 238000000638 solvent extraction Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 238000012986 modification Methods 0.000 description 3
- 230000004048 modification Effects 0.000 description 3
- 230000000694 effects Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
- 230000002123 temporal effect Effects 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Landscapes
- Picture Signal Circuits (AREA)
- Image Analysis (AREA)
Abstract
本发明实施例提供了图像的噪声估计方法、视频图像去噪方法及装置,确定待处理图像的两个方向的拉普拉斯算子,其中,对两个方向的拉普拉斯算子进行预设线性运算得到的差异算子中,各行数值之和以及各列数值之和均为零,且位于中心的数值的绝对值大于位于非中心位置的数值的绝对值;利用差异算子对待处理图像进行平滑处理,得到平滑图像;基于第一预设尺寸,将平滑图像进行不重叠分块,得到各图像块;分别确定各图像块的方差,并从各方差中,按照方差值从小到大的顺序确定预设数量个方差,并将预设数量个方差的加权平均值确定为估计的待处理图像的噪声方差。解决现有技术中对图像的噪声估计不准确的问题,本发明涉及图像处理领域。
Description
技术领域
本发明涉及图像处理领域,尤其涉及图像的噪声估计方法、视频图像去噪方法及装置。
背景技术
图像或者视频在传输的过程中,会受到各种噪声的污染,而导致接收方接收到的图像或视频与原始的图像或视频相比分辨率降低,不仅影响视觉效果,对于需要从中获取或识别运动目标的图像或视频,更是影响了获取或识别工作的准确性。因此,现有技术中需要对图像或者视频帧图像中的噪声进行估计,以便为后续的去噪处理做准备。
现有技术中,对图像或者视频帧图像中的噪声的估计方法主要包括如下步骤:
步骤1、对待估计图像或者视频帧图像进行不重叠分块处理;
步骤2、对步骤1中分块处理得到的各图像块分别进行方差计算;
步骤3、从步骤2计算出的各方差中,选取多个值较小方差进行噪声估计。
但是,由于图像或者视频帧图像中的目标存在边缘,而目标的边缘影响了图像的平滑性,因为平滑块的像素值波动可认为是由噪声引起的,因此,通过估计平滑块的方差,就可以获得原始图像中噪声方差的估计,也就是说,经过步骤1的分块处理之后,不存在目标边缘的图像块会比较平滑,计算出的方差与噪声方差很相近,而存在目标边缘的图像块不够平滑,导致计算出的方差与噪声方差相比相差很大,对于低噪声图像,由于噪声较小,边缘块估计出的噪声方差会大于实际噪声方差,导致过高地估计其噪声,对于高噪声图像,由于噪声很大,边缘块估计出的噪声会小于实际噪声方差,导致过低地估计其噪声,使得噪声估计不准确。
发明内容
本发明实施例提供了图像的噪声估计方法、视频图像去噪方法及装置,用以解决现有技术中对图像的噪声估计不准确的问题。
基于上述问题,本发明实施例提供的一种图像的噪声估计方法,包括:
确定待处理图像的两个方向的拉普拉斯算子,其中,所述两个方向的拉普拉斯算子满足如下条件:对所述两个方向的拉普拉斯算子进行预设线性运算得到的差异算子中,各行数值之和以及各列数值之和均为零,且位于中心的数值的绝对值大于位于非中心位置的数值的绝对值;
利用所述差异算子对所述待处理图像进行平滑处理,得到平滑图像;
基于第一预设尺寸,将所述平滑图像进行不重叠分块,得到各图像块;
分别确定所述各图像块的方差,并
从各方差中,按照方差值从小到大的顺序确定预设数量个方差,并将所述预设数量个方差的加权平均值确定为估计的所述待处理图像的噪声方差。
本发明实施例提供的一种视频图像去噪方法,包括:
基于上述一种图像的噪声估计方法,确定待处理视频的当前帧图像的噪声方差,并根据所述噪声方差确定运动检测阈值;
基于第二预设尺寸,分别将所述当前帧图像和前一帧图像进行不重叠分块,得到所述当前帧图像的各图像块和所述前一帧图像的各图像块;
针对所述当前帧图像的每个图像块,确定该图像块与所述前一帧图像中与该图像块位置对应的图像块之间的相似度;并
将所述相似度与所述运动检测阈值进行比较,并根据比较结果确定该图像块为运动图像块或者静止图像块;
分别对所述当前帧图像中的运动图像块和静止图像块进行去噪处理。
本发明实施例提供的一种图像的噪声估计装置,包括:
算子确定模块,用于确定待处理图像的两个方向的拉普拉斯算子,其中,所述两个方向的拉普拉斯算子满足如下条件:对所述两个方向的拉普拉斯算子进行预设线性运算得到的差异算子中,各行数值之和以及各列数值之和均为零,且位于中心的数值的绝对值大于位于非中心位置的数值的绝对值;
平滑处理模块,用于利用所述差异算子对所述待处理图像进行平滑处理,得到平滑图像;
第一分块模块,用于基于第一预设尺寸,将所述平滑图像进行不重叠分块,得到各图像块;
方差确定模块,用于分别确定所述各图像块的方差,并
噪声估计模块,用于从各方差中,按照方差值从小到大的顺序确定预设数量个方差,并将所述预设数量个方差的加权平均值确定为估计的所述待处理图像的噪声方差。
本发明实施例提供的一种视频图像去噪装置,包括:
阈值确定模块,用于通过上述图像的噪声估计装置,确定待处理视频的当前帧图像的噪声方差,并根据所述噪声方差确定运动检测阈值;
第二分块模块,用于基于第二预设尺寸,分别将所述当前帧图像和前一帧图像进行不重叠分块,得到所述当前帧图像的各图像块和所述前一帧图像的各图像块;
相似度确定模块,用于针对所述当前帧图像的每个图像块,确定该图像块与所述前一帧图像中与该图像块位置对应的图像块之间的相似度;并
运动检测模块,用于将所述相似度与所述运动检测阈值进行比较,并根据比较结果确定该图像块为运动图像块或者静止图像块;
去噪处理模块,用于分别对所述当前帧图像中的运动图像块和静止图像块进行去噪处理。
本发明实施例的有益效果包括:
本发明实施例提供的图像的噪声估计方法、视频图像去噪方法及装置,确定待处理图像的两个方向的拉普拉斯算子,其中,两个方向的拉普拉斯算子满足如下条件:对两个方向的拉普拉斯算子进行预设线性运算得到的差异算子中,各行数值之和以及各列数值之和均为零,且位于中心的数值的绝对值大于位于非中心位置的数值的绝对值;利用差异算子对待处理图像进行平滑处理,得到平滑图像;基于第一预设尺寸,将平滑图像进行不重叠分块,得到各图像块;分别确定各图像块的方差,并从各方差中,按照方差值从小到大的顺序确定预设数量个方差,并将预设数量个方差的加权平均值确定为估计的待处理图像的噪声方差。在对待处理图像进行分块之前,对图像进行平滑处理,使得图像中的目标边缘变得平滑,那么,划分的各图像块中的图像变得平滑,各图像块计算的方差可以较为正常,利用方差进行图像的噪声估计更加准确;
进一步地,利用上述噪声估计方法确定待处理视频的当前帧图像的噪声方差,并根据噪声方差确定运动检测阈值;基于第二预设尺寸,分别将当前帧图像和前一帧图像进行不重叠分块,得到当前帧图像的各图像块和前一帧图像的各图像块;针对当前帧图像的每个图像块,确定该图像块与前一帧图像中与该图像块位置对应的图像块之间的相似度;并将相似度与运动检测阈值进行比较,并根据比较结果确定该图像块为运动图像块或者静止图像块;分别对当前帧图像中的运动图像块和静止图像块进行去噪处理。在上述对视频图像进行去噪处理的过程中,基于准确的噪声估计方法对图像进行噪声估计,再进行去噪处理,使得去噪处理效果更好。
附图说明
图1为本发明实施例提供的一种图像的噪声估计方法的流程图;
图2为本发明实施例提供的一种视频图像去噪方法的流程图;
图3为本发明实施例1提供的一种图像的噪声估计方法的流程图;
图4a-图4c为本发明实施例1提供的拉普拉斯算子模板举例;
图5为本发明实施例2提供的一种视频图像去噪方法的流程图;
图6为本发明实施例3提供的一种视频图像去噪方法的流程图;
图7为本发明实施例提供的一种图像的噪声估计装置的结构示意图;
图8为本发明实施例提供的一种视频图像去噪装置的结构示意图。
具体实施方式
本发明实施例提供了图像的噪声估计方法、视频图像去噪方法及装置,以下结合说明书附图对本发明的优选实施例进行说明,应当理解,此处所描述的优选实施例仅用于说明和解释本发明,并不用于限定本发明。并且在不冲突的情况下,本申请中的实施例及实施例中的特征可以相互组合。
本发明实施例提供一种图像的噪声估计方法,如图1所示,包括:
S101、确定待处理图像的两个方向的拉普拉斯算子,其中,该两个方向的拉普拉斯算子满足如下条件:对该两个方向的拉普拉斯算子进行预设线性运算得到的差异算子中,各行数值之和以及各列数值之和均为零,且位于中心的数值的绝对值大于位于非中心位置的数值的绝对值。
S102、利用S101中确定的差异算子对待处理图像进行平滑处理,得到平滑图像。
S103、基于第一预设尺寸,将S102中得到的平滑图像进行不重叠分块,得到各图像块。
S104、分别确定S103中确定的各图像块的方差。
S105、从S104中确定的各方差中,按照方差值从小到大的顺序确定预设数量个方差。
S106、并将S105中确定的预设数量个方差的加权平均值确定为估计的待处理图像的噪声方差。
相应的,基于上述图像的噪声估计方法,本发明实施例还提供一种视频图像去噪方法,如图2所示,包括:
S201、基于上述图像的噪声估计方法,确定待处理视频的当前帧图像的噪声方差,并根据噪声方差确定运动检测阈值。
S202、基于第二预设尺寸,分别将当前帧图像和前一帧图像进行不重叠分块,得到当前帧图像的各图像块和前一帧图像的各图像块。
S203、针对当前帧图像的每个图像块,确定该图像块与前一帧图像中与该图像块位置对应的图像块之间的相似度。
S204、将S203中确定的相似度与S201中确定的运动检测阈值进行比较,并根据比较结果确定该图像块为运动图像块或者静止图像块。
S205、分别对当前帧图像中的运动图像块和静止图像块进行去噪处理。
进一步地,上述步骤S201的执行与步骤S202-步骤S203没有先后顺序。
下面结合附图,用具体实施例对本发明提供的方法及相关设备进行详细描述。
实施例1:
本发明实施例1中,提供了一种图像的噪声估计方法,如图3所示,具体包括如下步骤:
S301、确定待处理图像的两个方向的拉普拉斯算子,其中,该两个方向的拉普拉斯算子满足如下条件:对该两个方向的拉普拉斯算子进行预设线性运算得到的差异算子中,各行数值之和以及各列数值之和均为零,且位于中心的数值的绝对值大于位于非中心位置的数值的绝对值。
进一步地,假设确定的待处理图像的两个方向的拉普拉斯算子进行如式(1)的线性运算,得到差异算子dL:
dL=aL1-bL2;(1)
其中,L1和L2分别表示确定的两个方向的拉普拉斯算子;且差异算子dL满足各行数值之和以及各列数值之和均为零,且位于中心的数值的绝对值大于位于其他位置的数值的绝对值。
下面举例说明:
假设a和b分别取2和1,图4a-图4c分别为L1、L2和dL对应的模板。如图4a-图4c所示,由于dL对应的模板中,各行数值之和均为0,各列数值之和均为0,且位于中间的数值的绝对值4,大于位于其他位置的数值的绝对值,因此,L1、L2能够满足条件,可以将L1、L2确定为选择的两个方向的拉普拉斯算子。
进一步地,现有技术中可以利用拉普拉斯算子的零交叉性质进行图像中的目标边缘定位,而本实施例中,由于差异算子的性质,可以使用差异算子对待处理图像进行平滑处理。
S302、利用S301中确定的差异算子采用如下公式对所述待处理图像进行平滑处理:
其中,(x,y)表示所述待处理图像中的像素点;v(x,y)表示对像素点(x,y)进行平滑处理之前,像素点(x,y)的灰度值;n(x,y)表示对像素点(x,y)进行平滑处理之后,像素点(x,y)的灰度值;l的值为拉普拉斯算子模板的行数或者列数减1除以2。
进一步地,(2)式中以(x,y)为中心像素点的(2l+1)×(2l+1)图像块与拉普拉斯算子模板相对应,对该中心像素点(x,y)进行平滑处理之后的灰度值n(x,y)为使用拉普拉斯模板dL(x+i,y+j)中的每个值分别与图像块中与该值位置对应的像素点的灰度值相乘,再将得到的各乘积相加的得到的值来代替中心像素(x,y)的灰度值即为处理后的灰度值n(x,y)。
S303、基于第一预设尺寸,将S302中得到的平滑图像进行不重叠分块,得到各图像块。
S304、分别确定S303中确定的各图像块的方差。
S305、从S304中确定的各方差中,按照方差值从小到大的顺序确定预设数量个方差。
较佳地,本步骤中可以将各方差由小到大进行排序,并确定前6%的方差。
S306、并将S305中确定的预设数量个方差的加权平均值确定为估计的待处理图像的噪声方差。
本步骤中,确定的各个方差的权值可以相等,即将预设数量个方差的平均值确定为估计的待处理图像的噪声方差。
实施例2:
本发明实施例2中,提供了一种视频图像去噪方法,如图5所示,具体包括如下步骤:
S501、基于上述图像的噪声估计方法,确定待处理视频的当前帧图像的噪声方差,并根据噪声方差确定运动检测阈值。
进一步地,可以根据实际需要使用预设系数与噪声方差相乘,并将乘积确定为运动检测阈值。
进一步地,本步骤与步骤S502-S503没有严格的执行顺序,也就是说,只要在步骤S504之前确定出运动检测阈值即可。
S502、基于第二预设尺寸,分别将当前帧图像和前一帧图像进行不重叠分块,得到当前帧图像的各图像块和前一帧图像的各图像块。
S503、针对当前帧图像的每个图像块,确定该图像块与前一帧图像中与该图像块位置对应的图像块之间的相似度。
S504、将S503中确定的相似度与S501中确定的运动检测阈值进行比较,并根据比较结果确定该图像块为运动图像块或者静止图像块。
S505、分别对当前帧图像中的运动图像区域和静止图像区域进行去噪处理。
S506、确定当前帧图像中每个像素点的梯度值。
进一步地,本步骤与步骤S501-步骤S505的执行没有严格的先后顺序,也就是说,只要在执行步骤S507之前,确定出当前帧图像中每个像素点的梯度值即可。
进一步地,也可以将本步骤修改为:确定当前帧图像中的运动图像块中每个像素点的梯度值。
S507、从运动图像块中的各像素点中,确定梯度值大于预设梯度阈值的像素点。
进一步地,本步骤中,梯度值大于预设梯度阈值的像素点可以认为是作为运动目标的边缘的像素点。
S508、针对确定的像素点中的每个像素点,将对该像素点进行去噪处理之后得到的灰度值,与预设比例的对该像素点进行去噪处理之前的灰度值的和,进一步确定为对该像素点进行去噪处理之后的灰度值。
进一步地,对于当前帧图像中的运动图像块中和静止图像块中,作为非运动目标边缘的像素点,可以将去噪处理之后,得到的像素点的灰度值直接作为去噪处理之后该像素点的灰度值,而对于作为运动目标边缘的像素点,还需要按照本步骤进行进一步的处理。经过进一步地处理,能够使得图像中运动目标的细小边缘完整地保留,使得处理后的图像不会出现拖影现象。
本步骤中,假设作为运动目标边缘的像素点,进行去噪处理之前灰度值为Y,进行去噪处理之后灰度值为Y',预设比例为α,则该像素点进行去噪处理之后,最终灰度值可以为:Y'+αY。
实施例3:
本发明实施例3中,提供了一种视频图像去噪方法,如图6所示,具体包括如下步骤:
S601、基于上述图像的噪声估计方法,确定待处理视频的当前帧图像的噪声方差,并根据噪声方差确定运动检测阈值。
进一步地,本步骤与步骤S602-S604没有严格的执行顺序,也就是说,只要在步骤S605之前确定出运动检测阈值即可。
S602、基于第二预设尺寸,分别将当前帧图像和前一帧图像进行不重叠分块,得到当前帧图像的各图像块和前一帧图像的各图像块。
S603、针对当前帧图像和前一帧图像的每个图像块,确定该图像块所在帧中,以该图像块为中心的s×s个图像块构成的图像区域。
S604、针对当前帧图像的每个图像块,采用如下公式确定该图像块与前一帧图像中与该图像块位置对应的图像块之间的相似度D(i,j,k):
其中,(i,j)表示该图像块的中心坐标;h的值为(s-1)/2;k和k-1表示帧号;(m,n)表示以该图像块为中心的图像区域中除该图像块之外的图像块相对于该图像块的坐标;SB(m,n,k)表示第k帧中中心坐标为(m,n)的图像块中所有像素点灰度值的和;SB(m,n,k-1)表示第k-1帧中中心坐标为(m,n)的图像块中所有像素点灰度值的和。
进一步地,假设s的值为3,也就是说,针对当前帧图像和前一帧图像的每个图像块,确定该图像块以及该图像块所在帧中,以该图像块为中心的3×3个图像块构成的图像区域,则以该图像块为中心的图像区域中,包括9个图像块。那么式(3)中,分别计算当前帧的9个图像块中每个图像块中所有像素点灰度值的和,以及前一帧的9个图像块中每个图像块中所有像素点灰度值的和,并分别确定当前帧中以及前一帧中各自的9个图像块中,每两个位置对应的图像块中所有像素点灰度值的和的差值的平方,并将得到的9个差值的平方相加,得到的和作为该图像块的相似度。
S605、当S604中确定的相似度大于S601中确定的运动检测阈值时,确定当前帧的该图像块相对上一帧中与该图像块位置对应的图像块为运动图像块;
当S604中确定的相似度不大于S601中确定的运动检测阈值时,确定当前帧的该图像块相对上一帧中与该图像块位置对应的图像块为静止图像块。
进一步地,本步骤中,当D(i,j,k)大于运动检测阈值时,可以将(i,j)表征的图像块中的每个像素点标记为1,否则,标记为0,那么被标记为1的像素点构成的图像块表征运动图像块,被标记为0的像素点构成的图像块表征静止图像块。
S606、分别对当前帧图像中的运动图像块和静止图像块进行去噪处理。
本步骤中,可以采用改进的3D-Rational滤波公式对当前帧图像中的静止图像块进行去噪处理,如式(4)所示:
针对静止图像块中的每个像素点,假设该像素点进行去噪处理之前的灰度值为x0,该像素点进行去噪处理之后的灰度值为y0,t和t-1分别用来标识帧号(不表示幂值),Ω表示与该像素点相邻的九个像素点的索引集合,i和j分别表示Ω中的任一像素点,xi和xj分别表示i像素点和j像素点的灰度值,kt,At均表示时域去噪调节因子。
本步骤中,可以采用改进的快速NLM方法,对当前帧图像中的运动图像块进行去噪处理,如式(5)所示:
其中,
在式(5)中,NLv(i)表示对像素点i进行去噪处理之后像素点i的灰度值,像素点j为与像素点i相邻的像素点,v(j)表示对像素点j进行去噪处理之前像素点j的灰度值,w(i,j)表示利用像素点j进行去噪处理之前像素点j的灰度值确定像素点i进行去噪处理之后像素点i的灰度值时,使用的各个权值。
式(6)为权值w(i,j)的计算公式,其中,||N(i)-N(j)||2表示像素点i和像素点j之间的欧式距离,h为衰减因子,▽v(i)表示像素点i的梯度值,▽v(j)表示像素点j的梯度值,σ▽表示另一个梯度阈值,表示像素点i的均值,表示像素点j的均值,η1和η2分别表示第一均值阈值和第二均值阈值,较佳地,η1和η2可以分别取值为0.9和1.2。
基于同一发明构思,本发明实施例还提供了一种图像的噪声估计装置,由于这些装置所解决问题的原理与前述一种图像的噪声估计方法相似,因此该装置的实施可以参见前述方法的实施,重复之处不再赘述。
本发明实施例提供的一种图像的噪声估计装置,如图7所示,包括如下模块:
算子确定模块701,用于确定待处理图像的两个方向的拉普拉斯算子,其中,所述两个方向的拉普拉斯算子满足如下条件:对所述两个方向的拉普拉斯算子进行预设线性运算得到的差异算子中,各行数值之和以及各列数值之和均为零,且位于中心的数值的绝对值大于位于非中心位置的数值的绝对值;
平滑处理模块702,用于利用所述差异算子对所述待处理图像进行平滑处理,得到平滑图像;
第一分块模块703,用于基于第一预设尺寸,将所述平滑图像进行不重叠分块,得到各图像块;
方差确定模块704,用于分别确定所述各图像块的方差,并
噪声估计模块705,用于从各方差中,按照方差值从小到大的顺序确定预设数量个方差,并将所述预设数量个方差的加权平均值确定为估计的所述待处理图像的噪声方差。
进一步地,所述平滑处理模块702,具体用于利用所述差异算子,采用如下公式对所述待处理图像进行平滑处理:
其中,(x,y)表示所述待处理图像中的像素点;v(x,y)表示对像素点(x,y)进行平滑处理之前,像素点(x,y)的灰度值;n(x,y)表示对像素点(x,y)进行平滑处理之后,像素点(x,y)的灰度值;l的值为所述差异算子模板的行数或者列数减1除以2。
基于同一发明构思,本发明实施例还提供了一种基于上述图像的噪声估计装置的视频图像去噪装置,由于这些装置所解决问题的原理与前述视频图像去噪方法相似,因此该装置的实施可以参见前述方法的实施,重复之处不再赘述。
本发明实施例提供的一种视频图像去噪装置,如图8所示,包括如下模块:
阈值确定模块801,用于通过上述图像的噪声估计装置,确定待处理视频的当前帧图像的噪声方差,并根据所述噪声方差确定运动检测阈值;
第二分块模块802,用于基于第二预设尺寸,分别将所述当前帧图像和前一帧图像进行不重叠分块,得到所述当前帧图像的各图像块和所述前一帧图像的各图像块;
相似度确定模块803,用于针对所述当前帧图像的每个图像块,确定该图像块与所述前一帧图像中与该图像块位置对应的图像块之间的相似度;并
运动检测模块804,用于将所述相似度与所述运动检测阈值进行比较,并根据比较结果确定该图像块为运动图像块或者静止图像块;
去噪处理模块805,用于分别对所述当前帧图像中的运动图像块和静止图像块进行去噪处理。
进一步地,上述视频图像去噪装置,还包括:
梯度值确定模块806,用于确定所述当前帧图像中每个像素点的梯度值;
像素点确定模块807,用于在所述去噪处理模块分别对所述当前帧图像中的运动图像块和静止图像块进行去噪处理之后,从所述运动图像块中的各像素点中,确定梯度值大于预设梯度阈值的像素点;
边缘处理模块808,用于针对确定的像素点中的每个像素点,将对该像素点进行去噪处理之后得到的灰度值,与预设比例的对该像素点进行去噪处理之前的灰度值的和,进一步确定为对该像素点进行去噪处理之后的灰度值。
进一步地,所述相似度确定模块803,具体用于针对所述当前帧图像和所述前一帧图像的每个图像块,确定该图像块以及该图像块所在帧中,以该图像块为中心的s×s个图像块构成的图像区域;
针对所述当前帧的每个图像块,采用如下公式确定该图像块与所述前一帧中与该图像块位置对应的图像块之间的相似度D(i,j,k):
其中,(i,j)表示该图像块的中心坐标;h的值为(s-1)/2;k和k-1表示帧号;(m,n)表示以该图像块为中心的图像区域中除该图像块之外的图像块相对于该图像块的坐标;SB(m,n,k)表示第k帧中中心坐标为(m,n)的图像块中所有像素点灰度值的和;SB(m,n,k-1)表示第k-1帧中中心坐标为(m,n)的图像块中所有像素点灰度值的和。
进一步地,所述运动检测模块804,具体用于当所述相似度大于所述运动检测阈值时,确定当前帧的该图像块相对上一帧中与该图像块位置对应的图像块为运动图像块;当所述相似度不大于所述运动检测阈值时,确定当前帧的该图像块相对上一帧中与该图像块位置对应的图像块为静止图像块。
上述各单元的功能可对应于图1至图3、图5-图6所示流程中的相应处理步骤,在此不再赘述。
本发明实施例提供的图像的噪声估计方法、视频图像去噪方法及装置,确定待处理图像的两个方向的拉普拉斯算子,其中,两个方向的拉普拉斯算子满足如下条件:对两个方向的拉普拉斯算子进行预设线性运算得到的差异算子中,各行数值之和以及各列数值之和均为零,且位于中心的数值的绝对值大于位于非中心位置的数值的绝对值;利用差异算子对待处理图像进行平滑处理,得到平滑图像;基于第一预设尺寸,将平滑图像进行不重叠分块,得到各图像块;分别确定各图像块的方差,并从各方差中,按照方差值从小到大的顺序确定预设数量个方差,并将预设数量个方差的加权平均值确定为估计的待处理图像的噪声方差。在对待处理图像进行分块之前,对图像进行平滑处理,使得图像中的目标边缘变得平滑,那么,划分的各图像块中的图像变得平滑,各图像块计算的方差可以较为正常,利用方差进行图像的噪声估计更加准确;
进一步地,利用上述噪声估计方法确定待处理视频的当前帧图像的噪声方差,并根据噪声方差确定运动检测阈值;基于第二预设尺寸,分别将当前帧图像和前一帧图像进行不重叠分块,得到当前帧图像的各图像块和前一帧图像的各图像块;针对当前帧图像的每个图像块,确定该图像块与前一帧图像中与该图像块位置对应的图像块之间的相似度;并将相似度与运动检测阈值进行比较,并根据比较结果确定该图像块为运动图像块或者静止图像块;分别对当前帧图像中的运动图像块和静止图像块进行去噪处理。在上述对视频图像进行去噪处理的过程中,基于准确的噪声估计方法对图像进行噪声估计,再进行去噪处理,使得去噪处理效果更好。
通过以上的实施方式的描述,本领域的技术人员可以清楚地了解到本发明实施例可以通过硬件实现,也可以借助软件加必要的通用硬件平台的方式来实现。基于这样的理解,本发明实施例的技术方案可以以软件产品的形式体现出来,该软件产品可以存储在一个非易失性存储介质(可以是CD-ROM,U盘,移动硬盘等)中,包括若干指令用以使得一台计算机设备(可以是个人计算机,服务器,或者网络设备等)执行本发明各个实施例所述的方法。
本领域技术人员可以理解附图只是一个优选实施例的示意图,附图中的模块或流程并不一定是实施本发明所必须的。
本领域技术人员可以理解实施例中的装置中的模块可以按照实施例描述进行分布于实施例的装置中,也可以进行相应变化位于不同于本实施例的一个或多个装置中。上述实施例的模块可以合并为一个模块,也可以进一步拆分成多个子模块。
上述本发明实施例序号仅仅为了描述,不代表实施例的优劣。
显然,本领域的技术人员可以对本发明进行各种改动和变型而不脱离本发明的精神和范围。这样,倘若本发明的这些修改和变型属于本发明权利要求及其等同技术的范围之内,则本发明也意图包含这些改动和变型在内。
Claims (12)
1.一种图像的噪声估计方法,其特征在于,包括:
确定待处理图像的两个方向的拉普拉斯算子,其中,所述两个方向的拉普拉斯算子满足如下条件:对所述两个方向的拉普拉斯算子进行预设线性运算得到的差异算子中,各行数值之和以及各列数值之和均为零,且位于中心的数值的绝对值大于位于非中心位置的数值的绝对值,预设线性运算公式为:dL=aL1-bL2,其中,L1和L2分别表示确定的两个方向的拉普拉斯算子,dL表示所述差异算子,a和b分别为常数,所述确定的两个方向的拉普拉斯算子的行数和列数相等且为奇数;
利用所述差异算子对所述待处理图像进行平滑处理,得到平滑图像;
基于第一预设尺寸,将所述平滑图像进行不重叠分块,得到各图像块;
分别确定所述各图像块的方差,并
从各方差中,按照方差值从小到大的顺序确定预设数量个方差,并将所述预设数量个方差的加权平均值确定为估计的所述待处理图像的噪声方差。
2.如权利要求1所述的方法,其特征在于,利用所述差异算子,采用如下公式对所述待处理图像进行平滑处理:
其中,(x,y)表示所述待处理图像中的像素点;v(x,y)表示对像素点(x,y)进行平滑处理之前,像素点(x,y)的灰度值;n(x,y)表示对像素点(x,y)进行平滑处理之后,像素点(x,y)的灰度值;l的值为拉普拉斯算子模板的行数或者列数减1除以2。
3.一种视频图像去噪方法,其特征在于,包括:
基于权利要求1-2任一所述图像的噪声估计方法,确定待处理视频的当前帧图像的噪声方差,并根据所述噪声方差确定运动检测阈值;
基于第二预设尺寸,分别将所述当前帧图像和前一帧图像进行不重叠分块,得到所述当前帧图像的各图像块和所述前一帧图像的各图像块;
针对所述当前帧图像的每个图像块,确定该图像块与所述前一帧图像中与该图像块位置对应的图像块之间的相似度;并
将所述相似度与所述运动检测阈值进行比较,并根据比较结果确定该图像块为运动图像块或者静止图像块;
分别对所述当前帧图像中的运动图像块和静止图像块进行去噪处理。
4.如权利要求3所述的方法,其特征在于,还包括:
确定所述当前帧图像中每个像素点的梯度值;
分别对所述当前帧图像中的运动图像块和静止图像块进行去噪处理之后,还包括:
从所述运动图像块中的各像素点中,确定梯度值大于预设梯度阈值的像素点;
针对确定的像素点中的每个像素点,将对该像素点进行去噪处理之后得到的灰度值,与预设比例的对该像素点进行去噪处理之前的灰度值的和,进一步确定为对该像素点进行去噪处理之后的灰度值。
5.如权利要求3所述的方法,其特征在于,针对所述当前帧图像的每个图像块,确定该图像块与所述前一帧图像中与该图像块位置对应的图像块之间的相似度,具体包括:
针对所述当前帧图像和所述前一帧图像的每个图像块,确定该图像块所在帧中,以该图像块为中心的s×s个图像块构成的图像区域;
针对所述当前帧的每个图像块,采用如下公式确定该图像块与所述前一帧中与该图像块位置对应的图像块之间的相似度D(i,j,k):
其中,(i,j)表示该图像块的中心坐标;h的值为(s-1)/2;k和k-1表示帧号;(m,n)表示以该图像块为中心的图像区域中除该图像块之外的图像块相对于该图像块的坐标;SB(m,n,k)表示第k帧中中心坐标为(m,n)的图像块中所有像素点灰度值的和;SB(m,n,k-1)表示第k-1帧中中心坐标为(m,n)的图像块中所有像素点灰度值的和。
6.如权利要求3所述的方法,其特征在于,将所述相似度与所述运动检测阈值进行比较,并根据比较结果确定该图像块为运动图像块或者静止图像块,具体包括:
当所述相似度大于所述运动检测阈值时,确定当前帧的该图像块相对上一帧中与该图像块位置对应的图像块为运动图像块;
当所述相似度不大于所述运动检测阈值时,确定当前帧的该图像块相对上一帧中与该图像块位置对应的图像块为静止图像块。
7.一种图像的噪声估计装置,其特征在于,包括:
算子确定模块,用于确定待处理图像的两个方向的拉普拉斯算子,其中,所述两个方向的拉普拉斯算子满足如下条件:对所述两个方向的拉普拉斯算子进行预设线性运算得到的差异算子中,各行数值之和以及各列数值之和均为零,且位于中心的数值的绝对值大于位于非中心位置的数值的绝对值,预设线性运算公式为:dL=aL1-bL2,其中,L1和L2分别表示确定的两个方向的拉普拉斯算子,dL表示所述差异算子,a和b分别为常数,所述确定的两个方向的拉普拉斯算子的行数和列数相等且为奇数;
平滑处理模块,用于利用所述差异算子对所述待处理图像进行平滑处理,得到平滑图像;
第一分块模块,用于基于第一预设尺寸,将所述平滑图像进行不重叠分块,得到各图像块;
方差确定模块,用于分别确定所述各图像块的方差,并
噪声估计模块,用于从各方差中,按照方差值从小到大的顺序确定预设数量个方差,并将所述预设数量个方差的加权平均值确定为估计的所述待处理图像的噪声方差。
8.如权利要求7所述的装置,其特征在于,所述平滑处理模块,具体用于利用所述差异算子,采用如下公式对所述待处理图像进行平滑处理:
其中,(x,y)表示所述待处理图像中的像素点;v(x,y)表示对像素点(x,y)进行平滑处理之前,像素点(x,y)的灰度值;n(x,y)表示对像素点(x,y)进行平滑处理之后,像素点(x,y)的灰度值;l的值为拉普拉斯算子模板的行数或者列数减1除以2。
9.一种视频图像去噪装置,其特征在于,包括:
阈值确定模块,用于通过权利要求7-8任一所述图像的噪声估计装置,确定待处理视频的当前帧图像的噪声方差,并根据所述噪声方差确定运动检测阈值;
第二分块模块,用于基于第二预设尺寸,分别将所述当前帧图像和前一帧图像进行不重叠分块,得到所述当前帧图像的各图像块和所述前一帧图像的各图像块;
相似度确定模块,用于针对所述当前帧图像的每个图像块,确定该图像块与所述前一帧图像中与该图像块位置对应的图像块之间的相似度;并
运动检测模块,用于将所述相似度与所述运动检测阈值进行比较,并根据比较结果确定该图像块为运动图像块或者静止图像块;
去噪处理模块,用于分别对所述当前帧图像中的运动图像块和静止图像块进行去噪处理。
10.如权利要求9所述的装置,其特征在于,还包括:
梯度值确定模块,用于确定所述当前帧图像中每个像素点的梯度值;
像素点确定模块,用于在所述去噪处理模块分别对所述当前帧图像中的运动图像块和静止图像块进行去噪处理之后,从所述运动图像块中的各像素点中,确定梯度值大于预设梯度阈值的像素点;
边缘处理模块,用于针对确定的像素点中的每个像素点,将对该像素点进行去噪处理之后得到的灰度值,与预设比例的对该像素点进行去噪处理之前的灰度值的和,进一步确定为对该像素点进行去噪处理之后的灰度值。
11.如权利要求9所述的装置,其特征在于,所述相似度确定模块,具体用于针对所述当前帧图像和所述前一帧图像的每个图像块,确定该图像块所在帧中,以该图像块为中心的s×s个图像块构成的图像区域;
针对所述当前帧的每个图像块,采用如下公式确定该图像块与所述前一帧中与该图像块位置对应的图像块之间的相似度D(i,j,k):
其中,(i,j)表示该图像块的中心坐标;h的值为(s-1)/2;k和k-1表示帧号;(m,n)表示以该图像块为中心的图像区域中除该图像块之外的图像块相对于该图像块的坐标;SB(m,n,k)表示第k帧中中心坐标为(m,n)的图像块中所有像素点灰度值的和;SB(m,n,k-1)表示第k-1帧中中心坐标为(m,n)的图像块中所有像素点灰度值的和。
12.如权利要求9所述的装置,其特征在于,所述运动检测模块,具体用于当所述相似度大于所述运动检测阈值时,确定当前帧的该图像块相对上一帧中与该图像块位置对应的图像块为运动图像块;当所述相似度不大于所述运动检测阈值时,确定当前帧的该图像块相对上一帧中与该图像块位置对应的图像块为静止图像块。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310603614.4A CN104680483B (zh) | 2013-11-25 | 2013-11-25 | 图像的噪声估计方法、视频图像去噪方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310603614.4A CN104680483B (zh) | 2013-11-25 | 2013-11-25 | 图像的噪声估计方法、视频图像去噪方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104680483A CN104680483A (zh) | 2015-06-03 |
CN104680483B true CN104680483B (zh) | 2016-03-02 |
Family
ID=53315482
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310603614.4A Active CN104680483B (zh) | 2013-11-25 | 2013-11-25 | 图像的噪声估计方法、视频图像去噪方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104680483B (zh) |
Families Citing this family (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105894457A (zh) * | 2015-12-04 | 2016-08-24 | 乐视致新电子科技(天津)有限公司 | 一种去除图像噪声的方法及装置 |
CN106127729A (zh) * | 2016-06-08 | 2016-11-16 | 浙江传媒学院 | 一种基于梯度的图像噪声水平估计方法 |
CN110175995B (zh) | 2019-05-29 | 2021-04-30 | 腾讯科技(深圳)有限公司 | 一种基于病理图像的图像状态确定方法、装置以及*** |
CN110363726A (zh) * | 2019-07-24 | 2019-10-22 | 陈泰杉 | 一种基于正则图拉普拉斯变换的视觉降噪处理方法 |
CN112311962B (zh) * | 2019-07-29 | 2023-11-24 | 深圳市中兴微电子技术有限公司 | 一种视频去噪方法和装置、计算机可读存储介质 |
CN113129221B (zh) * | 2019-12-31 | 2023-08-18 | 杭州海康威视数字技术股份有限公司 | 图像处理方法、装置、设备及存储介质 |
CN111556227B (zh) * | 2020-05-19 | 2022-04-15 | 广州市百果园信息技术有限公司 | 一种视频去噪方法、装置、移动终端和存储介质 |
CN111724421B (zh) * | 2020-06-29 | 2024-01-09 | 深圳市慧鲤科技有限公司 | 图像处理方法及装置、电子设备及存储介质 |
CN115330628B (zh) * | 2022-08-18 | 2023-09-12 | 盐城众拓视觉创意有限公司 | 基于图像处理的视频逐帧去噪方法 |
CN118037584B (zh) * | 2024-04-12 | 2024-06-14 | 禹创半导体(深圳)有限公司 | 一种基于过驱动的图像边界平滑方法、设备及介质 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102722874A (zh) * | 2012-05-29 | 2012-10-10 | 飞依诺科技(苏州)有限公司 | 超声图像斑点噪声抑制方法 |
CN103095967A (zh) * | 2011-10-28 | 2013-05-08 | 浙江大华技术股份有限公司 | 一种视频噪声量化计算方法及*** |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8310548B2 (en) * | 2009-10-23 | 2012-11-13 | Contec Llc | System and method for video quality parametric tests |
-
2013
- 2013-11-25 CN CN201310603614.4A patent/CN104680483B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103095967A (zh) * | 2011-10-28 | 2013-05-08 | 浙江大华技术股份有限公司 | 一种视频噪声量化计算方法及*** |
CN102722874A (zh) * | 2012-05-29 | 2012-10-10 | 飞依诺科技(苏州)有限公司 | 超声图像斑点噪声抑制方法 |
Also Published As
Publication number | Publication date |
---|---|
CN104680483A (zh) | 2015-06-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104680483B (zh) | 图像的噪声估计方法、视频图像去噪方法及装置 | |
CN104408707B (zh) | 一种快速数字成像模糊鉴别与复原图像质量评估方法 | |
Wang et al. | Improved scheme of estimating motion blur parameters for image restoration | |
US9798951B2 (en) | Apparatus, method, and processor for measuring change in distance between a camera and an object | |
CN104504120B (zh) | 一种相似图片检测方法及装置 | |
US10255673B2 (en) | Apparatus and method for detecting object in image, and apparatus and method for computer-aided diagnosis | |
US20150254814A1 (en) | Globally dominant point spread function estimation | |
CN104599286B (zh) | 一种基于光流的特征跟踪方法及装置 | |
CN105890768B (zh) | 一种红外图像非均匀性校正的方法及装置 | |
CN108960012B (zh) | 特征点检测方法、装置及电子设备 | |
CN105447828B (zh) | 沿运动模糊路径进行一维反卷积的单视点图像去模糊方法 | |
KR20150032822A (ko) | 이미지를 필터링하기 위한 방법 및 장치 | |
US20150187051A1 (en) | Method and apparatus for estimating image noise | |
Pok et al. | Efficient block matching for removing impulse noise | |
JP2015106290A5 (zh) | ||
US9569866B2 (en) | Flexible video object boundary tracking | |
Kryjak et al. | Real-time implementation of foreground object detection from a moving camera using the vibe algorithm | |
CN108062523B (zh) | 一种红外远小目标检测方法 | |
Bansal et al. | A Comparative Analysis of Iterative and Ostu's Thresholding Techniques | |
CN104376578A (zh) | 一种应用于直录播***中的运动目标检测处理方法及设备 | |
Tian et al. | A fractional-order Laplacian operator for image edge detection | |
CN105469369B (zh) | 基于分割图的数字图像滤波方法和*** | |
Haque et al. | Robust feature-preserving denoising of 3D point clouds | |
CN104182971B (zh) | 一种高精度图像矩定位方法 | |
Grimaldi et al. | Detection and parameters estimation of locally motion blurred objects |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |