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
本站声明： 1. iCode9 技术分享网（下文简称本站）提供的所有内容，仅供技术学习、探讨和分享； 2. 关于本站的所有留言、评论、转载及引用，纯属内容发起人的个人观点，与本站观点和立场无关； 3. 关于本站的所有言论和文字，纯属内容发起人的个人观点，与本站观点和立场无关； 4. 本站文章均是网友提供，不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属；如您发现该文章侵犯了您的权益，可联系我们第一时间进行删除； 5. 本站为非盈利性的个人网站，所有内容不会用来进行牟利，也不会利用任何形式的广告来间接获益，纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。