遥感图像的监督分类与处理-赵文彪

1、遥感图像的监督分类与处理赵文彪杭州师范大学 遥感图像的监督分类和处理实验杭州师范大学遥感原理与应用实验报告题目：遥感图像的监督分类与处理实验姓名：赵文彪学号：2014212425班级：地信141学院：理学院1杭州师范大学 遥感图像的监督分类和处理实验1实验目的运用envi软件对自己家乡的遥感影像经行分类和分类后操作。2概述分类方法：监督分类和非监督分类监督分类一一从遥感数据中找到能够代表已知地面覆盖类型的均质 样本区域（训练样区），然后用这些已知区域的光谱特征（包括均值、标 准差、协方差矩阵和相关矩阵等）来训练分类算法，完成影像剩余部分的 地面覆盖制图（将训练样区外的每个像元划分到具有最大相似

2、性的类别 中）。非监督分类一一依据一些统计判别准则将具有相似光谱特征的像元 组分分为特定的光谱类；然后，再对这些光谱类进行标识并合并成信息类光谱特征空间同名地物点在JT同波段图像中亮度的观测量将构成一个多维的随机 向量X,称为光谱特征向量。而这些向量在直角坐标系中分布的情况为光 谱特征空间。同类地物在光谱特征空间中不可能是一个点，而是形成一个相对聚集 的点群。开同地物的点群在特征空间内一般具有不同的分布。特征点集群的分布情况：?理想情况：至少在一个子空间中可以相互区分?典型情况：任一子空间都有相互重叠，总的特征空间可以区分?一般情况：任一子空间都存在重叠现象监督分类，又称训练分类法，用被确认类

3、别的样本像元去识别其他未 知类别像元的过程。在分类前通过目视判读和野外调查，对遥感图像上 某些样区中影像地物的类别属性有了先验知识，对每一种类别选取一定数 量的训练样本，计算机计算每种训练样区的统计或其他信息，同时用这些 种子类别对判决西数逆行训练，使其符合于对各种子类别分类的要求，随 后用训练好的判决曲数去对其他待分数据逆行分类。使每个像元和训练样 本做比较，按不同的规则将其划分到和其最相似的样本类，以此完成对整 个图像的分类。3实验步骤遥感影像图的剪切用envi打开下载的遥感影像图，剪切出一个地貌信息丰富的区域（因为一景遥感影像太大，分类时间较长，故而采用剪切的方法，剪切一个地 貌丰富的遥

4、感影像图。既便于分类也使得分类种数不至于减小的太多）以下为剪切出来的遥感影像2杭州师范大学遥感图像的监督分类和处理实验类别定义/特征判别用envi打开剪切后的影像图，判断可分辨地物为：裸土地、村庄、林 地、河流湖泊、城市及城市道路和其他六类。样本选择在管理图层layer manager中、剪切区.dat图层右键选择 New region of interest (新建感兴趣区)打开 Region of Interest (ROI)Tool面板，开始选择 样本。在Region of Interest (ROI)Tool面板上，设置以下参数：ROI Name (感兴趣区名字)ROI Color (

5、感兴趣区颜色)默认ROIs绘制类型为多边形，在影像上辨别林地区域并单击鼠标左键 开始绘制多边形样本，一个多边形绘制结束后，双击鼠标左键或者点击鼠标右键，选择 Completeand Accept Polygon ,完成一个多边形样本的选择；完成裸土地、村庄、林地、河流湖泊、城市及城市道路和其他六类感 兴趣区的分类如图：3杭州师范大学遥感图像的监督分类和处理实验计算样本的可分离性。计算样本的可分离性。在 Region of Interest (ROI)Tool面板上，选择 OptionCompute ROI Separability ,在 Choose ROIsH板，将几类样本都打勾，点击OK;

6、表示各个样本类型的可分离性，用 Jeffries-Matusita,Transformed Divergence参数表示，这两个参数的值在 02.0之间，大于1.9 说明样本之间可分离性好，属于合格样本；小于 1.8,需要编辑样本或者 重新选择样本；小于1,考虑将两类样本合成一类样本。从如下图可知，分类良好，不需要合并。4杭州师范大学遥感图像的监督分类和处理实验在图层管理器中，选择 Region of interest ,点击右键，save as保存 为.xml格式的样本文件。5杭州师范大学遥感图像的监督分类和处理实验分类器选择平行管道法(Parallelepiped)平行管道法聚类分析：以地

