（一）Arcpy 提取面要素文件的质心

Lucky_云佳

已于 2024-05-04 10:41:06 修改

阅读量535

点赞数 21

文章标签：学习经验分享 arcgis

于 2024-04-19 11:15:23 首次发布

本文链接：https://blog.csdn.net/weixin_51924669/article/details/137957511

版权

文章目录

前言
01-Arcpy 提取面要素文件的质心
结语
- 总结
- 下文预告：

前言

本系列Arcpy实验包含以下内容：

提取江苏省县区级行政区域的质心，形成一个点类型shape文件，至少含有县区名称属性；
将2021年各区数据整理为excel数据，两列，一列为县区名称，一列为2021年县区常住人口；
编写一个脚本工具，将县区常住人口数据作为属性值附加到县区质心点数据上；
编写脚本工具，对江苏2021年常住人口数据进行插值，插值方法任选；
制作一个绘图模版，绘制江苏各地级市2021年人口分布图（注意是地级市，需要用地级数据分别裁剪人口分布数据）。

01-Arcpy 提取面要素文件的质心

【目的要求】：

提取江苏省县区级行政区域的质心，形成一个点类型shape文件，至少含有县区名称属性。

【主要思路】：

通过arcpy.da.SearchCursor()函数定位输入文件的属性字段，利用令牌"SHAPE@XY"获取质心坐标，并利用arcpy.CreateFeatureclass_management()函数创建点要素。

用到的重要函数：

arcpy.CreateFeatureclass_management()；
arcpy.AddField_management()；
arcpy.da.SearchCursor()；
arcpy.InsertCursor()

【实现步骤】

首先进行路径设置，设置输入输出文件的路径，输入文件是江苏省县区级行政区域country.shp。
接着使用arcpy.CreateFeatureclass_management()函数，创建点要素。
使用arcpy.AddField_management()函数，给创建的点要素添加质心坐标X、Y以及Name的字段。
使用arcpy.da.SearchCursor()函数，利用令牌"SHAPE@XY"获取质心坐标，然后将获取到的输入文件（江苏省县区级行政区域country.shp）的质心坐标与地区名称字段的属性值，使用for循环依次分别添加到提前创建好的3个空列表中。
使用arcpy.Point()函数，创建普通点。这里是给新生成的文件添加字段属性值，因此游标应使用arcpy.InsertCursor()函数依次进行后续的生成地理点并添加相关字段属性的操作。
使用for循环，并根据前面提取到的质心坐标依次生成质心点并添加字段属性值。
添加自定义工具箱Geo_Tools，在ArcToolbox上右键-【添加工具箱(A)】。
将上述代码稍作改进，主要是将路径设置等变量用arcpy.GetParameterAsText()函数来替换，然后在Geo_Tools自定义工具箱上右键-【添加(A)】-【脚本(S)】，创建自定义脚本工具。
在该脚本工具的【常规】中填写名称、标签、描述并勾选存储相对路径名；在【源】中填写.py的脚本代码路径；在【参数】中设置相关参数名称、数据类型、参数属性等。属性设置界面如下：
提取质心点脚本工具的界面如下所示：

【实验结果】

【脚本使用说明】：提取质心点脚本工具，需要输入待提取质心的面矢量要素，以及形成的质心文件需要的从面矢量要素的赋值字段属性，设置输出文件夹以及文件名（不写后缀），结果的质心文件可自由选择是否添加质心X和Y坐标的字段。
填写相关参数，在Arcmap中运行该脚本，所得结果如下：
在这里插入图片描述

源代码（仅供参考）

# -*- coding: GB2312 -*-
import arcpy

arcpy.env.overwriteOutput = True
# 00-路径设置
inputFile = arcpy.GetParameterAsText(0)  # 待求质心的面状矢量文件
Filed_Name = arcpy.GetParameterAsText(1)  # 将输入文件的某一字段属性赋值给生成的质心文件中
out_path = arcpy.GetParameterAsText(2)  # 生成的质心文件的保存路径（文件夹）
out_name = arcpy.GetParameterAsText(3)  # 生成的质心文件的文件名（不加后缀名）
Centroid_Flag = arcpy.GetParameterAsText(4)  # 质心标志位，1表示生成质心点，0表示不生成
pointFile = out_path + "/" + out_name + ".shp"  # 生成的质心文件的绝对路径

# 01-生成输入文件的质心坐标与地区名称并添加到创建的列表中
arcpy.AddMessage("读取输入文件……")
point_X = []
point_Y = []
point_Name = []
for row in arcpy.da.SearchCursor(inputFile, ["SHAPE@XY", Filed_Name]):
    point_X.append(row[0][0])  # 质心X坐标
    point_Y.append(row[0][1])  # 质心Y坐标
    point_Name.append(row[1])  # 名称属性

# 02-创建点要素
arcpy.AddMessage("创建点要素")
arcpy.CreateFeatureclass_management(out_path, out_name, "POINT")

# 03-给创建的点shp要素添加字段
arcpy.AddMessage("点要素添加字段")
fieldName = ['X', 'Y', 'Name']
if Centroid_Flag == 'true':
    arcpy.AddField_management(pointFile, fieldName[0], "DOUBLE")
    arcpy.AddField_management(pointFile, fieldName[1], "DOUBLE")
arcpy.AddField_management(pointFile, fieldName[2], "TEXT")
# arcpy.AddMessage(Filed_Name)

# 04-生成点要素，并添加属性值
arcpy.AddMessage("点要素添加属性值中……")
point = arcpy.Point()
rows = arcpy.InsertCursor(pointFile)
for i in range(len(point_X)):
    print(point_Name[i])
    row = rows.newRow()
    # 添加点坐标
    point.X = float(point_X[i])
    point.Y = float(point_Y[i])
    pointGeometry = arcpy.PointGeometry(point)
    row.shape = pointGeometry
    # 添加字段属性
    if Centroid_Flag == 'true':
        row.X = float(point_X[i])
        row.Y = float(point_Y[i])
    row.Name = point_Name[i]
    rows.insertRow(row)
del rows
arcpy.AddMessage("质心点生成成功！")

# 05-添加生成的点要素到当前图层
mxd = arcpy.mapping.MapDocument("CURRENT")
df = arcpy.mapping.ListDataFrames(mxd)[0]
addLayer = arcpy.mapping.Layer(pointFile)
arcpy.mapping.AddLayer(df, addLayer, "AUTO_ARRANGE")
# 图层刷新
arcpy.RefreshActiveView()
arcpy.RefreshTOC()
arcpy.AddWarning("生成文件已添加至当前图层！")