【图形图像的C++ 实现 01/20】 2D 和 3D 贝塞尔曲线

目录

  • 一、说明
  • 二、贝塞尔曲线特征
  • 三、模拟
  • 四、全部代码如下

一、说明

   以下文章介绍了用 C++ 计算和绘制的贝塞尔曲线（2D 和 3D）。
   贝塞尔曲线具有出色的数学能力来计算路径（从起点到目的地点的曲线）。曲线的形状由“控制点”决定。所讨论的曲线最重要的特征是平滑度。
   在许多应用和领域中，平滑度是不可或缺的。我们可以考虑机器人或其他机器的运动，其中运动必须是可预测的，以确保人员和硬件的安全（低磨损系数）。当机器人关节的轨迹被计算为平滑路径时，我们可以假设机器人将按照规划的路径平滑地移动，不会出现急动或意外移动。请注意，在我们考虑的机器人技术中，除了路径之外，还有速度、加速度、冲击力和电机扭矩。所有这些参数主要影响最终路径。
   除了机器人技术之外，贝塞尔曲线还用于动画、游戏和设计。

   为了绘图的目的，我将使用我之前的文章中讨论过的 C++ 的 matplotlib 库。
   头文件（用于绘图库）必须与您的 cpp 位于同一文件夹中。您的程序可以按如下方式编译，

​//compile
g++ my_prog.cpp -o my_prog -I/usr/include/python3.8 -lpython3.8// 
//run
./my_prog
//folder tree
├── my_prog
├── my_prog.cpp
├── matplotlibcpp.h

​

二、贝塞尔曲线特征

   可以计算点集的贝塞尔曲线： { P0, P1, P2 …Pn}，其中n定义我们建模的曲线（多项式）的阶数。在每种情况下，第一个点和最后一个点定义曲线的起点和终点的位置。其他点 - 控制点通常不属于计算的曲线，而是影响贝塞尔曲线的形状。

   2D中的每个点P都有两个{x,y}笛卡尔坐标，但在3D中，点P按预期由三个{x, y, z}定义。

   贝塞尔曲线的显式定义可以指定如下（我们将在模拟中使用这个公式）。

在这里插入图片描述
这里
在这里插入图片描述

   是二项式系数。

   在我们的例子中，二项式系数的计算如下（如果您查看维基百科，您会发现递归实现，但这是最简单的版本或更直观）。

   C++ 中的实现可以如下所示，

​double computeBinominal(int n, int k)
{
    double value = 1.0;
    for (int i = 1; i <= k; i++)
    {
        value = value * ((n + 1 - i) / i);
    }
    if (n == k){
        value = 1;
    }
    return value;
}

​
平面空间中的四个点P 0 、P 1 、P 2 和P 3 定义三次贝塞尔曲线。该曲线可以建模为三阶多项式。
在这里插入图片描述

当提供六个点P 0、P 1、P 2、P 3、P4和P5时，贝塞尔曲线被计算为五阶多项式。

在这里插入图片描述

三、模拟

   现在我们将显示上面定义的曲线的 2D 和 3D 模拟（针对 4 点和 6 点）。下面的代码为您提供了计算和绘制您想要的任何数字点P 的贝塞尔曲线的绝佳机会。

x{2.5, 1.5, 6.0, 10.0}; 
y{0.5, 5.0, 5.0, 0.5};

x{2.5, 1.5, 6.0, 10.0}; 
//与 2D y{0.5, 5.0, 5.0, 0.5}相同；
//与 2D z{1.0, 2.0, 3.0, 4.0}相同；

X{2.5, 1.5, 6, 10, 7, 3}; 
Y{0.5, 5.0, 5.0, 0.5, 1.0, 2.0};

X{2.5, 1.5, 6.0, 10.0, 7.0, 3.0}; // 对于 2D 
Y{0.5, 5.0, 5.0, 0.5, 1.0 , 2.0}; // 对于 2D 
Z{1.0, 2.0, 3.0, 4.0, 5.0, 0.1};

