1. 数据准备
创建SLC、ORBITS、AUX、Orbits四个文件夹
-s SLC_DIRNAME #A folder with downloaded Sentinel-1 SLC’s.
-o ORBIT_DIRNAME #A folder containing the Sentinel-1 orbits. Missing orbit files will be downloaded automatically
-a AUX_DIRNAME #A folder containing the Sentinel-1 Auxiliary files
-d DEM_FILENAME #A DEM (Digital Elevation Model) referenced to wgs84
文件夹图片如下:
1.1 创建DEM文件夹
在相应的官网下载每度dem,在ISCE数据处理中可以从USGS的网站上下载这些地形数据dem数据官网(需要账号,注册一下即可)或者可以直接在网址后面接上要下载的地形文件名称https://e4ftl01.cr.usgs.gov/MEASURES/SRTMGL1.003/2000.02.1/N10E031.SRTMGL1.hgt.zip 官网之前也可以python自动下载(最近拒绝链接),然后,进行dem每度数据拼接,代码如下:
mkdir DEM; cd DEM
dem.py -a stitch -b 36 39 100 104 -r -s 1 -c
rm demLat*.dem demLat*.dem.xml demLat*.dem.vrt
cd ..
1.2 创建SLC文件夹
在欧空局官网自动下载Sential-1A/B数据(账号需要自己注册),网址两个1.哨兵数据 2.asf官网,放入SLC文件夹。
1.3 创建Orbits文件夹
方法1.设置轨道文件夹为空,在后续的stackSentinel.py命令下自动下载轨道文件。
方法2.在对应的轨道数据网站下载,网址1.批量轨道数据 2.欧空局轨道数据
下载完放入Orbits文件夹
1.4 创建AxuDir文件夹
一般为空,2015年3月前的需要下载。具体内容官网是这么说的:下面的校准辅助(AUX_CAL)文件用于天线模式校正,以补偿IPF版本002.36的SAFE产品的距离相位偏移(主要用于2015年3月之前获取的图像)。如果您所有的SAFE产品都来自另一个IPF版本,那么就不需要AUX文件。
2.开始预处理
创建Process文件夹,进行数据的预处理,代码如下,-b后面是形变范围区域,由于ESA的哨兵轨道有时候会变动偏移,同一景影像几年以后会偏移,所以,我们需要选择其中对应的形变即可,另外可以减少数据处理,提高效率。
stackSentinel.py -s ../SLC/ -d ../DEM/demLat_N36_N39_Lon_E100_E104.dem.wgs84 -a ../AuxDir/ -o ../Orbits -b '37 38.5 100.3 102' -W slc
stackSentinel.py -s ../SLC/ -d ../DEM/demLat_N36_N39_Lon_E099_E103.dem.wgs84 -a ../AuxDir/ -o ../Orbits -b '36.43 38.4 99.43 102.59' -W slc
stackSentinel.py -s ../SLC/ -d ../DEM/demLat_N36_N39_Lon_E099_E103.dem.wgs84 -a ../AuxDir/ -o ../Orbits -b '36.44 38.4 99.44 102.57' -W slc
最终生成run_files文件夹,里面有13个文件夹,图片如下:
给run_files文件夹赋于权限:
sudo chmod 777 -R run_files/
然后,对每个文件进行依次处理,这个处理的时间很长,慢慢等待,命令如下:
./run_01_unpack_topo_reference
最终,处理完生成merge文件夹,是我们所需的文件,图片如下:
这就是预处理的所有步骤。
参考的文章、视频、官网如下:
InSAR_ISCE软件数据处理流程之:利用哨兵数据获取2015年尼泊尔地震同震形变场