自然资源遥感, 2022, 34(3): 17-26 doi: 10.6046/zrzyyg.2021319

技术方法

多特征准则融合的遥感图像脉冲噪声的识别处理

马晓剑,1, 赵法舜1, 刘艳宾2

1.东北林业大学理学院,哈尔滨 150040

2.中国地质大学(北京)地球科学与资源学院,北京 100083

Multi-feature fusion-based recognition and processing of impulse noise in remote sensing images

MA Xiaojian,1, ZHAO Fashun1, LIU Yanbin2

1. College of Science, Northeast Forestry University, Harbin 150040, China

2. School of Earth Sciences and Resources, China University of Geosciences(Beijing), Beijing 100083, China

责任编辑: 李瑜

收稿日期: 2021-09-30   修回日期: 2022-03-6  

基金资助: 中央高校基本科研业务费专项资金项目“证据理论融合算法在图像处理中的研究与应用”(2572018BC21)

Received: 2021-09-30   Revised: 2022-03-6  

作者简介 About authors

马晓剑(1977-),女,副教授,主要从事图像处理研究。Email: mxjzy@nefu.edu.cn

摘要

消除脉冲噪声,获取高质量的遥感图像对应用研究有着重要意义。消除高密度脉冲噪声的同时,保持原有遥感图像的边缘细节信息一直是这一领域中的难题。该文认为被脉冲噪声冲击后的图像会出现不确定性突变,为了解决这种不确定性问题,基于证据理论,利用脉冲噪声的多个特征进行了不确定性建模; 融合了BJS散度和信度熵,给出新的权重分配,得到了新的概率指派; 再根据融合规则和概率转换,给出噪声与信号点的分类依据,从而有效降低了高度冲突发生的可能性。实验结果表明,在噪声密度达到90%以上时,该文提出的方法仍然有效,且在消噪后的遥感图像中对不同地物信息的细节保持良好。

关键词: 证据理论; 不确定性建模; 融合规则; 高度冲突; 遥感图像

Abstract

Eliminating impulse noise of high-quality remote sensing images is of great significance for applied research. It has always been a challenge to eliminate high-density impulse noise while remaining detailed information on edges in original remote sensing images. This study concluded that uncertain changes will appear when a remote sensing image is corrupted by impulse noise. Given this, an uncertainty model based on the evidence theory was constructed using multiple features of impulse noise. The BJS divergence and the reliability entropy were fused into the model to obtain new weights and a new probability assignment. Then, the classification between noise and signals was given according to fusion rules and probability transformation, thus effectively reducing the possibility of high-level conflicts. The experimental results show that the classification method proposed in this study is effective even when the noise density is up to over 90% and can well maintain detailed information on different ground objects in the denoised remote sensing images.

Keywords: evidence theory; uncertainty modeling; fusion rules; highly conflict; remote sensing image

PDF (5845KB) 元数据 多维度评价 相关文章 导出 EndNote| Ris| Bibtex  收藏本文

本文引用格式

马晓剑, 赵法舜, 刘艳宾. 多特征准则融合的遥感图像脉冲噪声的识别处理[J]. 自然资源遥感, 2022, 34(3): 17-26 doi:10.6046/zrzyyg.2021319

MA Xiaojian, ZHAO Fashun, LIU Yanbin. Multi-feature fusion-based recognition and processing of impulse noise in remote sensing images[J]. Remote Sensing for Land & Resources, 2022, 34(3): 17-26 doi:10.6046/zrzyyg.2021319

0 引言

遥感图像是用来描述地表信息的重要数据源,能否获取清晰、高质量的遥感图像直接关系到有用信息的提取、分析和应用。但是在采集、切换和传送遥感图像时,常常因为仪器性能、操作手段或者成像环境,导致图像被脉冲噪声所冲击,造成图像质量下降。因此,在利用遥感图像解决深层问题之前,选择合适的方法在去除脉冲噪声的同时,有效保留遥感图像原有的有用信息(如边缘和细节等)就显得至关重要。

遥感图像去噪的方法可以基于不同的理论模型,比如概率统计、偏微分方程、稀疏表示以及多尺度分析等[1-6]。DBA算法[5]、NSBMF[6]等非线性方法是在统计的基础上添加噪声判断,而文献[2]提出了一种基于梯度倒数的自适应开关算法对噪声进行判断并消除。此外,为了解决遥感图像在滤波过程中丢失细节信息的问题,文献[3]提出使用组合滤波方法,在利用多尺度分析和中值滤波后,提取边缘信息,将提取结果再融合到消噪后的遥感影像中,以此提升细节信息。以上方法的共同特点是先对随机脉冲噪声与信号点进行分类[4],再对噪声点进行消除,那么分类的准确性会直接影响着滤波的效果。

