采用PS-InSAR或SBAS-InSAR对同一个地区得到不同时间段内的形变结果,有时候我们需要将这些不同时间段内形变结果进行合并,得到以最早时间点为基准、覆盖整个时间跨度的形变结果。
SARscape从5.5.3版本开始提供栅格或矢量地表形变结果合成工具,计算得到时间序列累计沉降量和平均形变速率,平均速度的估算是通过对合并后的时间序列进行线性回归计算得出的。
需要合并的形变结果需要满足以下条件:
(1)PS或者SBAS生产的矢量或者栅格结果;
(2)最低输入数据要求为 2个升序数据或2个降序数据。该合并过程可基于由不同卫星传感器获取的结果进行计算;
(3)为了获得可靠的结果,输入形变结果的视线入射角度必须具有可比性,角度差异在 5 至 10 度范围内;
注:如果矢量地表形变结果由多个shp文件组成,这些文件需要放在时序处理过程中创建的原始文件夹中,软件会自动将它们重新组合一个shp文件。
本文使用的例子为矢量格式,如下图所示。
图1:例子数据(左边为2017.3~2018.10年形变结果,右边为2018.11~2020.6年形变结果)
(1)在Toolbox中,选择/SARscape/Interferometric Stacking/Post Processing Tools/Shape Combination工具。打开Shape Combination工具面板。
注:如果是栅格结果文件,选择/SARscape/Interferometric Stacking/Post Processing Tools/Meta Combination工具。
(2)数据输入(Input Files):点击按钮,输入需要合并的矢量结果文件,如SI_XX*_Y.shp 或者shp文件。
(3)参数设置(Parameters):主要参数(Principal Parameters)。
- 组合类型(Combination Type):Merge。
- 匹配圆圈半径(Radius (m) ):匹配矢量点最大半径,在这个半径内匹配输入的点,输出的点的位置为输入圆圈的重心。
- 半径搜索纬度(Radius Search Dimension):它定义了有效半径约束的维度(2D或3D)。如果选择3D,建议输入的数据为PS处理结果。
- 矢量文件最大点个数(Max Point in Shape):设置单个shp矢量文件最多储存点个数,超出最大数会自动输出多个shp文件。
(4)数据输出(Output Files):选择输入目录和设置输出文件名。
(5)单击Exec按钮开始执行。
图2:参数设置
处理过程中可能会出现以下提示,这个提示意味输入的结果文件中有些点在设置的圆圈半径内没有匹配到点。这里可以忽略这个提示。
图3:提示信息
如下图为合并后的结果。
图4:合并后的结果
分别从三个矢量上获取相同点的累计沉降量。
图5:相同点合并前后的累积沉降量曲线