ISCE2.3+StaMPS4.1处理哨兵数据(三)——PS处理

首先根据上篇博客继续进行永久散射体处理（ps处理）哨兵数据预处理
ps处理流程如下：
1、在INSAR_master_date文件下运行mt_prep_isce 0.4 3 2 50 200命令。
按照设定的幅度离差阈值、距离向方位向分块和重叠范围等参数进行处理，选定初始PS候选点主要作用是提取所需信息，幅度校正，建立分块，选取初始PS候选点。生成的PATCH文件夹内容如下：

2、修改参数
PS处理的一些参数是默认的，如果需要修改则可以运行下列命令：

>>matlab
>>getparm
>>setparm(‘param_name’,param_value)

param_name参数名称，param_value参数值，正式进行PS处理
3、第一步

>>stamps(1,1)

把数据加载到matlab工作区
4、第二步

>>stamps(2，2)

估计相位噪声，这是一个迭代过程，用于估计每个干涉图中每个候选像素的相位噪声值，计算每个候选点的 值。
5、第三步

>>stamps(3，3)

PS点的选择，基于像素的噪声特性来进行PS点的选择，初选点方法：幅度离差指数阈值法，并且估计随机像素点（非PS点）的百分占比。
6、第四步

>>stamps(4，4)

去除那些只在一些干涉图中相位稳定的PS点及其受临近PS点影响而表示为PS点特征的像元。同时根据噪声标准差去除噪声点。将选择的新的PS点存储到新的工作空间

7、第五步

>>stamps(5，5)

相位校正，主要是DEM引起的相位误差。
8、第六步

>>stamps(6，6)

对选取的PS像元进行相位解缠
9、第七步

>>stamps(7，7)

空间失相关的视角(SULA)误差在步骤3中计算，并在步骤5中去除。在步骤7中，计算空间相关的视角(SCLA)误差，这几乎完全是由于空间相关的数字高程模型误差(这包括数字高程模型本身的误差，以及数字高程模型到雷达坐标的不正确映射)。同时估计主大气和轨道误差(AOE)相位
10、第八步

>>stamps(8，8)

11、生成Google地球的kml文件

ps_plot(‘v-d’,-1)
load ps_plot_v-d ph_disp
ps_gescatter(‘project_velo.kml’,ph_disp,10,0.4)

其他的一些可视化可在ps_plot的matlab代码中查看并进行可视化。
注意注意！！！后期写小基线集的处理流程以及将ps与小基线结合处理的过程，感兴趣的朋友可以评论进行讨论，有问题可指出！