在 SARscape PS(永久散射体)处理期间消除速度跳跃的策略
-
示例
使用 PS 处理SAR 数据时,可能会在地理编码步骤后观察到 PS 组之间的速度跳跃。如下图所示,其中蓝色 PS 组位于绿色 PS 组中间。
这种跳跃在地理编码步骤之后很明显,但也可以在第一次反演后的结果中看到。
-
造成了这种影响的原因
这种效应与 PS 过程有关,如果感兴趣区域很大(默认情况下超过 25 平方公里),则在第一次反演期间会将整个感兴趣区域切分为多块子区域。
可以在第一次反演的参数中检查此阈值:该参数名为“Subarea for Single Reference Point”。其默认值为 25 平方公里,即25平方公里一个参考点。
在第一次反演步骤中,每个分区会自动设置一个参考点。该 GCP 点通常对应于给定子区域中的最高相干性和最低速度点。然后,每个分区中的 PS 速度将根据其自己的参考点进行计算,这可能会导致最终输出速度的跳跃。
-
如何可视化分区和参考点
可以使用以下方法显示每个分区的参考点
- 在 ENVI 中加载 Velocity_first 文件。该文件位于PS进程的输出\work\work_first_inversion\目录下
- 加载 SubAreas.shp 文件以显示第一次反演期间定义的子区域。该文件位于PS进程的输出\work\work_first_inversion\目录下
- 加载 Ref_GCP.shp 文件以显示每个分区的参考 GCP。该文件位于PS进程的输出\work\work_first_inversion\目录下
- 每个分区的ID在SubAreas.shp文件的属性表中定义,它对应于多边形ID
-
消除速度跳跃的策略
选项 1: 通过在第一次反演步骤中增加“Subarea for Single Reference Point”的参数值来增大创建子区域的最小尺寸。这将减少用于计算的分区数量,从而一定程度限制速度跳跃现象。
然而,它可能有副作用:
- 检测到的 PS 数量减少;
- 降低大气校正效率,这将导致远离参考点的区域的速度估计出现错误;
选项2: 在子区域之间的重叠区域中选择新的参考点,可以使用以下工具来完成:
/SARscape/Interferometric Stacking/Stacking Tools/PS Compute Specific Sub-Area
选择分区ID后,在分区重叠的区域中选择一个新的GCP点。
然后,需要重新运行第一次反演和第二次反演,以及之后的地理编码步骤。在第一次反演期间需要执行合并步骤,并且需要将第一次反演中"Other Parameters"的"SubArea Merging Only"参数设置为真。
最后,可以通过显示速度结果来检查 GCP 更新是否有效。