Arcpy基础入门-5、读取自定义格式xml

2024-01-27 05:18

本文主要是介绍Arcpy基础入门-5、读取自定义格式xml,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

arcpy作为一种脚本开发语言,其功能的目的是为用户提供强有力的可扩展性。

这一节通过示例,介绍比较下使用python和arcpy读取自定义格式的xml并构建要素的方法。

先简单描述下自定义的xml格式:

j


自定义了一种北京54坐标系的航迹数据,每一个航迹包括了属性和对应的几何点,由于通用的航迹格式GPS只支持经纬度坐标的读取,所以需要自己写脚本,对要素进行处理


将xml格式的数据转化为要素有两种方式:一种是较为常见的创建要素,赋属性的方式,

使用的是arcpy.insertcursor,这种方式和ArcEngine开发有点相似,代码如下:

#coding=utf-8
import arcpy, os, time
import xml.dom.minidom
import random
from arcpy import env
arcpy.env.overwriteOutput = True #inGPXFile = arcpy.GetParameterAsText(0)
inpath = 'F:\\gpx\\航迹-settingDB-20150605031316.gpx'
uipath = unicode(inpath , "utf8")#outputFC = arcpy.GetParameterAsText(1)
outputFC = r'F:\gpx\20150605031316.shp'
uoutputFC = unicode(outputFC , "utf8")# 打开xml文档
dom = xml.dom.minidom.parse(uipath)
#得到文档元素对象
root = dom.documentElementdef get_attrvalue(node, attrname):return node.getAttribute(attrname) if node else ''def get_nodevalue(node, index = 0):if node:if node.childNodes:return node.childNodes[index].nodeValue   else:return Nonedef get_xmlnode(node,name):return node.getElementsByTagName(name) if node else []trk_nodes = get_xmlnode(root,'trk')
#获得节点测试
print trk_nodes#创建shp
path = os.path.split(outputFC)[0]
name = os.path.split(outputFC)[1]
print path
print name
arcpy.CreateFeatureclass_management(path,name,"POLYLINE")fieldName1 = "name"
fieldAlias = "轨迹名称"
fieldLength = 50
arcpy.AddField_management(outputFC, fieldName1, "TEXT", "","",fieldLength,fieldAlias, "NULLABLE")fieldName2 = "time"
fieldAlias = "时间"
fieldLength = 50
arcpy.AddField_management(outputFC, fieldName2, "TEXT", "","",fieldLength,fieldAlias, "NULLABLE")fieldName3 = "length"
fieldPrecision = 9
fieldAlias = "长度"
arcpy.AddField_management(outputFC, fieldName3, "DOUBLE", fieldPrecision, "", "",fieldAlias, "NULLABLE")fieldName4 = "jiange"
fieldAlias = "间隔"
arcpy.AddField_management(outputFC, fieldName4, "LONG", "","","",fieldAlias, "NULLABLE")fieldName5 = "tpkname"
fieldAlias = "图幅名字"
fieldLength = 50
arcpy.AddField_management(outputFC, fieldName5, "TEXT", "","",fieldLength,fieldAlias, "NULLABLE")fieldName6 = "daino"
fieldAlias = "带号"
arcpy.AddField_management(outputFC, fieldName6, "LONG", "","","",fieldAlias, "NULLABLE")cur=arcpy.InsertCursor(outputFC)
for trk_node in trk_nodes:trk_node_time = get_xmlnode(trk_node,'time')trk_node_length = get_xmlnode(trk_node,'length')trk_node_jiange = get_xmlnode(trk_node,'jiange')trk_node_name = get_xmlnode(trk_node,'name')trk_node_daino = get_xmlnode(trk_node,'daino')trk_node_tpkname = get_xmlnode(trk_node,'tpkname')trk_time =get_nodevalue(trk_node_time[0]).encode('utf-8','ignore')trk_length =get_nodevalue(trk_node_length[0]).encode('utf-8','ignore')trk_jiange =get_nodevalue(trk_node_jiange[0]).encode('utf-8','ignore')trk_name =get_nodevalue(trk_node_name[0]).encode('utf-8','ignore')trk_daino =get_nodevalue(trk_node_daino[0]).encode('utf-8','ignore')trk_tpkname =get_nodevalue(trk_node_tpkname[0]).encode('utf-8','ignore')array = arcpy.Array()trkseg_node = get_xmlnode(trk_node,'trkseg')print trkseg_nodetrkpt_nodes = get_xmlnode(trkseg_node[0],'trkpt')print trkpt_nodesfor trkpt_node in trkpt_nodes:trkpt_node_x = get_xmlnode(trkpt_node,'X')trkpt_node_y = get_xmlnode(trkpt_node,'Y')trkpt_x =get_nodevalue(trkpt_node_x[0]).encode('utf-8','ignore')trkpt_y =get_nodevalue(trkpt_node_y[0]).encode('utf-8','ignore')x=float(trkpt_x)y=float(trkpt_y)pnt=arcpy.Point()pnt.X = trkpt_ypnt.Y = trkpt_xarray.add(pnt)print xprint ypolyline = arcpy.Polyline(array)array.removeAll()feat= cur.newRow()feat.shape = polylinefeat.setValue("shape",polyline)feat.setValue("time",trk_time)feat.setValue("length",float(trk_length))feat.setValue("jiange",int(trk_jiange))feat.setValue("name",trk_name)feat.setValue("daino",int(trk_daino))feat.setValue("tpkname",trk_tpkname)cur.insertRow(feat)del  feat
del cur
还有另外一种方式使用arcpy.da模块下面的插入游标,根据arcgis的帮助文档,这种方式速度会更快, 详情见这里,代码如下:

