Bezier曲线

已于 2024-03-31 13:45:09 修改
阅读量359 收藏 9
点赞数 4
分类专栏： 计算机图形学 文章标签： python 开发语言
于 2024-03-31 13:33:08 首次发布
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/qq_55552561/article/details/137047482
版权

1. 实验要求

2. Bezier曲线的原理 以及 公式推导

参考贝塞尔曲线（Bezier Curve）原理及公式推导_bezier曲线-CSDN博客

Bezier曲线的一些特性：

使用n个控制点\left \{ P_{1},P_{2},P_{3}......P_{n} \right \}来控制曲线形状

曲线通过起始点P_{1}和终止点P_{n}，接近但不通过中间点 P_{2}-P_{n-1}

2.1 直观理解

Step 1. 在二维平面内选三个不同的点并依次用线段连接

Step 2. 在线段AB 和BC 上找到 D、E 两点，使得 \frac{AD}{DB}=\frac{BE}{EC}

Step 3. 连接DE ，并在 DE 上找到 F 点，使其满足 \frac{DF}{FE}=\frac{AD}{DB}=\frac{BE}{EC}    抛物线的三切定理

Step 4. 找出符合上述条件的所有点

2.2 bezier曲线公式

一次贝塞尔曲线：B(t) = (1-t)P_{0} + tP_{1} ,t\epsilon [0,1]

二次贝塞尔曲线：B(t) = (1-t)^{2} P_{0} + 2t(1-t) P_{1} + t^{2} P_{2}, t\epsilon [0,1]

三次贝塞尔曲线：B(t) = (1-t)^{3} P_{0} + 3t(1-t)^{2} P_{1} + 3t^{2}(1-t) P_{2}+t^{3}P_{3}, t\epsilon [0,1]

3. 简单的三次贝塞尔曲线的绘制

使用PyOpenGL（环境的配置可以查看PyOpenGL环境配置教程-CSDN博客

from OpenGL.GL import *
from OpenGL.GLUT import *
import numpy as np

# 控制点
control_points = [(-0.5, -0.5), (-0.6, 0.25), (0.25, 0.25), (0.75, -0.5)]


# 计算Bezier曲线上的点
def bezier_curve(control_points, t):
    n = len(control_points) - 1
    result = [0, 0]
    for i in range(n + 1):
        result[0] += control_points[i][0] * binomial_coefficient(n, i) * ((1 - t) ** (n - i)) * (t ** i)
        result[1] += control_points[i][1] * binomial_coefficient(n, i) * ((1 - t) ** (n - i)) * (t ** i)
    return result


# 计算二项式系数
def binomial_coefficient(n, k):
    result = 1
    for i in range(1, k + 1):
        result *= (n - i + 1) / i
    return result


# 绘制Bezier曲线
def draw_bezier_curve():
    glBegin(GL_LINE_STRIP)
    for t in np.arange(0, 1.01, 0.01):
        point = bezier_curve(control_points, t)
        glVertex2fv(point)
    glEnd()

    # 绘制控制点
    glColor3f(1.0, 0.0, 0.0)  # 设置控制点颜色为红色
    glPointSize(5.0)  # 设置控制点大小
    glBegin(GL_POINTS)
    for point in control_points:
        glVertex2fv(point)
    glEnd()

    # 将相邻的控制点连线
    glColor3f(0.0, 1.0, 0.0)  # 设置连线颜色为绿色
    glBegin(GL_LINES)
    for i in range(len(control_points) - 1):
        glVertex2fv(control_points[i])
        glVertex2fv(control_points[i + 1])
    glEnd()


def display():
    glClear(GL_COLOR_BUFFER_BIT)
    glColor3f(1.0, 1.0, 1.0)
    draw_bezier_curve()
    glFlush()


def main():
    glutInit()
    glutInitDisplayMode(GLUT_SINGLE | GLUT_RGB)
    glutInitWindowSize(400, 400)
    glutCreateWindow(b"Bezier Curve")

    glClearColor(0.0, 0.0, 0.0, 0.0)
    glMatrixMode(GL_PROJECTION)
    glLoadIdentity()
    glOrtho(-1.0, 1.0, -1.0, 1.0, -1.0, 1.0)
    glMatrixMode(GL_MODELVIEW)

    glutDisplayFunc(display)

    glutMainLoop()


