python arcpy练习之面要素重叠拓扑检查
更新时间:2021年9月18日 08:06 点击:1509
需求
有多个文件地理数据库(gdb),数据库内有多个面要素类图层,每个图层不能有自重叠,也不能和其他图层重叠。所以,需要为每个文件地理数据库(gdb)进行拓扑检查。
思路
1.在输出文件夹下新建与输入文件夹下同名的gdb。
2.在gdb内新建要素类数据集。
3.在数据集内创建拓扑。
4.向拓扑中添加要素类。
5.添加拓扑规则。
6.拓扑验证。
7.导出拓扑错误。
将需要拓扑的gdb放入同一个文件夹下,用python代码串连起来利用循环遍历gdb便可以进行批量拓扑了。
代码
# -*- coding: gbk -*- import os import arcpy import sys reload(sys) sys.setdefaultencoding('utf-8') arcpy.AddMessage(sys.getdefaultencoding()) arcpy.env.XYResolution = "0.0001 Meters" #设置XY分辨率 arcpy.env.XYTolerance = "0.001 Meters" #设置XY容差 def checkTopology(in_path, out_path): cout = 0 fail = 0 warning = 0 fail_list = [] warning_list = [] arcpy.env.workspace = in_path workspaces = arcpy.ListWorkspaces("*", "ALL") for workspace in workspaces: arcpy.AddMessage("=" * 60) cout += 1 arcpy.AddMessage(" (" + str(cout) + ") " + os.path.basename(workspace)) arcpy.AddMessage(" ") try: # 将每个gdb设为工作区 arcpy.env.workspace = workspace fc_list = arcpy.ListFeatureClasses() #列出gdb内的要素类图层 in_fc_path_list = [] for fc in fc_list: in_fc_path_list.append(os.path.join(workspace,fc)) # 在输出路径out_path下新建gdb gdb_name1 = os.path.basename(workspace) result_gdb = os.path.join(out_path, gdb_name1) if not arcpy.Exists(result_gdb): arcpy.AddMessage(" 新建 " + result_gdb + " ...") arcpy.CreateFileGDB_management(out_path, gdb_name1) else: arcpy.AddWarning(result_gdb + " 已存在!") warning_list.append(result_gdb) warning += 1 continue arcpy.AddMessage(" 在 " + result_gdb + " 内新建要素数据集...") dataset_name = "dataset" dataset_path = os.path.join(result_gdb,dataset_name) # 引用包含要应用的空间参考的要素类或要素数据集, # 将in_fc_path_list[0]的坐标作为要素数据集的坐标,所以gdb内的所有要素类应该是统一坐标 arcpy.CreateFeatureDataset_management(result_gdb, dataset_name, in_fc_path_list[0]) #导入要素类到数据集 arcpy.AddMessage(" 向 " + dataset_path + " 导入要素类图层...") arcpy.FeatureClassToGeodatabase_conversion(in_fc_path_list,dataset_path) #在dataset数据内创建拓扑 arcpy.AddMessage(" 创建拓扑...") topology_name = "topology" topology_path = os.path.join(dataset_path, topology_name) arcpy.CreateTopology_management(dataset_path, topology_name) # 将导入dataset的要素类添加到拓扑中 arcpy.AddMessage(" 向拓扑中添加要素类...") dataset_fc_path_lsit = [] #数据集中的要素类绝对路径列表 for ifc_name in fc_list: # 拼接数据集中的要素类绝对路径 dataset_fc_path = os.path.join(dataset_path,ifc_name) dataset_fc_path_lsit.append(dataset_fc_path) arcpy.AddFeatureClassToTopology_management(topology_path, dataset_fc_path, "1", "1") arcpy.AddMessage(" 添加拓扑规则...") for i in range(len(dataset_fc_path_lsit)): fc_path1 = dataset_fc_path_lsit[i] #规则1:不能重叠 arcpy.AddRuleToTopology_management(topology_path, "Must Not Overlap (Area)", fc_path1, "", "", "") for j in range(i + 1, len(dataset_fc_path_lsit)): fc_path2 = dataset_fc_path_lsit[j] # 规则2:不能重叠与其他要素类重叠 arcpy.AddRuleToTopology_management(topology_path, "Must Not Overlap With (Area-Area)", fc_path1, "", fc_path2, "") arcpy.AddMessage(" 拓扑验证...") arcpy.ValidateTopology_management(topology_path, "Full_Extent") arcpy.AddMessage(" 导出拓扑错误...") arcpy.ExportTopologyErrors_management(topology_path, dataset_path, "topoError") except Exception as e: arcpy.AddError(e.message) arcpy.AddError(" " + os.path.basename(workspace) + " 失败!") fail_list.append(os.path.basename(workspace)) fail += 1 arcpy.AddMessage('+' * 60) arcpy.AddMessage(u" 成功:" + str(cout - fail - warning) + u" 个!") if warning > 0: arcpy.AddWarning(u" 警告:" + str(warning) + u" 个! 如下:") arcpy.AddWarning("####" + '*' * 20) for fff in warning_list: arcpy.AddWarning(" " + fff) arcpy.AddWarning("####" + '*' * 20) if fail > 0: arcpy.AddError(u" 失败:" + str(fail) + u" 个! 如下:") arcpy.AddError("####" + '*' * 20) for ff in fail_list: arcpy.AddError(" " + ff) arcpy.AddError("####" + '*' * 20) arcpy.AddMessage('+' * 60) if __name__ == "__main__": in_path_globle = arcpy.GetParameterAsText(0) out_path_globe = arcpy.GetParameterAsText(1) checkTopology(in_path_globle,out_path_globe)
效果图
总结
本篇文章就到这里了,希望能够给你带来帮助,也希望您能够多多关注猪先飞的更多内容!
下一篇: Python logging简介详解
相关文章
- 这篇文章主要介绍了python-opencv-画外接矩形框的实例代码,代码简单易懂,对大家的学习或工作具有一定的参考借鉴价值,需要的朋友可以参考下...2021-09-04
Python astype(np.float)函数使用方法解析
这篇文章主要介绍了Python astype(np.float)函数使用方法解析,文中通过示例代码介绍的非常详细,对大家的学习或者工作具有一定的参考学习价值,需要的朋友可以参考下...2020-06-08- 2022虎年新年即将来临,小编为大家带来了一个利用Python编写的虎年烟花特效,堪称全网最绚烂,文中的示例代码简洁易懂,感兴趣的同学可以动手试一试...2022-02-14
- 在本篇文章里小编给大家分享的是一篇关于python中numpy.empty()函数实例讲解内容,对此有兴趣的朋友们可以学习下。...2021-02-06
python-for x in range的用法(注意要点、细节)
这篇文章主要介绍了python-for x in range的用法,具有很好的参考价值,希望对大家有所帮助。一起跟随小编过来看看吧...2021-05-10- 这篇文章主要介绍了Python 图片转数组,二进制互转操作,具有很好的参考价值,希望对大家有所帮助。一起跟随小编过来看看吧...2021-03-09
- 这篇文章主要介绍了Python中的imread()函数用法说明,具有很好的参考价值,希望对大家有所帮助。一起跟随小编过来看看吧...2021-03-16
- 这篇文章主要介绍了python如何实现b站直播自动发送弹幕,帮助大家更好的理解和学习使用python,感兴趣的朋友可以了解下...2021-02-20
python Matplotlib基础--如何添加文本和标注
这篇文章主要介绍了python Matplotlib基础--如何添加文本和标注,帮助大家更好的利用Matplotlib绘制图表,感兴趣的朋友可以了解下...2021-01-26- 这篇文章主要介绍了解决python 使用openpyxl读写大文件的坑,具有很好的参考价值,希望对大家有所帮助。一起跟随小编过来看看吧...2021-03-13
- 今天小编就为大家分享一篇python 计算方位角实例(根据两点的坐标计算),具有很好的参考价值,希望对大家有所帮助。一起跟随小编过来看看吧...2020-04-27
- 这篇文章主要为大家详细介绍了python实现双色球随机选号,文中示例代码介绍的非常详细,具有一定的参考价值,感兴趣的小伙伴们可以参考一下...2020-05-02
- 在本篇文章里小编给大家整理的是一篇关于python中使用np.delete()的实例方法,对此有兴趣的朋友们可以学习参考下。...2021-02-01
- 这篇文章主要介绍了使用Python的pencolor函数实现渐变色功能,本文通过实例代码给大家介绍的非常详细,对大家的学习或工作具有一定的参考借鉴价值,需要的朋友可以参考下...2021-03-09
Python getsizeof()和getsize()区分详解
这篇文章主要介绍了Python getsizeof()和getsize()区分详解,文中通过示例代码介绍的非常详细,对大家的学习或者工作具有一定的参考学习价值,需要的朋友们下面随着小编来一起学习学习吧...2020-11-20- 这篇文章主要介绍了python自动化办公操作PPT的实现,文中通过示例代码介绍的非常详细,对大家的学习或者工作具有一定的参考学习价值,需要的朋友们下面随着小编来一起学习学习吧...2021-02-05
- 这篇文章主要介绍了解决python 两个时间戳相减出现结果错误的问题,具有很好的参考价值,希望对大家有所帮助。一起跟随小编过来看看吧...2021-03-12
- 这篇文章主要为大家详细介绍了python实现学生通讯录管理系统,文中示例代码介绍的非常详细,具有一定的参考价值,感兴趣的小伙伴们可以参考一下...2021-02-25
- 这篇文章主要介绍了PyTorch一小时掌握之迁移学习篇,本文给大家介绍的非常详细,对大家的学习或工作具有一定的参考借鉴价值,需要的朋友可以参考下...2021-09-08
- 这篇文章主要介绍了Python绘制的爱心树与表白代码,本文通过实例代码给大家介绍的非常详细,对大家的学习或工作具有一定的参考借鉴价值,需要的朋友可以参考下...2021-04-06