利用 MODIS 的 L1B 数据反演气溶胶的流程如下：该方案中的数据和算法介绍：1、MODIS 数据是采用的 MODIS L1B 1KM 数据应严格按照说明进行操作，例如数据是1km 的，数据的合成是反射率在上，发射率在下等，下面将详细介绍各个步骤2、这是在 ENVI 5.0 下做的北京市气溶胶反演，具体包括 MODIS 影像（HDF）的辐射校正、几何校正、云检测、气溶胶反演3、气溶胶反演算法采用经典的暗像元法（DDV）也叫浓密植被法，因此对于冬季反演的气溶胶效果不好4、气溶胶反演的查找表是利用 IDL 调用 6S 辐射模型得到的，采用的是一般的参数，因此 3-9 月期间都可以用这一个查找表进行气溶胶反演，也可以自己制定查找表5、七纬查找表，从左向右，依次为太阳天顶角，卫星天顶角，相对方位角，P0 、T、 S（辐射传输方程参数），最后一列为气溶胶光学厚度（AOD）6.感兴趣的可以提供用到的 modis 云检测工具和气溶胶反演工具以及详细的 pdf 文档MODIS 数据下载地址：http://ladsweb.nascom.nasa.gov/data/search.html一、MODIS 影像的辐射校正在 ENVI 5.0 中打开 MODIS 影像的 HDF 文件就已经做了辐射校正，打开 HDF 文件的方法是 File->Open As->EOS->MODIS,打开后在数据列表中可以看到三个文件，第一个是发射率 Emissive（band20-band36），第二个是辐射率 Radiance（band1-band26）,第三个是反射率 Reflectance（band1-band26），如图 1 所示。
图 1 ENVI5.0 打开 HDF 文件二、几何校正2.1 发射率文件的几何校正（1）Georeference MODIS 工具ENVI5.0 下有对特定传感器进行几何校正的工具，其中就有专门针对 MODIS 数据的几何校正，如图 2 Georeference MODIS 工具的位置图 2 Georeference MODIS 工具的位置（2）选择发射率文件双击 Georeference MODIS 工具，打开输入 MODIS 数据对话框，选中发射率文件，点击 OK，如图 3 所示 图 3 输入发射率文件（3）Georeference MODIS 参数设置Georeference MODIS Parameters 对话框中，选择投影信息见图 4，其中要注意保持GCP 控制点以及对 MODIS 影像做双眼皮去除，然后点击 OK 图 4 MODIS 几何校正参数设置（4）发射率几何校正参数设置Registration Parameters 对话框中，注意分辨率是 1000（默认）选择保持路径点击OK 即可 图 5 保存几何校正文件（5）几何校正方法建立 GCP 控制点图 6，最后开始进行几何校正和双眼皮去除注意几何校正方法为Triangulation 见图 7。
图 6 建立 GCP 控制点图 7 进行几何校正（6）发射率几何校正结果图 8 发射率几何校正结果2.2 反射率文件的几何校正（1）Wrap from GCPs：Image to Map Registration 工具发射率文件几何校正之后，可以利用导出的 GCP 控制点来校正其他文件在工具箱中找到 Wrap from GCPs：Image to Map Registration 工具，如图 9 所示该工具的位置 图 9 Wrap from GCPs：Image to Map Registration 工具（2）选择 GCP 控制点文件双击 Wrap from GCPs：Image to Map Registration 工具后打开选择 GCP 控制点对话框，选择之前保持的 GCP 控制点文件，如图 10 图 10 选择 GCP 控制点文件（3）Image to Map 投影设置接着打开 Image to Map Registration 对话框，调整参数如图，注意更改分辨率为1000，然后点击 OK，如图 11图 11 参数调整（4）选择待校正影像，调整几何校正方法选择待校正的反射率文件图 12。
打开 Registration Parameters 对话框参数设置如图，注意更改几何校正方法为 Triangulation 和重采样方法 Bilinear，这样才能与之前的发射率校正结果相匹配，如图 13图 12 选择反射率文件图 13 调整几何校正参数（5）反射率几何校正进度和结果然后进行几何校正如图 14，反射率几何校正结果如图 15.图 14 进行几何校正图 15 反射率几何校正结果（6）反射率校正结果与发射率校正结果对比ENVI 5.0 中提供了更多的图像间的对比查看，可以点击 Views 菜单选择创建多个视图，然后进行 link 查看，也可以将两个图像添加进一个视图利用影像查看功能查看，改功能包括四个，一般在工具条上下面我们就以 Portal 来查看几何校正效果将反射率和发射率的第一波段添加进内容列表，然后点击工具条上的 Portal 按钮即可建立一个新图层效果图如图 16，从图上可以看出两个图像几何校正之后完全重合说明该方法可行图 16 Portal 两个几何校正图像2.3 角度数据集的几何校正该教程数据采用的是 MODIS L1B 1KM 数据，其中包括了角度信息，但是角度数据集的行列号与科学数据集的行列号不同，因此如果想用校正发射率的 GCP 控制点来校正角度数据集必须在校正角度数据集之前进行重采样到行列号为 1354 ,2030。
1）打开角度数据集角度数据集的打开不能和科学数据集一样直接打开，打开路径为 File -> Open As ->Generic Formats -> HDF,点击后选择 HDF 文件，接着弹出 HDF 数据集选择对话框选择四个角度数据集，这四个角度数据集分别为：卫星天顶角，卫星方位角，太阳天顶角，太阳方位角，如图 17 图 17 选择角度数据集（2）查看角度数据集选择之后可以再 Data Manager 窗口中看到打开的数据集图 18,如果打开的角度数据集显示的是科学数据集，这是由于 ENVI 5.0 下建立金字塔之后产生.ecp 文件，删除和 HDF同名的.ecp 文件即可可以右击查看元数据中的 auxiliary URL图 18 打开的角度数据集（3）重采样工具从图 18 中可以看到，角度数据集的行列数是 271*406 的，而发射率的行列数是1354*2030 的，因此要用之前的 GCP 文件来校正角度数据集必须重采样重采样工具Resize Data 的位置如图 19 所示图 19 Resize Data 工具（4）选择重采样的角度数据双击 Resize Data 工具，打开输入重采样数据的对话框，选择其中一个角度数据，点击 OK，如图 20 所示。
图 20 选择重采样角度数据（5）重采样参数设置在重采样参数设置窗口中，设置 X Y 的采样比例：X=4.9963 Y=5 点击 Enter，即可采样到 1354*2030，选择重采样方法为 Blinear，将采样结果暂时保存到内存即可，如图 21所示依次将四个角度数据进行重采样之后就可以进行几何校正了图 21 重采样参数设（6）角度数据集几何校正角度数据的几何校正和发射率的几何校正是一样的，根据校正发射率产生的 GCP 控制点，利用 Wrap from GCPs：Image to Map Registration 工具进行校正三、波段合成和裁剪波段合成工具的位置如图 22. 图 22 Layer Satacking 工具3.1 反射率和发射率的合成（1）Layer Satacking 工具在 ENVI5.0 下打开之前校正好的反射率文件和发射率文件，利用 Layer Stacking 工具进行合成，双击 Layer Satacking 工具，打开参数设置窗口如图 23 所示图 23 合成窗口参数设置（2）空间裁剪点击 Import File…选择几何校正后的反射率文件和发射率文件，然后利用北京市的行政范围来裁剪，如图 24 所示。
图 24 选择反射率和发射率文件并选择裁剪文件（3）调整合成文件的顺序点击 Recorder File 调整顺序，必须是反射率在上，发射率在下，如图 25. 图 25 调整文件顺序（4）合成点击 OK 进行合成，合成进度见图 26图 26 合成进度3.2 角度数据的合成角度数据的合成和反射率发射率的合成是一样的，只是要注意角度数据的合成时的顺序是卫星天顶角（Sensor Zenith）、卫星方位角（Sensor Azimuth）、太阳天顶角（Solar Zenith）、太阳方位角（Solar Zenith），顺序如图 27 所示，合成结果如图 28所示图 27 角度数据的顺序图 28 角度合成结果四、合成后处理4.1 云检测云检测工具是扩展工具，该工具实现对反射率和发射率的合成文件进行去云处理，将modis_cloud.sav 文件放在 ENVI 5.0 安装目录下的 Extensions 文件夹下，重启 ENVI 即可看到 modis_cloud 工具1）modis_cloud 工具双击 modis_cloud 工具，选择几何校正结果，就是反射率和发射率合成裁剪后的结果，如图 29 所示。
图 29 选择数据合成结果（2）保存云检测结果选择保存路径，保存云检测结果，如图 30 所示图 30 云检测结果保存4.2 角度数据HDF 中的角度数据时扩大了 100 倍的，所以在进行气溶胶反演之前要将角度合成数据乘以 0.011）band math 工具在工具箱中找到 band math 工具，如图 31 所示 图 31 band math 工具（2）角度数据波段运算双击打开 band math 工具后，写入公式 b1*0.01,如图 32 所示 图 32 写入公式（3）选择角度合成文件图 33 选择角度合成文件（4）保存角度数据处理结果五、气溶胶反演气溶胶反演工具（modis_aerosol_inversion.sav）也是扩展工具，放在 ENVI5.0 安装目录下的 Extensions 文件夹下重启 ENVI 即可5.1 选择云检测结果图 34 选择云检测结果5.2 选择角度合成波段运算结果图 35 选择角度合成结果5.3 选择查找表文件查找表文件是个通用的文本文件，该查找表适合 3-9 月的 MODIS 影像图 36 选择查找表文件5.4 保存气溶胶反演结果 图 37 保存气溶胶反演结果5.5 查看反演结果图 38 北京市气溶胶反演结果。