CN101719268B - 基于改进Directionlet域的广义高斯模型图像去噪方法 - Google Patents
基于改进Directionlet域的广义高斯模型图像去噪方法 Download PDFInfo
- Publication number
- CN101719268B CN101719268B CN2009102193489A CN200910219348A CN101719268B CN 101719268 B CN101719268 B CN 101719268B CN 2009102193489 A CN2009102193489 A CN 2009102193489A CN 200910219348 A CN200910219348 A CN 200910219348A CN 101719268 B CN101719268 B CN 101719268B
- Authority
- CN
- China
- Prior art keywords
- coefficient
- noise
- subgraph
- band
- frequency sub
- 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
Images
Landscapes
- Image Processing (AREA)
Abstract
本发明公开了一种基于改进Directionlet域的广义高斯模型图像去噪方法,主要解决现有去噪方法边缘细节丢失严重、均匀区域过平滑的问题。其实现步骤为:(1)选测试图像,加入高斯噪声,得到噪声图像;(2)对噪声图像进行子图分割,确定各子图的变换矩阵;(3)对子图采样,得到陪集;(4)对各陪集进行各向异性小波变换;(5)估计高频子带广义高斯模型的形状参数和局部标准差;(6)由含噪系数估计无噪系数;(7)对无噪系数进行各向异性小波逆变换;(8)根据变换矩阵加权综合,重构各子图;(9)将重构的子图合成,得到去噪结果。本发明具有边缘细节保持好、均匀区域失真少和峰值信噪比高的优点,可用于去除自然图像中的高斯噪声。
Description
技术领域
本发明属于数字图像处理技术领域,具体地说是一种图像去噪方法,该方法可用于去除被污染的自然图像中的噪声。
背景技术
图像在获取和传输过程中,常会受到噪声的污染。如何对图像进行去噪并尽可能多地保留图像的细节信息以改善图像的质量,成为图像处理中的一个重要任务。人们根据图像的特点、噪声的统计特征和频谱的分布规律,发展了多种去噪方法,其中最为直观的是根据噪声一般集中在高频,而图像频谱则分布在一个有限区间这一特点,采用低通滤波来进行去噪,如滑动平均窗滤波,Wiener滤波等。其他的去噪方法还有基于秩-阶滤波的方法,基于马尔可夫场模型的方法等,这些方法仅具有空域或频域的局部分析性能。
近年来,由于小波具有良好的时频特性和多分辨特性,小波阈值方法被广泛应用到各种去噪处理中,但对经过正交小波变换后的高频系数进行阈值去噪时,会出现“过扼杀”,去噪图像边缘处会出现振荡现象。在小波基础上发展而来的平移不变的平稳小波变换,克服了正交小波变换去噪时的不足,但在阈值去噪时均匀区域出现过平滑现象。为解决二维或更高维奇异性而出现的多尺度几何工具Contourlet,对图像中的奇异曲线的逼近接近最优,Contourlet阈值去噪时边缘信息丢失较少,但在均匀区域会产生虚假成分,即蚊状噪声很明显。
Vladan Velisavljevi′提出的Directionlet变换,是一种新的多尺度几何分析工具。当Directionlet基函数的方向与图像中各向异性目标的方向匹配时,对图像的逼近效果较好,可以稀疏的表示图像的边缘等细节信息;不匹配时Directionlet则退化为小波。但现有的Directionlet变换,变换方向和队列方向都是任意选取的,比如0度,90度,45度和-45度,并没有根据图像的特征进行方向自适应的变换,导致在图像去噪时边缘信息丢失严重,出现模糊失真现象。
发明内容
本发明的目的在于克服上述已有技术的不足,提供一种基于改进Directionlet域的广义高斯模型图像去噪方法,以根据图像的特征进行方向自适应的变换,实现在去除噪声的同时,尽可能保持图像边缘等细节特征的目的,提高图像的质量。
实现本发明目的的技术方案是,首先自适应确定Directionlet变换的变换矩阵,使它根据图像的主要方向进行变换,在去噪处理时,用广义高斯模型对Directionlet系数建模,根据广义高斯模型形状参数的大小,结合局部标准差采取不同的无噪系数估计策略,对图像进行去噪处理,其具体去噪步骤包括如下:
(1)选取测试图像,加入零均值的高斯噪声,得到噪声图像;
(2)对噪声图像进行64×64的子图分割,并用二进小波变换自适应确定各分割子图的Directionlet变换矩阵MΛ;
(3)利用变换矩阵MΛ,对各分割子图进行采样,得到分割子图的|det(MΛ)|个陪集,|det(MΛ)|是矩阵MΛ行列式的绝对值;
(4)对各分割子图的每个陪集沿Directionlet变换矩阵MΛ的变换方向和队列方向分别进行n1=2与n2=1次的一维小波变换,得到Directionlet变换的高频和低频子带系数;
(5)对各个高频子带,利用该子带所有的变换系数,估计广义高斯模型的形状参数υ和局部标准差σx;
(6)对各高频子带的广义高斯模型形状参数υ进行判断:
若0<υ<0.5,按下式对噪声图像的无噪系数进行估计,
其中是无噪系数的估计,y是含噪系数,Tυ=Cυσ2/(2-υ)σx -(υ/(2-υ)),Cυ=(2-υ)(2-2υ)-(1-υ/2-υ)η(υ)υ/(2-υ), Γ是Gamma函数, o(y2(υ-1))是y2(υ-1)的高阶无穷小量,σ是噪声标准差,Y=(y1,y2,...,yN)是含噪高频子带系数,N是高频子带系数的数目;
若0.5≤υ<1,则采用阈值 进行软阈值处理;
(7)对低频子带和估计得到的无噪高频子带沿矩阵MΛ的变换方向和队列方向分别进行n1=2与n2=1次的一维小波逆变换;
(8)根据变换矩阵MΛ的变换方向和队列方向加权综合,重构各个分割子图;
(9)将重构的各分割子图按其在原图像中的位置合成,得到去噪后的图像。
本发明与现有技术相比具有如下优点:
(1)本发明采用Directionlet变换的矩阵是根据图像的特征自适应确定的,由于Directionlet变换系数对图像中的边缘各向异性特征逼近较好,对图像去噪时边缘细节信息保持较好,因而解决了现有Directionlet变换方法由于任意选取变换矩阵而退化成小波变换,在图像去噪处理中因对边缘等各向异性特征逼近不够精确,边缘信息丢失严重而引起的边缘模糊问题;
(2)本发明由于在无噪系数估计时充分考虑了广义高斯模型形状参数的影响,采用的局部标准差具有局部适应性,能更精确地估计高频子带的各个系数,解决了现有小波阈值方法对含噪系数的“过扼杀”问题和平稳小波阈值方法均匀区域过平滑问题,且本发明方法去噪后图像的均匀区域更清晰,失真较少。
仿真实验结果表明,本发明在对加有高斯噪声的自然图像去噪中,边缘等细节信息保持较好,均匀区域更清晰,得到了更高的峰值信噪比,提高了图像的质量。
附图说明
图1是本发明的主要操作过程示意图;
图2是用本发明方法与已有去噪方法对测试图像Lena的去噪效果对比图;
图3是用本发明方法与已有去噪方法对测试图像Barbara的去噪效果对比图;
图4是用本发明方法与已有去噪方法对测试图像Peppers的去噪效果对比图。
具体实施方式
参照图1,本发明的具体实现步骤如下:
步骤1:选取测试图像,加入高斯噪声,得到噪声图像,该噪声图像的数学表达式为:
f=g+n
其中g={g(i1,j1)|i1,j1=1,2,...N}表示测试图像,n={n(i1,j1)|i1,j1=1,2,...N}表示均值为零方差为σ2的高斯噪声,噪声图像记为f={f(i1,j1)|i1,j1=1,2,...N},N表示图像大小。
步骤2:把噪声图像分割成互不重叠,且大小为64×64的子图,寻找每个分割子图的两个主要方向θ1和θ2,根据主要方向θ1和θ2确定各分割子图的Directionlet变换矩阵MΛ,具体步骤如下:
2a)对分割子图进行二进小波变换,得到其水平细节图h(i,j)和垂直细节图v(i,j),其中(i,j)是二进小波变换系数的位置,i,j=1,2,...64;
2b)根据h(i,j)和v(i,j),计算分割子图在(i,j)处的方向θ(i,j):
若h(i,j)<<v(i,j),即
若h(i,j)>>v(i,j),即 θ(i,j)=0;
若 或
2c)统计分割子图的方向θ(i,j)的分布,找出出现次数最多的两个方向θ1和θ2;
2d)分别求两个主要方向θ1和θ2的反正切值,得到两个近似有理斜率r1和r2,r1≈arctanθ1=b1/a1,r2≈arctanθ2=b2/a2,根据有理斜率r1和r2,构造变换矩阵 其中沿r1的方向称为变换矩阵MΛ的变换方向,沿r2的方向称为队列方向,a1,a2,b1,b2都是整数。
步骤3:利用变换矩阵MΛ,对各分割子图分别沿r1和r2方向进行采样,得到分割子图的|det(MΛ)|个陪集,|det(MΛ)|是矩阵MΛ行列式的绝对值。
步骤4:对各分割子图的每个陪集沿r1方向进行n1=2次的一维小波变换,沿r2方向进行n2=1次的一维小波变换,得到Directionlet变换的高频和低频子带系数。
步骤5:由含噪高频子带系数估计广义高斯模型的形状参数υ和局部标准差σx按如下步骤进行:
3a)设高频子带的大小为l1×l2,由含噪的高频子带系数Y按下式计算出含噪高频子带系数的二阶矩σY 2和含噪高频子带系数的峰度kY:
σY 4=(σY 2)2
式中是(i2,j2)处含噪高频系数的平方,是(i2,j2)处含噪高频系数的4次方;
3b)由含噪高频子带系数Y,利用中值估计法求出噪声标准差: 再根据噪声标准差σ求出无噪系数的标准差:
3c)由以上得到的kY,σY 2和σ值,利用数值计算法按照如下公式求出形状参数υ:
3d)由含噪高频子带系数Y,利用最小均方误差估计法得到无噪系数的初始估计值
其中,是(i2,j2)处的含噪高频系数,i2=1,2,...l1,j2=1,2,...l2;
3e)由形状参数υ和无噪系数的初始估计值X,计算出局部标准差:σx=SE
其中S是只与形状参数υ有关的因子,
E是大小为5×5的局部邻域内无噪系数初始估计值的绝对均值,
i3=2,...l1-1,j3=2,...l2-1。
步骤6:将形状参数υ分成两个区间0<υ<0.5和0.5≤υ<1,使用以下估计策略由含噪系数估计无噪系数:
若0<υ<0.5,按下式对噪声图像的无噪系数进行估计,
其中是无噪系数的估计,y是含噪系数,Tυ=Cυσ2/(2-υ)σx -(υ/(2-υ)),Cυ=(2-υ)(2-2υ)-(1-υ/2-υ)η(υ)υ/(2-υ), Γ是Gamma函数, o(y2(υ-1))是y2(υ-1)的高阶无穷小量,σ是噪声标准差;
若0.5≤υ<1,则采用阈值 对含噪系数进行软阈值处理,得到无噪系数的估计。
步骤7:对低频子带和估计得到的无噪高频子带沿r1方向进行n1=2次的一维小波逆变换,沿r2方向进行n2=1次的一维小波逆变换。
步骤8:对步骤7得到的结果沿r1方向和r2方向进行加权综合,重构出各个分割子图。
步骤9:将重构的各分割子图按其在原图像中的位置合成,得到去噪后的图像。
本发明的效果通过以下仿真进一步说明。
一、仿真条件
采用图像去噪中常用的标准图像Lena 512×512、Barbara 512×512和Peppers512×512图像,并对这三幅图像分别加入标准差σ为10,15,20,25,30,35,40,45和50的零均值加性高斯白噪声,分别用本发明去噪方法、小波软阈值去噪方法、平稳小波硬阈值去噪方法、Contourlet硬阈值去噪方法和Contourlet软阈值去噪方法进行仿真去噪处理。
二、仿真结果分析
仿真1,用本发明方法与现有去噪方法对测试图像Lena的去噪效果对比图,结果如图2,其中:
图2(a)是测试图像Lena;图2(b)是σ=20的Lena含噪图像,峰值信噪比PSNR=22.13dB;图2(c)是小波软阈值的去噪结果图,PSNR=26.74dB;图2(d)是平稳小波硬阈值去噪结果图,PSNR=28.38dB;图2(e)是Contourlet硬阈值去噪结果图,PSNR=28.93dB;图2(f)Contourlet软阈值去噪结果图,PSNR=26.76dB;图2(g)是本发明方法去噪结果图,PSNR=30.19dB;
从图2(a)和图2(c)可以看出,现有小波软阈值方法去噪后Lena的帽沿处出现震荡现象,脸部和其它均匀区域失真严重。
从图2(d)可以看出,现有平稳小波硬阈值方法去噪后虽然Lena的帽沿比较清晰,但帽子上的纹理有些模糊,均匀区域出现过平滑。
从图2(e)和图2(f)可以看出,现有Contourlet阈值方法去噪后,产生了虚假成分,即蚊状噪声很明显。
从图2(g)可以看出本发明方法去噪后Lena的帽沿及帽子上的纹理比较清晰,脸部和其它均匀区域较平滑,失真较少,并且有较高的峰值信噪比,有效去除了图2(b)的噪声。
仿真2,用本发明方法与现有去噪方法对测试图像Barbara的去噪效果对比图,结果如图3,其中:
图3(a)是测试图像Barbara;图3(b)是σ=20的Barbara含噪图像,PSNR=22.19dB;图3(c)是小波软阈值方法的去噪结果图,PSNR=23.67dB;图3(d)是平稳小波硬阈值方法去噪结果图,PSNR=25.39dB;图3(e)是Contourlet硬阈值方法去噪结果图,PSNR=25.86dB;图3(f)Contourlet软阈值方法去噪结果图,PSNR=23.87dB;图3(g)本发明方法去噪结果图,PSNR=27.87dB。
从图3(a)和图3(c)可以看出,现有小波软阈值方法去噪后Barbara裤子和围巾上的纹理信息丢失严重。
从图5(d)可以看出,现有平稳小波硬阈值方法去噪后Barbara裤子和围巾上的纹理信息有些丢失,桌子腿部分出现震荡。
从图3(e)和图3(f)可以看出,现有Contourlet阈值去噪方法去噪后,Barbara纹理保持较好,图像较清晰,但产生了严重的虚假成分,即蚊状噪声很明显。
从图3(g)可以看出,本发明方法在对图3(b)去除噪声的同时更多的保持了Barbara裤子,围巾等处的纹理细节,可见本发明在细节保持方面有相当大的优势,均匀区域也更清晰。
仿真3,用本发明方法与现有去噪方法对测试图像Peppers的去噪效果对比图,结果如图4,其中:
图4(a)是测试图像Peppers;图4(b)是σ=20的Peppers含噪图像,PSNR=22.21dB;图4(c)是小波软阈值方法的去噪结果图,PSNR=27.21dB;图4(d)是平稳小波硬阈值方法去噪结果图,PSNR=28.59dB;图4(e)是Contourlet硬阈值方法去噪结果图,PSNR=28.65dB;图4(f)Contourlet软阈值方法去噪结果图,PSNR=26.51dB;图4(g)本发明方法去噪结果图,PSNR=30.82dB。
从图4(a)和图4(c)可以看出,现有小波软阈值方法去噪后Peppers边缘处出现震荡现象。
从图4(d)可以看出,现有平稳小波硬阈值方法去噪后,Peppers均匀区域有过平滑现象。
从图4(e)和图4(f)可以看出,现有Contourlet阈值去噪方法去噪较彻底,但蚊状噪声很明显。
从图4(g)可以看出,本发明方法对图4(b)去噪后各种蔬菜的边缘很清晰,均匀区域较清晰,和测试图像图4(a)相比差别更小,明显比其他方法去噪后的效果要好。
综上,本发明方法不仅视觉效果好,而且峰值信噪比也比较高,在边缘保持方面有较大优势。
Claims (2)
1.一种基于改进Directionlet域的广义高斯模型图像去噪方法,包括如下步骤:
(1)选取测试图像,加入零均值的高斯噪声,得到噪声图像;
(2)对噪声图像进行64×64的子图分割,用二进小波变换自适应确定各分割子图的Directionlet变换矩阵MΛ:
2a)对分割子图进行二进小波变换,得到其水平细节图h(i,j)和垂直细节图v(i,j),其中(i,j)是二进小波变换系数的位置,i,j=1,2,...64;
2b)根据h(i,j)和v(i,j),计算分割子图在(i,j)处的方向θ(i,j):
若h(i,j)>>v(i,j),即θ(i,j)=0;
2c)统计分割子图的θ值,找出出现次数最多的两个方向θ1和θ2;
2d)由θ1和θ2得到近似有理斜率r1和r2,r1≈arctanθ1=b1/a1,r2≈arctanθ2=b2/a2,根据有理斜率r1和r2,构造变换矩阵其中沿r1的方向称为变换矩阵MΛ的变换方向,沿r2的方向称为队列方向,a1,a2,b1,b2都是整数;
(3)利用变换矩阵MΛ,对各分割子图进行采样,得到分割子图的|det(MΛ)|个陪集,|det(MΛ)|是矩阵MΛ行列式的绝对值;
(4)对各分割子图的每个陪集沿Directionlet变换矩阵MΛ的变换方向和队列方向分别进行n1=2与n2=1次的一维小波变换,得到Directionlet变换的高频和低频子带系数;
(5)对各个高频子带,利用该子带所有的变换系数,估计广义高斯模型的形状参数υ和局部标准差σx;
(6)对各高频子带的广义高斯模型形状参数υ进行判断:
若0<υ<0.5,按下式对噪声图像的无噪系数进行估计,
其中 是无噪系数的估计,y是含噪系数,Tυ=Cυσ2/(2-υ)σx -(υ/(2-υ)),Cυ=(2-υ)(2-2υ)-(1-υ/2-υ)η(υ)υ/(2-υ), Γ是Gamma函数, o(y2(υ-1))是y2(υ-1)的高阶无穷小量,σ是噪声标准差;
(7)对低频子带和估计得到的无噪高频子带沿矩阵MΛ的变换方向和队列方向分别进行n1=2与n2=1次的一维小波逆变换;
(8)根据变换矩阵MΛ的变换方向和队列方向加权综合,重构各个分割子图;
(9)将重构的各分割子图按其在原图像中的位置合成,得到去噪后的图像。
2.根据权利要求1所述的去噪方法,其中步骤(5)所说的估计广义高斯模型的形状参数υ和局部标准差σx,按如下步骤进行:
5a)设高频子带的大小为l1×l2,由含噪的高频子带系数Y计算出含噪高频子带系数的二阶矩σY 2和含噪高频子带系数的峰度kY:
σY 4=(σY 2)2
5c)由以上得到的ky,σY 2和σ值,利用数值计算法求出形状参数υ:
5e)由形状参数υ和无噪系数的初始估计值X,计算出局部标准差:σx=SE
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2009102193489A CN101719268B (zh) | 2009-12-04 | 2009-12-04 | 基于改进Directionlet域的广义高斯模型图像去噪方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2009102193489A CN101719268B (zh) | 2009-12-04 | 2009-12-04 | 基于改进Directionlet域的广义高斯模型图像去噪方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101719268A CN101719268A (zh) | 2010-06-02 |
CN101719268B true CN101719268B (zh) | 2011-10-19 |
Family
ID=42433838
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2009102193489A Active CN101719268B (zh) | 2009-12-04 | 2009-12-04 | 基于改进Directionlet域的广义高斯模型图像去噪方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101719268B (zh) |
Families Citing this family (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101984461A (zh) * | 2010-11-05 | 2011-03-09 | 西安电子科技大学 | 基于可操纵金字塔的统计模型图像去噪方法 |
CN102073999B (zh) * | 2011-01-20 | 2012-08-29 | 西安电子科技大学 | 基于双冗余字典学习的自然图像去噪方法 |
CN102142133B (zh) * | 2011-04-19 | 2013-01-23 | 西安电子科技大学 | 基于非下采样Directionlet变换和压缩感知的乳腺X线图像增强方法 |
CN102903099A (zh) * | 2012-09-05 | 2013-01-30 | 西安电子科技大学 | 基于directionlet变换的彩色图像边缘检测方法 |
CN102867293A (zh) * | 2012-09-05 | 2013-01-09 | 西安电子科技大学 | 基于广义高斯最大似然估计的表面波变换视频去噪方法 |
CN103067671B (zh) * | 2012-12-31 | 2015-09-23 | 华为技术有限公司 | 一种显示图像的方法及装置 |
CN110097132B (zh) * | 2019-05-07 | 2020-12-08 | 电子科技大学 | 一种识别数字照片与拍摄相机的方法 |
CN110458773B (zh) * | 2019-08-01 | 2023-04-21 | 南京信息工程大学 | 一种基于边缘增强算子的各向异性扩散噪声处理方法 |
CN113436093A (zh) * | 2021-06-16 | 2021-09-24 | 中国电子科技集团公司第五十四研究所 | 一种基于能量域变换的遥感影像去噪方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1917577A (zh) * | 2006-09-01 | 2007-02-21 | 上海大学 | 一种图像组合降噪方法 |
CN101477679A (zh) * | 2009-01-16 | 2009-07-08 | 西安电子科技大学 | 基于轮廓波Contourlet变换的图像去噪方法 |
-
2009
- 2009-12-04 CN CN2009102193489A patent/CN101719268B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1917577A (zh) * | 2006-09-01 | 2007-02-21 | 上海大学 | 一种图像组合降噪方法 |
CN101477679A (zh) * | 2009-01-16 | 2009-07-08 | 西安电子科技大学 | 基于轮廓波Contourlet变换的图像去噪方法 |
Non-Patent Citations (4)
Title |
---|
S.Grace Chang,Bin Yu,Martin Vetterli.Adaptive Wavelet Thresholding for Image Denoising and Compression.《IEEE TRANSACTIONS ON IMAGE PROCESSING》.2000,第9卷(第9期),P1532-P1546. * |
Wentao Huang,Licheng Jiao,Jianhua Jia,Hang Yu.A neural contextual model for detecting perceptually salient contours.《Pattern Recognition Letters》.2009,P985-P993. * |
侯建华,熊承义,田晓梅.广义高斯分布及其在图像去噪中的应用.《中南民族大学学报自然科学版》.2005,第24卷(第3期),P45-P47. * |
白静,侯彪,王爽,焦立成.基于提升Directionlet域高斯混合尺度模型的SAR图像噪声抑制.《计算机学报》.2008,第31卷(第7期),P1234-P1241. * |
Also Published As
Publication number | Publication date |
---|---|
CN101719268A (zh) | 2010-06-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101719268B (zh) | 基于改进Directionlet域的广义高斯模型图像去噪方法 | |
CN101482617B (zh) | 基于非下采样轮廓波的合成孔径雷达图像去噪方法 | |
CN103854264B (zh) | 一种基于改进型阈值函数的小波变换图像去噪方法 | |
CN103761719B (zh) | 一种基于邻域相关性的自适应小波阈值去噪方法 | |
CN101739671B (zh) | 一种含噪图像的降噪处理方法及系统 | |
CN102509263B (zh) | 基于sar图像局部统计特性的k-svd相干斑抑制方法 | |
CN101477680A (zh) | 基于滑窗邻域数据选择的小波图像降噪方法 | |
CN101980286B (zh) | 结合双树复小波与双变量模型的sar图像降斑方法 | |
CN101477679A (zh) | 基于轮廓波Contourlet变换的图像去噪方法 | |
CN101719267B (zh) | 一种带噪图像降噪处理方法和系统 | |
CN101639537A (zh) | 基于方向波域混合高斯模型的sar图像噪声抑制方法 | |
CN101566688B (zh) | 基于邻域方向性信息的sar图像降斑方法 | |
CN101477681A (zh) | 基于自适应滑窗邻域阈值的小波图像降噪方法 | |
CN107451981A (zh) | 基于dct和梯度协方差矩阵的图像噪声水平估计方法 | |
CN107133938A (zh) | 基于小波系数统计模型的鲁棒图像融合方法 | |
CN103679663A (zh) | Tetrolet变换域结合PDE和GCV理论的图像去噪方法 | |
CN106023103B (zh) | 一种基于精确局部方差先验建模的自适应正交小波图像去噪方法 | |
CN102547073B (zh) | 基于表面波变换的自适应阈值视频去噪方法 | |
CN103955893B (zh) | 基于可分离全变差模型的图像去噪方法 | |
CN103426145A (zh) | 一种基于多分辨分析的合成孔径声纳相干斑噪声抑制方法 | |
CN103024377B (zh) | 基于目标区域提取与方向波的sar图像压缩方法 | |
CN104809714A (zh) | 基于多形态稀疏表示的图像融合方法 | |
CN101984461A (zh) | 基于可操纵金字塔的统计模型图像去噪方法 | |
CN104182944A (zh) | 基于曲波与小波变换相串联的光学图像去噪方法 | |
CN102509268B (zh) | 基于免疫克隆选择的非下采样轮廓波域图像去噪方法 |
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 |