CN107330863B - 一种基于噪声估计的图像去噪方法 - Google Patents
一种基于噪声估计的图像去噪方法 Download PDFInfo
- Publication number
- CN107330863B CN107330863B CN201710390337.1A CN201710390337A CN107330863B CN 107330863 B CN107330863 B CN 107330863B CN 201710390337 A CN201710390337 A CN 201710390337A CN 107330863 B CN107330863 B CN 107330863B
- Authority
- CN
- China
- Prior art keywords
- pixel
- image
- super
- denoising
- smooth
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 25
- 230000011218 segmentation Effects 0.000 claims abstract description 10
- 238000009499 grossing Methods 0.000 claims description 4
- 229910010888 LiIn Inorganic materials 0.000 claims description 2
- 230000000694 effects Effects 0.000 abstract description 6
- 239000000463 material Substances 0.000 abstract description 5
- 238000010586 diagram Methods 0.000 description 3
- 230000003044 adaptive effect Effects 0.000 description 1
- 239000000654 additive Substances 0.000 description 1
- 230000000996 additive effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 239000003550 marker Substances 0.000 description 1
- 238000012634 optical imaging Methods 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 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
- 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
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20021—Dividing image into blocks, subimages or windows
-
- 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/20182—Noise reduction or smoothing in the temporal domain; Spatio-temporal 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/20172—Image enhancement details
- G06T2207/20192—Edge enhancement; Edge preservation
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
本发明公开了一种基于噪声估计的图像去噪方法,利用超像素分割根据图像内容将其分为若干同源区域,利用平坦信息更能代表图像噪声污染水平的先验知识,以图像信息熵为依据找出图像中较为平滑的同源区域,并估计平滑区域的噪声标准差,以这些平滑区域的噪声标准差作为整幅图像的噪声水平,达到更加准确估计噪声水平的目的,从而根据噪声水平对非局部均值(NLM)去噪方法进行修正,实现根据噪声程度合理地控制去噪程度,能够自适应地处理带噪图像,使去噪后的图像整体效果较传统的NLM算法大幅度提升,在去噪的同时更好地保留了细节。全过程可实现自动化与自能化,无需人工干预。
Description
技术领域
本发明涉及图像处理领域,尤其涉及一种噪声估计的图像去噪方法。
背景技术
随着数字图像和多媒体技术的快速发展,各种类型的光学成像***也越来越多,人们对成像***的画面质量要求也越来越高,但在图像的传输过程中,难免会受到各种噪声的污染,从而导致所获取的图像与原始的图像质量有所降低,不仅会影像图像的效果,噪声严重时甚至会影响图像的可读性,从而导致后续获取信息出现错误。因此,在图像处理领域中,对数字图像进行去噪是十分有意义的。
目前一种常用的去噪方法是性能优异且边缘保持能力较强的非局部均值(NLM)图像去噪算法。其基本思想是在整幅图像中寻找与待去噪像素相似的像素并以这些相似像素的权值平均值作为其去噪结果,相似度越高的像素对应的权值也越高。其优点在于引入非局部思想和基于图像块的相似度的概念。在自然图像中,每个像素对应的相似像素往往不仅限于其周围较小范围内,且去噪效果也一般与参与去噪的相似像素数量成正比,因此引入非局部搜索方式对提升去噪性能具有重大意义。相似度的准确性是非局部搜索发挥作用的基本保证,若相似度的准确性无法得到保证,使用图像块向量间的高斯加权欧式距作为像素间的相似度,显著提升了噪声环境中相似度的鲁棒性。
带噪图像中的噪声通常被认为是均值为零标准差位置的加性白噪声,因此,合理地估计噪声水平对于去噪可以起到一定的引导作用,而图像中平滑区域的噪声水平往往比细节丰富的区域更接近图像的噪声水平,因此利用超像素分割得到的图像内容相关的同源块,避免了边缘的影响,同时,结合信息熵来评价超像素同源区域的平滑程度,信息熵值越小代表图像中的纹理信息越少,则图像越平滑。因此,选出较为平滑的超像素同源区域用来估计图像的噪声水平,以此为依据修正NLM去噪算法,更有效地对图像进行去噪,在对图像进行去噪的同时很好地保留了细节和边缘信息。
发明内容
本发明的目的在于提出一种自适应的,易于实现,鲁棒性强的图像去噪方法,从而解决现有图像去噪方法由于未知噪声水平而引起的去噪程度不足或去噪过度的问题,并在去噪的同时有效地保留了边缘和纹理信息。
本发明的目的是通过以下技术方案实现的:一种图像去噪的方法,该方法包括以下步骤:
(1)对带噪图J进行基于熵率的超像素分割,得到若干超像素同源区域Li;
(2)根据图像信息熵理论,计算每个超像素同源区域的信息熵:
Li代表第i个超像素同源区域,l代表像素的灰度级,代表第i个超像素同源区域中灰度级为l的概率,代表Li中的信息熵;
(3)根据计算得到的各个超像素同源区域的信息熵的大小进行排序,选出值最小的t个超像素同源区域(占总数的10~20%),即为图像中最为平滑的部分区域,并标记为1,其他区域标记为0,得到平滑区域的二值图;
(4)计算各个平滑同源区域的标准差
是Li中像素x的强度,代表Li中平均强度,为超像素块中像素的个数;
(5)计算上述平滑同源区域的平均标准差,以此作为估计得到的整幅图像的噪声水平:
(6)利用估计得到的噪声水平修正NLM去噪算法对带噪图像进行去噪,得到去噪图F,其中,去噪图像F中的像素x的强度F(x)为:
I代表带噪图J中像素x的邻域,y表示邻域I中的像素,J(y)表示像素y在带噪图J中的强度;w(x,y)代表像素x与像素y之间的相似性:
Z(x)是像素x处的归一化参数,h代表控制去噪程度的平滑参数,
C为常数,C∈[300,500] (16)。
本发明的有益效果在于:利用超像素分割得到的图像内容相关的同源块,避免了边缘的影响,同时,结合信息熵来评价超像素同源区域的平滑程度,信息熵值越小代表图像中的纹理信息越少,则图像越平滑。因此,选出较为平滑的超像素同源区域用来估计图像的噪声水平,以此为依据修正NLM去噪算法,更有效地对图像进行去噪,在对图像进行去噪的同时很好地保留了细节和边缘信息。
附图说明
图1为发明方法的流程框图。
图2为任意带噪图像。
图3为带噪图像进行超像素分割得到的分割图。
图4为通过计算同源区域的信息熵得到的二值图。
图5为图2的最终去噪效果图。
图6~9为带噪图、使用NLM方法去噪图、以及本发明方法(S+NLM)去噪图的三个比较示例图。
具体实施方式
本发明一种基于噪声估计的图像去噪方法,先利用超像素分割根据图像内容将其分为若干同源区域,以图像信息熵为依据找出图像中较为平滑的同源区域,并估计平滑区域的噪声标准差,从而根据噪声水平对非局部均值(NLM)去噪方法进行修正,实现根据噪声程度合理地控制去噪程度,在去噪的同时更好地保留了细节。
下面结合附图和实例进行详细说明:
图1为本发明方法的简易流程框图。下面结合实施例对本发明作进一步说明。
(1)获取一张带噪图J,如图2所示;
(2)对图像进行基于熵率的超像素分割,得到若干超像素同源区域,本实施例中,共分割出200个同源区域,如图3所示;
(3)根据图像信息熵理论,计算每个超像素同源块的信息熵:
Li代表第i个超像素同源区域,l代表像素的灰度级,代表第i个超像素同源区域Li中灰度级为l的概率,代表Li的信息熵;
(4)根据计算得到的各个超像素同源区域的信息熵的大小进行排序,选出值最小的20个超像素同源区域,即为图像中最为平滑的部分区域,本实施例中以选出百分之十为例,并将这部分同源区域标记为1,其他区域标记为0,得到平滑区域的二值图,如图4所示;
(5)计算各个平滑同源区域的标准差:
是Li超像素块中像素x的强度,代表Li超像素块中平均强度,为超像素块中像素的个数;
(6)计算上述平滑同源区域的平均标准差,以此作为估计得到的整幅图像的噪声水平:
t代表所选出的超像素同源区域的个数,本实施例中为20个;
(7)利用估计得到的噪声水平修正NLM去噪算法对带噪图像进行去噪,得到图5所示的最终去噪图。计算像素x在去噪图像F中的强度:
I代表带噪图J中像素x的邻域,y表示邻域I中的像素,J(y)表示像素y在带噪图J中的强度;w(x,y)代表像素x与像素y之间的相似性,其值决定于以x与y为中心的矩形区域的距离,计算公式如下:
Z(x)是归一化参数,h代表控制去噪程度地平滑参数项,这里根据所估计的噪声水平来修正平滑参数项,
C为常数,本实施例中取值为300,通过估计得到的噪声水平来控制去噪程度,以得到更好的去噪效果。
本发明一种基于噪声估计的去噪方法,先利用超像素分割根据图像内容将其分为若干同源区域,以图像信息熵为依据找出图像中较为平滑的同源区域,并估计较为平滑区域的噪声标准差,以这些平滑区域的噪声标准差作为整幅图像的噪声水平,从而根据噪声水平对非局部均值(NLM)去噪方法进行修正,实现根据噪声程度合理地控制去噪程度,能够自适应地处理带噪图像。分别按照现有的NLM方法进行去噪以及本发明的去噪方法对如图6~9所示的4幅带噪图像进行去噪处理,从图中可以看出,去噪后的图像整体效果较传统的NLM算法大幅度提升,在去噪的同时更好地保留了细节。
Claims (1)
1.一种基于噪声估计的图像去噪方法,其特征在于,该方法包括以下步骤:
(1)对带噪图J进行基于熵率的超像素分割,得到若干超像素同源区域Li;
(2)根据图像信息熵理论,计算每个超像素同源区域的信息熵:
Li代表第i个超像素同源区域,l代表像素的灰度级,代表第i个超像素同源区域中灰度级为l的概率,代表Li中的信息熵;
(3)根据计算得到的各个超像素同源区域的信息熵的大小进行排序,选出值最小的t个超像素同源区域占总数的10~20%,即为图像中最为平滑的部分区域,并标记为1,其他区域标记为0,得到平滑区域的二值图;
(4)计算各个平滑同源区域的标准差
是Li中像素x的强度,代表Li中平均强度,为超像素块中像素的个数;
(5)计算上述平滑同源区域的平均标准差,以此作为估计得到的整幅图像的噪声水平:
(6)利用估计得到的噪声水平修正NLM去噪算法对带噪图像进行去噪,得到去噪图F,其中,去噪图像F中的像素x的强度F(x)为:
I代表带噪图J中像素x的邻域,y表示邻域I中的像素,J(y)表示像素y在带噪图J中的强度;w(x,y)代表像素x与像素y之间的相似性:
Z(x)是像素x处的归一化参数,h代表控制去噪程度的平滑参数,
C为常数,C∈[300,500] (8)。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710390337.1A CN107330863B (zh) | 2017-05-27 | 2017-05-27 | 一种基于噪声估计的图像去噪方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710390337.1A CN107330863B (zh) | 2017-05-27 | 2017-05-27 | 一种基于噪声估计的图像去噪方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107330863A CN107330863A (zh) | 2017-11-07 |
CN107330863B true CN107330863B (zh) | 2019-10-18 |
Family
ID=60192961
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710390337.1A Expired - Fee Related CN107330863B (zh) | 2017-05-27 | 2017-05-27 | 一种基于噪声估计的图像去噪方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107330863B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109064418B (zh) * | 2018-07-11 | 2022-03-08 | 成都信息工程大学 | 一种基于非局部均值的非均匀噪声图像去噪方法 |
CN109544468B (zh) * | 2018-10-26 | 2023-10-13 | 浙江师范大学 | 一种图像数据扩增方法 |
WO2022246663A1 (zh) * | 2021-05-25 | 2022-12-01 | 深圳市大疆创新科技有限公司 | 图像处理方法、设备、***和存储介质 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101661611A (zh) * | 2009-09-25 | 2010-03-03 | 西安电子科技大学 | 基于贝叶斯非局部均值滤波器的实现方法 |
CN102298773B (zh) * | 2011-09-19 | 2012-12-26 | 西安电子科技大学 | 一种形状自适应的非局部均值去噪方法 |
CN103116873A (zh) * | 2013-03-20 | 2013-05-22 | 中国科学院深圳先进技术研究院 | 图像降噪方法 |
CN103208105A (zh) * | 2013-05-02 | 2013-07-17 | 中国电子科技集团公司第二十八研究所 | 一种红外图像细节增强与噪声自适应抑制方法 |
CN105913396A (zh) * | 2016-04-11 | 2016-08-31 | 湖南源信光电科技有限公司 | 一种噪声估计的图像边缘保持混合去噪方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9196017B2 (en) * | 2013-11-15 | 2015-11-24 | Linear Algebra Technologies Limited | Apparatus, systems, and methods for removing noise from an image |
-
2017
- 2017-05-27 CN CN201710390337.1A patent/CN107330863B/zh not_active Expired - Fee Related
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101661611A (zh) * | 2009-09-25 | 2010-03-03 | 西安电子科技大学 | 基于贝叶斯非局部均值滤波器的实现方法 |
CN102298773B (zh) * | 2011-09-19 | 2012-12-26 | 西安电子科技大学 | 一种形状自适应的非局部均值去噪方法 |
CN103116873A (zh) * | 2013-03-20 | 2013-05-22 | 中国科学院深圳先进技术研究院 | 图像降噪方法 |
CN103208105A (zh) * | 2013-05-02 | 2013-07-17 | 中国电子科技集团公司第二十八研究所 | 一种红外图像细节增强与噪声自适应抑制方法 |
CN105913396A (zh) * | 2016-04-11 | 2016-08-31 | 湖南源信光电科技有限公司 | 一种噪声估计的图像边缘保持混合去噪方法 |
Also Published As
Publication number | Publication date |
---|---|
CN107330863A (zh) | 2017-11-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Zheng et al. | Single-image vignetting correction | |
CN103679173B (zh) | 图像显著区域检测方法 | |
CN107330863B (zh) | 一种基于噪声估计的图像去噪方法 | |
CN109377450B (zh) | 一种边缘保护的去噪方法 | |
Wang et al. | Variational single nighttime image haze removal with a gray haze-line prior | |
CN110428371A (zh) | 基于超像素分割的图像去雾方法、***、存储介质及电子设备 | |
CN103914820B (zh) | 一种基于图像分层增强的图像去雾方法及*** | |
Hua et al. | Extended guided filtering for depth map upsampling | |
CN108133462B (zh) | 一种基于梯度场区域分割的单幅图像的复原方法 | |
CN112488955B (zh) | 一种基于波长补偿的水下图像复原方法 | |
CN108182671B (zh) | 一种基于天空区域识别的单幅图像去雾方法 | |
CN112200746B (zh) | 一种雾天交通场景图像的去雾方法和设备 | |
CN107403414B (zh) | 一种利于模糊核估计的图像区域选择方法和*** | |
CN102208101A (zh) | 一种红外图像的自适应线性变换增强方法 | |
CN102306307B (zh) | 一种彩色显微图像序列中固点噪声的定位方法 | |
CN109410235B (zh) | 融合边缘特征的目标跟踪方法 | |
CN108229247A (zh) | 一种移动车辆检测方法 | |
CN111476744B (zh) | 一种基于分类及大气成像模型的水下图像增强方法 | |
CN105023246B (zh) | 一种基于对比度和结构相似度的图像增强方法 | |
Liu et al. | Texture filtering based physically plausible image dehazing | |
CN110111342B (zh) | 一种抠图算法的优化选取方法及装置 | |
CN112053298A (zh) | 一种图像去雾方法 | |
CN111738931A (zh) | 光伏阵列无人机航拍图像的阴影去除算法 | |
CN117541582A (zh) | 一种高频变频器igbt绝缘质量检测方法 | |
CN117092647A (zh) | 一种制作区域星载光学和sar影像dom的方法和*** |
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: 20191018 Termination date: 20200527 |
|
CF01 | Termination of patent right due to non-payment of annual fee |