基于后验概率空间变化向量分析的NSCT高分辨率遥感影像变化检测

doi:10.6046/zrzyyg.2023079

基于后验概率空间变化向量分析的NSCT高分辨率遥感影像变化检测

宋嘉鑫^,¹, 李轶鲲^,¹^,²^,³, 杨树文¹^,²^,³, 李小军¹^,²^,³

1.兰州交通大学测绘与地理信息学院,兰州 730070

2.地理国情监测技术应用国家地方联合工程研究中心,兰州 730070

3.甘肃省地理国情监测工程实验室,兰州 730070

NSCT-based change detection for high-resolution remote sensing images under the framework of change vector analysis in posterior probability space

SONG Jiaxin^,¹, LI Yikun^,¹^,²^,³, YANG Shuwen¹^,²^,³, LI Xiaojun¹^,²^,³

1. Faculty of Geomatics,Lanzhou Jiaotong University, Lanzhou 730070, China

2. National-Local Joint Engineering Research Center of Technologies and Applications for National Geographic State Monitoring, Lanzhou 730070, China

3. Gansu Provincial Engineering Laboratory for National Geographic State Monitoring, Lanzhou 730070, China

通讯作者: 李轶鲲(1978-),男,博士,副教授,主要从事影像处理方面的研究。Email:liyikun2003@hotmail.com。

责任编辑: 张仙

收稿日期: 2023-03-21 修回日期: 2023-06-20

基金资助:

国家重点研发计划项目“边海重点区域安全态势异常感知与互联互通分析技术”(2022YFB3903604)
国家自然科学基金项目“西北重点城市彩钢板建筑群与产业园区时空关联关系”(42161069)

Received: 2023-03-21 Revised: 2023-06-20

作者简介 About authors

宋嘉鑫(1997-),男,硕士研究生,研究方向为遥感图像处理。Email: 3012571995@qq.com。

摘要

非下采样轮廓波变换(non-subsampled contourlet transform,NSCT)和变化向量分析法(change vector analysis,CVA)在高分辨率遥感影像变化检测中,会因不同地物的变化幅度有显著差异,而在单一阈值下无法保证较高的检测精度。为此,文章在后验概率空间变化向量分析(change vector analysis in posterior probability space,CVAPS)的框架下,提出了一种基于模糊C均值聚类(fuzzy C-means, FCM)和简单贝叶斯网络(simple Bayesian network,SBN)的NSCT变化检测方法(FCM-SBN-CVAPS-NSCT)。该方法首先将FCM与SBN耦合,计算出后验概率变化强度图; 之后,通过NSCT将后验概率变化强度图分解为不同尺度和方向的子图,通过保留高频子图中的细节并消除噪声,优化了重构后的后验概率变化强度图,实现了后验概率空间下的多尺度、多方向的变化检测,最终提高了变化检测的精度。实验结果表明,所提方法在3个研究区中得到的Kappa系数比FCM-SBN-CVAPS分别高出了0.100 9,0.056 6和0.067 4,具有一定的优越性。

关键词： 简单贝叶斯网络; 模糊C均值聚类; 后验概率空间变化向量分析; 非下采样轮廓波变换

Abstract

In the change detection for high-resolution remote sensing images, non-subsampled contourlet transform (NSCT) and change vector analysis (CVA) cannot ensure high detection accuracies under single thresholds due to significantly different changes in surface features. Hence, under the framework of change vector analysis in posterior probability space (CVAPS), this study proposed a NSCT-based change detection method combining fuzzy C-means (FCM) clustering and a simple Bayesian network (SBN): the FCM-SBN-CVAPS-NSCT method. First, the proposed method coupled FCM with an SBN to generate a change intensity map in posterior probability space. Then, the change intensity map was decomposed into submaps of different scales and directions through NSCT. The reconstructed change intensity map was optimized by preserving the details and eliminating noise in the high-frequency submaps. Finally, the multi-scale and multi-directional change detection in posterior probability space was achieved, enhancing the change detection accuracy. As indicated by the experimental results, the Kappa values obtained by the proposed method for three study areas were 0.100 9, 0.056 6, and 0.067 4 higher than those derived from the FCM-SBN-CVAPS method, demonstrating certain superiority.

