使用geopandas写入PostGIS后无法识别geometry类型如何解决

1. 自定义查询报ValueError: Query missing geometry column原因 #

你执行的SQL查询语句仅查询了st_cnty_fips字段，没有包含geometry列。gpd.read_postgis()会将查询结果转换为GeoDataFrame，要求返回结果必须包含你通过geom_col参数指定的空间列，和数据库内是否存在该列无关，只要查询结果里没有就会触发报错。

2. 读全表时报SAWarning: Did not recognize type 'geometry'原因 #

你创建的SQLAlchemy引擎没有注册PostGIS的geometry空间类型，在自动反射表结构时无法识别该字段类型，因此抛出警告，通常是缺少GeoAlchemy2依赖、驱动声明不明确导致的。

• 解决查询报错：
  • 如果你需要返回带空间信息的GeoDataFrame，修改SQL语句，在SELECT中加入geometry列：

    gpd.read_postgis('SELECT DISTINCT st_cnty_fips, geometry FROM table_name', con=self.engine, geom_col='geometry')
    

  • 如果你仅需要查询属性字段不需要空间信息，直接用pandas的read_sql()方法即可，无需调用geopandas的空间读取接口：

    import pandas as pd
    df = pd.read_sql('SELECT DISTINCT(st_cnty_fips) FROM table_name', con=engine)
    

• 解决类型识别警告：
  1. 先安装必要依赖：

     pip install psycopg2-binary geoalchemy2
     

  2. 创建引擎时明确指定psycopg2驱动，SQLAlchemy会自动关联GeoAlchemy2的空间类型注册逻辑：

     engine = create_engine(f'postgresql+psycopg2://user:password@localhost:5432/{dbname}')
     

  3. 如果仍有警告，读取时可手动指定geometry列的类型，对应填写你数据的SRID和几何类型即可：

     from geoalchemy2 import Geometry
     gdf = gpd.read_postgis(
         'table_name',
         con=engine,
         geom_col='geometry',
         dtype={'geometry': Geometry(srid=4326, geometry_type='POLYGON')}
     )
     

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