使用技巧-大数据量高分辨率影像的分块处理
标签:
分块规则集cart决策树高分辨率影像 |
分类: 易康整理 |
http://s3/bmiddle/002KgeW3ty71U73YbBM42&690
第一步,在Developer中对影像进行分块,将分块结果单独输出为影像。这里简单地将实验数据划分为左右两块,考虑到X方向是0-1206,Y方向是0-1130,因此只需要将X分为两段,这里设置的是0-600和600-1206,即有重复的像素,经过实验,如果设置为0-600和601-1206,那么分块的影像之间会出现1个像素的界线,如下图所示(导入到ArcGIS里即可看到)。
http://s11/mw690/002KgeW3ty71UjAsYJA7a&690
对影像的分块和输出,需要在新建工程的时候进行,如下图所示,首先将原始影像添加进来,然后点击【Subset
Selection】,弹出右边的【Subset
Selection】对话框;1表示影像预览界面中的红色边框,即分块的边界,2表示X方向设置当前分块影像的范围是0-600,当然如果获取右边的分块影像,X方向设置应该为600和1206,这样设置完毕分块影像的具体范围后,再选中3里面的【Use
geocoding for subset】和【Store subset in own
files】,前者使分块后影像继承原始影像的坐标信息,后者表示输出当前分块影像,点击OK就可以把当前的分块影像输出到一个文件夹,同时以当前分块影像新建一个Developer工程。
http://s3/mw690/002KgeW3ty71UkZh7GOc2&690
第二步,在其中一个工程里,建立分类体系和分类规则集。考虑到分块影像,分割的算法和参数设置一样的,分类方式一样的,分类体系是一样的,最终的结果输出格式等也是一样的,因此可以将一个工程里建立好的规则集ruleset,直接加载到另外分块影像所在的工程使用。这里以左边块状影像为例,对于大数据的影像分类而言,建议大家选用决策树分类的方式,而不是最邻近分类:(1)最邻近分类和决策树分类都需要选择样本,只不过是选择样本的方式不同,前者用【Select
Samples】,后者用【Manual
Editing】;(2)最邻近分类还需要去做分类特征集的优化,当分类样本数量较多,初始的分类特征集里面包含很多纹理特征时,速度会非常慢,而决策树分类的特征优选是放在训练阶段即train阶段后台完成的;(3)最邻近分类所有的类别使用相同的特征空间,特征空间如果包含若干个纹理特征,那么后面执行监督分类的速度会非常非常慢,十几甚至二十多个小时都很常见,而且分类效果也不一定满足要求,决策树分类是采用的CART算法,根据初始的特征空间能够去构建一棵二叉树,从而决定了类别提取的先后顺序、使用的特征和特征对应的阈值,从而避免了我们手动去逐个实验究竟哪个特征究竟采用哪个阈值来提取哪个类别的效果较好,而且分类精度明显好于最邻近分类。因此如下图所示,在一个工程的Process
Tree里建立的规则集,包含了分割、决策树分类、输出当前分类结果影像三个简单的规则,这个规则集直接可以导入到另外一个工程里进行执行,只不过需要重新选择一下样本而已。
http://s3/mw690/002KgeW3ty71UlRDAkO92&690
第三步,输出分类结果无非两种格式,一种是tif影像(栅格格式),一种是shp数据(矢量格式)。输出分类结果影像,可以在Process
Tree里添加结果输出的算法【export classification
view】,结果自动输出到工程所在文件夹里自动生成的【results】文件夹中;输出分类结果矢量数据,可以点击菜单项【Export】下面的【Export
Results】,弹出的界面如下图所示,选择要输出的类别,以及这些类别所带的属性(这里仅使用了类别名称,以方便在ArcGIS里做专题显示)。
http://s7/mw690/002KgeW3ty71Umfrg9M76&690
第四步,将输出的结果数据导入到ArcGIS的ArcMap软件打开,首先来看结果影像的叠加,分块边界处因为基本都属于同一个类别,因此并没有明显的分块线,这也是上面提到分块要考虑几个因素的原因。

加载中…