Arcpy读取excel生成shp文件

标签: Arcgis 分类: Gis 创建时间:2019-10-24 06:21:13 更新时间:2025-01-17 10:39:21

excel表中存储的内容格式如下:

读取excel文件使用了openpyxl库,定义点arcpy.Point,将点存储在数组中,然后将点数组构建成线。

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
#!/usr/bin/python
# -*- coding: UTF-8 -*-

import openpyxl as excel
import arcpy
import os

## 读取excel生成shp文件
def excelToShp(path):

arcpy.env.workspace = os.getcwd()
outName = r"polyline.shp"
polylineGeometryList = [] # a list to hold the PointGeometry objects
spRef = arcpy.SpatialReference(4326)


## 读取excel,获取站点关联线段
excelFile=path
wb=excel.load_workbook(excelFile)
## 读取第一行,获取列头
sheet=wb["Sheet1"]
rows=sheet.max_row
columns=sheet.max_column
dic={}
## 存储列名与列号的对应
for i in range(1,columns+1):
dic[sheet.cell(row=1,column=i).value]=i

## 遍历剩余的行,根据对应关系求坐标
polylinGeometryList=[]
for i in range(2,rows+1):
start_n=str(sheet.cell(i,dic['start_n']).value)

if start_n is None:
continue

## 开始坐标和结束坐标
start_lon=sheet.cell(i,dic['start_lon']).value
start_lat=sheet.cell(i,dic['start_lat']).value
end_lon=sheet.cell(i,dic['end_lon']).value
end_lat=sheet.cell(i,dic['end_lat']).value
## 构建线段
startPoint=arcpy.Point(float(start_lon), float(start_lat))
endPoint=arcpy.Point(float(end_lon), float(end_lat))
array=arcpy.Array()
array.append(startPoint)
array.append(endPoint)
polylinGeometry = arcpy.Polyline(array,spRef)
polylinGeometryList.append(polylinGeometry)

## 保存内容
arcpy.CopyFeatures_management(polylinGeometryList, outName)

print 'done'



if __name__=='__main__':

excelPath=r"D:\zlc\hdxs\filedata\creat_test.xlsx"
excelToShp(excelPath)

其他知识点
(1) 官方文档中,有这么一个高级的python用法,叫列表推导式,加了星号(*)的变量名会存放所有未命名的变量参数,称为不定长参数。

我觉得就类似于javascript中的map函数,前面是对每一个对象所作的事,后面是在哪里遍历对象。

(2) 关于python的路径问题

1
2
3
4
5
6
7
8
9
10
import os
## 创建路径
if os.path.isdir(path + '/' + x):
print(x + ' exists.')
else:
os.mkdir(path + '/' + x)
print(x + ' created.')

## 拼接路径
os.path.join('/my/root/directory', 'in', 'here')
小额赞助
本人提供免费与付费咨询服务,感谢您的支持!赞助请发邮件通知,方便公布您的善意!
**光 3.01 元
Sun 3.00 元
bibichuan 3.00 元
微信公众号
广告位
诚心邀请广大金主爸爸洽谈合作
每日一省
isNaN 和 Number.isNaN 函数的区别?

1.函数 isNaN 接收参数后,会尝试将这个参数转换为数值,任何不能被转换为数值的的值都会返回 true,因此非数字值传入也会返回 true ,会影响 NaN 的判断。

2.函数 Number.isNaN 会首先判断传入参数是否为数字,如果是数字再继续判断是否为 NaN ,不会进行数据类型的转换,这种方法对于 NaN 的判断更为准确。

每日二省
为什么0.1+0.2 ! == 0.3,如何让其相等?

一个直接的解决方法就是设置一个误差范围,通常称为“机器精度”。对JavaScript来说,这个值通常为2-52,在ES6中,提供了Number.EPSILON属性,而它的值就是2-52,只要判断0.1+0.2-0.3是否小于Number.EPSILON,如果小于,就可以判断为0.1+0.2 ===0.3。

每日三省
== 操作符的强制类型转换规则?

1.首先会判断两者类型是否**相同,**相同的话就比较两者的大小。

2.类型不相同的话,就会进行类型转换。

3.会先判断是否在对比 null 和 undefined,是的话就会返回 true。

4.判断两者类型是否为 string 和 number,是的话就会将字符串转换为 number。

5.判断其中一方是否为 boolean,是的话就会把 boolean 转为 number 再进行判断。

6.判断其中一方是否为 object 且另一方为 string、number 或者 symbol,是的话就会把 object 转为原始类型再进行判断。

每日英语
Happiness is time precipitation, smile is the lonely sad.
幸福是年华的沉淀,微笑是寂寞的悲伤。