如何用Python（Geopandas等）创建覆盖所有建筑多边形的边界多边形

要生成覆盖所有建筑的范围多边形，你可以用两种简单高效的方法：轴对齐边界框或凸包，完全满足你"无需贴合顶点、覆盖所有范围"的需求。

核心思路 #

• 边界框：生成一个轴对齐的矩形，直接包裹所有建筑的最小/最大坐标，操作最简单，速度最快。
• 凸包：生成一个包围所有建筑的最小凸多边形，比边界框更紧凑，依然不需要贴合建筑的外围顶点。

代码实现 #

1. 安装依赖 #

先确保安装了必要的库：

pip install geopandas matplotlib shapely

2. 完整代码示例 #

import geopandas as gpd
from shapely.geometry import Polygon
import matplotlib.pyplot as plt

# 读取建筑多边形图层（支持shp、geojson等格式）
buildings = gpd.read_file("你的建筑数据文件.shp")

# 预处理：过滤无效几何（如果有的话）
buildings = buildings[buildings.is_valid]

# --- 方法1：生成轴对齐边界框 ---
# 获取所有建筑的整体边界范围
total_bounds = buildings.total_bounds  # 返回格式：[minx, miny, maxx, maxy]
# 将边界范围转为多边形对象
bbox_poly = Polygon([
    (total_bounds[0], total_bounds[1]),
    (total_bounds[2], total_bounds[1]),
    (total_bounds[2], total_bounds[3]),
    (total_bounds[0], total_bounds[3]),
    (total_bounds[0], total_bounds[1])
])
# 转为GeoDataFrame，保持和原数据一致的坐标系统
bbox_gdf = gpd.GeoDataFrame([{"geometry": bbox_poly}], crs=buildings.crs)

# --- 方法2：生成凸包多边形 ---
# 合并所有建筑多边形，再生成凸包
convex_hull = buildings.unary_union.convex_hull
convex_hull_gdf = gpd.GeoDataFrame([{"geometry": convex_hull}], crs=buildings.crs)

# --- 可视化对比效果 ---
fig, (ax1, ax2, ax3) = plt.subplots(1, 3, figsize=(18, 6))
# 原始建筑图层
buildings.plot(ax=ax1, color="#87CEEB", edgecolor="black", linewidth=0.5)
ax1.set_title("原始建筑多边形")
# 边界框覆盖效果
bbox_gdf.plot(ax=ax2, color="none", edgecolor="#FF4500", linewidth=2)
buildings.plot(ax=ax2, color="#87CEEB", edgecolor="black", linewidth=0.5)
ax2.set_title("轴对齐边界框覆盖")
# 凸包覆盖效果
convex_hull_gdf.plot(ax=ax3, color="none", edgecolor="#32CD32", linewidth=2)
buildings.plot(ax=ax3, color="#87CEEB", edgecolor="black", linewidth=0.5)
ax3.set_title("凸包覆盖")

plt.tight_layout()
plt.show()

# 保存结果到文件
bbox_gdf.to_file("建筑边界框.shp")
convex_hull_gdf.to_file("建筑凸包.shp")

选择建议 #

• 如果只需要快速获取一个矩形范围，选边界框即可，代码少、效率高。
• 如果想要更紧凑的非矩形范围，选凸包，它会自动贴合建筑群的整体轮廓，但依然是凸多边形，不会凹进去贴合单个建筑的边角。

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