上一篇《RSD的规则法遥感数据分类》介绍了分类方法,这里介绍一个用于识别乡村土路的例子。注意这里不是一个道路识别的项目,只是通过分类的方法提取得到乡村土路对象。
示例数据来自 //pan.baidu.com/s/1nvIJekT/RSD_Torial_Sample_Data\GFDM_0.5m融合数据土路分类.rsd。这是一个高分多模LiIFT融合后数据,仍然保留着8个多光谱波段,空间分辨率0.5m。见图1。
图1 GFDM 0.5m融合示例数据
乡村土路基本都是裸露的土壤自然形成或者由沙石堆填而成,可见光波段反射率较高目视轮廓明显,属于比较容易辨识的目标,分类应该比较容易。
1. 对象分割
在主窗口滚动鼠标滚轮直至出现黑色光标,右击,在弹出菜单选择 对象分析->基于图像特征的分割。分割得到图2a,再 对象分析->全局合并小对象 得到图2b。
图2 对象分割结果
全局合并小对象较慢,这一步不是必须的。道路的介质通常都比较均匀,具有很好的合并特性,合并能够得到一条条狭长的对象,容易分类。当然土路的对象特征有很多,不一定就非需要是狭长形状。
这里有保存的分割结果,可以从对象分析->打开对象分割结果文件 直接打开。 //pan.baidu.com/s/1nvIJekT\RSD_Torial_Sample_Data\GFDM_0.5m融合数据土路分类.seg
2. 建立规则
在规则文件中写入
规则文件 规则数 3 :
01. 平均值 大于 700.0
02. 欧氏距离 "8 1108.1252 739.6690 837.1032 601.7449 529.6680 522.1005 677.4884 564.7296 " 小于 0.2 与前述结论进行逻辑 AND
03. 形状指数 大于 3.0000 与前述结论进行逻辑 AND
三条规则中,“平均值”、“欧氏距离”和“形状指数”均为特征名,是不可以随便写的。从 面向对象分类->创建规则文件 菜单命令可以创建规范的规则文件。
第1条规则,通道计数平均值大于700是基于乡村土路反射率普遍偏高。各波段均值达700以上可能是道路。
第2条是和双引号中聚类核心的欧氏距离要小于0.2。括号中第1个数据8表示后面有8个数据,各自代表聚类核心每个波段的均值。这个是从乡村土路的样本数据计算而来的。
样本数据采集可以按下述步骤进行:
黑箭头光标右击待采集样本对象位置,在弹出菜单选则 新建分类样本类别。在对话框选确定后,打开 样本 浮动窗口,可见采集了一个对象样本。图3。
图3 采集样本
点击 样本 窗口的 1新建类别 使其点亮,然后在主窗口双击道路的对象,对象被点亮并添加到样本窗口里面,反复采集更多的对象。
采集足够的对象后,在弹出菜单选择 面向对象分类->设置节点规则 。然后在 选择样本类别 选中 新建样本。点击确定。出现一个对话框 ,见图4。在图4中勾选 欧氏距离 ,最右边的按钮被激活,点击,变成了 1 。
图4 选择特征
然后在本对话框点击确定。出现一个提示。图5
图5 保存规则
点击是,就将本条规则保存为一个规则文件。里面就有图5中的样本聚类核心数据。可以直接在这里编辑规则文件,也可以将聚类核心数据复制到上述规则文件中。
第3条规则是形状指数(新形状指数),可以理解为对象的长宽比的近似值。这个值越大表示对象越扁长。1最小,表示对象接近正方形或者圆形。
3. 分类
点击 面向对象分类->基于规则分类,提示需要打开规则文件,打开上面准备的有3条规则的规则文件。然后,在出现的对话框里选择基于先验知识和规则的类划分。还可以点击 分类颜色 按钮更改类别的颜色。点击确定,分类结果见图6。
图6 乡村土路分类结果
4. 分类后处理
分类都有一定误差,有的对象应该归类没有归类产生漏分,还有的把不是本类的对象分到本类产生误分,分类后编辑尽量修正这些漏分和误分。
上述分类创建了一个新层(图6中层2),在层列表窗口勾选本层,右击该层,在弹出菜单选择 设置为分割目标层 此时这一层的名字变成红色,这就是对象分析一章里面说的红层。这里红层作为接受来自其它层裁剪数据的缺省层。
勾选层1,选中漏分的对象,选择复制或者移动对象到目标层(红层),这样就能修正漏分。选择层2也就是分类结果,删除错分的对象,这样就修正了错分。这一部分参阅《RSD用户手册第15章对象分析》对象编辑内容。
5. 导出分类结果
勾选上述红层,主菜单选择 GRID层管理->GRID层转换为矢量数据->GRID转Shape文件 。出现一个文件保存的对话框,保存即可。再将该Shape文件添加到任务,见图7。
图7 乡村土路分类导出的Shape格式加载后结果
6. 小结
① 这是分类结果,进一步的道路结果可以继续修改分类结果,栅格或者矢量数据(Shape)。
② 基于对象的编辑将数据一块一块的移动,相比对像素操作,能缓解操作人员的焦虑感。
③ 分类方法简单,操作不落地。
请下载RSD专业版 SetupRsd3.2.5。下载地址:https://pan.baidu.com/s/1T-LBvaD_zVCwJsGf_hCyCg
详情加企鹅群136965427,在这里解答和讨论有关遥感数据处理和RSD平台的有关技术问题。