利用MODIS的L1B数据反演气溶胶的流程如下.docx
- 文档编号:12269454
- 上传时间:2023-04-17
- 格式:DOCX
- 页数:28
- 大小:1.28MB
利用MODIS的L1B数据反演气溶胶的流程如下.docx
《利用MODIS的L1B数据反演气溶胶的流程如下.docx》由会员分享,可在线阅读,更多相关《利用MODIS的L1B数据反演气溶胶的流程如下.docx(28页珍藏版)》请在冰豆网上搜索。
利用MODIS的L1B数据反演气溶胶的流程如下
利用MODIS的L1B数据反演气溶胶的流程如下:
该方案中的数据和算法介绍:
1、MODIS数据是采用的MODISL1B1KM数据。
应严格按照说明进行操作,例如数据是1km的,数据的合成是反射率在上,发射率在下等,下面将详细介绍各个步骤。
2、这是在ENVI5.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影像的辐射校正
在ENVI5.0中打开MODIS影像的HDF文件就已经做了辐射校正,打开HDF文件的方法是File->OpenAs->EOS->MODIS,打开后在数据列表中可以看到三个文件,第一个是发射率Emissive(band20-band36),第二个是辐射率Radiance(band1-band26),第三个是反射率Reflectance(band1-band26),如图1所示。
图1ENVI5.0打开HDF文件
二、几何校正
2.1发射率文件的几何校正
(1)GeoreferenceMODIS工具
ENVI5.0下有对特定传感器进行几何校正的工具,其中就有专门针对MODIS数据的几何校正,如图2GeoreferenceMODIS工具的位置。
图2GeoreferenceMODIS工具的位置
(2)选择发射率文件
双击GeoreferenceMODIS工具,打开输入MODIS数据对话框,选中发射率文件,点击OK,如图3所示。
图3输入发射率文件
(3)GeoreferenceMODIS参数设置
GeoreferenceMODISParameters对话框中,选择投影信息见图4,其中要注意保持GCP控制点以及对MODIS影像做双眼皮去除,然后点击OK。
图4MODIS几何校正参数设置
(4)发射率几何校正参数设置
RegistrationParameters对话框中,注意分辨率是1000(默认)选择保持路径点击OK即可。
图5保存几何校正文件
(5)几何校正方法
建立GCP控制点图6,最后开始进行几何校正和双眼皮去除注意几何校正方法为Triangulation见图7。
图6建立GCP控制点
图7进行几何校正
(6)发射率几何校正结果
图8发射率几何校正结果
2.2反射率文件的几何校正
(1)WrapfromGCPs:
ImagetoMapRegistration工具
发射率文件几何校正之后,可以利用导出的GCP控制点来校正其他文件。
在工具箱中找到WrapfromGCPs:
ImagetoMapRegistration工具,如图9所示该工具的位置。
图9WrapfromGCPs:
ImagetoMapRegistration工具
(2)选择GCP控制点文件
双击WrapfromGCPs:
ImagetoMapRegistration工具后打开选择GCP控制点对话框,选择之前保持的GCP控制点文件,如图10。
图10选择GCP控制点文件
(3)ImagetoMap投影设置
接着打开ImagetoMapRegistration对话框,调整参数如图,注意更改分辨率为1000,然后点击OK,如图11。
图11参数调整
(4)选择待校正影像,调整几何校正方法
选择待校正的反射率文件图12。
打开RegistrationParameters对话框参数设置如图,注意更改几何校正方法为Triangulation和重采样方法Bilinear,这样才能与之前的发射率校正结果相匹配,如图13。
图12选择反射率文件
图13调整几何校正参数
(5)反射率几何校正进度和结果
然后进行几何校正如图14,反射率几何校正结果如图15.
图14进行几何校正
图15反射率几何校正结果
(6)反射率校正结果与发射率校正结果对比
ENVI5.0中提供了更多的图像间的对比查看,可以点击Views菜单选择创建多个视图,然后进行link查看,也可以将两个图像添加进一个视图利用影像查看功能查看,改功能包括四个,一般在工具条上
。
下面我们就以Portal来查看几何校正效果。
将反射率和发射率的第一波段添加进内容列表,然后点击工具条上的Portal按钮即可建立一个新图层。
效果图如图16,从图上可以看出两个图像几何校正之后完全重合。
说明该方法可行。
图16Portal两个几何校正图像
2.3角度数据集的几何校正
该教程数据采用的是MODISL1B1KM数据,其中包括了角度信息,但是角度数据集的行列号与科学数据集的行列号不同,因此如果想用校正发射率的GCP控制点来校正角度数据集必须在校正角度数据集之前进行重采样到行列号为1354,2030。
(1)打开角度数据集
角度数据集的打开不能和科学数据集一样直接打开,打开路径为File ->OpenAs->GenericFormats->HDF,点击后选择HDF文件,接着弹出HDF数据集选择对话框。
选择四个角度数据集,这四个角度数据集分别为:
卫星天顶角,卫星方位角,太阳天顶角,太阳方位角,如图17。
图17选择角度数据集
(2)查看角度数据集
选择之后可以再DataManager窗口中看到打开的数据集图18,如果打开的角度数据集显示的是科学数据集,这是由于ENVI5.0下建立金字塔之后产生.ecp文件,删除和HDF同名的.ecp文件即可。
可以右击查看元数据中的auxiliaryURL。
图18打开的角度数据集
(3)重采样工具
从图18中可以看到,角度数据集的行列数是271*406的,而发射率的行列数是1354*2030的,因此要用之前的GCP文件来校正角度数据集必须重采样。
重采样工具ResizeData的位置如图19所示。
图19ResizeData工具
(4)选择重采样的角度数据
双击ResizeData工具,打开输入重采样数据的对话框,选择其中一个角度数据,点击OK,如图20所示。
图20选择重采样角度数据
(5)重采样参数设置
在重采样参数设置窗口中,设置XY的采样比例:
X=4.9963Y=5点击Enter,即可采样到1354*2030,选择重采样方法为Blinear,将采样结果暂时保存到内存即可,如图21所示。
依次将四个角度数据进行重采样之后就可以进行几何校正了。
图21重采样参数设
(6)角度数据集几何校正
角度数据的几何校正和发射率的几何校正是一样的,根据校正发射率产生的GCP控制点,利用WrapfromGCPs:
ImagetoMapRegistration工具进行校正。
三、波段合成和裁剪
波段合成工具的位置如图22.
图22LayerSatacking工具
3.1反射率和发射率的合成
(1)LayerSatacking工具
在ENVI5.0下打开之前校正好的反射率文件和发射率文件,利用LayerStacking工具进行合成,双击LayerSatacking工具,打开参数设置窗口。
如图23所示。
图23合成窗口参数设置
(2)空间裁剪
点击ImportFile…选择几何校正后的反射率文件和发射率文件,然后利用北京市的行政范围来裁剪,如图24所示。
图24选择反射率和发射率文件并选择裁剪文件
(3)调整合成文件的顺序
点击RecorderFile调整顺序,必须是反射率在上,发射率在下,如图25.
图25调整文件顺序
(4)合成
点击OK进行合成,合成进度见图26。
图26合成进度
3.2角度数据的合成
角度数据的合成和反射率发射率的合成是一样的,只是要注意角度数据的合成时的顺序是卫星天顶角(SensorZenith)、卫星方位角(SensorAzimuth)、太阳天顶角(SolarZenith)、太阳方位角(SolarZenith),顺序如图27所示,合成结果如图28所示。
图27角度数据的顺序
图28角度合成结果
四、合成后处理
4.1云检测
云检测工具是扩展工具,该工具实现对反射率和发射率的合成文件进行去云处理,将modis_cloud.sav文件放在ENVI5.0安装目录下的Extensions文件夹下,重启ENVI即可看到modis_cloud工具。
(1)modis_cloud工具
双击modis_cloud工具,选择几何校正结果,就是反射率和发射率合成裁剪后的结果,如图29所示。
图29选择数据合成结果
(2)保存云检测结果
选择保存路径,保存云检测结果,如图30所示。
图30云检测结果保存
4.2角度数据
HDF中的角度数据时扩大了100倍的,所以在进行气溶胶反演之前要将角度合成数据乘以0.01。
(1)bandmath工具
在工具箱中找到bandmath工具,如图31所示。
图31bandmath工具
(2)角度数据波段运算
双击打开bandmath工具后,写入公式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北京市气溶胶反演结果
(此文档部分内容来源于网络,如有侵权请告知删除,文档可自行编辑修改内容,供参考,感谢您的配合和支持)
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 利用 MODIS L1B 数据 反演 气溶胶 流程 如下