如何从GeoPandas MultiPolygon几何列中去除SRID=4326

我用.shp文件中的空间数据，通过GeoPandas转换为MultiPolygon，要导入PostgreSQL表中类型为geometry(multipolygon, 4326) NULL的列。目前遇到的问题是生成的几何数据带有SRID=4326;前缀，比如格式为SRID=4326;MULTIPOLYGON(((-113.82162 31.50802........，但我只需要后面的MULTIPOLYGON(((-113.82162 31.50802........部分。刚接触GeoPandas，附上我的代码：

gdf = gpd.read_file(shapefile_path).to_crs(epsg=4326)
gdf['geometry'] = gdf['geometry'].apply(lambda x: MultiPolygon([x]) if isinstance(x, Polygon) else x)
gdf['geometry'] = gdf['geometry'].apply(lambda x: WKTElement(x.wkt, srid=4326))
gdf.to_sql(table_name, engine, schema='shapes', if_exists='append', index=False, method='multi', dtype={'geometry': Geometry('MULTIPOLYGON', srid=4326)})

问题出在你用WKTElement包装WKT字符串时，它默认会带上SRID前缀。其实GeoPandas配合sqlalchemy的Geometry类型导入PostgreSQL时，完全不需要手动转成WKTElement，直接用原生几何对象即可。

修改后的代码如下：

from shapely.geometry import Polygon, MultiPolygon
import geopandas as gpd
from sqlalchemy import create_engine
from geoalchemy2 import Geometry

# 读取shp文件并转换到4326坐标系
gdf = gpd.read_file(shapefile_path).to_crs(epsg=4326)

# 将单个Polygon转为MultiPolygon
def convert_to_multipolygon(geom):
    if isinstance(geom, Polygon):
        return MultiPolygon([geom])
    return geom

gdf['geometry'] = gdf['geometry'].apply(convert_to_multipolygon)

# 直接导入数据库，无需手动处理WKT
engine = create_engine('postgresql://user:password@host:port/dbname')
gdf.to_sql(
    table_name,
    engine,
    schema='shapes',
    if_exists='append',
    index=False,
    method='multi',
    dtype={'geometry': Geometry('MULTIPOLYGON', srid=4326)}
)

关键说明 #

• 移除WKTElement转换步骤：GeoPandas会自动处理几何对象到PostGIS兼容格式的转换，不需要手动生成带SRID前缀的WKT。
• Geometry dtype已经指定了srid=4326，PostgreSQL会自动识别坐标系，无需在WKT中额外标注。

如果一定要手动处理WKT字符串（不推荐，易出错），可以用以下方式去掉前缀：

gdf['geometry'] = gdf['geometry'].apply(lambda x: x.wkt.split(';')[-1] if 'SRID=' in x.wkt else x.wkt)

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