利用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。

1.2、dem数据准备

参考上一篇博文

二、干涉像对生成

干涉像对的生成步骤为：首先计算所有满足时间基线的干涉对，计算这些干涉对的空间基线长度，再从中挑选空间基线长度满足要求的像对进行干涉处理，生成干涉图。

2.1、时间基线与空间基线选取

生成干涉像对之前要确定最长时间间隔，即干涉像对的最大时间基线，一味追求尽量多的时间基线会导致干涉像对的相干性降低，影响结果质量。最大时间基线的选取要考虑研究区域的地形地貌特征以及使用的影像数据的波长，一般在寒冷干燥的区域时间基线可以适量放大，在气候湿润、草木茂盛的地区时间基线要适量减小，本例中最大时间基线选择是200，即生成所有时间基线不超过200天的干涉对（其实这一步不需要计算完成的干涉图，只需要计算这些干涉对的空间基线即可）。

首先需要更改SAR影像文件名，将数据获取时间作为影像的文件名。如图，需要将文件名从左边更改为右边所示。

mv S1A_IW_SLC__1SDV_20170425T050421_20170425T050448_016296_01AF5C_3FD0.SAFE  20170425

像这样一个一个改比较麻烦，我们可以将safe文件写入一个文本文件，再用foreach循环读取批量更改，脚本如下：

#!/bin/csh 
rm -f list_all
ls > safe.list
foreach safe (`cat safe.list`)
        echo $safe
        set slc=`echo $safe | awk '{print substr($1,18,8)}'`
        e