《遥感实验ENVI.doc》由会员分享，可在线阅读，更多相关《遥感实验ENVI.doc（17页珍藏版）》请在三一文库上搜索。
1、遥感实验 ENVI在envi5.0中打开modis影像文件HDF-file-open as-eosmodis早打开的图像中 1、MODIS影像的辐射校正 在ENVI 5.0中打开MODIS影像的HDF文件就已经做了辐射校正，先将文件里的exp文件删掉，仅仅留下hdf，打开HDF文件的方法是File-Open As-MODIS，打开后在数据列表中可以看到三个文件，第一个是发射率Emissive（band20-band36），第二个是辐射率Radiance（band1-band26），第三个是反射率Reflectance（band1-band26）， 图1 ENVI5.0打开HDF文件 2、几何校
2、正 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对话框中，选择投
3、影信息，见图4，其中要注意保持GCP控制点以及对MODIS影像做双眼皮去除，choose，然后点击OK保存。 图4 MODIS几何校正参数设置 （4）发射率几何校正参数设置 Registration Parameters对话框中，注意分辨率是1000（默认）点击choose选择保持路径点击OK即可。图5 保存几何校正文件 （5）几何校正方法 建立GCP控制点，如图6，最后开始进行几何校正和双眼皮去除，注意几何校正方法为Triangulation，见图7. 图6 建立GCP控制点 图7 进行几何校正 （6）发射率几何校正结果 图8 发射率几何校正结果 2.2反射率文件的几何校正 （1）Wrap
4、from GCP：Image to Map Registration工具 发射率文件几何校正之后，可以利用导出的GCP控制点来校正其他文件。在工具箱中找到Wrap from GCP：Image to Map Registration工具，如图9所示该工具的位置。 图9 Wrap from GCPs：Image to Map Registration工具 （2）选择GCP控制点文件 双击Wrap from GCP：Image to Map Registration工具后打开选择GCP控制点对话框，选择之前保持的GCP控制点文件，如图10. 图10 选择GCP控制点文件 （3）Image to M
5、ap投影设置 接着打开Image to Map Registration对话框，调整参数如图，注意更改分辨率为1000，然后点击OK，如图11。 图11 参数调整 （4）选择待校正影像，调整几何校正方法 在上一步ok后，选择待校正的反射率文件，如图12,。打开Registration Parameters对话框参数设置如图，注意更改几何校正方法为Triangulation和重采样方法为Bilinear，这样才能与之前的发射率校正结果相匹配，如图13. 图12 选择发射率文件 图13 调整几何校正参数 （5）反射率几何校正进度和结果 然后进行几何校正如图14，反射率几何校正结果如图15。 图14
6、 进行几何校正 图15 反射率几何校正结果 （6）反射率校正结果与发射率校正结果对比 ENVI5.0中提供了更多的图像间的对比查看，可以点击Views菜单选择创建多个视图，然后进行link查看，也可以将两个图像添加进一个视图利用影像查看功能查看，该功能括四个，一般在工具条上。下面我们就以Portal来查看几何校正效果。 将反射率和发射率的第一波段添加进内容列表，然后点击工具条上的Portal按钮即可建立一个新图层。效果图如图16，从图上可以看出两个图像几何校正之后完全重合。说明该方法可行。 图16 Portal两个几何校正图像 2.3角度数据集的几何校正 本实验采用的是MODIS L1B 1K
7、M数据，其中包括了角度信息，但是角度数据集的行列号与科学数据集的行列号不同，因此如果想用校正发射率的GCP控制点来校正角度数据集必须在校正角度数据 集之前进行重采样到行列号为1354,2030. （1）打开角度数据集 角度数据集的打开不能和科学数据集一样直接打开，打开路径为File-Open As-Generic Formats-HDF，点击后选择HDF文件，接着弹出HDF数据集选择对话框。选择四个角度数据集，这四个角度数据集分别为：卫星天顶角，卫星方位角，太阳天顶角，太阳方位角，如图17。 图17 选择角度数据集（2）查看角度数据集 选择之后可以在Data Manager窗口中看到打开的数据
8、集图18，如果打开的角度数据集显示的是科学数据集，这是由于ENVI5.0下建立金字塔之后产生.ecp文件，删除和HDF同名的.ecp文件即可。可以右击查看元数据中的auxiliary URL。 将envi5.0关闭，在文件夹里删除enp文件，然后打开数据，操作如上面的步骤，再重新来过然后点击菜单第二个Data Manager得到如下的图像。 图18 打开的角度数据集 （3）重采样工具 从图18中可以看到，角度数据集的行列数是271*406的，而发射率的行列数是1354*2030的，因此要用之前的GCP文件来校正角度数据集必须重采样。重采样工具Resize Data的位置如图19所示。 图19
9、Resize Data工具 （4）选择重采样的角度数据 双击Resize Data工具，打开输入重采样数据的对话框，选择其中一个角度数据，点击OK，如图20所示。 图20 选择重采样角度数据 （5）重采样参数设置 在重采样参数设置窗口中，设置XY的采样比例：X=4.9963 Y=5点击Enter，即可采样到1354*2030，选择重采样方法为Blinear，将采样结果暂时保存到内存即可，如图21所示。依次将四个角度数据进行重采样之后就可以进行几何校正了。 图21 重采样参数设置 （6）角度数据集几何校正 角度数据的几何校正和发射率的几何校正是一样的（应该是反射率的步骤一样的，要用到发射率生成的
10、控制点来完成），根据校正发射率产生的GCP控制点，利用Wrap From GCPs：Image to Map Registration工具进行校正。双击Image to Map Registration，添加发射率生成的控制点GCP文件，打开，在弹出的对话框image to map registration中修改相关的坐标为1000，ok。在弹出的input warp image中选中重采样好的数据，其他如上面的反射率步骤和相关数据的修改，并将以此经四个几何校正。3、波段合成和裁剪 波段合成工具的位置如图22。 图22 Layer Satacking工具 3.1反射率和发射率的合成 （1）La
11、yer Satacking工具 在ENVI5.0下打开之前矫正好的反射率文件和发射率文件，利用Layer Stacking工具进行合成，双击Layer Satacking工具，打开参数设置窗口。如图23所示。 图23 合成窗口参数设置 （2）空间裁剪 点击Import File选择几何校正后的反射率文件和发射率文件，然后利用行政区范围来裁剪，如图24所示。在点击ROI/EVF,选择气溶胶数据-辅助数据北京市范围（矢量）-beijing（一定要选中后再点击ok），打开，ok。 图24 选择反射率和发射率文件并选择裁剪文件 （3）调整合成文件的顺序 点击Recorder File调整顺序，必须是反
12、射率在上，发射率在下，如图25。 图25 调整文件顺序 （4）合成 点击OK进行合成，合成进度见图26。 图26 合成进度 3.2角度数据的合成 角度数据的合成和反射率发射率的合成是一样的，只是要注意角度数据的合成时的顺序是卫星天顶角、卫星方位角、太阳天顶角、太阳方位角，顺序如图27所示，合成结果如图28所示。 点击layer stackingimport file在打开的窗口中选择几何校正好的文件点击，出现变化后，点击spatial subset-再出现的对话框中点击ROI/EVF-在弹出的enterROI/EVF filename中选择辅助数据-背景矢量数据中的-beijing打开，在打开
13、的对话框中仍要选中beijing，选中变成蓝色的时候就可以点击ok。然后继续操作。 图27 角度数据的顺序 图28 角度合成结果 4合成后处理 4.1云检测 云检测工具是扩展工具，该工具实现对反射率和发射率的合成文件进行去云处理，将modis_cloud.sav文件放在ENVI5.0安装目录下的Extensions文件夹下，重启ENVI即可看到modis_cloud工具。 （1）modis_cloud工具 双击modis_cloud工具，选择几何校正结果，就是反射率和发射率合成裁剪后的结果，如图29所示。 图29 选择数据合成结果 （2）保存云检测结果 选择保存路径，保存云检测结果，如图30所
14、示。 图30 云检测结果保存 4.2角度数据 HDF中的角度数据是扩大了100倍的，所以在进行气溶胶反演之前要将角度合成数据乘以0.01。 （1）band math工具 在工具箱中找到band math工具，如图31所示。 图31 band math工具 （2）角度数据波段运算 双击打开band math工具后，写入公式b1*0.01，如图32所示。 图32 写入公式 （3）选择角度合成文件 图33 选择角度合成文件 （4）保存角度数据处理结果 5气溶胶反演 气溶胶反演工具（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 最终反演结果