遥感图像处理二（ENVI5.6 Classic）

1 实验目的和内容

1.1 实验目的

本次上机旨在继续深入了解ENVI软件的基本使用，并对提供的实验数据进行基本的图像分割和地物分类等操作并分析结果。

1.2 实验内容

1.2.1 图像分割

对教材示例数据“C7图像分割”中的风景图、兰花图和娃娃图分别进行图像分割操作。

1.2.2 图像分类

选影像数据进行图像分类操作，要求至少使用两种非监督分类方法五种监督分类方法，对比不同方法的分析效果并作精度评价。

1.3 实验流程

通过对实验内容进行分析，不难发现，本次实习主要包括图像分割和图像分类这两大部分，在该条件下，笔者总结了本次实验的操作逻辑，如图1.3-1所示。

图1.3-1 实习2操作流程图

2 图像操作处理过程

2.1 图像分割

2.1.1 风景图操作

首先，依次点击【File】-【Open Image File】在ENVI Classic中添加数据“地物与直方图DSCF0153”，并将其按照其原本的RGB三波段加载出来。如图2.1-1所示。

图2.1-1 添加图像数据

点击【Enhance】-【Interactive Stretching】即可观察RGB三波段的阈值，不难发现其直方图呈现双峰分布的特点，而这也表明了图像中主要有两类地图的样式，结合图2.1-1也不难发现，风景图主要是以颜色偏深色的地面和偏亮一些的天空为主。RGB波段直方图如图2.1-2所示。

图2.1-2 三波段直方图

接下来可以使用【Basic Tool】-【Band Math】进行波段的相关计算，这里输入公式为“(b1 gt 150)*(b2 gt 160)*(b3 gt 150)”，如图2.1-3所示。

图2.1-3 波段计算

由于这里的运算只是针对b1、b2和b3，所以需要将这里的b1、b2和b3指定到相应的风景图的RGB波段上。操作如图2.1-4所示。

图2.1-4 指定波段

得到波段计算的结果如图2.1-5所示，不难发现经过初步的波段运算，图像的前景和地面被很好的分离了，这样就形成了一个良好的掩膜文件。

图2.1-5 初步波段运算结果

接下来可以继续使用【Basic Tool】-【Band Math】进行波段计算，这里输入公式为“b4*(1-(b1 gt 150)*(b2 gt 160)*(b3 gt 150))”，如图2.1-6所示。

图2.1-6 波段运算

同上理，这里需要设置一下b1、b2、b3和b4波段分别为风景图的R、G、B三个波段以及通过【Map Variable to Input File】选择的风景图多波段影像。操作如图2.1-7所示。

图2.1-7 波段选择

最后，即可得到相应的分割结果，主要就是将地物给单独的提取出来，和GIS操作中的掩膜提取十分类似，不过这里是直接通过对栅格数据的直接操作而实现的。分割结果如图2.1-8所示。

图2.1-8 风景图分割结果

2.1.2 兰花图操作

兰花图同样位于“C7图像分割”文件夹中，依次点击【File】-【Open Image File】在ENVI Classic中添加数据“兰花图”，并将其按照其原本的RGB三波段加载出来。如图2.1-9所示。

图2.1-9 添加兰花图数据

通过观察兰花图像的直方图，不难发现其蓝光波段的峰值远大于红光和绿光波段，因此可以利用这种特性来针对兰花进行提取。兰花图像直方图如图2.1-10所示。

图2.1-10 兰花图像三波段直方图

接下来可以继续使用【Basic Tool】-【Band Math】进行波段计算，这里输入公式为“(b1 gt b2)*(b1 gt b3)”，如图2.1-11所示。

图2.1-11 波段运算

同上理，这里需要设置一下b1、b2、b3波段分别为风景图的B、R、G三个波段。操作如图2.1-12所示。

Tips：因为这里针对蓝光波段进行了一步操作，所有要将特征关注放在蓝光波段上，即b1、b2、b3波段不是直接赋RGB三个波段，而有一定顺序的调整。

