从 CA 状态平面坐标到 Lat、Long 的 PyProj 转换不符合预期

Question

我正在尝试使用 pyproj 将加利福尼亚州平面区域 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 假定您以米为单位工作。要继续使用英制系统，您必须包含选项 preserve_units=True.

从 CA 状态平面坐标到 Lat、Long 的 PyProj 转换不符合预期

PyProj transformation from CA state plane coordinates to Lat, Long not as expected

python

gis

latitude-longitude

python-3.x

pyproj