if __name__ == "__main__":
    main()

结果

修改控制点

control_points = [(-0.5, -0.5), (0, 0.5), (0.25, -0.25), (0.5, 0.25), (0.75, 0),(0.6,-0.1)]

得到5阶贝塞尔曲线

 

 4.增加交互功能

4.1 实现鼠标可以拖动控制点

添加了两个函数，鼠标的控制

from OpenGL.GL import *
from OpenGL.GLUT import *
import numpy as np

# 控制点
control_points = [(-0.5, -0.5), (0, 0.5), (0.25, -0.25), (0.5, 0.25), (0.75, 0),(0.6,-0.1)]

# 选中的控制点索引
selected_point_index = -1

# 计算Bezier曲线上的点
def bezier_curve(control_points, t):
    n = len(control_points) - 1
    result = [0, 0]
    for i in range(n + 1):
        result[0] += control_points[i][0] * binomial_coefficient(n, i) * ((1 - t) ** (n - i)) * (t ** i)
        result[1] += control_points[i][1] * binomial_coefficient(n, i) * ((1 - t) ** (n - i)) * (t ** i)
    return result

# 计算二项式系数
def binomial_coefficient(n, k):
    result = 1
    for i in range(1, k + 1):
        result *= (n - i + 1) / i
    return result

# 鼠标点击事件处理函数
def mouse(button, state, x, y):
    global selected_point_index
    if button == GLUT_LEFT_BUTTON and state == GLUT_DOWN:
        # 将鼠标点击位置的窗口坐标转换为归一化坐标
        x_normalized = 2 * (x / glutGet(GLUT_WINDOW_WIDTH)) - 1
        y_normalized = 1 - 2 * (y / glutGet(GLUT_WINDOW_HEIGHT))

        # 寻找离鼠标点击位置最近的控制点
        min_distance = float('inf')
        for i, point in enumerate(control_points):
            distance = np.sqrt((x_normalized - point[0]) ** 2 + (y_normalized - point[1]) ** 2)
            if distance < min_distance:
                min_distance = distance
                selected_point_index = i



# 鼠标移动事件处理函数
def motion(x, y):
    global selected_point_index
    if selected_point_index != -1:
        # 将鼠标位置的窗口坐标转换为归一化坐标
        x_normalized = 2 * (x / glutGet(GLUT_WINDOW_WIDTH)) - 1
        y_normalized = 1 - 2 * (y / glutGet(GLUT_WINDOW_HEIGHT))

        # 更新选中控制点的坐标
        control_points[selected_point_index] = (x_normalized, y_normalized)

    glutPostRedisplay()  # 通知OpenGL窗口需要重新绘制


# 绘制函数
def draw():
    glClear(GL_COLOR_BUFFER_BIT)
    glColor3f(1.0, 1.0, 1.0)

    # 绘制Bezier曲线
    glBegin(GL_LINE_STRIP)
    for t in np.arange(0, 1.01, 0.01):
        point = bezier_curve(control_points, t)
        glVertex2fv(point)
    glEnd()

    # 绘制控制点
    glColor3f(1.0, 0.0, 0.0)
    glPointSize(5.0)
    glBegin(GL_POINTS)
    for point in control_points:
        glVertex2fv(point)
    glEnd()

    # 连接相邻的控制点
    glColor3f(0.0, 1.0, 0.0)
    glBegin(GL_LINES)
    for i in range(len(control_points) - 1):
        glVertex2fv(control_points[i])
        glVertex2fv(control_points[i + 1])
    glEnd()

    glFlush()

# 主函数
def main():
    glutInit()
    glutInitDisplayMode(GLUT_SINGLE | GLUT_RGB)
    glutInitWindowSize(400, 400)
    glutCreateWindow(b"Bezier Curve")

    glClearColor(0.0, 0.0, 0.0, 0.0)
    glMatrixMode(GL_PROJECTION)
    glLoadIdentity()
    glOrtho(-1.0, 1.0, -1.0, 1.0, -1.0, 1.0)
    glMatrixMode(GL_MODELVIEW)

    glutMouseFunc(mouse)
    glutMotionFunc(motion)
    glutDisplayFunc(draw)

    glutMainLoop()

