利用Python轻松在ArcGIS中进行图斑分类

已知数据：图斑数据文件

使用软件：ArcMap

要求：国有土地与待确认土地，国有土地与待确认重叠部分超过待确认土地面积的一半则判定待确认土地为国有。

下面开始操作：

1、 打开ArcMap，将国有土地与待确认添加到数据框内，点击地理处理—相交

输入要素待确认和国有土地，输出要素类位置选择我们自建的数据库，点击确定。

2、 打开ArcToolbox—分析工具—统计分析—汇总统计数据

输入表选择intersect，统计字段选择shape_area，统计类型选择SUM,分组字段选择fid_待确认。确认后我们得到一个表

3、 将表intersect_statistics与待确认属性表进行关联。右击待确认图层点击连接，连接选择某一表的属性，连接字段选择fid_待确认，

然后点击验证连接

验证无误后，右击待确认图层打开属性表就能看到连接到属性表的内容。

我们添加字段，名称写djmj_1（叠加面积），文件类型选择双精度

然后我们在表头右击—字段计算器djmj=sum_shape_area

接下来在图框中右击待确认图层，移除连接，移除后的属性表如下

4、 操作完成后接下来我们右击权属—字段计算器

开始判定地块权属

解析程序选择Python，在显示代码块前打勾，调出预逻辑脚本代码。

def qs(a,b):

if a >= 0.5* b:

return "国有"

else:

return "集体"

权属= qs( !djmj_1!, !SHAPE_Area!)

If前空格数为4个，return前空格数为8个

这里我们要判断的是待确认的叠加部分是否超过原有面积的0.5。

点击确定后权属的属性即可自动填入。

- END -