Keywords： simple Bayesian network; fuzzy C-means clustering; change vector analysis in posterior probability space; non-subsampled contourlet transform

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

本文引用格式

宋嘉鑫, 李轶鲲, 杨树文, 李小军. 基于后验概率空间变化向量分析的NSCT高分辨率遥感影像变化检测[J]. 自然资源遥感, 2024, 36(3): 128-136 doi:10.6046/zrzyyg.2023079

SONG Jiaxin, LI Yikun, YANG Shuwen, LI Xiaojun. NSCT-based change detection for high-resolution remote sensing images under the framework of change vector analysis in posterior probability space[J]. Remote Sensing for Land & Resources, 2024, 36(3): 128-136 doi:10.6046/zrzyyg.2023079

0 引言

遥感影像变化检测技术通过比较同一地理区域不同时相的遥感影像,以识别该地理区域发生的变化。由于其具备高时效和低成本的优点,遥感影像变化检测技术近年来应用非常广泛,在土地利用和土地覆盖变化检测、森林植被变化检测、灾害治理和国土资源调查方面发挥了极大作用^{[1⇓⇓⇓-5]}。

目前,多数变化向量分析(change vector analysis,CVA)法^[6⇓-8]基于变化强度图进行变化检测,如Li等 ^[9]提出了基于CVA和非下采样轮廓波变换(nonsubsampled-contourlet transform,NSCT)的遥感图像变化检测方法(CVA-NSCT),该方法通过NSCT分解CVA变化强度图,实现了不同方向和尺度的变化检测^[10⇓-12]。然而,不同地物在CVA变化强度图中的变化幅度往往不同,使用单一阈值时CVA-NSCT法难以保证检测精度。由此,Chen等 ^[13]通过研究分类后比较(post classification comparison,PCC)法和CVA的特点^[14-15],提出了一种后验概率空间变化向量分析(change vector analysis in posterior probability space,CVAPS)法。在CVAPS框架下,由于每个像素属于某种地物的后验概率在0~1范围内,使得基于后验概率变化向量计算得到的不同类型变化幅度在相同的范围内,使其更适于使用单一阈值进行变化检测。因此,与CVA-NSCT相比,将CVAPS与NSCT结合很可能提高变化检测的精度。

然而,原始的CVAPS法利用支持向量机(support vector machine,SVM)估计像素的后验概率向量,会受到图像中混合像元等问题的影响^[16]。因此,李轶鲲等^[17]提出了耦合模糊C均值聚类(fuzzy C-means, FCM)和简单贝叶斯网络(simple Bayesian network,SBN)的CVAPS方法(FCM-SBN-CVAPS)。该方法通过FCM算法以不同的隶属度将遥感影像中混合像元分解为不同的信号类,有效地缓解了混合像元问题。之后,通过SBN^[18]建立信号类和地物之间的多对多的随机链接,计算出遥感影像中各像元的后验概率空间向量。但是,该模型仅适用于中低分辨率的图像,在高分辨率遥感影像变化检测中无法取得较高的精度。

综上所述,本文将FCM-SBN-CVAPS与NSCT相结合,提出了一种新的变化检测方法FCM-SBN-CVAPS-NSCT。该方法用NSCT将FCM-SBN-CVAPS法生成的变化强度图在不同方向和尺度上,分解成噪声方差不同的高频子图。由于CVAPS变化强度图中不同类型的变化幅度范围一致,因此每个高频子图可以通过单一阈值较好地检测出某一尺度和方向的变化,从而提高最终变化检测的精度。本文提出的变化检测方法的贡献在于: 使用NSCT分解CVAPS变化强度图,在后验概率空间,实现多尺度和多方向的变化检测,以提高高分辨率遥感影像变化检测的精度。

1 FCM-SBN-CVAPS-NSCT模型构建

本文提出了一种多模式结合的变化检测算法FCM-SBN-CVAPS-NSCT,总体流程如图1所示。

图1