事实上这种分类问题可以视作是一种不确定性问题,因为被脉冲噪声冲击的遥感图像在某些像素点上会发生不确定性的突变,这种不确定性主要体现在突变像素强度的不确定性、突变位置的不确定性和突变数量的不确定性。D-S证据理论具有表达“不确定”和“不知道”的能力,在没有先验信息的情况下它对不确定性和非精确性进行建模有着不错的灵活性和有效性[7]。因此,它比一般概率理论的适用范围更加广泛,作为解决不确定性问题的有力工具,它已经被应用到许多实际领域,如目标识别、遥感分类和图像分割等[8-10]。D-S证据理论已经成为不确定性信息处理与信息融合领域的重要理论,但是很少有人将该理论与图像脉冲噪声和信号点的分类问题结合起来。

直到2018年Zhang等[11]提出了ASMF-DBER算法,同时考虑脉冲噪声间断特征和极值特征2种准则进行了不确定性建模,创新性地将D-S证据理论应用到了噪声与信号的分类上,并取得了不错的结果。但是在证据高度冲突的情况下,仅使用2个特征和冲突系数无法准确显示融合的效果。比如,在脉冲噪声密度较高时,图像由于受到大量脉冲噪声的冲击而遭到严重损坏,ASMF-DBER算法中的2个特征准则融合后的结果出现了违背直觉或者难以分类的情况,本文考虑添加第3种特征来降低遥感图像中出现的证据高度冲突发生的可能性。所谓“高度冲突”的情况在本文中主要表现为以下几种情况: 一是遥感图像脉冲噪声密度过大; 二是待分类的像素处在遥感图像的边缘位置或者处在图像明暗区域的交界处; 三是待分类的像素恰好体现遥感图像的细节信息。一直以来这些情况下的噪声消除都是难以解决的。

为了更加有效地刻画遥感图像的固有结构、边缘以及纹理等细节特征,并尽可能地将脉冲噪声从遥感图像中去除掉,本文把脉冲噪声的3个特征准则: 极值性、间断性和相似性作为证据,提出多特征准则模型,再根据BJS散度和信度熵,构建新的权重分配以解决证据融合出现的高度冲突的问题,然后综合评估每个像素被判断为噪声点的概率大小,即通过概率决策完成分类。另外,在消除噪声的阶段,采用的是自适应中值滤波方法。需要指出的是本文的重点是脉冲噪声与信号的分类,分类的准确性是本文重点研究的内容。实验结果证实本文算法在能够提高图像滤波后的信噪比的同时,还做到了有效保留遥感图像原有的边缘细节等信息,即使是在高度冲突的情况下也有理想的结果。

1 不确定性噪声模型及证据理论

1.1 脉冲噪声模型

为了描述脉冲噪声的不确定性,本文选择文献[12]中的噪声模型,当图像(以8位灰色图像为研究对象)被脉冲噪声破坏后,噪声随机出现在像素集合N1={0,1,2,,l}N2={255-l,255-(l-1),,255}中,其概率密度函数如下:

f(xi,j)=p2   (xi,jN1)1-p  (xi,j=Si,j)p2   (xi,jN2)

式中:xi,jSi,j分别为噪声图像和原始图像中处在第i行,第j列的位置的像素点;p为图像被脉冲噪声破坏的概率。被脉冲噪声污染的遥感图像的像素明显偏离正常值,主要集中在0或者255附近,具有一定的不确定性,因此,本文选取l=2进行研究,即噪声像素强度为0,1,2,253,254或255。与l=0时,脉冲噪声强度仅为0或者255的椒盐噪声相比,本文给出的噪声模型增大了噪声强度的不确定性。此外,由脉冲噪声产生的机制不难发现,噪声出现位置和数量也具有不确定性。可见,式(1)确定的脉冲噪声模型更能反映遥感图像被脉冲噪声污染的实际情况。

1.2 D-S证据理论

