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

置顶菲林可乐

已于 2025-02-10 21:20:11 修改

阅读量4.8k

点赞数 5

分类专栏： InSAR isce sentinel 文章标签： ubuntu

于 2022-10-03 09:09:51 首次发布

本文链接：https://blog.csdn.net/m0_60086135/article/details/127149180

版权

InSAR 同时被 3 个专栏收录

23 篇文章

订阅专栏

isce

5 篇文章

订阅专栏

sentinel

1 篇文章

订阅专栏

本文详细介绍使用哨兵-1A/B数据进行InSAR处理的准备工作，包括创建必要的文件夹结构、下载DEM数据及哨兵数据等。此外，还介绍了如何通过stackSentinel.py脚本进行数据预处理，以提高后续处理效率。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

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数据，请引用我的文献，非常感谢
Wang X, Li S, Tao T, et al. Interseismic strain accumulation across the Tuolaishan–Lenglongling segment of the Qilian–Haiyuan fault zone prior to the 2022 Mw 6.7 Menyuan earthquake from Sentinel-1 InSAR time series[J]. Acta Geophysica, 2025, 73(1): 143-161.

参考的文章、视频、官网如下：

InSAR_ISCE软件数据处理流程之：利用哨兵数据获取2015年尼泊尔地震同震形变场

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

isce_github官网