在这里插入图片描述
对于相同阶的多项式（三阶），我们可以计算 3D 贝塞尔曲线。

x{2.5, 1.5, 6.0, 10.0}; //same as 2D
y{0.5, 5.0, 5.0, 0.5}; //same as 2D
z{1.0, 2.0, 3.0, 4.0};

在这里插入图片描述
这是一条 2D 贝塞尔曲线，它是针对五阶多项式（六点）计算的。

X{2.5, 1.5, 6, 10, 7, 3};
Y{0.5, 5.0, 5.0, 0.5, 1.0 , 2.0};

在这里插入图片描述
和以前一样，我们可以绘制 3D 贝塞尔曲线。

X{2.5, 1.5, 6.0, 10.0, 7.0, 3.0}; //as for 2D
Y{0.5, 5.0, 5.0, 0.5, 1.0 , 2.0}; //as for 2D
Z{1.0, 2.0, 3.0, 4.0, 5.0, 0.1};

在这里插入图片描述

四、全部代码如下

/// g++ bezier_curve.cpp -o t -I/usr/include/python3.8 -lpython3.8

#include <iostream>
#include <vector>
#include <tuple>
#include <math.h>

#include "matplotlibcpp.h"

namespace plt = matplotlibcpp;

//-----------------------------------------------------------

std::tuple<std::vector<double>, std::vector<double>> computeBesierCurve2D(std::vector<double> xX, std::vector<double> yY)
{

    std::vector<double> bCurveX;
    std::vector<double> bCurveY;
    double bCurveXt;
    double bCurveYt;

    for (double t = 0.01; t <= 1; t += 0.01)
    {

        bCurveXt = std::pow((1 - t), 3) * xX[0] + 3 * std::pow((1 - t), 2) * t * xX[1] + 3 * std::pow((1 - t), 1) * std::pow(t, 2) * xX[2] + std::pow(t, 3) * xX[3];
        bCurveYt = std::pow((1 - t), 3) * yY[0] + 3 * std::pow((1 - t), 2) * t * yY[1] + 3 * std::pow((1 - t), 1) * std::pow(t, 2) * yY[2] + std::pow(t, 3) * yY[3];

        bCurveX.push_back(bCurveXt);
        bCurveY.push_back(bCurveYt);
    }

    return std::make_tuple(bCurveX, bCurveY);
}

//-----------------------------------------------------------

void plot2D(std::tuple<std::vector<double>, std::vector<double>> data)
{

    std::vector<double> xX = std::get<0>(data);
    std::vector<double> yY = std::get<1>(data);


    plt::plot(xX, yY);
    plt::show();
}

//-----------------------------------------------------------

std::tuple<std::vector<double>, std::vector<double>, std::vector<double>> computeBesierCurve3D(std::vector<double> xX, std::vector<double> yY, std::vector<double> zZ)
{

    std::vector<double> bCurveX;
    std::vector<double> bCurveY;
    std::vector<double> bCurveZ;
    double bCurveXt;
    double bCurveYt;
    double bCurveZt;

    for (double t = 0.01; t <= 1; t += 0.01)
    {

        bCurveXt = std::pow((1 - t), 3) * xX[0] + 3 * std::pow((1 - t), 2) * t * xX[1] + 3 * std::pow((1 - t), 1) * std::pow(t, 2) * xX[2] + std::pow(t, 3) * xX[3];
        bCurveYt = std::pow((1 - t), 3) * yY[0] + 3 * std::pow((1 - t), 2) * t * yY[1] + 3 * std::pow((1 - t), 1) * std::pow(t, 2) * yY[2] + std::pow(t, 3) * yY[3];
        bCurveZt = std::pow((1 - t), 3) * yY[0] + 3 * std::pow((1 - t), 2) * t * yY[1] + 3 * std::pow((1 - t), 1) * std::pow(t, 2) * yY[2] + std::pow(t, 3) * yY[3];

        bCurveX.push_back(bCurveXt);
        bCurveY.push_back(bCurveYt);
        bCurveZ.push_back(bCurveZt);
    }

    return std::make_tuple(bCurveX, bCurveY, bCurveZ);
}