n个两两互斥的元素组成的有限的完备集合称为辨识框架,记为Θ。本文中的辨识框架Θ={N,S},N为脉冲噪声,S为非脉冲噪声(以下称为信号)。对于辨识框架Θ,如果它的幂集2Θ到[0,1]的映射m,即m:2Θ[0,1],满足以下条件:m(Ø)=0,AΘm(A)=1,称该映射m为基本概率指派,记作BPA。m(A)被视为准确分配给焦元A的信度。

在辨识框架Θ上,设有n个BPA函数m1,m2,,mn,焦元分别为A1,A2,,An,并假设这n个BPA函数是相互独立的,则BPA函数之间的Dempster融合规则定义如下:

m(A)=1Ki=1nAi=Ai=1nmi(Ai) (AØ)0        (A=Ø)

式中冲突系数K=inAiØinmi(Ai)是一个常数,表示BPA函数之间的冲突,称其为冲突系数。Dempster的融合规则仅适用于冲突系数K<1的情况。

概率转换函数[13]是指假设m(A)是辨识框架Θ下的BPA函数,Pignistic概率转换函数将一个BPA转换为一个概率测度Betp,公式为:

Betp(θi)=θiAm(A)A,θiΘ

式中A为焦元A的基数; θi为脉冲噪声N或非脉冲噪声S。一般地,通过辨识框架ΘBetp的最高值来进行最终的决策。

2 多特征准则融合的分类算法

2.1 脉冲噪声的多特征准则不确定性建模

与相邻信号相比,被脉冲噪声污染的像素有显著的特征,比如: 极值性、相似性以及间断性。本文根据这些特征利用区间距离、相似度函数以及像素的相对差异分别构造特征准则证据m1,m2m3

2.1.1 极值性特征准则BPA的构造

由噪声模型可以看出,脉冲噪声像素的强度主要分布在0或者255附近,接近或者达到极值。为了将噪声和信号点进行分类,本文利用区间距离来描述像素的强度和极值之间的接近程度,定义为[14]:

d(a,b)=3(a1+a2)-(b1+b2)22+(a2-a1)-(b2-b1)22

式中ab为区间,即区间a=(a1,a2)和区间b=(b1,b2),则像素x与图像之间的区间距离可以表示为:

d(I,[x,x])=(Imax-Imin-2x)2+(Imax-Imin)2

这里I=[Imin,Imax]是图像的强度区间,IminImax分别为图像的最小强度和最大强度,x[Imin,Imax]

图1所示,当x取[0,255]的中值时,d(I,[x,x])达到最小值; 当x接近极值0或255时,区间距离变大。因此,d(I,[x,x])可以用来描述像素强度x与极值IminImax之间的接近程度。根据区间距离的这一性质,本文构造BPA函数m1,公式如下:

m1(N)=dc-d0+εdext-d0+ε    m1(S)=1-dc-dmeddext-dmed   m1(Θ)=1-m1(N)-m1(S)

式中:dc为待分类像素的强度和区间I之间的距离;d0为在当前窗口中的像素和I之间的最小距离;dextI和极值(IminImax)的距离;dmedII的中值的距离。在噪声密度高的情况下,m1(N)应该接近于1,但是文献[11]中m1(N)趋近于0,这样就导致了噪声点被判断成信号的错误分类结果,本文将文献[11]的m1函数做了改进,用ε(ε=0.001)去调整BPA函数,避免了这种情形下的高度冲突结果的发生。

图1

图1   不同像素的d(I,[x,x])变化图

Fig.1   d(I,[x,x]) of pixels


2.1.2 相似性特征准则BPA的构造

往往图像的细节信息与周围信号既有差异性又有相似性,只用极值性特征准则进行分类还不够准确,本文用待分类像素与窗口内其他信号的差异,作为分类的又一依据来构建相似性特征准则BPA函数。

若以(i,j)为中心的窗口是WF(i,j)=xi-s,j-t-(WF-1)2s,t(WF-1)2,s,tZ,其窗口大小为WF×WF,则

rg(xi,j)=xi,j-xi-s,j-t

式中: rg(xi,j)为像素差,g=1,2,,WF×WF;xi,j(i,j)处的待分类像素强度;xi-s,j-t(i-s,j-t)处的像素强度。令SIMijxi,j与窗口内像素的相似度[15],公式为:

SIMij=g=1WF×WFexp-rg(xi,j)22σμ2WF×WF

式中σμ取当前窗口区域像素的均值。则SIMij表征了xi,j与周围像素的相似程度。

利用相似度构造特征准则BPA如下:

