监督分类中常用的具体分类方法.doc

《监督分类中常用的具体分类方法.doc》由会员分享，可在线阅读，更多相关《监督分类中常用的具体分类方法.doc（22页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、.word.监视分类中常用的具体分类方法包括：最小距离分类法minimum distance classifier：最小距离分类法是用特征空间中的距离作为像元分类依据的。最小距离分类包括最小距离判别法和最近邻域分类法。最小距离判别法要求对遥感图像中每一个类别选一个具有代表意义的统计特征量均值，首先计算待分象元与类别之间的距离，然后将其归属于距离最小的一类。最近邻域分类法是上述方法在多波段遥感图像分类的推广。在多波段遥感图像分类中，每一类别具有多个统计特征量。最近邻域分类法首先计算待分象元到每一类中每一个统计特征量间的距离，这样，该象元到每一类都有几个距离值，取其中最小的一个距离作为该象元到该类

2、别的距离，最后比拟该待分象元到所有类别间的距离，将其归属于距离最小的一类。最小距离分类法原理简单，分类精度不高，但计算速度快，它可以在快速浏览分类概况中使用。多级切割分类法multi-level slice classifier：是根据设定在各轴上值域分割多维特征空间的分类方法。通过分割得到的多维长方体对应各分类类别。经过反复对定义的这些长方体的值域进展内外判断而完成各象元的分类。这种方法要求通过选取训练区详细了解分类类别总体的特征，并以较高的精度设定每个分类类别的光谱特征上限值和下限值，以便构成特征子空间。多级切割分类法要求训练区样本选择必须覆盖所有.word.的类型，在分类过程中，需要利用

3、待分类像元光谱特征值与各个类别特征子空间在每一维上的值域进展内外判断，检查其落入哪个类别特征子空间中，直到完成各像元的分类。多级分割法分类便于直观理解如何分割特征空间，以及待分类像元如何与分类类别相对应。由于分类中不需要复杂的计算，与其它监视分类方法比拟，具有速度快的特点。但多级分割法要求分割面总是与各特征轴正交，如果各类别在特征空间中呈现倾斜分布，就会产生分类误差。因此运用多级分割法分类前，需要先进展主成分分析，或采用其它方法对各轴进展相互独立的正交变换，然后进展多级分割。最大似然分类法(maximum likelihood classifier)：最大似然分类法是经常使用的监视分类方法之一

4、，它是通过求出每个像元对于各类别归属概率似然度likelihood，把该像元分到归属概率似然度最大的类别中去的方法。最大似然法假定训练区地物的光谱特征和自然界大局部随机现象一样，近似服从正态分布，利用训练区可求出均值、方差以及协方差等特征参数，从而可求出总体的先验概率密度函数。当总体分布不符合正态分布时，其分类可靠性将下降，这种情况下不宜采用最大似然分类法。最大似然分类法在多类别分类时，常采用统计学方法建立起一个判别函数集，然后根据这个判别函数集计算各待分象元的归.word.属概率似然度。这里，归属概率似然度是指：对于待分象元 x，它附属于分类类别 k 的后验概率。设从类别 k 中观测到 x

5、的条件概率为 P(x|k)，那么归属概率Lk 可表示为如下形式的判别函数：式中 P(k)为类别 k 的先验概率，它可以通过训练区来决定。此外，由于上式中分母和类别无关，在类别间比拟的时候可以忽略。最大似然分类必须知道总体的概率密度函数 P(x|k)。由于假定训练区地物的光谱特征和自然界大局部随机现象一样，近似服从正态分布(对一些非正态分布可以通过数学方法化为正态问题来处理)，因此通常可以假设总体的概率密率函数为多维正态分布，通过训练区，按最大似然度测定其平均值及方差、协方差。此时，像元 X 归为类别 k 的归属概率 Lk 表示如下这里省略了和类别无关的数据项。(6-9)式中：n:特征空间的维数

