卢丽君，张继贤，王 腾

1.中国测绘科学研究院，北京100830；2.武汉大学测绘遥感信息工程国家重点实验室，湖北武汉430079

1 引 言

干涉合成孔径雷达（interferometric synthetic aperture radar，InSAR）作为一种主动的微波遥感技术在地形测绘已经显现出越来越大的优势［1－2］。针对国家一些重大的测绘生产任务，如西部测图任务，在坡度陡、高差大的西部地区进行高精度的DEM制图已经成为一项极具挑战性的任务。同时，由于TerraSAR－X、COSMO－SkyMed等高分率雷达遥感卫星的相继发射，这些具有强干涉能力的雷达传感器使高精度的DEM地形制图已经成为可能。

在应用InSAR技术生成DEM过程中，主要的误差源来自于相位误差［3－5］。相位误差主要由三部分组成：① 轨道误差引起的相位趋势误差；② 由于时间、空间去相关和热噪声造成的误差；③由于主从影像成像时间产生的大气差的扰动所引起的误差。为了更加精确地去除相位误差，外部DEM已经被应用到雷达干涉测量中，文献［6］提出了一种用粗的DEM方法减少干涉图的局部相位带宽，较低的相位带宽减少了相位残余，从而使相位解缠更加容易；文献［7］则有用外部的DEM减少高程的搜索范围；文献［8］使用SRTM通过构建线性模型最大程度地去除了相位误差。然而，针对坡度陡、高差大的复杂地形地区，单一的线性模型已无法拟合整个地形趋势，同时，在高分辨率雷达数据处理步骤中，如建立外部DEM模拟相位和干涉解缠相位的对应关系，在轨道参数精度低的条件限制下也严重影响最终的DEM精度。

本文针对坡度陡、高……