isce D-InSAR

一、数据准备

二、配置文件insarApp.xml

我设置的解缠常用参数grass和snaphu

<topsApp>
  <component name="topsinsar">
    <property name="Sensor name">SENTINEL1</property>
      <component name="reference">
      <property name="orbit directory">/media/wx/Elements/menyuan1.17-1.05/raw</property>
      <property name="output directory">reference</property>
      <property name="safe">/media/wx/Elements/menyuan1.17-1.05/raw/S1A_IW_SLC__1SDV_20220105T111059_20220105T111126_041325_04E9B1_F077.SAFE</property>
    </component>
      <component name="secondary">
      <property name="orbit directory">/media/wx/Elements/menyuan1.17-1.05/raw</property>
      <property name="output directory">secondary</property>
      <property name="safe">/media/wx/Elements/menyuan1.17-1.05/raw/S1A_IW_SLC__1SDV_20220117T111059_20220117T111126_041500_04EF67_CDCA.SAFE</property>
    </component>
    <property name="demFilename">/media/wx/Elements/menyuan1.17-1.05/dem/demLat_N36_N39_Lon_E100_E104.dem.wgs84</property>
<!--    <property name="swaths">[1,2,3]</property>  -->
 <!--         <property name="region of interest">[37.1,38.6,99.5,102.7]</property> -->
         <property name="do unwrap">True</property>
         <property name="unwrapper name">grass</property>
         <property name="filter strength">0.4</property>
   <!--        <property name="geocode bounding box">[37.1,38.6,99.5,102.7]</property> -->
         <property name="correct phase error caused by ionosphere azimuth shift">1</property> 
</component>
</topsApp>

三、干涉处理

文件配置总的干涉步骤

topsApp.py insarApp.xml

四、干涉结果格式转化

4.1形变信号提取

sce干涉生成的形变信号包含在filt_topophase.unw.geo文件的波段2中，可以使用gdal_translate 转化成gmt绘图软件使用的网格文件。

gdal_translate -of GMT -b 2 filt_topophase.unw.geo filt_topophase.unw.geo.grd

生成的filt_topophase.unw.geo.grd是解缠后的相位信号，需要乘以（波长/４π)转化为形变，本例中采用的是哨兵数据，波长为5.56cm，所以los向形变为：

gmt grdmath filt_topophase.unw.geo.grd -0.44 MUL = los.grd

加负号将方向转换成由地面指向卫星。

4.2、相干性文件

isce干涉生成的相干性文件为:phsig.cor.geo

gdal_translate -of GMT -b 1 phsig.cor.geo cor.grd

可以通过cor.grd文件对形变信号做一个掩膜处理，去掉一些低相干性部分。

4.3 LOS方向向量

los向方向向量包含在los.rdr.geo中：

imageMath.py --eval='sin(rad(a_0))*cos(rad(a_1+90));sin(rad(a_0)) * sin(rad(a_1+90));cos(rad(a_0))' --a=los.rdr.geo -t FLOAT -s BIL -o enu.rdr.geo

enu.rdr.geo中波段１、２、３分别为e n u方向:

gdal_translate -of GMT -b 1 enu.rdr.geo e.grd
gdal_translate -of GMT -b 2 enu.rdr.geo n.grd
gdal_translate -of GMT -b 3 enu.rdr.geo u.grd

4.4 DEM转换为grd

gdal_translate -of GMT -b 1 demLat_N36_N37_Lon_E120_E121.dem.wgs84 dem.grd

4.5

sh mask.sh filt_topophase.unw.geo.grd 0.3

mdx.py topophase.flat.geo