CN103593827A - 一种基于证据理论的图像去噪方法 - Google Patents
一种基于证据理论的图像去噪方法 Download PDFInfo
- Publication number
- CN103593827A CN103593827A CN201310559241.5A CN201310559241A CN103593827A CN 103593827 A CN103593827 A CN 103593827A CN 201310559241 A CN201310559241 A CN 201310559241A CN 103593827 A CN103593827 A CN 103593827A
- Authority
- CN
- China
- Prior art keywords
- row
- evidence
- gray scale
- synthetic
- current pixel
- 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
Images
Landscapes
- Image Processing (AREA)
Abstract
Description
技术领域
本发明涉及图像处理技术领域,具体涉及一种基于证据理论的图像去噪方法。
背景技术
随着各种数字仪器和数码产品的普及,图像和视频已成为人类活动中最常用的信息载体,它们包含着物体的大量信息,成为人们获取外界原始信息的主要途径.然而在图像的获取、传输和存贮过程中常常会受到各种噪声的干扰和影响而使图像降质,并且图像预处理算法的好坏又直接关系到后续图像处理的效果,如图像分割、目标识别、边缘提取等,所以为了获取高质量数字图像,很有必要对图像进行降噪处理,尽可能的保持原始信息完整性(即主要特征)的同时,又能够去除信号中无用的信息。所以,去噪处理一直是图像处理和计算机视觉研究的热点。
图像视频去噪的最终目的是改善给定的图像,解决实际图像由于噪声干扰而导致图像质量下降的问题.通过去噪技术可以有效地提高图像质量,增大信噪比,更好的体现原来图像所携带的信息,作为一种重要的预处理手段,人们对图像去噪算法进行了广泛的研究,在现有的去噪算法中,有的去噪算法在低维信号图像处理中取得较好的效果,却不适用于高维信号图像处理;或者去噪效果较好,却丢失部分图像边缘信息,或者致力于研究检测图像边缘信息,保留图像细节。如何在抵制噪音和保留细节上找到一个较好的平衡点,成为近年来研究的重点。
传统去噪方法中,均值滤波是典型的线性滤波,对高斯噪声具有较好的抑制效果;中值滤波是常用的非线性滤波方法,对椒盐噪声特别有效;维纳滤波对高斯噪声有明显的抑制作用;对小波系数进行阈值处理可以在小波变换域中去除低幅值的噪声和不期望的信号。均值滤波是一种典型的线性去噪方法,因为其运算简单快速, 同时又能够较为有效地去除高斯噪声,因而适用面较广,至今仍是一种常用的去噪方法, 许多滤除噪声方法都是在此基础上发展而来,其缺点是严重破坏了图像的边缘,模糊了图像。中值滤波去除脉冲噪声的效果很好,然而抑制噪声和保持图像中的细节往往是一对矛盾,也是图像处理中尚未很好解决的问题。所以保护细节的中值滤波的研究成为非线性滤波器研究的一个重要方面。小波变换对突变信号和非平稳信号处理具有优异性能,能较好地模拟视觉模型,使得它在数字图像处理、计算机视觉等方面引起了广泛的关注。但是,由于小波变换缺乏方向性, 仅具有水平、垂直、对角方向的信息, 不能很好地捕获二维图像中的线和面奇异,不能最优地表示含线或面奇异的二维图像,从而使得传统小波变换在处理二维图像时表现出一定的局限性。维纳滤波的优点是适应面较广,无论平稳随机过程是连续的还是离散的,是标量的还是向量的,都可应用,对某些问题,还可求出滤波器传递函数的显式解,并进而采用由简单的物理元件组成的网络构成维纳滤波器。 维纳滤波器的缺点是,要求得到半无限时间区间内的全部观察数据的条件很难满足,同时它也不能用于噪声为非平稳的随机过程的情况,对于向量情况应用也不方便。因此,维纳滤波在实际问题中应用不多。
发明内容
本发明所要解决的技术问题是提供一种基于证据理论的图像去噪方法。
本发明解决其技术问题所采用的技术方案:
一种基于证据理论的图像去噪方法,其步骤如下:
(1) 读入一帧图像,令第一个像素点为当前像素点;
(2) 判断所述当前像素点是否为噪声点,若是则按照以下步骤进行处理,若不是则置下一个像素点为当前像素点,接着判断,直到遇到噪声点或算法结束;上述判断噪声点的方法采用区间法,即灰度值位于0-5或250-255之间的为噪声点,否则为信号点;
(3) 以噪声点为中心,在所述图像上对目标像素构造一个n×n正方形模板,该模板为在横轴和纵轴上均包含n个像素点,n为除“1”之外的奇数;
如果图像边界的像素点为噪声,则对于以该像素点为中心构造的n×n正方形模板中超出图像区域的部分忽略不计,防止超出边界时发生溢出问题;
式中:
i为行标,i=-N, ……,-2,-1,0,1,2,……N;
j为列标,j=-N, ……,-2,-1,0,1,2,……N;
所述的i为n×n正方形模板中的行标,其排列顺序为从上至下依次排列,其中-N为第1行,N为第n行;
所述的j为n×n正方形模板中的列标,其排列顺序为从左至右依次排列,其中-N为第1列,N为第n列;
式中:
i为行标,i=-N, ……,-2,-1,0,1,2,……N;
j为列标,j=-N, ……,-2,-1,0,1,2,……N;
a. 第一轮灰度证据合成
b.第二轮灰度证据合成:
将上述第一轮灰度证据合成的新的各行的灰度证据再进行第二轮灰度证据合成,第二轮合成后,得到第二轮合成后的各行的灰度证据;
c.再进行第三轮合成及更多轮的灰度证据合成,直至只剩下两行时为止,将最后剩下的两行进行证据合成得到灰度证据合成值m 0 ;
所述灰度证据合成值m 0 按下述公式计算:
(7)噪音去除:
(8) 如果所述图像中所有的像素点处理完毕,则算法结束;否则,置下一个像素点为当前像素点,返回上述步骤(2),继续处理。
本发明采用D-S证据理论,其结构的最大特点是在证据中引入了不确定性,建立了基本概率赋值函数、信任函数、似真度函数等满足比概率论弱的公理,放松了传统理论需要完整的先验概率和条件概率知识以及需要有统一的辨识框架要求,可对相互重叠、非互不相容的命题进行证据组合等特点。
以下介绍本发明的理论基础:
①;
似真函数定义为
, (4)
函数称为下限函数,表示了对命题的总的信任程度,其中,.函数也称为上限函数或不否定函数,表示不否定的信任度,是所有与相交的集合的之和.当证据拒绝时,等于零,当没有证据反对时,为1,容易证明。就表示对的信任区间.如果辨识框架的一个子集为,且,则称的子集为信任函数的焦元。
;
本发明的有益效果为:本发明的去噪方法是根据证据理论计算的像素点的灰度证据,再利用证据合成公式进行灰度证据合成,能有效减少计算复杂度,得到更为精确的去噪效果。本发明更适合处理含不确定灰度证据的图像,在处理全局噪声判断方面优于传统的去噪方法。
附图说明
图1 为证据理论中的证据区间示意图;
图2 为3×3正方形模板及其各像素的灰度证据示意图;
图3 为5×5正方形模板及其各像素的灰度证据示意图;
图4-1至图4-7为实施例1的琳达(Linda)图的去噪效果图。
具体实施方式
实施例1(采用3×3正方形模板,参见图2):
采用3×3正方形模板时的去噪步骤如下:
(1) 读入一帧图像,令第一个像素点为当前像素点;
(2) 判断所述当前像素点是否为噪声点,若是则按照以下步骤进行处理,若不是则置下一个像素点为当前像素点,接着判断,直到遇到噪声点或算法结束;上述判断噪声点的方法采用区间法,即灰度值位于0-5或250-255之间的为噪声点,否则为信号点;
如果图像边界的像素点为噪声,则对于以该像素点为中心构造的3×3正方形模板中超出图像区域的部分忽略不计,防止超出边界时发生溢出问题;
式中:
i为行标,i=-1,0,1;
j为列标,j=-1, 0, 1 ;
式中:
i=-1,0,1;j=-1, 0, 1 ;
(7)噪音去除:
(8) 如果所述图像中所有的像素点处理完毕,则算法结束;否则,置下一个像素点为当前像素点,返回上述步骤(2),继续处理。
本实施例1的去噪效果(参见图4-1至图4-7及附表1):
由图4-1至4-7的图像显示可知,本实施例1在去除噪音(尤其是局部变化急剧的范围内噪声)有较好的效果。数据显示,实施例1较标准中值滤波均方误差小而峰值信噪比大,尤其是均方误差在50%噪声情况下差别达到27倍(见表1),能够达到除噪的一个较理想效果。
表1 峰值信噪比
实施例2 (采用5×5正方形模板,参见图3):
实施例2与实施例1不同之处如下:
i= -2,-1, 0,1,2;
j= -2,-1, 0,1,2;
为5×5正方形模板中像素点(i,j)的灰度;
其中:i= -2,-1, 0,1,2;
j= -2,-1,0, 1,2 ;
a. 第一轮灰度证据合成:
首先将第1行与第5行灰度证据进行合成:
式中:j= -2,-1, 0,1,2 ;
合成后的灰度证据为5个灰度证据,构成新1行灰度证据;
然后,将第2行与第4行灰度证据进行合成:
b.第二轮灰度证据合成(将新构成的两行灰度证据进行合成):
因为只剩这两行灰度证据,所以直接计算灰度证据合成值m 0:
Claims (1)
1.一种基于证据理论的图像去噪方法,其特征在于步骤如下:
(1) 读入一帧图像,令第一个像素点为当前像素点;
(2) 判断所述当前像素点是否为噪声点,若是则按照以下步骤进行处理,若不是则置下一个像素点为当前像素点,接着判断,直到遇到噪声点或算法结束;上述判断噪声点的方法采用区间法,即灰度值位于0-5或250-255之间的为噪声点,否则为信号点;
(3) 以噪声点为中心,在所述图像上对目标像素构造一个n×n正方形模板,该模板为在横轴和纵轴上均包含n个像素点,n为除“1”之外的奇数;
如果图像边界的像素点为噪声,则对于以该像素点为中心构造的n×n正方形模板中超出图像区域的部分忽略不计,防止超出边界时发生溢出问题;
式中:
i为行标,i=-N, ……,-2,-1,0,1,2,……N;
j为列标,j=-N, ……,-2,-1,0,1,2,……N;
所述的i为n×n正方形模板中的行标,其排列顺序为从上至下依次排列,其中-N为第1行,N为第n行;
所述的j为n×n正方形模板中的列标,其排列顺序为从左至右依次排列,其中-N为第1列,N为第n列;
(2)
式中:
i为行标,i=-N, ……,-2,-1,0,1,2,……N;
j为列标,j=-N, ……,-2,-1,0,1,2,……N;
第一轮灰度证据合成
(3)
式中:为第一行与第n行的灰度证据合成(即新1行的灰度证据);
b.第二轮灰度证据合成:
将上述第一轮灰度证据合成的新的各行的灰度证据再进行第二轮灰度证据合成,第二轮合成后,得到第二轮合成后的各行的灰度证据;
c.再进行第三轮合成及更多轮的灰度证据合成,直至只剩下两行时为止,将最后剩下的两行进行证据合成得到灰度证据合成值m 0 ;
所述灰度证据合成值m 0 按下述公式计算:
;
(7)噪音去除:
(8) 如果所述图像中所有的像素点处理完毕,则算法结束;否则,置下一个像素点为当前像素点,返回上述步骤(2),继续处理。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310559241.5A CN103593827B (zh) | 2013-11-12 | 2013-11-12 | 一种基于证据理论的图像去噪方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310559241.5A CN103593827B (zh) | 2013-11-12 | 2013-11-12 | 一种基于证据理论的图像去噪方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103593827A true CN103593827A (zh) | 2014-02-19 |
CN103593827B CN103593827B (zh) | 2017-01-11 |
Family
ID=50083952
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310559241.5A Expired - Fee Related CN103593827B (zh) | 2013-11-12 | 2013-11-12 | 一种基于证据理论的图像去噪方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103593827B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105931197A (zh) * | 2016-04-13 | 2016-09-07 | 石家庄铁道大学 | 一种基于模糊度理论的图像去噪方法 |
CN106408548A (zh) * | 2016-04-13 | 2017-02-15 | 石家庄铁道大学 | 一种图像调和平均去噪方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101976435A (zh) * | 2010-10-07 | 2011-02-16 | 西安电子科技大学 | 基于对偶约束的联合学习超分辨方法 |
JP2013065302A (ja) * | 2011-09-16 | 2013-04-11 | Toshiba Corp | 画像処理装置及び方法 |
CN103077505A (zh) * | 2013-01-25 | 2013-05-01 | 西安电子科技大学 | 基于字典学习和结构聚类的图像超分辨率重建方法 |
CN103136727A (zh) * | 2012-12-14 | 2013-06-05 | 西安电子科技大学 | 基于梯度一致性与各向异性正则化的超分辨图像重建方法 |
-
2013
- 2013-11-12 CN CN201310559241.5A patent/CN103593827B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101976435A (zh) * | 2010-10-07 | 2011-02-16 | 西安电子科技大学 | 基于对偶约束的联合学习超分辨方法 |
JP2013065302A (ja) * | 2011-09-16 | 2013-04-11 | Toshiba Corp | 画像処理装置及び方法 |
CN103136727A (zh) * | 2012-12-14 | 2013-06-05 | 西安电子科技大学 | 基于梯度一致性与各向异性正则化的超分辨图像重建方法 |
CN103077505A (zh) * | 2013-01-25 | 2013-05-01 | 西安电子科技大学 | 基于字典学习和结构聚类的图像超分辨率重建方法 |
Non-Patent Citations (1)
Title |
---|
杨海峰 等: "基于证据理论的小波萎缩图像去噪", 《光学技术》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105931197A (zh) * | 2016-04-13 | 2016-09-07 | 石家庄铁道大学 | 一种基于模糊度理论的图像去噪方法 |
CN106408548A (zh) * | 2016-04-13 | 2017-02-15 | 石家庄铁道大学 | 一种图像调和平均去噪方法 |
Also Published As
Publication number | Publication date |
---|---|
CN103593827B (zh) | 2017-01-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107016642B (zh) | 用于对有噪输入图像进行分辨率上调的方法和装置 | |
CN106910169B (zh) | 一种防止边缘模糊的图像椒盐噪声去除方法 | |
CN110349103A (zh) | 一种基于深度神经网络和跳跃连接的无干净标签图像去噪方法 | |
CN103455991B (zh) | 一种多聚焦图像融合方法 | |
CN108564597B (zh) | 一种融合高斯混合模型和h-s光流法的视频前景目标提取方法 | |
CN109064418B (zh) | 一种基于非局部均值的非均匀噪声图像去噪方法 | |
CN101930597B (zh) | 一种基于数学形态学的图像边缘检测方法 | |
CN106408522A (zh) | 一种基于卷积对神经网络的图像去噪方法 | |
CN109377450B (zh) | 一种边缘保护的去噪方法 | |
CN104680483B (zh) | 图像的噪声估计方法、视频图像去噪方法及装置 | |
CN103578092A (zh) | 一种多聚焦图像融合方法 | |
CN107730536B (zh) | 一种基于深度特征的高速相关滤波物体追踪方法 | |
CN103116875A (zh) | 自适应双边滤波图像去噪方法 | |
CN103533214A (zh) | 一种基于卡尔曼滤波和双边滤波的视频实时去噪方法 | |
CN110992288B (zh) | 一种用于矿山井筒环境下的视频图像盲去噪方法 | |
CN103886553A (zh) | 图像非局部均值去噪方法和系统 | |
CN103226824B (zh) | 维持视觉显著性的视频重定向系统 | |
CN107169932A (zh) | 一种适用于中子成像系统图像的基于高斯‑泊松混合噪声模型的图像复原方法 | |
CN103493482A (zh) | 一种提取及优化图像深度图的方法与装置 | |
CN111899200A (zh) | 一种基于3d滤波的红外图像增强方法 | |
CN105608689A (zh) | 一种全景拼接消除图像特征误匹配方法及装置 | |
CN103593827A (zh) | 一种基于证据理论的图像去噪方法 | |
CN105279742B (zh) | 一种快速的基于分块噪声能量估计的图像去噪方法 | |
CN103489201A (zh) | 基于运动模糊信息的目标跟踪方法 | |
CN109242797A (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 | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20170111 Termination date: 20181112 |