图2.1-12 波段选择

将特征波段进行提取即可得到二值化的波段图像数据，这里就是兰花的掩膜图像，如图2.1-13所示。

图2.1-13 兰花掩膜图像

同2.1.1理，这里可以继续对兰花图像进行掩膜提取的操作，这里输入的波段运算公式为“b1*b2”，并将b1选择为上一步的掩膜图像，将b2选择为整个兰花图像的多波段图像。操作如图2.1-14所示。

图2.1-14 波段运算及赋值操作

最终即可得到提取后的兰花影像，如图2.1-15所示。

图2.1-15 兰花图图像分割结果

2.1.3 娃娃图操作

娃娃图也是位于“C7图像分割”文件夹中，依次点击【File】-【Open Image File】在ENVI Classic中添加数据“娃娃图”，并将其按照其原本的RGB三波段加载出来。如图2.1-16所示。

图2.1-16 加载娃娃图图像

同样，这里也依次点击【Enhance】-【Interactive Stretching】观察图像的RGB三波段的直方图分布。和2.1.1中的风景图与2.1.2中的兰花图不同的是，娃娃图中没有特定的特征颜色，其直方图分布的多峰特征也恰恰印证了这一点，娃娃图图像直方图如图2.1-16所示。

图2.1-16 娃娃图三波段直方图

这里以特征提取中心的娃娃抱桃子为例，和前两节中采用相同的方式，先提取掩膜在进行感兴趣区域的提取。使用【Basic Tool】-【Band Math】进行波段运算，将公式输入为“float(b1)/float(b2)”，并将相应的b1、b2分别设置为娃娃图的R和G波段。操作如图2.1-17所示。

图2.1-17 娃娃图波段运算

据此运算可以实现通过放大红光波段的影响从而实现掩膜提取，提取的结果如图2.1-18所示。

图2.1-18 黑白图像提取结果

通过如图2.1-18的提取结果可以很明显的观察到，当前的二值化提取并不理想，现在呈现的灰度值过于多样，比较难辨别，因此这里可以尝试继续通过设定阈值来提取符合预期的掩膜图像。提取后的娃娃图的图像直方图如图2.1-19所示。

图2.1-19 提取后直方图

接下来，继续使用【Band Math】工具，将公式设置为“b1 gt 0,98”，同时将这里的b1设定为刚刚波段计算得到的黑白图像。操作如图2.1-20所示。

图2.1-20 波段运算

此次运算即可得到较为理想的二值化的掩膜图像，如图2.1-21所示。

图2.1-21 掩膜图像结果

最后，和2.1.1节与2.1.2节中一样，将原始图像和掩膜图像叠加运算求得掩膜提取结果。操作如图2.1-22所示。

图2.1-22 波段运算

最终，可以得到娃娃图提取的结果，如图2.1-23所示。

图2.1-23 掩膜文件提取

2.1.4 图像分割操作结果分析

在本次实验中，通过对风景图、兰花图和娃娃图的图像分割操作，我们得到了一系列的分割结果。这些结果展示了不同图像分割技术在处理不同类型图像时的有效性和局限性。

1.风景图分割：通过波段运算和掩膜提取，我们成功地将天空和地面分离。这一过程表明，通过分析直方图和选择合适的波段阈值，可以有效地对图像中的主要地物进行分割。

2.兰花图分割：利用兰花图像中蓝光波段的峰值特性，我们通过波段计算和掩膜提取得到了兰花的掩膜图像。这一结果说明了特定波段的特征可以作为图像分割的重要依据。

3.娃娃图分割：娃娃图的分割过程较为复杂，因为没有特定的特征颜色。通过波段运算和阈值设定，我们得到了较为理想的二值化掩膜图像。这一过程展示了在缺乏明显特征时，如何通过多次迭代和参数调整来优化分割结果。

总体而言，图像分割操作结果分析表明，选择合适的分割方法和参数对于获得准确的分割结果至关重要。同时，这也强调了在实际应用中，需要根据图像的具体特征和需求来灵活选择和调整分割策略。