if __name__ == "__main__":
    main()

结果5阶bezier鼠标拖动演示-CSDN直播

5阶bezier鼠标拖动演示

4.2 采用鼠标添加控制点
from OpenGL.GL import *
from OpenGL.GLUT import *
import numpy as np

# 控制点
control_points = []

# 是否允许添加新的控制点
allow_adding_points = True

# 是否处于控制点移动模式
is_move_mode = False

# 选中的控制点索引
selected_point_index = -1


# 计算Bezier曲线上的点
def bezier_curve(control_points, t):
    n = len(control_points) - 1
    result = [0, 0]
    for i in range(n + 1):
        result[0] += control_points[i][0] * binomial_coefficient(n, i) * ((1 - t) ** (n - i)) * (t ** i)
        result[1] += control_points[i][1] * binomial_coefficient(n, i) * ((1 - t) ** (n - i)) * (t ** i)
    return result


# 计算二项式系数
def binomial_coefficient(n, k):
    result = 1
    for i in range(1, k + 1):
        result *= (n - i + 1) / i
    return result


# 鼠标点击事件处理函数
# 鼠标点击事件处理函数
def mouse(button, state, x, y):
    global allow_adding_points, selected_point_index, is_move_mode
    if allow_adding_points and button == GLUT_LEFT_BUTTON and state == GLUT_DOWN:
        # 将鼠标点击位置的窗口坐标转换为归一化坐标
        x_normalized = 2 * (x / glutGet(GLUT_WINDOW_WIDTH)) - 1
        y_normalized = 1 - 2 * (y / glutGet(GLUT_WINDOW_HEIGHT))

        # 添加新的控制点
        control_points.append((x_normalized, y_normalized))

        # 通知OpenGL窗口需要重新绘制
        glutPostRedisplay()
    elif button == GLUT_LEFT_BUTTON and state == GLUT_DOWN:
        # 寻找离鼠标点击位置最近的控制点
        min_distance = float('inf')
        for i, point in enumerate(control_points):
            distance = np.sqrt((2 * x / glutGet(GLUT_WINDOW_WIDTH) - 1 - point[0]) ** 2 + (
                        1 - 2 * y / glutGet(GLUT_WINDOW_HEIGHT) - point[1]) ** 2)
            if distance < min_distance:
                min_distance = distance
                selected_point_index = i
        if min_distance < 0.05:  # 如果鼠标点击到了控制点附近，则进入控制点移动模式
            is_move_mode = True
    elif button == GLUT_LEFT_BUTTON and state == GLUT_UP:
        is_move_mode = False
        selected_point_index = -1


# 鼠标移动事件处理函数
def motion(x, y):
    global selected_point_index, is_move_mode
    if is_move_mode:
        # 将鼠标位置的窗口坐标转换为归一化坐标
        x_normalized = 2 * (x / glutGet(GLUT_WINDOW_WIDTH)) - 1
        y_normalized = 1 - 2 * (y / glutGet(GLUT_WINDOW_HEIGHT))

        # 更新选中控制点的坐标
        control_points[selected_point_index] = (x_normalized, y_normalized)

        # 通知OpenGL窗口需要重新绘制
        glutPostRedisplay()


# 键盘事件处理函数
def keyboard(key, x, y):
    global allow_adding_points
    if key == b'\r':  # 按下回车键
        allow_adding_points = False
    elif key == b'\x1b':  # 按下ESC键
        allow_adding_points = True


# 绘制函数
def draw():
    glClear(GL_COLOR_BUFFER_BIT)
    glColor3f(1.0, 1.0, 1.0)

    # 绘制Bezier曲线
    glBegin(GL_LINE_STRIP)
    for t in np.arange(0, 1.01, 0.01):
        point = bezier_curve(control_points, t)
        glVertex2fv(point)
    glEnd()

    # 绘制控制点
    glColor3f(1.0, 0.0, 0.0)
    glPointSize(5.0)
    glBegin(GL_POINTS)
    for point in control_points:
        glVertex2fv(point)
    glEnd()

    # 连接相邻的控制点
    glColor3f(0.0, 1.0, 0.0)
    glBegin(GL_LINES)
    for i in range(len(control_points) - 1):
        glVertex2fv(control_points[i])
        glVertex2fv(control_points[i + 1])
    glEnd()

    glFlush()


