ENVI实习-高光谱遥感

1、文档大全文档大全高光谱遥感第三次实习实习任务：运用变换后的波段以及散点图工具提取端元运用变换后的波段以及纯净像元指数工具以及维可视化仪提取端元运用提取的端元进行分类和制图实习目标以及用时：学习运用软件进行纯净像元的提取方法教学方式：依据实习指导书进行实验，并完成实习报告使用器材：美国内华达的赤铜矿遥感数据，该数据已经经过大气校正，遥感软件具体实习过程本次实习主要内容：本章选用的实验数据是一幅经过校准的图像，处理的结果用于地质学应用，这主要是考虑到，到目前为止地质学研究仍然是高光谱遥感的主要应用领域之一。在主菜单下选择：，在打开的文件选择窗口中选择图像文件，点击打开图像：这是一幅经过校准的有个波

2、段的图像，图中显示的是将第18、319、320波7段分别赋红、绿、蓝合成的彩色图像。我们可以打开它的散点图观察一下。在主图像窗口中选择：o在随即弹出的波段选择窗口中任意选择两个波段，点击构成散点图。这里选择的是第17、217波3段。口回GflSca-tlerPict口回GflSca-tlerPict在这幅三点图上我们可以观察到，在由和波段组成的光谱特征空间中图像上的点明显地呈线状点云分布，说明这两个波段的相关性极强。遥感图像的某些波段之间往往存在着很高的相关性，直观上波段图像彼此很相似，从提取有用信息的角度考虑，有相当一部分数据是多余和重复的，解决这一问题的有效方法是进行特征提取和特征选择，去

3、相关和分离噪声。在多光谱遥感图像处理中，我们会采取旋转，但是相比之下，变化更适用于高光谱遥感数据。下面我们就用变换对图像进行处理。最低噪声分数（）变换用以确定图像数据的内在维度、隔离噪声以及降低后处理的计算要求。变换的本质就是两个叠置的主成分变换。第一次变换（基于估计的噪声协方差矩阵）用于分离和重新调节数据中的噪声。第一步导致了转换数据的噪声个体的变异和波段与波段的不相关。第二步是标准主成分变换。在主菜单下选择：变换完成后得到如下特征值曲线，其横坐标为变换后的波段数（之所以只有个波段是因为我在之前的参数设置中更改了输出波段数），纵坐标为特征值。把鼠标移到曲线上并点击左键，曲线上会出现一条以点击

4、位置为交叉点的十字，同时在窗口左下角显示当前的波段数和其特征值。特征值越高说明信息量越丰富回TJENFEigenvaluesFileEditOp+icnsFlot_FuiLcti-jnHelp10回TJENFEigenvaluesFileEditOp+icnsFlot_FuiLcti-jnHelp102030405C3,22,9653EigenvalueNumber另外，我们还可以利用散点图检查去相关的情况。以同样的方法打开散点图，但需要注意的是我们要用图像的第、波段，也就是信息最集中的两个波段构成散点图：26,21回G#1Sca-tlerPictFil已C1=seOptiotle26,21回

5、G#1Sca-tlerPictFil已C1=seOptiotleHelp可见，经过变换后的图像波段之间的相关性有效地降低了，并且出现了多个拐点，这些拐点就是我们要找的端元m在打开的散点图上点击鼠标中键会出现一个红色的小方块，在主图像窗口中对应这个小方块区域中的点同时呈现红色；在主图像窗口中点击鼠标左键，在散点图上对应点击位置的像素同时呈现红色。这一功能方便我们察看图上像素点与散点图上的像素点的相互对应位置。端元对应图像上的“纯”像元，是否能很好地提取它对于我们的分类是十分重要的。下面的操作就是利用处理后图像散点图选择端元生成样本区用以分类的过程。在打开的散点图上利用制图功能将点云拐角零散的几个

6、点圈起来，并且不同的用不同的颜色表示。同时在图上这些也显示了出来。回桀2ScatterPlotP匚口口宝一2.56fileCltEE回桀2ScatterPlotP匚口口宝一2.56fileCltEE在散点图窗口中选择：示如下窗口：将选择的区域输出为的，显在散点图窗口中选择：示如下窗口：将选择的区域输出为的，显在此为了便于区分，我用将每个样本区更名为其颜色的名称。选好了样本区后还要对样本区中的点进行训练。我接下来要进行的操作就是要通过链接散点图和剖面图察看样本区中点的光谱特征曲线，删除差别显著的点，达到训练样本区的目的。在散点图窗口中选择，在文件列表中选择反射数据，随即显示一个空白波段图。在散点

7、图窗口中点击鼠标右键，在空白的散点图窗口中就会出现当前的像素剖面图及其坐标值，将差别显著的点删去（用白色画样本区即可删除）。训练完成后，在散点图窗口中选择:提取各样本区的平均训练完成后，在散点图窗口中选择:提取各样本区的平均波谱曲线我们现在已经得到了七类地物，但还没有确定它们的种类。在此，我们可以运用的波谱分析功能来解决。波谱分析首先需要打开一个波谱库，然后将未知波谱与波谱库中的波谱进行匹配处理，并得出一系列匹配系数，系数越大就说明与这种地物越匹配。下面就进行具体的说明。在主菜单下选择：l在弹出的窗口中选择波谱库，在此我们选择（美国地质调查局）波谱库，点击，在弹出的“”窗口中设定计算参数后点击

