arcpy批量处理矢量面要素,合并面中被包围的小面

您好!根据您的要求,您可以使用arcpy模块中的适当工具来批量处理矢量面要素。以下是可能实现此目标的一些步骤:

1. 导入必要的模块和设置工作环境


import arcpy
from arcpy import env
env.workspace = r"path/to/workspace"
 

2. 循环遍历每个面要素,并使用SelectLayerByLocation工具选择被包含在其中的面要素


for fc in arcpy.ListFeatureClasses(feature_type="Polygon"):
    lyr = arcpy.MakeFeatureLayer_management(fc, "lyr")
    arcpy.SelectLayerByLocation_management(lyr, "COMPLETELY_WITHIN", fc)
 

3. 对于选择的面要素,计算其面积与包含要素面积的比例,并对小于1%的要素进行删除,并将其属性写入到包含要素的属性中


    for sub_fc in arcpy.ListFeatureClasses(feature_type="Polygon"):
        sub_lyr = arcpy.MakeFeatureLayer_management(sub_fc, "sub_lyr")
        arcpy.SelectLayerByLocation_management(sub_lyr, "COMPLETELY_WITHIN", lyr)
        if int(arcpy.GetCount_management(sub_lyr).getOutput(0)) > 0:
            for row in arcpy.da.SearchCursor(sub_lyr, ["SHAPE@AREA", "*"]):
                area_ratio = float(row[0]) / float(arcpy.Describe(lyr).shapeArea)
                if area_ratio < 0.01:
                    arcpy.DeleteRows_management(sub_lyr)
                    with arcpy.da.UpdateCursor(lyr, ["*"]) as upd_cur:
                        for upd_row in upd_cur:
                            if upd_row[0] == row[1]:
                                upd_row[-1] = "该面要素被删掉合并"
                                upd_cur.updateRow(upd_row)
 

4. 清理选择的要素图层,以便进行下一次循环


    arcpy.ClearWorkspaceCache_management()
 

希望这可以帮助您实现您的目标!

猜你喜欢

转载自blog.csdn.net/weixin_58851039/article/details/131161399