SAR图像变化检测综述

SAR变化检测研究综述1闻彩焕1,2,3，刘君臣1,3，蔡红奎41中国矿业大学环境与测绘学院，江苏徐州 (221116)2中国测绘科学研究院摄影测量与遥感所，北京 (100830)3江苏省资源环境信息工程重点实验室，江苏徐州（221116）4河南理工大学土木工程学院，河南焦作 (454000)E-mail:wencah@摘 要：SAR由于其具有的传统光学遥感不可比拟的全天候、全天时的优势而逐渐成为变化检测的重要数据源，在国土资源监测、灾害评估与监测等各个方面起到日益重要的作用。

本文介绍了SAR变化检测的主要研究内容，常用的SAR变化检测算法，并在广泛阅读主要遥感刊物上发表的文献的基础上总结了国内外的研究现状，最后对SAR变化检测研究的前景做了展望。

关键词：遥感，SAR，变换检测中图分类号：TP791 引言遥感变化检测是指通过对不同时期同一区域的遥感图像进行比较分析，根据图像之间的差异得到我们所需要的地物或目标的变化信息。

现代遥感技术的飞速发展为变化检测提供了一种便捷的途径，遥感数据成为变化检测的主要数据源。

与可见光和红外遥感相比，微波遥感具有无可比拟的优点：微波能穿透云雾、雨雪，具有全天候、全天时的工作能力。

二，微波对地物有一定穿透能力。

三，采用侧视方式成像，覆盖面积大。

正是这些优点，使得SAR图像日益成为变化检测的重要数据源。

SAR变化检测技术的需求日益广泛。

目前，全球坏境变化加剧，城市急速发展，洪水、地震等自然灾害时有发生，这些都需要及时掌握相关动态信息，为相关决策部门提供支持，而SAR的种种优点为快速响应提供了技术支持和应急保障。

2 变化检测的研究内容2.1图像配准图像的配准精度对变化检测的精度影响很大，如果图像没有较高精度的配准，则图像的错位会带来大量的虚警。

因此，针对这一方面，Dai和Khorran[1]得出的结论为：好的变化检测结果需要高精度的配准算法。

Jeannie Moulton [2]等针对多时相SAR图像的平移及偏转误差提出一种解决方法。

意大利的Paolo [3]针对高分辨SAR影像配准的困难，提出基于象元与特征相结合的方法进行未精确配准的变化检测研究。

2.2噪声抑制基于SAR图像的变化检测对斑点的抑制要求很高：降班算法要既能有效降斑，还能保持图像的细节信息。

因此，需要采用性能良好的自适应降斑方法来完成。

相干斑抑制方法可以分为两类：一类是成像前的多视平滑处理，但同时也降低了图像的空间分辨率，所以目前常用的方法为成像后的去相干斑噪声的滤波技术。

滤波方法主要有：（1）传统方法，如均1项目来源：973课题“高性能的遥感影像网格计算与信息解译的智能方法”（编号：2006CB701303）子课题“SAR影像目标提取与信息解译的智能化处理”。

值滤波、中值滤波等；（2）模型方法，假定静态的噪声模型，采用相应的滤波器进行处理，如Kalman滤波和Lee滤波；（3）几何滤波方法，如Gamma MAP滤波；（4）基于局域统计自适应滤波，这种方法能在平滑噪声的同时较好的保持边缘信息，是目前的SAR变化检测常用的方法。

2.3变化信息提取及后处理进行变化信息提取时，可以选择适当的变化检测方法对两个时相的遥感影像进行变化信息的提取。

根据实际野外观测提供目标区变化的真实信息，并借助相关统计资料和专题资料对变化信息进行后处理，归并小图斑，消除提取出的变化信息中的伪变化信息[4]。

2.4精度评价通常采取误差矩阵法进行精度评价。

根据真实变化图与变化信息图比较，归纳到误差矩阵中。

变化误差矩阵如表1。

表1误差矩阵变化检测结果实际的变化像元Ci 实际的非变化像元Cj 检测出的变化像元 C0 Cid Cjd检测出的非变化像元C1 Cir Cjr根据这一变化误差矩阵，变化检测的性能可以通过下面的参数来定量分析：虚警率Pjd为实际未变化但检测为发生变化的比例，表示了检测出的变化点中误判点的比例。

Pjd= Cjd/ C0 (1) 漏检率Pir为实际发生变化但检测为未变化的比例，表示了未检测出的变化点在所有变化点中的比例。

Pir= Cir/ C1 (2) 检测概率Pid为实际发生变化或非变化且检测为变化或非变化的比例，表示了对所有变化区域的检测百分比。

Pid= （Cid+ Cir ）/ （C0+ C1） (3) 3 常用的变化检测算法SAR图像变化检测是研究不同时期的同一场景图像之间发生的变化。

归结起来，近些年逐步形成的SAR图像变化检测方法可以分为：(1)基于简单代数运算的变化检测，经典的变化检测方法包括图像差值法、图像比值法、对数比值法；(2)基于图像变换的变换检测，经典的变化检测方法包括主成分分析、变化向量分析法、相关分析法图像变换法；(3)基于图像分类的变化检测方法。

3.1 代数运算法早期变化检测的主要方法是利用图像差值法，主要利用两幅图像对应像素点的灰度差值作为变化判别依据，该方法简单、直接，但是容易受成像质量、噪声等的影响。

图像差值法应用于海岸线环境、森林变化、沙漠化等。

与差值法相比，比值差异图像与原始SAR图像的平均强度值无关，对于SAR图像的乘性噪声不敏感，而且减少了辐射定标误差的影响。

