第一作者简介: 孙亚飞(1986-),男,研究生,主要从事TanDEM-X双站InSAR数据处理与应用。Email: yfsun2013@126.com。
较全面、系统地介绍了TanDEM-X/TerraSAR-X双站SAR科学计划,重点涉及其科学目标、TanDEM-X卫星参数、轨道结构以及干涉数据获取模式等相关内容,并讨论了双站SAR成像、极化InSAR和数字波束成形等干涉测量新技术及其研究进展。这些双站SAR新技术的实现将大大地推动SAR干涉测量在全球地形测绘、冰川学、海洋学及地质学等领域中的应用。
This paper comprehensively describes the scientific research plan of the TanDEM-X/TerraSAR-X bistatic SAR mission, with emphasis placed on its scientific objectives, TanDEM-X satellite parameters, orbital configuration and data acquisition modes. Then, the advantages of the new InSAR techniques including bistatic SAR, Pol-InSAR and digital beam forming are briefly discussed and the progress of the techniques of TanDEM-X bistatic SAR interferometry is analyzed. These new SAR techniques will greatly promote the application potential of SAR interferometry in Earth sciences such as global topography mapping, glaciology, oceanography and geology.
近几十年来, 星载SAR技术与应用取得了突破性进展, SAR影像已成为十分重要的对地观测数据源[1]。干涉SAR(InSAR)是雷达遥感的重要发展方向。2000年美国航天飞机雷达地形测图计划(SRTM)展示了InSAR技术在地形测绘方面的卓越性能[2]。作为X-SAR/SRTM计划的延续, 德国宇航局分别于2007年6月15日和2010年6月21日成功发射了TerraSAR-X(TSX)卫星和TanDEM-X(TDX)卫星, 2颗卫星编队飞行, 构建双站SAR干涉系统, 获取全球范围的高分辨、高精度DEM数据[3]。这些数据对于地形测绘、冰川学、水文学和海洋学等地球科学领域具有重要意义。
本文较为详细地阐述了TanDEM-X双站SAR科学计划, 重点分析了双站SAR成像、极化InSAR和数字波束成形等TDX/TSX双站SAR干涉测量新技术, 并评述了该技术在地球科学领域的应用潜力及国内外最新研究成果。
2010年6月21日德国宇航局成功发射TanDEM-X卫星, 与之前的TerraSAR-X卫星进行双星星座模式运行, 开创了真正意义的星载双站SAR干涉测量[4]。与其他卫星星座运行方式相比, TanDEM-X双站SAR干涉测量采用近距离螺旋结构运行, 消除了时间去相干和大气去相干等去相干源, TDX/TSX卫星飞行效果如图1所示。
与2000年NASA获取的SRTM数据相比, TanDEM-X DEM由南北纬60° 范围的覆盖扩展到全球覆盖, 包括南北极的全部地区, 且分辨率和高程精度分别提高到12 m和2 m[5], 根据用户需求, 甚至可以分别达到6 m和0.8 m。未来几年内, TanDEM-X DEM将会替代SRTM, 成为新一代全球DEM基础数据。此外, TanDEM-X科学计划将会探索一些新SAR技术, 如双站SAR成像、极化InSAR、数字光束成形和超高分辨率图像等[6]。
TanDEM-X卫星的系统参数和成像模式与TerraSAR-X卫星的基本一致, 但TanDEM-X卫星安装了协调星座间运行的冷气推进系统, 以及接收TerraSAR-X卫星GPS位置信息和状态信息的天线装置[7]。2颗卫星的基本参数如表1所示。
| 表1 TanDEM-X/TerraSAR-X卫星基本参数 Tab.1 Basic parameters of TanDEM-X/TerraSAR-X satellites |
TanDEM-X双站SAR科学计划的首要目标是获得一个全球的、连续的和高分辨率的DEM。TanDEM-X卫星设计寿命为5 a, 与TerraSAR-X卫星至少有3 a的编队飞行时间[8]。在3 a的重叠期中, 首先进行6个月的试运行阶段, 之后2 a内将会调整空间基线范围, 获取全球DEM数据。其中, 在重叠期的最后半年, 将获取地形复杂地区的DEM数据[9]; 重叠期之后, TanDEM-X计划将专注于获取大基线实验数据。
由表2可知, 第一年的垂直基线范围为200 ~ 450 m, 第二年的垂直基线将会扩大到第一年的1.3倍, 范围为450 ~ 600 m。之后半年内, TanDEM-X计划将调整垂直基线范围为500 ~ 3 000 m, 以减小阴影、叠掩和透视收缩所引起的误差, 获取复杂地形的DEM。3 a后, TerraSAR-X卫星如果仍能与TanDEM-X卫星保持编队飞行, TanDEM-X计划将进行各种基线的数据获取实验。欧空局为使其保持长期编队飞行, 在TerraSAR-X卫星停止运行后, 将会发射TerraSAR-X2继续维持编队飞行[10]。
| 表2 TanDEM-X任务时间表 Tab.2 TanDEM-X mission timeline |
TanDEM-X和TerraSAR-X卫星在非常接近的轨道面上协同运行[6], 通过调整双星编队飞行结构的节点线角度、近地点角度、轨道扁率和双星相位, 结合在不同升降轨节点偏离轨道平面的位移和在不同扁率矢量垂直方向放射状的偏离, 这种双星编队形成了一种沿轨道螺旋状卫星之间的相对飞行[11], 形成螺旋(helix)飞行结构, 如图2(a)所示。该轨道结构不存在相交情况, 保证了2颗卫星安全运行。2颗卫星空间基线的顺轨方向分量小于1 km, 甚至在预定纬度上可以小于250 m, 交叉轨方向分量在200 ~ 600 m之间, 以满足不同应用的需求。Helix编队可保持稳定的高程模糊度(height of ambiguity), 且使南北半球同一纬度能以相同的有效基线观测, 如图2(b)所示[12]。
TDX/TSX卫星编队采用3种协作模式获取干涉数据(图3), 并且这3种模式与聚焦、条带和宽幅SAR成像模式自由组合[6]。其中双站模式是TanDEM-X计划中生成DEM的最常用模式, 而其他数据获取模式只进行系统标定和验证。
1)双站模式。如图3(a)所示, 在该模式下, 任何一颗卫星都可以作为雷达脉冲发射卫星, 2颗卫星同时接收回波信号, 消除了时间去相干源和大气去相干源, 但必须保证2颗卫星间PRF(pulse repetition frequency)同步。双站模式是TanDEM-X计划中获取DEM的标准模式。
2)单站追踪模式。如图3(b)所示, 单站追踪模式不需要保持2颗卫星同步, 可以单独运行。为避免2颗卫星雷达信号的干扰, 2颗卫星空间基线在顺轨方向分量必须大于10 km, 除了在植被、中高风速水面等地区外, 大部分地物都能保证较高的相干性。该模式与传统重复轨道干涉模式相似, 但在时间去相干方面远远优于重复轨道干涉模式。
3)双站交替模式。如图3(c)所示, 双站交替模式要求TanDEM-X和TerraSAR-X轮流发射雷达脉冲, 同时接收回波信号。2颗卫星轮流发射完一次脉冲称为一次交替, 每次发射脉冲的卫星获得的图像称为主图像, 相应另一颗卫星获得的图像称为辅图像。双站交替模式一次交替获取2幅主图像和2幅辅图像, 同一颗卫星的主图像和辅图像用来测定振荡器引起的相位误差, 进行双站SAR干涉测量精确相位校正。相位校正后, 2幅辅图像可以合为一个双倍PRF的辅图像。同时, 一次交替生成2幅具有不同高程模糊度的干涉图:
①一颗卫星的主图像与另一颗卫星的辅图像生成一幅干涉图, 干涉图的高程模糊度为
hamb=(λ r sin θ i)/B⊥ , (1)
式中: λ 为脉冲波长; r为斜距; θ i为入射角; B⊥为视线向垂直基线。
同样地, 另一对主图像和辅图像也生成一幅干涉图。通过合并2幅干涉图来提高相位校正精度。
②2幅主图像也可以生成一幅干涉图, 且相位模糊度为第一幅干涉图的2倍。该干涉图的高程模糊度为
hamb=(λ r sin θ i)/2B⊥ , (2)
双站交替模式采用交替发射脉冲方式, 获取2种不同相位模糊度的干涉图。这2种干涉图既能够提高干涉相位解缠质量和效率, 又能够区分一次散射和二次散射。
TanDEM/TerraSAR双站SAR新技术主要包括双站SAR成像、极化InSAR、数字波束成形、超高分辨率及多基线InSAR数据处理等[7]。
1)双站SAR成像。双站SAR成像技术是将主图像和辅图像融合为多角度观测图像的技术。基于雷达信号特征的双站雷达散射面对目标进行探测和识别, 通过比较双站SAR图像的后向散射系数进行图像分割和图像分类。双站SAR成像还可应用于海洋状态参数反演、地表粗糙度估算和地形坡度计算, 以及林业、立体航空摄影测量、气象学和大气学等领域。
2)极化InSAR[13]。采用多种极化数据减小干涉相位变化, 可以获取高质量干涉图和高分辨率DEM。结合星载或机载L波段极化SAR数据, 反演出森林树高和树冠层等农林参数, 生成数字表面模型(digital surface model, DSM)。
3)数字波束成形。数字波束成形是结合一些无方向的小天线模拟成一条大且带方向的天线。TerraSAR-X卫星和TanDEM-X卫星具有双天线接收设备, 能够获取4个中心的相位数据。数字波束成形能够获取一个高分辨率宽刈幅图像, 且能解决交替双星模式的时间冲突问题。
4)超高分辨率。2颗卫星接收同一地面点的回波信号具有不同的角度, 通过合成这2种信号生成超高分辨率SAR图像。
5)多基线InSAR数据处理技术。TanDEM-X可以获取宽幅SAR、聚焦干涉、双基线双站交替干涉和分离带宽干涉SAR等各种SAR数据, 扩大DEM的覆盖区域, 提高其分辨率。此外, TanDEM-X将进行多基线数据处理, 提高DEM精度, 获取高阶差分相位和双差分InSAR相干值。与之前ERS开创的重复轨道干涉测量一样, TanDEM-X将开创多样SAR干涉处理新时代。
根据实际需求的不同, 按照交叉轨干涉SAR技术和顺轨干涉SAR技术分别对TanDEM-X应用领域进行分析与展望。
交叉轨干涉SAR技术主要用来生成高精度DEM。对2条SAR天线获取的图像进行相位差值, 并将相位转换成高程生成DEM数据。精确的DEM或DEM副产品广泛应用于地形测绘、冰川学、水文学、海洋学、地质学、城区制图、导航、灾害应急管理、土地覆盖和植被等领域。
1)全球地形测绘。目前, 获取局部地区高分辨率和高精度的DEM相对容易, 但获取一个全球的、高分辨率、高精度和有统一数据格式及投影的DEM比较困难[14]。2000年NASA航天飞机获取SRTM DEM覆盖了全球南北纬60° 之间的范围, 但存在年数久远、分辨率低、精度低以及未能全球覆盖和有数据空洞等问题。TanDEM-X计划作为SRTM项目的延续, 其主要目的仍是全球地形测绘, 将生成精度满足HRTI-3(high resolution terrain information)标准的全球覆盖、12m高分辨率的DEM用于科学研究和商业用途[15]。
2)冰川学。冰川物质平衡、冰川气候相互影响和冰川盆地径流的监测及建模, 都需要精确的DEM, 且多期DEM可以进行时序物质平衡探测和冰川退缩监测[16]。而TanDEM-X脉冲信号为低穿透性的X波段, 可以提供冰川物质平衡变化和冰川退缩引起海平面上升方面的精确数据。DEM还可以作为一个重要输入参数[17], 用于极地冰川变化模拟。尽管CryoSAT-2数据也能提供冰川地形[18], 但分辨率较低, 而TanDEM-X可提供高分辨的DEM。
3)水文学。高分辨率DEM可用于区域洪水淹没风险评估和与水文有关的地貌特征提取(如水系网、集水区和地表湿润系数等)[19]。TanDEM-X提供的高精度DEM可缩小现有水文模型的差距, 开辟遥感在水文学方面新的应用领域。
4)海洋学。目前只有ERS SAR和ENVISAT ASAR星载单天线SAR能够进行洋流波形测量[20]。而TanDEM-X干涉数据有望克服传统测量方法的缺陷, 提高海洋波谱的估计精度。另外, SAR图像显示的精细洋流信息也可开辟新的应用领域, 譬如洋流-波形相互作用参数。TanDEM-X还可以提供海洋学中的2个重要参数, 即二维海洋波谱和风场参数。且TanDEM-X的SAR数据能够计算出二维海洋波谱的波高、波传播方向和波长参数以及高分辨率风场。
5)地质学。高精度的DEM是地质科学应用的基础数据。地面沉降、火山、地震和海啸等地质灾害都需要最新且精确的DEM, 用以计算灾害发生前后的变化信息[21, 22, 23]。目前, 只有TanDEM-X计划能够获取全球高精度的DEM数据。
6)城区制图。城市土地利用、城市规划、城市环境和城市安全等领域都需要精确的DSM[24, 25]。TanDEM-X将提供高分辨率且适合提取城市地区DSM的干涉对数据, 解决了阴影、叠掩和多路径效应导致SAR图像及干涉相位解译的困难。
顺轨干涉SAR技术是通过计算2幅同一地区、不同时刻的SAR图像相位差, 探测移动目标并估计其速度。与交叉轨干涉SAR相比, 顺轨干涉的2幅SAR数据获取时间方面存在短时间的差异。这种差异可以通过调整2颗卫星在顺轨方向的距离实现, 也可采用单颗卫星的双天线接收(DRA)设备实现。2颗卫星间距可以在0~1 000 m范围内调整, 以实现不同速率的动态目标探测和估算。顺轨干涉SAR可应用于监测海洋洋流、冰川流速和交通流量等领域。
1)洋流和河流流速提取。海岸带洋流地图集的制作、深海变化监测、潮汐发电机位置确定和不能现场实施的洋流和河流流速提取等, 都需要高分辨的流场数据[26]。顺轨干涉测量技术能够提供该类数据。2000年的SRTM项目已成功获取顺轨干涉测量数据, 而TanDEM-X顺轨干涉测量技术的优势在于能够提供更高质量的数据, 并可根据用户需求提供实验区的时序数据。
2)冰川流速监测。除冰川地形和接地线位置外, 冰川流速也是计算冰川物质平衡变化的重要参数[27]。TanDEM-X能够提供适合的时间基线和空间基线, 有效监测各种流速的冰川变化, 成为监测冰川流速的最佳手段之一。
3)交通流量监测。交通流量监测需要探测移动目标及其速度等参数[28]。TanDEM-X计划的双站顺轨干涉, 尤其2颗卫星上的DRA设备, 在没有外部参考数据情况下能够精确测量动态目标的前后偏移量、速度和加速度等参数。
自TanDEM-X卫星成功发射以来, 国内外对TanDEM-X双站SAR干涉测量研究越来越重视, 已在城区DEM提取、森林制图、生物量估算、冰川地形提取及质量平衡估算、火山监测、洋流和河流流速提取等领域取得了突出的研究成果。
1)城区DEM提取。城市扩张监测及城市测绘等城区管理方面均需要高精度DEM。由于传统重复轨道数据获取模式在时间去相干和时间基线等方面存在缺陷, 难以在城区DEM提取方面得到广泛推广。Rossi等[29]采用TanDEM-X双站SAR数据, 利用自适应多视技术, 结合TanDEM-X双站SAR高相干特性, 提取了德国柏林地区10 m分辨率DEM。
2)森林制图。碳排放量是影响气候变化的一个重要因素, 拥有有效的森林覆盖率监测手段尤为重要。TanDEM-X获取了全球覆盖、高分辨率的双站SAR数据, 克服了时间去相干、体散射去相关等相关源, 成为森林土地覆盖、森林分类的重要数据源。Schlund等[30]采用TanDEM-X双站SAR数据的相干特性, 有效地进行了森林分类及土地覆盖监测。
3)植被生物量估算。目前, 利用航空激光雷达、光学遥感等手段可以进行生物量的估算, 但受数据精度低、连续性差及云雨天气等的影响, 在大范围、高精度生物量估算方面受到限制。Solberg等[31]通过TanDEM-X双站SAR干涉测量, 提取了挪威森林的时间序列DSM, 并对森林随时间变化的DSM进行分析, 得到了森林DSM与生物量之间的关系。Caicoya[32], Lange[33]和Treuhaft[34]等也采用TanDEM-X数据对Boreal, Lardal及Brazil地区进行了植物生物量估算和植被分类。
4)冰川地形提取及质量平衡估算。山岳冰川和冰原的消融对全球气候变化和区域水资源、生态环境具有显著影响。GRACE重力卫星资料在监测南北极冰盖和冰架质量平衡方面发挥了重要的作用, 但受400 km分辨率限制, 难以准确估算山岳冰川的质量变化。利用TanDEM-X双站数据, Jaber W[35]和Jiang[36]等已成功提取了南美Patagonia冰原、青藏高原Puruogangri冰原、格陵兰Petermann冰川以及喜马拉雅Chhota Shigri冰川、阿尔卑斯Otztal冰川等的高分辨率、高精度DEM, 并与SRTM等DEM进行差值分析, 获取了较为可靠的冰川高程变化及质量平衡参数[35, 36, 37, 38, 39]。
5)火山监测。火山地形、体积及其增长率是火山研究的重要参数, 传统重复轨道InSAR数据重访周期为11~35 d, 较低的相干性很难监测火山变化。Charbonnier[40], Froger[41]和Kubanek[42]等利用TanDEM-X双站SAR提取了Merapi和Piton de la等火山的地形及喷发前后的体积变化等信息, 并对喷发火山进行了灾害评估。
6)洋流和河流流速提取。TanDEM-X与TerraSAR-X卫星能够构成单站追踪模式获取顺轨干涉数据, 且每颗卫星安装有双接收天线设备(DRA)。这种灵活多样的顺轨数据获取模式能够用来提取海洋洋流和河流流速信息。Runge[43]和Desroches[44]等利用TanDEM-X顺轨数据(ATI)和双接收天线数据(DRA)成功提取了洋流和河流流速。
TanDEM-X科学计划运行了全球首个星载双站InSAR系统, 这是高分辨率星载SAR发展的重要里程碑。本文详细介绍了TanDEM-X科学计划, 分析了该计划的科学目标、TanDEM-X卫星参数、Helix轨道结构和干涉数据获取模式。双站InSAR系统实现了双站SAR成像、极化InSAR、超高分辨率、数字波束成形和多基线SAR干涉等新SAR技术, 大大促进了星载InSAR技术的发展, 并在全球地形测绘、冰川学、水文学、海洋学及地质学等领域显示出了巨大的应用前景。
志谢: 本研究得到了德国宇航局(DLR)Michael Eineder博士和Dana Floricioiu博士的帮助与指导。本文所采用的TanDEM-X双站SAR数据由DLRTanDEM-X任务AO项目(No. XTI_LAND0413)提供,在此一并志谢。
The authors have declared that no competing interests exist.
| [1] |
|
| [2] |
|
| [3] |
|
| [4] |
|
| [5] |
|
| [6] |
|
| [7] |
|
| [8] |
|
| [9] |
|
| [10] |
|
| [11] |
|
| [12] |
|
| [13] |
|
| [14] |
|
| [15] |
|
| [16] |
|
| [17] |
|
| [18] |
|
| [19] |
|
| [20] |
|
| [21] |
|
| [22] |
|
| [23] |
|
| [24] |
|
| [25] |
|
| [26] |
|
| [27] |
|
| [28] |
|
| [29] |
|
| [30] |
|
| [31] |
|
| [32] |
|
| [33] |
|
| [34] |
|
| [35] |
|
| [36] |
|
| [37] |
|
| [38] |
|
| [39] |
|
| [40] |
|
| [41] |
|
| [42] |
|
| [43] |
|
| [44] |
|

