【ArcGIS技巧】用地块生成界址点去重、顺时针编号挂接DKBM属性

"农经权二轮延包上传到网签数据包，矢量数据不止地块矢量数据。还需要界址点、界址线等要素图层。这个数据不是很重要，但是得有，小编就写了个arxtools的工具来生成此图层，并且进行去重与编号字段、地块编码字段的生成，后续补充些字段即可，感兴趣的可以下载来使用。"

1、工具使用

在开始使用前，建议先对面要素图层进行简化，这样就防止拐点隔太近导致后续质检报错。打开arcmap，依次找到工具箱------系统工具箱------cartography tools.tbx------制图综合------简化面。（简化容差填1米）

插件界面很简单，打开工具箱脚本，输入的要素图层：1、地块图层（简化后的面要素图层），2、输出位置（选择一个文件夹目录）。注意：地块要素图层要有字段"DKBM"。

执行程序，在文件夹里会生成一个"JZD.shp"的要素图层，里面会包含字段"DKBM"、"JZDH"，如果一个点属于多个地块共有，用"/"隔开。

下载后记得加载py文件，不会的参考:【ArcGIS技巧】分享个判断是否基本农田的工具。

2、代码脚本

脚本的代码分享给大家，可以在此基础上修改(写代码的过程中，在pycharm和其他python解析器都能运行，在arcmap报错，后来找了一圈是有些变量需要设置成全局变量global)，有的编号会在前面加字母"J"：

python 复制代码

# -*- coding: utf-8 -*-
import arcpy
import math
from collections import OrderedDict
import sys
import os
# 设置默认编码为UTF-8
reload(sys)
sys.setdefaultencoding('utf-8')
# 输入输出路径
dk_shp = arcpy.GetParameterAsText(0)
jzd_folder = arcpy.GetParameterAsText(1)
output_fc = os.path.join(jzd_folder, u"JZD.shp")
# 创建空字典存储界址点信息（使用OrderedDict保持顺序）
global point_dict1
point_dict1 = OrderedDict()
# 处理地块要素
with arcpy.da.SearchCursor(dk_shp, ["DKBM", "SHAPE@"]) as cursor:
    global point_dict1
    for row_idx, row in enumerate(cursor):
        # try:  # 添加错误处理外层
        dkbm = row[0]
        geom = row[1]
        if geom.isMultipart or geom.type != "polygon":
            arcpy.AddMessage(u'地块 {} 包含多部件或非面要素，已跳过'.format(dkbm))
            continue
        extent = geom.extent
        northwest = arcpy.Point(extent.XMin, extent.YMax)
        # 优化顶点处理逻辑
        vertices = []
        for part in geom.getPart():
            part_vertices = []
            for pnt in part:
                if pnt:
                    part_vertices.append((pnt.X, pnt.Y))
            # 移除每个part的闭合点
            if len(part_vertices) > 1 and part_vertices[-1] == part_vertices[0]:
                part_vertices.pop()
            vertices.extend(part_vertices)
        if not vertices:
            arcpy.AddMessage(dkbm)
            continue
        # 顶点有效性检查
        if vertices == []:
            arcpy.AddMessage('CUO')
            continue
        # 添加计算保护
        try:
            distances = [
                math.hypot(v[0] - northwest.X, northwest.Y - v[1])
                for v in vertices
            ]
            min_dist = min(distances)
            start_idx = distances.index(min_dist)
        except:
            arcpy.AddMessage(dkbm)
            continue
        n = len(vertices)
        ordered_indices = [(start_idx + i) % n for i in range(n)]
        if not len(ordered_indices)==len(vertices):
            arcpy.AddMessage(dkbm)
            continue
        for point_number, vertex_idx in enumerate(ordered_indices, start=1):
            raw_x,raw_y= vertices[vertex_idx]
            hashed_coord = (round(raw_x,3),round(raw_y,3))
            if hashed_coord not in point_dict1:
                point_dict1[hashed_coord] = {'point':arcpy.Point(raw_x, raw_y),'dkbms':[],'jzdhs':[]}
            current_entry = point_dict1[hashed_coord]
            if dkbm not in current_entry['dkbms']:
                current_entry['dkbms'].append(dkbm)
                current_entry['jzdhs'].append(str(point_number))
# 创建界址点要素
global sr
sr = arcpy.Describe(dk_shp).spatialReference
arcpy.CreateFeatureclass_management(jzd_folder, "JZD.shp", "POINT", spatial_reference=sr)
arcpy.AddField_management(output_fc, "DKBM", "TEXT", field_length=255)
arcpy.AddField_management(output_fc, "JZDH", "TEXT", field_length=255)
# 写入要素
with arcpy.da.InsertCursor(output_fc, ["SHAPE@", "DKBM", "JZDH"]) as cursor:
    for key in point_dict1.values():
        dkbm_str = "/".join(key['dkbms'])
        jzdh_str = "/".join(key['jzdhs'])
        cursor.insertRow((
            arcpy.PointGeometry(key['point'], sr),
            dkbm_str,
            jzdh_str
        ))
arcpy.AddMessage(u"finish")