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
0
投稿
猜你喜欢
- array_unique() 定义和用法 array_unique() 函数移除数组中的重复的值,并返回结果数组。 当几个数组元素的值相等时
- 循环是我们经常用到的一个概念,比如,循环计算数字叠加、循环输出文字内容等。循环是运行重复内容的一个最简单的方法,简化了代码流程,增加了时效性
- 本文主要介绍Python3.6及TensorFlow的安装和配置流程。一、Python官网下载自己电脑和系统对应的Python安装包。网址:
- 机器学习可应用在各个方面,本篇将在系统性进入机器学习方向前,初步认识机器学习,利用线性回归预测波士顿房价;原理简介利用线性回归最简单的形式预
- 一、常见模型分类1.1、循环服务器模型循环接收客户端请求,处理请求。同一时刻只能处理一个请求,处理完毕后再处理下一个。优点:实现简单,占用资
- 本例详细介绍了如何在wiondws XP下安装与配置MySQL5.0.37 ,图文并茂,相信对初学mysql的朋友有所帮助。1 点击MySQ
- 对文本类文件(*.txt;*.html;*.doc;等等),图片类文件(*.jpg;*.gif等等)直接点击链接时会在浏览器打开,而无法出现
- 很多人喜欢把一个网站中相同的部分象是统一的页面logo,版权声明等做成一个过程,然后放到一个include文件中,这样所有的页面就都可以使用
- 新年钟声刚过,淘宝新版首页全“心”上线了,这次设计大胆的将布局从 960px 伸展至 1000px,页面更通透,新首页更大范围的实践了 HT
- 安装pyinstallerpip install pyinstaller制作项目的.spec文件 进入django项目所在路径,
- 目录一、前言二、基础命令三、正常显示中文:四、设置图样或子图五、设置x轴或y轴相关属性:六、设置标题:七、设置图例:八、进行标注:九、保存图
- python实现超市扫码仪计费的程序主要是使用超市扫码仪扫商品的条形码,读取商品信息,实现计费功能。主要用到的技术是串口通信,数据库的操作,
- 用过软件的朋友都知道,进度条是一个优秀软件的重要组成部分。它的存在能够使用户及时掌握程序的运行进度,确认应用程序正常工作。可是ASP中似乎没
- 本文实例讲述了Golang排列组合算法问题之全排列实现方法。分享给大家供大家参考,具体如下:【排列组合问题】一共N辆火车(0<N<
- 阅读上一篇教程:WEB2.0网页制作标准教程(9)第一个CSS布局实例如果我们想在3列布局的最后加一行页脚,放版权之类的信息。就遇到必须对齐
- ul设置浮动后不能自适应高度,也就是不能撑开父容器,不能自适应内容的高度。解决方法是在ul结束标签前加个清除浮动。 &
- 先看下面例子的效果:<INPUT TYPE="text" NAME=""&
- 本文实例讲述了php测试kafka项目。分享给大家供大家参考,具体如下:概述Kafka是最初由Linkedin公司开发,是一个分布式、分区的
- 字符串是日常开发中用到最多的数据类型之一,了解字符串首先要明白定界符。定界符有4种,单引号、双引号、heredoc、nowdoc,常使用的单
- 效果展示爬取目标网址:酷我音乐工具使用开发工具:pycharm开发环境:python3.7, Windows10使用工具包:requests