ArcPy快速检查相邻面要素多个属性取值是否一致

需求

  • 规则名称:地类图斑属性字段重复性检查
  • 检查内容:检查相邻地类图斑属性字段内容是否重复
  • 检查要求:地类图斑层空间位置为相邻,且地类编码相同的图斑,权属单位代码、坐落单位代码、耕地类型、耕地坡度级别、图斑细化代码、耕地种植属性代码、耕地等别等内容不能全部相同。

难点

对一个区县来说,检查的图斑个数可达数万甚至十万之多。如果使用ArcGIS Engine,逐个图斑去遍历、查询、对比,性能势必成问题。需要引入ArcPy,使用空间分析的手段来解决问题。

思路

  1. 既然要求相邻图斑多个属性字段取值不能完全一致,那么,我们将所有属性字段取值完全一致的相邻图斑进行融合,便能得到要检查不通过的图斑。
  2. 对不通过检查的融合后的图斑,需要找到其对应原图斑的位置,进行报错。

步骤

  1. 将多个属性完全相同的相邻图斑融合到一个新的图层。注意,需融合为单部件要素。原图层与融合后的图层分别记为source_name、dissolve_name。
arcpy.Dissolve_management(source_name, dissolve_name, dissolve_fields, None, "SINGLE_PART", "DISSOLVE_LINES")
  1. 将融合前后两个图层分别转为多段线,做几何差集处理,得到相异的部分,即为被融合的线段,记为sym_diff_name图层。
arcpy.PolygonToLine_management(source_name, source_line_name)
arcpy.PolygonToLine_management(dissolve_name, dissolve_line_name)
arcpy.SymDiff_analysis(dissolve_line_name, source_line_name, sym_diff_name, "ONLY_FID")
  1. 将dissolve_name构建为图层,使用sym_diff_name选择被包含的要素并导出为dissolve_selection_name图层。
arcpy.MakeFeatureLayer_management(dissolve_name, dissolve_layer_name)
arcpy.SelectLayerByLocation_management(dissolve_layer_name, "CONTAINS", sym_diff_name, None, "NEW_SELECTION")
arcpy.FeatureClassToFeatureClass_conversion(dissolve_layer_name, target_gdb_path, dissolve_selection_name)
  1. 将source_name构建为图层,使用dissolve_selection_name选择被包含的要素并导出为source_selection_name图层,则所有相邻且属性一致的要素均在该图层中,且可使用dissolve_selection_name逐个要素包含查询查到。
arcpy.MakeFeatureLayer_management(source_name, source_layer_name)
arcpy.SelectLayerByLocation_management(source_layer_name, "WITHIN", dissolve_selection_name, None,
                                       "NEW_SELECTION")
arcpy.FeatureClassToFeatureClass_conversion(source_layer_name, target_gdb_path, source_selection_name)
  • 3
    点赞
  • 12
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值