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

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

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

3、类法要求训练区样本选择必须覆盖所有 的类型， 在分类过程中， 需要利用待分类像元光谱特征值与各个 类别特征子空间在每一维上的值域进行内外判断， 检查其落入哪 个类别特征子空间中，直到完成各像元的分类。多级分割法分类便于直观理解如何分割特征空间，以及待分 类像元如何与分类类别相对应。由于分类中不需要复杂的计算， 与其它监督分类方法比较， 具有速度快的特点。 但多级分割法要 求分割面总是与各特征轴正交， 如果各类别在特征空间中呈现倾 斜分布，就会产生分类误差。因此运用多级分割法分类前，需要 先进行主成分分析， 或采用其它方法对各轴进行相互独立的正交 变换，然后进行多级分割。最大似然分类法 (max

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

5、象元的归属概率（似然度）。这里，归属概率（似然度）是指：对于待分 象元x，它从属于分类类别k的（后验）概率。设从类别k中观测到x的条件概率为P（x|k），则归属概率Lk 可表示为如下形式的判别函数：S =叫）=松）山伞）/迟砂）X（4）式中P（k）为类别k的先验概率，它可以通过训练区来 决定。此外，由于上式中分母和类别无关，在类别间比较的时候 可以忽略。最大似然分类必须知道总体的概率密度函数P（x|k）。由于假定训练区地物的光谱特征和自然界大部分随机现象一样，近似服从正态分布（对一些非正态分布可以通过数学方法化为正态问题来 处理），因此通常可以假设总体的概率密率函数为多维正态分布， 通过训练区

6、，按最大似然度测定其平均值及方差、协方差。此时，像元X归为类别k的归属概率Lk表示如下（这里省略了和类别 无关的数据项）。厶“2丹严阿-1x优-颅）辽：（兀-如（6-9）式中：n:特征空间的维数;P（k）:类别k的先验概率;Lk（x）:像元X归并到类别k的归属概率；X:像元向量；a k 类别k的平均向量（n维列向量）；det :矩阵A的行列式刀k :类别k的方差、协方差矩（n x n矩阵）.这里注意:各个类别的训练数据至少要为特征维数的2到3倍以上这样才能测定具有较高精度的均值及方差、协方差；如果2个以上的波段相关性强，那么方差协方差矩阵的逆矩阵可能不存 在，或非常不稳定，在训练样本几乎都取相

7、同值的均质性数据组 时这种情况也会出现。此时，最好采用主成分变换，把维数压缩 成仅剩下相互独立的波段，然后再求方差协方差矩阵；当总体分布不符合正态分布时，不适于采用正态分布的假设为基础的最大 似然分类法。当各类别的方差、协方差矩阵相等时，归属概率变成线性判别函数，如果类别的先验概率也相同， 此时是根据欧氏距离建立 的的线性判别函数，特别当协方差矩阵取为单位矩阵时，最大似然判别函数退化为采用欧氏距离建立的最小距离判别法。监督分类流程图（Erdas环境）在专业遥感图像处理软件 Erdas环境下，监督分类的流程图可以 表示如下：原始遥感彫椽1训经区1 (aoi)训练区2 (aoi)训练区-&

8、D训练区n Caoi)持征文件1(我Q特征文件2 Cig)持征文件持征交件总特征袁件(sig)选择分真器分类结果不符合裳求图2-1监督分类流程图监督分类注意事项(1) 分类应从下往上，即每一地类应先细分为若干小类，然后 再依需要自下而上合并成大类。(2) 每一类的训练区文件 aoi与特征文件sig应该一一对应， 即每一类对应的训练区和特征文件都应该保存为一个单独的文件，以方便在调整训练区的时候进行修改。(3) 精度检验后若精度不符合要求，需要重新调整训练区，再 次分类，直到精度满足要求为止。监督分类过程示例1.图2-2为某地TM遥感影像，432波段假彩色合成。图2-2TM影像（432波段合成）

9、2.确定分类类别通过色调、纹理等图像特征，确定该区域分类类别为水体，植被 和滩涂。各类分类特征如表 2-1所示。类别判别特從邑调呈暗蓝鱼，纹理变化较卜植祓色调呈现红色滩徐色调为白色.分布在水域附近表2-1分类特征AOI丸 o凶Z '罠+吐口£團令/&&b Clo$e ?图2-3 AOI浮动工具栏3.为每一类选择训练区及特征文件（1）AOI操作工具简介在Viewer窗口中选择“ AOI”f“Tools ”，调出 AOI （Area OfIn terest ，感兴趣区）浮动工具栏（如 图2-3所示）。其中较为常用的工具按钮为：按 钮功能简穴lx点选按钮,按下Shi

