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 

的回答你的問題的第二部分：更容易，但並非如此精確，將使用單位球被突出3210