(三)连接图编辑
对上一步生成的每一对干涉结果进行查看和筛选,筛选出效果差的,同时记录两个最好的。
文件存放的位置:…\SBAS_processing\interferogram_stacking\interf_tiff
每一对结果包括相干系数图(_cc)、相位解缠图(_upha)、滤波后的干涉图(_fint)等,最快的筛选方法只需要看upha图就好,相位解缠图是灰度图像,黑色代表没有解缠结果的地方,所以黑色越多,图像是越差的;灰色部分越多,效果是很好的,保留。(凭你的主观印象进行删除和保留)
(图1 效果很差,不要)
(图2 效果较好,保留)
(图3 效果最好的,保留并记录在.txt中)
使用工具:SBAS Edit Connection Graph ,进行连接图编辑。可以同时删掉多个,ctrl 多选选好后,点击remove pairs 就可以。
(四)第一步反演
SARscape5.6.2版本跑SBAS时,不再单独进行轨道精炼与重去平,而是在第一次、第二次反演和地理编码这三步中进行轨道精炼与重去平。
输入auxiliary.sml文件,控制点可以自己添加,也可以空着(系统自动选点)。这一步运行也很久,半天到一天时间。
(四)第二步反演
(五)地理编码
打开反演第一步的结果中的两个精度文件:FI_H_precision 和 FI_V_precision,这两个文件在 \SBAS_processing\first_inversion\inversion 路径中。
右击图层打开quick stats,查看直方图,大致呈正态分布就是对的。然后看高程和速率的阈值。(比如下图阈值是28)
‘
(六)结果分析
结果保存在 SBAS_processing\geocoding\raster 下,要看的主要有两个数据结果。
① G_SI_velocity_geo 数据,该数据的像元值是平均位移速率,单位是mm/y。data值的正负代表累积量沉降量。
② G_disp_geo_meta 时序形变结果,即从开始到结束时间的累计形变量。打开SARscape/General Tools/Time Series Analyzer/Raster工具,定位到目标像元,点击plot,就可以得到该像元的时间序列上的变化情况。