新窗口打开| 下载原图ZIP| 生成PPT

图1 FCM-SBN-CVAPS-NSCT算法流程

Fig.1 Change detection flow chart of FCM-SBN-CVAPS-NSCT

首先,通过FCM-SBN-CVAPS计算后验概率变化强度图; 其次,利用NSCT分解变化强度图,实现多尺度和多方向的高分辨率遥感影像变化检测。在NSCT分解过程中,先通过NSCT将后验概率图像分解到尺度s,生成1个低频信号图像和3组高频信号图像,每组有8个方向的高频子图。由于受输入图像噪声的干扰,要对高频图像中每个子带图像进行阈值处理,阈值 $T_{s, l}$ 的计算公式为:

(1)

T_{s, l} = C δ_{s, l}

式中: $δ_{s, l}$ 为s尺度上l子带上图像的噪声方差; C为平衡细节信息和噪声的参数。

鉴于NSCT分解后第三尺度的低频图像中信息丢失较为严重,因此使用FCM-SBN-CVAPS计算出的原始后验概率强度图作为重构的第三尺度低频成分,结合去噪后的高频子带图像,重构成新的变化强度图。最后,使用Otsu法^[19]计算变化阈值,生成变化二值图,经过形态学滤波,去除细小相斑和孔洞,生成最终的变化结果。

1.1 SBN模型

本文所使用的FCM-SBN-CVAPS方法模型为像素、信号和地物构成的3层SBN网络(图2)。

图2

新窗口打开| 下载原图ZIP| 生成PPT

图2 简单贝叶斯网络

Fig.2 Simple Bayesian network

SBN基于像素( $p_{i}$ )与信号类( $S_{k}$ )和信号类(S_k)与地物( $O_{v}$ )的随机链接 $P (S_{k} | p_{i})$ 及 $P (O_{v} | S_{k}),$ 计算得到像素属于不同地物的后验概率 $P (O_{v} | p_{i}),$ 其计算公式为:

(2)$P\left( {{O}_{v}}|{{p}_{i}} \right)=\overset{n}{\mathop{\underset{k=1}{\mathop \sum }\,}}\,P({{O}_{v}}|{{S}_{k}})P({{S}_{k}}|{{p}_{i}})$,

式中: $P (S_{k} | p_{i})$ 可以通过模糊隶属度 $u_{k} (i)$ 估计得到; $P (O_{v} | S_{k})$ 可以通过专家提供的训练样本估计得到。根据贝叶斯公式,可表示为:

(3)$P\left( {{O}_{v}}|{{p}_{i}} \right)=P\left( {{O}_{v}} \right)\overset{n}{\mathop{\underset{k=1}{\mathop \sum }\,}}\,\frac{P({{S}_{k}}|{{O}_{v}})P({{S}_{k}}|{{p}_{i}})}{P\left( {{S}_{k}} \right)}$,

式中: $P (S_{k} | p_{i})$ 为像素 $p_{i}$ 属于信号类 $S_{k}$ 的程度,可由FCM估计得到; $P (S_{k} | O_{v})$ 为信号类 $S_{k}$ 发生在地物 $O_{v}$ 中的概率; $P (O_{v})$ 为地物 $O_{v}$ 的先验概率; $P (S_{k})$ 为信号类 $S_{k}$ 的概率。由此可计算出像素 $p_{i}$ 属于不同地物的后验概率向量 $α = (α_{1}, \dots, α_{v}, \dots, α_{m})$ ,其中 $α_{v} = P (O_{v} | p_{i})$ 。

1.2 FCM算法

FCM聚类源自于普通C均值聚类,能够将单一像素以不同隶属度分解成不同的信号类,从而解决遥感图像中常存在的混合像元问题。因此,本文使用FCM计算模糊隶属度建立像元与信号类间的随机链接,FCM聚类采用误差平方和准则函数,公式为:

(4)$J\left( U,\Psi \right)=\overset{N\times M}{\mathop{\underset{i=1}{\mathop \sum }\,}}\,\overset{n}{\mathop{\underset{k=1}{\mathop \sum }\,}}\,{{({{u}_{k}}\left( i \right))}^{q}}{{\left\| {{p}_{i}}-{{\psi }_{k}} \right\|}^{2}}$,

