[go: up one dir, main page]

CN104680483B - Image noise estimation method, video image denoising method and device - Google Patents

Image noise estimation method, video image denoising method and device Download PDF

Info

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
Application number
CN201310603614.4A
Other languages
Chinese (zh)
Other versions
CN104680483A (en
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.)
Zhejiang Dahua Technology Co Ltd
Original Assignee
Zhejiang Dahua 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 Zhejiang Dahua Technology Co Ltd filed Critical Zhejiang Dahua Technology Co Ltd
Priority to CN201310603614.4A priority Critical patent/CN104680483B/en
Publication of CN104680483A publication Critical patent/CN104680483A/en
Application granted granted Critical
Publication of CN104680483B publication Critical patent/CN104680483B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Picture Signal Circuits (AREA)
  • Image Analysis (AREA)

Abstract

The embodiment of the invention provides a noise estimation method of an image, a denoising method of a video image and a device, which are used for determining Laplacian operators in two directions of the image to be processed, wherein in a difference operator obtained by carrying out preset linear operation on the Laplacian operators in the two directions, the sum of numerical values of each row and the sum of numerical values of each column are zero, and the absolute value of the numerical value in the center is greater than the absolute value of the numerical value in a non-center position; smoothing the image to be processed by using a difference operator to obtain a smooth image; based on the first preset size, carrying out non-overlapping blocking on the smooth image to obtain each image block; and respectively determining the variances of each image block, determining a preset number of variances according to the sequence of variance values from small to large from each variance, and determining the weighted average of the preset number of variances as the estimated noise variance of the image to be processed. The invention solves the problem of inaccurate noise estimation of an image in the prior art, and relates to the field of image processing.

Description

图像的噪声估计方法、视频图像去噪方法及装置Image noise estimation method, video image denoising method and device

技术领域technical field

本发明涉及图像处理领域,尤其涉及图像的噪声估计方法、视频图像去噪方法及装置。The invention relates to the field of image processing, in particular to an image noise estimation method, a video image denoising method and a device.

背景技术Background technique

图像或者视频在传输的过程中,会受到各种噪声的污染,而导致接收方接收到的图像或视频与原始的图像或视频相比分辨率降低,不仅影响视觉效果,对于需要从中获取或识别运动目标的图像或视频,更是影响了获取或识别工作的准确性。因此,现有技术中需要对图像或者视频帧图像中的噪声进行估计,以便为后续的去噪处理做准备。During the transmission of images or videos, they will be polluted by various noises, which will cause the resolution of the images or videos received by the receiver to be reduced compared with the original images or videos, which not only affects the visual effect, but also affects the image or video that needs to be obtained or recognized. Images or videos of moving targets affect the accuracy of acquisition or recognition. Therefore, in the prior art, it is necessary to estimate noise in an image or a video frame image, so as to prepare for subsequent denoising processing.

现有技术中,对图像或者视频帧图像中的噪声的估计方法主要包括如下步骤:In the prior art, methods for estimating noise in images or video frame images mainly include the following steps:

步骤1、对待估计图像或者视频帧图像进行不重叠分块处理;Step 1, performing non-overlapping block processing on the image to be estimated or the video frame image;

步骤2、对步骤1中分块处理得到的各图像块分别进行方差计算;Step 2, performing variance calculation on each image block obtained by block processing in step 1;

步骤3、从步骤2计算出的各方差中,选取多个值较小方差进行噪声估计。Step 3. From the variances calculated in step 2, select multiple values with smaller variances for noise estimation.

但是,由于图像或者视频帧图像中的目标存在边缘,而目标的边缘影响了图像的平滑性,因为平滑块的像素值波动可认为是由噪声引起的,因此,通过估计平滑块的方差,就可以获得原始图像中噪声方差的估计,也就是说,经过步骤1的分块处理之后,不存在目标边缘的图像块会比较平滑,计算出的方差与噪声方差很相近,而存在目标边缘的图像块不够平滑,导致计算出的方差与噪声方差相比相差很大,对于低噪声图像,由于噪声较小,边缘块估计出的噪声方差会大于实际噪声方差,导致过高地估计其噪声,对于高噪声图像,由于噪声很大,边缘块估计出的噪声会小于实际噪声方差,导致过低地估计其噪声,使得噪声估计不准确。However, since there are edges in the target in the image or video frame image, and the target's edge affects the smoothness of the image, because the pixel value fluctuation of the smooth block can be considered to be caused by noise, therefore, by estimating the variance of the smooth block, it is An estimate of the noise variance in the original image can be obtained, that is, after the block processing in step 1, the image block without the target edge will be relatively smooth, and the calculated variance is very similar to the noise variance, while the image with the target edge The block is not smooth enough, resulting in a large difference between the calculated variance and the noise variance. For low-noise images, due to the small noise, the estimated noise variance of the edge block will be greater than the actual noise variance, resulting in overestimation of its noise. For high For noisy images, due to the large noise, the estimated noise of the edge block will be smaller than the actual noise variance, resulting in an underestimation of its noise, making the noise estimation inaccurate.

发明内容Contents of the invention

本发明实施例提供了图像的噪声估计方法、视频图像去噪方法及装置,用以解决现有技术中对图像的噪声估计不准确的问题。Embodiments of the present invention provide an image noise estimation method, a video image denoising method and device, to solve the problem of inaccurate image noise estimation in the prior art.

基于上述问题,本发明实施例提供的一种图像的噪声估计方法,包括:Based on the above problems, an image noise estimation method provided by an embodiment of the present invention includes:

确定待处理图像的两个方向的拉普拉斯算子,其中,所述两个方向的拉普拉斯算子满足如下条件:对所述两个方向的拉普拉斯算子进行预设线性运算得到的差异算子中,各行数值之和以及各列数值之和均为零,且位于中心的数值的绝对值大于位于非中心位置的数值的绝对值;Determining Laplacian operators in two directions of the image to be processed, wherein the Laplacian operators in the two directions satisfy the following conditions: preset the Laplacian operators in the two directions In the difference operator obtained by linear operation, the sum of the values of each row and the sum of the values of each column are both zero, and the absolute value of the value at the center is greater than the absolute value of the value at the non-center position;

利用所述差异算子对所述待处理图像进行平滑处理,得到平滑图像;smoothing the image to be processed by using the difference operator to obtain a smooth image;

基于第一预设尺寸,将所述平滑图像进行不重叠分块,得到各图像块;Based on the first preset size, the smooth image is divided into non-overlapping blocks to obtain each image block;

分别确定所述各图像块的方差,并Determine the variance of each image block respectively, and

从各方差中,按照方差值从小到大的顺序确定预设数量个方差,并将所述预设数量个方差的加权平均值确定为估计的所述待处理图像的噪声方差。From each variance, a preset number of variances are determined in ascending order of variance values, and a weighted average of the preset number of variances is determined as the estimated noise variance of the image to be processed.

本发明实施例提供的一种视频图像去噪方法,包括:A video image denoising method provided by an embodiment of the present invention includes:

基于上述一种图像的噪声估计方法,确定待处理视频的当前帧图像的噪声方差,并根据所述噪声方差确定运动检测阈值;Based on the above-mentioned noise estimation method of an image, determine the noise variance of the current frame image of the video to be processed, and determine the motion detection threshold according to the noise variance;

基于第二预设尺寸,分别将所述当前帧图像和前一帧图像进行不重叠分块,得到所述当前帧图像的各图像块和所述前一帧图像的各图像块;Based on the second preset size, the current frame image and the previous frame image are respectively divided into non-overlapping blocks to obtain each image block of the current frame image and each image block of the previous frame image;

针对所述当前帧图像的每个图像块,确定该图像块与所述前一帧图像中与该图像块位置对应的图像块之间的相似度;并For each image block of the current frame image, determine the similarity between the image block and the image block corresponding to the position of the image block in the previous frame image; and

将所述相似度与所述运动检测阈值进行比较,并根据比较结果确定该图像块为运动图像块或者静止图像块;Comparing the similarity with the motion detection threshold, and determining that the image block is a moving image block or a still image block according to the comparison result;

分别对所述当前帧图像中的运动图像块和静止图像块进行去噪处理。Denoising processing is performed on moving image blocks and still image blocks in the current frame image respectively.

本发明实施例提供的一种图像的噪声估计装置,包括:An image noise estimation device provided by an embodiment of the present invention includes:

算子确定模块,用于确定待处理图像的两个方向的拉普拉斯算子,其中,所述两个方向的拉普拉斯算子满足如下条件:对所述两个方向的拉普拉斯算子进行预设线性运算得到的差异算子中,各行数值之和以及各列数值之和均为零,且位于中心的数值的绝对值大于位于非中心位置的数值的绝对值;An operator determination module, configured to determine Laplacian operators in two directions of the image to be processed, wherein the Laplacian operators in the two directions satisfy the following conditions: the Laplacian operators in the two directions In the difference operator obtained by the preset linear operation of the Las operator, the sum of the values of each row and the sum of the values of each column are both zero, and the absolute value of the value at the center is greater than the absolute value of the value at the non-center position;