m2(N)=μ·(1-SIMij)m2(S)=μ·SIMij   m2(Θ)=1-μ    

式中μSIMij的权重,介于0~1之间。从图2中可以看出: 当rg的绝对值越小时,SIMij越大,说明xi,j的纹理、强度等特征与周围像素越相似,即它没有被脉冲破坏的可能性越大,则被判断成是信号点的概率就越大; 当rg绝对值增大时,SIMij减小,xi,j被判断成是脉冲噪声的概率变大。可见用SIMij构造BPA函数来对脉冲噪声与信号分类是合理的。通过实验,本文取μ=0.9效果最佳,容易证明m2满足BPA函数的定义。

图2

图2   rgSIMij关系曲线

Fig.2   The curve of relation between rg and SIMij


2.1.3 间断性特征准则BPA的构造

处在明暗交界处的信号往往比其相邻信号的强度高或低很多,这体现了间断性的特征,同时也是脉冲噪声的特点,因此这些信号往往被错误地分类成噪声。本文构建了间断性特征准则BPA,记作m3函数,以减小分类误判的概率,公式如下:

m3(N)=γg=1nrg(xi,j)n(IWmax-IWmin)m3(S)=γ(1-g=1nrg(xi,j)n(IWmax-IWmin))m3(Θ)=1-γ

式中:IWmax为当前窗口像素的最大值,IWmin为当前窗口像素的最小值;n为当前窗口WF(i,j)中除去中心像素的像素个数;γ为权重,介于0~1之间。

当窗口固定时,式(10)中的g=1nrg(xi,j)表示待分类像素强度与周围其他像素强度差异的总和,g=1nrg(xi,j)的值越小则说明其间断程度越小,相反地,g=1nrg(xi,j)的值越大则说明它与其他周围像素强度间断程度较大,被判定为脉冲噪声的可能性更大; 反之,根据待分类像素与其他周围像素强度的局部连续特征,可以更大概率地将其判断为信号点。又注意到rg(xi,j)<(IWmax-IWmin),g=1nrg(xi,j)n(IWmax-IWmin)是介于[0,1]的数,容易证明m3满足BPA函数的定义,本文选取γ=0.9

2.2 基于信度熵和BJS散度的证据融合

虽然根据脉冲噪声的极值性、相似性、间断性这3个特征进行了不确定性证据建模,但是证据之间的冲突可能还会出现,为了减少证据冲突的可能性,本文提出在利用证据融合规则之前,对特征准则BPA函数进行修正处理,利用BJS散度与信度熵对证据进行重新分配,以校正高度冲突情形下噪声与信号分类的错误结果。

Ai为信度函数m的假设,m1m2为在同一辨识框架Θ上的两个BPA函数,则m1m2之间的BJS散度[16]可表示为:

BJS(m1,m2)=12Sm1,m1+m22+Sm2,m1+m22

式中:S(m1,m2)=im1(Ai)logm1(Ai)m2(Ai);imj(Ai)=1;i=1,2,,M;j=1,2。当BPA函数的值为0时,将使用10-6来替换0值。

根据BJS散度建立证据之间的差异度量矩阵M,即:

M=0BJS12BJS210BJS1nBJS2nBJSn1BJSn20

则证据可信度Rei可由BJSij来表示,记Rei=1Acdi,i=1,2,,M。其中Acdi=j=1,jiMBJSijM-1,i=1,2,,M

事实上,BJS散度的显著特征之一是可以为每种概率分布分配不同的权重,它代表了证据间可靠性的差异大小,2个证据间越相似,BJS散度就越大,尤其在多于2个证据时,它对分类结果会产生至关重要的作用。

此外,当分类的不确定性程度很高时,信度熵越大,得到证据的支持度越少,可以用它来进一步表示信号与噪声分类的不确定性。所以令Ai为信度函数m的假设,则集合Ai的信息量为:

ICi=e-Ed

式中Ed=-im(Ai)logm(Ai)2Ai-1,i=1,2,,M为集合Ai的信度熵[17]

综上,为降低高度冲突发生的可能性,提高分类效果,本文利用BJS散度建立证据可信度Rei,再根据信度熵,给出证据信息量ICi,最后融合两者作为证据的新权重,实现对证据的重新分配,具体步骤如下:

1)BJS散度计算3个证据的差异度量矩阵Mij,1i3,1j3

2)计算证据可信度Rei,i=1,2,3

