基于python语言的坐标信息转换为shapefile

# -*- coding: cp936 -*-
import numpy
import arcpy
from arcpy import env
env.workspace = “C:\\Users\\kathy\\Desktop\\China“

env.overwriteOutput = True #为了避免输出文件名重复出现异常，设置文件可覆盖
input_file = open（“c:\\Users\\kathy\\Desktop\\China\\China_city.txt““r“）
IDs = []
names = []
Lons = []
Lats = []
PGList = []
pointchange = []
spatialRef1 = arcpy.SpatialReference（“C:\\Users\\kathy\\Desktop\\China\\China_geo.prj“）
spatialRef2 = arcpy.SpatialReference（“C:\\Users\\kathy\\Desktop\\China\\China_prj.prj“）
for s in input_file:
    IDnameLonLat = s.split（““）#用，将txt中的的信息隔开
    IDs.append（ID）
    names.append（name）
    Lons.append（Lon）
    Lats.append（Lat）
    point = arcpy.Point（LonLat）
    pointGeometry = arcpy.PointGeometry（pointspatialRef1）#建立几何对象，使其坐标系统为spatialRef1
    pointGeometry2 = pointGeometry.projectAs（spatialRef2）#转换几何