平滑处理模块,用于利用所述差异算子对所述待处理图像进行平滑处理,得到平滑图像;a smoothing processing module, configured to use the difference operator to perform smoothing processing on the image to be processed to obtain a smoothed image;

第一分块模块,用于基于第一预设尺寸,将所述平滑图像进行不重叠分块,得到各图像块;The first block module is used to perform non-overlapping block on the smooth image based on the first preset size to obtain each image block;

方差确定模块,用于分别确定所述各图像块的方差,并a variance determination module, configured to determine the variance of each image block respectively, and

噪声估计模块,用于从各方差中,按照方差值从小到大的顺序确定预设数量个方差,并将所述预设数量个方差的加权平均值确定为估计的所述待处理图像的噪声方差。A noise estimation module, configured to determine a preset number of variances from the variances in ascending order of variance values, and determine a weighted average of the preset number of variances as the estimated value of the image to be processed noise variance.

本发明实施例提供的一种视频图像去噪装置,包括:A video image denoising device provided in an embodiment of the present invention includes:

阈值确定模块,用于通过上述图像的噪声估计装置,确定待处理视频的当前帧图像的噪声方差,并根据所述噪声方差确定运动检测阈值;A threshold determination module, configured to determine the noise variance of the current frame image of the video to be processed through the above image noise estimation device, and determine the motion detection threshold according to the noise variance;

第二分块模块,用于基于第二预设尺寸,分别将所述当前帧图像和前一帧图像进行不重叠分块,得到所述当前帧图像的各图像块和所述前一帧图像的各图像块;The second block module is used to divide the current frame image and the previous frame image into non-overlapping blocks based on a second preset size, so as to obtain each image block of the current frame image and the previous frame image Each image block of

相似度确定模块,用于针对所述当前帧图像的每个图像块,确定该图像块与所述前一帧图像中与该图像块位置对应的图像块之间的相似度;并A similarity determination module, configured to, for each image block of the current frame image, determine the similarity between the image block and the image block corresponding to the position of the image block in the previous frame image; and

运动检测模块,用于将所述相似度与所述运动检测阈值进行比较,并根据比较结果确定该图像块为运动图像块或者静止图像块;A motion detection module, configured to compare the similarity with the motion detection threshold, and determine that the image block is a motion image block or a still image block according to the comparison result;

去噪处理模块,用于分别对所述当前帧图像中的运动图像块和静止图像块进行去噪处理。The denoising processing module is used to respectively perform denoising processing on the moving image block and the still image block in the current frame image.

本发明实施例的有益效果包括:The beneficial effects of the embodiments of the present invention include:

本发明实施例提供的图像的噪声估计方法、视频图像去噪方法及装置,确定待处理图像的两个方向的拉普拉斯算子,其中,两个方向的拉普拉斯算子满足如下条件:对两个方向的拉普拉斯算子进行预设线性运算得到的差异算子中,各行数值之和以及各列数值之和均为零,且位于中心的数值的绝对值大于位于非中心位置的数值的绝对值;利用差异算子对待处理图像进行平滑处理,得到平滑图像;基于第一预设尺寸,将平滑图像进行不重叠分块,得到各图像块;分别确定各图像块的方差,并从各方差中,按照方差值从小到大的顺序确定预设数量个方差,并将预设数量个方差的加权平均值确定为估计的待处理图像的噪声方差。在对待处理图像进行分块之前,对图像进行平滑处理,使得图像中的目标边缘变得平滑,那么,划分的各图像块中的图像变得平滑,各图像块计算的方差可以较为正常,利用方差进行图像的噪声估计更加准确;The image noise estimation method, video image denoising method and device provided by the embodiments of the present invention determine the Laplacian operators in two directions of the image to be processed, wherein the Laplacian operators in the two directions satisfy the following Condition: In the difference operator obtained by performing preset linear operations on Laplacian operators in two directions, the sum of the values in each row and the sum of the values in each column are both zero, and the absolute value of the value at the center is greater than that at the non- The absolute value of the numerical value of the center position; use the difference operator to smooth the image to be processed to obtain a smooth image; based on the first preset size, the smooth image is divided into non-overlapping blocks to obtain each image block; determine each image block respectively variance, and from each variance, determine a preset number of variances in ascending order of variance values, and determine the weighted average of the preset number of variances as the estimated noise variance of the image to be processed. Before the image to be processed is divided into blocks, the image is smoothed so that the target edge in the image becomes smooth, then the image in each divided image block becomes smooth, and the variance calculated by each image block can be relatively normal. Using It is more accurate to estimate the noise of the image with the variance;

进一步地,利用上述噪声估计方法确定待处理视频的当前帧图像的噪声方差,并根据噪声方差确定运动检测阈值;基于第二预设尺寸,分别将当前帧图像和前一帧图像进行不重叠分块,得到当前帧图像的各图像块和前一帧图像的各图像块;针对当前帧图像的每个图像块,确定该图像块与前一帧图像中与该图像块位置对应的图像块之间的相似度;并将相似度与运动检测阈值进行比较,并根据比较结果确定该图像块为运动图像块或者静止图像块;分别对当前帧图像中的运动图像块和静止图像块进行去噪处理。在上述对视频图像进行去噪处理的过程中,基于准确的噪声估计方法对图像进行噪声估计,再进行去噪处理,使得去噪处理效果更好。Further, the noise variance of the current frame image of the video to be processed is determined using the above noise estimation method, and the motion detection threshold is determined according to the noise variance; based on the second preset size, the current frame image and the previous frame image are respectively non-overlapping blocks to obtain each image block of the current frame image and each image block of the previous frame image; for each image block of the current frame image, determine the difference between the image block and the image block corresponding to the image block position in the previous frame image The similarity between them; and the similarity is compared with the motion detection threshold, and according to the comparison result, it is determined that the image block is a moving image block or a still image block; respectively denoising the moving image block and the still image block in the current frame image deal with. In the above process of denoising the video image, noise estimation is performed on the image based on an accurate noise estimation method, and then denoising processing is performed, so that the denoising processing effect is better.

附图说明Description of drawings

图1为本发明实施例提供的一种图像的噪声估计方法的流程图;FIG. 1 is a flowchart of an image noise estimation method provided by an embodiment of the present invention;

图2为本发明实施例提供的一种视频图像去噪方法的流程图;FIG. 2 is a flow chart of a video image denoising method provided by an embodiment of the present invention;

图3为本发明实施例1提供的一种图像的噪声估计方法的流程图;FIG. 3 is a flowchart of an image noise estimation method provided in Embodiment 1 of the present invention;

图4a-图4c为本发明实施例1提供的拉普拉斯算子模板举例;Figure 4a-Figure 4c is an example of the Laplacian operator template provided by Embodiment 1 of the present invention;

图5为本发明实施例2提供的一种视频图像去噪方法的流程图;FIG. 5 is a flow chart of a video image denoising method provided in Embodiment 2 of the present invention;

图6为本发明实施例3提供的一种视频图像去噪方法的流程图;FIG. 6 is a flowchart of a video image denoising method provided by Embodiment 3 of the present invention;

图7为本发明实施例提供的一种图像的噪声估计装置的结构示意图;FIG. 7 is a schematic structural diagram of an image noise estimation device provided by an embodiment of the present invention;

图8为本发明实施例提供的一种视频图像去噪装置的结构示意图。FIG. 8 is a schematic structural diagram of a video image denoising device provided by an embodiment of the present invention.

具体实施方式detailed description

本发明实施例提供了图像的噪声估计方法、视频图像去噪方法及装置,以下结合说明书附图对本发明的优选实施例进行说明,应当理解,此处所描述的优选实施例仅用于说明和解释本发明,并不用于限定本发明。并且在不冲突的情况下,本申请中的实施例及实施例中的特征可以相互组合。The embodiment of the present invention provides an image noise estimation method, a video image denoising method and device. The preferred embodiments of the present invention will be described below in conjunction with the accompanying drawings. It should be understood that the preferred embodiments described here are only for illustration and explanation The present invention is not intended to limit the present invention. And in the case of no conflict, the embodiments in the present application and the features in the embodiments can be combined with each other.

本发明实施例提供一种图像的噪声估计方法,如图1所示,包括:An embodiment of the present invention provides an image noise estimation method, as shown in FIG. 1 , including:

