SARscape中导入外部GCP点用于轨道精炼

在 SAR 处理时，有时会加入 GCP 点文件， SAR 处理中用到的控制点分为两类：用于校正地理位置的几何控制点（ Geometry GCP ）和用于轨道精炼的控制点（ Refinement GCP ）。

几何控制点的详情请参考： http://blog.sina.com.cn/s/blog_764b1e9d0102xd72.html 。

轨道精炼控制点在以下处理步骤使用：

l   InSAR/DInSAR 处理工作流： /SARscape/Interferometry/Phase Processing/4 - Refinement and Re-flattening

l   立体 SAR 处理工作流： /SARscape/Interferometry/Stereo-Radargrammetry/2 - Shift Refinement and Re-flattening

l   MAI 处理工作流： /SARscape/Interferometry/MAI Processing/2 - MAI Refinement and Re-flattening

l   SBAS 处理工作流： /SARscape/Interferometric Stacking/SBAS/3 - Refinement and Re-Flattening

轨道精炼 GCP 点的选择常规方法在此不再赘述，选择依据等可参考 http://blog.sina.com.cn/s/blog_764b1e9d0102xdfj.html 。

有些情况下，用户希望能从已有的底图上（光学底图或 GoogleEarth 上）选择 GCP ，比如做地表形变测量时，用户基于对工程区已有的经验，在光学底图上能确定稳定的点的位置，那么在光学底图上找到的 GCP 点如何导入 SARscape 中用于轨道精炼呢？本文以 SBAS 工作流，从 Google Earth 上选择 GCP 为例，介绍从已有底图上选的 GCP 点用于轨道精炼方法。

以下的操作是在 SBAS 工作流已经完成了连接图生成、干涉工作流之后进行：

第一步，打开工程区的 kml 文件，在 Google Earth 上，使用 工具，在工程区内定位若干个 GCP 点（点要位于没有形变的位置、最好是人工地物上，相干性高），将位置另存为 kmz 文件。

图 在 GoogleEarth 上选择 GCP 点

第二步：在 ArcGIS 中，将 kmz 文件转为 shapefile 文件，工具是 kmz to layer ，再将 layer 另存为 shapefile 文件。

注：此方法仅供参考，把 kmz 转成 shapefile 即可。

图 在 ArcGIS 中将 kmz 文件转为 shp

注：以上两步是在 Google Earth 上选 GCP 点的操作，如果用户要在地理坐标的光学底图上选 GCP 点，直接在 ENVI 或 ArcGIS 中生成 shape 矢量点图层即可。 如果有已知的GCP点的xy坐标，直接用点坐标文件生成shape文件即可。

第三步：将上一步得到的地理坐标系的 shp ，在 SARscape 转为 SAR 坐标系的 shapefile 文件。打开 /SARscape/Basic/Intensity Processing/Geocoding/Map to SAR Shape Conversion 工具。

l   Input Files ：

-      Input File ：选择上一步准备好的地理坐标系的 shapefile 文件；

-      Input Reference File ：选择 SAR 坐标系的参考数据，（在接下来的轨道精炼处理中， GCP 点要在哪个图像上使用，就选哪个），此处选择干涉工作流中生成的某一对像对的 fint 文件作为参考数据；

图 输入文件面板

l   Optional Files ：输入用于几何校正的 GCP 文件，此处不输入。

l   DEM/Cartographic System ：输入参考 DEM 或者选择坐标系，作为地理坐标系的参考，此处输入参考 DEM 文件。

图 输入参考 DEM 文件

图 输出斜距的 shapefile 文件

第四步：在轨道精炼一步，选择 GCP 的界面上，导入上一步生成的斜距 shapefile 文件，即可进行轨道精炼。

    在轨道精炼数据输入面板上，输入工程文件，点击 Create GCP ，选择 input file 和参考 DEM 文件，点击 Next ，

图 点击 Create GCP

在 GCP 生成面板上上，点击 Load GCPs ，选择上一步生成的 SAR 坐标系的 shapefile 文件，导入 GCP 点，可检查点在参考图像上的位置，如果有不在图像范围的 GCP 点，删除即可。导入的 GCP 点号有乱码的情况，忽略即可。也可以手动再添加新的 GCP 点。

图 点击 Load GCPs 导入 SAR 坐标系的 GCP 点 shapefile 文件

切换到 Cartographic System ，选择地理坐标系： GEO-GLOBAL-GEO-WGS84

图 选择地理坐标系

切换到 Export 面板，设置输出 xml 点文件的路径，点击 Finish 。

图 SARscape 中控制点文件生成

在面板上确认输入生成的控制点文件之后，进行轨道精炼操作即可。

图 轨道精炼数据输入面板