1、中国地质大学（武汉）遥感图像处理实习报告 中国地质大学（武汉）姓名：学号：联系方式：日期：中国地质大学（武汉）第 1 页1、实习要求（一）掌握使用 ENVI 进行各种图像基本操作；（二）熟练运用 ENVI 中工具进行图像图像校正、裁剪拼接、融合及图像增强处理；2、实习操作过程与实现结果（1）辐射校正及大气校正1、辐射校正（1）选择 File-open，选择 Landset8 武汉数据中的LC81230392016205LGN00_MTL.txt文件。（2）选择 Toolbox-Radiometric Correction-Radiometric Calibration 工具，选择要校正的LC8
2、1230392016205LGN00_MTL_MultiSpectral多光谱数据，设置中国地质大学（武汉）第 2 页定标参数（存储格式：BIL；单位转换“Scale Factor”的设置，单击 Apply FLAASH Settings 得到相应的参数） ，得到辐射定标后的结果。2、大气校正（1）选择 Toolbox-Radiometric Correction-Atmospheric Correction Module-FLAASH Atmospheric Correction 工具；打开工具后设置参数：在 FLAASH Atmospheric Correction Module Inpu
3、t Parameters 面板中如图设置各项参数；点击 apply 运行大气校正。中国地质大学（武汉）第 3 页（2）大气校正运行结果（2）图像裁剪与拼接1、15 米全色波段图像裁剪拼接（1）选择 File-open，选择县界.shp中国地质大学（武汉）第 4 页LC81230382016205LGN00_MTL及LC81230392016205LGN00_MTL文件。（2）选择 Toolbox-Regions of Interest-Subset Date from ROIs 工具；双击打开后 input file 面板选择 38 区段 15 米分辨率文件，input ROIs 面板选择县界
4、文件。点击OK ，38 区段文件裁剪后如图。中国地质大学（武汉）第 5 页（3）重复（2）中工具选择步骤；双击打开后在 input file 面板选择 39 区段 15 米分辨率文件，在 input ROIs 面板选择县界文件。点击OK ，39 区段裁剪后如图。（4）选择 Toolbox-Mosaicking-Seamless Mosaic 工具，双击打开，选择裁剪后的 15 米分辨率 38、39 区段文件；设置Date Ignore Value为 0， Feathering Distance为 10， Color Correction为overlap area only ， Resembin
5、g Methods设为Cubic Convolutions 。中国地质大学（武汉）第 6 页（5）点击 Finish，得到拼接后武汉 15 米边界图。（6）依据（1） （2） （3）步，分别得到裁剪后的 30 米分辨率图像。注：在裁剪时与上面三步不同的是打开文件为 30 米分辨率的中国地质大学（武汉）第 7 页38、39 区段文件。（38 区段裁剪 30m）（39 区段裁剪 30m）（7）类似第（4）步，设置拼接参数，选择裁剪好的 30 米分辨率38、39 区段图像文件，得到拼接后武汉 30 米边界图。中国地质大学（武汉）第 8 页（3）图像融合1、全色波段与 MODIS09A1-H27V05
6、 融合（1）选择 File-open,打开武汉 MODIS-UTM-剪切边界和武汉15 米边界图（2）选择 Toolbox-Image Sharping-Gram-Schmidt Pan Sharping,双击打开后选择武汉 MODIS-UTM-剪切边界和武汉15 米边界图文件。（3）设置Pan Sharping Parameters面板：Sensor-Landset 8 oil，Resampling-Cubic Convention，out put-ENVI。点击OK得到融合后图像。中国地质大学（武汉）第 9 页（全色与 MODIS 融合）（全色与 MODIS 融合细节图）2、全色波段与多光
7、谱波段图像融合（1）选择 File-open,打开武汉 15 米边界图和武汉边界 30米文件。中国地质大学（武汉）第 10 页（2）选择 Toolbox-Image Sharping-Gram-Schmidt Pan Sharping,双击打开后选择武汉边界图 30 米和武汉 15 米边界图文件。（3）设置Pan Sharping Parameters面板：Sensor-Landset 8 oil，Resampling-Cubic Convention，out put-ENVI。点击OK得到融合后图像。（全色与多光谱融合图像）中国地质大学（武汉）第 11 页（全色与多光谱融合图像细节图）3、两
8、次融合对比图像（右为全色与 MODIS 融合，左为全色与多光谱融合）中国地质大学（武汉）第 12 页（对比细节图）对比结果：左图共有 7 个波段，右图只有一个波段。左图的亮度低，对比度高，对地物细节还原较好，对于肉眼分辨地物更实用，可以明显分别出房屋，道路，植被，水域等；而右图的分辨率则较低，但也可以分辨出河流，道路等地物。（4）图像增强1、辐射增强处理（1）交互式直方图拉伸选择 File-open,打开 30 米武汉边界图文件。选择菜单栏中Display-Custom Stretch，显示 R G B 三个波段直方图。中国地质大学（武汉）第 13 页选择 R 波段，移动边界，得到前后两幅不同
9、图像（前）中国地质大学（武汉）第 14 页（后）选择 G 波段，同样操作，得到两幅图像。（前）中国地质大学（武汉）第 15 页（后）选择 B 波段，同样操作得到两幅图像中国地质大学（武汉）第 16 页（前）中国地质大学（武汉）第 17 页（后）Equalization 直方图均衡化（以 R 波段为例）（前）中国地质大学（武汉）第 18 页（后）Gaussian 高斯拉伸（以 R 波段为例）（前）中国地质大学（武汉）第 19 页（后）Square Root 平方根（以 R 波段为例）中国地质大学（武汉）第 20 页（前）（后）中国地质大学（武汉）第 21 页Logarithmic 对数拉伸（以
10、R 波段为例）（前）（后）中国地质大学（武汉）第 22 页（2）去条带选择 Toolbox-Raster-Destripe,双击工具在 Destriping input file 面板中选择未处理的武汉边界 30 米文件。设置 Detectors：条带周期Landsat TM 值为 16。点击OK得到处理后图像。中国地质大学（武汉）第 23 页2、假彩色合成由上一个去条带得到文件destrip文件，右键点击文件，选择Change RGB Bands ，分别是红绿蓝波段对应 543 波段，点击ok加载得到假彩色图像。中国地质大学（武汉）第 24 页3、评价及结论通过这次上机实习操作，我掌握了基本的 ENVI 操作，熟悉了基本的图像处理方法，但在上机时存在的一些问题需要总结经验，不断学习。一是对于 ENVI 工具箱 Toolbox 中基本工具和遥感课程中一些基本专业词汇的英文名称不熟悉，这需要发挥主观能动性，在操作学习中不断积累，也希望老师在课堂上讲解时能够涉及到相关概念的英文名称；二是在进行图像融合时，由于时间原因没有进行同一分辨率图像不同计算方式的融合对比，但我认为这种不同计算方式融合对比更有利于我直观地掌握各种计算方式的优、缺点；三是在保存文件过程中应该更加注重层次和条理化，这样有利于后期的再处理，也方便查询。
