写一个arcgis工具箱 Python 2.7 版本脚本:实现arcgis空间连接多对一,可以根据设置面积大于多少的值进行多对一空间连接吗,多对一的赋值的字段用“\”隔开。
时间: 2023-12-06 09:45:18 浏览: 113
下面是一个实现多对一空间连接的 ArcGIS 工具箱 Python 脚本,可以根据设置面积大小的值进行多对一空间连接,并且可以设置多对一赋值的字段用“\”隔开。
```python
# -*- coding: utf-8 -*-
import arcpy
# 定义工具箱参数
input_features = arcpy.Parameter(
name="input_features",
displayName="输入要素",
datatype="DEFeatureClass",
parameterType="Required",
direction="Input")
join_features = arcpy.Parameter(
name="join_features",
displayName="连接要素",
datatype="DEFeatureClass",
parameterType="Required",
direction="Input")
join_type = arcpy.Parameter(
name="join_type",
displayName="连接类型",
datatype="GPString",
parameterType="Required",
direction="Input")
join_type.filter.type = "ValueList"
join_type.filter.list = ["ONE_TO_ONE", "ONE_TO_MANY", "MANY_TO_ONE", "MANY_TO_MANY"]
join_fields = arcpy.Parameter(
name="join_fields",
displayName="连接字段",
datatype="Field",
parameterType="Required",
direction="Input",
multiValue=True)
output_features = arcpy.Parameter(
name="output_features",
displayName="输出要素",
datatype="DEFeatureClass",
parameterType="Required",
direction="Output")
area_threshold = arcpy.Parameter(
name="area_threshold",
displayName="面积阈值",
datatype="GPLong",
parameterType="Required",
direction="Input")
join_value_fields = arcpy.Parameter(
name="join_value_fields",
displayName="连接赋值字段",
datatype="GPString",
parameterType="Required",
direction="Input")
# 设置连接字段
join_fields.parameterDependencies = [join_features]
join_fields.filter.list = [field.name for field in arcpy.ListFields(join_features)]
# 设置连接赋值字段
join_value_fields.value = "field1\\field2"
# 将工具箱参数添加到参数列表
parameters = [input_features, join_features, join_type, join_fields, output_features, area_threshold, join_value_fields]
# 定义脚本函数
def multiple_join(input_features, join_features, join_type, join_fields, output_features, area_threshold, join_value_fields):
# 多对一空间连接
arcpy.SpatialJoin_analysis(input_features, join_features, output_features, join_type, "KEEP_ALL", join_fields)
# 计算面积
arcpy.AddField_management(output_features, "SHAPE_Area", "DOUBLE")
arcpy.CalculateGeometryAttributes_management(output_features, "SHAPE_Area")
# 根据面积阈值进行筛选
expression = "SHAPE_Area >= {}".format(area_threshold)
arcpy.Select_analysis(output_features, "selected_features", expression)
# 对选择的要素进行多对一赋值
join_value_fields_list = join_value_fields.split("\\")
join_value_fields_count = len(join_value_fields_list)
if join_value_fields_count == 1:
# 多对一赋值字段只有一个
arcpy.CalculateField_management("selected_features", join_value_fields_list[0], "!{}!".format(join_fields[0].name), "PYTHON")
else:
# 多对一赋值字段有多个
for i in range(join_value_fields_count):
arcpy.CalculateField_management("selected_features", join_value_fields_list[i], "!{}!".format(join_fields[i].name), "PYTHON")
# 执行脚本函数
if __name__ == '__main__':
arcpy.AddMessage("开始执行多对一空间连接...")
multiple_join(*parameters)
arcpy.AddMessage("多对一空间连接执行完成!")
```
使用时,将上述脚本保存为 .py 文件,并将其添加到 ArcGIS 工具箱中,即可在 ArcGIS 中使用该工具箱进行多对一空间连接操作。
阅读全文