应用envi进行高光谱波谱分析和地物识别.pdf

航天星图科技（北京）有限公司 应用应用 ENVI 进行高光谱波谱分析和地物识别进行高光谱波谱分析和地物识别 在本章中， 主要介绍这些方法用 ENVI 如何实现， 以求在此过程中更好地熟悉和理解高光谱遥感图像的处理方法和步骤本章选用的实验数据是一幅经过校准的 AVIRIS 图像，处理的结果用于地质学应用，这主要是考虑到，到目前为止地质学研究仍然是高光谱遥感的主要应用领域之一最后，我对一幅相比之下空间分辨率更高的用于军事的高光谱图像进行了部分改进的分析操作，以便比较分类效果 在ENVI主菜单下选择： File > Open Image File， 在打开的文件选择窗口中选择图像文件cup95eff，点击 OK 打开图像： 这是一幅经过校准的有 50 个波段的 AVIRIS 图像，图中显示的是将第 183、193、207 波段分别赋红、绿、蓝合成的彩色图像 我们可以打开它的 2-D 散点图观察一下在主图像窗口中选择：Tools > 2-D Scatter Plots，在随即弹出的波段选择窗口中任意选择两个波段，点击 OK 构成 2-D 散点图这里选择的是第 172、173波段 在这幅 2-D 三点图上我们可以观察到，在由 172 和 173 波段组成的光谱特征空间中图像上的点明显地呈线状点云分布，说明这两个波段的相关性极强。

遥感图像的某些波段之间往往存在着很高的相关性，直观上波段图像彼此很相似，从提取有用信息的角度考虑，有相当一部分数据是多余和重复的，解决这一问题的有效方法是进行特征提取和特征选择，去相关和分离噪声 在多光谱遥感图像处理中，我们会采取 PC 旋转，但是相比之下，MNF 变化更适用于高光谱遥感数据下面我们就用 MNF 变换对图像进行处理 最低噪声分数（MNF）变换用以确定图像数据的内在维度、隔离噪声以及降低后处理的计算要求MNF 变换的本质就是两个叠置的主成分变换第一次变换（基于估计的噪声协方差矩阵）用于分离和重新调节数据中的噪声第一步导致了转换数据的噪声个体的变异和波段与波段的不相关第二步是标准主成分变换 在 ENVI 主菜单下选择：Transform > MNF Rotation > Forward MNF > Estimate Noise Statistics 航天星图科技（北京）有限公司 from Data 变换完成后得到如下 MNF 特征值曲线，其横坐标为变换后的波段数（之所以只有 25 个波段是因为我在之前的 MNF 参数设置中更改了输出波段数） ，纵坐标为特征值 把鼠标移到曲线上并点击左键，曲线上会出现一条以点击位置为交叉点的十字，同时在窗口左下角显示当前的波段数和其特征值。

特征值越高说明信息量越丰富 另外，我们还可以利用 2-D 散点图检查去相关的情况 以同样的方法打开 2-D 散点图，但需要注意的是我们要用 MNF 图像的第 1、2 波段，也就是信息最集中的两个波段构成散点图： 可见，经过 MNF 变换后的图像波段之间的相关性有效地降低了，并且出现了多个拐点，这些拐点就是我们要找的端元——Endmembers 在打开的散点图上点击鼠标中键会出现一个红色的小方块，在主图像窗口中对应这个小方块区域中的点同时呈现红色；在主图像窗口中点击鼠标左键，在散点图上对应点击位置的像素同时呈现红色这一功能方便我们察看图上像素点与散点图上的像素点的相互对应位置 端元对应图像上的“纯”像元，是否能很好地提取它对于我们的分类是十分重要的下面的操航天星图科技（北京）有限公司 作就是利用 MNF 处理后图像 2-D 散点图选择端元生成样本区用以分类的过程 在打开的 2-D 散点图上利用 ROI 制图功能将点云拐角零散的几个点圈起来，并且不同的 ROI 用不同的颜色表示同时在图上这些 ROIs 也显示了出来 在散点图窗口中选择：Options > Export All 将选择的区域输出为 ENVI 的 ROIs，显示如下窗口： 在此为了便于区分，我用 Edit 将每个样本区更名为其颜色的名称。

