arcgis 经纬度转大地坐标_大地经纬度转UTM

在日常GNSS应用中，常需要应用到UTM坐标，而python提供了较为完备坐标转换工具，无需再去用arcgis了。首先贴出几个帖子，值得一读，

https://blog.csdn.net/NobodyWu/article/details/81158298​blog.csdn.net https://blog.csdn.net/hqs_09/article/details/100624829​blog.csdn.net

其次，说一下要用到的工具，pyproj，其源于proj4库：

安装 : pip install pyproj

调用：import pyproj

再者，库函数使用，先上代码：

    def WGS2UTM(lon,lat):
    #WGS84转UTM
    p1 = pyproj.Proj(init='epsg:4326')
    p2 = pyproj.Proj(init='epsg:32748')
    x1, y1 = pyproj.transform(p1, p2, lon, lat)
    return x1, y1

解释：p1 为现有坐标系统，p2为目标坐标系统

init值代表坐标系统标号，这个可以去Coordinate Systems Worldwide上查询。

transform函数可以直接转出目标坐标，同时lon,lat可以是列表值。

epsg中有"epsg:3857"与常规的UTM不一致，他是Web UTM，其xy也与正常地理坐标不一致，其经度是x轴，纬度是y轴。