3)由信度熵计算每个证据的信息量ICi,i=1,2,3

4)由可信度和信息量融合生成新权重Wfi=Rei×ICi,i=1,2,3

5)将权重归一化分配给原始证据得到新的证据m'(Ai)=i=13(Wfi'×mi(Ai)),其中Wfi'=Wfii=13Wfi,i=1,2,3

于是再将修改后的证据采用Dempster的融合规则进行自身间的融合,并使用Pignistic概率转换(式(3))得到概率测度Betp,并将其作为最终脉冲噪声与信号点分类的依据。如果Betp>0.5,则判断该像素点为脉冲噪声,反之则为信号点。本文提出的多特征准则的融合分类算法具体过程见图3

图3

图3   多特征准则的建模与融合

Fig.3   Modeling and fusion process of multi-feature criterion


本文算法采用先对脉冲噪声与信号分类,再进行噪声消除的思路,提出了基于IBDND[18]的消噪改进方法,步骤如下:

1)将过滤窗口的初始大小WF设置为3,并将最大窗口大小Wmax设置为11。

2)设置当前窗口大小WF×WF,并以(i,j)的目标像素为中心设置滤波窗口WF(i,j)

3)判断在当前滤波窗口中是否满足条件:Nu>14PNt。如果满足,则转到5),否则进行4)。

4)将过滤窗口的大小扩展为(WF+1)×(WF+1),然后重复2)和3),直到当前过滤窗口大小达到最大窗口11×11。

5)对当前过滤窗口应用近似中值像素替换,输出像素值为:

Zij=Yij+1Dk=1NuVu(k)-Yijdist(k)

式中:Nu为当前窗口脉冲噪声像素个数;P为脉冲噪声密度;Nt为当前窗口像素总数;Zij为在窗口信号点中值Yij的基础上增添了基于距离索引的调整项;Vu为当前窗口信号点像素;dist(k)Vu中第k个值的信号与位置(i,j)处的噪声像素之间的空间距离,公式为:

dist(k)=(a1+a2)-(b1+b2)22+(a2-a1)-(b2-b1)22

式中:a1,a2为第k个像素的行与列的索引;b1,b2为中心像素在当前窗口的行与列的索引;D=k=1Nu1dist(k)

3 实验与分析

为了验证本文算法的有效性,文章从SIRI-WHU遥感影像中的公园、住宅区、工业区、和河坝等数据库[19-21]里随机选取了6幅含有不同类型地物信息的遥感图像(见图4)进行实验。这些图像具有区域平滑、色调对比明显或者背景复杂、纹理信息丰富等不同特征。采用的不确定性噪声模型见式(1)。为说明算法的优势,本文还与DBA,NSBMF,IBDND,SAMF[22]以及ASMF-DBER这5种算法进行了比较。实验主要从“高冲突”角度出发,对高噪声密度、边缘、纹理细节等情况进行视觉分析和定量分析。定量指标选择常用的结构相似度(structural similarity index,SSIM)、准确率(accuracy rate,AR)和峰值信噪比(peak signal to noise ratio,PSNR),指标公式如下:

AR=NS-FP-FNNS
PSNR=10lg2552MNi=1Mj=1N(xij︿-xij)2
SSIM=(2μx︿μx+CN1)(2σx︿x+CN2)(μx︿2+μx2+CN1)(σx︿2+σx2+CN2)

式中: NS为像素总数; FP为误检数; FN为漏检数;MN为图像尺寸;x为初始图像;x︿为过滤图像;μx︿μx分别为xx︿的均值;σx︿σx分别为xx︿的方差;σx︿xxx︿的协方差;CN1CN2为常数。

图4

图4   不同地物信息的遥感图像

Fig.4   Remote sensing images for different ground object


由于本文算法分2步进行,即先将遥感图像中的脉冲噪声与信号分类,再消除噪声,因而先选择了同为2步消噪的IBDND和ASMF-DBER算法进行对比。由式(16)可以看出AR指标可以衡量分类的准确性,其值越大说明遥感图像中信号与噪声的分类越准确;SSIM指标则表示消噪后的遥感图像与原图的接近程度。因此,以图4(a)为例,表1给出这2个指标下的对比结果。

表1   不同算法的ARSSIM

Tab.1  AR and SSIM of different algorithms