6、；P(k)：类别 k 的先验概率；Lk(x)：像元 X 归并到类别 k 的归属概率；.word.X：像元向量；k 类别 k 的平均向量(n 维列向量)；det：矩阵 A 的行列式k：类别 k 的方差、协方差矩(nn 矩阵).这里注意:各个类别的训练数据至少要为特征维数的 2 到 3 倍以上这样才能测定具有较高精度的均值及方差、协方差；如果 2 个以上的波段相关性强，那么方差协方差矩阵的逆矩阵可能不存在，或非常不稳定，在训练样本几乎都取一样值的均质性数据组时这种情况也会出现。此时，最好采用主成分变换，把维数压缩成仅剩下相互独立的波段，然后再求方差协方差矩阵；当总体分布不符合正态分布时，不适于采用

7、正态分布的假设为根底的最大似然分类法。当各类别的方差、协方差矩阵相等时，归属概率变成线性判别函数，如果类别的先验概率也一样，此时是根据欧氏距离建立的的线性判别函数，特别当协方差矩阵取为单位矩阵时，最大似然判别函数退化为采用欧氏距离建立的最小距离判别法。监视分类流程图监视分类流程图ErdasErdas 环境环境在专业遥感图像处理软件 Erdas 环境下，监视分类的流程图可以表示如下：.word.图 2-1 监视分类流程图监视分类考前须知监视分类考前须知1分类应从下往上，即每一地类应先细分为假设干小类，然后再依需要自下而上合并成大类。2每一类的训练区文件 aoi 与特征文件 sig 应该一一对应，

8、即每一类对应的训练区和特征文件都应该保存为一个单独的文件，以方便在调整训练区的时候进展修改。3精度检验后假设精度不符合要求，需要重新调整训练区，再次分类，直到精度满足要求为止。监视分类过程例如监视分类过程例如.word.1图 2-2 为某地 TM 遥感影像，432 波段假彩色合成。图 2-2 TM 影像432 波段合成2确定分类类别通过色调、纹理等图像特征，确定该区域分类类别为水体，植被和滩涂。各类分类特征如表 2-1 所示。表 2-1 分类特征.word.3为每一类选择训练区及特征文件1AOI 操作工具简介在 Viewer 窗口中选择“AOI“Tools，调出 AOIArea Of Inte

9、rest，感兴趣区浮开工具栏如图 2-3 所示。图 2-3 AOI 浮开工具栏其中较为常用的工具按钮为：.word.2特征文件操作工具简介特征文件从 AOI 区域中获得。使用“Erdas“Classifier“Signature Editor，调出特征文件编辑器，如图 2-4 所示。图 2-4 特征文件编辑器其中较为常用的工具为：翻开一个特征文件。新建一个特征文件/翻开新的特征文件编辑器。添加选中的 AOI 的特征到特征文件中。使用选中的 AOI特征替换当前特征。合并选中的特征文件中的特征到一个特征。.word.一般建立特征文件的步骤是，在Viewer 窗口中使用AOI 工具勾画感兴趣区，使用

10、把该AOI 区域中的特征添加到特征文件中。也可以选中多个 AOI 批量添加到特征文件中。2为各类别建立训练区文件和特征文件。把遥感影像放大到像元级，选择矩形AOI选择工具，根据建立的判读标识，在遥感影像上选择 AOI 区域，然后使用 依次添加图 2-5图 2-6.word.特征到特征文件中。注：作为例如，本例选择 3 个AOI 区域，且没有细分小类。选择完成的 AOI 区域和特征文件如图 2-5 和图 2-6 所示。分别保存为“水体.aoi和“水体.sig。在 Viewer 窗口中使用去除已经保存完毕的 AOI图层，重新选择其他类别的训练区，并建立新的特征文件。分别保存为“植被.aoi和“植被

