如何将OpenCV检测到的轮廓转换为平滑样条曲线或类贝塞尔曲线

我正在对彩色位图进行矢量化处理：大致流程为统计颜色数量，将每种颜色提取为黑白（b/w）掩码，再检测轮廓，最后将所有内容导出为SVG文件。但目前输出效果较差，颜色边缘的轮廓过于生硬粗糙，同时还存在大量“椒盐噪声”，导致输出中出现过多无用小轮廓。

为了消除椒盐噪声，我尝试在轮廓检测前应用不同的模糊效果（median blur、gaussian），以及形态学操作（erode/dilate/open/close），但这些方法都不适用于我的场景，会模糊甚至删除需要保留清晰边缘的微小对象/轮廓。

最终我打算将检测到的轮廓点集转换为样条曲线或贝塞尔曲线，但还没找到利用cv2.findContour(...)返回的数组实现该需求的方法。

我曾试验过对轮廓应用c = cv2.approxPolyDP(c, simplification_param, False)，但处理后的轮廓反而更棱角分明、存在断裂，没有得到平滑连续的效果。

请问是否有方法可以将cv2.approxPolyDP(...)的输出结果转换为样条曲线，或者有没有其他库可以实现类似approxPolyDP的效果，但输出为贝塞尔曲线？

以下是我用于通过cv2.findContours(...)检测轮廓的黑白掩码示例，最终需要将这些轮廓转换为平滑的连续曲线（无需一定是连通或闭合轮廓）。

谢谢！

解决方案 #

1. 高效消除椒盐噪声（不损失有效微小轮廓） #

不用在掩码预处理阶段做模糊或形态学操作，直接在轮廓检测完成后过滤小面积轮廓即可，仅丢弃面积低于设定阈值的无效噪声轮廓，完全不会影响需要保留的有效轮廓：

import cv2
import numpy as np

# 读取单通道黑白掩码
mask = cv2.imread("你的掩码路径.png", 0)
# 检测轮廓，CHAIN_APPROX_NONE保留所有轮廓点，适合后续平滑拟合
contours, _ = cv2.findContours(mask, cv2.RETR_EXTERNAL, cv2.CHAIN_APPROX_NONE)

# 按面积过滤噪声轮廓，min_area阈值根据你的图像分辨率自行调整
min_area = 3
filtered_contours = [c for c in contours if cv2.contourArea(c) > min_area]

2. 轮廓转平滑样条/贝塞尔曲线 #

不需要依赖cv2.approxPolyDP的输出做转换，直接用原始轮廓点集即可做平滑拟合，以下两种方案都可实现：

方案1：用Scipy做B样条拟合（灵活度高）

通过调整平滑参数可以平衡边缘保真度和顺滑度，拟合后的样条可以直接采样导出，也可以转换为贝塞尔曲线段写入SVG：

from scipy.interpolate import splprep, splev

def get_smooth_contour(contour, smoothness=2, sample_num=200, is_closed=True):
    # 转换cv2轮廓的格式为独立的x、y数组
    x = contour[:, 0, 0]
    y = contour[:, 0, 1]
    # 拟合B样条，s为平滑权重，数值越大越平滑，per=1对应闭合轮廓
    tck, _ = splprep([x, y], s=smoothness, per=1 if is_closed else 0)
    # 采样生成平滑后的坐标点
    u = np.linspace(0, 1, sample_num)
    smooth_x, smooth_y = splev(u, tck)
    return np.column_stack((smooth_x, smooth_y)).astype(np.int32)

smoothness参数取值越小，拟合结果越贴近原始轮廓边缘，取值越大平滑效果越强，可根据你的需求灵活调整。

方案2：用OpenCV扩展模块直接输出贝塞尔控制点

如果你安装了opencv-contrib-python包，可以直接调用ximgproc模块的轮廓拟合接口，直接输出三次贝塞尔曲线的控制点，无需额外转换：

import cv2
from cv2 import ximgproc

# 初始化拟合器，指定贝塞尔曲线阶数，3为三次贝塞尔
fitter = ximgproc.createContourFitting(ctr=3)
# 输入原始轮廓点，输出贝塞尔控制点和对应系数
ctrl_points, coeffs = fitter.fit(contour, 0.01)

拟合完成后你可以直接将曲线数据按SVG规范拼接，或者用svgpathtools库直接生成SVG文件即可。

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