前言
国产系列卫星的L1级数据因为没有进行集合校正，所以相对于实际位置有偏移，这种情况在GF-1的WFV传感器数据表现得尤为强烈，但是同时在数据下载的时候，会带有RPC参数文件，可以获取比较准确的位置，相对来说，国产卫星的RPC参数比较准，但是这里介绍一下，基于RPC参数，再基于landsat影像进行自动匹配同名点来进行正射校正，效果更佳。landsat影像提供的是L1T数据，经过全球控制点的几何校正，相对来说位置比较准，尤其是，最近USGS正在生产L1TP数据，位置更加准确，我们来看一下他的介绍。Tier 1 (T1) – Contains the highest quality Level-1 Precision Terrain (L1TP) data considered suitable for time-series analysis. The georegistration is consistent and within prescribed tolerances [<12m root mean square error (RMSE)].反正意思就是说误差很低.RMSE只有１２Ｍ。
链接：https://landsat.usgs.gov/landsat-processing-details
然后再看一下GF-4的偏移，相对来说很小。
我们再来看下GF-1的WFV数据：
简直不忍直视！！！
步骤
此次以截取的一小段GF-4数据演示，软件ENVI5.3.1，在这里下载：https://www.ixxin.cn/software.html
加载影像数据
值得注意的是，landsat影像要完全覆盖你要校正的数据。
打开RPC Orthor......Using REference ....
打开这个工具后，按要求输入参数
Input Raster：选择待校正GF4影像
Input Reference Raster：选择拼接好的Landsat8-OLI全色影像
Input DEM Raster：选择DEM数据（GDEM V2 30M）
DEM Is Height Above Ellipsoid：选择No
Requested Number of GCPs：输入100
Image Resampling：选择三次卷积
Grid Spacing：保持默认，10
Output Raster：设置输出路径和文件名
Output GCPs：设置自动采集的控制点文件输出路径和文件名
Display result：勾选，校正结束后加载校正结果
前后对比
与landsat8叠合情况：
使用生成的控制点再次纠正
上一步除了正射校正外还生成了一些控制点，我们把这些控制点导入到正射校正模块，进行有控制点的正射校正。
在左侧控制点列表框中，选中某一控制点，当前视图窗口会自动跳转到该控制点所在位置，同时在右侧GCP Properties面板中，可以看到该控制点的详细信息，包括控制点名称、经纬度坐标（X，Y，Z）、像素坐标（X，Y）、控制点是否参与RPC模型的调整（Adjustment参与，Independent不参与）以及该控制点在X和Y方向的误差。图5为选中GCP1后在视图中的位置，可以借助工具栏上放大和移动按钮进一步查看。
另外，在RPC Refinement面板右下角，经过一番调整后（删除误差较大的点），可以看到在当前控制点参与调整后平面误差为15.33，单位为米，基本为0.3个像元误差，达到预期，点击finish。
结果，高分四：
我们再看一下偏移很大的高分一号也回到的正确的位置：
参考文章
http://blog.sina.com.cn/s/blog_764b1e9d0102ww9y.html
后语
以参考影像进行高精度正射校正，达到了比较好的效果，但是Erdas中有根据立体像对进行正射校正，效果也很好，有机会写一下。
参考网址;https://landsat.usgs.gov/landsat-collections
打赏直接扫描二维码，打赏名单请看：https://www.ixxin.cn/about.html