11、.sig；“滩涂.aoi和“滩.word.涂.sig。3合并特征文件在各个类别的特征文件建立完毕后，需要合并成为一个总体特征文件。新建一个特征文件编辑器，选择翻开保存的“水体.sig文件。注意选择“Append添加把特征文件添加进来，而非“Replace替换。如图 2-7 所示。图 2-7 添加特征文件把水体特征文件添加进来之后，全部选中所有的特征，如图 2-8所示。.word.图 2-8 选中所有特征使用 工具，把选中的水体的所有特征合并为一个总体的水体特征，右单击“Class#列表，选择“Delete Selection删除原有特征如图 2-9 所示。图 2-9 删除原有特征重命名总体水体

12、特征的“Signature Name为“水体。如图2-10 所示。.word.图 2-10 总体水体特征如此添加其他两类进入，并合并成各自的总体特征，分别命名为“植被、“滩涂。并更改 Value 值为 1，2，3，并另存为 SaveAs“结果特征文件.sig如图 2-11 所示。图 2-11 结果特征文件4分类选择“Erdas“Classifier“Supervised Classification，在分类设置对话框中如图 2-12 设置。.word.图 2-12 监视分类设置在该对话框中，使用 输入待分类的图像“subset.img、分类特征文件“结果特征文件.sig并指定分类结果的保存路径

13、及名称，如“分类结果.img。分类方法选择“Maximum Likelihood最大似然，其余可以默认。点击“OK，系统将对原始影像依据指定的特征文件进展分类。运算完毕界面如图 2-13 示。.word.图 2-13 运算完成5分类结果分类的结果如图 2-14 所示。图 2-14 分类结果为了更好的表达分类结果，可以使用 Viewer 窗口中的“Raster“Attributes，更改“水体和“植被的显示颜色为蓝色RGB 为 0 0 1和绿色RGB 为 0 1 0，如图 2-15 示。.word.图 2-15 调整颜色调整颜色后的分类结果如图 2-16 所示。精度检验精度检验 同时翻开原始影像

14、和分类结果图，在任一幅图中单击右键，在弹出的菜单中选择“Geo.Link/Unlink，然后在另一幅图中单击左键，关联两幅影像。.word.使用“Erdas“Classifier“Accuracy Assessment，调出精度检验设置窗口。图 2-17 精度检验窗口 使用该窗口中“File“Open，翻开原始影像“Subset.img，调入内存。使用“View“Select Viewer，选择已经翻开的分类图，用以显示将要读取的点位信息。读入 GPS 测量的点。格式为标准的 txt 文本。文件格式化为 3 列，第一列存储 x 坐标，第二列存储 y 坐标，第三列存储类别代码即分类时指定的 Va

15、lue 值。如本例中存储的 GPS 点文件如表 2-3 所示表 2-3 GPS 点位491355.864486746.252.word.491070.454487008.523490754.194484941.222486997.564485905.443486797.004486707.683492096.394486615.111489118.864486815.963486434.454483151.611486920.424483028.193487375.534482665.641488069.784482449.652491070.454482657.931492266.09448

16、4439.821489604.834483167.041486881.854487216.791487984.934487085.661490769.614487116.511489905.674483483.303.word.489080.294483606.733487074.704483120.752文件中存储的坐标投影应与影像投影完全一致，如本例中影像投影为 UTM/Clark1866 N50。使用“Edit“Import User-defined Points，读入 GPS 点位文件。选项如图 2-18 所示。图 2-18 导入选项读取的结果如图 2-19 所示。.word.图 2-19 导入结果 输入各点位分类类别使用“Viewer“Show All，把读入的 GPS 点位在分类图中全部显示出来，逐一对照，在 Reference 列中输入分类影像的类别代码图 2-20。.word.图 2-20 输入代码 精度检验使用“Report“Accuracy report，分析分类精度情况图2-21。图 2-21 精度检验结果.word.假设精度符合要求，那么承受分类结果，假设不符合要求，那么重新分类。