2.2 图像分类

2.2.1 非监督分类 – K-Means聚类

ENVI的所有操作基础都是基于多级窗口来实现的，这个原理其实应用过二次开发就很明了，因为数据都需要通过存储到内存才能进行调用处理，因此在进行图像分类处理也不例外，都需要先进行数据的加载。

这里先点击【File】-【Open Image File】在ENVI Classic中添加数据“”，并将其按照3、2、1波段赋值RGB真彩色加载出来。如图2.2-1所示。

图2.2-1 加载影像

接下里，依次点击【Classification】-【Unsupervised】-【K-Means】使用K-Means聚类对如上的Landsat7影像进行一步无监督分类。通过观察图2.2-1的影像地物的大致类别可以将其区分为五类进行分类。操作如图2.2-2所示。

图2.2-2 聚类分析操作

得到的分类结果如图2.2-3所示，不难发现，这里分类得到的影像呈现了一个比较明显的特征就是大体上可以将特征地物给区分开，但由于区分密度不一致，导致特征光谱曲线的结果不一致，同样由于植被的电磁波反射率与水体类似，导致了植被和水体呈现了相同的分类类别。

图2.2-3 K-Means聚类分析结果

显然，在如上的分析中，相当于是纯粹只通过特征光谱进行区分，而没有任何的附加区分条件，因此这里可以选择叠加一些条件从而辅助其更好的进行区分，这里的策略就是通过计算NDVI来辅助进行精度更高的分类。

Tips：这里使用归一化植被指数的主要原因就是因为在直接使用光谱进行分类时没有很好的分类植被和水体。

这里依次点击【Transform】-【NDVI】进行NDVI的计算，操作如图2.2-4所示。

图2.2-4 NDVI计算

NDVI计算得到的结果如图2.2-5所示。

图2.2-5 NDVI计算结果

接下来，可以使用【Layer Stacking】工具将NDVI计算的结果与原始影像进行合并处理，操作如图2.2-6所示。

图2.2-6 合并影像

接下来，就可以继续使用【K-Means】进行聚类分析了，最终得到的分类结果如图2.2-7所示。

图2.2-7 二次聚类分析结果

对比单纯利用光谱进行聚类分析（左边影像）和叠加NDVI后再进行据分类分析的结果（右边影像），可以很明显的察觉到，当叠加一些参量后，可以显著提高聚类分析的结果，尤其是在已知分类较为不理想的地物时。

Tips：这里笔者举例时叠加了NDVI，这主要是由于前面单一借助光谱进行分类时可以看到的植被分类不明显所致。如果时针对水体分类不明显则可以使用NDWI指数来辅助分类，这里的选择是弹性的，主要还是针对初步结果分析后所得的。

2.2.2 非监督分类 – ISODATA

除去K-Means，ISODATA也是一种优良的非监督分类方法。这里依次点击【Classification】-【Unsupervised】-【ISODATA】对2.2.1节中已经添加的原始影像进行分类。通过对地物进行观察，可以适当的将这里的动态分类设置为为“5-7”类，并且将主要依据像素的分类参数进行适当的调整。操作如图2.2-8所示。

图2.2-8 ISODATA分类

得到的分类结果如图2.2-9所示。

图2.2-9 ISODATA分类结果

这里可以将2.2.1节中的K-Means初步分类的结果与ISODATA分类的结果进行一步对比，如图2.2-10所示。

图2.2-10 两种非监督分类结果对比

通过对只基于光谱反射率分类的ISODATA分类结果（左边影像）和K-Means分类结果（右边影像）进行对比分析，很明显，基于ISODATA分类的效果是更优良的，其对于植被和水体的分类效果显然更加显著且贴合实际，减少了因为反射率相同而产生的问题。

2.2.3 监督分类 – 平行管道法

监督分类和无监督分类一个显著的区别就是无监督分类可以自行分类而无须干预，这在2.2.1节和2.2.2节中的操作中便可看到。而监督分类的特征就是需要人为进行干预，而这种干预手段就是预设训练集。

