积分商城

最新搜索： HDMI2.1 环形队列 ZXDN10 定位 ADS1115 TM1723
您现在的位置是：首页 > 技术资料 > Python ArcPy在地质信息处理中的基本应用
推荐星级：
  • 1
  • 2
  • 3
  • 4
  • 5

Python ArcPy在地质信息处理中的基本应用

更新时间：2019-12-26 10:48:50 大小：1M 上传用户：songhuahua查看TA发布的资源 标签：python地质信息处理 下载积分：1分 评价赚积分 （如何评价?） 打赏 收藏 评论(0) 举报

资料介绍

Python语言作为一种开源语言,在GIS软件中,常作为辅助工具使用,其为Arc GIS提供了大量的处理模块,能更高效地完成地质信息处理工作。利用Python语言,调用集成Arc GIS中Arc Py软件包,能实现地质数据的输入、投影变换以及标准图幅查询等功能,同时其代码简单易读,便于使用,可为地质数据处理人员提供更方便的工具。


部分文件列表

文件名 大小
Python_ArcPy在地质信息处理中的基本应用.pdf 1M

【关注视频号领20积分】   【关注公众号立即送20积分】

部分页面预览

（完整内容请下载后查看）
·15·  
2017，33( 1) : 15 ～ 17  
CN 53 － 1124 /TD ISSN 1007 － 9394  
地矿测绘  
Surveying and Mapping of Geology and Mineral Resources  
*
Python ArcPy  
在地质信息处理中的基本应用  
叶宝莹  
(
，
中国地质大学 地质调查研究院 北京  
100083)  
: Python  
，
语言作为一种开源语言 在  
GIS  
，
，
ArcGIS ，  
提供了大量的处理模块 能更高效地完  
摘要  
软件中 常作为辅助工具使用 其为  
ArcGIS ArcPy ， 、  
软件包 能实现地质数据的输入 投影变换以及标准图幅查询等  
。
成地质信息处理工作 利用  
Python  
，
语言 调用集成  
中
， ， ，  
功能 同时其代码简单易读 便于使用 可为地质数据处理人员提供更方便的工具  
。
: Python; ArcGIS; ArcPy;  
关键词  
中图分类号  
地质信息处理  
: P 209; TP 311  
: A  
: 1007 － 9394( 2017) 01 － 0015 － 03  
文章编号  
文献标识码  
Basic Application of Python ArcPy in Geological Information Processing  
YE Bao-ying  
( Institute of Geological Survey，China University of Geosciences，Beijing 100083，China)  
Abstract: Python as an open source language，in the GIS software，often has been used as an auxiliary tool． Especial-  
ly in ArcGIS，it provides a large number of processing modules，more efficient completion of geological information pro-  
cessing． The use of Python language and ArcPy software package integrated in ArcGIS is so easy to realize input，map pro-  
jection transformation and query work of geological data． The code is simple，readable and easy to use，which provides a  
more convenient tool for geological data processing personnel．  
Key words: Python; ArcGIS; ArcPy; geological information processing  
0
1
引言  
数据输入  
［1］  
， 、 、  
在实际地质应用中 主要涉及矿点 矿权 地质路线的生成  
。
Python  
C
是基于 语言的高级动态解释性语言  
。Python  
具
，
野外采集的数据多为坐标对 转入  
GIS  
，
时 需要生成相应的点线  
: 、 ; ，  
有以下特点 简单易学 免费开源 跨平台的解释性高级语言 具  
。
面文件 点文件  
，ArcGIS 10． x desktop  
ad-  
软件中有现成的工具  
Excel，dbf 。Excel  
;
有很强的健壮性的面向对象编程语言 可扩性  
。
dxydata  
。
工具可用 点文件可以为  
等格式  
文件  
ArcPy ArcGIS  
。
处理脚本的软件包 其目的是利用  
是一套  
建立实用而且高效的地理数据分析 转换 数据管理以及  
Python ArcPy  
， ，  
无法设置属性字段的长度 类型 会引起错误  
。dbf  
文件则需要  
Python  
自动制图的处理工具  
，Python  
、
、
Access  
， ，  
中进行字段设置 然后导出才可以使用 过程有些  
在
。
在
中使用  
。Python  
是一种通用的编程语言  
的另一个主要原  
是一种支持动态  
输入的解释型语言 适用于交互式操作以及为称为脚本的一次  
。
繁琐  
因是  
ArcGIS  
———ArcInfo Workstaion  
，
版 本 提 供 了  
早 期 产 品  
，
point、line、polygon  
，generate 。Work-  
命令  
DOS  
等文本文件导入命令  
，
性程序快速制作原型 同时其具有编写大型应用程序的强大功  
station  
，
为数据批量处理的绝佳工具 只是其是基于  
命令  
。
ArcPy ArcGIS  
，
能
利用  
来自多个不同领域的  
小群体开发的附加模块  
在进行地质信息管理与计算分析时 对于各种复杂的地质  
编写的  
应用程序的优势在于 可以使用  
， ， ，  
的 需要使用者熟记数百个命令 并熟悉操作流程 而且过程不  
GIS  
Python  
专业人员和程序员组成的众多  
［2］  
， ， 。  
可逆 若出现操作错误 会引起无法挽回的损失 其是专业数据  
。
，
处理人员才使用的工具 但  
ESRI  
Workstaion 10  
公司在 版本并  
，
。
， ， ，  
数据 进行程序化处理 可以大大减少冗余工作 并保障数据的  
不提供更新版本  
，
基于此 若有大量文本数据需要转入  
ArcGIS  
，
中时 工作比  
。ESRI ArcPy  
公司提供的  
处理质量  
软件包正是针对该类问题  
Python ，  
程序 实现地质数据  
。
应运而生的 本文通过几个简单的  
。ESRI  
ArcPy  
， Python  
模块 可利用 语言可以  
较繁琐  
。
完成这类工作  
公司提供的  
、
的输入 投影变换以及图幅查询工作  
。
*
: 2017 － 01 － 05  
收稿日期  
·16·  
2017  
3
年 月  
地矿测绘  
1． 1  
row =［pnt． num，pnt． name，points］  
cursor． insertRow( row)  
f． close( )  
del cursor  
点文件生成  
首先确定点数据文本文件格式 最好是  
Excel csv shapefile  
，
txt  
，
格式 也可以是  
。
。
格式  
格式 生成新的文件格式为  
点文本文件格式  
点文件可以带有坐标及相应的属性 属性可以添加逗号分  
1． 1． 1  
。
1． 2  
、
线 面文件生成  
， 1  
开 如图 所示  
。
， ， 、  
处理地质路线以及矿区 工作区等信息时 需要生成线 面  
。
等图形文件  
， ，  
线和面文件生成方法几乎相同 文本数据格式也相同 线和  
。
面的区别在于面文件要求首尾坐标重合  
1． 2． 1  
线由一组有序坐标对组成 并且具有属性 为  
1
图
点文件格式  
Fig． 1 Point file format  
Python  
线面文本文件格式  
，
，
txt  
，
格式 具体  
3
如图 所示  
。
，
，
语言处理成该格式 处  
如果不是该格式 可以利用  
。 : ; 、  
理方法简单有效 要求 利用坐标生成点文件 并将序号 名称  
， 2  
作为属性写入到属性字段中 如图 所示  
。
3
图
线文件格式  
Fig． 3 Line file format  
， ( ，  
面由一组有序坐标对组成 且首尾坐标相同 若不同 首尾  
) ， ，  
自动连线形成封闭 并且具有属性 为  
txt  
，
格式 具体如图  
4
。
所示  
4
图
面文件格式  
Fig． 4 Polygon file format  
ESRI shapefile  
，
面文件 并且属性写入属性字  
要求生成矢量  
2
图
点文件生成属性表  
。
段中  
Fig． 2 Generating attribute table from point file  
1． 2． 2  
关键代码  
1． 1． 2  
关键代码  
， ‘polygon’  
生成面文件的程序代码 线文件只需将参数 替换  
#
读取点文件  
‘polyline’。  
为
f = file( " D: \\soft \\yePython \\test \\points． txt" )  
f = file( " D: \\soft \\yePython \\test \\polys． txt" )  
fc = " D: \\soft \\yePython \\test \\newpoly． shp"  
#
shapefile  
确定新  
文件的名称及位置  
fc = " D: \\soft \\yePython \\test \\newpoint． shp"  
# define projection ，3 methods ，code ，name or * ． prj file  
spatiref = ArcPy． SpatialReference( " Beijing 1954" )  
# create a new shapefiles  
#
， ， 、  
定义投影 有三种方法可以选择 投影代码 投影名称以  
． prj  
及投影文件名称  
# spatiref = ArcPy． SpatialReference( 4214)  
ArcPy． CreateFeatureclass_management( " D: \\ soft \\ yePython  
\\test" ，" newpoly． shp" ，" Polygon" ，" " ，" " ，" " ，spatiref)  
# create fields，num，name  
# attention，spatial reference name MUST be ＇show in arcMap＇．  
spatiref = ArcPy． SpatialReference( " Beijing 1954" )  
#
shapefile  
点文件  
创建新的  
ArcPy． CreateFeatureclass_management( " D: \\ soft \\ yePython  
\\test" ，" newpoint． shp" ，" POINT" ，" " ，" " ，" " ，spatiref)  
ArcPy． AddField_management ( fc，" name" ，" TEXT" )  
# insert data records  
cursor = ArcPy． da． InsertCursor ( fc，( " num" ，" name" ，"  
SHAPE@ " ) )  
#
， ， ，  
添加需要的属性字段名称 类型 长度 别名  
ArcPy． AddField_management ( fc，" num"，" SHORT"，""，""，"  
arrays = ArcPy． Array( )  
")  
序号  
pnt = ArcPy． Point( )  
ArcPy． AddField_management ( fc，" name"，" TEXT"，""，""，"  
for line in f． readlines( ) :  
")  
名称  
j = line． split( ＇，＇)  
#
插入一条记录  
cursor = ArcPy． da． InsertCursor ( fc，( " num" ，" name" ，"  
SHAPE@ " ) )  
pnt = ArcPy． Point( )  
for px in range( 2，len( j) － 1，2) :  
pnt． num，pnt． name，pnt． X，pnt． Y = j［0］，j［1］，  
j［px］，j［px + 1］  
printpnt． num，pnt． X，pnt． Y  
for line in f． readlines( ) :  
arrays． append( pnt)  
j = line． split( ＇，＇)  
polys = ArcPy． Polygon( arrays)  
pnt． num，pnt． name，pnt． X，pnt． Y = j［0］，j［1］，float( j  
［2］) ，float( j［3］)  
row =［pnt． num，pnt． name，polys］  
cursor． insertRow( row)  
points = ArcPy． Point( pnt． X，pnt． Y)  

相关下载

全部评论(0)

暂无评论
评论赚积分>>

上传资源 上传优质资源有赏金

热门标签

更多>>

最新上传

  • 打赏
  • 30日榜单

热门下载

推荐下载

本站上的所有资源均为源于网上收集或者由用户自行上传，仅供学习和研究使用，无任何商业目的，版权归原作如有侵权，请 来信指出，本站将立即改正。

ICP许可证号：京ICP证070360号     21IC电子网 2000- 版权所有

京ICP备11013301号

京公网安备 11010802024343号