基于小波基线误差估计相干目标InSAR方法研究及应用.pdfVIP

中国地球科学联合学术年会 2014 ·1535 · 基于小波基线误差估计相干目标 InSAR 方法研究及应用 杨怀宁 1* 闫世勇2 杨 凌 1 1.中国地震应急搜救中心 北京 100049 ； 2. 中国矿业大学 徐州 221116 摘要：本文介绍了基线误差估计的常用方法，详细说明了小波基线误差估计的原理。以青海玉树地 震作为研究区，选用 ALOS1 PALSAR 数据，分离出基线误差，并得到研究区的时序形变特征。说明小 波基线误差估计方法的有效性。 一、引言 基线是描述 SAR 传感器平台间的空间位置，根据干涉的几何关系，基线是 InSAR 形变计算中一个 非常重要的参数。对于星载的 SAR 系统，可以利用卫星轨道的星历或状态向量计算基线，主要是基于 多普勒方程、距离方程、椭球方程的三维基线估计[1] 。由于 SAR 卫星提供的轨道数据精度有限，基线 参数的误差就会传播到干涉产品中。目前，研究人员提出了多个算法致力于减少由基线估计引起的误差， 主要可以分为两类:基线的精确估计和干涉产品校正。 基线的精确估计方法不但可以提高基线成分的精度，还可以分离出解缠相位中的轨道误差。Small et al.[2]、Joughin et al.[3]提出类似的算法，选取一定数量已知高程和解缠相位的 GCP，基于最小二乘求解。 Kohlhase et al.[4]建议通过调整卫星轨道迹线达到改正轨道误差的目的。该方法首先较好地估计轨道迹 线，粗略计算干涉图中的平行条纹数，从而得到最初迹线在径向和跨轨道方向的改正量。该方法适用于 相关性较好的干涉图，从而轨道误差条纹可数。刘广[5]提出采用残余平地相位校正垂直基线误差的方法， 并以 New Orleans 地区为例进行永久散射体方法处理。 [6] [7] 干涉产品校正方法主要包括基于干涉图条纹频率 ，或曲面拟合 等。基于干涉图自身信息的基线 估计主要是利用干涉条纹图中蕴含的信息如干涉条纹、干涉相位差等来估算基线分量。GAMMA 公司 采用傅里叶变换的方法获得垂直基线，并且还采用地面控制点的方法对基线进行了估计。Singh et al.[8] [9] 提出利用从 SAR 干涉测量数据中提取的干涉相干条纹精确估测基线的方法，廖静娟等 对该方法进行 了试验验证。另一个广泛应用的方法是通过估计方位向和距离向的相位斜坡来解决轨道误差。 本文详细说明了小波基线误差估计的原理。以青海玉树地震作为研究区，选用 ALOS1 PALSAR 数 据，分离出基线误差，并得到研究区的时序形变特征。 二、小波基线误差估计算法 小波基线误差估计方法属于干涉产品校正方法中的曲面拟合。该方法基于小波多尺度分解和稳健回 归，实现重轨差分干涉中轨道残差的校正[10] 。小波多尺度分解将信号分解成短、长波成分，其中长波 成分包括有轨道残差信息。而且，稳健回归可以有效地估计出轨道误差。 1 小波多尺度分解 二进制小波变换将原数据分解成 1 个低频成分（数据的平缓变化和趋势部分）和 3 个高频成分（数 据的较快变化和局部特征部分）,各成分占总数据量的 1/ 2。低频成分再经过小波变换又可以生成新的高 频和低频成分。通过小波多尺度分解,不仅可以从具有细部特征的空间数据集派生较概略程度的任意尺 度的空间数据集, 而且可以完成其逆过程, 实现完备的重构。 对解缠相位图进行小波分解，这种分解的目的是区分出解缠相位的不同波长成分。轨道误差的波 长相对较长，会影响到整幅干涉图，然而 DEM 误差和失相关噪声的波长相对于轨道误差的要短很多。 地表形变和大气延迟则表现出近似连续的相位梯度变化，与轨道误差大体类似。为解决这个问题，所研 究干涉图的区域应大于受影响的地表形变和大气延迟的区域尺度。 2 稳健回归 对解缠相位进行小波分解，去除短波成分，剩余相位主要