噪声
密度/%
ASMF-DBMRIBDND本文算法
SSIMAR/%SSIMAR/%SSIMAR/%
100.98099.610.98199.430.98199.65
200.95899.590.93099.380.96099.33
300.93299.570.78499.200.92699.38
400.89999.520.54898.480.89699.48
500.86399.510.31596.930.85199.62
600.81799.480.16694.120.80499.76
700.75099.590.08989.590.75699.76
800.65599.620.04584.140.67899.73
900.53199.820.01875.810.55399.84

新窗口打开| 下载CSV


表1中可以看出,在低噪声密度时,2种指标下的IBDND算法稍逊一筹,当噪声密度达到50%后,即“高度冲突”体现越明显时,本文算法的分类准确率AR表现就越突出,这意味着更多的噪声被准确识别出来。另外本文算法的SSIM在“高度冲突”下表现稳定,说明消噪后的遥感图像结构变化小于其他算法,细节保持更优。

表1不难发现,在噪声密度较低时,ASMF-DBER算法的结构相似度略微优于本文算法,但是从定性视觉直观的角度(图5)能够发现,当噪声密度为50%时,消噪后选定明暗交界处,放大横跨河水两岸的建筑物的部分细节(图5(e)),可以明显看出ASMF-DBER算法对梯度变化较大的信息保留要稍差,而本文算法符合地理要素的空间自相关性,消噪后遥感图像中的建筑物仍能保持良好的连续性,失真度小,边缘更流畅,对比度损失小。为了进一步说明本文算法对图像边缘或明暗交替的位置具有较高的还原度,本文对河坝图的局部信息(20×20的像素区域)进行误差分析,将未污染噪声的细节图与消噪后的细节图做差,生成像素误差图,如果误差分析图越接近0附近的颜色,那么消噪后的误差就会越小。由图6不难发现相比于ASMF-DBER算法,本文算法误差波动幅度较小,图像颜色显示误差值分布更接近0附近,说明本文算法处理后的细节结果更接近原图像,从而再次说明本文算法对遥感图像的细节保持好,后期应用处理会更具优势。

图5

图5   噪声密度为50%时ASMF-DBER算法与本文算法的对比图

Fig.5   Comparison between ASMF-DBER algorithm and the proposed algorithm at the noise level of 50%


图6

图6   50%噪声密度下的局部信息对比图

Fig.6   Comparison of local information at the noise level of 50%


为了进一步从视觉角度说明本文算法适合不同地物类型的遥感图像的去噪,本文将图4(b),(c),(d),(e)以及(f)添加10%~90%不同密度的脉冲噪声,考察遥感图像被不同噪声密度干扰时的恢复能力。这5幅不同类型的遥感图像中住宅区和立交桥的遥感图像背景相对复杂,像素反差小,细节信息较多; 公园和工业区的遥感图像类型像素强度高低相间,背景相对平滑,图像的连续性较为复杂; 云层图像色调反差大,纹理信息明显。原图及各算法处理后的效果如表2所示。从中可以看出,SAMF算法由于分类能力最差,明显有大量噪声残留,故仅适用于低噪声密度的情况; DBA算法处理后的图像过于模糊,影响后续的地物信息的提取; NSBMF算法对于细节信息多的遥感图像来说,该方法的弱点就暴露出来,比如边缘连续性差、细节清晰度低、明暗交界的轮廓位置出现了毛糙; 而本文算法对于“高冲突”处理结果理想,噪声提取准确,对于遥感图像中不同类型的地物信息还原性强。同时从表2中还可以发现,公园和云层的遥感图像中的白色高亮部分经过不同恢复算法后损失严重,这是因为白色高亮部分的像素强度接近脉冲噪声的像素强度,在脉冲噪声的识别过程中容易产生误判。本文利用多特征准则融合算法提取脉冲噪声,有效降低了误判发生的可能性,白色高亮的细节部分得到更好的保留。

表2   噪声密度依次为10%,30%,50%,70%和90%的滤波效果对比图

Tab.2  Comparison of filtering effect under 10%,30%,50%,70% and 90% noise density

新窗口打开| 下载CSV


为验证不同噪声密度条件下本文算法均有良好的还原能力,当噪声密度由10%变化到90%时,图7为应用5种算法对图4不同类型遥感图像分别处理后的定量指标PSNR的结果。很显然,随着噪声密度的升高,PSNR指标普遍下降,IBDND算法下降速度最快,本文算法的滤波优势明显,即使噪声密度达到90%,相比其他算法,本文算法处理过的遥感图像对噪声的分类仍然有效。

