從CA狀態平面坐標到Lat的PyProj轉換，Long不是預期的

Question

我正在嘗試使用pyproj將“加利福尼亞州立平面區域5”的XY坐標轉換為“經度”。

為了進行驗證，我知道CA 5區XY坐標（6559361.78613，1834842.95456）是地址13422 Ankerton St，Bassett，CA。 緯度，經度應為（34.0342403°，-118.0076074°）。

inProj = Proj(init='epsg:2229')
outProj = Proj(init='epsg:4326')
x1,y1 = x1,y1=6559361.78613, 1834842.95456
LONGITUDE,LATITUDE = transform(inProj,outProj,x1,y1)
print(LATITUDE,LONGITUDE)

但是，輸出為(34.65142393815357°, -65.96879755500356°) ，它位於大西洋上某個遠離加利福尼亞5區的地方。

Answer 1

inProj = Proj(init='epsg:2229', preserve_units=True)
outProj = Proj(init='epsg:4326')
x1,y1 = x1,y1=la_df['X_COORDINATE'][0], la_df['Y_COORDINATE'][0]
LONGITUDE,LATITUDE = transform(inProj,outProj,x1,y1)
print(LATITUDE,LONGITUDE)

事實證明pyproj假設您是以米為單位工作。 要繼續使用英制，您必須包括選項preserve_units=True 。