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

這篇文章概述了問題第二部分的答案。