8、，随即弹出波谱分析窗口，点击p在弹出的输入波段列表中选择我们要鉴别的波谱。选择一个待鉴别的波谱后波谱分析窗口中就会显示这一波谱与波谱库中的波谱匹配分析的结果，如下图：从结果来看，这一类显著匹配明矾石，所以我们基本上可以将这类定为“明矾石”。双击列表中的第一个波段名称绘制未知波段，波谱库波谱也会同时显示以作比较：依此方法将其余的几类逐一进行分析，确定其特征，并更改类别名称。通过分析我们发现最后两类显示了极其相似的特征，即中长石特征，于是我把这两类合并。在窗口中选择g弹出如下窗口，选择基准类及合并类后点击即完成合并。在保存了新的样本区文件后，就可以利用适当的监督分类方法进行分类了。以上介绍的是通过

9、散点图定义样本区的方法，对于高光谱遥感来说由于其数据特点，考虑到运用散点图分析效果可能会更好。此外在得到处理图像以后通常还用到其它一些处理方法，下面就一一介绍。像素纯度指数是一种在多光谱和高光谱图像中寻找波谱最纯的像元的方法。波谱最纯像元与混合的端元相对应。像素纯度指数通过迭代将散点图影射为一个随机单位向量。每次影射的极值像元被记录下来，并且每个像元被标记为极值的总次数也记下来。一幅“像素纯度图像”被建立，在这幅图像上，每个像素的值与像元被标记为极值的次数相对应。下图概括中的使用过程：在主菜单下选择：，选择前面得到的图像进行处理。经过次迭代后得到的图像如下所示：越亮的像素说明它被标记为极值的次

10、数越多相应地也越纯；相反，暗一些的图像纯度就低。在主图像窗口中选择：尝试不同的交互式拉伸以理解图像的直方图和数据分类。上图显示的是一个输入和一个输出直方图的比较窗口，在图中显示了当前的输入数据和各自拉伸的结果。两条垂直的星布线标志着当前拉伸的最小值和最大值。在窗口的底部列出了拉伸类型和直方图的来源。拖拽星布线的最大值和最小值，然后点击“”l拉伸将自动执行。这幅图像是前面的图像经过一万次迭代得到的结果，图像上像素点的值表明了它在迭代过程中有多少次作为极值像元被记录下来。这些数值显示了每个像素周围的数据云的局部突面程度以及每个像素和数据的突起外壳的亲近程度。简言之，值越高越接近散点图的数据源拐角，

11、拥有这些值的数据，其纯度比值低的像素高。零值像素是从未被作为极值的像素。下面要执行的操作是由图像生成样本区。在对话框中选择建立一个只包含拥有高值像素的。选择输入的文件，在弹出的对话框中输入最小极限值：。从下。从下点击“”，即生成包含在迭代过程中次以上作为极值的最纯像素的图中可以看出有298个符9合条件的点被提取出来，生成了红色的样本区。生成的包含最纯像素的位置，但却没有区分它们相应的端元。维空间观察仪可以帮我们解决这个问题。波谱可以被认为是散点图中的点（其中是波段数）。空间中的点坐标由个值组成，它们只是一个给定像元的每个波段中波谱辐射或反射值。这些点在空间中的分布可以用来估计波谱的端元数以及它

12、们的纯波谱信号数。维观察仪为维空间中选择端元提供了一个交互式工具。观察仪用于连接最小噪声分数转换和要定位、识别的纯净像元指数，并收集数据集中最纯的像元和极值波谱反应。观察仪允许数据在维空间中交互式旋转，选择像元组进行分类，以及聚集类，使其它类的选择更容易。选择的类可以输出到，并用作分类、不混溶和匹配的滤波技术的输入。在主菜单中选择,在弹出的对话框中选择前面处理好的文件，选择其前十个波段进行观察。这里默认使用刚才生成的只拥有高值的。如果有多个可选择的，它会让你选择用哪一个。点击后将弹出可以选择到波段的维控制对话框和维散点图窗口。选择前五个波段构成散点图。并选择控制对话框中的选项。随后在控制窗口中

13、点击进行旋转：这是一个维空间散点图的任意位置的动态显示功能，在这个模型中，小于输入波段数的任意多的波段都能同时被检查。在这一过程中我们可以切实感受到维数据在空间中的分布，可以确定数据真的是高维的，可见对于处理高光谱图像数据二维散点图是无法满足要求的。运用维观察器进行动态旋转观察，当出现感兴趣的位置时停止旋转，选择一种颜色将点云的拐角处圈起来作为一个，点击鼠标左键对进行定义，点击右键完成的建立。继续旋转，根据需要对类进行修改。以下的关于定义样本区、训练样本区、波谱分析法确定地物类别（测试波谱和波谱库中的波谱进行比较）等操作均与前面介绍的用散点图生成样本区的一系列操作类似，在此就不一一赘述了。下面只给出操作结果。通过散点图定义样本区并输出：回运用剖面图进行样本区的训练，之后输出样本区:回11ROIToolFileR口工_TypuOp+iorisHelpImageScrollE00mOffAvailableRegionsOtIntermit：n_DnDn_DDnDClassClas