10、ft可以选择多个AOI区域。F9+1LJ框选按孙可以选择多个AOIE域.矩形AOI工具，用以绘制犯形的AOI区域。o箭圆丸01工島 用以绘制荊圆型AOI区域。任意AOIZ具，绘制任意形狀的AOI区域。AOI剪切删除工具，用于剪切/删除选中的AOI区域。(2) 特征文件操作工具简介 特征文件从AOI区域中获得。使用“ Erdas” f“ Classifierf“ Sig nature Editor ”，调出特征文件编辑器，如图2-4所示Signature Editor (No File)冋皀 Edit View Evaluate Feature Classify Help上1乂由E j已S A

11、(A A图2-4特征文件编辑器其中较为常用的工具为:打开一个特征文件。征文件编辑器。新建一个特征文件/打开新的特S3添加选中的aoi的特征到特征文件中。I使用选中的AOI特征替换当前特征。=L合并选中的特征文件中的特征到一个特征比f三1口1划般建立特 征文件的步骤 是，在Viewer窗 口中使用AOI工 具勾画感兴趣 区，使用03把 该AOI区域中的 特征添加到特征 文件中。也可以 选中多个AOI批 量添加到特征文 件中。（2）为各类别建 立训练区文件和 特征文件。把遥感影像放大到像元级，窃凸'ED 岡雪0 蚩兼凶凸 +乓图2-5啦 Signature Editor (No File)

12、-inlxiF3e Edit Vi&w Evaluate FeatureClassify Help汾口 4 I斗s V LA A图2-6选择矩形AOI选 择工具，根据建立的判读标识， 在遥感影像上选 择AOI区域，然 后使用依次添 加特征到特征文 件中。（注：作 为示例，本例选 择3个AOI区域， 且没有细分小 类。）选择完成 的AOI区域和特 征文件如图2-5 和图2-6所示。分别保存为“水体.aoi ”和 “水体.sig ”。在Viewer窗 口中使用去 除已经保存完毕 的AOI图层，重 新选择其他类别的训练区，并建立新的特征文件。分别保存为“植被.aoi ”和“植被.sig ”；

13、“滩涂.aoi ”和“滩涂.sig ”。（3）合并特征文件在各个类别的特征文件建立完毕后，需要合并成为一个总体特征 文件。新建一个特征文件编辑器，选择曰打开保存的“水体.sig ”文件。注意选择“ Append'（添加）把特征文件添加进来，而非 “Replace” （替换）。如图2-7所示。图2-7 添加特征文件把水体特征文件添加进来之后，全部选中所有的特征，如图2-8 ID所示。楊口斗丸送0 Lt Class itl?JSiqntuie NameColorFeeGreenEIufValueOrder02770.3120.3C611活EG.2770.305CJI2ZClasi 30.2

14、810.3130.J1433Class 40.2730.30B0.OT44Class： 50.2840.331C.4355Class 60.2B30.349C.36SG&Chss 70.2910.70 3S2?7Cla«£： 90.307o.doe0.454日32弓ignatur亡 E由tm*水File Edit View Evaluate Feature Classfy Help图2-8选中所有特征使用 工具，把选中的水体的所有特征合并为一个总体的水体特征，右单击“ Class# ”列表，选择“ Delete Selection ”删除原有特征如图2-9所示Sig

15、nature Editor （水作,sig）冋箕1Fde Eck 决賊 EvaluateFeature Gawify He4pD *U +-* =USLAACl*$s tt !| >Swa lure NameColorr RedSreftfi8lue=ValueOfdn _I0.27703120.3061P 1ROrW SewctiGn1 |0.2770 3050.3C1225ele:t rjone|Q 2610 3130l3U33Select Al02?80.30644 1Irwwihi I|必40.3310l34655in戶jlAr kinn0.2880 343D.3CC66CZI0

16、.2S10 357a 38277厂士 jmaIh,QX70 408DJ54pH«eLrlCftris. >A11D.28503340.345g上1tr 112J0t0.图2-9删除原有特征重命名总体水体特征的“ Sig nature Name ”为“水体”。如图2-10所示。Eb Siqn-ature EditorFile Edit View Eyaiuate Feature Classify Help 金口 7 * 三S A as售 It H、| Sjqnatuie Mare | 5lor | RedEi 匕tiBlue | x/atue | dider | Count |

17、氏兰图2-10 总体水体特征如此添加其他两类进入，并合并成各自的总体特征，分别命名为“植被”、“滩涂”。并更改Value值为1,2, 3,并另存为(SaveAs)“结果特征文件.sig ”如图2-11所示。粉.Signature Edit or （结果特征直件予ig） ID刈FIs EditView Evaluate feature 匚Help詣0卜 L, =US 7遍1T Class itSignature N-anneColorRedGreenBlueValuerder IJ1028513340.345TiQ8481.0001 00016I>0.0770.477 555§2

18、4图2-11结果特征文件(4)分类f“ Supervised选择“ Erdas ” f“ ClassifierClassification ”，在分类设置对话框中如图2-12设置2S1Inpi* R-aster File' p.img)subsel.mnginput Sonatuie Fite： (Ksig)皓果特征立件.比Cl-assied 冋 ecC-img)Distance FileAthibute Options.Alena me: A imgFussy Classifies buntiet Clares % Piwe>Decision Rule;H(?n-p-3rame

19、tric Rule;parametric Fiu!everlap) Huie:IJrcJasiled Rule:|Paiametiic Rule21Parametric Ruh：曲imuiTi Likeliho«jK.Classify zeiosUse ProbabilitiesBatchA3I.CancelHelp I图2-12监督分类设置在该对话框中，使用 输入待分类的图像“”、分类特征文件“结 果特征文件.sig ”并指定分类结果的保存路径及名称，如“分类结果.img ”。分类方法选择“ Maximum Likelihood ” （最大似 然），其余可以默认。点击“ OK