S101、确定待处理图像的两个方向的拉普拉斯算子,其中,该两个方向的拉普拉斯算子满足如下条件:对该两个方向的拉普拉斯算子进行预设线性运算得到的差异算子中,各行数值之和以及各列数值之和均为零,且位于中心的数值的绝对值大于位于非中心位置的数值的绝对值。S101. Determine the Laplacian operators in two directions of the image to be processed, wherein the Laplacian operators in the two directions meet the following conditions: preset the Laplacian operators in the two directions In the difference operator obtained by linear operation, the sum of the values in each row and the sum of the values in each column are both zero, and the absolute value of the value at the center is greater than the absolute value of the value at the non-center position.

S102、利用S101中确定的差异算子对待处理图像进行平滑处理,得到平滑图像。S102. Using the difference operator determined in S101 to perform smoothing processing on the image to be processed to obtain a smoothed image.

S103、基于第一预设尺寸,将S102中得到的平滑图像进行不重叠分块,得到各图像块。S103. Based on the first preset size, divide the smooth image obtained in S102 into non-overlapping blocks to obtain image blocks.

S104、分别确定S103中确定的各图像块的方差。S104. Determine respectively the variance of each image block determined in S103.

S105、从S104中确定的各方差中,按照方差值从小到大的顺序确定预设数量个方差。S105. From the variances determined in S104, determine a preset number of variances in ascending order of variance values.

S106、并将S105中确定的预设数量个方差的加权平均值确定为估计的待处理图像的噪声方差。S106. Determine the weighted average of the preset number of variances determined in S105 as the estimated noise variance of the image to be processed.

相应的,基于上述图像的噪声估计方法,本发明实施例还提供一种视频图像去噪方法,如图2所示,包括:Correspondingly, based on the above image noise estimation method, an embodiment of the present invention also provides a video image denoising method, as shown in FIG. 2 , including:

S201、基于上述图像的噪声估计方法,确定待处理视频的当前帧图像的噪声方差,并根据噪声方差确定运动检测阈值。S201. Determine the noise variance of the current frame image of the video to be processed based on the above image noise estimation method, and determine a motion detection threshold according to the noise variance.

S202、基于第二预设尺寸,分别将当前帧图像和前一帧图像进行不重叠分块,得到当前帧图像的各图像块和前一帧图像的各图像块。S202. Based on the second preset size, respectively divide the current frame image and the previous frame image into non-overlapping blocks to obtain each image block of the current frame image and each image block of the previous frame image.

S203、针对当前帧图像的每个图像块,确定该图像块与前一帧图像中与该图像块位置对应的图像块之间的相似度。S203. For each image block of the current frame image, determine the similarity between the image block and the image block corresponding to the position of the image block in the previous frame image.

S204、将S203中确定的相似度与S201中确定的运动检测阈值进行比较,并根据比较结果确定该图像块为运动图像块或者静止图像块。S204. Compare the similarity determined in S203 with the motion detection threshold determined in S201, and determine that the image block is a moving image block or a still image block according to the comparison result.

S205、分别对当前帧图像中的运动图像块和静止图像块进行去噪处理。S205. Perform denoising processing on moving image blocks and still image blocks in the current frame image respectively.

进一步地,上述步骤S201的执行与步骤S202-步骤S203没有先后顺序。Further, the execution of step S201 and step S202-step S203 are not performed sequentially.

下面结合附图,用具体实施例对本发明提供的方法及相关设备进行详细描述。The method and related equipment provided by the present invention will be described in detail below with specific embodiments in conjunction with the accompanying drawings.

实施例1:Example 1:

本发明实施例1中,提供了一种图像的噪声估计方法,如图3所示,具体包括如下步骤:In Embodiment 1 of the present invention, an image noise estimation method is provided, as shown in FIG. 3 , which specifically includes the following steps:

S301、确定待处理图像的两个方向的拉普拉斯算子,其中,该两个方向的拉普拉斯算子满足如下条件:对该两个方向的拉普拉斯算子进行预设线性运算得到的差异算子中,各行数值之和以及各列数值之和均为零,且位于中心的数值的绝对值大于位于非中心位置的数值的绝对值。S301. Determine the Laplacian operators in two directions of the image to be processed, wherein the Laplacian operators in the two directions meet the following conditions: preset the Laplacian operators in the two directions In the difference operator obtained by linear operation, the sum of the values in each row and the sum of the values in each column are both zero, and the absolute value of the value at the center is greater than the absolute value of the value at the non-center position.

进一步地,假设确定的待处理图像的两个方向的拉普拉斯算子进行如式(1)的线性运算,得到差异算子dL:Further, assuming that the determined Laplacian operators in two directions of the image to be processed perform the linear operation as in formula (1), the difference operator dL is obtained:

dL=aL1-bL2;(1)dL=aL 1 -bL 2 ; (1)

其中,L1和L2分别表示确定的两个方向的拉普拉斯算子;且差异算子dL满足各行数值之和以及各列数值之和均为零,且位于中心的数值的绝对值大于位于其他位置的数值的绝对值。Among them, L 1 and L 2 respectively represent the Laplacian operators in the two determined directions; and the difference operator dL satisfies the sum of the values of each row and the sum of the values of each column are zero, and the absolute value of the value at the center Absolute value that is greater than a numeric value elsewhere.

下面举例说明:The following example illustrates:

假设a和b分别取2和1,图4a-图4c分别为L1、L2和dL对应的模板。如图4a-图4c所示,由于dL对应的模板中,各行数值之和均为0,各列数值之和均为0,且位于中间的数值的绝对值4,大于位于其他位置的数值的绝对值,因此,L1、L2能够满足条件,可以将L1、L2确定为选择的两个方向的拉普拉斯算子。Assuming that a and b are 2 and 1 respectively, Figure 4a-Figure 4c are templates corresponding to L 1 , L 2 and dL respectively. As shown in Figure 4a-Figure 4c, in the template corresponding to dL, the sum of the values of each row is 0, the sum of the values of each column is 0, and the absolute value of the value in the middle is 4, which is greater than the value of the value in other positions Therefore, L 1 and L 2 can satisfy the conditions, and L 1 and L 2 can be determined as Laplacian operators in the two selected directions.

进一步地,现有技术中可以利用拉普拉斯算子的零交叉性质进行图像中的目标边缘定位,而本实施例中,由于差异算子的性质,可以使用差异算子对待处理图像进行平滑处理。Further, in the prior art, the zero-crossing property of the Laplacian operator can be used to locate the target edge in the image, but in this embodiment, due to the nature of the difference operator, the difference operator can be used to smooth the image to be processed deal with.

S302、利用S301中确定的差异算子采用如下公式对所述待处理图像进行平滑处理:S302. Use the difference operator determined in S301 to perform smoothing processing on the image to be processed by using the following formula:

nno (( xx ,, ythe y )) == ΣΣ jj == -- ll ll ΣΣ ii == -- ll ll dLL (( xx ++ ii ,, ythe y ++ jj )) vv (( xx ++ ii ,, ythe y ++ jj )) ;; -- -- -- (( 22 ))

其中,(x,y)表示所述待处理图像中的像素点;v(x,y)表示对像素点(x,y)进行平滑处理之前,像素点(x,y)的灰度值;n(x,y)表示对像素点(x,y)进行平滑处理之后,像素点(x,y)的灰度值;l的值为拉普拉斯算子模板的行数或者列数减1除以2。Wherein, (x, y) represents the pixel in the image to be processed; v(x, y) represents the gray value of the pixel (x, y) before smoothing the pixel (x, y); n(x, y) represents the gray value of the pixel point (x, y) after smoothing the pixel point (x, y); the value of l is the number of rows or columns of the Laplacian operator template minus 1 divided by 2.

进一步地,(2)式中以(x,y)为中心像素点的(2l+1)×(2l+1)图像块与拉普拉斯算子模板相对应,对该中心像素点(x,y)进行平滑处理之后的灰度值n(x,y)为使用拉普拉斯模板dL(x+i,y+j)中的每个值分别与图像块中与该值位置对应的像素点的灰度值相乘,再将得到的各乘积相加的得到的值来代替中心像素(x,y)的灰度值即为处理后的灰度值n(x,y)。Furthermore, the (2l+1)×(2l+1) image block with (x, y) as the center pixel in formula (2) corresponds to the Laplacian template, and the center pixel (x , y) after smoothing, the gray value n(x, y) is each value in the Laplacian template dL(x+i, y+j) corresponding to the position of the value in the image block Multiply the gray values of the pixels, and then add the obtained values to replace the gray value of the central pixel (x, y), which is the processed gray value n(x, y).

S303、基于第一预设尺寸,将S302中得到的平滑图像进行不重叠分块,得到各图像块。S303. Based on the first preset size, divide the smooth image obtained in S302 into non-overlapping blocks to obtain image blocks.