这里还是先加载Landsat7影像并以真彩色显示到窗口，接下来依次点击【Overlay】-【Region of Intrest】设置训练集，这里笔者依据对图像的理解，选择设置了五类地物，并依据“少量多次”的原则勾选处了训练集样本，如图2.2-11所示。

图2.2-11 设置训练集

Tips：这里笔者只是按照个人理解和个人感兴趣区域进行地物划分，显然可以继续细分类别设置更加精细的训练集，不过这样划分的训练时间成本一定会相应增加。

为了方便后续其它监督分类，可以使用【Save ROIs to File】将勾画的训练集进行保存，操作如图2.2-12所示。

图2.2-12 文件保存操作

同理，为了后续方便进行监督分类的精度评价，这里也勾画测试数据集并保存，操作如图2.2-13所示。

图2.2-13 保存测试集

Tips：这里需要注意的一个细节是测试集所勾选的区域不要和训练集相重合，以免对精度评价产生影响。

进行完数据集的处理工作后，就可以进行监督分类了，这里先以平行管道法进行示例。在确保训练集兴趣点是打开的情况下，依次点击【Classification】-【Supervised】-【Parallelepiped】进行监督分类，将五种分类类别全部勾选，即可进行监督分类，操作如图2.2-14所示。

图2.2-14 平行管道法监督分类

最后可得训练结果如图2.2-15所示。

图2.2-15 平行管道法监督分类结果

不难发现，和此前的无监督分类的两种方法相比，分类结果还是比较理想的。

2.2.4 监督分类 – 最小距离法

最小距离法也是监督分类方法中的一种，其主要是基于距离度量，将像元分配给与其最近的训练样本类别，而这里通常使用欧氏距离。

依次点击【Classification】-【Supervised】-【Minimum Distance】进行监督分类，将五种分类类别全部勾选，即可进行监督分类，操作如图2.2-16所示。

图2.2-16 最小距离法分类

分类结果如图2.2-17所示。

图2.2-17 最小距离法监督分类结果

2.2.5 监督分类 – 马氏距离法

同理，接下来使用马氏距离法进行分类。依次点击【Classification】-【Supervised】-【Mahalanobis Distance】进行监督分类，将五种分类类别全部勾选，即可进行监督分类，操作如图2.2-18所示。

图2.2-18 马氏距离法监督分类

分类结果如图2.2-19所示。

图2.2-19 马氏距离法监督分类结果

2.2.6 监督分类 – 最大似然法

接下来使用最大似然法进行分类。依次点击【Classification】-【Supervised】-【Maximum Likelihood】进行监督分类，将五种分类类别全部勾选，即可进行监督分类，操作如图2.2-20所示。

图2.2-20 最大似然法监督分类

分类结果如图2.2-21所示。

图2.2-21 最大似然法监督分类结果

2.2.7 监督分类 – 二进制编码法

接下来使用二进制编码法进行分类。依次点击【Classification】-【Supervised】-【Binary Encoding】进行监督分类，将五种分类类别全部勾选，即可进行监督分类，操作如图2.2-22所示。

图2.2-22 二进制编码法监督分类

分类结果如图2.2-23所示。

图2.2-23 二进制编码法监督分类结果

2.2.8 监督分类 – 支持向量机

接下来使用支持向量机进行分类。依次点击【Classification】-【Supervised】-【Support Vector Machine】进行监督分类，将五种分类类别全部勾选，即可进行监督分类，操作如图2.2-24所示。

图2.2-24 支持向量机监督分类

分类结果如图2.2-25所示。

图2.2-25 支持向量机监督分类结果

2.2.9 监督分类 – 神经网络

接下来使用神经网络进行分类。依次点击【Classification】-【Supervised】-【Neural Net】进行监督分类，将五种分类类别全部勾选，即可进行监督分类，操作如图2.2-26所示。

图2.2-26 神经网络监督分类

分类结果如图2.2-27所示。

图2.2-27 神经网络监督分类结果

