高分遥感数据预处理

高分遥感数据预处理

以高分一号为例

1. 安装GF-1插件

下载网址：http://www.enviidl.com/appstore/
安装：放在目录C:\Program Files\Exelis\ENVI53\extensions\下面

2. 以GF-1插件加载文件
在ENVI5.3中使用GF1专用插件（Open GF1）打开xml后缀文件

3. 正射校正（附带几何校正）
在Toolbox中打开RPC Orthorectification Using Reference Image工具：

1. 在Input Raster中输入待正射校正的影像；
2. 在Input Reference Raster中输入基准影像（一般为Landsat8-OLI影像）；
3. 在Input DEM Raster中输入正射校正所需的地形数据，使用ENVI自带的DEM数据，文件路径为"D:\Program Files (x86)\Exelis\ENVI53\data\GMTED2010.jp2"
4. 重采样方法选择最邻近法；
5. Output Raster为输出路径，选择指定位置，输出的GCPs也可以指定路径。
   
   4. 研究区裁剪
   在工具栏中选择File-Save as-Save as…(ENVI, NITF, TIFF)，在File Selection中选择正射校正完的影像,并点击下方Spatial Subset按钮,进行研究区裁剪，拉选红框至研究区，点击ok后选择输出路径。
   
   5. 辐射定标
   在Toolbox 中选择Radiometric Calibration进行辐射定标，在File Selection中选择裁剪过后的数据，点击ok。
   根据之后所进行大气校正的方法选择定标后的物理量。
   （1）若是需要进行6S大气校正，则需要将Calibration Type选择为Reflectance；
   （2）若是需要进行FLAASH大气校正，则需要点击下方按钮Apply FLAASH Settings；
   （3）若是需要进行QUAC大气校正，则保持默认设置就好。
   选择完定标的方式后，确定输出路径，并点击ok。
   
   6. 大气校正
   （1）待辐射定标完毕，现进行大气校正，以FLAASH为例。
   在Toolbox中打开FLAASH Atmospheric Correction工具，并进行一系列的参数输入。（若以中国为研究区时间要-8）
   

在Multispectral Settings中的Filter Function File中加入光谱响应函数（sli文件，插件安装后会自动载入高分1的sli文件）

设置完毕后，点击Save…将所有的参数进行保存，最后点击Apply运行。
（2）以QUAC为例，在Toolbox中输入QUAC，打开选择进行QUAC的影像，并选择Genric方式进行大气校正。

补充：有全色影像的需要进行融合

7. 转换为遥感反射率

此处需要使用band math工具，在band math中输入公式
1.0 * b1 / (!PI * 10000)
输入完后点击Add to list然后点击ok

点击Map Variable to Input File，选择大气校正后影像

8. 参数反演

:同样需要在BandMath中完成，以SSD反演为例，在Bandmath中输入公式:

9. 使用边界shp对反演得到的影像进行裁剪

将shp文件拖入ENVI窗口，或者在File-Open中打开shp文件。在Toolbox中选择Subset Data From ROIs，选择要裁剪的反演影像。

Mask pixels outside of ROI选择Yes，并在Mask Background Values输入nan输入输出文件路径。
PS：从实验室前辈中学习的