图7

图7   遥感图像的PSNR对比图

Fig.7   Comparison of the PSNR of remote sensing images


4 结论

脉冲噪声的存在极大影响着遥感图像的信息判读与分析,为了实现对遥感图像的高质量复原,扩展遥感图像的应用范围和应用价值,本文提出了多特征准则融合的脉冲噪声识别方法。在文献[11]的基础上,本文增加了脉冲噪声的相似性特征,根据脉冲噪声的相似性、极值性和间断性特征进行了新的不确定性建模; 同时利用BJS散度和信度熵对证据进行了新的权重分配,根据D-S证据理论进行了证据融合,避免冲突信息的发生。通过实验得到以下结论:

1)针对被脉冲噪声损坏后的SIRI-WHU影像中工业区、公园等不同类别的遥感图片,本文算法均表现出很强的恢复能力,还原后的遥感影像符合地理要素的空间自相关性,图像清晰且对比度损失小,说明本文算法可以广泛适用于不同类型遥感图像的去噪处理。

2)即使在信息“高冲突”的情况下,本文算法仍然能保证脉冲噪声被准确的识别,完成图像的去噪。在脉冲噪声密度达到90%时,本文算法处理后的遥感图像仍能保持原有的纹理、边缘以及明暗交界等细节信息,图像的连续性仍能保持较佳的水平,说明本文提出的多特征准则信息融合算法适合对细节信息丰富或者被高密度噪声干扰的图像进行恢复处理。

参考文献

赵洪臣, 周兴华, 彭聪, .

一种去除遥感影像混合噪声的集成BM3D方法

[J]. 武汉大学学报(信息科学版), 2019, 44(6):925-932.

[本文引用: 1]

Zhao H C, Zhou X H, Peng C et al.

An integrated BM3D method for removing mixed noise in remoting sensing image

[J]. Geomatics and Information Science of Wuhan University, 2019, 44(6):925-932.

[本文引用: 1]

汪贵平, 杜晶晶, 宋京, .

基于梯度倒数的无人机遥感图像融合滤波方法

[J]. 科学技术与工程, 2018, 18(31):190-194.

[本文引用: 2]

Wang G P, Du J J, Song J, et al.

A fusion filter method for unmanned aerial vehicle remote sensing image based on gradient inverse

[J]. Science Technology and Engineering, 2018, 18(31):190-194.

[本文引用: 2]

朱建军, 周靖鸿, 周璀, .

一种新的去除遥感影像混合噪声组合滤波方法

[J]. 武汉大学学报(信息科学版), 2017, 42(3):348-354.

[本文引用: 2]

Zhu J J, Zhou J H, Zhou C, et al.

A new combination filtering method to removing mixed noise of remote sensing images

[J]. Geomatics and Information Science of Wuhan University, 2017, 42(3):348-354.

[本文引用: 2]

刘帅. 基于分层稀疏学习和协同表示的高光谱图像去噪和分类[D]. 西安: 西安电子科技大学, 2016.

[本文引用: 2]

Liu S. Hierarchical sparse learning and collaborative representation for hyperspectral imagery restoration and classification[D]. Xi’an: Xidian University, 2016.

[本文引用: 2]

Srinivasan K S, Ebenezer D.

A new fast and efficient decision-based algorithm for removal of high-density impulse noises

[J]. IEEE Signal Processing Letters, 2007, 14(3):189-192.

DOI:10.1109/LSP.2006.884018      URL     [本文引用: 2]

Jayaraj V, Ebenezer D.

A new switching-based median filtering scheme and algorithmfor removal of high density salt and pepper noise in images

[J]. Journal on Advances in Signal Processing, 2010(1):409-413.

[本文引用: 2]

Dempster A P.

Upper and lower probabilities induced by a multivalued mapping

[J]. The Annals of Mathematical Statistics, 1967, 38(2):325-339.

DOI:10.1214/aoms/1177698950      URL     [本文引用: 1]

蒋雯. 邓鑫洋. D-S证据理论信息建模与应用[M]. 北京: 科学出版社, 2018.

[本文引用: 1]

Jiang W, Deng X Y. D-S evidence theory information modeling and application[M]. Beijing: Science Press, 2018.

[本文引用: 1]

童涛, 杨桄, 李昕, .

基于D-S证据理论的多特征融合SAR图像目标识别方法