图像比值法对于SAR图像上的乘性噪声是不敏感的，被广泛用于城区变化检测。

对数比值法，在得到对应像素的比值后，再取其对数，能压缩图像的变化范围，将乘性噪声转换为加性噪声。

代数运算法的优点是直接、简单，主要缺点在于(1)没有考虑多时相图像之间的相关性，简单计算可能出现虚警率；(2)变化阈值确定困难。

3.2图像变换法Lillesand[5]提出的主分量法，又称K-L变换，是一种经典的数学变换方法，把原来多个波段中的有用信息集中到互不相关的新成分图像中，达到冗余压缩和信息集中的目的。

在进行主成分变换时，利用协方差矩阵和相关矩阵得到的主分量是不同的，由相关矩阵推导的主分量变换对于多时相分析是尤其有用的[6][7]。

Malila提出了变化向量分析法(Change Vector Analysis,CVA)[6],Johnson[8]等详细阐述了这种方法的具体内容、优缺点。

变化向量法是一种多变量的方法，描述从第一时间到第二时间的光谱变化的强度和方向。

如果变化向量的幅值超过给定的门限，则判定该像素发生变化，变化向量的方向包含变化类型信息。

该方法可用在多通道极化SAR图像的变化检测，或者用于多特征，如空间结构特征，纹理特征等分析。

相关分析法，斑点的时相去相关提供了地表在结构或介电特性上可能的变化信息，并且这种变化检测不依赖于定标精度。

实际应用中，通常使用一般的强度图像来计算斑点相关系数，不考虑相位信息。

3.3分类法分类后比较法，在70年代末开始应用于Landsat卫星影像的变化检测。

分类后比较法是对每幅图像单独进行分类，然后对图像的分类结果图进行比较，以确定变化的类别和区域。

如果对应像素的分类类别相同，则认为该像素没有发生变化，否则认为该像素发生了变化。

使用分类法的优点主要有：(1)可以给出区域变化的类型、位置、数量等有关地物性质信息；(2)受图像配准的影响要小。

此方法的主要缺点是它受分类器的影响很大，不同的分类器能够决定变化检测的精度高低。

同时分类法，将多时相图像构成一副复合影像进行分类，得到的分类结果图中的每一类代表一种变化类型。

此方法能够减少分类时间，但是难于标记变化类别，而且分类更加复杂，对训练样本的选取相当高。

随着计算机技术的发展，解决复杂问题的各种数学方法与专家知识被逐渐加入。

人工神经网络、向量机、专家智能等方法成为热点。

4 国内外研究现状近几十年来，SAR变化检测技术被越来越多的国内外学者所关注。

国外，SAR变化检测技术方兴未艾，国内，该技术正处于起步阶段[9]。

4.1 国外现状国外，Vilasenor [10]利用两幅重轨ERS-1 SAR影像对阿拉斯加北坡地区进行了变化检测，验证了雷达后向散射强度的差值比较可以作为变化检测的一种方法。

Rignot[11]等人根据多视SAR强度影像服从Gamma分布的假设，从理论上证明了比值法更适合多时相SAR影像的变化检测。

Bazi[12] [13] [14][15]、Bruzzone[16]等在多时相SAR比值图像的阈值变化检测算法上展开一系列研究。

针对直接代数运算法没有考虑到象元的空间上下文信息，误检和漏检现象严重，近几年来出现了在此基础上的基于上下文空间关系的变化检测，主要从两个方面进行改进。

一方面，单尺度上利用差异图像的上下文关系进行进一步的优化得到变化图, Bruzzone[17]首先利用 Markov随机场对差值图像进行分析，并求得最优的变化检测结果。

文献[18][19][20][21][22]都是基于马尔科夫随机场进行上下文分析。

文献[23][24]利用神经网络进行上下文分析；另一方面，利用差异图像的小波多尺度信息进行优化得到变化图，如Bovolo[25],Inglada[26],Celik.T[27][28]利用多尺度小波信息进行变化信息的优化。

Fung [29],Gong [30]等先后将主成份变换应用于变化检测中。

Qiu[31]等人将分块主分量分析(Multi-Block PCA)法应用于图像变化检测。

Rignot[11]等人利用ERS-1的重轨影像进行了不同方法的变换检测，得出，对于单视影像或者是视数较小的多视影像，基于去相关的变换检测方法效果最好。

Corri[32]通过计算三幅多时相影像两两之间的相关系数进行变化检测，并提出一种对建筑物变化敏感的检测方法。

Hyung Sup [33]等利用多时相JERS-1SAR数据，进行了多种方法的变化检测，认为去相关法对于检测物体属性的变化是非常有用的方法。

Takashi[34]等用相关系数法以及标准偏差系数法进行了比较，发现，相关系数法提取的区域要小于标准差系数的区域，对于大面积区域，更适合用标准差系数法来进行变化检测。

Roger[35]等利用Seasat和SIR-B不同时相图像，用最小距离分类器得到了森林主要变化类型的大小和范围。

White [36]对基于模型与基于神经网络方法进行了比较，认为神经网络训练得当，则此方法优于基于模型的方法。

Fabio [37]等提取出后向强度、相关参数以及纹理参数特征，进行神经网络监督分类，进行变化检测。

Welismiller[38]等采用同时分类法进行了德克萨斯地区海岸带环境变化检测。

White[36]针对SAR图像的斑点噪声以及配准问题，提出基于图像分割的方法，并利用神经元网络进行对象识别，在此基础上进行变化检测。

Craves[39]对多时相ERS-SAR数据进行了分割，分割后的图像叠加起来进行检测。

日本L．Zhu [40]利用分水岭分割实现建筑物提取，在此基础上进行变化检测。