Convert to and from UTM coordinates with pyproj

最新推荐文章于 2024-09-06 22:47:21 发布

x5675602

最新推荐文章于 2024-09-06 22:47:21 发布

阅读量515

点赞数

分类专栏：遥感

本文链接：https://blog.csdn.net/x5675602/article/details/87992112

版权

遥感专栏收录该内容

31 篇文章 9 订阅

订阅专栏

相应的python的代码如下所示：

from pyproj import Proj
import numpy as np
from pandas import DataFrame

x = (-909419.00, 529114.00, 545227.00, 545582.00)
y = (4536905, 7422343.00, 7435702.00, 7435741.00)

df = DataFrame(np.c_[x, y], columns=['Meters East', 'Meters South'])
print 'The original points:\n', df
myProj = Proj("+proj=utm +zone=13N, +north +ellps=WGS84 +datum=WGS84 +units=m +no_defs")
lon, lat = myProj(df['Meters East'].values, df['Meters South'].values, inverse=True)
df = DataFrame(np.c_[x, y, lon, lat], columns=['UTMx', 'UTMy', 'Lon', 'Lat'])

print 'The result (x, y) --> (lon, lat)\n', df
UTMx, UTMy = myProj(lon, lat)
print(UTMx, UTMy)
df = DataFrame(np.c_[UTMx, UTMy, lon, lat], columns=['UTMx', 'UTMy', 'Lon', 'Lat'])
print 'The result (lon, lat) --> (x, y)\n', df