如何在大型GeoPandas DataFrame中将Shapely几何对象转换为GeoJSON字符串并避免复制警告

核心解决方案 #

你遇到的SettingWithCopyWarning是因为直接对切片副本赋值导致的，而且没必要绕GeoSeries来转换单个几何对象——shapely本身就提供了直接转换的方法，同时我们可以通过at/loc直接操作原DataFrame，或者用apply批量处理，完全不需要创建新的DataFrame。

1. 单个元素的安全赋值

如果你只需要修改某一行的geometry，用at直接定位原DataFrame的单元格，结合shapely的to_json函数转换：

from shapely import to_json

# 直接修改原DataFrame的第0行geometry列
file.at[0, 'geometry'] = to_json(file.loc[0, 'geometry'])

at是pandas中定位单个单元格最高效的方式，它直接操作原DataFrame，不会产生副本，因此不会触发警告。

2. 批量处理整个Geometry列（适合大规模数据集）

如果要把整个列的shapely对象都转成GeoJSON字符串，用apply直接修改原列即可，效率很高：

from shapely import to_json

# 直接替换原geometry列，无需创建新DataFrame
file['geometry'] = file['geometry'].apply(to_json)

这样处理后，geometry列的每个元素都会变成对应的GeoJSON几何字符串（比如{"type": "Polygon", "coordinates": [...]}），完全符合你的需求。

为什么你的原代码会触发警告？ #

file.iloc[0]返回的是原DataFrame的切片副本（pandas无法确定是视图还是副本，所以抛出警告），你在这个副本上修改geometry，可能不会同步到原DataFrame，这就是警告的原因。而用at/loc直接定位原DataFrame的单元格，或者直接对整列用apply操作，都是直接修改原数据，不存在副本的问题。

额外说明 #

如果你需要的是包含属性的Feature格式GeoJSON（而不是单纯的几何），可以用geopandas.GeoSeries的to_json方法指定输出格式，但针对单个元素的话，也可以手动构造：

# 单个元素转Feature格式GeoJSON
feature_json = geopandas.GeoSeries(file.loc[0]).to_json(na='drop')
file.at[0, 'geometry'] = feature_json

不过根据你的需求描述，单纯的几何GeoJSON字符串应该就足够了。

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