这里可以统一将七种监督分析进行简单的对比，如图2.2-28所示。

图2.2-28 七种监督分类结果对比

通过对其中监督分类的结果进行对比，显然，这里的最小距离法的分类效果最为优良，而支持向量机、神经网络的分类则比较一般。

2.2.10 监督分类精度评估

在2.2.1节中笔者便已讲述过，这里可以使用测试数据集对监督分类后的结果进行精度评估，这里可以通过构建混淆矩阵来进行预测。在确保打开了测试数据集的情况下，依次点击【Classification】-【Post Classification】-【Generate Random Sample】-【Using Ground Truth ROIs】构建矩阵，并将选择分类图像选择为“Unclassified”，操作如图2.2-29所示。

图2.2-29 构建混淆矩阵

最后即可得到相应的混淆矩阵，不难发现应用二进制编码法的监督分类精度较低，只有66.7793%，显然这种监督分类是不理想的。

图2.2-30 混淆矩阵

2.2.11 图像分类结果分析

在图像分类部分，笔者使用了多种监督和非监督分类方法，并对结果进行了精度评估。以下是对分类结果的分析：

（1）非监督分类：

1.K-Means聚类：虽然能够大致区分地物类别，但由于植被和水体的光谱反射率相似，导致分类结果中这两类地物难以区分。

2.ISODATA：相较于K-Means，ISODATA分类结果更为优良，尤其是在植被和水体的分类上。

（2）监督分类：

平行管道法、最小距离法、马氏距离法、最大似然法、二进制编码法、支持向量机和神经网络：这些方法都需要预先设定的训练集，分类结果的准确性依赖于训练集的质量和代表性。其中，最小距离法的分类效果最为优良。

通过构建混淆矩阵，笔者对分类结果进行了精度评估。结果显示，二进制编码法的监督分类精度较低，这可能与其对地物特征的复杂性和多样性的敏感性有关。

综上所述，图像分类结果分析表明，不同的分类方法适用于不同的应用场景，且分类结果的准确性受到多种因素的影响，包括地物特征、训练集的选择和分类方法本身的特性。

3 实验总结

通过本次《遥感原理与应用》的上机实习二，我深入学习了ENVI软件在图像分割和分类方面的应用。实习内容包括对风景图、兰花图和娃娃图的图像分割操作，以及使用多种非监督和监督分类方法对影像数据进行分类，并进行了精度评估。

在图像分割方面，我掌握了如何利用直方图分析、波段运算和掩膜提取等技术来分离图像中的不同地物。这些技术在处理具有明显特征的图像时效果显著，但在处理特征不明显的图像时则需要更多的迭代和参数调整。

在图像分类方面，我学习了K-Means聚类、ISODATA、平行管道法、最小距离法、马氏距离法、最大似然法、二进制编码法、支持向量机和神经网络等多种分类方法。通过对比不同方法的分类结果，我认识到了每种方法的优缺点，并了解了如何根据具体的应用需求选择合适的分类方法。此外，通过构建混淆矩阵对分类结果进行精度评估，我学会了如何量化分类结果的准确性，并据此优化分类策略。

总体而言，本次实习不仅加深了我对遥感图像处理技术的理解，还提升了我的实际操作能力。通过亲身实践，我更加明确了在实际工作中如何根据不同的遥感影像和应用需求，选择合适的图像分割和分类方法，以及如何评估和优化分类结果。这些知识和技能将为我未来的学习和工作奠定坚实的基础。

4 参考资料

[1]   最全综述 | 图像分割算法-CSDN博客

[2]   机器学习入门（九）：非监督学习：5种聚类算法+2种评估模型_聚类算法 c#-CSDN博客

[3]   ENVI操作：监督分类_envi 按类别删除面roi-CSDN博客

[4]   基于ENVI的遥感影像的非监督分类_envi isodata-CSDN博客

图像分割操作文件位置：【免费】遥感图像处理分割源文件资源-CSDN文库

图像分类操作文件位置：：【免费】遥感图像处理分类源文件资源-CSDN文库