//-----------------------------------------------------------

void plot3Dexample()
{

    std::vector<double> xX;
    std::vector<double> yY;
    std::vector<double> zZ;
    double theta;
    double r;
    double z_inc = 4.0 / 99.0;
    double theta_inc = (8.0 * M_PI) / 99.0;

    for (double i = 0; i < 100; i += 1)
    {
        theta = -4.0 * M_PI + theta_inc * i;
        zZ.push_back(-2.0 + z_inc * i);
        r = zZ[i] * zZ[i] + 1;
        xX.push_back(r * std::sin(theta));
        yY.push_back(r * std::cos(theta));
    }

    plt::plot3(xX, yY, zZ);
    plt::show();
}

//-----------------------------------------------------------

void plot3D(std::tuple<std::vector<double>, std::vector<double>, std::vector<double>> data)
{

    std::vector<double> xX = std::get<0>(data);
    std::vector<double> yY = std::get<1>(data);
    std::vector<double> zZ = std::get<2>(data);

    plt::plot3(xX, yY, zZ);
    plt::xlabel("x");
    plt::ylabel("y");
    plt::set_zlabel("z");
    plt::show();
}

//-----------------------------------------------------------

double computeBinominal(int n, int k)
{

    double value = 1.0;

    for (int i = 1; i <= k; i++)
    {

        value = value * ((n + 1 - i) / i);
    }

    if (n == k){
        value = 1;
    }
    
    return value;
}

//-----------------------------------------------------------

std::tuple<std::vector<double>, std::vector<double>> computeNVertexBasierCurve2D(std::vector<double> xX, std::vector<double> yY)
{

    std::vector<double> bCurveX;
    std::vector<double> bCurveY;

    int n = xX.size() - 1;
    std::cout << "n :" << n << "\n";

    for (double t = 0.0; t <= 1.0; t += 0.01)
    {

        double bCurveXt{0};
        double bCurveYt{0};

        for (int i = 0; i <= n; ++i)
        {

            bCurveXt += computeBinominal(n, i) * std::pow((1 - t), (n - i)) * std::pow(t, i) * xX[i];
            bCurveYt += computeBinominal(n, i) * std::pow((1 - t), (n - i)) * std::pow(t, i) * yY[i];
            //std::cout << " t= "<< t<< " i=" << i << " bCurveXt=" << bCurveXt << " = " << computeBinominal(n, i)  << " * " << std::pow((1 - t), (n - i))  << " * " << std::pow(t, i) << " * " << xX[i] << std::endl;
        }

        bCurveX.push_back(bCurveXt);
        bCurveY.push_back(bCurveYt);
    }

    return std::make_tuple(bCurveX, bCurveY);
}

std::tuple<std::vector<double>, std::vector<double>, std::vector<double>> computeNVertexBasierCurve3D(std::vector<double> xX, std::vector<double> yY, std::vector<double> zZ)
{

    std::vector<double> bCurveX;
    std::vector<double> bCurveY;
    std::vector<double> bCurveZ;

    int n = xX.size() - 1;
    std::cout << "n :" << n << "\n";

    for (double t = 0.0; t <= 1.0; t += 0.01)
    {

        double bCurveXt{0};
        double bCurveYt{0};
        double bCurveZt{0};

        for (int i = 0; i <= n; ++i)
        {

            bCurveXt += computeBinominal(n, i) * std::pow((1 - t), (n - i)) * std::pow(t, i) * xX[i];
            bCurveYt += computeBinominal(n, i) * std::pow((1 - t), (n - i)) * std::pow(t, i) * yY[i];
            bCurveZt += computeBinominal(n, i) * std::pow((1 - t), (n - i)) * std::pow(t, i) * zZ[i];
            //std::cout << " t= "<< t<< " i=" << i << " bCurveXt=" << bCurveXt << " = " << computeBinominal(n, i)  << " * " << std::pow((1 - t), (n - i))  << " * " << std::pow(t, i) << " * " << xX[i] << std::endl;
        }

        bCurveX.push_back(bCurveXt);
        bCurveY.push_back(bCurveYt);
        bCurveZ.push_back(bCurveZt);
    }

    return std::make_tuple(bCurveX, bCurveY, bCurveZ);
}


