StaMPS软件操作流程

注意：冒号后面都是解释，如果步骤。

没说更改目录就一直是前一步的目录。

每个命令都有一个log对应
1成像处理（ROI_PAC）
1.1 SAR影像原始数据准备
建立SLC文件夹，将SAR原数据拷贝到SLC目录里，为每景影像按照获取时间分别建立yyyymmdd文件夹，每景SAR影像原始数据对应的数据和头文件分别命名为IMAGERYyymmdd和SARLEADERyyyymmdd。

1.2 主影像成像(对应的是0级)
按照主影像选取原则选取主影像，进入到主影像文件夹master_date，对于ERS-1/2数据，运行step_slc_ers命令成像。

step_slc_ers主要通过调用ROI_PAC软件中的make_raw.pl, roi_prep.pl, roi.pl命令，及其StaMPS本身的roipac2doris命令来完成。

1.3 选择主影像兴趣区域重新成像I yao
观测主影像成像图像，确定研究区域，按照研究区域边界扩展1000像元的范围来编辑roi.proc。

重新运行step_slc_ers，按照兴趣区域重新成像。

step_slc_ers命令用法同1.2.//粗裁。

直接在SLC文件夹下生成有这个文件，在里面修改，别忘记去掉#号。

1.4 建立主影像
精裁，观察主影像成像图像，编辑master_crop.in文件，再次确定裁剪区域。

运行step_master_setup命令建立主影像。

//从ROI_PAC_SCR中复制master_crop.in到主影像成像的文件夹中。

1.5 辅影像成像
在SLC目录下，运行make_slcs_ers命令，对所有辅影像按照兴趣范围成像。

运行过程为逐个进行辅影像目录，执行step_slc_ers命令成像。

step_slc_ers用法同1.2.
//注意其他辅影像成像的时候是按照粗裁的范围裁剪的，他不进行精裁，由于主影像又进行了精裁，一般辅影像比较大。

对应1级产品不用成像，直接仿照说明书，读取即可
1建立连接并读取主影像
link_slcs_路径（初始数据的）
cd master_date（即是主影像日期文件夹）//进入主影像日期文件夹
step read whole XXX (where XXX is ‘ERS’, ‘Envisat’, ‘RSAT’, or ‘TSX’)第一次对主影像读取，全部读取，不裁剪
2裁剪主影像
cp $MY SCR/master_crop.in 或者从安装文件里考
并编辑
3裁剪后主影像再次读取
step_master_read
4读取各个辅影像
返回到SLC目录
make_read
2 差分干涉处理（DORIS）
2.1 提取精密轨道信息（冒号后面的是此步骤的解释说明，不是让你
这么操作）
在insarmaster文件夹下
step_master_orbit_ODR：主影像精密轨道信息提取。

//没有精轨不运行这步，直接下步make_orbits：为每幅辅影像建立相应的文件夹，在每个文件夹内执行step_orbit命令，读取相应辅影像精密轨道信息。

2.2 主辅影像配准
2.2.1 主辅影像粗配准
make_coarse：进入每幅辅影像相应的文件夹，在每个文件夹内执行step_coarse命令，（coarse.dorisin）进行主辅影像一级轨道配准和二级像素级配准。

2.2.2 主辅影像精配准
make_coreg：精配准，精配准过程中所有基线距小于设定的bperp_max值的直接与主影像配准，基线距比较大的和离它最近的三个辅影像配准。

然后利用Hooper提出的配准方法，利用中间过渡的方法实现所有主辅影像的配准。

make_coreg命令运行Doris软件中的FINE和COREGPM步骤，即三级子像素精配准和估计配
2.3 主影像强度图模拟与地形相位计算(平坦地区地形可不执行)
cd INSAR_master_date//在INSAR_master_date文件夹，修改timing.dorisin文件参数
step_master_timing：模拟主影像强度图，计算DEM与主影像的Timing Error。

make_dems：进入到每个辅影像目录，执行step_dem命令，计算每个干涉对的地形相位。

2.4 辅影像重采样
make_resample：对辅影像按照精配准算出的多项式来重采样。

进入到每个辅影像目录，执行step_resample。

(resample.dorisin,自动加入RS_DBOW 1000 10999 1500 3019)
2.5差分干涉处理
make_ifgs：主辅影像进行干涉处理，去除平地相位和地形相位。

进入到每个辅影像目录，运行step_ifg命令。

（integram.dorisin）（make_ifgs_nodem对于平坦地区不用去地形）
step_ifg运行Doris软件中的SUBTRREFDEM、COMPREFPHA和SUBTRREFPHA步骤，即生成干涉图，利用轨道来计算平地相位，减去平地相位，减去由粗DEM反算得到的地形相位：
2.6 地理编码
step_geo：从雷达坐标系转换到地图坐标系//只在INSAR_masterdate文件夹下任意一个辅文件（就是日期那个文件夹）夹下运行一次就行。

