PyProj从CA状态平面坐标转换为Lat，Long不如预期

Question

我正在尝试使用pyproj将Calfornia State Plane Zone 5 X-Y坐标转换为纬度，经度。

为了验证，我知道CA Zone 5 X-Y坐标（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假设你在米里工作。要继续使用Imperial系统，您必须包含选项preserve_units=True。