杨成生，张 勤，张双成，赵超英

(长安大学地质工程与测绘学院，西安 710054)

0 引言

InSAR已被证实为一种强大的空间对地观测技术，可以用来研究许多种地球物理现象[1－4]。然而，如同其他微波遥感一样，当雷达信号2次穿越大气层时，会因大气水汽含量在时间和空间上的不确定性而产生不同程度的传播延迟，从而影响了差分InSAR的监测精度[5－8]。因此，大气延迟成为影响差分InSAR高精度获取缓慢微小地表形变的重要因素。作为一种新型的大气探测手段，20世纪90年代发展起来的地基GPS可以实时、连续、全天候、高精度地获取大气可降水汽值，且精度可以达到1～2 mm[9－10]。随着越来越多的区域连续运行GPS网的建立，为利用GPS水汽观测值建立区域水汽延迟改正模型，从而对SAR影像中的大气延迟进行改正提供了可能性。但是，GPS在区域上是以离散点形式分布的，即使相距最近的区域性GPS连续跟踪站间距一般也有几km到几十km不等，远远大于SAR影像干涉像元的实际大小。因此，在利用GPS观测数据改正InSAR影像中的大气延迟时，必须要对GPS水汽观测值进行空间插值，以建立与干涉图像元的空间分辨率相当的水汽延迟改正模型。以往的研究主要是用逆距离权(inverse distance weighted，IDW)倒数和克里格(Kriging)等插值方法直接对GPS水汽值进行空间插值[11]，但这些方法并未考虑到地形对大气延迟的影响，效果并不理想。

本文提出了一种改进的克里格插值(IKriging)方法。该方法在分别考虑大气在高程和距离上变异特性的基础上，按照双因素对插值结果的影响，建立一个综合的样本变异模型，并根据该变异模型建立已知数据与待插值点之间的协方差矩阵，求解已知数据对待插点的贡献权值，进而确定待插值点值。……