中国矿业大学遥感实习报告参考模板

《遥感原理与应用》实习报告学院：班级：学号：姓名：中国矿业大学环境与测绘学院遥感实习2010年12月26日遥感原理与应用-1-目录1、实验一高光谱数据分析………………………………………..22、实验二影像镶嵌……...……………………………………….153、实验三影像配准………………………………………………264、实验四非监督分类……………………………………………385、实验五监督分类………………………………………………416、实习总结………………………………………………………..49遥感原理与应用-2-实验一高光谱数据分析一、实验目的本专题旨在向用户介绍波谱库的概念，并描述如何从感兴趣区中提取波谱信息，然后还将进行彩色合成，并使用二维散点图进行简单的分类。让学生学会如何使用ENVI先进的高光谱工具对多光谱数据进行分析。更好地理解高光谱处理的概念及其工具。本专题将从特定矿物质的感兴趣区中提取其波谱曲线，并与波谱库中的波谱曲线进行比较，找出显示波谱信息的最佳RGB彩色组合。使用二维散点图定位独特的像元，探究其数据的分布特点，然后进行简单的分类。二、实验数据介绍文件描述zhumuqian95_at.intCuprite地区ATREM校正后的反射率数据，50个波段（整型）zhumuqian95_at.hdrENVI相应的头文件Jpl1.sliENVI格式的JPL波谱库Jpl1.hdrENVI相应的头文件zhumuqian_min.sliENVI格式的USGS波谱库zhumuqian_min.hdrENVI相应的头文件zhumuqian95_av.roi保存的感兴趣区文件数据：1995年的航空可见光/红外成像光谱仪（AirborneVisible/InfraredImagingSpectrometer，简称为AVIRIS）所采集的表观反射率数据，该数据是美国内华达州（Nevada）Cuprite地区的表观反射率数据，它使用ATREM大气纠正建模软件进行了校正。这个数据子集共包含50个波段，波谱分辨率近似为10nm宽，其波长范围为1.99~2.48μm。三、实验内容本部分将介绍以下内容：波谱库操作、浏览和提取影像反射率波谱、ENVI中感兴趣区（ROI）的定义及进行彩色合成影像的选取，其目的是为了鉴别波谱类型。1加载AVIRIS影像数据并显示灰阶影像在ENVI主菜单中，选择File→OpenImageFile，然后选择进入envi47/data目录。选择zhumuqian95_at.int文件作为输入文件名，点击Open弹出可用波段列表，它将列出50个波段的名字。在可用波段列表对话框中，选择Band193（2.2008um），点击GrayScale单选按钮，然后点击LoadBand。将灰度影像加载到显示窗口中。遥感原理与应用-3-从主影像窗口菜单中选择ToolsProfilesZProfile(Spectrum)，提取表观反射率波谱曲线。2浏览影像波谱并同波谱库进行比较在影像上移动缩放指示矩形框，同时查看#1SpectralProfile窗口中的波谱曲线，浏览整个影像的表观反射率波谱曲线。在主影像窗口中，使用鼠标左键点击并拖动缩放指示矩形框或者直接点击鼠标左键，将缩放指示矩形框移动到以所选像素点为中心的区域中。将从影像中获取的表观反射率波谱曲线同所选波谱库中的波谱曲线进行比较。ENVI提供了几个不同的波谱库，根据本专题的目的，我们将会使用JPL波谱库（Groves等，1992）以及USGS波谱库（Clarke等，1993）。从ENVI主菜单中选择SpectralSpectralLibrariesSpectralLibraryViewer。在SpectralLibraryInputFile对话框中，点击OpenFile按钮，从spec_lib/jpl_lib子目录中，选择jpl1.sli波谱库文件，点击OK。选择SelectInputFile区域中的jpl1.sli，点击OK。在SpectralLibraryViewer对话框中，选择Options→Edit(x,y)ScaleFactors，并在YDataMultiplier文本框中，输入值1.000，以匹配影像表观反射率范围（1-1000），点击OK。遥感原理与应用-4-在SpectralLibraryViewer对话框中，选择下列波谱名称，绘制它们的波谱曲线：ALUNITESO-4ABUDDINGTONITEFELDSTS-11ACALCITEC-3DKAOLINITEWELLORDEREDPS-1A得到如下波谱曲线绘制图：从绘制（plot）窗口菜单中，选择EditPlotParameters，自定义波谱曲线的绘制图。在PlotParameters对话框中，按下面的步骤进行：将Charsize减少为0.50。选择X-Axis单选按钮，然后将Range调整为1.90到2.45。还是选中X-Axis单选按钮，点击Left/RightMargins的箭头增量按钮，直到达到所需X方向的页边距。选择Y-Axis单选按钮，将AxisTitle改为“Reflectance”。还是选中Y-Axis单选按钮，点击Top/BottomMargins的箭头增量按钮，直到所需的Y方向的页边距。点击Apply，然后再点击Cancel。要显示波谱名称的图例，可以在绘制窗口中点击鼠标右键，从弹出的快捷菜单中选择PlotKey。将绘制窗口拖动到所需的大小，以容纳下波谱名称。在绘制窗口中，选择OptionsStackPlots，分别查看绘制的波谱曲线。绘制的遥感原理与应用-5-波谱曲线如下图所示：在#1SpectralProfile绘图窗口中点击鼠标右键，从弹出的快捷菜单中选择PlotKey，这将显示出波谱曲线的图例，该图例指出了所选像素的X和Y的像素坐标。在#1SpectralProfile绘图窗口中，选择OptionsNewWindow:Blank，打开一个新的绘图窗口。然后重新放置#1SpectralProfile绘图窗口和新绘图窗口的位置，使能够同时看到这两个绘图窗口。从主影像窗口菜单中，选择ToolsPixelLocator。使用PixelLocator对话框，定位到下列各点精确的像素位置上：位置点名称列（带偏移）行（带偏移）StonewallPlaya590570VarnishedTuff435555SilicaCap494514OpaliteZonewithAlunite531541StronglyArgillizedZonewithKaolinite502589BuddingtoniteZone448505Calcite260613在PixelLocator对话框中，输入像素的坐标，列（sample）590、行（line）570，使缩放指示矩形框移动到以这个像素为中心的影像地区，即StonewallPlaya地区，然后点击Apply，将矩形框移动到这个位置上。同时#1SpectralProfile绘图窗口将更新显示所选点的波谱曲线，其所对应的图例为：“X：590Y：570”。在新的绘图窗口中，点击鼠标右键，在弹出的快捷菜单中，选择PlotKey，打开显示了X和Y坐标位置的图例。在#1SpectralProfile绘图窗口中，使用鼠标左键，点击并按住图例“X：590Y：570”，将这个波谱曲线图例拖到新的绘图窗口中。对上表所列的每一个像素点重复上面的步骤，直到新的绘图窗口包含了所有的7种波谱曲线。在新的绘图窗口中，选择OptionsStackPlots。新的绘图窗口如下图所示：遥感原理与应用-6-将这些波谱曲线同已经提取的波谱库中的波谱曲线进行目视比较。注意到，实验室所测的波谱曲线同某些影像表观反射率波谱曲线之间存在着相似的形状及吸收特性。基于这些相似的特性，我们可以得到这样的结论：与alunite、buddingtonite、calcite和kaolinite实验室所测波谱曲线相似的影像波谱曲线是由上面相应的矿物质所构成的。从波谱库波谱曲线的绘图窗口中，将波谱曲线拖动到#1SpectralProfile绘图窗口中，进行直接比较。3鉴别波谱曲线使用SpectralAnalystTM来鉴别波谱曲线：ENVI提供了一个波谱匹配工具，它根据波谱库中的波谱曲线对影像中的波谱曲线进行评分。波谱分析使用多种方法产生一个在0到1之间的分数值，其中分数值1相当于完全匹配。从ENVI主菜单中，选择Spectral→SpectralAnalyst。点击SpectralAnalystInputSpectralLibrary对话框底部的OpenSpecLib按钮。选择进入usgs_min波谱库目录，选择usgs_min.sli波谱库文件，点击Open。usgs_min.sli文件出现在SpectralAnalystInputSpectralLibrary对话框中，选中该文件，点击OK。在EditIdentifyMethodsWeighting对话框中，点击OK。从主影像窗口菜单中，选择ToolsProfilesZProfile(Spectrum)。然后在#1SpectralProfile绘图窗口中，点击鼠标右键，从弹出的快捷菜单中，选择PlotKey，显示波谱曲线名称的图例。从主影像窗口菜单中，选择Tools→PixelLocator。在PixelLocator对话框中，输入像素点的坐标列502，行589，点击Apply。遥感原理与应用-7-在SpectralAnalyst对话框中，选择Options→EditMethodWeights。在EditIdentifyMethodsWeighting对话框中，为每一个Weight文本框输入值0.33，然后点击OK。在SpectralAnalyst对话框中，点击Apply。如果在#1SpectralProfile绘图窗口中显示了多条波谱曲线，那么将会出现一个波谱曲线列表。如果出现了该波谱曲线列表，那么就选择像素（列502，行589）所对应的那条波谱曲线。SpectralAnalyst对话框将如下图所示：SpectralAnalyst对话框，显示出了矿物质kaolinite的波谱曲线与像素（列502，行589）的波谱曲线最匹配(左图)SpectralAnalyst将根据波谱库中的波谱曲线对未知地物的波谱曲线进行评分。上图显示了对像素（列502，行589）的波谱曲线进行鉴别的结果。注意到列表的第一行显示kaolinite的波谱曲线评分最高。这个相对较高的分数值表明了该像素对应的地物与kaolinite最相似。用鼠标双击列表中的第一条波谱曲线。在同一绘图窗口中，绘制出未知地物的波谱曲线以及波谱库中的波谱曲线，以进行比较。该绘制图如下图所示：未知地物与波谱库中最匹配的kaolinite的波谱曲线的比较图，表明两者高度匹配遥感原理与应用-8-使用波谱分析（SpectralAnalyst），比较波谱曲线图，对刚才提取的影像波谱曲线进行鉴别，查看其为哪一种矿物质。当已经鉴别了几种矿物质后，就可以进入下一部分了。也可以选择将USGS波谱库中的波谱曲线同影像波谱曲线以及JPL波谱库中的波谱曲线进行比较。USGS波谱库所对应的波谱数据文件为usgs_min.sli。5定义感兴趣区感兴趣区（ROIs）用来提取像素集合的统计信息及其平均波谱曲线。可以在任何显示的影像上定义足够多的感兴趣区。从主影像窗口菜单中，选择OverlayRegionofInterest，打开ROITool对话框。5．1创建新的感兴趣区：在影像中，点击鼠标左键。在多边形顶点处，点击鼠标左键，绘制感兴趣区，或者点击并按住鼠标左键，移动鼠标，连续绘制感兴趣区。点击鼠标右键，封闭该多边形，完成感兴趣区的定义。再次点击鼠标右键锁定感兴趣区的位置。点击Stats按钮，计算感兴趣区的统计信息，绘制均值波谱曲线（白色），均值波谱上下各有一条的标准差曲线（绿色），以及最小和最大值的包络波谱曲线（红色），其包含了感