D-INSAR中参考DEM误差和地表形变对相位贡献的灵敏度研究.pdfVIP

《测绘通报》测绘科学前沿技术论坛论文集 1 D-INSAR中参考DEM误差 与地表形变对相位贡献的灵敏度研究 1.2 1 1 1 龙四春 ，李陶 ，刘经南 ，张诗玉 （1.武汉大学 卫星导航定位技术研究中心，湖北武汉，430079） （2.湖南科技大学 建筑与城乡规划学院，湖南湘潭，411201 ） 摘 要：从多参数相位贡献的 D-INSAR 基本原理出发,分析讨论了地形相位和地表形变相位对干涉相位的影响,详细导出了参考 DEM 误差和地表形变量对差分干涉相位贡献灵敏度的影响因子、推理步骤和计算公式，得出雷达差分干涉测量中地表形变对相位贡 献灵敏度要远远高于参考DEM 误差对相位贡献灵敏度的结论。最后通过Bam 地震区的实例验证了该结论的正确性。 关键词：差分干涉测量 DEM 误差 形变量 灵敏度 一、引言 卫星合成孔径雷达差分干涉测量(Differential Interferometric Synthetic Aperture Radar , D-InSAR)凭借其独特的技术优势，如全天时、全天候、广覆盖、高空间分辨率和高形 变敏感度, 在区域地表形变测量中已得到广泛应用[1-9]。其基本思想是将覆盖同一地区的两幅 SAR 影像进行配准并取样到相同空间, 对同名像素进行相位差分得到干涉相位图, 从干涉相位 图中扣除参考相位和大气相位等因素的影响 , 最终提取形变信息。 初始干涉图每一像元上的初始干涉相位是由多个参数分量贡献的总和，可用如下函数模型 表示： Φ Φ +Φ +Φ +Φ +Φ （1） ref topo defo atmo noise 其中Φ表示干涉对生成的干涉相位，Φ 为参考椭球面引起的相位，Φ 为地形相位， ref topo Φ 为两次成像期间因地表位移引起的形变相位，Φ 为两次成像时大气状态不一致引起的 defo atmo 相位延迟，Φ 为残余噪声相位。本文假设其他相位贡献全部已经消除，只讨论差分干涉测量 noise 中参考DEM 误差及地表形变的相位贡献灵敏度。 二、地形相位和地表形变相位对干涉相位的影响 采用雷达干涉测量方法进行地表形变监测时必须满足两个条件：一、两次观测期间的变化 不能太大，像元之间变化梯度必须小于某一数值；二、两次观测期间，一个像元内的所有地物 散射体的特性应具有相似性，特别是表面位置变化中误差应只是雷达波长的10%-20%的量级[1]。 SAR 干涉测量原理示意图1 中，1、2 两雷达天线从同一目标表面单元接收到的相位差ϕ为： 4π λϕ 1 ϕ δ ⇒δ （2） λ r r 4π 作者简介：龙四春，湖南科技大学讲师，武汉大学卫星导航定位技术研究中心博士生，现主要从事INSAR 和大地测 量方面的研究。 基金项目：国家自然科学基金40604002）；武汉大学地球空间环境与大地测量教育部重点实验室基金 （905152533-05-08） 2 《测绘通报》测绘科学前沿技术论坛论文集 式中λ为波长，r ，r +δ 为两天线距散射体的距离，δ 可以看成地表形变沿雷达视线方 1 1 r