S304、分别确定S303中确定的各图像块的方差。S304. Determine respectively the variance of each image block determined in S303.

S305、从S304中确定的各方差中,按照方差值从小到大的顺序确定预设数量个方差。S305. From the variances determined in S304, determine a preset number of variances in ascending order of variance values.

较佳地,本步骤中可以将各方差由小到大进行排序,并确定前6%的方差。Preferably, in this step, the variances can be sorted from small to large, and the top 6% variances can be determined.

S306、并将S305中确定的预设数量个方差的加权平均值确定为估计的待处理图像的噪声方差。S306. Determine the weighted average of the preset number of variances determined in S305 as the estimated noise variance of the image to be processed.

本步骤中,确定的各个方差的权值可以相等,即将预设数量个方差的平均值确定为估计的待处理图像的噪声方差。In this step, the weights of the determined variances may be equal, that is, the average value of a preset number of variances is determined as the estimated noise variance of the image to be processed.

实施例2:Example 2:

本发明实施例2中,提供了一种视频图像去噪方法,如图5所示,具体包括如下步骤:In Embodiment 2 of the present invention, a video image denoising method is provided, as shown in FIG. 5 , which specifically includes the following steps:

S501、基于上述图像的噪声估计方法,确定待处理视频的当前帧图像的噪声方差,并根据噪声方差确定运动检测阈值。S501. Determine the noise variance of the current frame image of the video to be processed based on the above image noise estimation method, and determine a motion detection threshold according to the noise variance.

进一步地,可以根据实际需要使用预设系数与噪声方差相乘,并将乘积确定为运动检测阈值。Further, a preset coefficient may be used to multiply the noise variance according to actual needs, and the product may be determined as the motion detection threshold.

进一步地,本步骤与步骤S502-S503没有严格的执行顺序,也就是说,只要在步骤S504之前确定出运动检测阈值即可。Further, there is no strict execution sequence between this step and steps S502-S503, that is, as long as the motion detection threshold is determined before step S504.

S502、基于第二预设尺寸,分别将当前帧图像和前一帧图像进行不重叠分块,得到当前帧图像的各图像块和前一帧图像的各图像块。S502. Based on the second preset size, respectively divide the current frame image and the previous frame image into non-overlapping blocks to obtain each image block of the current frame image and each image block of the previous frame image.

S503、针对当前帧图像的每个图像块,确定该图像块与前一帧图像中与该图像块位置对应的图像块之间的相似度。S503. For each image block of the current frame image, determine the similarity between the image block and the image block corresponding to the position of the image block in the previous frame image.

S504、将S503中确定的相似度与S501中确定的运动检测阈值进行比较,并根据比较结果确定该图像块为运动图像块或者静止图像块。S504. Compare the similarity determined in S503 with the motion detection threshold determined in S501, and determine that the image block is a moving image block or a still image block according to the comparison result.

S505、分别对当前帧图像中的运动图像区域和静止图像区域进行去噪处理。S505. Perform denoising processing on the moving image area and the still image area in the current frame image respectively.

S506、确定当前帧图像中每个像素点的梯度值。S506. Determine the gradient value of each pixel in the current frame image.

进一步地,本步骤与步骤S501-步骤S505的执行没有严格的先后顺序,也就是说,只要在执行步骤S507之前,确定出当前帧图像中每个像素点的梯度值即可。Further, there is no strict sequence between this step and the execution of step S501-step S505, that is, as long as the gradient value of each pixel in the current frame image is determined before step S507 is executed.

进一步地,也可以将本步骤修改为:确定当前帧图像中的运动图像块中每个像素点的梯度值。Further, this step can also be modified to: determine the gradient value of each pixel in the moving image block in the current frame image.

S507、从运动图像块中的各像素点中,确定梯度值大于预设梯度阈值的像素点。S507. From each pixel in the moving image block, determine a pixel whose gradient value is greater than a preset gradient threshold.

进一步地,本步骤中,梯度值大于预设梯度阈值的像素点可以认为是作为运动目标的边缘的像素点。Further, in this step, the pixel points whose gradient values are greater than the preset gradient threshold may be regarded as the pixels of the edge of the moving object.

S508、针对确定的像素点中的每个像素点,将对该像素点进行去噪处理之后得到的灰度值,与预设比例的对该像素点进行去噪处理之前的灰度值的和,进一步确定为对该像素点进行去噪处理之后的灰度值。S508. For each pixel in the determined pixels, the sum of the gray value obtained after denoising the pixel and the gray value before denoising the pixel in a preset ratio , which is further determined as the gray value of the pixel after denoising processing.

进一步地,对于当前帧图像中的运动图像块中和静止图像块中,作为非运动目标边缘的像素点,可以将去噪处理之后,得到的像素点的灰度值直接作为去噪处理之后该像素点的灰度值,而对于作为运动目标边缘的像素点,还需要按照本步骤进行进一步的处理。经过进一步地处理,能够使得图像中运动目标的细小边缘完整地保留,使得处理后的图像不会出现拖影现象。Further, for the pixels in the moving image block and the still image block in the current frame image as the edge of the non-moving object, the gray value of the pixel obtained after the denoising process can be directly used as the pixel point after the denoising process. The gray value of the pixel point, and for the pixel point as the edge of the moving target, it needs to be further processed according to this step. After further processing, the fine edge of the moving object in the image can be completely preserved, so that no smear phenomenon will appear in the processed image.

本步骤中,假设作为运动目标边缘的像素点,进行去噪处理之前灰度值为Y,进行去噪处理之后灰度值为Y',预设比例为α,则该像素点进行去噪处理之后,最终灰度值可以为:Y'+αY。In this step, assuming that the pixel as the edge of the moving target has a gray value of Y before the denoising process, and a gray value of Y' after the denoising process, and the preset ratio is α, then the pixel is denoised. After that, the final gray value can be: Y'+αY.

实施例3:Example 3:

本发明实施例3中,提供了一种视频图像去噪方法,如图6所示,具体包括如下步骤:In Embodiment 3 of the present invention, a video image denoising method is provided, as shown in FIG. 6 , which specifically includes the following steps:

S601、基于上述图像的噪声估计方法,确定待处理视频的当前帧图像的噪声方差,并根据噪声方差确定运动检测阈值。S601. Determine the noise variance of the current frame image of the video to be processed based on the above image noise estimation method, and determine a motion detection threshold according to the noise variance.

进一步地,本步骤与步骤S602-S604没有严格的执行顺序,也就是说,只要在步骤S605之前确定出运动检测阈值即可。Further, there is no strict execution sequence between this step and steps S602-S604, that is, as long as the motion detection threshold is determined before step S605.

S602、基于第二预设尺寸,分别将当前帧图像和前一帧图像进行不重叠分块,得到当前帧图像的各图像块和前一帧图像的各图像块。S602. Based on the second preset size, respectively divide the current frame image and the previous frame image into non-overlapping blocks to obtain each image block of the current frame image and each image block of the previous frame image.

S603、针对当前帧图像和前一帧图像的每个图像块,确定该图像块所在帧中,以该图像块为中心的s×s个图像块构成的图像区域。S603. For each image block of the current frame image and the previous frame image, determine an image area composed of s×s image blocks centered on the image block in the frame where the image block is located.

S604、针对当前帧图像的每个图像块,采用如下公式确定该图像块与前一帧图像中与该图像块位置对应的图像块之间的相似度D(i,j,k):S604. For each image block of the current frame image, the following formula is used to determine the similarity D(i, j, k) between the image block and the image block corresponding to the position of the image block in the previous frame image:

DD. (( ii ,, jj ,, kk )) == ΣΣ mm == ii -- hh ii ++ hh ΣΣ nno == jj -- hh jj ++ hh [[ SS BB (( mm ,, nno ,, kk )) -- SS BB (( mm ,, nno ,, kk -- 11 )) ]] 22 -- -- -- (( 33 ))

其中,(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)的图像块中所有像素点灰度值的和。Among them, (i, j) represents the center coordinates of the image block; the value of h is (s-1)/2; k and k-1 represent the frame number; (m, n) represents the image centered on the image block The coordinates of the image blocks other than the image block in the area relative to the image block; S B (m, n, k) represents the grayscale of all pixels in the image block whose center coordinates are (m, n) in the kth frame The sum of values; S B (m,n,k-1) represents the sum of the gray values of all pixels in the image block whose center coordinates are (m,n) in the k-1th frame.

