因此，我只是为PostGIS数据库实现了这一点，可以在这里粘贴我的方法。对于MySQL，您必须修改代码。
第一步是在地理编码列转换为WKB十六进制字符串，因为我使用SQLAlchemy的，基于发动机pyscopg，并且这两个包的不理解地理类型本身。下一步是照常将数据写入SQL
DB（请注意，所有几何列都应转换为包含WKB十六进制字符串的文本列），最后通过执行查询将列的类型更改为几何。请参考以下伪代码：
# importsimport sqlalchemy as salimport geopandas as gpd# Function to generate WKB hexdef wkb_hexer(line): return line.wkb_hex# Convert `'geom'` column in GeoDataframe `gdf` to hex # Note that following this step, the GeoDataframe is just a regular Dataframe # because it does not have a geometry column anymore. Also note that # it is assumed the `'geom'` column is correctly datatyped.gdf['geom'] = gdf['geom'].apply(wkb_hexer)# Create SQL connection engineengine = sal.create_engine('postgresql://username:password@host:socket/database')# Connect to database using a context managerwith engine.connect() as conn, conn.begin(): # Note use of regular Pandas `to_sql()` method. gdf.to_sql(table_name, con=conn, schema=schema_name, if_exists='append', index=False) # Convert the `'geom'` column back to Geometry datatype, from text sql = """ALTER TABLE schema_name.table_name ALTER COLUMN geom TYPE Geometry(LINESTRING, <SRID>) USING ST_SetSRID(geom::Geometry, <SRID>)""" conn.execute(sql)
