空间大地测量(sar和insar).docVIP

关于InSAR和D-InSAR的数据处理 一、合成孔径雷达干涉技术（InSAR） 合成孔径雷达干涉技术出现于20世纪60年代末.它是SAR与射电天文学干涉测量技术结合的产物。当SAR扫过地面同一目标区域时,利用成像几何关系,通过成像、一些特殊的数据处理和几何转换,即可提取地表目标区域的高程信息和形变信息。由于InSAR 技术有效利用了SAR的回波相位信息,测高精度为米级甚至亚米级,而一般雷达立体测量方法只利用灰度信息来实现三维制图,测高精度仅能达到数十米,因此该技术迅速引起了地学界及相关领域科研工作者的极大兴趣,现已成为微波遥感领域的研究热点. 干涉合成孔径雷达利用多个接收天线观测得到的回波数据进行干涉处理，可以对地面的高程进行估计，对海流进行测高和测速，对地面运动目标进行检测和定位。接收天线相位中心之间的连线称为基线，按照基线和航向的夹角，人们将InSAR分为基线垂直于航向的切轨迹干涉和沿航向的顺轨迹干涉。切轨迹干涉可以快速提取地面的三维信息，顺轨迹干涉主要用于动目标检测和海洋水流与波形测量。 二、 InSAR 基本原理 InSAR 测量模式主要有两种:一种是双天线单轨(Single Pass)模式,主要用来生成数字高程模型,一般用于机载SAR;另一种是双轨(Two Pass) 模式,主要用于获取地表变形,一般用于星载SAR. 下面以重复轨道干涉测量为例,简要介绍InSAR 技术的基本原理(见图1). 假设卫星以一定的时间间隔和轨道偏离(通常为几十米到1km左右)重复对某一区域成像,并在两次飞行过程中处于不同的空间位置和,则空间干涉基线向量为B,长度为B;基线向量B 与水平方向的夹角为基线倾角。和至地面点P的斜距分别为R 和R+△R;将基线沿视线方向分解,得到平行于和垂直于视线向的分量、;H为到参考面的高度;从发射波长为的信号经目标点P反射后被接收,得到测量相位, (1) 同样,另一空间位置上测量到相位， （2）式中,arg{}和arg{}表示不同散射特性造成的随机相位. 假设两幅图中随机相位的贡献相同,则和关于目标P点的相位差 （3） 也称为干涉相位,可由经过配准的两幅SAR SLC图共扼相乘得到.根据图1中的几何关系并利用余弦定理可得: （4） （5） 由于且，则（6） (4)、(5) 两式即为In SAR 确定高程的原理性公式. 三、合成孔径雷达差分干涉测量(D-InSAR) D-InSAR技术是在主动式微波合成孔径雷达 SAR 相干成像基础上发展起来的,它以合成孔径雷达复数据提供的相位信息为信息源,可从包含目标区域地形和形变等信息的一幅或多幅干涉纹图中提取地面目标的微小形变信息。它是迄今为止独一无二的基于面观测的形变遥感监测手段，可补充已有的基于点观测的低空间分辨率大地测量技术如全球定位系统、甚长基线干涉和精密水准等。它具有高灵敏度、高空间分辨率、宽覆盖率、全天候、几乎不受云雨天气制约和空中遥感等突出的技术优势，在对地震形变、地表沉陷及火山活动等大范围地表变形的测量研究中迅速得到了广泛的应用。 四、 D-InSAR原理 根据去除地形相位采用的数据和处理方法的不同，差分干涉测量技术可以分为二轨差分、三轨差分和四轨差分。下面以三轨差分为例简要说明D-InSAR原理。在下图中，设、分别组成形变前后的两对主辅影像，、所获取的影像为地形对，、为形变对。在三角形中有： 为地表形变量： 因此，由差分相位解缠之后的结果可求出地表形变量。同三轨差分类似，二轨D-InSAR 首先利用一对跨越形变期的SAR图像进行干涉处理，得到包含形变信息的干涉相位；然后利用已有的观测区数字高程模型DEM和SAR成像参数反演干涉相位；最后将干涉处理得到的相位与利用DEM反演得到的干涉相位进行差分处理，从而得到差分相位，解缠之后即可得到地表形变信息。 五、合成孔径雷达干涉测量数据处理流程 虽然 InSAR 的成像几何原理较为简单，但由于 SAR 本身成像的特殊性，在干涉测量过程中噪声及地形的不连续性使得具体的实现算法比较复杂。下面描述了具体的处理算法。 ①SAR 干涉数据的准备 针对不同的干涉应用选取合适的 SAR 干涉相对。成功的进行 InSAR 干涉处理要求 SAR 图像对必须相干，因此选取图像对时必须考虑两个方面，一是临界基线距的限制，二是时间去相干的影响。将获得的原始 SAR 信号采用 SAR 成像算法做成像处理，生成两个地面二维单视复图像。其目的是实现距离和方位向的二维成像，并得出图中各点的相位信息。此外，成像过程中对得出的多普勒中心、多普勒带宽和各次成像斜距 R1 和 R2 的估值，也是成像的重要参数。 ②SAR 影像的过采样与预滤波 为了避免在形成干涉条纹时出现频谱卷绕，需要对两幅 SAR 影像过采样。在进行干涉测量时，经常会出现两幅