//-----------------------------------------------------------

int main()
{

    std::vector<double> xX{2.5, 1.5, 6, 10};
    std::vector<double> yY{0.5, 5, 5, 0.5};
    std::vector<double> zZ{1.0, 2.0, 3.0, 4.0};

    std::tuple<std::vector<double>, std::vector<double>> bCurve2D = computeBesierCurve2D(xX, yY);
    plot2D(bCurve2D);
    
    std::tuple<std::vector<double>, std::vector<double>, std::vector<double>> bCurve3D = computeBesierCurve3D(xX, yY, zZ);
    plot3D(bCurve3D);
    
    std::vector<double> xXn{2.5, 1.5, 6, 10, 7, 3};
    std::vector<double> yYn{0.5, 5, 5, 0.5, 1.0 , 2.0};
    std::vector<double> zZn{1, 2, 3, 4, 5, 0.1};

    std::tuple<std::vector<double>, std::vector<double>> bCurve2DxN = computeNVertexBasierCurve2D(xXn, yYn);
    plot2D(bCurve2DxN);
    
    std::tuple<std::vector<double>, std::vector<double>, std::vector<double>> bCurve3DxN = computeNVertexBasierCurve3D(xXn, yYn, zZn);
    plot3D(bCurve3DxN);

}

​

​

​

​​

​

本文来自互联网用户投稿，该文观点仅代表作者本人，不代表本站立场。本站仅提供信息存储空间服务，不拥有所有权，不承担相关法律责任。如若转载，请注明出处：http://www.mfbz.cn/a/378777.html

如若内容造成侵权/违法违规/事实不符，请联系我们进行投诉反馈qq邮箱809451989@qq.com，一经查实，立即删除！

相关文章

可达鸭二月月赛——入门赛第四场T1题解

可达鸭二月月赛——入门赛第四场T1题解

姓名 王胤皓 AC 记录 题意 有一个圆桶&#xff0c;底面半径为 r r r &#xff0c;高为 h h h。 问&#xff1a;小可每天都需要喝水 20 20 20 升&#xff0c;请问小可至少需要用这个桶接几杯水呢&#xff1f; 思路 首先求出圆桶能装的水&#xff0c;也就是这个圆桶的体…
阅读更多...
上下固定中间自适应布局

上下固定中间自适应布局

实现上下固定中间自适应布局 1.通过position：absolute实现 定义如下结构 <body> <div class="container"> <div class="top"></div> <div class="center"></div> <div class="bottom"&…
阅读更多...
Unity BuffSystem buff系统

Unity BuffSystem buff系统

Unity BuffSystem buff系统 一、介绍二、buff系统架构三、架构讲解四、框架使用buff数据Json数据以及工具ShowTypeBuffTypeMountTypeBuffOverlapBuffShutDownTypeBuffCalculateType时间和层数这里也不过多说明了如何给生物添加buff 五、总结 一、介绍 现在基本做游戏都会需要些…
阅读更多...
springboot167基于springboot的医院后台管理系统的设计与实现

springboot167基于springboot的医院后台管理系统的设计与实现

简介 【毕设源码推荐 javaweb 项目】基于springbootvue 的 适用于计算机类毕业设计&#xff0c;课程设计参考与学习用途。仅供学习参考&#xff0c; 不得用于商业或者非法用途&#xff0c;否则&#xff0c;一切后果请用户自负。 看运行截图看 第五章 第四章 获取资料方式 **项…
阅读更多...
【JS逆向五】逆向模拟某网站的生成的【jsonKey】值 仅供学习

