高分遥感数据预处理
以高分一号为例
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工具:
- 在Input Raster中输入待正射校正的影像;
- 在Input Reference Raster中输入基准影像(一般为Landsat8-OLI影像);
- 在Input DEM Raster中输入正射校正所需的地形数据,使用ENVI自带的DEM数据,文件路径为"D:\Program Files (x86)\Exelis\ENVI53\data\GMTED2010.jp2"
- 重采样方法选择最邻近法;
- 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:从实验室前辈中学习的