arcpy检查图斑不合理分割面(即相邻面要素指定的多个字段值不能全部一样)
简单说明
- 不合理分割面定义:两个共用一条边的面要素的属性值不能完全一样,如果属性值完全一样应该合并为一个面,如图所示,选中的两个面要素我称之为不合理分割面。
- 目的:找出指定字段属性值完全一致且相邻的面要素,并导出以便我们查看和定位对应的面要素位置进行修改。
思路
-
将待检查的面要素进行融合得到一个融合后的面要素,注意,需要融合为单部件要素。
-
将待检查的面要素和融合的面要素分别转点要素,这里分别记为输入点和融合点,方便描述。注意,一定要勾选“内部”。
-
用融合点擦除输入点,得到擦除后的点要素。
-
用擦除后的点要素去按位置选择待检查的面要素,选中的面要素就是不合理分割面了,导出,完成!
arcpy脚本代码实现步骤
- 融合
#in_features 为输入的待检查面要素路径
#dissolve_features 为融合后的面要素路径
#dissolve_feild 要聚合要素的字段
arcpy.Dissolve_management(in_features, dissolve_features, dissolve_feild, "", "SINGLE_PART")
- 转点
#in_features_point 输入点要素路径
#dissolve_features_point 融合点要素路径
arcpy.FeatureToPoint_management(in_features, in_features_point, "INSIDE")
arcpy.FeatureToPoint_management(dissolve_features, dissolve_features_point, "INSIDE")
- 擦除
#erase_point 擦除后的点要素路径
arcpy.Erase_analysis(in_features_point, dissolve_features_point, erase_point, "0.001 Meters")
- 按位置选择并导出
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)
待续,有空再附上完整代码。谢谢!