CN106447616B - 一种实现小波去噪的方法和装置 - Google Patents
一种实现小波去噪的方法和装置 Download PDFInfo
- Publication number
- CN106447616B CN106447616B CN201510492608.5A CN201510492608A CN106447616B CN 106447616 B CN106447616 B CN 106447616B CN 201510492608 A CN201510492608 A CN 201510492608A CN 106447616 B CN106447616 B CN 106447616B
- Authority
- CN
- China
- Prior art keywords
- image
- channel
- images
- pixel point
- layer
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 24
- 238000005070 sampling Methods 0.000 claims abstract description 75
- 230000002194 synthesizing effect Effects 0.000 claims abstract description 12
- 238000001914 filtration Methods 0.000 claims description 32
- 230000009466 transformation Effects 0.000 claims description 27
- 239000000126 substance Substances 0.000 claims description 8
- 239000011159 matrix material Substances 0.000 claims description 7
- 230000000694 effects Effects 0.000 abstract description 9
- 230000009977 dual effect Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Images
Classifications
-
- G06T5/70—
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/10—Image enhancement or restoration by non-spatial domain filtering
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20048—Transform domain processing
- G06T2207/20064—Wavelet transform [DWT]
Abstract
本发明公开了一种实现小波去噪的方法和装置，包括：对图像进行降采样，分别对降采样后的各图像的各通道的图像进行小波去噪得到降采样后的各图像的各通道的去噪图像；将得到的各去噪图像合成完整的去噪图像。本发明提高了去噪效果。
Description
技术领域
本发明涉及图像处理技术，尤指一种实现小波去噪的方法和装置。
背景技术
现有的实现小波去噪的方法大致包括：
将图像进行对偶延拓；对对偶延拓后的图像进行小波变换得到第一层小波系数；根据第一层小波系数的高频部分计算图像的噪声方差，根据图像的噪声方差和图像的方差计算第一层图像小波系数；对第i层小波系数的低频部分进行小波变换得到第(i+1)层小波系数；根据图像的噪声方差和第i层小波系数的低频部分计算第(i+1)层图像小波系数；其中，i为大于或等于1的整数；对计算得到的各层图像小波系数进行小波逆变换得到去噪图像。
现有的实现去噪的方法中，由于小波去噪只适用于高斯分布的噪声，因此，在噪声完全符合高斯分布时，是根据第一层小波系数的高频部分计算图像的噪声方差，而实际上噪声在图像中不完全符合高斯分布，因此，去噪效果不佳。
发明内容
为了解决上述问题，本发明提出了一种实现小波去噪的方法和装置，能够提高去噪效果。
为了达到上述目的，本发明提出了一种实现小波去噪的方法，至少包括：
对图像进行降采样，分别对降采样后的各图像的各通道的图像进行小波去噪得到降采样后的各图像的各通道的去噪图像；
将得到的各去噪图像合成完整的去噪图像。
优选地，还包括：对所述完整的去噪图像进行灰度化处理。
优选地，所述对完整的去噪图像进行灰度化处理包括：
对所述完整的去噪图像的Y通道的图像进行低通滤波后进行高通滤波；
将高通滤波后的图像转换成二值图像；
判断出所述完整的去噪图像的某一个像素点为噪声造成的彩色，且所述二值图像中与所述完整的去噪图像的某一个像素点对应的像素点的像素值为1，减小所述完整的去噪图像的某一个像素点的U通道和V通道的像素值。
优选地，所述将高通滤波后的图像转换成二值图像包括：
判断出所述高通滤波后的图像的某一个像素点的像素值大于或等于第一预设值，将所述二值图像与所述高通滤波后的图像的某一个像素点对应的像素点的像素值设为1；
判断出所述高通滤波后的图像的某一个像素点的像素值小于所述第一预设值，将所述二值图像与所述高通滤波后的图像的某一个像素点对应的像素点的像素值设为0。
优选地，当所述高通滤波后的图像的所有像素点的像素值的平均值小于或等于60/255时，所述第一预设值为45/255；
当所述平均值大于60/255，且小于或等于70/255时，所述第一预设值为40/255；
当所述平均值大于70/255，且小于或等于80/255时，所述第一预设值为35/255；
当所述平均值大于80/255，且小于或等于90/255时，所述第一预设值为30/255；
当所述平均值大于90/255时，所述第一预设值为20/255。
优选地，所述判断出完整的去噪图像的某一个像素点为噪声造成的彩色包括：
判断出所述完整的去噪图像的某一个像素点满足abs(1.1398v)＜t且abs(0.3946u+0.5806v)＜t且abs(2.0321u)＜t；其中，v为所述完整的去噪图像的某一个像素点的V通道的像素值，u为所述完整的去噪图像的某一个像素点的U通道的像素值，t为第二预设值；abs()表示取绝对值。
优选地，所述减小完整的去噪图像的某一个像素点的U通道和V通道的像素值包括：
计算所述完整的去噪图像的某一个像素点的U通道的像素值和第二预设值之间的比值作为所述完整的去噪图像的某一个像素点的U通道的新的像素值，计算所述完整的去噪图像的某一个像素点的V通道的像素值和所述第二预设值之间的比值作为所述完整的去噪图像的某一个像素点的V通道的新的像素值。
优选地，当所述完整的去噪图像的某一个像素点的Y通道的像素值小于或等于30/255时，所述第二预设值为1.8；
当所述完整的去噪图像的某一个像素点的Y通道的像素值大于30/255且小于或等于60/255时，所述第二预设值为1.6；
当所述完整的去噪图像的某一个像素点的Y通道的像素值大于60/255且小于或等于90/255时，所述第二预设值为1.4；
当所述完整的去噪图像的某一个像素点的Y通道的像素值大于90/255且小于或等于255/255时，所述第二预设值为1.2。
优选地，所述分别对降采样后的各图像的各通道的图像进行小波去噪得到降采样后的各图像的各通道的去噪图像包括：
分别对降采样后的各图像的各通道的图像进行小波变换得到降采样后的各图像的各通道的图像的第一层小波系数；
根据降采样后的各图像的各通道的图像的第一层小波系数的高频部分计算降采样后的各图像的各通道的图像的噪声方差，根据降采样后的各图像的各通道的图像的噪声方差和降采样后的各图像的各通道的图像的方差计算降采样后的各图像的各通道的图像的第一层图像小波系数；
对降采样后的各图像的各通道的图像的第i层小波系数的低频部分进行小波变换得到降采样后的各图像的各通道的图像的第(i+1)层小波系数；其中，i为大于或等于1的整数；
根据降采样后的各图像的各通道的图像的第(i+1)层小波系数的高频部分计算降采样后的各图像的各通道的图像的第i层小波系数的低频部分的噪声方差，根据降采样后的各图像的各通道的图像的第i层小波系数的低频部分的噪声方差和降采样后的各图像的各通道的图像的第i层小波系数的低频部分的方差计算降采样后的各图像的各通道的图像的第(i+1)层图像小波系数；
对降采样后的各图像的各通道的图像的各层图像小波系数分别进行小波逆变换得到降采样后的各图像的各通道的去噪图像。
优选地，所述根据降采样后的各图像的各通道的图像的第一层小波系数/第(i+1)层小波系数的高频部分计算降采样后的各图像的各通道的图像/降采样后的各图像的各通道的图像的第i层小波系数的低频部分的噪声方差包括：
按照公式
计算所述降采样后的各图像的各通道的图像/降采样后的各图像的各通道的图像的第i层小波系数的低频部分的噪声方差；median()表示取中值；
其中，σni为所述降采样后的各图像的各通道的图像/降采样后的各图像的各通道的图像的第i层小波系数的低频部分的噪声方差，yni+1为所述降采样后的各图像的各通道的图像的第一层小波系数/第(i+1)层小波系数的高频部分的像素矩阵。
优选地，所述根据降采样后的各图像的各通道的图像/降采样后的各图像的各通道的图像的第i层小波系数的低频部分的噪声方差和降采样后的各图像的各通道的图像/降采样后的各图像的各通道的图像的第i层小波系数的低频部分的方差计算降采样后的各图像的各通道的图像的第一层图像小波系数/第(i+1)层图像小波系数包括：
按照公式
计算所述降采样后的各图像的各通道的图像/降采样后的各图像的各通道的图像的第i层小波系数的低频部分的方差，按照公式
计算所述降采样后的各图像的各通道的图像的第一层图像小波系数/第(i+1)层图像小波系数；
其中，按照公式
计算Tik；
其中，
为所述降采样后的各图像的各通道的图像/降采样后的各图像的各通道的图像的第i层小波系数的低频部分的第k个像素点的方差，N为邻域窗口的像素数，yijk为所述降采样后的各图像的各通道的图像/降采样后的各图像的各通道的图像的第i层小波系数的低频部分的第k个像素点的邻域窗口的第j个像素点的像素值；j为1到N的整数，k为1到所述降采样后的各图像的各通道的图像/降采样后的各图像的各通道的图像的第i层小波系数的低频部分的像素数的整数，wi+1k为所述降采样后的各图像的各通道的图像的第一层图像小波系数/第(i+1)层图像小波系数的第k个像素点的像素值，yi+1k1为所述降采样后的各图像的各通道的图像的第一层图像小波系数/第(i+1)层图像小波系数的第k个像素点的像素值的实部，yi+1k2为所述降采样后的各图像的各通道的图像的第一层图像小波系数/第(i+1)层图像小波系数的第k个像素点的像素值的虚部。
本发明还提出了一种实现小波去噪的装置，至少包括：
降采样模块，用于对图像进行降采样；
去噪模块，用于分别对降采样后的各图像的各通道的图像进行小波去噪得到降采样后的各图像的各通道的去噪图像。
优选地，还包括：
灰度化处理模块，用于对所述完整的去噪图像进行灰度化处理。
优选地，所述灰度化处理模块具体用于：
对所述完整的去噪图像的Y通道的图像进行低通滤波后进行高通滤波；将高通滤波后的图像转换成二值图像；判断出所述完整的去噪图像的某一个像素点为噪声造成的彩色，且所述二值图像中与所述完整的去噪图像的某一个像素点对应的像素点的像素值为1，减小所述完整的去噪图像的某一个像素点的U通道和V通道的像素值。
优选地，所述灰度化处理模块具体用于：
对所述完整的去噪图像的Y通道的图像进行低通滤波后进行高通滤波；判断出所述高通滤波后的图像的某一个像素点的像素值大于或等于第一预设值，将所述二值图像与所述高通滤波后的图像的某一个像素点对应的像素点的像素值设为1；判断出所述高通滤波后的图像的某一个像素点的像素值小于所述第一预设值，将所述二值图像与所述高通滤波后的图像的某一个像素点对应的像素点的像素值设为0；判断出所述完整的去噪图像的某一个像素点为噪声造成的彩色，且所述二值图像中与所述完整的去噪图像的某一个像素点对应的像素点的像素值为1，减小所述完整的去噪图像的某一个像素点的U通道和V通道的像素值。
优选地，所述灰度化处理模块具体用于：
对所述完整的去噪图像的Y通道的图像进行低通滤波后进行高通滤波；将高通滤波后的图像转换成二值图像；判断出所述完整的去噪图像的某一个像素点满足abs(1.1398v)＜t且abs(0.3946u+0.5806v)＜t且abs(2.0321u)＜t；其中，v为所述完整的去噪图像的某一个像素点的V通道的像素值；abs()表示取绝对值，u为所述完整的去噪图像的某一个像素点的U通道的像素值，t为第二预设值，且所述二值图像中与所述完整的去噪图像的某一个像素点对应的像素点的像素值为1，减小所述完整的去噪图像的某一个像素点的U通道和V通道的像素值。
优选地，所述灰度化处理模块具体用于：
对所述完整的去噪图像的Y通道的图像进行低通滤波后进行高通滤波；将高通滤波后的图像转换成二值图像；判断出所述完整的去噪图像的某一个像素点为噪声造成的彩色，且所述二值图像中与所述完整的去噪图像的某一个像素点对应的像素点的像素值为1，计算所述完整的去噪图像的某一个像素点的U通道的像素值和第二预设值之间的比值作为所述完整的去噪图像的某一个像素点的U通道的新的像素值，计算所述完整的去噪图像的某一个像素点的V通道的像素值和所述第二预设值之间的比值作为所述完整的去噪图像的某一个像素点的V通道的新的像素值。
优选地，所述去噪模块具体用于：
分别对降采样后的各图像的各通道的图像进行小波变换得到降采样后的各图像的各通道的图像的第一层小波系数；根据降采样后的各图像的各通道的图像的第一层小波系数的高频部分计算降采样后的各图像的各通道的图像的噪声方差，根据降采样后的各图像的各通道的图像的噪声方差和降采样后的各图像的各通道的图像的方差计算降采样后的各图像的各通道的图像的第一层图像小波系数；对降采样后的各图像的各通道的图像的第i层小波系数的低频部分进行小波变换得到降采样后的各图像的各通道的图像的第(i+1)层小波系数；其中，i为大于或等于1的整数；根据降采样后的各图像的各通道的图像的第(i+1)层小波系数的高频部分计算降采样后的各图像的各通道的图像的第i层小波系数的低频部分的噪声方差，根据降采样后的各图像的各通道的图像的第i层小波系数的低频部分的噪声方差和降采样后的各图像的各通道的图像的第i层小波系数的低频部分的方差计算降采样后的各图像的各通道的图像的第(i+1)层图像小波系数；对降采样后的各图像的各通道的图像的各层图像小波系数分别进行小波逆变换得到降采样后的各图像的各通道的去噪图像，将得到的各去噪图像合成完整的去噪图像。
优选地，所述根据降采样后的各图像的各通道的图像的第一层小波系数/第(i+1)层小波系数的高频部分计算降采样后的各图像的各通道的图像/降采样后的各图像的各通道的图像的第i层小波系数的低频部分的噪声方差为：
按照公式
计算所述降采样后的各图像的各通道的图像/降采样后的各图像的各通道的图像的第i层小波系数的低频部分的噪声方差；median()表示取中值；
其中，σni为所述降采样后的各图像的各通道的图像/降采样后的各图像的各通道的图像的第i层小波系数的低频部分的噪声方差，yni+1为所述降采样后的各图像的各通道的图像的第一层小波系数/第(i+1)层小波系数的高频部分的像素矩阵。
优选地，所述根据降采样后的各图像的各通道的图像/降采样后的各图像的各通道的图像的第i层小波系数的低频部分的噪声方差和降采样后的各图像的各通道的图像/降采样后的各图像的各通道的图像的第i层小波系数的低频部分的方差计算降采样后的各图像的各通道的图像的第一层图像小波系数/第(i+1)层图像小波系数为：
按照公式
计算所述降采样后的各图像的各通道的图像/降采样后的各图像的各通道的图像的第i层小波系数的低频部分的方差，按照公式
计算所述降采样后的各图像的各通道的图像的第一层图像小波系数/第(i+1)层图像小波系数；
其中，按照公式
计算Tik；
其中，
为所述降采样后的各图像的各通道的图像/降采样后的各图像的各通道的图像的第i层小波系数的低频部分的第k个像素点的方差，N为邻域窗口的像素数，yijk为所述降采样后的各图像的各通道的图像/降采样后的各图像的各通道的图像的第i层小波系数的低频部分的第k个像素点的邻域窗口的第j个像素点的像素值；j为1到N的整数，k为1到所述降采样后的各图像的各通道的图像/降采样后的各图像的各通道的图像的第i层小波系数的低频部分的像素数的整数，wi+1k为所述降采样后的各图像的各通道的图像的第一层图像小波系数/第(i+1)层图像小波系数的第k个像素点的像素值，yi+1k1为所述降采样后的各图像的各通道的图像的第一层图像小波系数/第(i+1)层图像小波系数的第k个像素点的像素值的实部，yi+1k2为所述降采样后的各图像的各通道的图像的第一层图像小波系数/第(i+1)层图像小波系数的第k个像素点的像素值的虚部。
与现有技术相比，本发明的技术方案包括：对图像进行降采样，分别对降采样后的各图像的各通道的图像进行小波去噪得到降采样后的各图像的各通道的去噪图像；将得到的各去噪图像合成完整的去噪图像。通过本发明的方案，将图像进行降采样后再进行小波去噪，由于降采样后的图像的噪声更符合高斯分布，提高了去噪效果。
进一步地，根据降采样后的各图像的各通道的图像的第(i+1)层小波系数的高频部分计算降采样后的各图像的各通道的图像的第i层小波系数的低频部分的噪声方差，进一步提高了去噪效果。
进一步地，对完整的去噪图像进行灰度化处理，去除了平坦区域中的噪声点，使得去噪效果更加优化。
附图说明
下面对本发明实施例中的附图进行说明，实施例中的附图是用于对本发明的进一步理解，与说明书一起用于解释本发明，并不构成对本发明保护范围的限制。
图1为本发明实现小波去噪的方法的流程图；
图2为本发明实现小波去噪的装置的结构组成示意图。
具体实施方式
为了便于本领域技术人员的理解，下面结合附图对本发明作进一步的描述，并不能用来限制本发明的保护范围。需要说明的是，在不冲突的情况下，本申请中的实施例及实施例中的各种方式可以相互组合。
参见图1，本发明提出了一种实现小波去噪的方法，包括：
步骤100、对图像进行降采样，分别对降采样后的各图像的各通道的图像进行小波去噪得到降采样后的各图像的各通道的去噪图像。
本步骤中，分别对降采样后的各图像的各通道的图像进行小波去噪得到降采样后的各图像的各通道的去噪图像包括：
分别对降采样后的各图像的各通道的图像进行小波变换得到降采样后的各图像的各通道的图像的第一层小波系数；根据降采样后的各图像的各通道的图像的第一层小波系数的高频部分计算降采样后的各图像的各通道的图像的噪声方差，根据降采样后的各图像的各通道的图像的噪声方差和降采样后的各图像的各通道的图像的方差计算降采样后的各图像的各通道的图像的第一层图像小波系数；对降采样后的各图像的各通道的图像的第i层小波系数的低频部分进行小波变换得到降采样后的各图像的各通道的图像的第(i+1)层小波系数；根据降采样后的各图像的各通道的图像的第(i+1)层小波系数的高频部分计算降采样后的各图像的各通道的图像的第i层小波系数的低频部分的噪声方差，根据降采样后的各图像的各通道的图像的第i层小波系数的低频部分的噪声方差和降采样后的各图像的各通道的图像的第i层小波系数的低频部分的方差计算降采样后的各图像的各通道的图像的第(i+1)层图像小波系数；对降采样后的各图像的各通道的图像的各层图像小波系数分别进行小波逆变换得到降采样后的各图像的各通道的去噪图像。
其中，i为大于或等于1的整数。i的最大取值可以根据实际需要设定。
本步骤中，图像可以是YUV格式的图像，如果图像是RGB格式的图像，则可以采用公
式
将其转换成YUV格式的图像。那么，图像具
有三个通道。
本步骤中，降采样后的各图像应为正方形图像，且边长为2n，n为大于或等于1的整数。n具体的取值大小根据实际需要来确定，n越小，精度越高，然而计算量越大。
另外，在进行降采样前，如果图像的边长不为2n，则需要预先将图像填充使得图像的边长为2n。
例如，图像大小为2304×4096，进行填充后变成2560×4096，这样，就能进行降采样得到512×512的图像，将512×512的图像的3各通道分别进行小波去噪。
本步骤中，如何分别对降采样后的各图像的各通道的图像进行小波变换得到第一层小波系数属于本领域技术人员的公知技术，并不用于限定本发明的保护范围，这里不再赘述。
本步骤中，在进行小波变换前，还可以对降采样后的各图像的各通道的图像进行对偶延拓，具体实现属于本领域技术人员的公知技术，并不用于限定本发明的保护范围，这里不再赘述。
其中，根据降采样后的各图像的各通道的图像的第一层小波系数的高频部分计算降采样后的各图像的各通道的图像的噪声方差包括：
按照公式
计算所述降采样后的各图像的各通道的图像的噪声方差；median()表示取中值。
其中，Vn0为降采样后的各图像的各通道的图像的噪声方差，yn1为降采样后的各图像的各通道的图像的第一层小波系统的高频部分的像素矩阵。
其中，根据降采样后的各图像的各通道的图像的噪声方差和降采样后的各图像的各通道的图像的方差计算降采样后的各图像的各通道的图像的第一层图像小波系数包括：
按照公式
计算降采样后的各图像的各通道的图像的方差，按照公式
计算降采样后的各图像的各通道的图像的第一层图像小波系数/第(i+1)层图像小波系数；
其中，按照公式
计算T0k；
其中，
为降采样后的各图像的各通道的图像的第k个像素点的方差，N为邻域窗口的像素数，y0jk为降采样后的各图像的各通道的图像的第k个像素点的邻域窗口的第j个像素点的像素值；j为1到N的整数，k为1到降采样后的各图像的各通道的图像的像素数的整数，w1k为降采样后的各图像的各通道的图像的第一层图像小波系数的第k个像素点的像素值，y1k1为降采样后的各图像的各通道的图像的第一层图像小波系数的第k个像素点的像素值的实部，y1k2为降采样后的各图像的各通道的图像的第一层图像小波系数的第k个像素点的像素值的虚部。
其中，如何对降采样后的各图像的各通道的图像的第i层小波系数的低频部分进行小波变换得到降采样后的各图像的各通道的图像的第(i+1)层小波系数属于本领域技术人员的公知技术，并不用于限定本发明的保护范围，这里不再赘述。
其中，根据降采样后的各图像的各通道的图像的第(i+1)层小波系数的高频部分计算降采样后的各图像的各通道的图像的第i层小波系数的低频部分的噪声方差包括：
按照公式
计算降采样后的各图像的各通道的图像的第i层小波系数的低频部分的噪声方差；
其中，σni为降采样后的各图像的各通道的图像的第i层小波系数的低频部分的噪声方差，yni+1为降采样后的各图像的各通道的图像的第(i+1)层小波系数的高频部分的像素矩阵。
其中，根据降采样后的各图像的各通道的图像的第i层小波系数的低频部分的噪声方差和降采样后的各图像的各通道的图像的第i层小波系数的低频部分的方差计算降采样后的各图像的各通道的图像的第(i+1)层图像小波系数包括：
按照公式
计算降采样后的各图像的各通道的图像的第i层小波系数的低频部分的方差，按照公式
计算降采样后的各图像的各通道的图像的第(i+1)层图像小波系数；
其中，按照公式
计算Tik；
其中，
为降采样后的各图像的各通道的图像的第i层小波系数的低频部分的第k个像素点的方差，N为邻域窗口的像素数，yijk为降采样后的各图像的各通道的图像的第i层小波系数的低频部分的第k个像素点的邻域窗口的第j个像素点的像素值；j为1到N的整数，k为1到降采样后的各图像的各通道的图像的第i层小波系数的低频部分的像素数的整数，wi+1k为降采样后的各图像的各通道的图像的第(i+1)层图像小波系数的第k个像素点的像素值，yi+1k1为降采样后的各图像的各通道的图像的第(i+1)层图像小波系数的第k个像素点的像素值的实部，yi+1k2为降采样后的各图像的各通道的图像的第(i+1)层图像小波系数的第k个像素点的像素值的虚部。
其中，如何对降采样后的各图像的各通道的图像的各层图像小波系数分别进行小波逆变换得到降采样后的各图像的各通道的去噪图像，将得到的各去噪图像合成完整的去噪图像属于本领域技术人员的公知技术，并不用于限定本发明的保护范围，这里不再赘述。
步骤101、将得到的各去噪图像合成完整的去噪图像。
通过本发明的方案，将图像进行降采样后再进行小波去噪，由于降采样后的图像的噪声更符合高斯分布，提高了去噪效果。
进一步地，根据降采样后的各图像的各通道的图像的第(i+1)层小波系数的高频部分计算降采样后的各图像的各通道的图像的第i层小波系数的低频部分的噪声方差，进一步提高了去噪效果。
进一步地，该方法还包括：
步骤102、对完整的去噪图像进行灰度化处理。具体包括：
对完整的去噪图像的Y通道的图像进行低通滤波后进行高通滤波；将高通滤波后的图像转换成二值图像；判断出完整的去噪图像的某一个像素点为噪声造成的彩色，且二值图像中与完整的去噪图像的某一个像素点对应的像素点的像素值为1，减小完整的去噪图像的某一个像素点的U通道和V通道的像素值。
其中，对完整的去噪图像的Y通道的图像进行低通滤波后进行高通滤波就是指先将完整的去噪图像的Y通道的图像与低通滤波器进行卷积，再与高通滤波器进行卷积。
其中，低通滤波器可以是[0.125，0.375，0.125]T。
高通滤波器可以是[0.125，0，0.375，0，0.375，0，0.125]T。
其中，将高通滤波后的图像转换成二值图像包括：
判断出高通滤波后的图像的某一个像素点的像素值大于或等于第一预设值，将二值图像与高通滤波后的图像的某一个像素点对应的像素点的像素值设为1；
判断出高通滤波后的图像的某一个像素点的像素值小于第一预设值，将二值图像与高通滤波后的图像的某一个像素点对应的像素点的像素值设为0。
其中，当高通滤波后的图像的所有像素点的像素值的平均值小于或等于60/255时，第一预设值为45/255；当平均值大于60/255，且小于或等于70/255时，第一预设值为40/255；当平均值大于70/255，且小于或等于80/255时，第一预设值为35/255；当平均值大于80/255，且小于或等于90/255时，第一预设值为30/255；当平均值大于90/255时，第一预设值为20/255。
其中，判断出完整的去噪图像的某一个像素点为噪声造成的彩色包括：
判断出完整的去噪图像的某一个像素点满足abs(1.1398v)＜t且abs(0.3946u+0.5806v)＜t且abs(2.0321u)＜t；其中，v为完整的去噪图像的某一个像素点的V通道的像素值，u为完整的去噪图像的某一个像素点的U通道的像素值，t为第二预设值；abs()表示取绝对值。
其中，t可以取值为25/255。
其中，减小完整的去噪图像的某一个像素点的U通道和V通道的像素值包括：
计算完整的去噪图像的某一个像素点的U通道的像素值和第二预设值之间的比值作为完整的去噪图像的某一个像素点的U通道的新的像素值，计算完整的去噪图像的某一个像素点的V通道的像素值和第二预设值之间的比值作为完整的去噪图像的某一个像素点的V通道的新的像素值。
其中，当完整的去噪图像的某一个像素点的Y通道的像素值小于或等于30/255时，第二预设值为1.8；当完整的去噪图像的某一个像素点的Y通道的像素值大于30/255且小于或等于60/255时，第二预设值为1.6；当完整的去噪图像的某一个像素点的Y通道的像素值大于60/255且小于或等于90/255时，第二预设值为1.4；当完整的去噪图像的某一个像素点的Y通道的像素值大于90/255且小于或等于255/255时，第二预设值为1.2。
本步骤中，对完整的去噪图像进行灰度化处理，去除了平坦区域中的噪声点，使得去噪效果更加优化。
参见图2，本发明还提出了一种实现小波去噪的装置，至少包括：
降采样模块，用于对图像进行降采样；
去噪模块，用于分别对降采样后的各图像的各通道的图像进行小波去噪得到降采样后的各图像的各通道的去噪图像。
本发明的装置中，还包括：
灰度化处理模块，用于对所述完整的去噪图像进行灰度化处理。
本发明的装置中，所述灰度化处理模块具体用于：
对所述完整的去噪图像的Y通道的图像进行低通滤波后进行高通滤波；将高通滤波后的图像转换成二值图像；判断出所述完整的去噪图像的某一个像素点为噪声造成的彩色，且所述二值图像中与所述完整的去噪图像的某一个像素点对应的像素点的像素值为1，减小所述完整的去噪图像的某一个像素点的U通道和V通道的像素值。
本发明的装置中，所述灰度化处理模块具体用于：
对所述完整的去噪图像的Y通道的图像进行低通滤波后进行高通滤波；判断出所述高通滤波后的图像的某一个像素点的像素值大于或等于第一预设值，将所述二值图像与所述高通滤波后的图像的某一个像素点对应的像素点的像素值设为1；判断出所述高通滤波后的图像的某一个像素点的像素值小于所述第一预设值，将所述二值图像与所述高通滤波后的图像的某一个像素点对应的像素点的像素值设为0；判断出所述完整的去噪图像的某一个像素点为噪声造成的彩色，且所述二值图像中与所述完整的去噪图像的某一个像素点对应的像素点的像素值为1，减小所述完整的去噪图像的某一个像素点的U通道和V通道的像素值。
本发明的装置中，所述灰度化处理模块具体用于：
对所述完整的去噪图像的Y通道的图像进行低通滤波后进行高通滤波；将高通滤波后的图像转换成二值图像；判断出所述完整的去噪图像的某一个像素点满足abs(1.1398v)＜t且abs(0.3946u+0.5806v)＜t且abs(2.0321u)＜t；其中，v为所述完整的去噪图像的某一个像素点的V通道的像素值；abs()表示取绝对值，u为所述完整的去噪图像的某一个像素点的U通道的像素值，t为第二预设值，且所述二值图像中与所述完整的去噪图像的某一个像素点对应的像素点的像素值为1，减小所述完整的去噪图像的某一个像素点的U通道和V通道的像素值。
本发明的装置中，所述灰度化处理模块具体用于：
对所述完整的去噪图像的Y通道的图像进行低通滤波后进行高通滤波；将高通滤波后的图像转换成二值图像；判断出所述完整的去噪图像的某一个像素点为噪声造成的彩色，且所述二值图像中与所述完整的去噪图像的某一个像素点对应的像素点的像素值为1，计算所述完整的去噪图像的某一个像素点的U通道的像素值和第二预设值之间的比值作为所述完整的去噪图像的某一个像素点的U通道的新的像素值，计算所述完整的去噪图像的某一个像素点的V通道的像素值和所述第二预设值之间的比值作为所述完整的去噪图像的某一个像素点的V通道的新的像素值。
本发明的装置中，去噪模块具体用于：
分别对降采样后的各图像的各通道的图像进行小波变换得到降采样后的各图像的各通道的图像的第一层小波系数；根据降采样后的各图像的各通道的图像的第一层小波系数的高频部分计算降采样后的各图像的各通道的图像的噪声方差，根据降采样后的各图像的各通道的图像的噪声方差和降采样后的各图像的各通道的图像的方差计算降采样后的各图像的各通道的图像的第一层图像小波系数；对降采样后的各图像的各通道的图像的第i层小波系数的低频部分进行小波变换得到降采样后的各图像的各通道的图像的第(i+1)层小波系数；其中，i为大于或等于1的整数；根据降采样后的各图像的各通道的图像的第(i+1)层小波系数的高频部分计算降采样后的各图像的各通道的图像的第i层小波系数的低频部分的噪声方差，根据降采样后的各图像的各通道的图像的第i层小波系数的低频部分的噪声方差和降采样后的各图像的各通道的图像的第i层小波系数的低频部分的方差计算降采样后的各图像的各通道的图像的第(i+1)层图像小波系数；对降采样后的各图像的各通道的图像的各层图像小波系数分别进行小波逆变换得到降采样后的各图像的各通道的去噪图像，将得到的各去噪图像合成完整的去噪图像。
本发明的装置中，去噪模块具体用于：
分别对降采样后的各图像的各通道的图像进行小波变换得到降采样后的各图像的各通道的图像的第一层小波系数；按照公式
计算所述降采样后的各图像的各通道的图像的噪声方差，根据降采样后的各图像的各通道的图像的噪声方差和降采样后的各图像的各通道的图像的方差计算降采样后的各图像的各通道的图像的第一层图像小波系数；对降采样后的各图像的各通道的图像的第i层小波系数的低频部分进行小波变换得到降采样后的各图像的各通道的图像的第(i+1)层小波系数；其中，i为大于或等于1的整数；按照公式
计算降采样后的各图像的各通道的图像的第i层小波系数的低频部分的噪声方差，根据降采样后的各图像的各通道的图像的第i层小波系数的低频部分的噪声方差和降采样后的各图像的各通道的图像的第i层小波系数的低频部分的方差计算降采样后的各图像的各通道的图像的第(i+1)层图像小波系数；对降采样后的各图像的各通道的图像的、的各层图像小波系数分别进行小波逆变换得到降采样后的各图像的各通道的去噪图像，将得到的各去噪图像合成完整的去噪图像。
本发明的装置中，去噪模块具体用于：
分别对降采样后的各图像的各通道的图像进行小波变换得到降采样后的各图像的各通道的图像的第一层小波系数；根据降采样后的各图像的各通道的图像的第一层小波系数的高频部分计算降采样后的各图像的各通道的图像的噪声方差，按照公式
计算降采样后的各图像的各通道的图像的方差，按照公式
计算所述降采样后的各图像的各通道的图像的第一层图像小波系数；其中，按照公式
计算T0k；；对降采样后的各图像的各通道的图像的第i层小波系数的低频部分进行小波变换得到降采样后的各图像的各通道的图像的第(i+1)层小波系数；其中，i为大于或等于1的整数；根据降采样后的各图像的各通道的图像的第(i+1)层小波系数的高频部分计算降采样后的各图像的各通道的图像的第i层小波系数的低频部分的噪声方差，按照公式
计算降采样后的各图像的各通道的图像的第i层小波系数的低频部分的方差，按照公式
计算降采样后的各图像的各通道的图像的第(i+1)层图像小波系数；其中，按照公式
计算Tik；对降采样后的各图像的各通道的图像的各层图像小波系数进行小波逆变换得到降采样后的各图像的各通道的去噪图像，将得到的各去噪图像合成完整的去噪图像。
本发明的所有模块的功能均可以是通过处理器执行存储在存储器中的程序/指令实现，还可以通过固件/逻辑电路/集成电路实现。
需要说明的是，以上所述的实施例仅是为了便于本领域的技术人员理解而已，并不用于限制本发明的保护范围，在不脱离本发明的发明构思的前提下，本领域技术人员对本发明所做出的任何显而易见的替换和改进等均在本发明的保护范围之内。
Claims (16)
1.一种实现小波去噪的方法，其特征在于，至少包括：
对图像进行降采样，分别对降采样后的各图像的各通道的图像进行小波去噪得到降采样后的各图像的各通道的去噪图像；
将得到的各去噪图像合成完整的去噪图像；
对所述完整的去噪图像的Y通道的图像进行低通滤波后进行高通滤波；
将高通滤波后的图像转换成二值图像；
判断出所述完整的去噪图像的某一个像素点为噪声造成的彩色，且所述二值图像中与所述完整的去噪图像的某一个像素点对应的像素点的像素值为1，减小所述完整的去噪图像的某一个像素点的U通道和V通道的像素值。
2.根据权利要求1所述的方法，其特征在于，所述将高通滤波后的图像转换成二值图像包括：
判断出所述高通滤波后的图像的某一个像素点的像素值大于或等于第一预设值，将所述二值图像与所述高通滤波后的图像的某一个像素点对应的像素点的像素值设为1；
判断出所述高通滤波后的图像的某一个像素点的像素值小于所述第一预设值，将所述二值图像与所述高通滤波后的图像的某一个像素点对应的像素点的像素值设为0。
3.根据权利要求2所述的方法，其特征在于，
当所述高通滤波后的图像的所有像素点的像素值的平均值小于或等于60/255时，所述第一预设值为45/255；
当所述平均值大于60/255，且小于或等于70/255时，所述第一预设值为40/255；
当所述平均值大于70/255，且小于或等于80/255时，所述第一预设值为35/255；
当所述平均值大于80/255，且小于或等于90/255时，所述第一预设值为30/255；
当所述平均值大于90/255时，所述第一预设值为20/255。
4.根据权利要求1所述的方法，其特征在于，所述判断出完整的去噪图像的某一个像素点为噪声造成的彩色包括：
判断出所述完整的去噪图像的某一个像素点满足abs(1.1398v)＜t且abs(0.3946u+0.5806v)＜t且abs(2.0321u)＜t；其中，v为所述完整的去噪图像的某一个像素点的V通道的像素值，u为所述完整的去噪图像的某一个像素点的U通道的像素值，t为第二预设值；abs()表示取绝对值。
5.根据权利要求1所述的方法，其特征在于，所述减小完整的去噪图像的某一个像素点的U通道和V通道的像素值包括：
计算所述完整的去噪图像的某一个像素点的U通道的像素值和第二预设值之间的比值作为所述完整的去噪图像的某一个像素点的U通道的新的像素值，计算所述完整的去噪图像的某一个像素点的V通道的像素值和所述第二预设值之间的比值作为所述完整的去噪图像的某一个像素点的V通道的新的像素值。
6.根据权利要求5所述的方法，其特征在于，当所述完整的去噪图像的某一个像素点的Y通道的像素值小于或等于30/255时，所述第二预设值为1.8；
当所述完整的去噪图像的某一个像素点的Y通道的像素值大于30/255且小于或等于60/255时，所述第二预设值为1.6；
当所述完整的去噪图像的某一个像素点的Y通道的像素值大于60/255且小于或等于90/255时，所述第二预设值为1.4；
当所述完整的去噪图像的某一个像素点的Y通道的像素值大于90/255且小于或等于255/255时，所述第二预设值为1.2。
7.根据权利要求1所述的方法，其特征在于，所述分别对降采样后的各图像的各通道的图像进行小波去噪得到降采样后的各图像的各通道的去噪图像包括：
分别对降采样后的各图像的各通道的图像进行小波变换得到降采样后的各图像的各通道的图像的第一层小波系数；
根据降采样后的各图像的各通道的图像的第一层小波系数的高频部分计算降采样后的各图像的各通道的图像的噪声方差，根据降采样后的各图像的各通道的图像的噪声方差和降采样后的各图像的各通道的图像的方差计算降采样后的各图像的各通道的图像的第一层图像小波系数；
对降采样后的各图像的各通道的图像的第i层小波系数的低频部分进行小波变换得到降采样后的各图像的各通道的图像的第(i+1)层小波系数；其中，i为大于或等于1的整数；
根据降采样后的各图像的各通道的图像的第(i+1)层小波系数的高频部分计算降采样后的各图像的各通道的图像的第i层小波系数的低频部分的噪声方差，根据降采样后的各图像的各通道的图像的第i层小波系数的低频部分的噪声方差和降采样后的各图像的各通道的图像的第i层小波系数的低频部分的方差计算降采样后的各图像的各通道的图像的第(i+1)层图像小波系数；
对降采样后的各图像的各通道的图像的各层图像小波系数分别进行小波逆变换得到降采样后的各图像的各通道的去噪图像。
8.根据权利要求7所述的方法，其特征在于，
所述根据降采样后的各图像的各通道的图像的第一层小波系数的高频部分计算降采样后的各图像的各通道的图像的噪声方差包括：
按照公式
计算降采样后的各图像的各通道的图像的噪声方差；median()表示取中值；
其中，σn0为降采样后的各图像的各通道的图像的噪声方差，yn1为降采样后的各图像的各通道的图像的第一层小波系统的高频部分的像素矩阵；
所述根据降采样后的各图像的各通道的图像的第(i+1)层小波系数的高频部分计算降采样后的各图像的各通道的图像的第i层小波系数的低频部分的噪声方差包括：
按照公式
计算降采样后的各图像的各通道的图像的第i层小波系数的低频部分的噪声方差；median()表示取中值；
其中，σni为降采样后的各图像的各通道的图像的第i层小波系数的低频部分的噪声方差，yni+1为所述降采样后的各图像的各通道的图像的第(i+1)层小波系数的高频部分的像素矩阵。
9.根据权利要求7所述的方法，其特征在于，
所述根据降采样后的各图像的各通道的图像的噪声方差和降采样后的各图像的各通道的图像的方差计算降采样后的各图像的各通道的图像的第一层图像小波系数包括：
按照公式
计算降采样后的各图像的各通道的图像的方差，按照公式
计算降采样后的各图像的各通道的图像的第一层图像小波系数；
其中，按照公式
计算T0k；
其中，
为降采样后的各图像的各通道的图像的第k个像素点的方差，N为邻域窗口的像素数，y0jk为降采样后的各图像的各通道的图像的第k个像素点的邻域窗口的第j个像素点的像素值；j为1到N的整数，k为1到降采样后的各图像的各通道的图像的像素数的整数，w1k为降采样后的各图像的各通道的图像的第一层图像小波系数的第k个像素点的像素值，y1k1为降采样后的各图像的各通道的图像的第一层图像小波系数的第k个像素点的像素值的实部，y1k2为降采样后的各图像的各通道的图像的第一层图像小波系数的第k个像素点的像素值的虚部；
所述根据降采样后的各图像的各通道的图像的第i层小波系数的低频部分的噪声方差和降采样后的各图像的各通道的图像的第i层小波系数的低频部分的方差计算降采样后的各图像的各通道的图像的第(i+1)层图像小波系数包括：
按照公式
计算所述降采样后的各图像的各通道的图像的第i层小波系数的低频部分的方差，按照公式
计算所述降采样后的各图像的各通道的图像的第(i+1)层图像小波系数；
其中，按照公式
计算Tik；
其中，
为所述降采样后的各图像的各通道的图像的第i层小波系数的低频部分的第k个像素点的方差，N为邻域窗口的像素数，yijk为所述降采样后的各图像的各通道的图像的第i层小波系数的低频部分的第k个像素点的邻域窗口的第j个像素点的像素值；j为1到N的整数，k为1到降采样后的各图像的各通道的图像的第i层小波系数的低频部分的像素数的整数，wi+1k为所述第(i+1)层图像小波系数的第k个像素点的像素值，yi+1k1为所述降采样后的各图像的各通道的图像的第(i+1)层图像小波系数的第k个像素点的像素值的实部，yi+1k2为所述降采样后的各图像的各通道的图像第(i+1)层图像小波系数的第k个像素点的像素值的虚部。
10.一种实现小波去噪的装置，其特征在于，至少包括：
降采样模块，用于对图像进行降采样；
去噪模块，用于分别对降采样后的各图像的各通道的图像进行小波去噪得到降采样后的各图像的各通道的去噪图像，将得到的各去噪图像合成完整的去噪图像；
灰度化处理模块，用于对所述完整的去噪图像的Y通道的图像进行低通滤波后进行高通滤波；将高通滤波后的图像转换成二值图像；判断出所述完整的去噪图像的某一个像素点为噪声造成的彩色，且所述二值图像中与所述完整的去噪图像的某一个像素点对应的像素点的像素值为1，减小所述完整的去噪图像的某一个像素点的U通道和V通道的像素值。
11.根据权利要求10所述的装置，其特征在于，所述灰度化处理模块具体用于：
对所述完整的去噪图像的Y通道的图像进行低通滤波后进行高通滤波；判断出所述高通滤波后的图像的某一个像素点的像素值大于或等于第一预设值，将所述二值图像与所述高通滤波后的图像的某一个像素点对应的像素点的像素值设为1；判断出所述高通滤波后的图像的某一个像素点的像素值小于所述第一预设值，将所述二值图像与所述高通滤波后的图像的某一个像素点对应的像素点的像素值设为0；判断出所述完整的去噪图像的某一个像素点为噪声造成的彩色，且所述二值图像中与所述完整的去噪图像的某一个像素点对应的像素点的像素值为1，减小所述完整的去噪图像的某一个像素点的U通道和V通道的像素值。
12.根据权利要求10所述的装置，其特征在于，所述灰度化处理模块具体用于：
对所述完整的去噪图像的Y通道的图像进行低通滤波后进行高通滤波；将高通滤波后的图像转换成二值图像；判断出所述完整的去噪图像的某一个像素点满足abs(1.1398v)＜t且abs(0.3946u+0.5806v)＜t且abs(2.0321u)＜t；其中，v为所述完整的去噪图像的某一个像素点的V通道的像素值；abs()表示取绝对值，u为所述完整的去噪图像的某一个像素点的U通道的像素值，t为第二预设值，且所述二值图像中与所述完整的去噪图像的某一个像素点对应的像素点的像素值为1，减小所述完整的去噪图像的某一个像素点的U通道和V通道的像素值。
13.根据权利要求10所述的装置，其特征在于，所述灰度化处理模块具体用于：
对所述完整的去噪图像的Y通道的图像进行低通滤波后进行高通滤波；将高通滤波后的图像转换成二值图像；判断出所述完整的去噪图像的某一个像素点为噪声造成的彩色，且所述二值图像中与所述完整的去噪图像的某一个像素点对应的像素点的像素值为1，计算所述完整的去噪图像的某一个像素点的U通道的像素值和第二预设值之间的比值作为所述完整的去噪图像的某一个像素点的U通道的新的像素值，计算所述完整的去噪图像的某一个像素点的V通道的像素值和所述第二预设值之间的比值作为所述完整的去噪图像的某一个像素点的V通道的新的像素值。
14.根据权利要求10所述的装置，其特征在于，所述去噪模块具体用于：
分别对降采样后的各图像的各通道的图像进行小波变换得到降采样后的各图像的各通道的图像的第一层小波系数；根据降采样后的各图像的各通道的图像的第一层小波系数的高频部分计算降采样后的各图像的各通道的图像的噪声方差，根据降采样后的各图像的各通道的图像的噪声方差和降采样后的各图像的各通道的图像的方差计算降采样后的各图像的各通道的图像的第一层图像小波系数；对降采样后的各图像的各通道的图像的第i层小波系数的低频部分进行小波变换得到降采样后的各图像的各通道的图像的第(i+1)层小波系数；其中，i为大于或等于1的整数；根据降采样后的各图像的各通道的图像的第(i+1)层小波系数的高频部分计算降采样后的各图像的各通道的图像的第i层小波系数的低频部分的噪声方差，根据降采样后的各图像的各通道的图像的第i层小波系数的低频部分的噪声方差和降采样后的各图像的各通道的图像的第i层小波系数的低频部分的方差计算降采样后的各图像的各通道的图像的第(i+1)层图像小波系数；对降采样后的各图像的各通道的图像的各层图像小波系数分别进行小波逆变换得到降采样后的各图像的各通道的去噪图像，将得到的各去噪图像合成完整的去噪图像。
15.根据权利要求14所述的装置，其特征在于，
所述根据降采样后的各图像的各通道的图像的第一层小波系数的高频部分计算降采样后的各图像的各通道的图像的噪声方差包括：
按照公式
计算降采样后的各图像的各通道的图像的噪声方差；median()表示取中值；
其中，σn0为降采样后的各图像的各通道的图像的噪声方差，yn1为降采样后的各图像的各通道的图像的第一层小波系统的高频部分的像素矩阵；
所述根据降采样后的各图像的各通道的图像的第(i+1)层小波系数的高频部分计算降采样后的各图像的各通道的图像的第i层小波系数的低频部分的噪声方差包括：
按照公式
计算降采样后的各图像的各通道的图像的第i层小波系数的低频部分的噪声方差；median()表示取中值；
其中，σni为降采样后的各图像的各通道的图像的第i层小波系数的低频部分的噪声方差，yni+1为所述降采样后的各图像的各通道的图像的第(i+1)层小波系数的高频部分的像素矩阵。
16.根据权利要求14所述的装置，其特征在于，
所述根据降采样后的各图像的各通道的图像的噪声方差和降采样后的各图像的各通道的图像的方差计算降采样后的各图像的各通道的图像的第一层图像小波系数包括：
按照公式
计算降采样后的各图像的各通道的图像的方差，按照公式
计算降采样后的各图像的各通道的图像的第一层图像小波系数；
其中，按照公式
计算T0k；
其中，
为降采样后的各图像的各通道的图像的第k个像素点的方差，N为邻域窗口的像素数，y0jk为降采样后的各图像的各通道的图像的第k个像素点的邻域窗口的第j个像素点的像素值；j为1到N的整数，k为1到降采样后的各图像的各通道的图像的像素数的整数，w1k为降采样后的各图像的各通道的图像的第一层图像小波系数的第k个像素点的像素值，y1k1为降采样后的各图像的各通道的图像的第一层图像小波系数的第k个像素点的像素值的实部，y1k2为降采样后的各图像的各通道的图像的第一层图像小波系数的第k个像素点的像素值的虚部；
所述根据降采样后的各图像的各通道的图像的第i层小波系数的低频部分的噪声方差和降采样后的各图像的各通道的图像的第i层小波系数的低频部分的方差计算降采样后的各图像的各通道的图像的第(i+1)层图像小波系数包括：
按照公式
计算所述降采样后的各图像的各通道的图像的第i层小波系数的低频部分的方差，按照公式
计算所述降采样后的各图像的各通道的图像的第(i+1)层图像小波系数；
其中，按照公式
计算Tik；
其中，
为所述降采样后的各图像的各通道的图像的第i层小波系数的低频部分的第k个像素点的方差，N为邻域窗口的像素数，yijk为所述降采样后的各图像的各通道的图像的第i层小波系数的低频部分的第k个像素点的邻域窗口的第j个像素点的像素值；j为1到N的整数，k为1到降采样后的各图像的各通道的图像的第i层小波系数的低频部分的像素数的整数，wi+1k为所述第(i+1)层图像小波系数的第k个像素点的像素值，yi+1k1为所述降采样后的各图像的各通道的图像的第(i+1)层图像小波系数的第k个像素点的像素值的实部，yi+1k2为所述降采样后的各图像的各通道的图像第(i+1)层图像小波系数的第k个像素点的像素值的虚部。
Priority Applications (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510492608.5A CN106447616B (zh) | 2015-08-12 | 2015-08-12 | 一种实现小波去噪的方法和装置 |
| PCT/CN2016/072312 WO2016131370A1 (zh) | 2015-08-12 | 2016-01-27 | 一种实现小波去噪的方法和装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510492608.5A CN106447616B (zh) | 2015-08-12 | 2015-08-12 | 一种实现小波去噪的方法和装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106447616A CN106447616A (zh) | 2017-02-22 |
| CN106447616B true CN106447616B (zh) | 2021-10-08 |
Family
ID=56692192
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510492608.5A Active CN106447616B (zh) | 2015-08-12 | 2015-08-12 | 一种实现小波去噪的方法和装置 |
Country Status (2)
| Country | Link |
|---|---|
| CN (1) | CN106447616B (zh) |
| WO (1) | WO2016131370A1 (zh) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110246087B (zh) * | 2018-03-07 | 2021-06-04 | 舜宇光学（浙江）研究院有限公司 | 参考多通道的多分辨率去除图像色度噪声的系统和方法 |
| CN109993707B (zh) * | 2019-03-01 | 2023-05-12 | 华为技术有限公司 | 图像去噪方法和装置 |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101201937A (zh) * | 2007-09-18 | 2008-06-18 | 上海医疗器械厂有限公司 | 基于小波重构与分解的数字图像增强方法及其装置 |
| CN102110289A (zh) * | 2011-03-29 | 2011-06-29 | 东南大学 | 基于变分框架的彩色图像对比度增强方法 |
| CN102496144A (zh) * | 2011-11-15 | 2012-06-13 | 哈尔滨工程大学 | 基于hsv色彩空间的nsct水声图像增强方法 |
| CN104715461A (zh) * | 2015-04-02 | 2015-06-17 | 哈尔滨理工大学 | 图像去噪方法 |
Family Cites Families (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2000307938A (ja) * | 1999-03-25 | 2000-11-02 | Texas Instr Inc <Ti> | カメラにおける画像を予見する方法及びディジタルカメラ |
| CN101442673B (zh) * | 2008-12-19 | 2010-08-11 | 太原理工大学 | 贝尔模板图像编码与解码方法 |
| CN102184523B (zh) * | 2011-04-13 | 2012-11-21 | 展讯通信（上海）有限公司 | 基于多分辨率和过完备变换的数字图像降噪方法及系统 |
| US9503054B2 (en) * | 2013-03-14 | 2016-11-22 | Avago Technolgies General IP (Singapore) Pte. Ltd. | Linear phase FIR biorthogonal wavelet filters with complementarity for image noise reduction |
| US9286662B2 (en) * | 2013-09-26 | 2016-03-15 | Siemens Aktiengesellschaft | Single image super resolution and denoising using multiple wavelet domain sparsity |
| CN104008538B (zh) * | 2014-05-12 | 2017-03-01 | 清华大学 | 基于单张图像超分辨率方法 |
-
2015
- 2015-08-12 CN CN201510492608.5A patent/CN106447616B/zh active Active
-
2016
- 2016-01-27 WO PCT/CN2016/072312 patent/WO2016131370A1/zh active Application Filing
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101201937A (zh) * | 2007-09-18 | 2008-06-18 | 上海医疗器械厂有限公司 | 基于小波重构与分解的数字图像增强方法及其装置 |
| CN102110289A (zh) * | 2011-03-29 | 2011-06-29 | 东南大学 | 基于变分框架的彩色图像对比度增强方法 |
| CN102496144A (zh) * | 2011-11-15 | 2012-06-13 | 哈尔滨工程大学 | 基于hsv色彩空间的nsct水声图像增强方法 |
| CN104715461A (zh) * | 2015-04-02 | 2015-06-17 | 哈尔滨理工大学 | 图像去噪方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106447616A (zh) | 2017-02-22 |
| WO2016131370A1 (zh) | 2016-08-25 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| EP3087730B1 (en) | Method for inverse tone mapping of an image | |
| CN108122213B (zh) | 一种基于YCrCb的低对比度图像增强方法 | |
| CN104240194B (zh) | 一种基于抛物线函数的低照度图像增强算法 | |
| US10672112B2 (en) | Method and system for real-time noise removal and image enhancement of high-dynamic range images | |
| WO2016206087A1 (zh) | 一种低照度图像处理方法和装置 | |
| Pang et al. | Improved single image dehazing using guided filter | |
| CN103186888B (zh) | 一种去除ct图像噪声的方法及装置 | |
| JP4862897B2 (ja) | 画像処理方法 | |
| CN103440623A (zh) | 基于成像模型提高雾天图像清晰度的方法 | |
| CN105225210A (zh) | 一种基于暗通道的自适应直方图增强去雾方法 | |
| CN106780417A (zh) | 一种光照不均图像的增强方法及系统 | |
| Manju et al. | Improved method for enhancing dark images based on CLAHE and morphological reconstruction | |
| CN101877123A (zh) | 一种图像增强方法和装置 | |
| CN110599418B (zh) | 一种变换域融合的全局色调映射方法 | |
| CN112508816B (zh) | 一种红外图像锐化方法、锐化处理系统及终端设备 | |
| CN106447616B (zh) | 一种实现小波去噪的方法和装置 | |
| CN107451986B (zh) | 一种基于融合技术的单幅红外图像增强方法 | |
| CN102903077A (zh) | 一种快速图像去模糊算法 | |
| CN109064413B (zh) | 图像对比度增强方法及采用其的图像采集医疗设备 | |
| US11915392B2 (en) | Image enhancement method and apparatus | |
| CN110175959B (zh) | 一种台风云图增强方法 | |
| CN106803236A (zh) | 基于模糊域奇异值分解的非均匀性校正方法 | |
| CN110706180B (zh) | 一种极暗图像视觉质量提升方法、系统、设备及介质 | |
| CN117011192A (zh) | 基于中性集合与双边滤波的nsct域声呐图像去噪方法及设备 | |
| CN110111268B (zh) | 基于暗通道和模糊宽度学习的单幅图像去雨方法及装置 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |