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自带有这个数据。打开方式如下图

在进行校正前，如下图所示，可以看到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. 校正的结果对比

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