利用MODIS的L1B数据反演气溶胶的流程如下

1、该方案中的数据和算法介绍：1、MODIS数据是采用的 MODIS L1B 1KM数据。应严格按照说明进行操作，例如数据是 1km的，数据的合成是反射率在上，发射率在下等，下面将详细介绍各个步骤。2、这是在 ENVI 下做的北京市气溶胶反演，具体包括MODIS影像（ HDF）的辐射校正、几何校正、云检测、气溶胶反演。3、气溶胶反演算法采用经典的暗像元法（DDV）也叫浓密植被法，因此对于冬季反演的气溶胶效果不好。4、气溶胶反演的查找表是利用 IDL 调用 6S 辐射模型得到的，采用的是一般的参数， 因此 3-9 月期间都可以用这一个查找表进行气溶胶反演，也可以自己制定查找表。5、七纬查找表， 从左

2、向右， 依次为太阳天顶角， 卫星天顶角， 相对方位角， P0?、T、?S （辐射传输方程参数），最后一列为气溶胶光学厚度（AOD）6. 感兴趣的可以提供用到的 modis云检测工具和气溶胶反演工具以及详细的 pdf 文档。MODIS数据下载地址 ： ?一、MODIS影像的辐射校正在 ENVI 中打开 MODIS影像的 HDF文件就已经做了辐射校正，打开 HDF文件的方法是 File-Open As-EOS-MODIS, 打开后在数据列表中可以看到三个文件，第一个是发射率 Emissive （band20-band36 ），第二个是辐射率 Radiance （ band1-band26 ）, 第

3、三个是反射 率 Reflectance （ band1-band26 ），如图 1 所示。?图1 打开 HDF文件二、几何校正发射率文件的几何校正（1）Georeference MODIS 工具下有对特定传感器进行几何校正的工具，其中就有专门针对MODIS数据的几何校正，如图 2 Georeference MODIS 工具的位置。图 2 Georeference MODIS 工具的位置（2）选择发射率文件双击 Georeference MODIS工具，打开输入 MODIS数据对话框，选中发射率文件，点击 OK，如图 3 所示。图 3 输入发射率文件（3）Georeference MODIS 参

4、数设置Georeference MODISP arameters 对话框中， 选择投影信息见图 4，其中要注意保持 GCP 控制点以及对 MODIS影像做双眼皮去除，然后点击 OK。图 4 MODIS 几何校正参数设置（4）发射率几何校正参数设置Registration Parameters 对话框中，注意分辨率是 1000 （默认）选择保持路径点击 OK即可。图5 保存几何校正文件（5）几何校正方法建立 GCP控制点图 6，最后开始进行几何校正和双眼皮去除注意几何校正方法为Triangulation 见图 7。图 6 建立 GCP控制点图7 进行几何校正6）发射率几何校正结果图8 发射率几何

5、校正结果反射率文件的几何校正（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控制点

6、文件，如图 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）反射率

7、几何校正进度和结果然后进行几何校正如图 14，反射率几何校正结果如图 15.图 14 进行几何校正图 15 反射率几何校正结果（6）反射率校正结果与发射率校正结果对比ENVI 中提供了更多的图像间的对比查看，可以点击 Views 菜单选择创建多个视图， 然 后进行 link 查看，也可以将两个图像添加进一个视图利用影像查看功能查看，改功能包括 四个，一般在工具条上面我们就以 Portal 来查看几何校正效果。将反射率和发射率的第一波段添加进内容列表， 然后点击工具条上的 Portal 按钮即可 建立一个新图层。效果图如图 16，从图上可以看出两个图像几何校正之后完全重合。说明 该方法可行。图

8、16 Portal 两个几何校正图像角度数据集的几何校正该教程数据采用的是 MODIS L1B 1KM数据，其中包括了角度信息，但是角度数据集的 行列号与科学数据集的行列号不同， 因此如果想用校正发射率的 GCP控制点来校正角度数据 集必须在校正角度数据集之前进行重采样到行列号为1354 ,2030 。（1）打开角度数据集角度数据集的打开不能和科学数据集一样直接打开，打开路径为File?- Open As-Generic Formats - HDF, 点击后选择 HDF文件，接着弹出 HDF数据集选择对话框。选择 四个角度数据集，这四个角度数据集分别为：卫星天顶角，卫星方位角，太阳天顶角，太阳

