在ArcGIS Pro中,可以使用Python编写代码来捕获与其他多边形重叠或包含的多边形。以下是一个解决方法的代码示例:
import arcpy
# 设置工作空间和输出要素类的路径
workspace = r"C:\data\workspace.gdb"
output_feature_class = "overlapping_polygons"
# 设置输入要素类的路径
input_feature_class = r"C:\data\input_polygons.shp"
# 创建一个空的要素类作为输出
arcpy.CreateFeatureclass_management(workspace, output_feature_class, "POLYGON")
# 定义要素类字段
fields = [
arcpy.FieldDefinitiion("ID", "LONG"),
arcpy.FieldDefinitiion("Overlap_ID", "LONG"),
]
arcpy.AddField_management(output_feature_class, fields[0])
arcpy.AddField_management(output_feature_class, fields[1])
# 使用游标遍历输入要素类中的所有要素
with arcpy.da.SearchCursor(input_feature_class, ["OID@", "SHAPE@"]) as cursor:
for row in cursor:
oid = row[0]
shape = row[1]
# 创建一个空的几何对象
output_shape = arcpy.Geometry()
# 使用游标遍历输入要素类中的所有要素(第二次遍历)
with arcpy.da.SearchCursor(input_feature_class, ["OID@", "SHAPE@"]) as cursor2:
for row2 in cursor2:
oid2 = row2[0]
shape2 = row2[1]
# 跳过自身
if oid == oid2:
continue
# 判断多边形之间的关系
if shape.contains(shape2) or shape2.contains(shape):
# 创建一个新的几何对象,将重叠的部分添加到其中
new_shape = shape.intersection(shape2)
# 将重叠的多边形添加到输出要素类中
with arcpy.da.InsertCursor(output_feature_class, ["ID", "Overlap_ID", "SHAPE@"]) as insert_cursor:
insert_row = (oid, oid2, new_shape)
insert_cursor.insertRow(insert_row)
请注意,以上代码仅提供了一个基本的解决方案框架,并且需要根据实际情况进行修改和适应。您需要根据自己的数据和需求修改工作空间、输入要素类和输出要素类的路径,以及字段名和字段类型。
上一篇:捕获所有以GUID结尾的URL