ISCE2StaMPS全流程

1.准备文件夹

准备四个文件夹：AxuDir 、DEM 、Orbits 、Process、 SLC

2.准备 dem

cd DEM
dem.py -a stitch -b 38 39 116 118 -r -s 1 -c

rm demLat*.dem demLat*dem.xml demLat*.dem.vrt
cd ..
1
2
3
4
5

3.预处理

#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


cd Process
stackSentinel.py -s ../SLC/ -d ../DEM/demLat_N38_N39_Lon_E116_E118.dem.wgs84 -a ../AuxDir/ -o ../Orbits -b '38.62 39.23 117.1 117.6' -W slc -z 5 -r 20
1
2
3
4
5
6
7

生成文件夹run_files 、configs 文件 SAFE_files.txt

4.给run_files赋予权限

sudo chmod 777 -R run_files/

cd run_files/
1
2
3

5.对每个文件依次进行处理(如果第一步报错，可能需要把DEM复制到该文件夹下)

./run_01_unpack_topo_reference
./run_02_unpack_secondary_slc
./run_03_average_baseline
./run_04_extract_burst_overlaps
./run_05_overlap_geo2rdr
./run_06_overlap_resample
./run_07_pairs_misreg
./run_08_timeseries_misreg
./run_09_fullBurst_geo2rdr
./run_10_fullBurst_resample
./run_11_extract_stack_valid_region
./run_12_merge_reference_secondary_slc
./run_13_grid_baseline
1
2
3
4
5
6
7
8
9
10
11
12
13

6.最终生成的所需文件都在merged文件夹下

包括baselines文件夹 geom_reference文件夹 SLC文件夹

7.生成input_file文本文件内容

source_data slc_stack
slc_stack_path /media/estar/es/INSAR/ISCE/Process/merged/SLC
slc_stack_reference 20170320
slc_stack_geom_path /media/estar/es/INSAR/ISCE/Process/merged/geom_reference
slc_stack_baseline_path /media/estar/es/INSAR/ISCE/Process/merged/baselines
range_looks 40
azimuth_looks 10
aspect_ratio 4
lambda 0.056
slc_suffix .full
geom_suffix .full
1
2
3
4
5
6
7
8
9
10
11

######## 如果做sbas，在isce2/contrib/timeseries/prepStackToStaMPS/bin/make_single_reference_stack_isce中第二行 setgenerate_ifgs=n,将y改成n

在主文件夹 （与SLC同级别）下运行

make_single_reference_stack_isce    # 这步操作会将input_file变为 INSAR_20*** 文件夹
1

SBAS 处理

一、在INSAR_20*** 文件夹下打开终端

mt_extract_info_isce
1

二、在此文件夹下打开matlab，在matlab中输入ps_load_info

ps_load_info
1

生成loading info into matlab
SERTPARM: heading = -10.3323
SERTPARM: lambda = 0.056

三、输入sb_find(0.4, 48, 150) 0.4代表相干性， 48代表最大时间基线48天 150代表最大空间基线150m

sb_find(0.4, 48, 150) 
1

RHO_MIN:最小相干性 （0.50）
DDIFF_MAX: 总去相关的天数（1500）
BDIFF_MAX:以m为单位的相关基线（1070）

四、在INSAR_20*** 文件夹下打开终端

make_small_baselines_isce
1

五、做多视，在INSAR_20*** 文件夹中放入多视文件 prep_ml_isce.m 然后在matlab中敲入下面的命令

# mt_prep_isce脚本下载地址`https://pan.baidu.com/s/1KWBl18m_xp2xFA8M-Pnehw?pwd=cmih`
1

prep_ml_isce('401')  #此处还存在疑问
1

修改mt_prep_isce的161行代码，加入SMALL_BASELINES(不做多视就不需要)

七、cd MultiProc文件夹，打开终端，输入

mt_prep_isce 0.6 3 3 50 200
1

八、在MultiProc文件夹中打开matlab，做sbas

setparm('small_baseline_flag','y');
setparm('unwrap_method','3D_QUICK');

stamps(1,8);
1
2
3
4