#coding=utf-8
import arcpy, os, time
import xml.dom.minidom
import random
from arcpy import env
arcpy.env.overwriteOutput = True inGPXFile = arcpy.GetParameterAsText(0)
#inpath = 'F:\\gpx\\航迹-settingDB-20150605031316.gpx'
#uipath = unicode(inpath , "utf8")outputFC = arcpy.GetParameterAsText(1)
#outputFC = r'F:\gpx\20150605031316.shp'
#uoutputFC = unicode(outputFC , "utf8")# 打开xml文档
dom = xml.dom.minidom.parse(inGPXFile)
#得到文档元素对象
root = dom.documentElementdef get_attrvalue(node, attrname):return node.getAttribute(attrname) if node else ''def get_nodevalue(node, index = 0):if node:if node.childNodes:return node.childNodes[index].nodeValue   else:return Nonedef get_xmlnode(node,name):return node.getElementsByTagName(name) if node else []trk_nodes = get_xmlnode(root,'trk')
#获得节点测试
print trk_nodes#创建shp
path = os.path.split(outputFC)[0]
name = os.path.split(outputFC)[1]
print path
print name
arcpy.CreateFeatureclass_management(path,name,"POLYLINE")fieldName1 = "name"
fieldAlias = "轨迹名称"
fieldLength = 50
arcpy.AddField_management(outputFC, fieldName1, "TEXT", "","",fieldLength,fieldAlias, "NULLABLE")fieldName2 = "time"
fieldAlias = "时间"
fieldLength = 50
arcpy.AddField_management(outputFC, fieldName2, "TEXT", "","",fieldLength,fieldAlias, "NULLABLE")fieldName3 = "length"
fieldPrecision = 9
fieldAlias = "长度"
arcpy.AddField_management(outputFC, fieldName3, "DOUBLE", fieldPrecision, "", "",fieldAlias, "NULLABLE")fieldName4 = "jiange"
fieldAlias = "间隔"
arcpy.AddField_management(outputFC, fieldName4, "LONG", "","","",fieldAlias, "NULLABLE")fieldName5 = "tpkname"
fieldAlias = "图幅名字"
fieldLength = 50
arcpy.AddField_management(outputFC, fieldName5, "TEXT", "","",fieldLength,fieldAlias, "NULLABLE")fieldName6 = "daino"
fieldAlias = "带号"
arcpy.AddField_management(outputFC, fieldName6, "LONG", "","","",fieldAlias, "NULLABLE")
cur = arcpy.da.InsertCursor(outputFC, ("time","length","jiange","name","daino","tpkname","SHAPE@"))
for trk_node in trk_nodes:trk_node_time = get_xmlnode(trk_node,'time')trk_node_length = get_xmlnode(trk_node,'length')trk_node_jiange = get_xmlnode(trk_node,'jiange')trk_node_name = get_xmlnode(trk_node,'name')trk_node_daino = get_xmlnode(trk_node,'daino')trk_node_tpkname = get_xmlnode(trk_node,'tpkname')trk_time =get_nodevalue(trk_node_time[0]).encode('utf-8','ignore')trk_length =get_nodevalue(trk_node_length[0]).encode('utf-8','ignore')trk_jiange =get_nodevalue(trk_node_jiange[0]).encode('utf-8','ignore')trk_name =get_nodevalue(trk_node_name[0]).encode('utf-8','ignore')trk_daino =get_nodevalue(trk_node_daino[0]).encode('utf-8','ignore')trk_tpkname =get_nodevalue(trk_node_tpkname[0]).encode('utf-8','ignore')array = arcpy.Array()trkseg_node = get_xmlnode(trk_node,'trkseg')print trkseg_nodetrkpt_nodes = get_xmlnode(trkseg_node[0],'trkpt')print trkpt_nodesfor trkpt_node in trkpt_nodes:trkpt_node_x = get_xmlnode(trkpt_node,'X')trkpt_node_y = get_xmlnode(trkpt_node,'Y')trkpt_x =get_nodevalue(trkpt_node_x[0]).encode('utf-8','ignore')trkpt_y =get_nodevalue(trkpt_node_y[0]).encode('utf-8','ignore')x=float(trkpt_x)y=float(trkpt_y)pnt=arcpy.Point()pnt.X = trkpt_ypnt.Y = trkpt_xarray.add(pnt)print xprint ypolyline = arcpy.Polyline(array)row =[trk_time,float(trk_length),int(trk_jiange),trk_name,int(trk_daino),trk_tpkname,polyline] cur.insertRow(row)
至于如何根据脚本创建模型工具,可以去看下基础入门之3和4



