SNAP实现DInSAR

0.snap下载

直接下载最新版本即可(官网下载,非收费软件)

1.snap导入文件与文件准备

0.导入文件类型

此处需导入SLC文件,需下载SLC文件,而不是其他类型,第一次导入的时候用了其他类型卡了这里卡住了很久,最后知道是文件类型不同。直接导入zip即可,不需要解压。
SNAP实现DInSAR_第1张图片

2.Coregistration(图像配准)

这一步操作,简单来说就是使两幅影像在空间位置上对齐。
在SNAP软件中,在Radar菜单下选择Coregistration>S1 TOPS Coregistration>S-1 TOPS Coregistration。
随后在弹出的页面中,进行参数设置,包括:(1)设定主(master)从(slave)影像;(2)设定subswath;(3)设定所使用的卫星轨道文件;(4)设定辅助的DEM数据

1.设定主(master)从(slave)影像

SNAP实现DInSAR_第2张图片

2设定subswath

Read1,Read2分别填入需要对比的两个sar数据,如填入的数据有效,格式正确,TOPSAR-split会自动识别,变成如下图所示。
SNAP实现DInSAR_第3张图片

3.设定所使用的卫星轨道文件

(3)步骤均使用默认参数,也可自行导入数据

4.设定辅助的DEM数据

此处使用外部dem数据,使用的是长宁县dem数据。
SNAP实现DInSAR_第4张图片

5.输出

选择保存的路径
SNAP实现DInSAR_第5张图片
点击run,完成后save
SNAP实现DInSAR_第6张图片

新生成的Orb_Stack

3.Form the Interferogram (产生干涉)

这一步的操作也比较明确,即在上一步产生Orb_Stack的基础上,产生相应的干涉结果。具体为Radar菜单下:Interferometric>Products>Interferogram Formation只需要弹出的页面中设定好生成的产品名称即可(这里相当于使用默认的干涉参数,没用进行地形的消除)。
SNAP实现DInSAR_第7张图片
计算的时候很浪费时间,需要较长的等待。
SNAP实现DInSAR_第8张图片

4.TOPS Deburst

对于Sentinel-1而言,每一个swath是由多个bursts构成,这一步其实就是将多个burst拼接成为一整幅影像。在Radar菜单下选择Sentinel-1 TOPS>S-1 TOPS Deburst。可以看到在进行了deburst操作之后,确实有所改变。
SNAP实现DInSAR_第9张图片

5.Topographic Phase Removel

对于干涉产品而言,其相位部分同时受到形变(deformation)与地形(topographic)两个方面的影响,如果我们关心的是形变的话,要通过已知的DEM数据,去除相位中由topographic带来的影响。在SNAP中的具体操作为Radar菜单下,Radar>Interferometric>Products>Topographic Phase Removal
SNAP实现DInSAR_第10张图片

6.Multi-looking

点开Radar>SAR Utilities>Multilooking,选择相应的产品i_ifg_VV,q_ifg_VV以及coh_IW1_VV三个产品。其算法的本质即是通过平均多个点的值来相应地减少噪声。
SNAP实现DInSAR_第11张图片

SNAP实现DInSAR_第12张图片

7.Phase Filtering

在Radar菜单下Interferometric>Filtering>Goldstein Phase Filtering,使用默认参数进行相位滤波,
SNAP实现DInSAR_第13张图片

8.Geocoding

对干涉结果进行地理编码(Radar>Geometric>Terrain Correction>Range-Doppler Terrain Correction),得到最终结果。
SNAP实现DInSAR_第14张图片

9.最后的结果

SNAP实现DInSAR_第15张图片

你可能感兴趣的:(遥感,图像处理,其他)