CN104680483B - Image noise estimation method, video image denoising method and device - Google Patents
Image noise estimation method, video image denoising method and device 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
Landscapes
- Picture Signal Circuits (AREA)
- Image Analysis (AREA)
Abstract
Description
技术领域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:
其中,(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:
其中,(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):
针对静止图像块中的每个像素点,假设该像素点进行去噪处理之前的灰度值为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):
其中,in,
在式(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:
其中,(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:
其中,(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)
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)
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)
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)
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/en active Active
Patent Citations (2)
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 |