选好了样本区后还要对样本区中的点进行训练我接下来要进行的操作就是要通过链接 2-D 散点图和 Z-剖面图察看样本区中点的光谱特征曲线，删除差别显著的点，达到训练样本区的目的 在 2-D 散点图窗口中选择 Options > Z-Profile ，在文件列表中选择反射数据 CUP95EFF.INT，随即显示一个空白波段图在 2-D 散点图窗口中点击鼠标右键，在空白的 2-D 散点图窗口中就会出现当前的像素剖面图及其坐标值，将差别显著的点删去（用白色画样本区即可删除） 航天星图科技（北京）有限公司 训练完成后，在 2-D 散点图窗口中选择：Options > Mean All 提取各样本区的平均波谱曲线： 我们现在已经得到了七类地物，但还没有确定它们的种类在此，我们可以运用 ENVI 的波谱分析功能来解决波谱分析首先需要打开一个波谱库，然后将未知波谱与波谱库中的波谱进行匹配处理，并得出一系列匹配系数，系数越大就说明与这种地物越匹配下面就进行具体的说明 在 ENVI 主菜单下选择：Spectral > Spectral Analyst，在弹出的窗口中选择波谱库，在此我们选择USGS（美国地质调查局）波谱库，点击 OK，在弹出的“Edit Identify Methods Weighting”窗口中设定计算参数后点击 OK，随即弹出波谱分析窗口，点击 Apply，在弹出的输入波段列表中选择我们要鉴别的波谱。

选择一个待鉴别的波谱后波谱分析窗口中就会显示这一波谱与波谱库中的波谱匹配分析的结果，如下图： 航天星图科技（北京）有限公司 从结果来看，这一类显著匹配明矾石，所以我们基本上可以将这类定为“明矾石” 双击列表中的第一个波段名称绘制未知波段，波谱库波谱也会同时显示以作比较： 依此方法将其余的几类逐一进行分析，确定其特征，并更改类别名称通过分析我们发现最后两类显示了极其相似的特征，即中长石特征，于是我把这两类合并在 ROIs Tools 窗口中选择 Option > Merge Regions弹出如下窗口，选择基准类及合并类后点击 OK 即完成合并 航天星图科技（北京）有限公司 在保存了新的样本区文件后，就可以利用适当的监督分类方法进行分类了 以上介绍的是通过 2-D 散点图定义样本区的方法，对于高光谱遥感来说由于其数据特点，考虑到运用 n-D 散点图分析效果可能会更好 此外在得到 MNF 处理图像以后通常还用到其它一些处理方法，下面就一一介绍 像素纯度指数(PPI)是一种在多光谱和高光谱图像中寻找波谱最纯的像元的方法 波谱最纯像元与混合的端元相对应像素纯度指数通过迭代将 N-D 散点图影射为一个随机单位向量。

每次影射的极值像元被记录下来，并且每个像元被标记为极值的总次数也记下来一幅“像素纯度图像”被建立，在这幅图像上，每个像素的 DN 值与像元被标记为极值的次数相对应 下图概括 ENVI 中 PPI 的使用过程： 在 ENVI 主菜单下选择：Pixel Purity Index > [FAST]New Output Band，选择前面得到的 MNF 图像进行 PPI 处理经过 10000 次迭代后得到的 PPI 图像如下所示： 航天星图科技（北京）有限公司 越亮的像素说明它被标记为极值的次数越多相应地也越纯；相反，暗一些的图像纯度就低在主图像窗口中选择： Enhance > Interactive Stretching 尝试不同的交互式拉伸以理解 PPI 图像的直方图和数据分类 上图显示的是一个输入和一个输出直方图的比较窗口，在图中显示了当前的输入数据和各自拉伸的结果两条垂直的星布线标志着当前拉伸的最小值和最大值在窗口的底部列出了拉伸类型和直方图的来源拖拽星布线的最大值和最小值，然后点击“Apply” ，拉伸将自动执行 航天星图科技（北京）有限公司 这幅 PPI 图像是前面的 MNF 图像经过一万次迭代得到的结果，图像上像素点的值表明了它在迭代过程中有多少次作为极值像元被记录下来。