进一步地,假设s的值为3,也就是说,针对当前帧图像和前一帧图像的每个图像块,确定该图像块以及该图像块所在帧中,以该图像块为中心的3×3个图像块构成的图像区域,则以该图像块为中心的图像区域中,包括9个图像块。那么式(3)中,分别计算当前帧的9个图像块中每个图像块中所有像素点灰度值的和,以及前一帧的9个图像块中每个图像块中所有像素点灰度值的和,并分别确定当前帧中以及前一帧中各自的9个图像块中,每两个位置对应的图像块中所有像素点灰度值的和的差值的平方,并将得到的9个差值的平方相加,得到的和作为该图像块的相似度。Further, assuming that the value of s is 3, that is, for each image block of the current frame image and the previous frame image, determine the image block and the frame in which the image block is located, and a 3× An image area composed of three image blocks includes nine image blocks in the image area centered on the image block. Then in formula (3), calculate the sum of the gray values of all pixels in each of the nine image blocks of the current frame, and the gray values of all pixels in each of the nine image blocks in the previous frame and determine the square of the difference between the sum of the gray values of all pixels in the image blocks corresponding to each two positions in each of the nine image blocks in the current frame and the previous frame, and will get The squares of the 9 differences are summed, and the sum is taken as the similarity of the image block.

S605、当S604中确定的相似度大于S601中确定的运动检测阈值时,确定当前帧的该图像块相对上一帧中与该图像块位置对应的图像块为运动图像块;S605. When the similarity determined in S604 is greater than the motion detection threshold determined in S601, determine that the image block in the current frame is a moving image block relative to the image block corresponding to the position of the image block in the previous frame;

当S604中确定的相似度不大于S601中确定的运动检测阈值时,确定当前帧的该图像块相对上一帧中与该图像块位置对应的图像块为静止图像块。When the similarity determined in S604 is not greater than the motion detection threshold determined in S601, it is determined that the image block in the current frame is a still image block relative to the image block corresponding to the position of the image block in the previous frame.

进一步地,本步骤中,当D(i,j,k)大于运动检测阈值时,可以将(i,j)表征的图像块中的每个像素点标记为1,否则,标记为0,那么被标记为1的像素点构成的图像块表征运动图像块,被标记为0的像素点构成的图像块表征静止图像块。Further, in this step, when D(i, j, k) is greater than the motion detection threshold, each pixel in the image block represented by (i, j) can be marked as 1, otherwise, marked as 0, then An image block composed of pixels marked as 1 represents a moving image block, and an image block composed of pixels marked 0 represents a still image block.

S606、分别对当前帧图像中的运动图像块和静止图像块进行去噪处理。S606. Perform denoising processing on moving image blocks and still image blocks in the current frame image respectively.

本步骤中,可以采用改进的3D-Rational滤波公式对当前帧图像中的静止图像块进行去噪处理,如式(4)所示:In this step, the improved 3D-Rational filtering formula can be used to denoise the still image block in the current frame image, as shown in formula (4):

ythe y 00 tt == xx 00 tt -- ΣΣ ii ,, jj ∈∈ ΩΩ -- xx ii tt -- 11 ++ 22 xx 00 tt -- xx jj tt -- 11 kk tt (( xx ii tt -- 11 -- xx jj tt )) 22 ++ AA tt -- -- -- (( 44 ))

针对静止图像块中的每个像素点,假设该像素点进行去噪处理之前的灰度值为x0,该像素点进行去噪处理之后的灰度值为y0,t和t-1分别用来标识帧号(不表示幂值),Ω表示与该像素点相邻的九个像素点的索引集合,i和j分别表示Ω中的任一像素点,xi和xj分别表示i像素点和j像素点的灰度值,kt,At均表示时域去噪调节因子。For each pixel in the still image block, it is assumed that the gray value of the pixel before the denoising process is x 0 , and the gray value of the pixel after the denoising process is y 0 , t and t-1 are respectively Used to identify the frame number (does not represent the power value), Ω represents the index set of nine pixels adjacent to the pixel, i and j represent any pixel in Ω, x i and x j represent i respectively The gray value of the pixel and the j pixel, k t , At both represent the time domain denoising adjustment factor.

本步骤中,可以采用改进的快速NLM方法,对当前帧图像中的运动图像块进行去噪处理,如式(5)所示:In this step, the improved fast NLM method can be used to denoise the moving image block in the current frame image, as shown in formula (5):

NLvNLv (( ii )) == ΣΣ jj ∈∈ II ww (( ii ,, jj )) vv (( jj )) -- -- -- (( 55 ))

其中,in,

ZZ (( ii )) == ΣΣ jj ww (( ii ,, jj )) ..

在式(5)中,NLv(i)表示对像素点i进行去噪处理之后像素点i的灰度值,像素点j为与像素点i相邻的像素点,v(j)表示对像素点j进行去噪处理之前像素点j的灰度值,w(i,j)表示利用像素点j进行去噪处理之前像素点j的灰度值确定像素点i进行去噪处理之后像素点i的灰度值时,使用的各个权值。In formula (5), NLv(i) represents the gray value of pixel i after denoising processing on pixel i, pixel j is the pixel point adjacent to pixel i, and v(j) represents the gray value of pixel i Point j is the gray value of pixel j before the denoising process, w(i, j) means that the gray value of pixel j before denoising process is used to determine pixel i after denoising process When the gray value of , the weights used.

