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)
读取googel earth 路径
最新推荐文章于 2023-08-18 15:55:58 发布