参照失效物理分析改进永久散射体提取过程.docVIP

参照失效物理分析改进永久散射体提取过程 关泽群 白俊武 （武汉大学遥感信息工程学院 ，武汉 430079） 摘 要：在雷达干涉测量中，永久散射体(PS)是一类不受时间、空间基线和大气延迟影响的目标，如何从SAR图像中自动识别出有效的永久散射体是PS干涉系统中关键的环节之一。本文从不确定性量化理论的角度考虑失效物理分析问题，利用失效物理分析方法研究PS影像特征与环境的关系，对多幅SAR图像进行独立分量分析（ICA）和小波包分解，然后考虑永久散射体的高信噪比特性，依据图像能量最大的原则获取分离后的稳定成分，并进行分割处理得到PS点的分布图。实验表明，该方法能够有效识别出可靠的PS点。 关键词：雷达干涉测量；永久散射体；失效分析； 独立分量分析； 小波包分解 中图分类号：P237 文献标识码：A 一、引言 本文研究了永久散射体（Permanent Scatterers，简写为PS）识别过程中的不确定性量化问题。由于在大部分的情形下，永久散射体的范围都很小，而且受到的干扰因素也很多，对其进行识别与弱目标检测是类似的。特定战场条件下的许多重要军事目标可以是永久散射体，诸如工业企业、城市小区、道路等人工环境下的许多物体也可以是永久散射体，因而不确定性量化问题研究，有助于这些永久散射体在信息比较弱的情形下的检测。 永久散射体与合成孔径雷达差分干涉测量技术关系密切。后者是一种具有低成本、大覆盖、高精度三大明显优势的地表形变监测手段[1]。但是在重复轨道雷达差分干涉的应用中，形变测量的精度和可靠性受到了相位观测误差、轨道误差、地形数据误差和解缠误差等的影响。其中相位观测误差主要受到时空失相关和大气失相关等因素的影响。为了解决这些问题，意大利Ferretti等人提出了永久散射体技术。这种技术通过同一地区获取的多幅多时相SAR图像，提取出经过长时间间隔仍具有较好相干性的（也就是相位稳定的）像元，这些像元不受时间、空间基线和大气延迟的影响，我们把这类像元称为永久散射体。可以利用这些永久散射体的插值拟合曲面，计算出DEM误差、视线方向的目标偏移值和大气效应贡献值，达到估计并去除大气效应相位贡献值、提高变形监测精度的目的[2]。如何从SAR图像中自动识别有效的永久散射体是PS干涉系统中的关键环节，尚有待深入研究。 目前，常用的PS探测方法主要有三种，即相干系数阈值法、振幅阈值法和相位离差阈值法[2, 3]。然而采用相干系数法探测的PS受邻域点信号的影响，会降低PS判别的准确性，相干系数的窗口大小的选取也影响探测结果的质量；而振幅离差指数阈值法只考虑了PS散射的稳定性，忽略了PS回波信号高信噪比的特性，也会产生误判。为改进PS识别的准确与合理性，罗小军等提出了探测永久散射体的振幅信息双阈值法，首先依据振幅阈值选取具有强反射特性点作为PS候选点，然后根据振幅离差阈值，选取散射特性稳定的候选点作为永久散射体[4]；陈强等提出了时序相关系数阈值法和振幅离差指数阈值法相结合的“双重阈值”探测法，即先采用较低的相关系数阈值剔除那些失相关比较严重的目标，为后续PS探测缩小搜索范围，然后基于时间序列振幅影像，采用振幅离差指数阈值法，进一步识别出可靠的PS点[5]。这些方法已在城区地面沉降、滑坡监测等研究中得到了应用，然而在某些情况下，特别是当SAR影像数量比较有限时，还不能保证PS目标的有效性。 本文拟从失效物理分析的角度，从数理统计和物理分析两个方面，对PS散射的稳定性进行考察，探索某些能容许较大的不确定性而其相干性不失效的PS的存在性[6,7]。基于这种想法，参照Ferretti等人的研究，本文提出了一种基于独立分量分析[8]和小波包分解[9]的PS点识别方法，该方法首先对同一地区的时间序列影像进行独立分量分析（Independent Component Analysis，简写为ICA），去除在较长时间内的各种干扰因素，然后根据图像能量最大的原则得到PS点的分布图，实验结果表明，这种方法能够更合理有效的提取永久散射体。 二、利用失效物理分析方法研究PS影像特征与环境的关系 建筑物、桥梁、工业用地中的集装箱、油罐等地物自身后向散射特性稳定，在没有发生明显变化的情况下，回波强度几乎不随时间变化，具有稳定的时间特性，是比较典型的几种永久散射体。一般来说它们的后向散射强度较强，在SAR图像上较亮，理想情况下，容易辨识。然而实际应用中，可以发现，受到环境因素的影响，它们的识别与提取变得复杂与困难。 通常在人工目标的四周都有背景目标，例如地表的草地、农田及山地等，这些背景目标属于分布目标类。表1给出了最常见的几种地表背景目标在SAR图像上的表现。当永久散射体处于不同的地表环境时，其在SAR影像上的表现不同，从而导致其辨识度也各不相同。根据永久