如何从x,y笛卡尔坐标tmerc投影得到Lat / Lon?

时间:2017-11-08 09:10:08

标签: python coordinates gis proj pyproj

我正在尝试将数据从OpenDrive笛卡尔坐标转换为Lat / Lon值。对于德国A9高速公路上的pt = [6.714150516498e + 05,5.434880530093e + 06],我想知道Lat / Lon坐标是什么。

我一直试图以这种方式解决这个问题,但是我在尼日利亚的沙漠中得到了拉脱/拉力坐标。

from pyproj import Proj, transform
inProj = Proj("+proj=tmerc +lat_0=0 +lon_0=9 +k=0.9996 +x_0=500000 +y_0=0 +datum=WGS84 +units=m +no_defs")
outProj = Proj(init='epsg:4326', preserve_units=True)

w = [6.714150516498e+05, 5.434880530093e+06]

lat, lon = transform(inProj, outProj, w[0], w[1])
print(lat, lon)

1 个答案:

答案 0 :(得分:4)

transform返回(x,y)。这意味着:x = lon, y = lat

>>> print transform(inProj, outProj, w[0], w[1])[::-1]
(49.04294637744738, 11.345684678187824)

德国的某个地方就在A9高速公路旁边。