[J]. 国土资源遥感, 2013, 25(2):37-41.doi: 10.6046/gtzyyg.2013.02.07.

[本文引用: 1]

Tong T, Yang G, Li X, et al.

Recognition method of multi-feature fusion based on D-S evidence theory in SAR image

[J]. Remote Sensing for Land and Resources, 2013, 25(2):37-41.doi: 10.6046/gtzyyg.2013.02.07.

[本文引用: 1]

李华朋, 张树清, 孙妍.

证据理论结合遥感分类数据能力定量评价研究

[J]. 国土资源遥感, 2011, 23(1):26-32.doi: 10.6046/gtzyyg.2011.01.05.

[本文引用: 1]

Li H P, Zhang S Q, Sun Y.

The quantitative evaluation of remoting sensing data for supervised evidential classification

[J]. Remote sensing for Land and Resources, 2011, 23(1):26-32.doi: 10.6046/gtzyyg.2011.01.05.

[本文引用: 1]

Zhang Z, Han D, Dezert J, et al.

A new adaptive switching median filter for impulse noise reduction with predetection based on evidential reasoning

[J]. Signal Processing, 2018(147):173-189.

[本文引用: 4]

Ng P E, Ma K.

A switching median filter with boundary discriminative noise detection for extremely corrupted images

[J]. IEEE Transactions on Image Processing, 2006, 15(6):1506-1516.

DOI:10.1109/TIP.2005.871129      URL     [本文引用: 1]

Han D, Dezert J, Duan Z.

Evaluation of probability transformations of belief functions fordecision making

[J]. IEEE Transactions on Systems,Man,and Cybernetics, 2016, 46(1):93-108.

[本文引用: 1]

Irpino R V.

Dynamic clustering of interval data using a wasserstein based distance

[J]. Pattern Recognition Letter, 2008, 29(11):1648-1658.

DOI:10.1016/j.patrec.2008.04.008      URL     [本文引用: 1]

钱晓亮, 郭雷, 余博.

基于目标尺度的自适应高斯滤波

[J]. 计算机工程与应用, 2010, 46(12):14-16.

[本文引用: 1]

Qian X L, Guo L, Yu B.

Adaptive Gaussian filter based on object scale

[J]. Computer Engineering and Applications, 2010, 46(12):14-16.

[本文引用: 1]

Xiao F.

Multi-sensor data fusion based on the belief divergence measure of evidencesand the belief entropy

[J]. Information Fusion, 2019,(46):23-32.

[本文引用: 1]

Deng Y.

Deng entropy

[J]. Chaos,Solitons & Fractals, 2016(91):549-553.

[本文引用: 1]

Jafar I F, AlNa’mneh R A, Darabkh K A.

Efficient improvements on the BDND filtering algorithm for the removal of high-density impulse noise

[J]. IEEE Transactios on Image Processing, 2013, 22(3):1223-1232.

[本文引用: 1]

Zhao B, Zhong Y, Xia G S, et al.

Dirichlet-derived multiple topic scene classification model fusing heterogeneous features for high spatial resolution remote sensing imagery

[J]. IEEE Transactions on Geoscience and Remote Sensing, 2016, 54(4):2108-2123.

DOI:10.1109/TGRS.2015.2496185      URL     [本文引用: 1]

Zhao B, Zhong Y, Zhang L, et al.

The fisher kernel coding framework for high spatial resolution scene classification

[J]. Remote Sensing, 2016, 8(2):157-176.

DOI:10.3390/rs8020157      URL     [本文引用: 1]

Zhu Q, Zhong Y, Zhao B, et al.

Bag-of-visual-words scene classifier with local and global features for high spatial resolution remote sensing imagery

[J]. IEEE Geoscience and Remote Sensing Letters, 2016, 13(6):747-751.

DOI:10.1109/LGRS.2015.2513443      URL     [本文引用: 1]

Haidi I, Nicholas S P K, Theam F N.

Simple adaptive median filter for the removal of impulse noise from highly corrupted images

[J]. IEEE Transactions on Consumer Electronics, 2008, 544(4):1920-1927.

[本文引用: 1]

/

京ICP备05055290号-2
版权所有 © 2015 《自然资源遥感》编辑部
地址:北京学院路31号中国国土资源航空物探遥感中心 邮编:100083
电话:010-62060291/62060292 E-mail:zrzyyg@163.com
本系统由北京玛格泰克科技发展有限公司设计开发