利用SBAS-InSAR技术获取Etna火山形变特征——基于ISCE+GIANT软件

时序InSAR数据处理过程——ISCE + GIANT 处理软件

ISCE 软件是一款常用的干涉处理软件，基于该软件可生成一系列干涉图，再通过时序分析软件（GIANT）即可完成时序InSAR数据处理过程。不过较为可惜的是目前GIANT软件的开发人员大部分被抽调去进行ISCE的维护与研发，GIANT软件已经很久没有更新过，本人暂时还没有在ubuntu18.04上成功安装过（只在ubuntu16.04上成功过）。如果无法正常安装giant，可以尝试利用stamps等软件进行时序分析。

一、数据准备

1.1、SAR数据准备

本文以意大利Enta火山为例，通过哨兵数据获取其喷发前的地表形变特征，以2017年与2018年的喷发事件作为研究时间窗口，下载得到覆盖火山口的降轨影像共41景（如图所示），轨道号为124。

全部更改完成后：

2.2 干涉处理

参考上一篇博文
本例中我们需要计算多个干涉对，但是原理上与计算单个干涉对是一样的，跟之前计算基线长度一样，加一个循环让计算过程实现自动化，能大大提高效率，需要脚本的可以邮件找我要。
本例中干涉结束一共得到358幅干涉图，需要将merged文件夹中生成的相位文件　filt_topophase.unw.geo和相干性　topophase.cor.geo　提取出来，并命名为类似20170612-20171103.unw　和　20170624-20170811.cor的格式（即主影像在前从影像在后的格式），方便后续时序处理。

注意：在进行干涉处理时一定要限制范围！就算是需要计算整幅影像的干涉图，同样需要限制范围，这个范围可以包含一整幅影像（可以从safe文件中的kml文件中查看影像四个角的经纬度），不然会导致时序分析那一步报错。

11行中data表示相干性文件与干涉图的储存路径
xlim 和 ylim表示干涉图的大小(可从data/20170419-20170425.unw.xml获取);因为在干涉处理中我们对干涉图的范围进行了约束，所以生成的干涉图的大小都是一致的。
rxlim 和 rylim是指研究范围内相干稳定的地块。
cohth　指相干性阈值

20170419 是指参考主影像，一般习惯将获取时间最早的一幅作为sbas-InSAR技术中的参考主影像。

3.1.2　userfn.py

userfn.py文件包含以下内容

3.1.3　ifg.list

ifg.list文件是指在生成空间基线步骤中生成的baseline.list文件，只需要将文件名更改为ifg.list就可以了。

这样我们就得到41个包含形变时间序列的文件，因为我们把获取时间最早的那幅设定为参考主影像，所以1.recons里面全部是0。
现在只需要将这些时间序列绘图制作出来即可

可以看出，2017年火山口并没有出现明显形变特征，到2018年接近火山喷发的时间，火山口才表现出强烈的地表隆升现象，最大视线向抬升超过40mm。

来源：CUG_施贺青