前言
国产系列卫星的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只有12M。
链接: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