ISCE2StaMPS全流程

软件的环境搭配可参考
ISCE安装
STAMPS安装

1.准备文件夹

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

2.准备 dem

切换到DEM文件夹下，运行dem.py下载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 ..

运行dem.py以后会生成以上文件，其中部分数据无用，可以用rm命令删除。

3.预处理

进入Process文件夹，其中38.62 39.23 117.1 117.6为研究区的范围，5，20为多视比

#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

运行stackSentinel.py 命令后会生成文件夹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

快捷方法

ls run* > list
cat run_list | while read rows
do
  echo "************************************** start run: $rows **************************************************"
  $rows
done 

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

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

7.生成input_file文本文件内容

在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变为input_file_20*** 并生成 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`

注意这个地方的40l要根据实际的多视比进行修改

prep_ml_isce('40l')  

修改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);

如有需要可通过主页的联系方式交流