【JS逆向五】逆向模拟某网站的生成的【jsonKey】值 仅供学习

逆向日期&#xff1a;2024.02.07 使用工具&#xff1a;Node.js 加密方法&#xff1a;AES 文章全程已做去敏处理&#xff01;&#xff01;&#xff01; 【需要做的可联系我】 可使用AES进行解密处理&#xff08;直接解密即可&#xff09;&#xff1a;在线AES加解密工具 1、打开…
阅读更多...
作业2.8

作业2.8

1、选择题 1.1、以下选项中,不能作为合法常量的是 ____B______ A&#xff09;1.234e04 B&#xff09;1.234e0.4 C&#xff09;1.234e4 D&#xff09;1.234e0 1.2、以下定义变量并初始化错误的是_____D________。 A) char c1 ‘H’ &#xff1b; B) char c…
阅读更多...
MPLS VPN功能组件

MPLS VPN功能组件

VPN实例 VPN实例即为VPN路由转发表VRF&#xff0c;不同VPN之间的路由隔离通过VPN实例实现&#xff0c;PE上存在多个路由转发表&#xff0c;包括一个公网路由转发表&#xff0c;以及一个或多个VPN路由转发表。 PE为每个直接相连的Site建立并维护专门的VPN实例&#xff0c;VPN实…
阅读更多...
JavaScript表单：事件和方法提交

JavaScript表单：事件和方法提交

&#x1f9d1;‍&#x1f393; 个人主页&#xff1a;《爱蹦跶的大A阿》 &#x1f525;当前正在更新专栏&#xff1a;《VUE》 、《JavaScript保姆级教程》、《krpano》、《krpano中文文档》 ​ ​ ✨ 前言 表单 是 web 开发中不可或缺的一部分&#xff0c;用于收集用户输入的…
阅读更多...
Bootstrap学习三

Bootstrap学习三

Bootstrap学习三 文章目录 前言四、Bootstrap插件4.1. 插件概览4.1.1. data属性4.1.2. 编程方式的API4.1.3. 避免命名空间冲突4.1.4. 事件 4.2. 模态框4.2.1. 引入4.2.2. 基本结构4.2.3. 基本使用4.2.4. 触发模态框的方法 4.3. 下拉菜单和滚动监听4.3.1. 下拉菜单4.3.2. 滚动监…
阅读更多...
创建个人的AI智能体agent

创建个人的AI智能体agent

我在魔塔空间创建了一个作息时间管理的agent >>>>>点我体验<<<<<< Agent并非ChatGPT升级版&#xff0c;它不仅告诉你“如何做”&#xff0c;更会帮你去做&#xff01; 尽管ChatGPT的知识储备接近“全知全能”&#xff0c;但当你试着以AI助手…
阅读更多...
Stable Diffusion 模型下载：ToonYou（平涂卡通）

Stable Diffusion 模型下载：ToonYou（平涂卡通）

文章目录 模型介绍生成案例案例一案例二案例三案例四案例五案例六案例七案例八案例九案例十 下载地址 模型介绍 ToonYou 是一个平涂风格的卡通模型&#xff0c;它的画风独特、光感强烈、画面表现力强、场景结构完整&#xff0c;一张张图犹如动画电影截图&#xff0c;非常值得推…
阅读更多...
ArcGIS学习（六）地理数据库

ArcGIS学习（六）地理数据库

ArcGIS学习（六）地理数据库 上个任务我们讲了一个非常重要的知识点一一坐标系。这个任务我们带来另外一个很重要的知识点一一地理数据库。 地理数据库的内容相比于坐标系简单很多！ 首先，先让我们来学习下地理数据库的理论。 ArcGIS 中的地理数据库（Geodatabase）是一个用…
阅读更多...
一文讲透Python函数中的形式参数和实际参数

