如何用Python+OpenCV从网格连通分量生成带洞的Shapely多边形？

完全可以通过点网格的连通分量标签生成支持孔洞的Shapely多边形，核心思路是结合OpenCV的轮廓提取能力与Shapely的多边形构建逻辑，具体步骤如下：

1. 从连通分量标签提取轮廓 #

假设你已经通过cv2.connectedComponentsWithStats得到了连通分量标签图label_img，先将目标分量对应的区域转为二值图像，再用cv2.findContours提取轮廓——必须使用cv2.RETR_CCOMP模式，它能区分外轮廓和孔洞的内轮廓：

import cv2
import numpy as np

# 目标连通分量的标签值
target_label = 1
# 生成目标分量的二值图像
binary_img = (label_img == target_label).astype(np.uint8) * 255
# 提取轮廓及层级结构（区分外轮廓/孔洞）
contours, hierarchy = cv2.findContours(binary_img, cv2.RETR_CCOMP, cv2.CHAIN_APPROX_SIMPLE)

2. 区分外轮廓与孔洞轮廓 #

hierarchy数组的每一行对应一个轮廓，其中第4个元素（索引3）表示当前轮廓的父轮廓索引：

• 外轮廓的父索引为-1
• 孔洞轮廓的父索引对应所属外轮廓的索引

outer_contours = []
hole_contours = []

for idx, cnt in enumerate(contours):
    # 筛选外轮廓
    if hierarchy[0][idx][3] == -1:
        outer_contours.append(cnt)
    # 筛选孔洞轮廓
    else:
        hole_contours.append(cnt)

3. 转换为Shapely多边形 #

OpenCV的轮廓是(N,1,2)形状的数组，需要先转为(N,2)的坐标格式，再用Shapely的Polygon类构建——孔洞需作为第二个参数传入，同时要确保内外轮廓的方向相反（外轮廓逆时针、孔洞顺时针，或反之）：

from shapely.geometry import Polygon

# 处理外轮廓坐标
outer_coords = outer_contours[0].reshape(-1, 2)
# 处理孔洞坐标，调整方向以符合Shapely要求
hole_coords_list = []
for hole_cnt in hole_contours:
    hole_coords = hole_cnt.reshape(-1, 2)
    # 检查并反转孔洞坐标方向，确保与外轮廓方向相反
    hole_poly = Polygon(hole_coords)
    if hole_poly.exterior.is_ccw == Polygon(outer_coords).exterior.is_ccw:
        hole_coords = hole_coords[::-1]
    hole_coords_list.append(hole_coords)

# 生成带孔洞的多边形
final_polygon = Polygon(outer_coords, hole_coords_list)

4. 验证与优化 #

• 若存在多个不连通的外轮廓，可使用shapely.geometry.MultiPolygon组合所有多边形
• 用final_polygon.is_valid检查多边形有效性，无效时可尝试final_polygon.buffer(0)修复
• 若点网格是基于索引的坐标，可通过乘以网格分辨率转换为实际空间坐标

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