公益林小班与林地一张图小班融合后碎片处理思路_肖玲_20150421.doc

公益林小班与林地一张图小班融合后碎片处理思路 肖玲 2015.1.12 总体流程图 公益林小班处理 目的：只保留公益林区划界线，通过合并方式清除公益林界定界线。 将公益林小班按森林类别、事权等级、林地所有权、林木所有权、国家级公益林保护等级、工程类别用disolve工具进行融合： 融合后的图层命名为“公益林小班面522628_2014_Dissolve”。 兑现面积处理： 将“公益林小班面522628_2014”图层生成点图层，图层命名为 “公益林小班面522628_2014_FeatureToP” ： 打开Spatial Join工具，“公益林小班面522628_2014_Dissolve”图层 Spatial Join“公益林小班面522628_2014_FeatureToP”图层，新图层命名为“公益林小班面522628_2014_融合”。 首先清除点图层字段： 按“+”键新增“兑现面积”字段，用于将融合后的公益林图斑的兑现面积值更新为其范围内原公益林小班面积之和： 按右键菜单“Add Input Field”给新增“兑现面积”字段赋值： “公益林小班面522628_2014_融合”图层增加“gyl_id”字段（Long integer），且gyl_id = ObjectID：（注：操作视频无此步骤，需补充！） 叠加林地一张图小班 利用union工具将“公益林小班面522628_2014_融合”叠加到林地一张图小班上： 林地年度变更小班面522628 union 公益林小班面522628_2014_融合。 新生成的图层命名为“林地年度变更小班面522628_union”。 利用Multipart to Singlepart工具打散“林地年度变更小班面522628_union”图层中的多部件，新图层命名为“林地年度变更小班面522628_Union_Multip”。 “林地年度变更小班面522628_Union_Multip”图层增加DJH（字段类型：TEXT，文本型；长度：17），并按以下代码赋值： [XIAN] [XIANG] [CUN] [XIAO_BAN] “林地年度变更小班面522628_Union_Multip”图层增加CUN_ID（字段类型：TEXT，文本型；长度：12），并按以下代码赋值： [XIAN] [XIANG] [CUN] “林地年度变更小班面522628_Union_Multip”图斑分成三类：一是DJH有值且为17位；二是DJH有值且为12位，这部分图斑为非林地图斑，无小班号；三是DJH为空字符串，这部分图斑为县界外补充的公益林图斑。 对“双眼皮”碎斑进行筛选 计算图斑外接矩形面积。 “林地年度变更小班面522628_union_Multip”图层增加maxX（X最大值）、minX（X最小值）、maxY（Y最大值）、minY（Y最小值）、area_envelope（外接矩形面积）、area_percent（图斑图形面积占外接矩形面积比率）字段，字段类型均为double： 其中： maxX、minX、maxY、minY用以下方法赋值： 在工具箱打开Calculate Field工具： 针对maxX、minX、maxY、minY在Expression输入框分别录入以下PYTHON_9.3语句： maxX：!shape.extent.XMax! minX：!shape.extent.XMin! maxY：!shape.extent.YMax! minY：!shape.extent.YMin! 然后用以下公式计算外接矩形面积和图斑图形面积占外接矩形面积比率： area_envelope ＝（[maxX] - [minX]）*（[maxY] - [minY]） area_percent ＝ [SHAPE_area] / [area_envelope] 通过属性查找工具筛选出area_percent（图斑图形面积占外接矩形面积比率）＜0.3（该阈值可根据实际情况调整）的图斑，这种图斑一般是细长条形的，其中就包括了“双眼皮”图斑。 “林地年度变更小班面522628_union_Multip”图层增加isSelected字段，字段类型Short Integer，将选中的图斑赋值为1，值为1的图斑为待合并的图斑（注：其中包含了合理区划的细长形小班，需在接下来的操作中进行剔除）。 对“林地年度变更小班面522628”重复（1）～（2）步骤，目的在于进一步排除正常的细长条形图斑。 筛选出“林地年度变更小班面522628”area_percent（图斑图形面积占外接矩形面积比率）＜0