一文讲透Python函数中的形式参数和实际参数

函数参数包括形式参数和实际参数&#xff0c;简称形参和实参。其中形式参数即是在定义函数时函数后面括号中的参数列表&#xff08;parameterlist&#xff09;&#xff0c;比如上一个帖子的示例中的width, length&#xff1b;实际参数则是调用函数时函数后面括号中的参数值&…
阅读更多...
前端架构: 从vue-cli探究脚手架原理

前端架构: 从vue-cli探究脚手架原理

从使用角度理解什么是脚手架 脚手架本质是一个操作系统的客户端 在终端中去执行一个命令&#xff0c;这个命令本身它就是一个客户端我们其实可以把脚手架理解为操作系统的一个客户端通过命令去执行它的时候&#xff0c;这个命令往往是这样的一个构造&#xff0c;如下 比如&…
阅读更多...
Ajax 详解及其使用

Ajax 详解及其使用

Ajax&#xff08;Asynchronous JavaScript and XML&#xff09;是一种在客户端与服务器之间进行异步通信的技术&#xff0c;它允许网页在不重新加载整个页面的情况下&#xff0c;与服务器交换数据并更新部分网页内容。Ajax 的核心是XMLHttpRequest&#xff08;XHR&#xff09;对…
阅读更多...
Go 语言中如何大小端字节序？int 转 byte 是如何进行的？

Go 语言中如何大小端字节序？int 转 byte 是如何进行的？

嗨&#xff0c;大家好&#xff01;我是波罗学。 本文是系列文章 Go 技巧第十五篇&#xff0c;系列文章查看&#xff1a;Go 语言技巧。 我们先看这样一个问题&#xff1a;“Go 语言中&#xff0c;将 byte 转换为 int 时是否涉及字节序&#xff08;endianness&#xff09;&#x…
阅读更多...
Java SE多态

Java SE多态

文章目录 1.多态&#xff1a;1.1.什么是多态&#xff1a;1.2.多态实现条件&#xff1a;1.2.1.重写&#xff1a;1.2.2.向上转型&#xff1a; 1.多态&#xff1a; 1.1.什么是多态&#xff1a; 多态的概念&#xff1a;通俗来说&#xff0c;就是多种形态&#xff0c;具体点就是去…
阅读更多...
软件应用实例分享，电玩计时计费怎么算，佳易王PS5游戏计时器系统程序教程

软件应用实例分享，电玩计时计费怎么算，佳易王PS5游戏计时器系统程序教程

软件应用实例分享&#xff0c;电玩计时计费怎么算&#xff0c;佳易王PS5游戏计时器系统程序教程 一、前言 以下软件教程以 佳易王电玩计时计费管理系统软件V17.9为例说明 软件文件下载可以点击最下方官网卡片——软件下载——试用版软件下载 点击开始计时后&#xff0c;图片…
阅读更多...
课时17：本地变量_命令变量

课时17：本地变量_命令变量

2.2.3 命令变量 学习目标 这一节&#xff0c;我们从 基础知识、简单实践、小结 三个方面来学习。 基础知识 基本格式 定义方式一&#xff1a;变量名命令注意&#xff1a; 是反引号定义方式二&#xff1a;变量名$(命令)执行流程&#xff1a;1、执行 或者 $() 范围内的命令…
阅读更多...
2021年通信工程师初级 实务 真题

2021年通信工程师初级 实务 真题

文章目录 一、第1章 现代通信网概述&#xff0c;通信网的定义。第10章 通信业务&#xff0c;普遍服务原则10.2.4 通信行业的发展趋势&#xff08;六化&#xff09; 二、第2章 传输网SDH帧结构SDH线路保护倒换&#xff0c;“11 保护”和“1:1保护”波长值λc/f&#xff0c;中心频…
阅读更多...
最新文章

Copyright @ 2022~2023

友情链接： 我的编程经验分享 一条长河网 尧图网站开发 尧图建网站