CN107292851A - 一种基于伪3d变换的bm3d图像降噪方法 - Google Patents
一种基于伪3d变换的bm3d图像降噪方法 Download PDFInfo
- Publication number
- CN107292851A CN107292851A CN201710576072.4A CN201710576072A CN107292851A CN 107292851 A CN107292851 A CN 107292851A CN 201710576072 A CN201710576072 A CN 201710576072A CN 107292851 A CN107292851 A CN 107292851A
- Authority
- CN
- China
- Prior art keywords
- image
- block
- coefficient
- bm3d
- pseudo
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 66
- 230000009467 reduction Effects 0.000 claims abstract description 66
- 238000006243 chemical reaction Methods 0.000 claims abstract description 54
- 238000012545 processing Methods 0.000 claims abstract description 44
- 230000009466 transformation Effects 0.000 claims description 19
- 238000001914 filtration Methods 0.000 claims description 17
- 238000011946 reduction process Methods 0.000 claims description 10
- 238000004364 calculation method Methods 0.000 abstract description 12
- 230000000694 effects Effects 0.000 description 6
- 238000005516 engineering process Methods 0.000 description 5
- 238000012360 testing method Methods 0.000 description 5
- 238000012544 monitoring process Methods 0.000 description 3
- 238000005457 optimization Methods 0.000 description 3
- 230000007812 deficiency Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000003384 imaging method Methods 0.000 description 2
- 230000007704 transition Effects 0.000 description 2
- 230000000007 visual effect Effects 0.000 description 2
- 241000287196 Asthenes Species 0.000 description 1
- 239000000654 additive Substances 0.000 description 1
- 230000000996 additive effect Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000003333 near-infrared imaging Methods 0.000 description 1
- 238000007781 pre-processing Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000012876 topography Methods 0.000 description 1
- 238000013519 translation Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/70—Denoising; Smoothing
-
- 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 using non-spatial domain filtering
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/11—Region-based segmentation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/30—Determination of transform parameters for the alignment of images, i.e. image registration
- G06T7/37—Determination of transform parameters for the alignment of images, i.e. image registration using transform domain methods
-
- 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]
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Image Processing (AREA)
- Image Analysis (AREA)
Abstract
本发明公开了一种基于伪3D变换的BM3D图像降噪方法,包括以下步骤:将待降噪图像分割成若干个图像块;以图像块为中心,在预设窗口范围内搜索相似图像块,并用图像块及其相似图像块构建三维图像块阵列;在阵列第三个维度方向进行频域变换处理,得到一维变换系数;对一维变换系数中的最低频系数进行二维正变换处理,得到伪变换系数;对伪变换系数进行滤波处理;对滤波后的伪变换系数进行二维逆变换处理,得到处理后的一维变换系数;在阵列第三个维度方向进行频域逆变换处理,得到降噪图像。该降噪方法在保证降噪效果的情况下,解决了BM3D图像降噪中计算量巨大、计算耗时长等问题,提高了计算效率。
Description
技术领域
本发明涉及BM3D图像降噪方法,具体涉及一种基于伪3D变换的BM3D图像降噪方法。
背景技术
随着智慧城市、大数据等领域的发展,基于视频监控的智能交通、车牌自动识别处理技术也逐渐成为研究的热门课题,并得到广泛应用。视频作为信息最为丰富的信息源,其质量直接影响到监控系统的效果,并且在交通领域,视频中的车辆、车牌等信息有时候直接作为司法取证的有效辅助证据。由于交通领域的视频监控通常处于室外环境,由于受到复杂多变的成像环境,以及成像设备硬件等因素的干扰,视频图像都存在一定的噪声。尤其是光线不足情况下、或者长焦摄像头的近红外成像,其视频图像噪声特别强,导致视频质量严重衰减。因此,监控视频通常需要对其进行图像降噪处理,以改善图像质量、尽可能的保留图像中的有用信息。
图像降噪是图像处理领域中一项基本而又是十分关键的技术,一直是图像处理领域的一个难题,属于图像预处理的关键步骤。图像去噪目的是尽可能恢复图像原貌,改善图像质量,突出图像本身的特征,从而为后续图像处理的其它工作打下良好的基础。能够有效去除噪声直接影响到图像的视觉效果和后续的图像处理结果。
图像降噪方法主要有经典空域滤波方法、频域滤波方法、多尺度几何分析降噪方法、非局部降噪方法以及由非局部均值降噪方法改进的各种降噪方法和盲图像去噪方法。空域滤波方法一般先假定噪声为加性模型,然后利用图像像素之间的的空间相关性对噪声进行滤波,如均值滤波和中值滤波。传统的空域滤波处理具有良好的去噪性能,并且方法简单易于实现,从而得到了很好的应用。然而这些经典去噪方法具有一定的局限性,它们都是基于局部统计的自适应滤波,因此对图像的同质区域有着较好的去噪效果,在消除图像噪声的同时会平滑图像的纹理或者不能有效的去除边缘周围的噪声,图像结构信息、点线等目标会在一定程度上被模糊或滤除。近来,由Buades等人提出的非局部均值滤波算法打破了局部邻域这一僵局,将传统的局部特征统计扩展到了非局部区域。局部算法通过计算像素点各自邻域间的欧氏距离,来判断像素点间的相似性,而没有对用于去除当前像素点噪声的大部分像素做任何假设。像素点间的欧氏距离可以用于判断邻域内的像素点与中心像素点的相似程度,并通过欧氏距离计算相似性权值,利用权值以及相似点对中心像素点加权平均得到去噪结果。非局部均值算法,将去噪从局部扩展到了非局部领域,提供了一个新的去噪思路。单个像素点的灰度值和空间距离不足以衡量像素点之间的相似性。非局部算法基本思想是以一个窗口内的图像块为处理单元,通过计算图像块之间的欧式距离,并以此衡量图像块之间的相似程度。由于大部分图像存在空间冗余,表现在图像块与图像块之间存在一定的相似性,以图像块的相似性代替像素点的相似性更科学,更具有可信性。基于非局部思想先后有非局部均值降噪、贝叶斯非局部降噪、三维块匹配算法(BM3D)降噪等方法被提出,并且取得了良好的降噪效果。其中BM3D算法是公认的降噪效果最佳的图像降噪算法。
基于三维变换的BM3D降噪方法如下:将待降噪图像分割成图像块,由图像块构建三维图像块阵列;对每个图像块作二维正变换,得到对应的二维变换系数;然后对图像块第三个维度作一维变换,得到对应的一维变换系数;由三维图像块阵列的二维变换系数和一维变换系数得到三维图像块阵列的三维变换系数;对三维变换系数滤波降噪处理;对滤波处理后的三维变换系数的第三个维度作一维逆变换,得到一维逆变换后的二维变换系数;对一维逆变换后的二维变换系数作二维逆变换,得到降噪后的图像块;最后得到降噪后的图像。BM3D图像降噪方法仍存在很多不足,其降噪的计算量巨大、计算耗时长,目前主要用于单帧图像降噪算法仿真,还无法对多帧视频进行实时降噪。
发明内容
针对现有技术存在的上述不足,本发明要解决的技术问题是:提供一种基于伪3D变换的BM3D图像降噪方法,在保证降噪效果的情况下,解决BM3D图像降噪过程中计算量巨大、计算耗时长等问题,可以提高计算效率、减少运算量。
为解决上述技术问题,本发明采用了如下的技术方案:
一种基于伪3D变换的BM3D图像降噪方法,包括以下步骤:
S1:对待降噪图像进行分割,得到若干个图像块。可以对图像块的大小进行预设,图像块大小通常设置为2的幂级数,以8*8为宜,图像块设置过小会导致频域变换效果不佳,图像块设置过大会极大增加图像块匹配计算量。
S2:以一个图像块为中心,在预设窗口范围内依据相似性准则搜索相似图像块,以图像块及其相似图像块构建三维图像块阵列;以三维图像块阵列中图像块的横向和纵向分别作为三维图像块阵列的第一个维度和第二个维度,以相似图像块对应点间的连接线作为第三个维度。根据背景技术中介绍的非局部算法基本思想,将一个窗口内的图像块作为一个处理单元,即所述的预设窗口范围。三维图像块阵列的构建,选择预设窗口的中心图像块作为三维图像块阵列中的第一个图像块,图像块的相似图像块按照搜索顺序排列。
S3:在三维图像块阵列第三个维度方向对每一组对应点分别进行频域变换处理,得到一维变换系数。在三维图像块阵列的第三个维度方向上,有多组相似图像块,对每一组相似图像块的对应点分别进行频域变换,得到多个频域系数,构成一维变换系数块阵列,即所述的一维变换系数。
S4:对一维变换系数块阵列中的最低频系数块进行二维正变换处理,得到BM3D伪变换系数;对BM3D伪变换系数进行滤波处理。一维变换系数块阵列中,第一个系数块对应为最低频系数块。仅对一维变换系数块阵列中的最低频系数块作二维正变换,对其他的频域系数不作变换处理,从而得到BM3D伪变换系数。通过对BM3D伪变换系数滤波处理来实现对图像块的降噪。
S5:对一维变换系数中最低频系数块对应的BM3D伪变换系数进行二维逆变换处理,得到处理后的一维变换系数。
S6:在三维图像块阵列第三个维度方向对每一组处理后的一维变换系数进行频域逆变换处理,得到降噪后的三维图像块阵列,从降噪后的三维图像块阵列中提取出降噪后的图像块,此处降噪后的图像块是指步骤S1中在预设窗口范围内处于中间位置的图像块经过上述降噪处理后的图像块;降噪后的三维图像块阵列中,除降噪后的图像块外,其余均为降噪后的图像块的相似图像块。
S7:循环步骤S2到步骤S6,对每个图像块分别进行降噪处理;对降噪处理后的若干图像块的像素点进行除权处理后,得到降噪图像。
进一步地,作为优化,步骤S3中的频域变换处理和步骤S6中的频域逆变换处理均采用哈尔小波变换(haar变换)、快速傅里叶变换(fft变换)或者变换编码(hadamard变换)。
进一步地,作为优化,步骤S4中的滤波处理采用硬阈值滤波和维纳滤波。
进一步地,作为优化,步骤S4的二维正变换处理和步骤S5中的二维逆变换处理均采用双正交小波或者离散余弦变换(DCT)进行二维变换。
本申请中伪3D变换与3D变换的分析如下:
3D变换是先对三维图像块阵列的三个维度均做正变换处理,得到三维变换系数,对三维变换系数进行滤波降噪,对滤波后的三维变换系数在三个维度均做逆变换处理,从而达到图像降噪的目的。本申请中的伪3D变换仅对三维图像块阵列的第三个维度方向做变换处理,先对三维图像块阵列的第三个维度做频域变换处理,得到一维变换系数;将得到的一维变换系数存放在原图像块位置,得到频域变换系数块;再对频域变换系数块中的最低频系数块作二维正变换,对其他的频域系数块不作变换处理,得到伪三维变换系数,对伪三维变换系数进行滤波降噪;对最低频系数块对应的滤波后伪三维变换系数进行二维逆变换处理,得到处理后的一维变换系数块;在三维图像块阵列第三个维度方向进行频域逆变换处理,得到降噪后的三维图像块阵列,从而达到图像降噪的目的。
伪3D变换仅对三维图像块阵列的第三个维度方向做变换处理,同时,也仅对一维变换系数中的最低频系数进行二维变换处理,对一维变换系数中的其他系数不做变换处理;逆变换亦只对三维图像块阵列的第三个维度方向做变换处理,只对一维变换系数中最低频系数对应的BM3D伪变换系数作变换;这样就大大地降低了三维变换的计算量,亦可以达到与3D变换几乎相同的降噪效果。
附图说明
图1为本发明实施例中基于伪3D变换的BM3D图像降噪方法的流程图;
图2为本发明实施例中基于伪3D变换的BM3D图像降噪方法的计算示意图;
图3为本发明实施例中不带噪声的原图;
图4为本发明实施例中带噪声的待降噪图像;
图5为本发明实施例中伪3D变换的BM3D图像降噪方法的第一次降噪效果图;
图6为本发明实施例中3D变换的BM3D图像降噪方法的第一次降噪效果图;
图7为本发明实施例中伪3D变换的BM3D图像降噪方法的第二次降噪效果图;
图8为本发明实施例中3D变换的BM3D图像降噪方法的第二次降噪效果图。
具体实施方式
下面将结合附图对本发明技术方案的实施例进行详细的描述。以下实施例仅用于更加清楚地说明本发明的技术方案,因此只是作为示例,而不能以此来限制本发明的保护范围。
实施例:
参照图1~图8,一种基于伪3D变换的BM3D图像降噪方法,包括以下步骤:
S1:对待降噪图像进行分割,得到若干个图像块。可以对图像块的大小进行预设,图像块大小通常设置为2的幂级数,以8*8为宜,图像块设置过小会导致频域变换效果不佳,图像块设置过大会极大增加图像块匹配计算量。
S2:以一个图像块为中心,在预设窗口范围内依据相似性准则搜索相似图像块,以图像块及其相似图像块构建三维图像块阵列;以三维图像块阵列中图像块的横向和纵向分别作为三维图像块阵列的第一个维度和第二个维度,以相似图像块对应点间的连接线作为第三个维度。根据背景技术中介绍的非局部算法基本思想,将一个窗口内的图像块作为一个处理单元,即所述的预设窗口范围。三维图像块阵列的构建,选择预设窗口的中心图像块作为三维图像块阵列中的第一个图像块,图像块的相似图像块按照搜索顺序排列。
S3:在三维图像块阵列第三个维度方向对每一组对应点分别进行频域变换处理,得到一维变换系数。在三维图像块阵列的第三个维度方向上,有多组相似图像块,对每一组相似图像块的对应点分别进行频域变换,得到多个频域系数,构成一维变换系数块阵列,即所述的一维变换系数。
S4:对一维变换系数块阵列中的最低频系数块进行二维正变换处理,得到BM3D伪变换系数;对BM3D伪变换系数进行滤波处理。一维变换系数块阵列中,第一个系数块对应为最低频系数块。仅对一维变换系数块阵列中的最低频系数块作二维正变换,对其他的频域系数不作变换处理,从而得到BM3D伪变换系数。通过对BM3D伪变换系数滤波处理来实现对图像块的降噪。
S5:对一维变换系数中最低频系数块对应的BM3D伪变换系数进行二维逆变换处理,得到处理后的一维变换系数。
S6:在三维图像块阵列第三个维度方向对每一组处理后的一维变换系数进行频域逆变换处理,得到降噪后的三维图像块阵列,从降噪后的三维图像块阵列中提取出降噪后的图像块,此处降噪后的图像块是指步骤S1中在预设窗口范围内处于中间位置的图像块经过上述降噪处理后的图像块;降噪后的三维图像块阵列中,除降噪后的图像块外,其余均为降噪后的图像块的相似图像块。
S7:循环步骤S2到步骤S6,对每个图像块分别进行降噪处理;对降噪处理后的若干图像块的像素点进行除权处理后,得到降噪图像。
本实施例中,步骤S3中的频域变换处理和步骤S6中的频域逆变换处理均采用哈尔小波变换(haar变换)、快速傅里叶变换(fft变换)或者变换编码(hadamard变换)。
本实施例中,步骤S4中的滤波处理采用硬阈值滤波和维纳滤波。
本实施例中,步骤S4的二维正变换处理和步骤S5中的二维逆变换处理均采用双正交小波或者离散余弦变换(DCT)进行二维变换。
参照图像2,为本发明基于伪3D变换的BM3D图像降噪方法的计算示意图。对待降噪图像进行分割,得到若干个图像块;以一个图像块为中心,在预设窗口范围内依据相似性准则搜索相似图像块,以图像块及其相似图像块构建三维图像块阵列;对三维图像块阵列的第三个维度做频域变换处理,得到一维变换系数,将得到的一维变换系数存放在原图像块位置,得到频域变换系数块;再对频域变换系数块中的最低频系数块作二维正变换,对其他的频域系数块不作变换处理,得到伪三维变换系数,对伪三维变换系数进行滤波降噪;对最低频系数块对应的伪三维变换系数进行二维逆变换处理,得到处理后的一维变换系数块;在三维图像块阵列第三个维度方向进行频域逆变换处理,得到降噪后的三维图像块阵列,从而达到图像降噪的目的。
参照图3到图8,本发明基于伪3D变换BM3D算法与基于3D变换BM3D算法降噪效果比较:采用ISO12233测试标板局部图像作为原始测试图像,图像大小为512*512,并添加标准差为25的高斯噪声;分别采用基于伪3D变换BM3D算法和基于3D变换BM3D算法进行降噪对比测试;所有参数设置均相同。
附图3是不含噪声的ISO12233测试标板图像;附图4是添加了高斯噪声的ISO12233测试标板图像,信噪比为20.18dB;附图5是通过伪3D变换BM3D算法的第一次降噪效果图,信噪比为33.37dB;附图6是通过3D变换BM3D算法的第一次降噪效果图,信噪比为33.36dB;附图7是通过伪3D变换BM3D算法的第二次降噪效果图,信噪比为34.12dB;附图8是通过3D变换BM3D算法的第二次降噪效果图,信噪比为34.38dB。无论是从视觉效果看,还是测量信噪比数值比较,本发明基于伪3D变换BM3D算法与基于3D变换BM3D算法的降噪效果几乎相同。
假设每一个图像块的大小为M×N,相似图像块的个数为K,相似图像块间的对应点组数为a。a组数据点的一维变换的计算量记作则一个M×N×K的三维图像块的3D变换的三维正变换和逆变换的计算量(p1)为:
采用伪3D变换的三维正变换和逆变换的计算量(p2)为:
假设M=8,N=8,K=16,则p1=20480,p2=8960
可以看出:基于伪3D变换BM3D算法与基于3D变换BM3D算法的降噪效果几乎相同,同时,每一个相似块的伪3D变换BM3D算法的计算量大约只有3D变换BM3D算法的三分之一,计算量降低了约三分之二,大幅度降低了BM3D算法的计算量,减少了运算量、提高运算效率。
本实施例中,基于伪3D变换的BM3D图像降噪方法的实施过程如下:
假设三维图像块阵列(T)为:
T∈Rm×n×k为三维图像块阵列,每个图像块为m行n列,共有k个相似图像块。
Ij表示第j个图像块,Ik表述第k个图像块,表示第j个图像块中的p行q列,表示k个相似图像块的p行q列像素点构成的向量。对三维图像块阵列T的第三个维度进行频域变换得到一维变换系数SZ:
其中,f(Γpq)表示对向量Γpq作频域变换,向量Γpq变换后的系数为频域变换系数Φpq,是一维变换后频域变换系数块阵列中的一个向量。
是一维变换后系数块阵列中的一个系数块。
对频域变换系数块阵列的最低频系数块采用Bior1.5双正交小波或者DCT进行二维变换,变换后的系数为
B(Φ1)表示对Φ1进行二维变换,其结果为:
为伪3D变换中的一个三维系数块,其它系数块Φj均为一维变换后对应的系数块,即伪3D变换后的系数并非是都是进行了三维变换的系数,故称为伪3D变换。
对伪3D变换后的系数进行滤波,结果表示为
表示对系数进行滤波,表示对系数进行滤波,滤波方式为硬阈值滤波或者维纳滤波。表示的滤波结果,表示对的滤波结果,表示滤波后的第一个系数块,表示滤波后的第j个系数块。
对滤波后的第一个系数块(即最低频系数块)进行二维逆变换,其结果为
表示对作二维逆变换,结果为:
表示二维逆变换后的系数块阵列中的一个向量。
对进行逆变换,其结果为
表示对作一维逆变换,其结果为表示一维逆变换后所得降噪图像块阵列中的一个向量,表示一维逆变换后所得降噪后的三维图像块阵列中的降噪后的图像块。
循环上述公式算法,对分割后的若干图像块进行降噪,得到若干个降噪后的图像块,对降噪处理后的若干图像块的像素点进行除权处理后,将降噪后的图像块拼接融合成降噪图像。进行除权处理的目的,是由于原图像中每个图像块的像素点在上述循环运算过程中可能被计算过多次,将每次的处理结果相加除以每次计算的权重系数即为除权处理,以获得每个像素点多次计算的最优估值。
现有技术中,BM3D图像降噪方法计算量大、计算时间长,无法进行硬件实时处理问题;本发明针对这一现状,提出了一种BM3D图像降噪的加速计算方法。针对三维正变换和逆变换计算量巨大问题,提出一种基于三维块匹配的不完全3D变换方法,该方法只需先对三维相似块的第三维度进行变换,并对上述变换后所得系数的最低频系数块进行二维变换,而无需对所有的系数块进行二维变换,从而将3D变换过程的计算量降低了约三分之二,降噪效果与完全3D变换降噪效果基本一致。
最后说明的是,以上实施例仅用以说明本发明的技术方案而非限制,尽管参照实施例对本发明进行了详细说明,本领域的普通技术人员应当理解,可以对本发明的技术方案进行修改或者等同替换,而不脱离本发明技术方案的宗旨和范围,其均应涵盖在本发明当中。
Claims (4)
1.一种基于伪3D变换的BM3D图像降噪方法,其特征在于,包括以下步骤:
S1:对待降噪图像进行分割,得到若干个图像块;
S2:以一个图像块为中心,在预设窗口范围内依据相似性准则搜索相似图像块,以图像块及其相似图像块构建三维图像块阵列;以三维图像块阵列中图像块的横向和纵向分别作为三维图像块阵列的第一个维度和第二个维度,以相似图像块对应点间的连接线作为第三个维度;
S3:在三维图像块阵列第三个维度方向对每一组对应点分别进行频域变换处理,得到一维变换系数;
S4:对一维变换系数中的最低频系数进行二维正变换处理,得到BM3D伪变换系数;对BM3D伪变换系数进行滤波处理;
S5:对一维变换系数中最低频系数对应的滤波后BM3D伪变换系数进行二维逆变换处理,得到处理后的一维变换系数;
S6:在三维图像块阵列第三个维度方向对每一组处理后的一维变换系数进行频域逆变换处理,得到降噪后的三维图像块阵列,从降噪后的三维图像块阵列中提取出降噪后的图像块;
S7:循环步骤S2到步骤S6,对每个图像块分别进行降噪处理;对降噪处理后的若干图像块的像素点进行除权处理后,得到降噪图像。
2.根据权利要求1所述的基于伪3D变换的BM3D图像降噪方法,其特征在于,步骤S3中的频域变换处理和步骤S6中的频域逆变换处理均采用哈尔小波变换、快速傅里叶变换或者变换编码。
3.根据权利要求1所述的基于伪3D变换的BM3D图像降噪方法,其特征在于,步骤S4中的滤波处理采用硬阈值滤波和维纳滤波。
4.根据权利要求1所述的基于伪3D变换的BM3D图像降噪方法,其特征在于,步骤S4的二维正变换处理和步骤S5中的二维逆变换处理均采用双正交小波或者离散余弦变换进行二维变换。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710576072.4A CN107292851B (zh) | 2017-07-14 | 2017-07-14 | 一种基于伪3d变换的bm3d图像降噪方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710576072.4A CN107292851B (zh) | 2017-07-14 | 2017-07-14 | 一种基于伪3d变换的bm3d图像降噪方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107292851A true CN107292851A (zh) | 2017-10-24 |
CN107292851B CN107292851B (zh) | 2019-09-27 |
Family
ID=60101804
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710576072.4A Expired - Fee Related CN107292851B (zh) | 2017-07-14 | 2017-07-14 | 一种基于伪3d变换的bm3d图像降噪方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107292851B (zh) |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109493295A (zh) * | 2018-10-31 | 2019-03-19 | 泰山学院 | 一种非局部哈尔变换图像去噪方法 |
CN109492652A (zh) * | 2018-11-12 | 2019-03-19 | 重庆理工大学 | 一种基于有序视觉特征单词库模型的相似图像判断方法 |
CN110060220A (zh) * | 2019-04-26 | 2019-07-26 | 中国科学院长春光学精密机械与物理研究所 | 基于改进bm3d算法的图像去噪方法及系统 |
WO2019191892A1 (zh) * | 2018-04-02 | 2019-10-10 | 北京大学 | 用于视频编解码的方法和设备 |
CN111754393A (zh) * | 2020-06-28 | 2020-10-09 | 展讯通信(上海)有限公司 | 图像处理方法、系统、电子设备和介质 |
CN112053295A (zh) * | 2020-08-21 | 2020-12-08 | 珠海市杰理科技股份有限公司 | 图像降噪方法、装置、计算机设备和存储介质 |
CN112188214A (zh) * | 2020-09-22 | 2021-01-05 | 展讯通信(上海)有限公司 | 图像处理方法、系统、电子设备和介质 |
CN112862723A (zh) * | 2021-03-04 | 2021-05-28 | 清华大学深圳国际研究生院 | 基于伪3d自相关网络的真实图像去噪方法 |
CN112991190A (zh) * | 2019-12-16 | 2021-06-18 | 荣耀终端有限公司 | 图像降噪方法及其介质和电子设备 |
WO2021135702A1 (zh) * | 2019-12-31 | 2021-07-08 | 百果园技术(新加坡)有限公司 | 一种视频去噪方法和电子设备 |
CN113643210A (zh) * | 2021-08-26 | 2021-11-12 | Oppo广东移动通信有限公司 | 图像处理方法、装置、电子设备和存储介质 |
CN114088658A (zh) * | 2021-10-09 | 2022-02-25 | 池明旻 | 用于近红外织物纤维成分无损清洁分析的降噪处理方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102184526A (zh) * | 2011-04-22 | 2011-09-14 | 西安电子科技大学 | 基于字典学习和块匹配的自然图像去噪方法 |
CN102663702A (zh) * | 2012-04-20 | 2012-09-12 | 西安电子科技大学 | 基于区域划分的自然图像去噪方法 |
CN104751418A (zh) * | 2013-12-26 | 2015-07-01 | 联发科技股份有限公司 | 3维区块匹配的图像去噪的方法与装置 |
-
2017
- 2017-07-14 CN CN201710576072.4A patent/CN107292851B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102184526A (zh) * | 2011-04-22 | 2011-09-14 | 西安电子科技大学 | 基于字典学习和块匹配的自然图像去噪方法 |
CN102663702A (zh) * | 2012-04-20 | 2012-09-12 | 西安电子科技大学 | 基于区域划分的自然图像去噪方法 |
CN104751418A (zh) * | 2013-12-26 | 2015-07-01 | 联发科技股份有限公司 | 3维区块匹配的图像去噪的方法与装置 |
Non-Patent Citations (3)
Title |
---|
AAJJAD AMANI ET AL.: "Seismic random noise attenuation via 3D block matchiing", 《JOURANL OF APPLIED GEOPHYSICS》 * |
KOSTADIN DABOV ET AL.: "Image denoising by sparse 3D transform-domain collaborative filtering", 《IEEE TRANSACTIONS ON IMAGE PROCESSING》 * |
YINGKUN HOU ET AL.: "Image Denoising by Block-matching and 1D Filtering", 《FOURTH INTERNATIONAL CONFERENCE ON MACHINE VISION(ICMV 2011)》 * |
Cited By (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2019191892A1 (zh) * | 2018-04-02 | 2019-10-10 | 北京大学 | 用于视频编解码的方法和设备 |
US11394966B2 (en) | 2018-04-02 | 2022-07-19 | SZ DJI Technology Co., Ltd. | Video encoding and decoding method and apparatus |
CN109493295A (zh) * | 2018-10-31 | 2019-03-19 | 泰山学院 | 一种非局部哈尔变换图像去噪方法 |
CN109493295B (zh) * | 2018-10-31 | 2022-02-11 | 泰山学院 | 一种非局部哈尔变换图像去噪方法 |
CN109492652A (zh) * | 2018-11-12 | 2019-03-19 | 重庆理工大学 | 一种基于有序视觉特征单词库模型的相似图像判断方法 |
CN109492652B (zh) * | 2018-11-12 | 2021-10-15 | 重庆理工大学 | 一种基于有序视觉特征单词库模型的相似图像判断方法 |
CN110060220A (zh) * | 2019-04-26 | 2019-07-26 | 中国科学院长春光学精密机械与物理研究所 | 基于改进bm3d算法的图像去噪方法及系统 |
CN112991190A (zh) * | 2019-12-16 | 2021-06-18 | 荣耀终端有限公司 | 图像降噪方法及其介质和电子设备 |
WO2021135702A1 (zh) * | 2019-12-31 | 2021-07-08 | 百果园技术(新加坡)有限公司 | 一种视频去噪方法和电子设备 |
CN111754393A (zh) * | 2020-06-28 | 2020-10-09 | 展讯通信(上海)有限公司 | 图像处理方法、系统、电子设备和介质 |
CN112053295A (zh) * | 2020-08-21 | 2020-12-08 | 珠海市杰理科技股份有限公司 | 图像降噪方法、装置、计算机设备和存储介质 |
CN112053295B (zh) * | 2020-08-21 | 2024-04-05 | 珠海市杰理科技股份有限公司 | 图像降噪方法、装置、计算机设备和存储介质 |
CN112188214A (zh) * | 2020-09-22 | 2021-01-05 | 展讯通信(上海)有限公司 | 图像处理方法、系统、电子设备和介质 |
CN112862723A (zh) * | 2021-03-04 | 2021-05-28 | 清华大学深圳国际研究生院 | 基于伪3d自相关网络的真实图像去噪方法 |
CN112862723B (zh) * | 2021-03-04 | 2022-08-02 | 清华大学深圳国际研究生院 | 基于伪3d自相关网络的真实图像去噪方法 |
CN113643210A (zh) * | 2021-08-26 | 2021-11-12 | Oppo广东移动通信有限公司 | 图像处理方法、装置、电子设备和存储介质 |
CN114088658A (zh) * | 2021-10-09 | 2022-02-25 | 池明旻 | 用于近红外织物纤维成分无损清洁分析的降噪处理方法 |
Also Published As
Publication number | Publication date |
---|---|
CN107292851B (zh) | 2019-09-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107292851B (zh) | 一种基于伪3d变换的bm3d图像降噪方法 | |
JP6759475B2 (ja) | シーンの多次元特徴に基づく船舶検出方法及びシステム | |
CN105741252B (zh) | 基于稀疏表示与字典学习的视频图像分级重建方法 | |
Yue et al. | Image denoising by exploring external and internal correlations | |
CN103049892B (zh) | 基于相似块矩阵秩最小化的非局部图像去噪方法 | |
CN110163818A (zh) | 一种用于海事无人机的低照度视频图像增强方法 | |
CN112215074A (zh) | 基于无人机视觉的实时目标识别与检测追踪系统及方法 | |
CN106485671A (zh) | 基于边缘的多方向加权tv和自相似性约束图像去模糊方法 | |
CN109509164A (zh) | 一种基于gdgf的多传感器图像融合方法及系统 | |
CN109509163B (zh) | 一种基于fgf的多聚焦图像融合方法及系统 | |
CN102243711A (zh) | 一种基于邻域嵌套的图像超分辨率重建方法 | |
CN110097617B (zh) | 基于卷积神经网络与显著性权重的图像融合方法 | |
CN107590779A (zh) | 一种基于图像块聚类字典训练的图像去噪去模糊方法 | |
CN103020933A (zh) | 一种基于仿生视觉机理的多源图像融合方法 | |
CN111899200B (zh) | 一种基于3d滤波的红外图像增强方法 | |
CN110889868B (zh) | 一种结合梯度和纹理特征的单目图像深度估计方法 | |
CN110415169A (zh) | 一种深度图超分辨率重建方法、系统及电子设备 | |
Qiu et al. | A GAN-based motion blurred image restoration algorithm | |
CN116433500A (zh) | 一种基于可逆网络的图像去雨方法、系统和可读存储介质 | |
CN109544477A (zh) | 基于自适应字典学习稀疏表示的图像去噪算法 | |
CN108198140A (zh) | 基于ncsr模型的三维协同滤波去噪方法 | |
Gang et al. | Infrared and visible dynamic image sequence fusion based on region target detection | |
CN101127120A (zh) | 一种有效抑制模板漂移的目标跟踪算法 | |
CN104601861B (zh) | 一种针对光纤监控视频序列的降噪方法及系统 | |
Liu et al. | Comparison and analysis of classical image denoising methods based on convolution neural network |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20190927 Termination date: 20200714 |
|
CF01 | Termination of patent right due to non-payment of annual fee |