【问题标题】:How to get Lat/Lon from x,y Cartesian Coordinates tmerc projection?如何从 x,y 笛卡尔坐标 tmerc 投影中获取纬度/经度?
【发布时间】:2018-04-20 21:31:52
【问题描述】:

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

我一直在尝试以这种方式解决它,但我得到了尼日利亚沙漠的纬度/经度坐标。

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)

【问题讨论】:

  • 我和他们核实了这个地方,我还找到了 u'Okha, Ovia South-West, Edo, Nigeria' 你能再检查一下纬度和经度吗?

标签: python coordinates gis proj pyproj


【解决方案1】:

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

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

位于德国 A9 高速公路旁的某处。

【讨论】:

    猜你喜欢
    • 2020-07-04
    • 2013-02-07
    • 2010-11-14
    • 2017-06-29
    • 1970-01-01
    • 1970-01-01
    • 1970-01-01
    • 1970-01-01
    • 1970-01-01
    相关资源
    最近更新 更多