标签：要素 features point 不合理 图斑 arcpy dissolve name

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)

待续，有空再附上完整代码。谢谢！

标签：要素,features,point,不合理,图斑,arcpy,dissolve,name
来源： https://blog.csdn.net/qz_CSDN/article/details/115256082

本站声明： 1. iCode9 技术分享网（下文简称本站）提供的所有内容，仅供技术学习、探讨和分享；
2. 关于本站的所有留言、评论、转载及引用，纯属内容发起人的个人观点，与本站观点和立场无关；
3. 关于本站的所有言论和文字，纯属内容发起人的个人观点，与本站观点和立场无关；
4. 本站文章均是网友提供，不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属；如您发现该文章侵犯了您的权益，可联系我们第一时间进行删除；
5. 本站为非盈利性的个人网站，所有内容不会用来进行牟利，也不会利用任何形式的广告来间接获益，纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。