ISCE-stamps的sbas时序insar基本流程(一):ISCE的预处理

5 篇文章 3 订阅
1 篇文章 0 订阅

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年尼泊尔地震同震形变场

ISCE2StaMPS小基线集时间序列处理(InSAR, SBAS)

isce_github官网

评论 14
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值