20、9;，系统将对原始影像依据指定的特征文件进行分类。 运算完毕界面如图2-13示。图2-13 运算完成（5）分类结果 分类的结果如图2-14所示噩File UtAty View AOI RasterEl H 3；Help图2-14分类结果为了更好的表达分类结果，可以使用Viewer窗口中的的显，如图“ Raster ” f“ Attributes ”，更改“水体”和“植被” 示颜色为蓝色（RGB为0 0 1 ）和绿色（RGB为 0 1 0 ） 2-15 示。|业 Raster Attribute Editor -分类结果»img(;L4yr.i) IFile Edit HelpD H

21、芍匚电IO)Layer Number卩日2厂 -ini ,1Eile Ut0ity Vieu A 囚 festerHelp爲If' 2滋皆乡出41425.00, 44B2913.00 (UTM ； Chrkf 抵图2-15 调整颜色调整颜色后的分类结果如图2-16所示。血 Viewer #2 :券精度检验 同时打开原始影像和分类结果图，在任一幅图中单击右键，在弹出的菜单中选择“ Geo. L ink/Un li nk ”，然后在另一幅图中单击左键，关联两幅影像。Assessment”，调出精度检验设置窗口Accuracy 使用“ Erdas ” “ Classifier图2-17 精度

22、检验窗口 使用该窗口中“ File ” “ Open”，打开原始影像“”，调 入内存。 使用“ View” f“ Select Viewer”，选择已经打开的分类 图，用以显示将要读取的点位信息。 读入GPS测量的点。格式为标准的txt文本。文件格式化为3列，第一列存储x坐标，第二列存储y坐标，第 三列存储类别代码（即分类时指定的 Value值）。如本例中存储 的GPS点文件如表2-3所示表2-3 GPS点位3233131312111111332文件中存储的坐标投影应与影像投影完全一致，如本例中影像投 影为 UTM/Clark1866 N50 使用“ Edit ”f“ Import User-

23、defined Points ”，读入 GPS点位文件。选项如图2-18所示。图2-18 导入选项读取的结果如图2-19所示File Edir View Report H»?lpP*i卅 ttNameXYClassR efareoce1|ID 41491355.060446674625022ID轮491070.4504487008.520331014343075419044G4341.22024lOtt+486997.5601405905.44035IDtt5486797.0004466707.66036IDttB492096.3904486615.11017IDI17489118.

24、SB04466S15.9603gID惬4S6434.450448311.610T9lDtt94869204204463028.1903101DB10487375.53044S26S5.640111iDttn4S9059.78044E2449.&50212IDW2491070.4504482G57.930113lDtt13492266.0904484439.320111IDti14489604.83044631 S7.04D115IDB1549Saai.85O448721 G.7901154B7994 93D4467035.560117IDH1749O709.G1O44G711G510119ID#1E489905.6