首先用isce进行预处理,首先在INSAR_master_date文件夹下执行。
(1)第一步,终端运行:
mt_extract_info
matlab
此时打开matlab,命令:
>>ps_load_info
读取基线信息加载到matlab工作空间。如果已经进行了PS处理,则可以省去第一步,直接在matlab中读取基线信息,即按照下面步骤进行。
(2)第二步,在matlab中
>>sb_find(1,1500,1070)
>```
其中第一个参数决定选择的小基线对的数量,数值越大选择的越少。每运行一次,选择的小基线对都会更新,注意不应该有孤立的小基线簇。也可以使用命令plot_sb_baselines调出基线图,之后处理会保存在SMALL_BASELINES文件夹下面。
(3)第三步,按照选择的小基线对,生成小基线的差分干涉图(终端在INSAR_master_date目录下运行)
make_small_baselines_isce
该指令将创建一个名为SMALL_BASELINES的新文件夹。
(4)第四步,在INSAR_master_date 文件夹的终端执行下列命令:
cd SMALL_BASELINES
mt_prep_isce 0.3 3 2 50 200
提取所需信息,幅度校正,建立分块,选取初始高相干候选点。(0.6这个值要比ps的高点。山区阈值选择0.3)其中3×2就是PATCH块的个数。
(5)同ps处理一样分别执行stamps1-5步骤。stamps(1,1)等步骤。注意matlab的工作空间是在SMALL_BASELINES文件夹下面。
(6)第六步
stamps(6,6)
一旦相位选择满意的话就可以使用第六步进行相位解缠,在解缠过程中,数据会暂时被重采样到网格中,网格的大小由参数“unwrap_grid_size”控制,并且至少应该与“merge_grid_size”一样大。继续将其设置为相同的值(200米),然后运行步骤6。
(7)第七步
stamps(7,7)
运行完成之后,便可得到估算的DEM误差相位、AOE误差相位、轨道误差相位和大气误差相位。
绘制结果同PS处理,可见上一篇博客。