GF1WFV数据预处理

GF1WFV数据,空间分辨率16m,在使用之前需要进行辐射校正,几何校正,下面详细说明一下预处理的过程。本文以中国山东青岛地区,2015年06月22日的数据为例。

工具
  1. ENVI 5.3.1

ENVI5.3.1带有自动正射校正的工具包,RPC Orthorectification Workflow。可以较为方便地对GF1WFV数据进行校正。

数据
  1. GF1WFV数据:GF1_WFV1_E119.6_N36.3_20150622_L1A0000876227
  2. Landsat8 全色波段数据:LC81200352013319_B8
  3. GMTED2010.jp2

Landsat8 全色波段数据空间分辨率15m,经过了精校正,因此可以用来作校正的参考图像。

GMTED2010.jp2为2010年全球高程数据。ENVI 5.3.1自带有这个数据。打开方式如下图

GMTED2010

在进行校正前,如下图所示,可以看到GF1宽幅数据失真是比较严重的。

校正前

1. 辐射校正

在ENVI中打开数据后。首先,进行大气校正。

工具箱
QUick Atmospheric Correction(QUAC)

Toolbox窗口,选择,可直接在搜索框中搜索。参数选择如图所示,选择unkown sensor

2. 生成控制点

自动生成校正的控制点。

工具箱
Gnerate GCPs from Reference Image

1). Toolbox窗口中选择对应工具箱
!

2). 设置参数

  • Input Raser: 待校正GF1WFV大气校正后的数据,20150622
  • Input Reference Raser: Landsat8 全色波段数据,LC81200352013319_B8
  • Input DEM Raser:GMTED2010.jp2

设置参数

3. 正射校正
工具箱
RPC Orthorectification Workflow
  1. Toolbox窗口中选择对应工具箱
  2. 输入待校正的数据20150622,和高程数据GMTED2010.jp2
  3. 打开之前保存的控制点文件
  4. 删除误差较大的点

一般来说,要保持到校正差小于0.5个像素,因此Horizontal Accuracy要小于8
这里有个技巧可以按F12切换到全图视图,点击工具箱的按钮,红色越深说明误差越大,鼠标可以选择控制点。
5. 设置参数

  1. Advanced->Output Pixel Size 输出的像素大小16
  2. Export->Output File 输出文件格式。我选择的是tiff便于后续matlab处理。

6.保存校正后的图像

保存图像

保存的时候可以根据需要,选择一部分感兴趣区域保存。由于原图较大,并且我需要多幅同一地区的多时相图像。保存方法是找到一些容易辨识的点,从这个点开始裁剪同样大小的图像。如下图所示:

裁剪保存图像
7. 校正的结果对比

可以看到校正结果是比较好的

猜你喜欢

转载自blog.csdn.net/qq_19531479/article/details/78978406