读取googel earth 点坐标

24 篇文章 0 订阅
23 篇文章 0 订阅
import re,arcpy
from arcpy import env
env.workspace = 'E:\\bao1'#%%工作空间
#reader=open('','r')#打开点的txt文件
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>.*</coordinates>')
while True:
    line=reader.readline()
    if len(line)==0:
        break
    line=line.strip()#删除首位空格
    m=pat.match(line)
    if m:
        corstr=line
        corstr=corstr.strip()
reader.close()
corstr=re.sub('<coordinates>(.*?)</coordinates>',r'\1',corstr)
coordilst=corstr.split(',')
coordilst=[float(a) for a in coordilst]
point=arcpy.Point()
point.X=coordilst[0]
point.Y=coordilst[1]
pointGeometry=arcpy.PointGeometry(point,prjFile)
outputname=inputname+'.shp'
arcpy.CopyFeatures_management(pointGeometry,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)
在这里插入代码片
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值