RSD的规则法遥感数据分类

这里笼统地说规则法可能会使大家一头雾水，那我再重新说一遍博文标题：一种基于先验知识的产生式规则遥感数据分类的一分类方法。这样说可能更糟糕了。 :)

先举个例子，我需要提取一个数据集的植被分布。你说植被区域NDVI值大于0.2，这就是先验知识。我处理数据，如果NDVI大于0.2就判定为是植被，这就是产生式规则的IF-THEN。然后将凡是NDVI大于0.2的区域提取出来，其它的一概不要，这就是一分类。通过这样一个简单的过程，我们就得到了植被的分类结果。

我们明确一下，上例这种分类方法在RSD里面重新赋予一个名字叫RSD规则法。或者就叫规则法，是为简便起见。不同于一些常用的最大似然函数、模糊隶属度等，规则法只得出二值结论。一分类只判断对象是或不是属于目标类别，而不是根据一个权值来判断对象到底归属哪一类别。

有同学可能会问，这不是一个基于阈值的图像分割吗？是的，不但如此，一分类还更像一个目标识别的问题。好在RSD是以解决问题为目标的，别纠结它是什么，还是多关注能做什么吧。

下面来举一个例子说明。

图1  GF2多光谱LiIFT光谱无损融合后0.8m数据

图1的示例数据来自：//pan.baidu.com/s/1nvIJekT/ RSD_Torial_Sample_Data\GF2_0.8m融合数据植被分类.rsd。这里的GF2多光谱数据是使用LiIFT光谱无损方式与PAN数据融合成0.8m数据。由于没有光谱失真问题，融合后数据仍能像融合前多光谱数据一样进行指数计算、分类、大气校正等处理。

这2个例子我们将密郁的山林和有裸地/受胁迫/作物种类不同的区域及进行分类。

1. 分割植被指数大于0.5的区域

第1步，先准备规则。编辑一个文本文件，里面写如下内容：

规则文件 规则数 1：

01 NDVI 大于 0.5

文件扩展名改为 .rule ，保存。RSD有专门编辑规则文件的菜单命令，不然用户不知道写什么特征名（这里的NDVI）才是正确的。这里用记事本写一个。

第2步，将上述示例数据进行对象分割。（这里有分割好的  ：//pan.baidu.com/s/1nvIJekT/ RSD_Torial_Sample_Data\GF2_0.8m融合数据植被分类.seg 可以滚动鼠标滚轮直至光标呈黑色，右击在弹出菜单选择 对象分析->打开对象分割结果文件）

第3步，从上一步的菜单选择 面向对象分类->基于规则分类 。点击打开第1步建立的规则文件，图2。

 图2 打开规则文件

第4步，分类。在图2点击确定。就开始分类，结果见图3。

图2 NDVI大于0.5的分类结果

可见不但有一个分类的掩膜，还新创建了一个分类结果的新层（层2）。隐藏（不显示）分类结果，隐藏层1（原始数据），可见层2（图2中蓝框）分类后数据如图3。

 图3 分类后的新层

2. 分割密郁的山林

勾选第2层也就是上一步分类的结果。打开对象浮动窗口，勾选NDVI和三个颜色贡献率的选项（图4中红框）。点击主窗口中山林的对象，可见都满足NDVI大于0.5，并且蓝贡献率都大于0.5。

图4 对密郁的林地对象特征求值

我们输入下列1行规则

01. 蓝贡献率 大于 0.4500

然后分类得到图5。

图5 分类得到的密郁的林地对象

3. 分割不同类型或者受胁迫的植被

稀疏的植被有裸露的土壤，或者受胁迫的作物红光反射率增加，可能表现在红光分量占比增加。如果我们直接对层1进行分类，规则可以写成

规则文件 规则数 2：

01 NDVI 大于 0.5

02. 红贡献率 大于 0.2400 与前述结论进行逻辑 AND

第2条规则表示与前述分类结果进行逻辑与，表示必须同时发生。得到的结果见图6。

 

图6 胁迫/不同类型植被（NDVI>0.5且高红分量）对象分割          左）分割前   右）分割后

这里的植被分类结果可能只是作物种类不同看起来偏黄。带有红边波段的数据对植被分类效果应该更好。

4. 小结

①. 先验知识可以是专家知识或者通过对象特征分析求值获得，并依此构造分类规则。对象特征的求算见《RSD用户使用手册-第16章》内容。

②. 基于光谱数据的特征可以使用公式求算，这样就可以构造更多几乎无限制的规则。例如：

...
;第2条规则是应用4波段和3波段计算NDVI，求得的NDVI值要求大于 0.5,与上一步是逻辑与关系
02.  公式 "(0.000334*_L1_B4-0.000246*_L1_B3)/(0.000334*_L1_B4+0.000246*_L1_B3)"  大于 0.5000 与前述结论进行逻辑 AND

这条规则是和计算NDVI同样的功能，其它还有一些基于对象位置关系形状大小等特征的计算手段。

③. 还是我常说的一句话，复杂的分类方法与分类精度相关不显著。和分类方法相比，高质量的对象分割可能更重要。

④. 见企鹅群 136965427 群文件《RSD用户手册16章-分类》。

详情加企鹅群136965427，在这里解答和讨论有关遥感数据处理和RSD平台的有关技术问题。

加企鹅758461012，原来的满了。