遥感图像处理实验报告

中国地质大学（武汉）遥感图像处理实习报告姓名：学号：联系方式：日期：中国地质大学（武汉）第一页一、实习要求（一）掌握使用ENVI进行各种图像基本操作；（二）熟练运用ENVI中工具进行图像图像校正、裁剪拼接、融合及图像增强处理；二、实习操作过程与实现结果（一）辐射校正及大气校正1、辐射校正（1）选择File-open，选择Landset8武汉数据中的‘LC81230392016205LGN00_MTL.txt’文件。（2）选择Toolbox-RadiometricCorrection-RadiometricCalibration工具，选择要校正的‘LC81230392016205LGN00_MTL_MultiSpectral’多光谱数据，设置定标参数（存储格式：BIL；单位转换“ScaleFactor”的设置，单击ApplyFLAASHSettings得到相应的参数），得到辐射定标后的结果。中国地质大学（武汉）第二页2、大气校正（1）选择Toolbox-RadiometricCorrection-AtmosphericCorrectionModule-FLAASHAtmosphericCorrection工具；打开工具后设置参数：在FLAASHAtmosphericCorrectionModuleInputParameters面板中如图设置各项参数；点击apply运行大气校正。中国地质大学（武汉）第三页（2）大气校正运行结果（二）图像裁剪与拼接1、15米全色波段图像裁剪拼接（1）选择File-open，选择‘县界.shp’‘LC81230382016205LGN00_MTL’及‘LC81230392016205LGN00_MTL’文件。中国地质大学（武汉）第四页（2）选择Toolbox-RegionsofInterest-SubsetDatefromROIs工具；双击打开后inputfile面板选择38区段15米分辨率文件，inputROIs面板选择‘县界’文件。点击‘OK’，38区段文件裁剪后如图。（3）重复（2）中工具选择步骤；双击打开后在inputfile面板选择39区段15米分辨率文件，在inputROIs面板选择‘县界’文件。点击‘OK’，39区段裁剪后如图。中国地质大学（武汉）第五页（4）选择Toolbox-Mosaicking-SeamlessMosaic工具，双击打开，选择裁剪后的15米分辨率38、39区段文件；设置‘DateIgnoreValue’为0，‘FeatheringDistance’为10，‘ColorCorrection’为‘overlapareaonly’，‘ResembingMethods’设为‘CubicConvolutions’。（5）点击Finish，得到拼接后武汉15米边界图。中国地质大学（武汉）第六页（6）依据（1）（2）（3）步，分别得到裁剪后的30米分辨率图像。注：在裁剪时与上面三步不同的是打开文件为30米分辨率的38、39区段文件。（38区段裁剪30m）（39区段裁剪30m）（7）类似第（4）步，设置拼接参数，选择裁剪好的30米分辨率38、39区段图像文件，得到拼接后武汉30米边界图。中国地质大学（武汉）第七页（三）图像融合1、全色波段与MODIS09A1-H27V05融合（1）选择File-open,打开‘武汉MODIS-UTM-剪切边界’和‘武汉15米边界图’（2）选择Toolbox-ImageSharping-Gram-SchmidtPanSharping,双击打开后选择‘武汉MODIS-UTM-剪切边界’和‘武汉15米边界图’文件。（3）设置‘PanSharpingParameters’面板：Sensor--Landset8oil，Resampling--CubicConvention，output--ENVI。点击‘OK’得到融合后图像。中国地质大学（武汉）第八页（全色与MODIS融合）（全色与MODIS融合细节图）2、全色波段与多光谱波段图像融合（1）选择File-open,打开‘武汉15米边界图’和‘武汉边界30米’文件。（2）选择Toolbox-ImageSharping-Gram-SchmidtPanSharping,双击打开后选择‘武汉边界图30米’和‘武汉15米边界图’文件。中国地质大学（武汉）第九页（3）设置‘PanSharpingParameters’面板：Sensor--Landset8oil，Resampling--CubicConvention，output--ENVI。点击‘OK’得到融合后图像。（全色与多光谱融合图像）（全色与多光谱融合图像细节图）中国地质大学（武汉）第十页3、两次融合对比图像（右为全色与MODIS融合，左为全色与多光谱融合）（对比细节图）对比结果：左图共有7个波段，右图只有一个波段。左图的亮度低，对比度高，对地物细节还原较好，对于肉眼分辨地物更实用，可以明显分别出房屋，道路，植被，水域等；而右图的分辨率则较低，但也可以分辨出河流，道路等地物。中国地质大学（武汉）第十一页（四）图像增强1、辐射增强处理（1）交互式直方图拉伸①选择File-open,打开30米武汉边界图文件。选择菜单栏中Display-CustomStretch，显示RGB三个波段直方图。②选择R波段，移动边界，得到前后两幅不同图像（前）中国地质大学（武汉）第十二页（后）③选择G波段，同样操作，得到两幅图像。（前）中国地质大学（武汉）第十三页（后）④选择B波段，同样操作得到两幅图像（前）中国地质大学（武汉）第十四页（后）⑤Equalization直方图均衡化（以R波段为例）（前）中国地质大学（武汉）第十五页（后）⑥Gaussian高斯拉伸（以R波段为例）（前）中国地质大学（武汉）第十六页（后）⑦SquareRoot平方根（以R波段为例）（前）中国地质大学（武汉）第十七页（后）⑧Logarithmic对数拉伸（以R波段为例）（前）中国地质大学（武汉）第十八页（后）（2）去条带①选择Toolbox-Raster-Destripe,双击工具在Destripinginputfile面板中选择未处理的武汉边界30米文件。②设置Detectors：条带周期——LandsatTM值为16。点击‘OK’得到处理后图像。中国地质大学（武汉）第十九页2、假彩色合成由上一个去条带得到文件‘destrip’文件，右键点击文件，选择‘ChangeRGBBands…’，分别是红绿蓝波段对应543波段，点击‘ok’加载得到假彩色图像。中国地质大学（武汉）第二十页三、评价及结论通过这次上机实习操作，我掌握了基本的ENVI操作，熟悉了基本的图像处理方法，但在上机时存在的一些问题需要总结经验，不断学习。一是对于ENVI工具箱Toolbox中基本工具和遥感课程中一些基本专业词汇的英文名称不熟悉，这需要发挥主观能动性，在操作学习中不断积累，也希望老师在课堂上讲解时能够涉及到相关概念的英文名称；二是在进行图像融合时，由于时间原因没有进行同一分辨率图像不同计算方式的融合对比，但我认为这种不同计算方式融合对比更有利于我直观地掌握各种计算方式的优、缺点；三是在保存文件过程中应该更加注重层次和条理化，这样有利于后期的再处理，也方便查询。