写一个arcgis工具箱脚本:连接要素中多个图斑叠加目标要素其中一个图斑时,取连接要素中多个图斑面积大于多少值的图斑利用空间连接赋值到目标要素图斑里,多个图斑相同字段值赋值到一个目标要素图斑字段时,用"/"隔开并赋值进去,最终导出成果
时间: 2023-12-03 09:44:36 浏览: 82
arcgis中根据excel坐标表批量生产shp图形
3星 · 编辑精心推荐
以下是一个实现上述功能的 Python 脚本:
```python
import arcpy
# 设置工作空间
arcpy.env.workspace = r"C:\path\to\workspace"
# 设置输入和输出要素类
input_fc = "input_features.shp"
target_fc = "target_features.shp"
output_fc = "output_features.shp"
# 设置连接要素中图斑面积的阈值
area_threshold = 1000
# 创建空间连接对象
spatial_join = arcpy.SpatialJoin_analysis(target_fc, input_fc, "in_memory\spatial_join", "JOIN_ONE_TO_MANY")
# 创建空字典用于存储目标要素和连接要素的信息
output_dict = {}
# 遍历连接要素
with arcpy.da.SearchCursor(spatial_join, ["TARGET_FID", "JOIN_FID", "Shape_Area"]) as cursor:
for row in cursor:
# 仅考虑连接要素中面积大于阈值的图斑
if row[2] > area_threshold:
target_fid = row[0]
join_fid = row[1]
# 如果目标要素已经在字典中,则将连接要素的信息添加到该目标要素的值中
if target_fid in output_dict:
output_dict[target_fid].append(join_fid)
# 否则创建新的目标要素和连接要素的键值对
else:
output_dict[target_fid] = [join_fid]
# 创建输出要素类
arcpy.CopyFeatures_management(target_fc, output_fc)
# 遍历输出要素类
with arcpy.da.UpdateCursor(output_fc, ["OID@", "FIELD_NAME"]) as cursor:
for row in cursor:
# 如果该目标要素在字典中,则将连接要素的字段值合并起来并赋值到该目标要素的字段中
if row[0] in output_dict:
field_value = "/".join(str(f) for f in output_dict[row[0]])
row[1] = field_value
cursor.updateRow(row)
```
该脚本首先进行空间连接,然后遍历连接要素,仅考虑面积大于阈值的图斑,并将其添加到一个字典中。然后遍历输出要素类,如果该目标要素在字典中,则将连接要素的字段值合并起来并赋值到该目标要素的字段中。最终将输出要素类导出为结果。请注意,该脚本仅供参考,具体实现可能需要根据实际情况进行修改。
阅读全文