这篇关于Arcpy基础入门-5、读取自定义格式xml的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/649161

相关文章

Python FastAPI入门安装使用

《PythonFastAPI入门安装使用》FastAPI是一个现代、快速的PythonWeb框架,用于构建API,它基于Python3.6+的类型提示特性,使得代码更加简洁且易于绶护,这篇文章主要介... 目录第一节:FastAPI入门一、FastAPI框架介绍什么是ASGI服务(WSGI)二、FastAP

Python中Windows和macOS文件路径格式不一致的解决方法

《Python中Windows和macOS文件路径格式不一致的解决方法》在Python中,Windows和macOS的文件路径字符串格式不一致主要体现在路径分隔符上,这种差异可能导致跨平台代码在处理文... 目录方法 1:使用 os.path 模块方法 2:使用 pathlib 模块(推荐)方法 3:统一使

Maven pom.xml文件中build,plugin标签的使用小结

《Mavenpom.xml文件中build,plugin标签的使用小结》本文主要介绍了Mavenpom.xml文件中build,plugin标签的使用小结,文中通过示例代码介绍的非常详细,对大家的学... 目录<build> 标签Plugins插件<build> 标签<build> 标签是 pom.XML

Java中使用注解校验手机号格式的详细指南

《Java中使用注解校验手机号格式的详细指南》在现代的Web应用开发中,数据校验是一个非常重要的环节,本文将详细介绍如何在Java中使用注解对手机号格式进行校验,感兴趣的小伙伴可以了解下... 目录1. 引言2. 数据校验的重要性3. Java中的数据校验框架4. 使用注解校验手机号格式4.1 @NotBl

SpringBoot自定义注解如何解决公共字段填充问题

《SpringBoot自定义注解如何解决公共字段填充问题》本文介绍了在系统开发中,如何使用AOP切面编程实现公共字段自动填充的功能,从而简化代码,通过自定义注解和切面类,可以统一处理创建时间和修改时间... 目录1.1 问题分析1.2 实现思路1.3 代码开发1.3.1 步骤一1.3.2 步骤二1.3.3

Python批量调整Word文档中的字体、段落间距及格式

《Python批量调整Word文档中的字体、段落间距及格式》这篇文章主要为大家详细介绍了如何使用Python的docx库来批量处理Word文档,包括设置首行缩进、字体、字号、行间距、段落对齐方式等,需... 目录关键代码一级标题设置  正文设置完整代码运行结果最近关于批处理格式的问题我查了很多资料,但是都没

基于Python开发PDF转Doc格式小程序

《基于Python开发PDF转Doc格式小程序》这篇文章主要为大家详细介绍了如何基于Python开发PDF转Doc格式小程序,文中的示例代码讲解详细,感兴趣的小伙伴可以跟随小编一起学习一下... 用python实现PDF转Doc格式小程序以下是一个使用Python实现PDF转DOC格式的GUI程序,采用T

dubbo3 filter(过滤器)如何自定义过滤器

《dubbo3filter(过滤器)如何自定义过滤器》dubbo3filter(过滤器)类似于javaweb中的filter和springmvc中的intercaptor,用于在请求发送前或到达前进... 目录dubbo3 filter(过滤器)简介dubbo 过滤器运行时机自定义 filter第一种 @A

Python如何实现读取csv文件时忽略文件的编码格式

《Python如何实现读取csv文件时忽略文件的编码格式》我们再日常读取csv文件的时候经常会发现csv文件的格式有多种,所以这篇文章为大家介绍了Python如何实现读取csv文件时忽略文件的编码格式... 目录1、背景介绍2、库的安装3、核心代码4、完整代码1、背景介绍我们再日常读取csv文件的时候经常

C#中读取XML文件的四种常用方法

《C#中读取XML文件的四种常用方法》Xml是Internet环境中跨平台的,依赖于内容的技术,是当前处理结构化文档信息的有力工具,下面我们就来看看C#中读取XML文件的方法都有哪些吧... 目录XML简介格式C#读取XML文件方法使用XmlDocument使用XmlTextReader/XmlTextWr