3 PS处理
3.1准备
注意小基线不是在此步初选点，是在确定小基线对后再初选点，别弄错了
//命令是mt_prep 0.4 3 2 50 200还是在INSAR master date目录下
mt_prep：按照设定的幅度离差阈值、距离向方位向分块和重叠范围等参数进行处理，选定初始PS候选点。

3.2 stamps处理
//在INSAR master date下先启动matlab
getparm(parmname),无参数时显示所有参数，默认参数从ps_parms_default.m提取；setparm(‘parmname’,value),修改某个参数值。

参数结果保存在parms.mat文件里。

//在INSAR master date下，运行stamps
stamps（start_step，end_step，patches_flag，est_gamma_parm）
Patches_flag 默认为‘y’，如果设置为‘n’把所有数据作为一块处理。

Est_gamma_parm是一个可选参数passed to ps_est_gamma_quick。

选y则利用ps_est_gamma_quick（est_gamma_parm），n利用ps_est_gamma（est_gamma_parm）。

stamps(3,3)：调用ps_select.m，基于计算出的每个候选点的γ值，同时考虑幅度离差选取
stamps(4,4):调用ps_weed.m，去除那些只在一些干涉图中相位稳定的PS点及其受临近PS点影响而表示为PS点特征的像元。

同时根据噪声标准差去除噪声点。

stamps(5,5):调用ps_correct_phase.m和ps_calc_ifg_std.m（计算每一个干涉图噪声标准差，作为权重），对选取的PS像元的缠绕相位去除spatiall-uncorrelated look angel (DEM) error.还要调用ps_merge_patches
stamps(7,7): 调用ps_calc_scla.m、ps_smooth_scla.m(减去SCLA后重新解缠时用)，估计
spatially-correlated look angle error ，同时估计出主影像大气影响相位和轨道误差相位（AOE ）
最后输出显示
ps_plot(‘v’，4) 4 小基线处理
如果没有运行前面的PS 处理，则需要先在INSAR
master date 目录下读取基线信息加
载到matlab 工作空间，命令如下
mt_extract_info
matlab
>>ps_load_info
4.1 确定小基线对
在matlab 里sb_find:寻找小基线干涉对//sb_find(rho_min,ddiff_max,bdiff_max)运行后显示基线图。

也可以使用命令plot_sb_baselines 调出基线图
//终端在INSAR master date 目录下运行
make_small_baselines ：
方位向、距离向滤波，按照选择的小基线对，生成差分干涉图。

运行Doris 软件中的m_filtazi ，s_filtazi ，FILTRANGE 与INTERFERO 步骤，及其Doris 的cpxsum 与cpxfiddle 命令。

并在INSAR master date 下生成SMALL BASELINES 文件夹
将原来配准后的文件slave_res.slc 考入到相应的干涉对组成的文件夹内分别进行主辅影像方位向滤波与距离向滤波，生成干涉图。

利用PS 处理中每个干涉对中的干涉图和去除地
形相位的差分干涉图，结合现有的干涉图，用Doris命令cpxsum处理，生成新的差分干涉
4.2 mt_prep
mt_prep：按照设定的幅度差分离差阈值、距离向方位向分块和重叠范围等参数进行处理，
选定初始高相干候选点。

//cd SMALL BASELINES进入这个文件夹
mt_prep 0.6 3 2 50 200
4.3 SB处理
运行过程同PS中stamps的使用，可参考上面。

先启动matlab还是在SMALL BASELINES下
getparm(parmname),无参数时显示所有参数，默认参数从ps_parms_default.m提取；setparm(‘parmname’,val ue),修改某个参数值。

参数结果保存在parms.mat文件里。

Stamps（start_step，end_step，patches_flag，est_gamma_parm）
Patches_flag 默认为‘y’，如果设置为‘n’把所有数据作为一块处理。

Est_gamma_parm是一个可选参数passed to ps_est_gamma_quick。

选y则利用ps_est_gamma_quick（est_gamma_parm），n利用ps_est_gamma（est_gamma_parm）。

stamps(3,3)：调用ps_select.m，基于计算出的每个候选点的γ值，同时考虑幅度离差选取高相干像元。

stamps(4,4):调用ps_weed.m，去除那些只在一些干涉图中相位稳定的高相干点及其受临近
stamps(5,5):调用ps_correct_phase.m和ps_calc_ifg_std.m，对选取的高相干像元的缠绕相位
stamps(6,6):调用ps_unwrap.m，sb_invert_uw.m, 对选取的高相干像元的缠绕相位进行相位
stamps(7,7)和stamps(8,8): 调用ps_calc_scla.m、ps_smooth_scla.m和ps_scn_filt.m，估计小
基线干涉对和PS处理干涉图中的spatially-correlated look angle error，同时估计出主影像大。