式中: $p_{i}$ 为图像里的第i个像素; $u_{k} (i)$ 为像素 $p_{i}$ 属于信号类 $S_{k}$ 的隶属度; q为模糊度参数; $Ψ = {ψ_{1}, \dots, ψ_{k}, \dots, ψ_{n}}$ 为信号类中心的集合; $ψ_{k}$ 为信号类 $S_{k}$ 的中心; U是由 $u_{k}$ 构成的隶属度矩阵; N和M分别为影像的高度和宽度。由于SBN中的 $P (S_{k} | p_{i})$ 和FCM中的 $u_{k} (i)$ 有着同样的约束条件,因此本文使用 $u_{k} (i)$ 估算 $P (S_{k} | p_{i})$ 。

FCM需要预先确定聚类数(信号类数)和模糊度参数q。其中,聚类数需要能够反映遥感影像光谱特点。过少的聚类数无法精确反映影像光谱信息,而过多的聚类数则会造成光谱信息冗余并增加FCM和SBN的运算量,两者都会降低变化检测精度。此外,过低的q会导致像素关于信号类的隶属度差异较大,产生较大的累积聚类误差; 过高的q会导致像素关于不同信号类的隶属度差异过低,无法有效反映像素的光谱特点,并增大FCM的运算量,降低检测效率,两者都会降低变化检测精度。本文通过实验确定最佳的聚类数和q。

1.3 基于模糊聚类隶属度的SBN学习

为了计算像素 $p_{i}$ 属于地物的后验概率 $P (O_{v} | p_{i})$ ,需要基于专家提供的地物 $O_{v}$ 的训练像素集合 $T_{v}$ 计算地物 $O_{v}$ 中信号类 $S_{k}$ 的频率 $F_{v} (k)$ 以计算 $P (S_{k} | O_{v})$ 。 $F_{v} (k)$ 的计算公式为:

(5)

F_{v} (k) = \sum u_{k} (x) \forall p_{x} \in T_{v}

式中 $p_{x}$ 为训练集 $T_{v}$ 里第x个像素。

依据式(5)中得到的 $F_{v} (k)$ ,近似估算出 $P (S_{k} | O_{v})$ ,公式为:

(6)$P({{S}_{k}}|{{O}_{v}},{{T}_{v}})=\frac{{{F}_{v}}\left( k \right)}{\underset{l}{\mathop \sum }\,{{F}_{v}}\left( l \right)}$。

先验概率 $P (S_{k})$ 是计算后验概率 $P (O_{v} | p_{i})$ 的归一化因子,可使用全概率公式^[20]计算得到。

1.4 后验概率空间向量

假设 $p_{i}$ 表示图像中第i个像素。由地面覆盖类型O₁和O₂组成。在时间 $t_{1}$ 和 $t_{2}$ 的遥感图像中,像素的光谱略有变化,导致在时间 $t_{1}$ ,后验概率 $P (O_{1} | p_{i})$ 和 $P (O_{2} | p_{i})$ 分别为51%和49%,而在时间t₂,后验概率 $P (O_{1} | p_{i})$ 和 $P (O_{2} | p_{i})$ 分别为49%和51%。因此,PCC基于最大后验概率原则判定像素 $p_{i}$ 是已变化的像素,而不考虑其在两时相图像中仅有轻微光谱差异。

由此可见,PCC对变化像元的过高估计,由像元的累计分类错误造成的。因此,为了避免PCC中累积分类误差对变化检测的不利影响^[12],Chen等^[13]提出了CVAPS。假设像素在 $t_{1}$ 时相的后验概率向量为 $α^{1} = (α_{1}^{1}, \dots, α_{v}^{1}, \dots, α_{m}^{1})$ ,在 $t_{2}$ 时相的后验概率向量为 $α^{2} = (α_{1}^{2}, \dots, α_{v}^{2}, \dots, α_{m}^{2})$ , $α_{v}^{1}$ 和 $α_{v}^{2}$ 分别为t₁和t₂时相的后验概率 $P (O_{v} | p_{i})$ ,m为遥感影像中地物的数量,则后验概率的变化向量为:

