geometry.py 文件源码

def compute_buffered_area_of_stops(lats, lons, buffer_meters, resolution=16):
    # geo_series = gp.GeoSeries([Point(lon, lat) for lon, lat in zip(lons, lats)])
    # geo_series.crs = {'init' :'epsg:4326'}
    # geo_series = geo_series.to_crs({'init':'epsg:3857'})

    # circles = geo_series.buffer(buffer_meters, resolution=resolution)
    # multi_points = circles.unary_union
    # return multi_points.area

    if len(lons) > 1:
        lon_meters, lat_meters = _get_lon_lat_meters(lons, lats)
    else:
        lon_meters = lons
        lat_meters = lats

    return MultiPoint(points=list(zip(lon_meters, lat_meters))).buffer(buffer_meters, resolution=resolution).area