公益林小班与林地一张图小班融合后碎片处理思路 - 肖玲 - 20150421

（5） “林地年度变更小班面522628_Union_Multip”图层增加CUN_ID（字段类型：

TEXT，文本型；长度：12），并按以下代码赋值： [XIAN]& [XIANG]& [CUN]

（6） “林地年度变更小班面522628_Union_Multip”图斑分成三类：一是DJH有值

且为17位；二是DJH有值且为12位，这部分图斑为非林地图斑，无小班号；三是DJH为空字符串，这部分图斑为县界外补充的公益林图斑。

4 对“双眼皮”碎斑进行筛选

（1） 计算图斑外接矩形面积。

“林地年度变更小班面522628_union_Multip”图层增加maxX（X最大值）、minX（X最小值）、maxY（Y最大值）、minY（Y最小值）、area_envelope（外接矩形面积）、area_percent（图斑图形面积占外接矩形面积比率）字段，字段类型均为double：

其中：

maxX、minX、maxY、minY用以下方法赋值：

a) 在工具箱打开Calculate Field工具：

b) 针对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]

（2） 通过属性查找工具筛选出area_percent（图斑图形面积占外接矩形面积比率）＜

0.3（该阈值可根据实际情况调整）的图斑，这种图斑一般是细长条形的，其中就包括了“双眼皮”图斑。

（3） “林地年度变更小班面522628_union_Multip”图层增加isSelected字段，字段

类型Short Integer，将选中的图斑赋值为1，值为1的图斑为待合并的图斑（注：其中包含了合理区划的细长形小班，需在接下来的操作中进行剔除）。