SARscape中的PS处理——修改第一次反演后子区域的结果

    SARscape中PS处理是流程化工具，包括连接图生成、干涉工作流、PS第一次反演、PS第二次反演、地理编码。每一步设置好参数，自动执行，不需要过多的人工干预就可以得到结果。如果得到的PS结果有明显的分块现象，如下图所示，这种情况需要检查和修改第一次反演的结果，本文介绍处理方法。

图 PS结果

第一次反演的结果存放在…\_PS_processing\first_inversion文件夹下，分别打开cc_first，Velocity_first、SubAreas.shp、Ref_GCP.shp文件，如下图所示：

图 查看第一次反演结果

图 PS第一次反演参数面板

PS第一步反演是分块进行的，默认的是25平方公里一个参考点（SubArea for Single Reference Point）。相邻分块重叠30%（SubArea Overlap）基于振幅离差指数自动选出一个或者多个PS参考点，计算出相位偏移。

    找到出现问题的那块区域，检查该区域的分块情况，找到该区域自动选出的参考点，发现参考点位于计算错误的位置。下面是强度图和cc图上该错误区域的位置，该区域包含了一块水体。

图 查看反演不正常的子区域

打开subarea和Ref_GCP，查看到ID为20、21的这两块cc明显和其他分块不连续，且和强度图对比，湖面左侧的区域不应该相干性很低，导致结果错误，在Ref_GCP中查看到两个子区域对应的参考点是同一个，可以看到这个参考点也位于结果不正常的区域。记下两个有问题的子区域的ID号：20、21。

图 查看子区域所用的参考点

打开像元定位工具，在这两个区域另外找一个像元，找点原则遵循PS点的原则，如屋顶、电杆、桥梁，通常是城市中的点，或其他人造地物，如温室、大坝、金属和混凝土的地物（如井盖周边，管道等）。除了这些之外，自然目标如裸露的岩石也可以是PS点。可以打开_PS_processing\interferogram_stacking路径下的mean和mu_sigma文件，结合cc文件，像元同时满足mean、musigma、cc值都比较高的位置，记录下所找到的点的行列号。本例中找到的点行列号是2173，1287，行列号即是像元在SAR图像的距离向和方位向的坐标。

图 寻找新的参考点

下面用新找到的这个参考点，对20、21子区域单独进行运算。

/SARscape/Interferometric Stacking/Stacking Tools/PS Compute Specific Sub-Area工具

图 子区域重新计算

运算完之后，重新运行PS反演第一步，速率范围和高程残差范围参数和之前做的时候保持一致。在Other Parameters下面，将SubArea Merging Only改为True，其他参数不变。运行第一次反演。

图 子区域重新计算后再做PS第一次反演

相同的方法，对21子区域也重新做，点可以用20子区域的同一个点，也可以重新选择一个，本例中用的是同一个点。

重新运算之后的第一次反演后的cc和vel如下图所示：

图 子区域重新计算之后的cc结果

图 子区域重新计算之后的vel结果

重新进行PS第一次反演之后，接着运行PS反演第二步和地理编码，rebuild all改为true，最终得到反演正常的PS结果。

图 子区域修改之后得到的最终的PS结果

标签：PS,结果,cc,反演,区域,参考点,后子
来源： https://www.cnblogs.com/enviidl/p/16275854.html