Step2 SBAS-InSAR处理流程

三、import sentinel

不同的SAR传感器有自己的数据格式，要先将数据导入成sarscape处理的标准数据格式。在处理之前还要把精密轨道数据准备好。5.6.2版本可以preference里面设置好路径，处理时不用再添加轨道数据，因为系统会自动添加。并同时进行子区域裁剪，就会很快，一景不超过十分钟吧。现在我们正式开始：

打开SARscape——import data——sar spaceborne——single sensor，我用的是哨兵一号数据，选择的是sentinel-1。

输入文件：sarscape5.6.2可以直接输入.zip文件，如果不是这个版本要先解压再输入manifest.safe文件。

感兴趣区域可输入地理坐标系的 shapefile 矢量文件或.kmz 文件，节省时间和硬盘空间。

极化方式根据数据而定（下载slc数据时候可以看见）

5.6.2 版本新增了“Skip Sample Selection”，默认是False，也就是说如果在Optional Files界面输入了子区域文件，得到的就是子区域裁剪的结果，不需要再单独进行子区域裁剪了（我用的地理坐标的.shp，但是报错invalid xml error，所以在Arcgis数据库里绘了新的feature，导出为shp文件，就搞定了。）如果要单独裁剪，用的是下图所示工具（可以用SAR坐标也可以地理坐标）。

对数据重命名（Rename the File Using Parameters）：True。软件会自动在文件名中增加“_VV_slc”。

​结果是这样：

四、SBAS正式处理

之后的工具都在这里：

(一)生成连接图

至少15幅影像结果才比较准确，所以我一次处理两年（17年3月-18年12月）的。“Super Master file”不用输入，5.6.2会自动选择一副主影像。时间基线最大设置120天（本人学艺不精，确实不知道如何设置，如果有朋友知道希望也能给我解惑，以普遍理性而言，时间基线越小，相关性越好），空间基线在一定阈值（最大空间基线的10%。最大基线可以从基线估算工具结果看到）内越大越好，我设为5。

结果如下：

（二）干涉处理流

输入上一步生成的 auxiliary.sml数据，Optional Files 不选择（5.6.2会自动选择控制点），输入_dem文件，参数设置如图。

补充：跑完一次sbas之后，效果不好，重新调整解缠方法，使用Delaunay MCF，如下图。

参考文献：施慧宇,王延霞,杨海燕等.相位解缠算法对比研究[J].黑龙江工程学院学报,2022,36(01):9-13.DOI:10.19352/j.cnki.issn1671-4679.2022.01.002.——结果表明相较于路径跟踪法和最小二乘法,Delaunay MCF算法在准确性和解缠运行时间上表现更为突出。

效果好了特别多。

这一步处理一到两年的InSAR数据需要1天，内存请准备好，否则会报错。