这些数值显示了每个像素周围的数据云的局部突面程度以及每个像素和数据的突起外壳的亲近程度简言之，值越高越接近 n-D 散点图的数据源拐角，拥有这些值的数据，其纯度比值低的像素高零值像素是从未被作为极值的像素 下面要执行的操作是由 PPI 图像生成样本区 在 ROIs Tool 对话框中选择 Options > Band Threshold to ROI 建立一个只包含拥有高 PPI 值像素的 ROI选择输入的 PPI 文件，在弹出的对话框中输入最小极限值： 点击“OK” ，即生成包含在迭代过程中 100 次以上作为极值的最纯像素的 ROI从下图中可以看出有2989 个符合条件的点被提取出来，生成了红色的样本区 航天星图科技（北京）有限公司 生成的 ROI 包含最纯像素的位置，但却没有区分它们相应的端元N 维空间观察仪可以帮我们解决这个问题 波谱可以被认为是n-D散点图中的点（其中n是波段数）n-D空间中的点坐标由n个值组成，它们只是一个给定像元的每个波段中波谱辐射或反射值这些点在n-D空间中的分布可以用来估计波谱的端元数以及它们的纯波谱信号数N维观察仪为N维空间中选择端元提供了一个交互式工具。

n-D观察仪用于连接最小噪声分数转换(MNF)和要定位、 识别的纯净像元指数，并收集数据集中最纯的像元和极值波谱反应 n-D观察仪允许数据在N维空间中交互式旋转，选择像元组进行分类，以及聚集类，使其它类的选择更容易选择的类可以输出到ROIs，并用作分类、不混溶和匹配的滤波技术的输入 在 ENVI 主菜单中选择 Spectral > n-Dimensional Visualizer > Visualize with New Data， 在弹出的对话框中选择前面处理好的 MNF 文件，选择其前十个波段进行观察这里默认使用刚才生成的只拥有高 PPI 值的 ROI如果有多个可选择的 ROIs，它会让你选择用哪一个点击 ROI 后将弹出可以选择1 到 10 波段的 N 维控制对话框和 N 维散点图窗口 选择前五个波段构成 n-D 散点图并选择 n-D 控制对话框中的 Options > Show Axes 选项随后在 n-D 控制窗口中点击 Start 进行旋转： 航天星图科技（北京）有限公司 这是一个 N 维空间散点图的任意位置的动态显示功能，在这个模型中，小于输入波段数的任意多的波段都能同时被检查。

在这一过程中我们可以切实感受到 N 维数据在空间中的分布，可以确定数据真的是高维的，可见对于处理高光谱图像数据二维散点图是无法满足要求的 运用 N 维观察器进行动态旋转观察，当出现感兴趣的位置时停止旋转，选择一种颜色将点云的拐角处圈起来作为一个 ROI，点击鼠标左键对 ROI 进行定义，点击右键完成 ROI 的建立继续旋转，根据需要对类进行修改 以下的关于定义样本区、训练样本区、波谱分析法确定地物类别（测试波谱和波谱库中的波谱进行比较）等操作均与前面介绍的用 2-D 散点图生成样本区的一系列操作类似，在此就不一一赘述了下面只给出操作结果 通过 n-D 散点图定义样本区并输出： 运用 Z-剖面图进行样本区的训练，之后输出样本区： 航天星图科技（北京）有限公司 利用波谱分析工具确定地物类别： 航天星图科技（北京）有限公司 在 n-D 控制窗口中选择 Options > Class Controls 弹出 n-D 分类控制窗口，在窗口中可以任意改变每一类的颜色，开启或关闭类，以及对类值片断的控制 点击窗口中的颜色块可以激活此类，并可对此类的显示符号进行修改，同时可以对其进行计算统计图表、平均波谱、划分类、清除或输出操作。

接下来我们要根据样本区进行分类在此我采用光谱角度制图法（SAM） 在 ENVI 主菜单下选择：。