Geopandas中如何基于多点几何生成可自定义参数的alpha shape？

方案1：使用alphashape第三方库（推荐，操作最简便） #

这个库专门用于生成alpha shape，和GeoPandas兼容性极强，调用逻辑和convex_hull几乎一致：

1. 先安装依赖：

pip install alphashape

2. 调整你的代码，在原有clusters生成后添加如下逻辑即可：

import alphashape

# 自定义alpha参数，单位和当前CRS的长度单位一致（当前EPSG:3311单位为米）
# alpha越小边界越紧凑，alpha足够大时输出结果等价于凸包
custom_alpha = 800

# 生成alpha shape，格式和convex_hull完全一致
clusters['alpha_shape'] = clusters['point'].apply(
    lambda multipoint: alphashape.alphashape(multipoint, alpha=custom_alpha)
)

# 如果不想手动调alpha，可自动计算最优边界，不需要传alpha参数即可
# clusters['alpha_shape_auto'] = clusters['point'].apply(alphashape.alphashape)

生成的alpha_shape列是标准的GeoPandas几何列，支持所有常规空间操作、可视化。

方案2：自定义实现（无额外第三方库依赖） #

如果不想安装额外库，可以基于Scipy的Delaunay三角剖分手动实现alpha shape逻辑，代码如下：

from scipy.spatial import Delaunay
from shapely.geometry import MultiPoint, Polygon, MultiPolygon

def generate_alpha_shape(multipoint: MultiPoint, alpha: float) -> Polygon|MultiPolygon:
    points = np.array([[geom.x, geom.y] for geom in multipoint.geoms])
    if len(points) < 3:
        return multipoint.convex_hull
    tri = Delaunay(points)
    triangles = points[tri.simplices]
    # 计算每个三角形的外接圆半径
    a = np.linalg.norm(triangles[:, 0] - triangles[:, 1], axis=1)
    b = np.linalg.norm(triangles[:, 1] - triangles[:, 2], axis=1)
    c = np.linalg.norm(triangles[:, 2] - triangles[:, 0], axis=1)
    s = (a + b + c) / 2.0
    area = np.sqrt(s * (s - a) * (s - b) * (s - c))
    circum_r = a * b * c / (4.0 * area + 1e-10)
    # 筛选外接圆半径小于alpha的三角形
    valid_tri = triangles[circum_r < alpha]
    if len(valid_tri) == 0:
        return multipoint.convex_hull
    # 合并有效三角形生成边界
    polygons = [Polygon(tri) for tri in valid_tri]
    return gpd.GeoSeries(polygons).unary_union

# 调用方法和之前一致
custom_alpha = 800
clusters['alpha_shape'] = clusters['point'].apply(lambda x: generate_alpha_shape(x, custom_alpha))

注意事项 #

• alpha的取值需要和你使用的CRS单位匹配，比如你当前用的EPSG:3311单位为米，设置alpha=1000代表仅保留外接圆半径小于1000米的三角形
• alpha取值过小时会出现边界破碎、生成多多边形的情况，取值过大时输出结果会和凸包完全一致，可以根据你的聚类点密度逐步调整到合适值

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