arcpy检查图斑不合理分割面

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

简单说明

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

思路

1. 将待检查的面要素进行融合得到一个融合后的面要素，注意，需要融合为单部件要素。
   

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

3. 用融合点擦除输入点，得到擦除后的点要素。
   

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

arcpy脚本代码实现步骤

1. 融合

#in_features 为输入的待检查面要素路径
#dissolve_features 为融合后的面要素路径
#dissolve_feild 要聚合要素的字段
arcpy.Dissolve_management(in_features, dissolve_features, dissolve_feild, "", "SINGLE_PART")

2. 转点

#in_features_point 输入点要素路径
#dissolve_features_point 融合点要素路径
arcpy.FeatureToPoint_management(in_features, in_features_point, "INSIDE")
arcpy.FeatureToPoint_management(dissolve_features, dissolve_features_point, "INSIDE")

3. 擦除

#erase_point 擦除后的点要素路径
arcpy.Erase_analysis(in_features_point, dissolve_features_point, erase_point, "0.001 Meters")

4. 按位置选择并导出

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)