SARscape处理ALOS数据实践

•  区域增长法(Region Growing)：选这种方法，则不要设置过高的相干性阈值（0.15-0.2是比较好的）以便留下足够的自由增长空间，相位突变部分在解缠后的图像上以解缠孤岛存在，这种方法降低了由相位突变引起的误差。
•  最小费用流(Minimum Cost Flow)：默认的解缠方法，当有大面积的低相干或是其他限制增长的因素而使解缠困难时，最小费用流算法可以取得比区域增长法更好的结果。这种方法采用正方形的格网，考虑了图像上所有的像元，对相干性小于阈值的像元做了掩膜处理。
• Delaunay MCF：和最小费用流法的不同在于，这种方法不是考虑了图像上所有的像元，而是仅考虑了相干性大于阈值的部分，而且不是用正方形的格网而是用了德罗尼三角形格网。只有对相干性高的部分进行解缠，不受低相干像元的影响。对于有大量相干性低的地物存在的时候，如影像上存在大量水体、浓密植被等，推荐使用该方法，最小化相位突变的影响。

 8.轨道精炼和重去平

当输入精确的轨道信息，为了矫正相位偏移，所以要进行轨道参数的修正，矫正的结果不会生成新的数据文件，而是将解缠图像头文件中的信息做了修正。

选择/SARscape/Interferometry/Phase Processing/4 - Refinement and Reflattening

用SARscape做干涉处理时（InSAR/DInSAR），在轨道精炼的一步，必须输入控制点文件（refinement GCP file），用于修正干涉相位和解缠后的相位。

点击按钮可启动控制点选择工具。在生成GCP工具的在文件选择面板，输入三个文件，分别是：_upha、_dem和_fint。

用于轨道精炼的控制点一般遵循以下几个原则：

• 不能位于高频的残余地形相位上
• 不能位于形变相位上，要远离形变区域，控制点的被认为是形变为0的点
• 不能位于解缠错误的相位跃变上（phase jump）
• 如果fint上存在规律的几何误差，如轨道条纹，GCP在满足上述条件的前提下，分布于整个图像，从而去除规律的相位误差

下面分不同情况来说明控制点选择的建议。

 如果是InSAR测量DEM，在整个图像中选择几个控制点用于去除系统的几何误差（轨道误差），GCP要远离陡峭的地形区域和残余地形相位，当地形起伏大的山区，最好是选择山谷底部的平地区域。

 如果DInSAR测量地表形变时，在没有轨道误差的情况下，远离形变区域的位置选择一个GCP点就可以。

 做InSAR测量地形或者DInSAR测量地表形变时，如果存在轨道误差的情况下，需要选择多个GCP点，作为稳定的参考点，程序从这些点中计算出误差相位从而去除。

 不论什么情况，在选择GCP之前，确定处理区域没有解缠错误的相位，也就是相位跃变的情况。

 有时候为了比较结果，会对同一区域用不同的观测几何的数据对进行处理，即升轨数据对和降轨数据对。这种情况下，GCP最好使用同一组，这样可以减少两个图像之间的相位误差，方法是：对两个数据对的去平后干涉图_fint进行地理编码，在地理坐标的图上选择GCP点。

 最后，有两个注意事项：（1）如果要选择子区域的话，最好选择比研究区大一些的区域，以保证除了形变区域以外还存在稳定区域用于选择控制点。（2）轨道精炼的方法有轨道优化和多项式优化两种。如果选择轨道优化的方法，至少需要7个控制点，否则，会自动改为多项式优化，控制点最少为多项式系数（Residual Phase Poly Degree），否则，多项式系数会自动减少。

 9.相位高程转换——生成 DEM

这一步是将经过绝对校准和解缠的实际相位，结合合成相位，转换为 DEM 并进行地理编码。生成进行了地理编码的 DEM 文件、相干图像，图像精度图像和分辨率图像。
选择/SARscape/Interferometry/Phase Processing/5A - Phase to Height Conversion and Geocoding

10.相位转形变以及地理编码

经过绝对校准和解缠的相位，结合合成相位，转换为形变数据以及地理编码到制图坐标系统，默认得到的是 LOS 方向的形变。
选择/SARscape/Interferometry/Phase Processing/5B - Phase to Displacement Conversion and Geocoding