ENVI下MODIS L1B数据大气校正（FLAASH&QUAC）

MODIS数据是著名的大尺度低空间分辨率数据源，一般用于大范围监测、宏观信息的快速提取。在MODIS数据预处理时，最常做的是几何校正。对MODIS数据进行几何校正的工具很多，诸如MCTK工具，或者直接使用ENVI提供的Georeference MODIS、Reproject GLT with Bowtie Correction等工具进行。在一些定量、半定量的监测中，诸如参数反演、植被指数计算时，需要做大气校正。

本文以一景L1B 级1km 的MODIS数据为例，介绍在ENVI5.6.2（其他版本类似）下MODIS L1B（即MOD02产品）数据预处理的方法和操作步骤，重点介绍MODIS L1B数据FLAASH大气校正的实现。

1 几何校正和定标处理

由于MODIS L1B数据本身自带几何定位信息（GLT），所以在进行大气校正之前，建议大家先做几何校正。几何校正建议采用ENVI扩展工具：MCTK进行，因为该工具不仅可以在几何校正的同时对数据进行定标处理，而且校正完毕后其结果图像会自动添加中心波长信息，便于下一步的FLAASH大气校正处理。

推荐使用 ENVI App Store 进行 ENVI 扩展工具的安装与管理。

App Store 下载地址：envi.geoscene.cn/appstore

启动ENVI。在Toolbox中，选择Extensions > modis_conversion_toolkit，弹出MODIS Conversion Toolkit面板。该面板大致可以分为左中右三部分，左侧部分用来设置输入输出信息，中间部分用来选择数据处理类型，右侧部分用来设置重采样方法和投影信息等。关于该工具的详细说明可以参考MCTK用户手册，本例中需要设置的参数项说明如下：

Input HDF：选择需要处理的MODIS数据（HDF格式）；
Select Output Type：选择输出类型，选择第二项Projected：Triangulation-based georeferencing；
Output Path：选择处理结果存储路径；
Rootname for output：设置输出文件前缀；
Processing Options：选择处理的数据类型，因为下面要进行FLAASH大气校正，所以这里我们选择Radiance/Emissivity；
Select Resampling Method：选择重采样方法，Nearest Neighbor即可；
Select Output Map Projection：选择输出投影类型，一般选择Geographic Lat/Lon经纬度投影；
Background Value to Use：背景值设置，设置为0，默认NaN；
其他参数默认。

注：本文用的是MCTK V2.1.11版本，其他参数有所区别。

具体参数设置如下图所示：

图1 MCTK工具参数设置面板

参数设置完毕后，点击Process进行处理。处理完毕后，数据自动加载到Data Manager中，如果默认没有自动加载到图层中，可以在Data Manager中，加载数据进行查看。

图2 经几何校正和定标处理后的MODIS数据

从上图可以看出，经几何校正后的图像具有投影坐标和中心波长信息，但其存储类型还是BSQ，而FLAASH大气校正工具对输入数据存储类型的要求是BIL或BIP。所以，在进行FLAASH大气校正之前，还需要进行存储类型的转换。

2 存储类型转换

在Toolbox中，选择Raster Management > Convert Interleave，弹出Convert File Input File对话框，选择上一步处理结果，点击OK；
在弹出的Convert File Parameters参数设置面板中， Interleave选择BIL。
选择输出文件名，单击OK执行。

3 FLAASH大气校正

ENVI中FLAASH大气校正工具对输入数据除了要求必须有中心波长信息（多光谱数据）、存储类型为BIL或BIP以及数据类型为整型（INT、UINT、LONG）或浮点型（FLOAT）之外，还要求输入数据是单位为μW/cm²*nm*sr的辐射亮度数据。但是，经过上述处理得到的辐射亮度数据单位为W/m²*μm*sr，二者相差10倍。所以，在下面选择输入数据时我们需要进行单位转换，这也就是下面Single scale factor需要输入10的原因。

注：如果不知道定标辐射亮度值的单位，可以查看像素值，一般μW/cm²*nm*sr的辐射亮度的数值在10以内。

数据准备无误后，就可以进行FLAASH大气校正。

在Toolbox中，选择Radiometric Correction > Atmospheric Correction Module > FLAASH Atmospheric Correction，弹出FLAASH Atmospheric Correction Model Input Parameters FLAASH大气校正参数设置面板。FLAASH大气校正参数设置面板可以大致分为基本参数设置、多光谱/高光谱参数设置以及高级参数设置三部分。

基本参数设置各项参数含义如下：

Input Radiance Image：选择辐射亮度数据，在弹出的Radiance Scale Factors中，选择Use single scale factor for all bands，输入10，点击OK；
Output Reflectance File：设置大气校正结果的输出路径和文件名，建议在原文件名的后面添加_flaash作为标识；
Output Directory for FLAASH Files：大气校正过程中生成其他文件的存储路径，默认在C盘用户临时文件夹下，如果C盘空间不足或没有读写权限，可以选择其他盘，负则会出现常见的102错误，关于该错误详细说明可参考博文：FLAASH 大气校正出现102错误解决方法。

图3 FLAASH大气校正辐射亮度数据单位转换

Sensor Type：modis；
Ground Elevation：地面平均高程（需要注意此处单位为km），可以通过Google Earth或对应区域DEM数据统计获得；
Pixel Size：选择传感器类型后自动获取；
Flight Date：数据获取时间，可以通过文件名或查看元数据信息获取；
Atmospheric Model：大气模型。参考帮助文档选择，具体选取时根据数据纬度和时间信息取交集。
Aerosol Model：气溶胶模型，根据实际情况选择；
Aerosol Retrieval：气溶胶反演方法，默认为2Band(K-T)；
Water Retrieval：是否进行水汽反演；
Water Column Multiplier：若Water Retrieval设置为Yes，则需要设置水汽反演波段，可选择940nm；
Initial Visibility：初始能见度设置，默认为40km，若选择了气溶胶反演方法，则这里的设置无效。

本例中基本参数设置面板具体参数设置如下图所示：

图4 FLAASH大气校正基本参数设置面板

在多光谱参数设置面板中，Water Retrieval水汽反演波段按照帮助中的表格进行设置。

本文选择的水汽吸收波段为940，水汽反演波段设置如下。

图5 水汽反演波段设置面板

K-T气溶胶反演波段选择Default > Over-Land Retrieval alternate(460:2100nm)，或者Over-Land Retrieval alternate(660:2100nm)，软件会自动选择Upper和Lower波段，MODIS都有对应的波段。其他参数默认

图6 K-T气溶胶反演参数设置面板

在高级参数设置面板中，需要将邻域校正Use Adjacency Correction关闭，否则会出现图8所示错误提示；其他参数可能需要修改的地方包括：Use Tiled Processing分块处理大小，根据内存情况设置；Output Reflectance Scale Factor反射率数据扩大的系数，默认扩大10000倍，如需要获取真实地表反射率数据，可在校正完毕后通过波段运算获取。