如何将GeoDataFrame插入MongoDB？insert_many方法无法生效

你的问题核心在于GeoDataFrame中的Shapely几何对象没有转换成MongoDB支持的GeoJSON格式。MongoDB的地理空间数据要求遵循GeoJSON规范，而默认的to_dict('records')会把几何对象转成字符串（比如WKT）或原生Python对象，Mongo无法识别为有效的地理空间类型，导致插入失败。

下面是具体的解决步骤和代码示例：

步骤1：将几何字段转换为GeoJSON格式 #

Shapely几何对象自带__geo_interface__属性，它正好是符合GeoJSON标准的字典结构，我们可以用这个属性来转换：

# 将geometry字段的Shapely对象转为GeoJSON字典
my_dataframe['geometry'] = my_dataframe['geometry'].apply(lambda geom: geom.__geo_interface__)

如果你的DataFrame中有多个几何字段，需要对每个字段都执行这个转换操作。

步骤2：转换为字典列表并插入MongoDB #

处理完几何字段后，再用to_dict('records')生成MongoDB可接受的字典列表，然后执行插入：

完整代码示例 #

import geopandas as gpd
from pymongo import MongoClient

# 初始化你的GeoDataFrame（假设df是原始DataFrame）
my_dataframe = gpd.GeoDataFrame(df)

# 关键：转换几何字段为GeoJSON格式
my_dataframe['geometry'] = my_dataframe['geometry'].apply(lambda geom: geom.__geo_interface__)

# 连接MongoDB（根据你的实际连接信息调整）
client = MongoClient('mongodb://localhost:27017/')
db = client['Marine_Database']

# 插入数据到集合
db.Theoritical_Coverage.insert_many(my_dataframe.to_dict('records'))

额外优化：创建地理空间索引 #

如果之后需要对地理数据进行查询，建议插入后创建2dsphere索引，这样MongoDB就能高效处理地理空间查询：

# 为geometry字段创建2dsphere索引
db.Theoritical_Coverage.create_index([('geometry', '2dsphere')])

补充说明 #

• 你的theoretica字段中存储的是WKT格式的多边形数据，如果不需要保留的话可以删除；如果需要保留，也可以用同样的方法将其转换为GeoJSON格式（先转成Shapely对象再取__geo_interface__）。
• 确保MongoDB版本支持地理空间特性（3.0以上版本都支持）。

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