7、物的光谱特性曲线为基础，假定同类地物 的光谱特性曲线相似作为判决的标准，设置一个相似阈值，这样，同类地 物在特征空间上表现为以特征曲线为中心，以相似阈值为半径的管子，此 即为所谓的“平行管道”最小距离(Minimum Distance ):利用训练样本数据计算出每一类的均值向量和标准差向量，然后以均 值向量作为该类在特征空间中的中心位置，计算输入图像中每个像元到各 类中心的距离，到哪一类中心的距离最小，该像元就归入到哪一类。最大似然分类算法(Maximum Likelihood )前面的分类器主要是根据距离测度进行判别而最大似然法的判别规则是基于概率最大似然法需要假设每个波段中各类训练数据都呈

8、正态分布，然后计算出概率密度函数神经网绚(Neural Net):指用计算机模拟人脑的结构，用许多小的处理单元模拟生物的神经元， 用算法实现人脑的识别、记忆、思考过程。支持向量机(Support Vector Machine):支持向量机分类 (Support Vector Machine或SVM)是一种建立在统计 学习理论(Statistical Learning Theory或SLT基础上的机器学习者斤法，可 以自动寻找那些对分类有较大区分能力的支持向量，由此构造出分类器， 有较好的推广性和较高的准确率。波谱角(Spectral Angle Mapper):它是在N维空间将像元不参照波谱逆

9、行匹配，通过计算波谱间的相似 度，之后对波谱之间相似度逆行角度的对比，较小的角度表示更大的相似 度。本实验选择支持向量机(Support Vector Machine)分类方法影像分类在 Toolbox/ Classification /Supervised Classification能找至U相应的分类方 法。6杭州师范大学遥感图像的监督分类和处理实验这里选择支持向量机分类方法。在 toolbox中选择/Classification/Supervised Classification/Support Vector Machine Classification, 选择待分类 影像，点击OK,按

10、照默认设置参数输出分类结果。开始分类7杭州师范大学遥感图像的监督分类和处理实验8杭州师范大学遥感图像的监督分类和处理实验分类还是比较精准的。分类后处理需要对初步的分类结果逆行一些处理，才能得到满足需求的分类结果， 这些处理过程就通常称为分类后处理。常用分类后处理通常包括：更改分类颜色、分类统计分析、小斑点处 理(类后处理)、栅矢转换等操作。包括更改类别颜色、分类后统计、小 斑坑处理、栅矢转换等。小斑块去除应用监督分类或者非监督分类以及决策树分类，分类结果中不可避免 地会产生一些面积较小的图斑，有必要对这些小图斑进行剔除或重新分类， 目前常用的方法有 Majority/Minority 分析、聚

11、类处理(clump)和过滤处理(Sieve)。以下对这三种方法融合处理Majority 和 Minority 分析Majority/Minority分析采用类似于卷积滤波的方法将较大类别的像元值归到该类中，定义一个变换核尺寸，主要分析( Majority Analysis)用变 换核中占主要地位（像元数最多）的像元类别代替中心像元的类别。如果 使用次要分析（Minority Analysis）,将用变换核中占次要地位的像元的类 别代替中心像元的类别。打开 Majority/Minority 分析工具，路径为 Toolbox/Classification/Post Classification/

12、Majority/Minority Analysis ,在弹出对话框中选择分类好的数 据，点击 OK;在 Majority/Minority Parameters 面板中，点击 Select All Items 选中所有的类别，其他参数按照默认即可，然后点击Choose按钮设置输出路径，点击OK执行操作。9杭州师范大学遥感图像的监督分类和处理实验10杭州师范大学遥感图像的监督分类和处理实验以下为处理后与处理前的对比，发现小斑块确实少了很多，图面变得光滑。聚类处理（Clump）聚类处理（clump）是将临近的类似分类区域聚类并逆行合并的一种 算法。分类图像经常缺少空间连续性（分类区域中斑点或洞的

13、存在）。低通滤波虽然可以用来平滑这些图像，但是类别信息常常会被临近类别的编 码干扰，聚类处理解决了这个问题。首先将被选的分类用一11杭州师范大学遥感图像的监督分类和处理实验个膨胀操作合并到一坑，然后用变换核对分类图像逆行腐蚀操作。打开聚类处理工具，路径为Toolbox/Classification/PostClassification/Clump Classes在弹出对话框中选择上一步处理好的数据，点击OK;在Clump Parameters面板中，点击Select All Items选中所有的类别， 其他参数按照默认即可，如下图所示。然后点击Choose按钮设置输出路径，点击OK执行操作。以

14、下为处理后与处理前的比较，发现图面变得更加光滑了。12杭州师范大学遥感图像的监督分类和处理实验过滤处理（Sieve）过滤处理（Sieve）解决分类图像中出现的孤岛问题。过滤处理使用斑 点分组方法来消除这些被隔离的分类像元。类别筛选方法通过分析周围的 4个或8个像元，判定一个像元是否不周围的像元同组。如果一类中被分 析的像元数少于输入的阈值，这些像元就会被从该类中删13杭州师范大学遥感图像的监督分类和处理实验除，删除的像元弼为未分类的像元（Unclassified） 0打开过滤处理工具，路径为Toolbox/Classification/PostClassification/Sieve Clas

15、ses在弹出对话框中选择要过滤处理的图像，点击 OK;在Sieve Parameters面板中，点击 Select All Items选中所有的类别， Group Min Threshold设置为5,其他参数按照默认即可，如下图所示。然 后点击Choose按钮设置输出路径，点击OK执行操作；以下为处理前与处理后的比较，发现小斑块更加少了，图面更加光滑 了。14杭州师范大学遥感图像的监督分类和处理实验精度验证对分类结果进行评价，确定分类的精度和可靠性。通常使用混淆矩阵用于精度验证，可以通过主菜单 -Classification-PostClassification-Confusion Matri

16、x来选择。真实参考源可以使用15杭州师范大学遥感图像的监督分类和处理实验两种方式：一是标准的分类图，二是选择的感兴趣区（验证样本区）。因为没有标准分类图和验证样本区，就用分类样本区来代替（如果用分类样本区来进行精度验证，毫无意义的）在 Toolbox 中， 选择 /Classification/Post Classification/Confusion MatrixUsing Ground Truth ROIs,选择分类结果，软件会根据分类代码自动匹配， 如不正确可以手动更改。点击OK后选择报表的表示方法（像素和百分比）， 点击OK,就可以得到精度报表。对混淆矩阵中的几项评价指标逆行说明：总体

17、分类精度：等于被正确分类的像元总和除以总像元数。被正确分类的像元数目沿 着混淆矩阵的对角线分布，总像元数等于所有真实参考源的像元总数。Kappa系数：它是通过把所有真实参考的像元总数 （N）乘以混淆矩阵对角线的和， 再减去某一类中真实参考像元数不该类中被分类像元总数之积之后，再除 以像元总数的平方减去某一类中真实参考像元总数不该类中被分类像元 总数之积对所有类别求和的结果。错分误差指被分为用户感兴趣的类，而实际属于另一类的像元，它显示在混淆 矩阵里面。漏分误差指本身属于地表真实分类，没有被分类器分到相应类别中的像元数。制图精度16杭州师范大学 遥感图像的监督分类和处理实验 是指分类器将整个 影

18、像的像元正确分为 A类的像元数（对角线值）不 A类真实参考总数（混 淆矩阵中A类列的总和）的比率。用户精度是指正确分到A类的像元总数（对角线值）不分类器将整个影像的像 元分为A类的像元总数（混淆矩阵中 A类行的总和）比率。分类统计分类统计（Class statistics可以基于分类结果计算源分类图像的统计信 息。基本统计包括：类别中的像元数、最小值、最大值、平均值以及类中 每个波段的标准差等。可以绘制每一类对应源分类图像像元值的最小值、 最大值、平均值以及标准差，还可以记每类的直方图，以及计算协方差矩 阵、相关矩阵、特征值和特征向量，并显示所有分类的总结记弓某。打开分类统计工具，路径为 To

19、olbox/Classification/PostClassification/Class Statistics，在弹出对话框中选择分类好的文件，点击OK；在Statistics Input File面板中，选择原始影像，点击OK;在弹出的ClassSelection面板中，点击Select All Items,统计所有分类的信息， 点击OK;在 Compute Statistics Parameters面板可以设置统计信息（全部勾选） 。Report Precision按钮可以设置输入精度，按默认即可。17杭州师范大学遥感图像的监督分类和处理实验18杭州师范大学遥感图像的监督分类和处理实验分

20、类叠加分类叠加（Overlay Classes功能，可以将分类结果的各种类别叠加在 一幅RGB彩色合成图或者灰度图像上，从而生成一幅RGB图像。打开分类结果和原始影像，（这里将原始影像的真彩色图像作为背景 图像）打开拉伸工具（Toolbox/Raster Management/Stretch Data）,在弹出的 对话框中选择原始影像，然后点击下方的Spectral Subset如下图所示）， 在弹出面板中选择波段1、2、3,点击OK;在 Data Stretching 面板中，Stretch Range 下的 min 为 2.0% max 为98.0%Output data range 下

21、min 为 0, max 为 255 data type 为 byte选择输出路径和文件名点击ok；19杭州师范大学遥感图像的监督分类和处理实验打开分类叠加工具，路径为Toolbox/Classification/PostClassification/Overlay Classes在打开的 Input OverlayRGB Image Input Bands面板中，R、G、B 分别 选择拉伸结果（上一步输出文件）的 and3、2、1,点击OK;20杭州师范大学遥感图像的监督分类和处理实验在ClassificationInput File面板中选择分类图像，点击 ok;在 Class Overlay to RGB Parameters面板中