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 ..
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
生成文件夹run_files 、configs 文件 SAFE_files.txt
4.给run_files赋予权限
sudo chmod 777 -R run_files/
cd run_files/
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
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
######## 如果做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*** 文件夹
SBAS 处理
一、在INSAR_20*** 文件夹下打开终端
mt_extract_info_isce
二、在此文件夹下打开matlab，在matlab中输入ps_load_info
ps_load_info
生成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)
RHO_MIN:最小相干性 （0.50）
DDIFF_MAX: 总去相关的天数（1500）
BDIFF_MAX:以m为单位的相关基线（1070）
四、在INSAR_20*** 文件夹下打开终端
make_small_baselines_isce
五、做多视，在INSAR_20*** 文件夹中放入多视文件 prep_ml_isce.m 然后在matlab中敲入下面的命令
# mt_prep_isce脚本下载地址`https://pan.baidu.com/s/1KWBl18m_xp2xFA8M-Pnehw?pwd=cmih`
prep_ml_isce('401') #此处还存在疑问
修改mt_prep_isce的161行代码，加入SMALL_BASELINES(不做多视就不需要)
七、cd MultiProc文件夹，打开终端，输入
mt_prep_isce 0.6 3 3 50 200
八、在MultiProc文件夹中打开matlab，做sbas
setparm('small_baseline_flag','y');
setparm('unwrap_method','3D_QUICK');
stamps(1,8);
