python - PyProj：使用 x,y 坐标将纬度/经度映射到矩形

Question

我有一系列纬度/经度坐标，我试图将它们作为 x,y 坐标投影到地图上。

我正在使用https://code.google.com/p/pyproj/库来转换坐标，但我似乎找不到任何解释函数输出含义的地方？

import PyProj

def ConvertToMapProjection(Coordinates):
    ''' Projects a mapped pair of coordinates onto a map '''

    # Define the projection
    RobinsonProjection = pyproj.Proj("+proj=robin +lon_0=0 +x_0=0 +y_0=0 +ellps=WGS84 +units=m +no_defs")

    East, North = RobinsonProjection(Coordinates[0],Coordinates[1])

    return [East,North]

ConvertToMapProjection([51.5072,0.1275])

返回：

[4866232.474090106, 13636.369990048854]

不清楚单位是什么？将其映射到 900 x 1100 矩形的最佳方法是什么？

Answer 1

proj的默认单位是米（您甚至可以在命令中指定它）。

您使用的参考椭球体是 WGS84，其平均半径为 6356752.3142 米（参见 参考资料proj -le）

您可以按此因子按比例缩小以获得“单位椭圆体”：

$ proj +proj=robin +ellps=WGS84 -m 1:6356752.3142

使用单位球体进行投影更容易，但不是那么精确：

$ proj +proj=robin +a=1

这篇文章概述了您问题第二部分的答案。