在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