Bezier曲线的公式推导及代码实现

Qmj0923

已于 2022-03-12 18:18:44 修改

阅读量621

点赞数

分类专栏： NJU计算机笔记文章标签：算法几何学 python

于 2022-03-12 17:57:53 首次发布

本文链接：https://blog.csdn.net/Qmj2333333/article/details/123445437

版权

NJU计算机笔记专栏收录该内容

3 篇文章 0 订阅

订阅专栏

本文仅简述Bezier曲线的公式推导，并给出了一种代码实现。在阅读本文之前，请确保你已经对Bezier曲线的背景知识有所了解。相关知识可以通过以下课程进行学习：MOOC-计算机图形学-中国农业大学-赵明或者观看B站搬运版。

算法原理

给定 $n + 1$ 个控制点 $P_i(i=0,1,\cdots,n)$ ，构造多项式函数
$\mathcal{P}(u)=\sum_{i=0}^{n}P_i\cdot BEZ_{i,n}(u), \quad 0\leqslant u\leqslant 1 \thinspace， \tag{1}$
来逼近曲线。其中 $BEZ_{i,n}(u)$ 称为 $n$ 次Bernstein基函数，且有
$BEZ_{i,n}(u)=C_n^i u^i(1-u)^{n-i}, \quad \left(\begin{array}{l} i=0,1,\cdots,n \\ u\in[0,1] \end{array}\right) \thinspace。 \tag{2}$

对于公式(1)，函数 $\mathcal{P}(u)$ 将实数区间 $[0, 1]$ 映射到实数域上的点集，该点集里的点组成了所求的逼近曲线。 $BEZ_{i,n}(u)$ 可以看作是控制点 $P_i$ 在 $u$ 下的权值。因为 $BEZ_{i,n}(u)$ 具有权性：
$\sum_{i=0}^{n} BEZ_{i,n}(u) \equiv 1, \quad \forall u\in(0,1) \thinspace。\tag{3}$

我们定义：
$\beta_{i}^{(r)} = \begin{cases} P_i, & r = 0 \\ (1-u)\beta_{i}^{(r-1)} + u\beta_{i+1}^{(r-1)}, & r = 1,2,\cdots,n \end{cases} \thinspace， \tag{4}$
其中 $i=0,1,\cdots,n-r$ 。

可以证明（证明过程参见维基百科De Casteljau算法）：
$\mathcal{P}(u) = \sum_{i=0}^{n} \beta_{i}^{(0)}\cdot BEZ_{i,n}(u) = \beta_{0}^{(n)} \thinspace。 \tag{5}$

于是得到算法逻辑如下：

在实数区间 $[0, 1]$ 上大量采样。对于每个样本 $u$ ，计算 $\mathcal{P}(u)$ 。
根据公式(5)，我们可以通过计算 $\beta_{0}^{(n)}$ 来计算 $\mathcal{P}(u)$ 。
$\beta_{0}^{(n)}$ 可以通过公式(4)进行递归计算。

代码实现

def draw_curve(p_list):
	"""
    :param p_list: (list of list of int: [[x0, y0], [x1, y1], [x2, y2], ...]) 曲线的控制点坐标列表
    :return: (list of list of int: [[x_0, y_0], [x_1, y_1], [x_2, y_2], ...]) 绘制结果的像素点坐标列表
    """
    result = []
    n = len(p_list) - 1
    freq = 1000
    for u in range(freq + 1):
        u /= freq
        tmp = [[p[0], p[1]] for p in p_list]
        for r in range(1, n + 1):
            tmp = [
                [(1 - u) * tmp[i][0] + u * tmp[i + 1][0],
                 (1 - u) * tmp[i][1] + u * tmp[i + 1][1]
                 ] for i in range(n - r + 1)
            ]
        result.append([round(tmp[0][0]), round(tmp[0][1])])
    return result