如何避免Shapely生成的多边形导入BigQuery后变为NULL？

核心原因 #

直接在WGS84经纬度（EPSG:4326）坐标系上调用Shapely的buffer()是平面运算，生成的多边形容易出现自相交、不符合球面几何规范的情况，即便BigQuery的SAFE.ST_GEOGFROMGEOJSON开启make_valid参数，也无法修复所有无效几何，最终导致部分条目返回NULL。

具体解决方案 #

1. 投影转换后再做Buffer（推荐） #

先将经纬度点转换为平面坐标系完成buffer运算，再转回WGS84，避免球面直接运算的偏差：

from shapely.geometry import Point
from shapely.ops import transform
import pyproj

# 定义投影转换：WGS84转全球通用平面投影（EPSG:3857，单位米）
proj_to_plane = pyproj.Transformer.from_crs(
    "EPSG:4326", "EPSG:3857", always_xy=True
).transform
# 定义逆转换：平面投影转回WGS84
proj_to_wgs84 = pyproj.Transformer.from_crs(
    "EPSG:3857", "EPSG:4326", always_xy=True
).transform

# 生成点并完成投影转换与buffer
point_wgs84 = Point(lng, lat)
point_plane = transform(proj_to_plane, point_wgs84)
polygon_plane = point_plane.buffer(r)  # r单位为米
polygon_wgs84 = transform(proj_to_wgs84, polygon_plane)

2. 提前校验并修复Shapely几何 #

在转换为GeoJSON前，主动修复无效几何，避免导入BigQuery时解析失败：

# 修复自相交等常见无效问题
if not polygon_wgs84.is_valid:
    polygon_wgs84 = polygon_wgs84.buffer(0)

# 统一输出为Polygon类型（避免buffer生成MultiPolygon）
if polygon_wgs84.geom_type == "MultiPolygon":
    polygon_wgs84 = list(polygon_wgs84.geoms)[0]

3. 直接用BigQuery原生函数生成缓冲区 #

跳过Shapely的本地运算，直接在BigQuery中基于经纬度点生成有效缓冲区：

SELECT ST_BUFFER(ST_GEOGPOINT(lng, lat), r) AS polygon
FROM your_table

注：BigQuery的ST_BUFFER默认单位为米，原生支持球面几何运算，生成的GEOGRAPHY类型100%有效。

4. 校验GeoJSON格式合法性 #

确保转换后的GeoJSON符合规范，避免格式错误导致解析失败：

import geojson

# 转换为GeoJSON并验证
geo_json_str = geojson.dumps(geojson.Feature(geometry=polygon_wgs84))
try:
    geojson.loads(geo_json_str)
except ValueError as e:
    print(f"无效GeoJSON: {e}")

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