(7)

Δ α = α^{2} - α^{1}

。

此时,后验概率空间向量变化幅度 $(Δ α)$ 采用欧氏距离可表示为:

(8)$\left\| \Delta \alpha \right\|=\sqrt{\overset{m}{\mathop{\underset{v=1}{\mathop \sum }\,}}\,{{(\alpha _{v}^{2}-\alpha _{v}^{1})}^{2}}}$。

根据式(8),后验概率向量的微小变化,不会造成变化幅度 $(Δ α)$ 的剧烈变化从而缓解了PCC对变化像元的过度估计问题。

1.5 NSCT

NSCT是一种多尺度、多方向的图像变换方法,它可以将图像分解成不同频率和方向的子带,提取出图像的局部特征,具有较好的局部特征表示和较好的多分辨率特性。

NSCT由正交对称的小波变换(undecimated wavelet transform,UWT)和轮廓波变换(contourlet transform,CT)2个步骤组成。UWT可以将图像分解成不同频率的子带。CT可以将图像分解成不同尺度和不同方向的子带,提取出图像中的局部特征。为了避免下采样的影响,NSCT使用了非下采样的方式进行CT,得到了更好的多尺度和多方向的特征表示。由于没有上采样和下采样的过程,所有子带的大小与原始图像相同,因而NSCT具有平移不变性的特点。

如图3所示,各层子带分解方向数为 $2^{n}$ ,图中所示为采用三层分解的NSCT结构在第一尺度上有8个方向的子带,在第二尺度上有4个方向的子带。信号完全重建的条件是非下采样金字塔滤波器组和非下采样方向滤波器组满足:

(9)

H_{0} (z) G_{0} (z) + H_{1} (z) G_{1} (z) = 1

式中: $G_{0} (z)$ 和 $G_{1} (z)$ 分别为低通和高通重构滤波器; $H_{0} (z)$ 和 $H_{1} (z)$ 分别为低通和高通分解滤波器。

图3

新窗口打开| 下载原图ZIP| 生成PPT

图3 NSCT三阶段分解结构

Fig.3 NSCT three-stage decomposition structure

1.6 变化强度图后处理

本文基于FCM-SBN-CVAPS-NSCT生成的变化强度图,采用Otsu阈值算法生成初始变化二值图。然而,由于影像中的噪音影响,初始变化二值图中有大量的细小相斑和孔洞,影响检测精度。本文在变化二值图中删除小于一定面积的相斑,并通过形态学膨胀填充孔洞。然而,经过形态学膨胀,有些孔洞依然无法被填充。因此,本文最终将小于一定面积的孔洞直接删除。相斑和孔洞的删除面积,及形态学膨胀的结构体尺寸均由实验确定。为了保证对比实验的客观性,本文所有对比算法均采用同样的后处理流程。

2 实验结果及分析

