地球信息科学学报 ›› 2014, Vol. 16 ›› Issue (5): 797-805.doi: 10.3724/SP.J.1047.2014.00797

• • 上一篇    下一篇

邻轨PS-InSAR地面沉降结果拼接处理方法与实验

熊思婷(), 曾琪明, 焦健*(), 章晓洁   

  1. 北京大学遥感与地理信息系统研究所,北京 100871
  • 收稿日期:2013-11-15 修回日期:2014-03-26 出版日期:2014-09-10 发布日期:2014-09-04
  • 通讯作者: 焦健 E-mail:xiongsiting90@gmail.com;jiaojian@pku.edu.cn
  • 作者简介:

    作者简介:熊思婷(1990-),女,硕士生,研究方向为PS-InSAR地表形变等研究。E-mail:xiongsiting90@gmail.com

  • 基金资助:
    国家自然科学基金项目(41171267);中国科技部—欧洲空间局合作“龙计划”三期项目(10665)

Research on Connecting PS-InSAR Results from Adjacent Tracks for Land Subsidence Monitoring

XIONG Siting(), ZENG Qiming, JIAO Jian*(), ZHANG Xiaojie   

  1. Institute of Remote Sensing and Geographical Information System, Peking University, Beijing 100871, China
  • Received:2013-11-15 Revised:2014-03-26 Online:2014-09-10 Published:2014-09-04
  • Contact: JIAO Jian E-mail:xiongsiting90@gmail.com;jiaojian@pku.edu.cn
  • About author:

    *The author: WANG Jiacheng, E-mail:shanqiangw@fync.edu.cn

摘要:

目前,我国许多经济发达地区都面临着地面沉降灾害的困扰,沉降范围扩大,程度日益加剧,逐渐成为城市发展中亟待解决的问题。永久散射体合成孔径雷达干涉测量(Persistent Scatterer Interferometric Synthetic Aperture Radar,PS-InSAR)作为地表形变测量的主要手段之一,在地面沉降监测中发挥着重要作用。而单轨星载SAR影像成像幅宽有限,在开展大范围地面沉降调查时需要将多轨道PS-InSAR沉降速率图进行拼接。本文重点讨论了入射角效应和参考点差异对PS-InSAR沉降结果的影响,分析了相邻轨道PS-InSAR沉降速率拼接中存在的PS点位置差异和沉降量偏移,鉴此,提出了采用区块法和插值法对异轨重叠区的形变结果求差的思路,以及基于现有软件PS-InSAR地面沉降速率的跨轨拼接处理流程,利用广东珠三角地区ENVISAT ASAR数据进行了实验分析。结果表明,相邻轨道入射角不同会造成沉降量的差异,在多轨道情况下对沉降量影响增大,因此,在拼接过程中需要进行入射角纠正。本文提出的区块法和插值法能有效地求解异轨重叠区的形变差,结果表明区块法优于插值法;相邻轨道参考点差异会造成沉降量偏移,通过区块法或插值法求差可以消除该偏移量。本文提出的拼接流程可将多轨道PS-InSAR地面沉降速率统一到同一基准下,从而获得大范围一致的地面沉降速率。

关键词: PS-InSAR, PS点, 地面沉降, 沉降速率, 跨轨拼接

Abstract:

At present, more and more economic developed areas are facing serious land subsidence problems with increasing coverage and intensity, which prevents the application of sustainable development in these areas. Thus, the land subsidence problems should be controlled urgently. PS-InSAR technology, as one of the most important tools to detect land subsidence, has played a significant role in land subsidence monitoring. However, multi-track PS-InSAR results must be connected to obtain the large scale subsidence map when an extended land subsidence survey area is required, since spaceborne SAR image from only one track is too limited in scan range to cover the whole area. This paper analyzes the processing chain of PS-InSAR and points out two key problems when connecting multi-track PS-InSAR results from adjacent tracks. First, many PS from adjacent tracks are physically different, which is mainly caused by different incidence angles and coherences of adjacent tracks. Therefore, it is hard to find the tie point in the overlapped area. Second, there is a subsidence offset between two adjacent tracks according to the overlapped area, which is attributed to the physical differences of PS and the processing parameters of PS-InSAR, such as the reference point, the temporal coverage and so on. Among these factors, impacts caused by different incidence angles and reference offset are emphasized. Two new methods, Block Method and Interpolation Method, are proposed in this paper to calculate the differences between subsidence velocities from adjacent tracks caused by the physical differences of PS in the overlapped area. This can avoid the difficulty of finding tie point in the overlapped area. Moreover, subsidence velocities are adjusted to mitigate the impacts from different incidence angles. Finally, based on the existing software, a complete procedure for multi-track PS-InSAR subsidence results connection is demonstrated. At the end of this paper, ENVISAT ASAR datasets of Pearl Delta River in Guangdong province are processed and analyzed in experiments. Experimental results show that different incidence angles could affect the land subsidence from adjacent tracks, so it is necessary to correct the differential subsidence caused by the incidence difference, especially when multiple tracks are involved. Block Method and Interpolation Method are all effective in calculating the differences of subsidence velocities, and Block Method was superior to Interpolation Method through statistical analyses. Meanwhile, different incidence angles and reference offset are effectively mitigated to obtain the conformed subsidence velocities using these two methods. Consequently, PS-InSAR subsidence velocities from two ENVISAT tracks are connected effectively, which shows the effectiveness of the processing procedure proposed in this paper.

Key words: PS-InSAR, PS point, land subsidence, subsidence velocity, multi-track connection