9、 方位角，如图 17 。图 17 选择角度数据集（2）查看角度数据集选择之后可以再 Data Manager 窗口中看到打开的数据集图 18, 如果打开的角度数据集 显示的是科学数据集，这是由于 ENVI 下建立金字塔之后产生 .ecp 文件，删除和 HDF同名 的 .ecp 文件即可。可以右击查看元数据中的auxiliary URL 。图 18 打开的角度数据集（3）重采样工具从图 18 中可以看到，角度数据集的行列数是 271*406 的，而发射率的行列数是 1354*2030 的，因此要用之前的 GCP文件来校正角度数据集必须重采样。 重采样工具 Resize Data 的位置如图 19

10、 所示。（4）选择重采样的角度数据图 19 Resize Data 工具双击 Resize Data 工具，打开输入重采样数据的对话框，选择其中一个角度数据，点 击 OK，如图 20 所示。图 20 选择重采样角度数据（5）重采样参数设置在重采样参数设置窗口中，设置 X Y 的采样比例： X= Y=5 点击 Enter ，即可采样到 1354*2030 ，选择重采样方法为 Blinear ，将采样结果暂时保存到内存即可，如图 21 所示。 依次将四个角度数据进行重采样之后就可以进行几何校正了。?图 21 重采样参数设（6）角度数据集几何校正角度数据的几何校正和发射率的几何校正是一样的，根据校正

11、发射率产生的GCP控制点，利用 Wrap from GCPs ： Image to Map Registration工具进行校正。三、波段合成和裁剪波段合成工具的位置如图 22.?图 22 Layer Satacking 工具反射率和发射率的合成（1）Layer Satacking 工具 在下打开之前校正好的反射率文件和发射率文件， 利用 Layer Stacking 工具进行合成， 双击 Layer Satacking 工具，打开参数设置窗口。如图 23 所示。图 23 合成窗口参数设置（2）空间裁剪点击 Import File 选择几何校正后的反射率文件和发射率文件， 然后利用北京市的行

12、政范围来裁剪，如图 24 所示。（3）调整合成文件的顺序点击 Recorder File 调整顺序，必须是反射率在上，发射率在下，如图图 24 选择反射率和发射率文件并选择裁剪文件25.图 25 调整文件顺序(4) 合成点击 OK进行合成，合成进度见图 26。图 26 合成进度角度数据的合成 角度数据的合成和反射率发射率的合成是一样的，只是要注意角度数据的合成时的顺 序是卫星天顶角( Sensor Zenith )、卫星方位角( Sensor Azimuth )、太阳天顶角( Solar Zenith )、太阳方位角( Solar Zenith )，顺序如图 27 所示，合成结果如图 28 所

13、示。图 27 角度数据的顺序图 28 角度合成结果四、合成后处理云检测 云检测工具是扩展工具，该工具实现对反射率和发射率的合成文件进行去云处理，将文件放在 ENVI 安装目录下的 Extensions 文件夹下，重启 ENVI 即可看到 modis_cloud 工具。(1) modis_cloud 工具双击 modis_cloud 工具， 选择几何校正结果，就是反射率和发射率合成裁剪后的结果，如图 29 所示。图 29 选择数据合成结果(2) 保存云检测结果选择保存路径，保存云检测结果，如图30 所示。?图 30 云检测结果保存角度数据HDF中的角度数据时扩大了 100 倍的，所以在进行气溶胶反演之前要将角度合成数据乘 以。(1) band math 工具在工具箱中找到 band math 工具，如图 31 所示。图 31 band math 工具 （2）角度数据波段运算 双击打开 band math 工具后，写入公式 b1*, 如图 32 所示。?图 32 写入公式（3）选择角度合成文件图 33 选择角度合成文件（4）保存角