CN111242856B - 一种基于剪切波的双重图像去噪方法 - Google Patents
一种基于剪切波的双重图像去噪方法 Download PDFInfo
- Publication number
- CN111242856B CN111242856B CN202010010833.1A CN202010010833A CN111242856B CN 111242856 B CN111242856 B CN 111242856B CN 202010010833 A CN202010010833 A CN 202010010833A CN 111242856 B CN111242856 B CN 111242856B
- Authority
- CN
- China
- Prior art keywords
- image
- frequency sequence
- denoising
- low
- frequency
- 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 46
- 230000009977 dual effect Effects 0.000 title claims abstract description 14
- 238000001914 filtration Methods 0.000 claims abstract description 30
- 230000009466 transformation Effects 0.000 claims abstract description 17
- 238000001035 drying Methods 0.000 claims abstract description 5
- 230000003044 adaptive effect Effects 0.000 claims description 6
- 238000000354 decomposition reaction Methods 0.000 abstract description 7
- 238000013519 translation Methods 0.000 abstract description 7
- 230000007547 defect Effects 0.000 abstract description 4
- 230000000694 effects Effects 0.000 description 13
- 238000010586 diagram Methods 0.000 description 5
- 238000012545 processing Methods 0.000 description 4
- 238000001514 detection method Methods 0.000 description 2
- 238000009499 grossing Methods 0.000 description 2
- 239000011159 matrix material Substances 0.000 description 2
- 239000006002 Pepper Substances 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 230000001105 regulatory effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000002441 reversible effect Effects 0.000 description 1
- 238000012163 sequencing technique Methods 0.000 description 1
- 238000010008 shearing Methods 0.000 description 1
- 230000001629 suppression Effects 0.000 description 1
Images
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
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10004—Still image; Photographic image
-
- 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/20024—Filtering details
-
- 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/20024—Filtering details
- G06T2207/20032—Median 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
-
- 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/20172—Image enhancement details
- G06T2207/20192—Edge enhancement; Edge preservation
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Image Processing (AREA)
Abstract
本发明公开了一种基于剪切波的双重图像去噪方法,包括:对噪声图像f(x,y)进行分解,得到不同尺度的低频序列和高频序列;分别对低频序列、高频序列进行去燥处理;将低频序列、高频序列进行逆变换,得到一重去噪后的图像F1;对图像F1进行二重去噪,得到图像F2。采用剪切波多尺度变换分解,克服小波分解不够精细的缺点,能够更好的对弱小信号进行区分,使该变换具有平移不变性和稳定性等优点;对三个高频序列采用自适应阈值法进行去噪处理,能够更好的消除图像中大部分的噪声,保持图像的边缘;对低频序列系数采用中值滤波去噪,能够消除图像中的孤立噪声,更好的维持图像的清晰度。
Description
技术领域
本发明属于图像处理方法技术领域,涉及一种基于剪切波的双重图像去噪方法。
背景技术
随着科技的不断进步,图像中包含的信息也越来越被人重视。图像处理是实现图像信息提取的重要前提,而在图像信息被提取或传输时容易受到不同类型噪声干扰,图像中的噪声会使图像质量大大降低,甚至模糊了图像的特征,掩盖了重要信息,造成图像的失真,因此对图像进行高效的去噪预处理异常重要。在众多学者的研究下,为了能够更好的减少或消除图像中的噪声,产生了许多的图像去噪方法,如小波阈值法、维纳滤波方法、中值滤波方法、高斯滤波方法等。这些方法各自去噪效果却各不相同。
小波阈值法因操作简单,受到了大家的青睐,此方法首先是把图像分为低频序列和高频序列,利用阈值化对高频序列的噪声进行处理,再对修正后的小波系数进行拼接,得到去噪后的图像。小波阈值法去噪效果的关键环节就是阈值的估量和阈值函数的选择,阈值选择过小,图像中的噪声会存有残留;选择太大则会覆盖掉图像中的细节特征,导致图像模糊。所以小波阈值函数去噪声的局限性在于阈值判定准则的选取。其他的一些方法如中值滤波对椒盐噪声去噪效果较好,空域低通滤波对高斯噪声去噪效果较好,但是这几种方法在去噪的同时模糊甚至是破坏了图像的边缘细节,影响了图像整体的去噪效果。
发明内容
本发明的目的是提供一种基于剪切波的双重图像去噪方法,解决了现有技术中存在的去噪时模糊甚至破坏图像边缘细节的问题。
本发明所采用的技术方案是,一种基于剪切波的双重图像去噪方法,包括:
步骤1、对噪声图像f(x,y)进行分解,得到不同尺度的低频序列和高频序列;
步骤2、分别对低频序列、高频序列进行去燥处理;
步骤3、将经步骤2处理后的低频序列、高频序列进行逆变换,得到一重去噪后的图像F1;
步骤4、对图像F1进行二重去噪,得到图像F2。
本发明的特点还在于,
步骤1具体过程为:采用剪切波变换对噪声图像f(x,y)进行分解,得到水平和竖直方向低频序列LL、水平低频和竖直高频序列LH、水平高频和竖直高频序列HL、水平高频和竖直高频序列HH。
步骤2中:
步骤4中采用维纳滤波对图像F1进行二重去噪处理。
本发明的有益效果是:
本发明的基于剪切波的双重图像去噪方法,采用剪切波多尺度变换分解,克服小波分解不够精细的缺点,能够更好的对弱小信号进行区分,使该变换具有平移不变性和稳定性等优点;特别是对于弱小目标检测,不会因为下采样操作而导致目标信息丢失;对三个高频序列采用自适应阈值法进行去噪处理,能够更好的消除图像中大部分的噪声,保持图像的边缘;对低频序列系数采用中值滤波去噪,能够消除图像中的孤立噪声,更好的维持图像的清晰度;对重构图像F1采用维纳滤波进行二重去噪处理,能更好的处理图像中的高斯噪声和乘性噪声,清除图像中的残留噪声,使图像变得更平滑,能保持图像的纹理结构等信息。
附图说明
图1是本发明一种基于剪切波的双重图像去噪方法的流程图;
图2是本发明一种基于剪切波的双重图像去噪方法的剪切波变换分解示意图;
图3是本发明一种基于剪切波的双重图像去噪方法的维纳滤波原理图;
图4是小波阈值去噪方法去除图像噪声的效果图;
图5是维纳滤波去噪方法去除图像噪声的效果图;
图6是本发明一种基于剪切波的双重图像去噪方法去除图像噪声的效果图。
具体实施方式
下面结合附图和具体实施方式对本发明进行详细说明。
一种基于剪切波的双重图像去噪方法,如图1所示,包括以下步骤:
步骤1、对噪声图像f(x,y)进行分解,得到不同尺度的低频序列和高频序列;
剪切波变换是对小波变换的继承和发展,它克服了小波变换的方向局限性和平移不变性的缺陷,其实现过程更加灵活有效。剪切波变换结合了小波的多尺度特点,能根据方向进行剖分,它构造简单,通过对函数进行伸缩、平移、旋转操作生成基函数,表达式为:
进一步的,如图2所示,采用剪切波变换对噪声图像f(x,y)进行分解,得到水平和竖直方向低频序列LL(体现图像的主体轮廓)、水平低频和竖直高频序列LH、水平高频和竖直高频序列HL、水平高频和竖直高频序列HH,三个高频序列反映图像的细节边缘,包含着大量的噪声。
步骤2、分别对低频序列、高频序列进行去噪处理;
具体原理为:在低频序列中选择奇数像素的窗口,将窗口内的像素按照灰度值排序,以中间位置像素的灰度值代替图像f(x,y)像素点的灰度值,作为窗口的灰度值g(x,y):
g(x,y)=median{f(x-k,y-l),(k,l∈M} (3);
中值滤波的输出像素是由邻域像素的中间值决定,而且中值滤波产生的模数较少,更适合消除低频序列的孤立噪声点,同时又能保持图像的细节。
自适应阈值法去噪的思想不是计算全局图像的阈值,而是根据图像不同区域亮度分布,计算其局部阈值,所以对于图像不同区域,能够自适应计算不同的阈值,因此被称为自适应阈值法。采用自适应阈值T分别对高频序列LH、HL、HH进行去噪处理:
其中,自适应阈值T的表达式为:
σx服从高斯分布的前提下有:
步骤3、将经步骤2处理后的低频序列、高频序列进行逆变换,得到一重去噪后的图像F1,重构图像的对比度得到了很大改善,大部分噪声也能够被滤除。
步骤4、采用维纳滤波对图像F1进行二重去噪,得到图像F2;
维纳滤波通常用于提取被噪声污染的图像,它是以最小均方误差准则进行滤波的,维纳滤波的主要思想分为三步:滤波、预测和平滑;滤波是利用直到当前时刻的随机过程的观测值,来得到当前信号的估计值;预测是利用直到当前时刻的随机过程的观测值,得到将来某个时刻信号的估计值;平滑是利用直到当前时刻的随机过程的观测值,得到过去某个时刻的估计值。具体原理如下:
如图3所示,将图像F1作为输入图像x(n),则x(n)=s(n)+v(n),s(n)表示无噪图像,v(n)表示噪声;对输入图像x(n)进行维纳滤波后输出y(n),即为图像F2,输出y(n)无限接近s(n),为期望值,则误差为:
均方误差为:
J=E[e2(n)] (10);
则均方误差J最小时滤波器系数为最优解。
分别对小波阈值去噪法、维纳滤波方法及本发明的去燥方法进行试验,去噪效果如下:
采用小波阈值去噪法虽然对噪声有一定的抑制效果,但会使图像变得模糊,低频成分通过,高频边缘细节则被滤去,见图4;维纳滤波方法保留图像结构信息,但去噪效果不太理想,存在残留噪声,见图5;本发明的去燥方法不仅能很好地抑制相关噪声,而且最大限度地保护去噪后图像的边缘信息,图像也较为清晰,取得更好的去噪效果,见图6。
表1为上述三种方法的去噪结果评价指标,峰值信噪比值与去噪效果成正比,值越大表明去噪效果越好,从表1可看出本发明去燥方法的峰值信噪比最高。
通过以上方式,本发明的基于剪切波的双重图像去噪方法,采用剪切波多尺度变换分解,克服小波分解不够精细的缺点,能够更好的对弱小信号进行区分,使该变换具有平移不变性和稳定性等优点;特别是对于弱小目标检测,不会因为下采样操作而导致目标信息丢失;对三个高频序列采用自适应阈值法进行去噪处理,能够更好的消除图像中大部分的噪声,保持图像的边缘;对低频序列系数采用中值滤波去噪,能够消除图像中的孤立噪声,更好的维持图像的清晰度;对重构图像F1采用维纳滤波进行二重去噪处理,能更好的处理图像中的高斯噪声和乘性噪声,清除图像中的残留噪声,使图像变得更平滑,能保持图像的纹理结构等信息。
Claims (2)
1.一种基于剪切波的双重图像去噪方法,其特征在于,包括:
步骤1、对噪声图像f(x,y)进行分解,得到不同尺度的低频序列和高频序列;
步骤2、分别对所述低频序列、高频序列进行去燥处理;
步骤3、将经步骤2处理后的低频序列、高频序列进行逆变换,得到一重去噪后的图像F1;
步骤4、对所述图像F1进行二重去噪,得到图像F2;
所述步骤1具体过程为:采用剪切波变换对噪声图像f(x,y)进行分解,得到水平和竖直方向低频序列LL、水平低频和竖直高频序列LH、水平高频和竖直高频序列HL、水平高频和竖直高频序列HH;
所述步骤2中:
2.根据权利要求1所述的一种基于剪切波的双重图像去噪方法,其特征在于,步骤4中采用维纳滤波对所述图像F1进行二重去噪处理。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010010833.1A CN111242856B (zh) | 2020-01-06 | 2020-01-06 | 一种基于剪切波的双重图像去噪方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010010833.1A CN111242856B (zh) | 2020-01-06 | 2020-01-06 | 一种基于剪切波的双重图像去噪方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111242856A CN111242856A (zh) | 2020-06-05 |
CN111242856B true CN111242856B (zh) | 2023-03-24 |
Family
ID=70870724
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010010833.1A Active CN111242856B (zh) | 2020-01-06 | 2020-01-06 | 一种基于剪切波的双重图像去噪方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111242856B (zh) |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102890820A (zh) * | 2012-09-18 | 2013-01-23 | 西安电子科技大学 | 基于shearlet变换和维纳滤波的图像去噪方法 |
EP3051485A1 (en) * | 2013-11-08 | 2016-08-03 | Huawei Device Co., Ltd. | Method and terminal used for image noise reduction |
CN106408532A (zh) * | 2016-09-09 | 2017-02-15 | 西北工业大学 | 基于剪切波域参数估计的合成孔径雷达sar图像去噪方法 |
CN109472756A (zh) * | 2018-11-15 | 2019-03-15 | 昆明理工大学 | 基于剪切波变换和带方向性局部维纳滤波的图像去噪方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP6207972B2 (ja) * | 2013-10-31 | 2017-10-04 | 東芝メディカルシステムズ株式会社 | 超音波診断装置、画像処理装置及び画像処理プログラム |
-
2020
- 2020-01-06 CN CN202010010833.1A patent/CN111242856B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102890820A (zh) * | 2012-09-18 | 2013-01-23 | 西安电子科技大学 | 基于shearlet变换和维纳滤波的图像去噪方法 |
EP3051485A1 (en) * | 2013-11-08 | 2016-08-03 | Huawei Device Co., Ltd. | Method and terminal used for image noise reduction |
CN106408532A (zh) * | 2016-09-09 | 2017-02-15 | 西北工业大学 | 基于剪切波域参数估计的合成孔径雷达sar图像去噪方法 |
CN109472756A (zh) * | 2018-11-15 | 2019-03-15 | 昆明理工大学 | 基于剪切波变换和带方向性局部维纳滤波的图像去噪方法 |
Non-Patent Citations (2)
Title |
---|
基于小波变换和均值滤波的图像去噪方法;梁利利等;《计算机与数字工程》(第05期);全文 * |
结合NSST和快速非局部均值滤波的刀具图像去噪;龙云淋等;《信号处理》(第11期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN111242856A (zh) | 2020-06-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112862716A (zh) | 基于双树复数小波和模糊理论的声呐图像目标增强方法 | |
Kumar et al. | Image Denoising by Wavelet Based Thresholding Method | |
Biswas et al. | An image denoising threshold estimation method | |
CN104616259A (zh) | 一种噪声强度自适应的非局部均值图像去噪方法 | |
CN111242856B (zh) | 一种基于剪切波的双重图像去噪方法 | |
Liu et al. | Wavelet image threshold denoising based on edge detection | |
Singh et al. | Noise reduction in ultrasound images using wavelet and spatial filtering techniques | |
Kaur et al. | Image de-noising techniques: a review paper | |
kumar Dass et al. | Improvising MSN and PSNR for finger-print image noised by GAUSSIAN and SALT & PEPPER | |
Bhonsle et al. | Medical Image De-Noising Using Combined Bayes Shrink and Total Variation Techniques | |
CN111147804A (zh) | 一种视频帧重建方法 | |
Liu | New method for image denoising while keeping edge information | |
Charde | A review on image denoising using wavelet transform and Median filter over AWGN channel | |
CN112465719A (zh) | 变换域图像去噪方法及系统 | |
Wu et al. | Research on medical image denoising based on wavelet pyramid fusion | |
Zhao et al. | An improved Roberts edge detection algorithm based on mean filter and wavelet denoising | |
Liu et al. | Bilateral Filtering for Image Processing Based on Pulse Coupled Neural Networks | |
Zhang et al. | Translation invariant RI-Spline wavelet and its application on de-noising | |
Cristobal et al. | Wavelet-based denoising methods: A comparative study with applications in microscopy | |
Santhanamari et al. | Adaptive bilateral filter and Bayesian threshold based image denoising | |
Khan et al. | A novel algorithm for removal of noise from X-Ray images | |
Dash et al. | Restoration of images corrupted with blur and impulse noise | |
Chen et al. | All Characteristics Preservation: Single Image Dehazing based on Hierarchical Detail Reconstruction Wavelet Decomposition Network | |
CN114820330A (zh) | 一种增强型显微图像降噪恢复方法 | |
JEBUR | Image Denoising using Mean Filter |
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 |