读取googel earth 路径

24 篇文章 0 订阅
23 篇文章 0 订阅
import re,arcpy
from arcpy import env
env.workspace = 'E:\\bao1'
s=arcpy.GetParameterAsText(0)
s=s.split('\\')[-1]
inputname=s.split('.')[0]
reader=open(arcpy.GetParameterAsText(0))
prjFile=arcpy.GetParameterAsText(1)
pat=re.compile('^<coordinates>$')
while True:
    line=reader.readline()
    if len(line)==0:
        break
    line=line.strip()
    m=pat.match(line)
    if m:
        corstr=reader.readline()
        corstr=corstr.strip()
reader.close()
coordilst=corstr.split(' ')
coordif=[]
for s in coordilst:
    xyz=s.split(',')
    xyz=[float(a) for a in xyz]
    coordif.append(xyz)
point=arcpy.Point()
array=arcpy.Array()
featureList=[]
for feature in coordif:
    point.X=feature[0]
    point.Y=feature[1]
    array.add(point)
polyline=arcpy.Polyline(array,prjFile)
featureList.append(polyline)
outputname=inputname+'.shp'
arcpy.CopyFeatures_management(featureList,outputname)
arcpy.AddField_management(outputname,'Name','TEXT',9,'','','Name','NULLABLE','REQUIRED')
rows=arcpy.UpdateCursor(outputname)
for row in rows:
    row.Name=inputname
    rows.updateRow(row)
del row,rows
arcpy.SetParameter(2,outputname)
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值