如何用PyMongo/GeoJSON定义圆/球形地理围栏？求Python实现方案

确实，GeoJSON标准本身没有原生支持圆形或球形地理围栏的类型，不过在Python生态里有不少成熟的工具可以解决这个问题，下面分享几种实用的方案：

1. 用Shapely + GeoPandas生成近似圆形的GeoJSON多边形 #

因为GeoJSON只支持多边形、线、点等几何类型，所以我们可以用高分辨率多边形来模拟圆形。如果需要考虑地球曲率（球形场景），关键是要处理好坐标系转换：

from shapely.geometry import Point
import geopandas as gpd
from pyproj import CRS

# 定义热点中心点（WGS84坐标系，即GPS常用的经纬度）
center_point = Point(-122.4194, 37.7749)  # 示例：旧金山坐标
radius_in_meters = 1000  # 热点半径，单位米

# 创建GeoSeries并设置坐标系
geo_series = gpd.GeoSeries([center_point], crs=CRS.from_epsg(4326))

# 关键步骤：先转成UTM投影（适配局部区域的平面坐标系，避免曲率误差），创建缓冲区，再转回WGS84
utm_crs = geo_series.estimate_utm_crs()
buffered_geo = geo_series.to_crs(utm_crs).buffer(radius_in_meters).to_crs(CRS.from_epsg(4326))

# 转换为GeoJSON格式
circle_geojson = buffered_geo.to_json()
print(circle_geojson)

这段代码生成的GeoJSON是一个多边形，顶点数量由Shapely的默认分辨率控制，你也可以在buffer()方法里传入resolution参数调整（值越大，顶点越多，形状越接近圆形）。

2. 直接判断点是否在球形热点内（无需生成GeoJSON） #

如果你的需求只是判断用户指定点是否在热点范围内，不需要生成围栏的GeoJSON，那么用geopy的距离计算会更直接高效：

from geopy.distance import distance

# 热点中心点（纬度在前，经度在后）
hotspot_center = (37.7749, -122.4194)
radius_meters = 1000

# 用户指定的目标点
user_target_point = (37.7760, -122.4200)

# 计算两点间的球面距离，判断是否在半径内
if distance(hotspot_center, user_target_point).meters <= radius_meters:
    print("该点位于热点范围内")
else:
    print("该点不在热点范围内")

geopy默认使用的是WGS84椭球计算球面距离，完全适配球形场景的需求。

3. 用PyProj直接计算球形缓冲区 #

如果你需要更底层的控制，可以用pyproj的Geod类来生成球形缓冲区的顶点，再手动构建GeoJSON多边形：

from pyproj import Geod
import json

geod = Geod(ellps="WGS84")
center_lon, center_lat = -122.4194, 37.7749
radius_meters = 1000
num_vertices = 36  # 生成36个顶点，足够接近圆形

# 生成缓冲区的顶点坐标
lon_list, lat_list, _ = geod.fwd(
    [center_lon] * num_vertices,
    [center_lat] * num_vertices,
    range(0, 360, int(360/num_vertices)),
    [radius_meters] * num_vertices
)

# 构建GeoJSON多边形（闭合，首尾顶点一致）
polygon_coords = list(zip(lon_list, lat_list))
polygon_coords.append(polygon_coords[0])
geojson = {
    "type": "Feature",
    "properties": {},
    "geometry": {
        "type": "Polygon",
        "coordinates": [polygon_coords]
    }
}

print(json.dumps(geojson, indent=2))

这些方案都能很好地替代你提到的JavaScript库的功能，根据你的具体需求（生成GeoJSON/仅判断点范围）选择最适合的即可。

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