式(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。Equation (6) is the calculation formula of weight w(i,j), where ||N(i)-N(j)|| 2 represents the Euclidean distance between pixel i and pixel j, and h is the attenuation Factor, ▽v(i) represents the gradient value of pixel i, ▽v(j) represents the gradient value of pixel j, σ represents another gradient threshold, Indicates the mean value of pixel i, Represents the average value of pixel point j, η1 and η2 represent the first average threshold and the second average threshold respectively, preferably, η1 and η2 can be 0.9 and 1.2 respectively.

基于同一发明构思,本发明实施例还提供了一种图像的噪声估计装置,由于这些装置所解决问题的原理与前述一种图像的噪声估计方法相似,因此该装置的实施可以参见前述方法的实施,重复之处不再赘述。Based on the same inventive concept, an embodiment of the present invention also provides an image noise estimation device. Since the principle of the problem solved by these devices is similar to the aforementioned image noise estimation method, the implementation of the device can refer to the implementation of the aforementioned method , the repetitions will not be repeated.

本发明实施例提供的一种图像的噪声估计装置,如图7所示,包括如下模块:An image noise estimation device provided by an embodiment of the present invention, as shown in FIG. 7 , includes the following modules:

算子确定模块701,用于确定待处理图像的两个方向的拉普拉斯算子,其中,所述两个方向的拉普拉斯算子满足如下条件:对所述两个方向的拉普拉斯算子进行预设线性运算得到的差异算子中,各行数值之和以及各列数值之和均为零,且位于中心的数值的绝对值大于位于非中心位置的数值的绝对值;The operator determination module 701 is configured to determine Laplacian operators in two directions of the image to be processed, wherein the Laplacian operators in the two directions satisfy the following conditions: the Laplacian operators in the two directions In the difference operator obtained by the preset linear operation of the Placian operator, the sum of the values of each row and the sum of the values of each column are both zero, and the absolute value of the value at the center is greater than the absolute value of the value at the non-center position;

平滑处理模块702,用于利用所述差异算子对所述待处理图像进行平滑处理,得到平滑图像;A smoothing processing module 702, configured to use the difference operator to perform smoothing processing on the image to be processed to obtain a smoothed image;

第一分块模块703,用于基于第一预设尺寸,将所述平滑图像进行不重叠分块,得到各图像块;The first block module 703 is configured to perform non-overlapping block on the smooth image based on a first preset size to obtain each image block;

方差确定模块704,用于分别确定所述各图像块的方差,并A variance determining module 704, configured to determine the variance of each image block respectively, and

噪声估计模块705,用于从各方差中,按照方差值从小到大的顺序确定预设数量个方差,并将所述预设数量个方差的加权平均值确定为估计的所述待处理图像的噪声方差。The noise estimation module 705 is configured to determine a preset number of variances from the variances in ascending order of variance values, and determine a weighted average of the preset number of variances as the estimated image to be processed noise variance.

进一步地,所述平滑处理模块702,具体用于利用所述差异算子,采用如下公式对所述待处理图像进行平滑处理:Further, the smoothing processing module 702 is specifically configured to use the difference operator to perform smoothing processing on the image to be processed by using the following formula:

nno (( xx ,, ythe y )) == ΣΣ jj == -- ll ll ΣΣ ii == -- ll ll dLL (( xx ++ ii ,, ythe y ++ jj )) vv (( xx ++ ii ,, ythe y ++ jj )) ;;

其中,(x,y)表示所述待处理图像中的像素点;v(x,y)表示对像素点(x,y)进行平滑处理之前,像素点(x,y)的灰度值;n(x,y)表示对像素点(x,y)进行平滑处理之后,像素点(x,y)的灰度值;l的值为所述差异算子模板的行数或者列数减1除以2。Wherein, (x, y) represents the pixel in the image to be processed; v(x, y) represents the gray value of the pixel (x, y) before smoothing the pixel (x, y); n(x, y) represents the gray value of the pixel point (x, y) after smoothing the pixel point (x, y); the value of l is the number of rows or columns of the difference operator template minus 1 Divide by 2.

基于同一发明构思,本发明实施例还提供了一种基于上述图像的噪声估计装置的视频图像去噪装置,由于这些装置所解决问题的原理与前述视频图像去噪方法相似,因此该装置的实施可以参见前述方法的实施,重复之处不再赘述。Based on the same inventive concept, the embodiment of the present invention also provides a video image denoising device based on the above image noise estimation device. Since the principle of the problem solved by these devices is similar to the aforementioned video image denoising method, the implementation of the device Reference may be made to the implementation of the aforementioned methods, and repeated descriptions will not be repeated here.

本发明实施例提供的一种视频图像去噪装置,如图8所示,包括如下模块:A video image denoising device provided in an embodiment of the present invention, as shown in FIG. 8 , includes the following modules:

阈值确定模块801,用于通过上述图像的噪声估计装置,确定待处理视频的当前帧图像的噪声方差,并根据所述噪声方差确定运动检测阈值;The threshold determination module 801 is used to determine the noise variance of the current frame image of the video to be processed through the above image noise estimation device, and determine the motion detection threshold according to the noise variance;

第二分块模块802,用于基于第二预设尺寸,分别将所述当前帧图像和前一帧图像进行不重叠分块,得到所述当前帧图像的各图像块和所述前一帧图像的各图像块;The second block module 802 is configured to perform non-overlapping block on the current frame image and the previous frame image based on the second preset size, so as to obtain each image block of the current frame image and the previous frame each image block of the image;

相似度确定模块803,用于针对所述当前帧图像的每个图像块,确定该图像块与所述前一帧图像中与该图像块位置对应的图像块之间的相似度;并A similarity determining module 803, configured to, for each image block of the current frame image, determine the similarity between the image block and the image block corresponding to the position of the image block in the previous frame image; and

运动检测模块804,用于将所述相似度与所述运动检测阈值进行比较,并根据比较结果确定该图像块为运动图像块或者静止图像块;A motion detection module 804, configured to compare the similarity with the motion detection threshold, and determine that the image block is a motion image block or a still image block according to the comparison result;

去噪处理模块805,用于分别对所述当前帧图像中的运动图像块和静止图像块进行去噪处理。The denoising processing module 805 is configured to perform denoising processing on moving image blocks and still image blocks in the current frame image respectively.

进一步地,上述视频图像去噪装置,还包括:Further, the above-mentioned video image denoising device also includes:

梯度值确定模块806,用于确定所述当前帧图像中每个像素点的梯度值;Gradient value determination module 806, configured to determine the gradient value of each pixel in the current frame image;

像素点确定模块807,用于在所述去噪处理模块分别对所述当前帧图像中的运动图像块和静止图像块进行去噪处理之后,从所述运动图像块中的各像素点中,确定梯度值大于预设梯度阈值的像素点;The pixel point determining module 807 is configured to, after the denoising processing module respectively performs denoising processing on the moving image block and the still image block in the current frame image, from each pixel in the moving image block, Determining pixels whose gradient values are greater than a preset gradient threshold;

边缘处理模块808,用于针对确定的像素点中的每个像素点,将对该像素点进行去噪处理之后得到的灰度值,与预设比例的对该像素点进行去噪处理之前的灰度值的和,进一步确定为对该像素点进行去噪处理之后的灰度值。The edge processing module 808 is configured to, for each pixel in the determined pixels, compare the gray value obtained after denoising the pixel with the preset ratio of the gray value before denoising the pixel The sum of the gray values is further determined as the gray value after the denoising process is performed on the pixel.

进一步地,所述相似度确定模块803,具体用于针对所述当前帧图像和所述前一帧图像的每个图像块,确定该图像块以及该图像块所在帧中,以该图像块为中心的s×s个图像块构成的图像区域;Further, the similarity determining module 803 is specifically configured to, for each image block of the current frame image and the previous frame image, determine the image block and the frame in which the image block is located, taking the image block as An image area composed of s×s image blocks in the center;

针对所述当前帧的每个图像块,采用如下公式确定该图像块与所述前一帧中与该图像块位置对应的图像块之间的相似度D(i,j,k):For each image block of the current frame, the following formula is used to determine the similarity D (i, j, k) between the image block and the image block corresponding to the position of the image block in the previous frame:

DD. (( ii ,, jj ,, kk )) == ΣΣ mm == ii -- hh ii ++ hh ΣΣ nno == jj -- hh jj ++ hh [[ SS BB (( mm ,, nno ,, kk )) -- SS BB (( mm ,, nno ,, kk -- 11 )) ]] 22

其中,(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)的图像块中所有像素点灰度值的和。Among them, (i, j) represents the center coordinates of the image block; the value of h is (s-1)/2; k and k-1 represent the frame number; (m, n) represents the image centered on the image block The coordinates of the image blocks other than the image block in the area relative to the image block; S B (m, n, k) represents the grayscale of all pixels in the image block whose center coordinates are (m, n) in the kth frame The sum of values; S B (m,n,k-1) represents the sum of the gray values of all pixels in the image block whose center coordinates are (m,n) in the k-1th frame.

进一步地,所述运动检测模块804,具体用于当所述相似度大于所述运动检测阈值时,确定当前帧的该图像块相对上一帧中与该图像块位置对应的图像块为运动图像块;当所述相似度不大于所述运动检测阈值时,确定当前帧的该图像块相对上一帧中与该图像块位置对应的图像块为静止图像块。Further, the motion detection module 804 is specifically configured to determine that the image block in the current frame is a moving image relative to the image block corresponding to the image block position in the previous frame when the similarity is greater than the motion detection threshold block; when the similarity is not greater than the motion detection threshold, determine that the image block in the current frame is a still image block relative to the image block corresponding to the position of the image block in the previous frame.

上述各单元的功能可对应于图1至图3、图5-图6所示流程中的相应处理步骤,在此不再赘述。The functions of the above-mentioned units may correspond to corresponding processing steps in the flow charts shown in FIG. 1 to FIG. 3 and FIG. 5-FIG.

本发明实施例提供的图像的噪声估计方法、视频图像去噪方法及装置,确定待处理图像的两个方向的拉普拉斯算子,其中,两个方向的拉普拉斯算子满足如下条件:对两个方向的拉普拉斯算子进行预设线性运算得到的差异算子中,各行数值之和以及各列数值之和均为零,且位于中心的数值的绝对值大于位于非中心位置的数值的绝对值;利用差异算子对待处理图像进行平滑处理,得到平滑图像;基于第一预设尺寸,将平滑图像进行不重叠分块,得到各图像块;分别确定各图像块的方差,并从各方差中,按照方差值从小到大的顺序确定预设数量个方差,并将预设数量个方差的加权平均值确定为估计的待处理图像的噪声方差。在对待处理图像进行分块之前,对图像进行平滑处理,使得图像中的目标边缘变得平滑,那么,划分的各图像块中的图像变得平滑,各图像块计算的方差可以较为正常,利用方差进行图像的噪声估计更加准确;The image noise estimation method, video image denoising method and device provided by the embodiments of the present invention determine the Laplacian operators in two directions of the image to be processed, wherein the Laplacian operators in the two directions satisfy the following Condition: In the difference operator obtained by performing preset linear operations on Laplacian operators in two directions, the sum of the values in each row and the sum of the values in each column are both zero, and the absolute value of the value at the center is greater than that at the non- The absolute value of the numerical value of the center position; use the difference operator to smooth the image to be processed to obtain a smooth image; based on the first preset size, the smooth image is divided into non-overlapping blocks to obtain each image block; determine each image block respectively variance, and from each variance, determine a preset number of variances in ascending order of variance values, and determine the weighted average of the preset number of variances as the estimated noise variance of the image to be processed. Before the image to be processed is divided into blocks, the image is smoothed so that the target edge in the image becomes smooth, then the image in each divided image block becomes smooth, and the variance calculated by each image block can be relatively normal. Using It is more accurate to estimate the noise of the image with the variance;

进一步地,利用上述噪声估计方法确定待处理视频的当前帧图像的噪声方差,并根据噪声方差确定运动检测阈值;基于第二预设尺寸,分别将当前帧图像和前一帧图像进行不重叠分块,得到当前帧图像的各图像块和前一帧图像的各图像块;针对当前帧图像的每个图像块,确定该图像块与前一帧图像中与该图像块位置对应的图像块之间的相似度;并将相似度与运动检测阈值进行比较,并根据比较结果确定该图像块为运动图像块或者静止图像块;分别对当前帧图像中的运动图像块和静止图像块进行去噪处理。在上述对视频图像进行去噪处理的过程中,基于准确的噪声估计方法对图像进行噪声估计,再进行去噪处理,使得去噪处理效果更好。Further, the noise variance of the current frame image of the video to be processed is determined using the above noise estimation method, and the motion detection threshold is determined according to the noise variance; based on the second preset size, the current frame image and the previous frame image are respectively non-overlapping blocks to obtain each image block of the current frame image and each image block of the previous frame image; for each image block of the current frame image, determine the difference between the image block and the image block corresponding to the image block position in the previous frame image The similarity between them; and the similarity is compared with the motion detection threshold, and according to the comparison result, it is determined that the image block is a moving image block or a still image block; respectively denoising the moving image block and the still image block in the current frame image deal with. In the above process of denoising the video image, noise estimation is performed on the image based on an accurate noise estimation method, and then denoising processing is performed, so that the denoising processing effect is better.

通过以上的实施方式的描述,本领域的技术人员可以清楚地了解到本发明实施例可以通过硬件实现,也可以借助软件加必要的通用硬件平台的方式来实现。基于这样的理解,本发明实施例的技术方案可以以软件产品的形式体现出来,该软件产品可以存储在一个非易失性存储介质(可以是CD-ROM,U盘,移动硬盘等)中,包括若干指令用以使得一台计算机设备(可以是个人计算机,服务器,或者网络设备等)执行本发明各个实施例所述的方法。Through the above description of the implementation manners, those skilled in the art can clearly understand that the embodiments of the present invention can be implemented by hardware, or by means of software plus a necessary general hardware platform. Based on this understanding, the technical solutions of the embodiments of the present invention can be embodied in the form of software products, which can be stored in a non-volatile storage medium (which can be CD-ROM, U disk, mobile hard disk, etc.), Several instructions are included to make a computer device (which may be a personal computer, a server, or a network device, etc.) execute the methods described in various embodiments of the present invention.

本领域技术人员可以理解附图只是一个优选实施例的示意图,附图中的模块或流程并不一定是实施本发明所必须的。Those skilled in the art can understand that the drawing is only a schematic diagram of a preferred embodiment, and the modules or processes in the drawing are not necessarily necessary for implementing the present invention.

本领域技术人员可以理解实施例中的装置中的模块可以按照实施例描述进行分布于实施例的装置中,也可以进行相应变化位于不同于本实施例的一个或多个装置中。上述实施例的模块可以合并为一个模块,也可以进一步拆分成多个子模块。Those skilled in the art can understand that the modules in the device in the embodiment can be distributed in the device in the embodiment according to the description in the embodiment, or can be located in one or more devices different from the embodiment according to corresponding changes. The modules in the above embodiments can be combined into one module, and can also be further split into multiple sub-modules.

上述本发明实施例序号仅仅为了描述,不代表实施例的优劣。The serial numbers of the above embodiments of the present invention are for description only, and do not represent the advantages and disadvantages of the embodiments.

显然,本领域的技术人员可以对本发明进行各种改动和变型而不脱离本发明的精神和范围。这样,倘若本发明的这些修改和变型属于本发明权利要求及其等同技术的范围之内,则本发明也意图包含这些改动和变型在内。Obviously, those skilled in the art can make various changes and modifications to the present invention without departing from the spirit and scope of the present invention. Thus, if these modifications and variations of the present invention fall within the scope of the claims of the present invention and their equivalent technologies, the present invention also intends to include these modifications and variations.

Claims (12)

1. A method of noise estimation of an image, comprising:
determining Laplacian operators of two directions of an image to be processed, wherein the Laplacian operators of the two directions meet the following conditions: in the difference operator obtained by performing preset linear operation on the laplacian operators in the two directions, the sum of numerical values of each row and the sum of numerical values of each column are all zero, the absolute value of the numerical value positioned in the center is greater than the absolute value of the numerical value positioned in a non-center position, and the preset linear operation formula is as follows: dL ═ aL1-bL2Which isIn, L1And L2Respectively representing Laplacian operators in two determined directions, dL represents the difference operator, a and b are respectively constants, and the number of rows and the number of columns of the Laplacian operators in the two determined directions are equal and are odd numbers;
smoothing the image to be processed by using the difference operator to obtain a smoothed image;
based on a first preset size, carrying out non-overlapping blocking on the smooth image to obtain each image block;
determining the variance of each image block, respectively, and
and determining a preset number of variances according to the sequence of variance values from small to large from each variance, and determining the weighted average value of the preset number of variances as the estimated noise variance of the image to be processed.
2. The method of claim 1, wherein the difference operator is used to smooth the image to be processed using the following formula:
n ( x , y ) = Σ j = - l l Σ i = - l l d L ( x + i , y + j ) v ( x + i , y + j ) ;
wherein, (x, y) represents a pixel point in the image to be processed; v (x, y) represents the gray value of the pixel point (x, y) before the smoothing processing is performed on the pixel point (x, y); n (x, y) represents the gray value of the pixel point (x, y) after the smoothing processing is performed on the pixel point (x, y); the value of l is the number of rows or columns of the Laplace operator template minus 1 divided by 2.
3. A video image denoising method is characterized by comprising the following steps:
the method for estimating noise based on the image according to any one of claims 1-2, determining the noise variance of the current frame image of the video to be processed, and determining the motion detection threshold according to the noise variance;
respectively carrying out non-overlapping blocking on the current frame image and the previous frame image based on a second preset size to obtain each image block of the current frame image and each image block of the previous frame image;
for each image block of the current frame image, determining the similarity between the image block and an image block corresponding to the position of the image block in the previous frame image; and are
Comparing the similarity with the motion detection threshold, and determining the image block to be a motion image block or a static image block according to the comparison result;
and denoising the moving image block and the static image block in the current frame image respectively.
4. The method of claim 3, further comprising:
determining the gradient value of each pixel point in the current frame image;
after the denoising processing is respectively performed on the moving image block and the static image block in the current frame image, the method further includes:
determining pixel points with gradient values larger than a preset gradient threshold value from all pixel points in the moving image blocks;
and aiming at each pixel point in the determined pixel points, further determining the sum of the gray value obtained after the pixel point is subjected to denoising processing and the gray value before the pixel point is subjected to denoising processing in a preset proportion as the gray value after the pixel point is subjected to denoising processing.
5. The method as claimed in claim 3, wherein determining, for each image block of the current frame image, a similarity between the image block and an image block corresponding to the position of the image block in the previous frame image comprises:
determining an image area formed by s multiplied by s image blocks taking the image block as a center in a frame where the image block is located aiming at each image block of the current frame image and the previous frame image;
for each image block of the current frame, determining the similarity D (i, j, k) between the image block and the image block corresponding to the position of the image block in the previous frame by adopting the following formula:
D ( i , j , k ) = Σ m = i - h i + h Σ n = j - h j + h [ S B ( m , n , k ) - S B ( m , n , k - 1 ) ] 2
wherein, (i, j) represents the center coordinates of the image block; h has a value of (s-1)/2; k and k-1 represent frame numbers; (m, n) denotes coordinates of image blocks other than the image block with respect to the image block in an image area centered on the image block; sB(m, n, k) represents the sum of gray values of all pixel points in the image block with the center coordinate (m, n) in the kth frame; sBAnd (m, n, k-1) represents the sum of gray values of all pixel points in the image block with the central coordinate of (m, n) in the (k-1) th frame.
6. The method according to claim 3, wherein comparing the similarity with the motion detection threshold and determining the image block as a moving image block or a stationary image block according to the comparison result specifically comprises:
when the similarity is larger than the motion detection threshold, determining that the image block of the current frame is a motion image block relative to the image block corresponding to the position of the image block in the previous frame;
and when the similarity is not greater than the motion detection threshold, determining that the image block of the current frame is a static image block relative to the image block corresponding to the position of the image block in the previous frame.
7. An apparatus for estimating noise of an image, comprising:
the operator determining module is used for determining laplacian operators in two directions of the image to be processed, wherein the laplacian operators in the two directions meet the following conditions: in the difference operator obtained by performing preset linear operation on the laplacian operators in the two directions, the sum of numerical values of each row and the sum of numerical values of each column are all zero, and the absolute value of the numerical value positioned in the center is greater than that positioned in the non-centerThe absolute value of the numerical value of the heart position is preset with a linear operation formula as follows: dL ═ aL1-bL2Wherein L is1And L2Respectively representing Laplacian operators in two determined directions, dL represents the difference operator, a and b are respectively constants, and the number of rows and the number of columns of the Laplacian operators in the two determined directions are equal and are odd numbers;
the smoothing processing module is used for smoothing the image to be processed by utilizing the difference operator to obtain a smooth image;
the first blocking module is used for carrying out non-overlapping blocking on the smooth image based on a first preset size to obtain each image block;
a variance determining module for determining the variance of each image block respectively and
and the noise estimation module is used for determining a preset number of variances from the variance values according to the sequence from small to large of the variance values, and determining the weighted average value of the preset number of variances as the estimated noise variance of the image to be processed.
8. The apparatus according to claim 7, wherein the smoothing module is specifically configured to perform smoothing on the image to be processed by using the difference operator according to the following formula:
n ( x , y ) = Σ j = - l l Σ i = - l l d L ( x + i , y + j ) v ( x + i , y + j ) ;
wherein, (x, y) represents a pixel point in the image to be processed; v (x, y) represents the gray value of the pixel point (x, y) before the smoothing processing is performed on the pixel point (x, y); n (x, y) represents the gray value of the pixel point (x, y) after the smoothing processing is performed on the pixel point (x, y); the value of l is the number of rows or columns of the Laplace operator template minus 1 divided by 2.
9. A video image denoising apparatus, comprising:
a threshold determination module, configured to determine, by the image noise estimation apparatus according to any one of claims 7 to 8, a noise variance of a current frame image of the video to be processed, and determine a motion detection threshold according to the noise variance;
the second blocking module is used for respectively carrying out non-overlapping blocking on the current frame image and the previous frame image based on a second preset size to obtain each image block of the current frame image and each image block of the previous frame image;
the similarity determining module is used for determining the similarity between each image block of the current frame image and the image block corresponding to the position of the image block in the previous frame image; and are
The motion detection module is used for comparing the similarity with the motion detection threshold and determining the image block to be a motion image block or a static image block according to the comparison result;
and the denoising processing module is used for respectively denoising the moving image block and the static image block in the current frame image.
10. The apparatus of claim 9, further comprising:
the gradient value determining module is used for determining the gradient value of each pixel point in the current frame image;
the pixel point determining module is used for determining pixel points with gradient values larger than a preset gradient threshold value from all pixel points in the moving image block after the denoising processing module denoises the moving image block and the static image block in the current frame image respectively;
and the edge processing module is used for further determining the sum of the gray value obtained after the pixel is subjected to denoising processing and the gray value before the pixel is subjected to denoising processing according to a preset proportion as the gray value after the pixel is subjected to denoising processing aiming at each pixel in the determined pixels.
11. The apparatus according to claim 9, wherein the similarity determining module is specifically configured to determine, for each image block of the current frame image and the previous frame image, an image area formed by s × s image blocks centered on the image block in a frame in which the image block is located;
for each image block of the current frame, determining the similarity D (i, j, k) between the image block and the image block corresponding to the position of the image block in the previous frame by adopting the following formula:
D ( i , j , k ) = Σ m = i - h i + h Σ n = j - h j + h [ S B ( m , n , k ) - S B ( m , n , k - 1 ) ] 2
wherein, (i, j) represents the center coordinates of the image block; h has a value of (s-1)/2; k and k-1 represent frame numbers; (m, n) denotes coordinates of image blocks other than the image block with respect to the image block in an image area centered on the image block; sB(m, n, k) represents the sum of gray values of all pixel points in the image block with the center coordinate (m, n) in the kth frame; sBAnd (m, n, k-1) represents the sum of gray values of all pixel points in the image block with the central coordinate of (m, n) in the (k-1) th frame.
12. The apparatus according to claim 9, wherein the motion detection module is specifically configured to determine that the image block of the current frame is a motion image block relative to an image block corresponding to the image block position in a previous frame when the similarity is greater than the motion detection threshold; and when the similarity is not greater than the motion detection threshold, determining that the image block of the current frame is a static image block relative to the image block corresponding to the position of the image block in the previous frame.
CN201310603614.4A 2013-11-25 2013-11-25 Image noise estimation method, video image denoising method and device Active CN104680483B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310603614.4A CN104680483B (en) 2013-11-25 2013-11-25 Image noise estimation method, video image denoising method and device

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310603614.4A CN104680483B (en) 2013-11-25 2013-11-25 Image noise estimation method, video image denoising method and device

Publications (2)

Publication Number Publication Date
CN104680483A CN104680483A (en) 2015-06-03
CN104680483B true CN104680483B (en) 2016-03-02

Family

ID=53315482

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310603614.4A Active CN104680483B (en) 2013-11-25 2013-11-25 Image noise estimation method, video image denoising method and device

Country Status (1)

Country Link
CN (1) CN104680483B (en)

Families Citing this family (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105894457A (en) * 2015-12-04 2016-08-24 乐视致新电子科技(天津)有限公司 Image noise removing method and device
CN106127729A (en) * 2016-06-08 2016-11-16 浙江传媒学院 A kind of picture noise level estimation method based on gradient
CN110443794B (en) 2019-05-29 2020-12-18 腾讯科技(深圳)有限公司 Pathological image-based image state determination method, device and system
CN110363726A (en) * 2019-07-24 2019-10-22 陈泰杉 A kind of vision noise reduction process method based on regular graph Laplace transform
CN112311962B (en) 2019-07-29 2023-11-24 深圳市中兴微电子技术有限公司 Video denoising method and device and computer readable storage medium
CN113129221B (en) * 2019-12-31 2023-08-18 杭州海康威视数字技术股份有限公司 Image processing method, device, equipment and storage medium
CN111556227B (en) * 2020-05-19 2022-04-15 广州市百果园信息技术有限公司 Video denoising method and device, mobile terminal and storage medium
CN111724421B (en) * 2020-06-29 2024-01-09 深圳市慧鲤科技有限公司 Image processing method and device, electronic equipment and storage medium
CN114331901A (en) * 2021-12-30 2022-04-12 北京超维景生物科技有限公司 Model training method and model training device
CN115330628B (en) * 2022-08-18 2023-09-12 盐城众拓视觉创意有限公司 Video frame-by-frame denoising method based on image processing
CN118037584B (en) * 2024-04-12 2024-06-14 禹创半导体(深圳)有限公司 Overdrive-based image boundary smoothing method, device and medium

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102722874A (en) * 2012-05-29 2012-10-10 飞依诺科技(苏州)有限公司 Ultrasonic image speckle noise suppression method
CN103095967A (en) * 2011-10-28 2013-05-08 浙江大华技术股份有限公司 Video noise quantization calculation method and video noise quantization calculation system

Family Cites Families (1)

* Cited by examiner, † Cited by third party
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

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103095967A (en) * 2011-10-28 2013-05-08 浙江大华技术股份有限公司 Video noise quantization calculation method and video noise quantization calculation system
CN102722874A (en) * 2012-05-29 2012-10-10 飞依诺科技(苏州)有限公司 Ultrasonic image speckle noise suppression method

Also Published As

Publication number Publication date
CN104680483A (en) 2015-06-03

Similar Documents

Publication Publication Date Title
CN104680483B (en) Image noise estimation method, video image denoising method and device
CN108230359B (en) Object detection method and apparatus, training method, electronic device, program, and medium
CN103955930B (en) Motion parameter estimation method based on gray integral projection cross-correlation function characteristics
CN109791695B (en) Determining a variance of an image block based on a motion vector of the block
CN104899830B (en) A kind of image super-resolution method
CN104599286B (en) A kind of characteristic tracking method and device based on light stream
US20150254814A1 (en) Globally dominant point spread function estimation
US11093778B2 (en) Method and system for selecting image region that facilitates blur kernel estimation
Zhang et al. Salt and pepper noise removal with image inpainting
CN106127729A (en) A kind of picture noise level estimation method based on gradient
CN103886553A (en) Method and system for non-local average value denoising of image
CN102298773A (en) Shape-adaptive non-local mean denoising method
KR20150032822A (en) Method and apparatus for filtering an image
CN110570379B (en) Non-local mean value CT image noise reduction method based on structure tensor
CN108960012B (en) Feature point detection method and device and electronic equipment
CN105590307A (en) Transparency-based matting method and apparatus
CN112308872A (en) Image edge detection method based on multi-scale Gabor first-order derivative
CN103150725B (en) Based on SUSAN edge detection method and the system of non-local mean
CN108921170B (en) An effective image noise detection and denoising method and system
CN106056577A (en) Hybrid cascaded SAR image change detection method based on MDS-SRM
CN108810319B (en) Image processing apparatus, image processing method, and program
CN105279742B (en) A kind of image de-noising method quickly based on piecemeal estimation of noise energy
US9679363B1 (en) System and method for reducing image noise
CN110738625B (en) Image resampling method, device, terminal and computer readable storage medium
CN106846366A (en) Use the TLD video frequency motion target trackings of GPU hardware

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