三、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天,内存请准备好,否则会报错。