arcpy检查图斑不合理分割面

arcpy检查图斑不合理分割面(即相邻面要素指定的多个字段值不能全部一样)

简单说明

  1. 不合理分割面定义:两个共用一条边的面要素的属性值不能完全一样,如果属性值完全一样应该合并为一个面,如图所示,选中的两个面要素我称之为不合理分割面。
    arcpy检查图斑不合理分割面
  2. 目的:找出指定字段属性值完全一致且相邻的面要素,并导出以便我们查看和定位对应的面要素位置进行修改。

思路

  1. 将待检查的面要素进行融合得到一个融合后的面要素,注意,需要融合为单部件要素。
    arcpy检查图斑不合理分割面

  2. 将待检查的面要素和融合的面要素分别转点要素,这里分别记为输入点和融合点,方便描述。注意,一定要勾选“内部”。
    arcpy检查图斑不合理分割面

  3. 用融合点擦除输入点,得到擦除后的点要素。
    arcpy检查图斑不合理分割面

  4. 用擦除后的点要素去按位置选择待检查的面要素,选中的面要素就是不合理分割面了,导出,完成!

arcpy脚本代码实现步骤

  1. 融合
#in_features 为输入的待检查面要素路径
#dissolve_features 为融合后的面要素路径
#dissolve_feild 要聚合要素的字段
arcpy.Dissolve_management(in_features, dissolve_features, dissolve_feild, "", "SINGLE_PART")
  1. 转点
#in_features_point 输入点要素路径
#dissolve_features_point 融合点要素路径
arcpy.FeatureToPoint_management(in_features, in_features_point, "INSIDE")
arcpy.FeatureToPoint_management(dissolve_features, dissolve_features_point, "INSIDE")
  1. 擦除
#erase_point 擦除后的点要素路径
arcpy.Erase_analysis(in_features_point, dissolve_features_point, erase_point, "0.001 Meters")
  1. 按位置选择并导出
layer_name = "a" #图层名称
out_path= "D:\\other\GISDATA\\test1data\\a.gdb" #导出位置
out_name = "检查结果" #导出要素名称
#添加图层
arcpy.MakeFeatureLayer_management(in_features, layer_name)
#按位置选择
arcpy.SelectLayerByLocation_management(layer_name, "INTERSECT", erase_point, "", "NEW_SELECTION")
#导出要素
arcpy.FeatureClassToFeatureClass_conversion(layer_name, out_path, out_name)
#移除图层
arcpy.Delete_management(layer_name)

待续,有空再附上完整代码。谢谢!

上一篇:2021-05-20


下一篇:20210322_23期_心跳检测_Task03_特征工程