SBAS实操

数据源：欧空局哨兵一号，升轨，范围：POLYGON((116.1326 40.0333,116.1877 40.0333,116.1877 40.064,116.1326 40.064,116.1326 40.0333))

---

1.Preference--参数设置

Load Preferences --- Sentinel TOPSAR

---

2.导入数据

2.1 SLC数据

2.2 研究区域--shp文件

2.3 导入

选Open Folder可以批量导入

参数默认

---

3.生成连接图

Input files：

Optional不用，图像多Delaunay 3D选True，注意最大时间基线填365，Normal Baseline一般填45

---

4.干涉工作流

1.Input files

输入上一步生成的_SBAS_processing文件夹下的.sml文件

2.DEM格式问题

下载的数据通常为.tif格式，需要导出为_dem格式

1.直接拖入tif

2.另存为.dat文件，保存时选择文件类型为.其实就是.dat，后缀会自动加上去

3.打开Orginal ENVI Format

4.导入.dat文件，然后设置参数

这里是直接利用的ALOS卫星的Ellipsdidal DEM（椭球DEM）。

关于DEM类型可以参考科普：PALSAR RTC DEM Information | Alaska Satellite Facility

通常下载到的DEM是基于大地水准面，而在进行干涉处理和高级InSAR时序分析时，都需要在参考椭球面下进行，所以需要对DEM进行转换，具体可以参考：SARscape5.6版本导入外部DEM数据 - ENVI-IDL技术殿堂 - 博客园 (cnblogs.com)

输出时手动将后缀改为_dem

*注：关于坐标系，WGS84 包括参考椭球体和大地水准面，但大地水准面已更新/替换为更精确的 EGM96。因此，人们在谈论WGS84时通常指的是参考椭球体

3.参数设置

解缠方式默认选择最小均匀流法

但是如果前面Delaunay 3D选择了True，这里解缠方式必须选择Delaunay MCF

---

5.轨道精炼和重去平

 注：生成GCP点

先采用Generate Ground Control Points模块结合干涉工作流生成的cc和uhpa文件生成GCP点，再进行轨道精炼，如果在轨道精炼模块下生成GCP点会卡死，估计是ENVI5.6的bug

注意输入文件和Reference文件不是此文件夹下的tif文件：

而是该文件夹下的无后缀文件：

总览：

这一步是估算和去除残余的恒定相位和解缠后还存在的相位坡道

均为默认操作：

---

6.第一次反演

这一步是 SBAS 反演的核心步骤，第一次估计位移速率和残余地形，用来对合成的干涉图进行去平，重新做相位解缠和精炼，生成更优化的结果，用于下一步计算。

均为默认操作：

---

7.第二次反演

这一步的核心是计算时间序列上的位移，在第一步得到的形变速率(_disp_first)基础上，进行定制的大气滤波， 从而估算和去除大气相位， 得到更加纯净的时间序列上的最终位移结果。在这里 “Displacement GCP file” 是用来去除还有残余的相位或相位斜坡，使用轨道精炼的时候所用的轨道控制点。

均为默认操作：

---

8.地理编码

均为默认操作：

结果存放在该文件夹下：

1.打开G_SI_velocity_geo文件，是平均位移速率图

进行彩色渲染并查看 DN 值，单位为 mm/year：

插入图例，修改图例格式：：

得到：

2.打开G_disp_geo_meta文件

双击/SARscape/General Tools/Time Series Analyzer/Raster 工具

点击 plot，绘制任一像元监测时间段内的形变历史图，单位为mm：

注：每一个形变值都是该时相相对于第一个时相的形变，负值沉降，正值抬升，相当于累计形变，不是一个周期的形变！参见《SARscape雷达图像处理入门与提高202302》:

各文件表示的涵义：

*9.关于误差校正、结果优化

可参考该文章里的：

6.2 连接图编辑的参数说明

6‘（11）单击 Exec 按钮执行处理。点击 Exec 进行处理

引用声明：该段内容为文档：[ENVI-SARscape入门教程-第8章短基线SBAS]
如果侵权请在评论区或私信联系我删除！谢谢