python arcpy练习之面要素重叠拓扑检查
作者:Ah·Week 发布时间:2021-12-24 12:41:42
标签:python,arcpy,练习
需求
有多个文件地理数据库(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)
效果图
来源:https://blog.csdn.net/qz_CSDN/article/details/120330887


猜你喜欢
- 一:需重定义神经网络继续训练的方法1.训练代码import numpy as npimport tensorflow as tfx_data
- 在我们学习python的过程中,学习序列是一门必修课。当我们掌握了序列过后,便会学习常用的两个排序函数sort()与sorted()。但很少
- 1.什么是多线程?多线程是为了同步完成多项任务,不是为了提高运行效率,而是为了提高资源使用效率来提高系统的效率。线程是在同一时间需要完成多项
- 本文实例讲述了Python Django框架模板渲染功能。分享给大家供大家参考,具体如下:项目名/settings.py(项目配置,配置模板
- MySQL Binary Log也就是常说的bin-log, ,是mysql执行改动产生的二进制日志文件,其主要作用有两个: * 数据回复
- 前言order by满足两种情况,会使用 index 方式排序:order by语句使用索引最左前列(最左匹配法则)where子句和orde
- 想用C++写项目,数据库是必须的,所以这两天学了一下C++操作MySQL数据库的方法。也没有什么教程,就是在网上搜的知识,下面汇总一下。 连
- 系统环境:Win10 64位MySQL版本:mysql-5.7.18-winX64部署的步骤就是按照网上说的:1:修改环境变量path,增加
- 1、什么是全局变量?在Python中,全局变量指的是可以作用于函数内部和外部的变量。在这里有两种情况:在函数的外部定义和内部定义添加glob
- Iterable – 可迭代对象能够逐一返回其成员项的对象。 可迭代对象的例子包括所有序列类型 (例如 list, str 和 tuple)
- 新闻、文章系统中经常会用到的一个功能,添加上一篇,下一篇或相关文章的功能可以增加访客停留的时间。也许新人在做这上一篇、下一篇功能时使用的是I
- 最新在学习Python的基础入门系列课程,今天学习到使用python 的内置库smtplib发送邮件内容。使用Python发送邮件步骤简单:
- 本实验中分别针对空库、脱机、联机三种方式,配置一主两从的mysql标准异步复制。只做整服务器级别的复制,不考虑对个别库表或使用过滤复制的情况
- 在日常的python编程中使用这几个函数来简化我们的编程工作,经常使用能使编程效率大大地提高。1. Map 函数map函数可以使用另外一个函
- javascript/js的ajax的GET请求代码如下所示:<script type="text/javascript&q
- pytorch自定义不可导激活函数今天自定义不可导函数的时候遇到了一个大坑。首先我需要自定义一个函数:sign_fimport torchf
- 先给大家看下Scratch3.0二次开发之windows环境下打包成exe的流程。1、需要先安装npm,安装过程不作过多介绍了。
- 搞了好几天的表格字体格式,一直想找一种能直接一次性修改表格所有字体格式的方法(函数),但是无论用什么方法都无法修改表格字体的格式,原因应该是
- 今天使用os.path.isdir()判断是否是文件夹的时候发现一个问题:lst = os.listdir(path) &nb
- 文件操作的一般内容:# 文件的操作# 打开文件 open 打开已存在文件 或者创建一个新文件open('./Test.txt'