# 主函数
def main():
    glutInit()
    glutInitDisplayMode(GLUT_SINGLE | GLUT_RGB)
    glutInitWindowSize(400, 400)
    glutCreateWindow(b"Bezier Curve")

    glClearColor(0.0, 0.0, 0.0, 0.0)  # 将窗体的背景色设置为黑色

    glMatrixMode(GL_PROJECTION)
    glLoadIdentity()
    glOrtho(-1.0, 1.0, -1.0, 1.0, -1.0, 1.0)
    glMatrixMode(GL_MODELVIEW)

    glutMouseFunc(mouse)
    glutMotionFunc(motion)
    glutDisplayFunc(draw)
    glutKeyboardFunc(keyboard)

    glutMainLoop()


if __name__ == "__main__":
    main()

结果

bezier曲线（鼠标添加控制点+鼠标拖动控制点

全是头发的羊羊羊
关注
  • 4
    点赞
  • 踩
  • 9
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
n次Bezier曲线
04-02
在vc上基于openGL编写的n次Bezier曲线，可以选择n的大小绘制指定次数的Bezier曲线，之后有较好的交互性，可以在n次Bezier曲线的基础上，将两点拖至重合，变成n-1次Bezier曲线
Bezier曲线的构建
sinat_19682279的博客
04-06 3686
Bezier曲线是由法国雷诺汽车公司的P.E.Bezier于20世纪70年代初为解决汽车外型设计而提出的一种新的参数表示法，这种方法的特点是：控制点的输入与曲线输出之间的关系明确，使设计人员比较直观地估计给定条件与设计出的曲线之间的关系。当设计人员（用户）使用交互手段改变输入控制点，就能很方便地在屏幕上改变拟合曲线的形状与代表它的多项式的次数以迎合设计要求。Bezier曲线是指用光滑参数曲线段逼近...
Bezier 曲线 2D
最新发布
言午日月
12-19 933
控制多边形的第一个顶点和最后一个顶点位于曲线上，多边形的第一条边和最后一条边表示了曲线在起点和终点的切矢量方向，其顶点则用于定义曲线的导数，阶次和形状。在工程项目中，最常见的是三次Bezier 曲线，其次数是二次Bezier 曲线，高次Bezier 曲线一般很少用。二次 Bezier 曲线，当 n=2 时，Bezier 曲线的控制多边形有 3 个控制点，三次 Bezier 曲线，当 n=3 时，Bezier 曲线控制多边形有 4 个控制点，， Bezier 曲线是二次多项式，称为二次 Bezier 曲线，
Bezier曲线及其性质
热门推荐
NEO's NOTEBOOK
08-10 1万+
本文为Bezier曲线mooc教程学习笔记。 目录Bezier曲线与曲面[1] Bezier曲线的背景[2] Bezier曲线[3] Bezier曲线详细定义[4] Bezier曲线举例[5] Bernstein 基函数性质[6] 贝塞尔曲线的性质[7] Bezier曲线的生成[8] Bezier曲线的拼接[9] Bezier曲线的升阶与降阶 Bezier曲线与曲面 [1] Bezier曲线的背景 如果要求曲线通过所有的数据点，则属于插值问题； 如果只要求曲线逼近这些数据点，则属于逼近.
Bezier曲线构造
hanmingjunv5的博客
07-10 3349
Bezier曲线构造 曲线公式 曲线：C(u)=∑i=0nBn,i(u)Pi基函数：Bn,i=n!i!(n−i)!ui(1−u)n−i 曲线：C(u) = \sum^n_{i=0}B_{n,i}(u)P_i\\ 基函数：B_{n,i}=\frac{n!}{i!(n-i)!}u^i(1-u)^{n-i} 曲线：C(u)=i=0∑n​Bn,i​(u)Pi​基函数：Bn,i​=i!(n−i)!n!​ui(1−u)n−i 曲线性质 1.Bezir曲线的阶数由控制点个数决定，n+1个控制点，则对应阶数为n. 例如3次
bezier曲线
04-12
bezier画曲线 opengl 计算机图形学作业 四个点构成一个曲线
OpenGL绘制Bezier曲线的方法
08-19
主要为大家详细介绍了OpenGL绘制Bezier曲线的方法，文中示例代码介绍的非常详细，具有一定的参考价值，感兴趣的小伙伴们可以参考一下
OpenGL画bezier曲线
01-20
Bezier Curve算法是根据参数曲线方程来得到光滑曲线的一种算法，曲线方程的参数由控制点决定。...如下图，t = AE:AB = BF:BC = CG:CD = EH:EF = FI:FG = HJ:HI，J即为Bezier曲线上的点   t取0到1的过程就把be
de Casteljau算法绘制Bezier曲线
03-22
利用de Casteljau算法绘制Bezier曲线，是利用了递归的思想
绘制三次 Bezier 曲线1
08-08
河北经贸大学计算机图形学课程设计课题名称：生成三次Bezier曲线学 号： 201672450124 姓 名： 刘振宇 软件开发环境操作系统：windows 1
二次Bezier曲线源代码
11-30
计算机图形学孔令德版第二版习题7.1 VC2010编译器编写，实现鼠标左键拖动控制多边形顶点，同时曲线随之发生变化
Matlab_画二次及三次Bezier曲线.rar_bezier利用MATLAB_三次曲线_二次及三次Bezier曲线
09-24
通过典型例子讲述如何利用matlab画二次及三次Bezier曲线，很经典。
Bezier曲线算法的实现
09-14
Bezier曲线算法的实现。 Bezier曲线的递推(de Casteljau)算法 Bezier曲线上的任一个点p(t)，都是其它相邻线段的同等比例(t)点处的连线，再取同等比例(t)的点再连线，一直取到最后那条线段的同等比例(t)处，该点就是...
Bezier曲线_Bezier曲线_
10-02
工程压缩包，可以直接打开运行，计算机图形学，共享
画bezier曲线代码
05-22
利用传入参数给定的控制顶点构成的顶点序列画出bezier曲线
MFC bezier曲线
05-13
计算机图形学实验 基于c++6.0的MFC类库实现的单界面程序源代码bezier曲线
bezier曲线算法
11-14
Bezier曲线是一种数学曲线，它可以用来创建光滑的曲线。Bezier曲线算法是通过一系列控制点来定义曲线的形状。以下是Bezier曲线算法的基本步骤： 1.定义控制点：定义一组控制点，这些点将决定曲线的形状。 2.计算Bezier曲线上的点：通过递归地计算相邻线段的同等比例点处的连线，再取同等比例的点再连线，一直取到最后那条线段的同等比例处，该点就是Bezier曲线上的点。 3.绘制Bezier曲线：将计算出的Bezier曲线上的点连接起来，就可以绘制出Bezier曲线。 下面是一个使用Python实现的例子： ```python import matplotlib.pyplot as plt import numpy as np def bezier_curve(control_points, num=1000): t = np.linspace(0, 1, num=num) n = len(control_points) result = [] for i in range(num): point = np.zeros((2,)) for j in range(n): point += control_points[j] * bernstein_poly(j, n-1, t[i]) result.append(point) return result def bernstein_poly(i, n, t): return comb(n, i) * t**i * (1-t)**(n-i) def comb(n, i): return np.math.factorial(n) / (np.math.factorial(i) * np.math.factorial(n-i)) control_points = np.array([[0, 0], [1, 2], [3, 1], [4, 3]]) curve_points = bezier_curve(control_points) plt.plot(control_points[:,0], control_points[:,1], 'ro--') plt.plot([p[0] for p in curve_points], [p[1] for p in curve_points], 'b-') plt.show() ``` 该例子中，我们定义了四个控制点，然后使用`bezier_curve`函数计算出Bezier曲线上的点，最后使用Matplotlib库绘制出Bezier曲线。你可以根据自己的需要修改控制点的坐标来创建不同形状的Bezier曲线。

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值