使用PyProj将NAD27+NGVD29高程转换为NAD27+NAVD88高程遇问题求助

问题分析 #

你遇到的核心问题是PyProj默认未加载垂直基准转换所需的地理格网数据，导致仅完成了水平坐标转换，垂直高程未发生变化。

修正方案 #

方法1：使用PROJ字符串明确指定格网 #

直接通过PROJ字符串定义包含垂直转换格网的CRS，确保PyProj能调用正确的转换数据：

from pyproj import Transformer, CRS

# 输入CRS：NAD27经纬度 + NGVD29高程（美制英尺）
crs_in = CRS.from_proj4("+proj=longlat +datum=NAD27 +vunits=us-ft +no_defs +geoidgrids=conus.gtx")
# 输出CRS：NAD27经纬度 + NAVD88高程（美制英尺）
crs_out = CRS.from_proj4("+proj=longlat +datum=NAD27 +vunits=us-ft +no_defs +geoidgrids=GEOID12B.gtx")

# 创建转换器并执行转换
transformer = Transformer.from_crs(crs_in, crs_out, always_xy=True)
lon, lat, z_in = -89.8469444, 31.23777778, 2.73
_, _, z_out = transformer.transform(lon, lat, z_in)

print(f"转换后高程：{z_out:.3f} ft")

方法2：组合水平与垂直CRS #

单独指定水平CRS和垂直CRS，再组合成3D坐标系进行转换：

from pyproj import Transformer, CRS

# 水平基准：NAD27经纬度（EPSG:4267）
horizontal_crs = CRS("epsg:4267")
# 垂直基准：NGVD29（美制英尺，EPSG:5702）、NAVD88（美制英尺，EPSG:6360）
vertical_in = CRS("epsg:5702")
vertical_out = CRS("epsg:6360")

# 组合3D输入输出CRS
crs_in = CRS(f"{horizontal_crs.to_wkt()} + {vertical_in.to_wkt()}")
crs_out = CRS(f"{horizontal_crs.to_wkt()} + {vertical_out.to_wkt()}")

transformer = Transformer.from_crs(crs_in, crs_out, always_xy=True)
lon, lat, z_in = -89.8469444, 31.23777778, 2.73
_, _, z_out = transformer.transform(lon, lat, z_in)

print(f"转换后高程：{z_out:.3f} ft")

关键注意事项 #

• 格网数据依赖：PyProj需要PROJ库的地理格网文件（如conus.gtx、GEOID12B.gtx）才能完成垂直转换。确保这些文件存在于PROJ的格网目录，或通过PROJ_DATA环境变量指定路径。
• 版本要求：使用最新版PyProj（>=3.0），旧版本对垂直CRS的支持存在缺陷。
• EPSG代码限制：部分EPSG代码默认未关联垂直转换格网，直接使用可能无法触发高程转换，需手动指定格网或组合CRS。

内容的提问来源于stack exchange，提问作者ZVY545