本文所使用的3组不同研究区影像来自商汤科技主办的AI瞰世界·2020人工智能遥感解译大赛所使用的SenseEarth平台数据集(https://rs.sensetime.com/),分辨率为3 m,影像尺寸为512像素×512像素。研究区1实验影像共包含4类地物,分别为: 建筑、荒地、林地和水体。研究区2实验影像共包括5类地物,分别为: 空地类型I,空地类型II,树木,草地,建筑物。研究区3实验影像共包括5类地物: 空地,道路,建筑物类型I,建筑物类型II和建筑物类型III。实验所用遥感影像经过了辐射校正、几何纠正和图像配准等预处理操作。算法运行环境为 Intel(R)Core(TM)i7-1165G7 2.80 GHz,24 G内存,4核8线程处理器。为了验证方法的有效性,本文实现了CVA-NSCT,CVA-NSCT-FCM,SVM-CVAPS和FCM-SBN-CVAPS算法,同本文提出的FCM-SBN-CVAPS-NSCT算法进行比较。

2.1 检测结果示例

经实验确定,FCM-SBN-CVAPS和FCM-SBN-CVAPS-NSCT方法的运行参数为: 模糊度参数q=3.5,聚类中心数为30,每类地物训练样本数为1 000。如图所示,图4—6分别为研究区1—3的遥感影像及其变化检测结果。其中黄色区域为人工检测的变化区域,红色标记为不同算法检测的变化区域,绿色标记为不同算法的错检区域,浅蓝色标记为不同算法漏检的区域。

图4

新窗口打开| 下载原图ZIP| 生成PPT

图4 研究区1变化检测算法比较

Fig.4 Comparison of change detection algorithms of study area 1

图5

新窗口打开| 下载原图ZIP| 生成PPT

图5 研究区2变化检测算法比较

Fig.5 Comparison of change detection algorithms of study area 2

图6

新窗口打开| 下载原图ZIP| 生成PPT

图6 研究区3变化检测算法比较

Fig.6 Comparison of change detection algorithms of study area 3

针对研究区1,从图4(d)—(e)中可以看出,CVA-NSCT没有明显的错检区域,而CVA-NSCT和CVA-NSCT-FCM两种算法的漏检区域面积均较大; 从图4(f)中可以看出相对其他算法,SVM-CVAPS算法错检区域范围明显更大。从图4(g)—(h)中可以发现,FCM-SBN-CVAPS-NSCT算法比FCM-SBN-CVAPS算法的错检范围小,且较好地保持了变化区域的边缘信息。

针对研究区2,从图5(d)—(f)中可以看出,CVA-NSCT,CVA-NSCT-FCM和SVM-CVAPS这3种算法有较大面积的漏检区域; 图5(e)—(h)中可以发现,FCM-SBN-CVAPS-NSCT算法的错检测区域和漏检区域相对较小。

针对研究区3,从图6(d)—(f)中可以看出,CVA-NSCT,CVA-NSCT-FCM和SVM-CVAPS这3种算法存在较多的错检区域; 从图6(g)—(h)中可以看出,FCM-SBN-CVAPS和FCM-SBN-CVAPS-NSCT算法只存在少量漏检区域,且FCM-SBN-CVAPS-NSCT的漏检区域更小。

2.2 实验参数影响

本节实验针对研究区1每种地物分别选择1 000,2 000,3 000,4 000和5 000个训练像素,检测训练样本对算法性能的影响。实验结果如图7所示。

图7

新窗口打开| 下载原图ZIP| 生成PPT

图7 不同数量训练像素数对算法Kappa系数的影响

Fig.7 Effect of different number of training pixels on the Kappa coefficient of the algorithm

在聚类数为10条件下,随训练像素数增加,本文算法Kappa系数最小为0.823 2,最大为0.845 9,差异为0.022 6。在聚类数为30条件下,随着训练像素数增加,本文算法Kappa系数最小为0.837 9,最大为0.845 4,差异为0.007 5,曲线变化最为稳定。在聚类数为50条件下,随着训练像素数增加,本文算法Kappa系数最小为0.834 3,最大为0.843 8,差异为0.009 5。另外,当训练像素数为1 000时,不同聚类数下本文算法Kappa系数分别为: 0.843 6(聚类数10),0.845 4(聚类数30)和0.839 4(聚类数50)。因此,当聚类数达到一定数量,基于少量的训练像素数,本文算法就可以获得较好的结果,并且大大降低运行的时间,提高了效率。然而,由于不同遥感影像的光谱特点有差异,本实验中小样本条件下的最佳聚类数(30)不一定适用于其他影像。因此,本课题组下一步计划设计能够根据影像光谱特点自动确定合适聚类数的FCM算法。另外,在所有训练样本数下,本文算法均明显高于SVM-CVAPS和FCM-SBN-CVAPS。总之,本文算法在不同聚类数和训练像素数下,精度较为稳定。因此,用户容易确定合适的聚类数,并且只需要耗费较少的工作量选择训练样本,即可取得满意的变化检测结果。

为了测试q和聚类数对本文算法的影响,q取6个值,即1.5,2,2.5,3,3.5,4,聚类数选取10,30和50进行实验。从图8可以看出,随模糊参数的增大,FCM-SBN-CVAPS-NSCT方法的Kappa系数随之上升。在q<2时,本文算法Kappa<0.75,在q>2.5时,3种聚类数下,本文算法Kappa均超过了0.8。本文算法Kappa值在q=3.5时最高,在q>3.5时下降。这是由于过低和过高的模糊度参数会降低混合像元分解的精度,从而降低之后的变化检测性能。另外,如图8所示,聚类数也对本文算法性能有一定影响。聚类数为30时,Kappa系数值在不同模糊度参数q下均较高。在不同模糊度参数下,FCM-SBN-CVAPS-NSCT方法的Kappa系数均高于FCM-SBN-CVAPS方法。与FCM-SBN-CVAPS相比,本文算法在不同q下,性能更为稳定。因此,用户较容易选择合适的q,以取得满意的变化检测结果。

图8

新窗口打开| 下载原图ZIP| 生成PPT

图8 不同聚类数和模糊参数q对Kappa系数的影响

Fig.8 Effect of different number of clusters and fuzzy parameter q on Kappa coefficient

NSCT对图像的尺度分解后得到的高频图像中含有大量的细节信息和噪声,而保留细节的同时有效去除噪声会对重构结果和变化检测有重大影响,因此首先根据式(1)在每个子带中独立地进行去噪处理。从图9可以看出,式(1)中的阈值系数C对检测精度有着重要的影响。在本文算法中系数C=10时效果达到最好,低于10时的Kappa系数均在0.77以上。系数C超过10后,过度消除了高频子带中有用的细节信息,导致Kappa系数快速下降。另外,从图9中可以看出,阈值系数C对CVA-NSCT法的Kappa值有类似的影响,说明有效保留变化检测细节的同时消除噪声对提高变化检测精度有重要意义。最后,从图中还可以看出,由于后验概率变化强度图中的不同类型变化强度在同一范围内,在不同的系数C下,本文算法的Kappa系数始终高于CVA-NSCT。因此,以上实验结果证明, 将NSCT与FCM-SBN-CVAPS相结合可以有效改善每个子带图像的去噪效果,从而提高变化检精度。

图9

新窗口打开| 下载原图ZIP| 生成PPT

图9 NSCT去噪阈值参数对Kappa系数的影响

Fig.9 Effect of NSCT denoising threshold parameters on Kappa coefficients

2.3 NSCT分解层数影响

通过分析表1可以得知,对于研究区1, NSCT的分解层数为4时,变化检测结果的Kappa系数为0.857 1,比分解层数为3时高出了0.011 7,而对于研究区2和研究区3,分解层数为3时,变化检测结果的Kappa比分解层数为4时分别高出了0.050 4和0.043 8。由此可见,不同研究区遥感影像的光谱纹理特点不同,其最佳分解层数也不同。因此NSCT的最佳分解层数需要通过相关实验确定。就本文的3个研究区而言,NSCT的分解层数为3时,本文算法均取得较好的检测精度,因此将NSCT分解层数统一设置为3。

表1 NSCT分解层数对Kappa系数的影响

Tab.1 The influence of the number of NSCT decomposition layers on the Kappa coefficients

研究区	分解层数
研究区	1	2	3	4	5
研究区1	0.758 6	0.827 1	0.845 4	0.857 1	0.823 1
研究区2	0.708 2	0.714 2	0.762 2	0.711 8	0.737 6
研究区3	0.562 4	0.706 9	0.889 5	0.845 7	0.839 3

新窗口打开| 下载CSV

2.4 算法综合性能比较

从研究区1—3的实验结果(表2—4)可知,本文算法虽然没有取得最低的错检率和漏检率,但总体精度和Kappa系数都明显高于其他算法。以表2为例,本文算法的总体精度和Kappa系数比FCM-SBN-CVAPS分别高出2.21百分点和0.100 9,证明了本文算法的有效性。同时,FCM-SBN-CVAPS-NSCT算法的性能也优于其他算法,其总体精度和Kappa系数比SVM-CVAPS分别高出5.67百分点和0.223 8,比 CVA-NSCT分别高出1.47百分点和0.139 2,比CVA-NSCT-FCM分别高出1.7百分点和0.118 0。本文算法性能优于CVA-NSCT和CVA-NSCT-FCM的原因是由于后验概率变化强度图中的不同类型变化强度在同一范围内,有效提高了单一阈值在变化检测中的效果。

表2 研究区1变化检测算法性能比较

Tab.2 Performance evaluation of change detection algorithms of study area 1

算法	参数	错检率/%	漏检率/%	总体精度/%	Kappa系数
FCM-SBN-CVAPS-NSCT	30聚类 q=3.5 1 000训练样本C=10	14.50	13.63	97.48	0.845 5
FCM-SBN-CVAPS	30聚类 q=3.5 1 000训练样本	32.13	10.92	95.27	0.744 6
SVM-CVAPS	cp=13 gamma=3 1 000训练样本	47.67	8.98	91.81	0.621 7
CVA-NSCT-FCM	30聚类 q=3.5	22.40	27.25	95.70	0.727 5
CVA-NSCT	C=10	6.48	40.43	96.01	0.706 3

新窗口打开| 下载CSV

表3 研究区2变化检测算法性能比较

Tab.3 Performance evaluation of change detection algorithms of study area 2

算法	参数	错检率/%	漏检率/%	总体精度/%	Kappa系数
FCM-SBN- CVAPS-NSCT	30聚类 q=3.5 1 000训练样本C=10	10.11	8.05	88.76	0.762 2
FCM-SBN-CVAPS	30聚类 q=3.5 1 000训练样本	13.98	3.73	86.31	0.705 6
SVM-CVAPS	cp=15 gamma=3 1 000训练样本	11.96	23.41	79.35	0.581 3
CVA-NSCT-FCM	30聚类 q=3.5	13.02	68.10	55.50	0.206 8
CVA-NSCT	C=10	5.47	68.59	56.95	0.238 9

新窗口打开| 下载CSV

表4 研究区3变化检测算法性能比较

Tab.4 Performance evaluation of change detection algorithms of study area 3

算法	参数	错检率/%	漏检率/%	总体精度/%	Kappa系数
FCM-SBN-CVAPS-NSCT	30聚类 q=3.5 1 000训练样本C=10	5.51	14.13	98.15	0.889 6
FCM-SBN-CVAPS	30聚类 q=3.5 1 000训练样本	0.87	27.66	97.27	0.822 2
SVM-CVAPS	cp=13 gamma=3 1 000训练样本	19.46	20.10	96.17	0.780 0
CVA-NSCT-FCM	30聚类 q=3.5	24.95	15.12	95.81	0.773 4
CVA-NSCT	C=10	13.31	15.85	97.22	0.838 6

新窗口打开| 下载CSV

3 结论与讨论

论文提出了一种有监督的遥感影像变化检测算法FCM-SBN-CVAPS-NSCT。实验结果表明,由于后验概率变化强度图中不同类型的变化幅度范围一致,使本文算法在单一阈值条件下,比CVA-NSCT精度更高。

另外,由于引入了NSCT变换,本文算法实现了后验概率空间下的多方向和多尺度的变化检测。与FCM-SBN-CVAPS算法相比,本文算法精度更高且不易受参数影响,具有更高鲁棒性。

虽然,所提方法取得了令人满意的效果,但它需要通过实验确定SBN的模糊度q、信号类数量和NSCT尺度分解层数。因此,在未来的工作中,我们计划通过训练影像,自动确定最优模糊度q、信号类别数和NSCT尺度分解层数,以增强所提出算法的实用性。此外,由于NSCT变化在分解和重构过程耗时较长,降低了本文算法的效率。因此,如何提高算法效率,是下一步研究的目标。

参考文献

原文顺序

文献年度倒序

文中引用次数倒序

被引期刊影响因子

[1]

李德仁

利用遥感影像进行变化检测

[J]. 武汉大学学报(信息科学版), 2003, 28(s1):7-12.