InSAR个人笔记（三）ISCE之DInSAR处理--topsApp.py

• 数据下载

• 1、哨兵1 slc数据下载

研究区的选定：首先确定研究区内的范围（经纬度），其次确定下载数据的时间，最后保证两张SAR影像是重复轨道拍摄的（哨兵单颗卫星时间周期为12天，星座为24天）

（1）下载网址：

SAF ：https://search.asf.alaska.edu/#/?zoom=5.000¢er=123.354,25.990https://search.asf.alaska.edu/#/?zoom=5.000¢er=123.354,25.990

欧空局：https://scihub.copernicus.eu/dhus/#/homehttps://scihub.copernicus.eu/dhus/#/home

2、轨道数据下载

（1）格式说明：

POD Precise Orbit Ephemerides（POD精密定轨星历数据）每个文件覆盖26个小时（一整天24小时加上一天开始前1小时和一天结束后的1小时）

POD Restituted Orbit（POD回归轨道数据），文件在接收到GNSS数据的3小时内产生。

2.下载网址

Index of /aux_poeorb/https://s1qc.asf.alaska.edu/aux_poeorb/  ctrl+f可根据日期检索

https://scihub.copernicus.eu/gnss/#/homehttps://scihub.copernicus.eu/gnss/#/home  用户名密码：gnssguest

3、DEM数据下载

DEM数据下载：

（1）STRM(官网):https://earthexplorer.usgs.gov/    

（2）地理空间数据云： http://www.gscloud.cn/search

（3）修改经纬度直接下载（左下角为起点，默认覆盖范围为1°×1°）：https://e4ftl01.cr.usgs.gov/MEASURES/SRTMGL1.003/2000.02.11/N31E094.SRTMGL1.hgt.zip

注：下载到dem数据是1°*1°的格式，范围大的区域需要将其进行拼接，用如下命令拼接（s n w e）：

dem.py -a stitch -b 27 30 83 86 -r -s 1 -c

• 二、干涉处理（topsApp.py）

• 1、参数文件（xx.xml）

• 

• 2、执行topsApp.py xx.xml

• 注：topsApp.py 各个步骤（https://gitlab.com/comet_licsar/licsar_documentation/-/wikis/ciw2019/isce）

• 

• 三、信号处理

  （1）格式转换

  形变信号包含在filt_topophase.unw.geo文件的中，可以使用gdal_translate 转化成gmt绘图软件使用的网格文件（grd）

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

• （2）形变相位与形变量之间的转换

• 相位信号需要乘以（波长/４π)转化为形变量，哨兵1的波长为5.6cm。命令如下：

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

• 四、GMT绘图（sh文件）

  sh xx.sh