geopandas.from_postgis在Python文件中触发Unknown WKB type 41错误

核心原因 #

WKB类型41对应PostGIS的三维(Z)几何类型，看似相同的包路径下，原生Python环境和Jupyter环境的shapely/GEOS库可能存在隐性差异——比如系统级GEOS库冲突、shapely安装方式（pip/conda）不同，导致原生环境的shapely无法解析三维WKB。

1. 强制转换几何类型兼容shapely #

直接在SQL查询里把三维几何转成二维，或者输出shapely兼容的WKB格式：

• 转二维（无需保留Z值）：

import geopandas as gpd
from sqlalchemy import create_engine

engine = create_engine("postgresql://user:pass@host:port/db")
# 用ST_Force2D把三维几何转成二维
sql = "SELECT ST_Force2D(geom_column) AS geom, other_columns FROM your_table"
gdf = gpd.GeoDataFrame.from_postgis(sql, engine, geom_col="geom")

• 保留Z值并输出兼容格式：

sql = "SELECT ST_AsEWKB(geom_column) AS geom, other_columns FROM your_table"
gdf = gpd.GeoDataFrame.from_postgis(sql, engine, geom_col="geom", crs=4326)

2. 核对shapely与GEOS的版本匹配 #

Jupyter可能用conda自带的高版本GEOS，而原生Python调用了系统低版本GEOS（低于3.8.0不支持三维WKB）。在两个环境分别执行以下代码对比版本：

import shapely
print(shapely.__version__)
print(shapely.geos_version)

如果原生环境GEOS版本过低，升级shapely到和geopandas 0.12.2兼容的稳定版：

# pip安装方式
pip install --upgrade shapely==1.8.5
# conda环境优先用conda升级，避免混合安装冲突
conda install shapely=1.8.5

3. 手动注册PostGIS适配器 #

原生Python的sqlalchemy连接可能未正确加载PostGIS类型适配器，导致WKB格式异常。手动注册适配器后再查询：

from sqlalchemy import create_engine
from geoalchemy2.types import _GEOBASE

engine = create_engine("postgresql://user:pass@host:port/db")
# 显式注册PostGIS类型解析器
_GEOBASE.registry.register(type(_GEOBASE), None)

sql = "SELECT geom_column, other_columns FROM your_table"
gdf = gpd.GeoDataFrame.from_postgis(sql, engine, geom_col="geom_column")

4. 排查环境变量差异 #

Jupyter可能加载了额外环境变量（比如Linux的LD_LIBRARY_PATH、Mac的DYLD_LIBRARY_PATH），导致调用不同的GEOS库。在原生脚本开头打印环境变量对比：

import os
print(os.environ.get("LD_LIBRARY_PATH"))

如果发现路径差异，在原生脚本开头手动指定正确的GEOS库路径（比如conda环境的lib目录）：

import os
# 替换为你的conda环境lib路径
os.environ["LD_LIBRARY_PATH"] = "/path/to/your/conda/env/lib:" + os.environ.get("LD_LIBRARY_PATH", "")
# 必须在设置环境变量后再导入shapely
import shapely

内容的提问来源于stack exchange，提问作者Kabrüggen