基于python实现地理坐标系与投影坐标系的转换
#-*-coding=utf-8 -*-
import pyproj
import pandas as pd
import numpy as np
data = pd.read_csv("test.csv",sep=',',header='infer')#将第一行作为表头
p1 = pyproj.Proj(init="epsg:4326")#wgs坐标系统的EPSG Code
p2 = pyproj.Proj(init="epsg:32651")#UTM 51N坐标系统的EPSG Code
lonprj, latprj = pyproj.transform(p1, p2,data['lon'], data['lat'])
xy = list(zip(lonprj,latprj))#封装数组,使每一个经度对应一个纬度
np.savetxt("UTM51.txt",xy,fmt="%.5f")#小数点后保留5位
本程序通过不同坐标系的EPSG代码实现坐标系的转换。