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

2024-02-09 11:04
文章标签 c++ 实现 3d 01 20 贝塞尔 曲线 2d 图形图像

本文主要是介绍【图形图像的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);}

​五、资源和下载

下面给出源代码资源下载链接地址：
https://download.csdn.net/download/gongdiwudu/88821722

​

​

​​

​

这篇关于【图形图像的C++ 实现 01/20】 2D 和 3D 贝塞尔曲线的文章就介绍到这儿，希望我们推荐的文章对编程师们有所帮助！


http://www.chinasem.cn/article/694041。

相关文章

Spring Boot 实现 IP 限流的原理、实践与利弊解析

Spring Boot 实现 IP 限流的原理、实践与利弊解析

《SpringBoot实现IP限流的原理、实践与利弊解析》在SpringBoot中实现IP限流是一种简单而有效的方式来保障系统的稳定性和可用性,本文给大家介绍SpringBoot实现IP限... 目录一、引言二、IP 限流原理2.1 令牌桶算法2.2 漏桶算法三、使用场景3.1 防止恶意攻击3.2 控制资源
阅读更多...
springboot下载接口限速功能实现

springboot下载接口限速功能实现

《springboot下载接口限速功能实现》通过Redis统计并发数动态调整每个用户带宽,核心逻辑为每秒读取并发送限定数据量,防止单用户占用过多资源,确保整体下载均衡且高效,本文给大家介绍spring... 目录 一、整体目标 二、涉及的主要类/方法✅ 三、核心流程图解（简化） 四、关键代码详解1️⃣ 设置
阅读更多...
Nginx 配置跨域的实现及常见问题解决

Nginx 配置跨域的实现及常见问题解决

《Nginx配置跨域的实现及常见问题解决》本文主要介绍了Nginx配置跨域的实现及常见问题解决,文中通过示例代码介绍的非常详细,对大家的学习或者工作具有一定的参考学习价值,需要的朋友们下面随着小编来... 目录1. 跨域1.1 同源策略1.2 跨域资源共享（CORS）2. Nginx 配置跨域的场景2.1
阅读更多...
Python中提取文件名扩展名的多种方法实现

Python中提取文件名扩展名的多种方法实现

《Python中提取文件名扩展名的多种方法实现》在Python编程中,经常会遇到需要从文件名中提取扩展名的场景,Python提供了多种方法来实现这一功能,不同方法适用于不同的场景和需求,包括os.pa... 目录技术背景实现步骤方法一：使用os.path.splitext方法二：使用pathlib模块方法三
阅读更多...
CSS实现元素撑满剩余空间的五种方法

CSS实现元素撑满剩余空间的五种方法

《CSS实现元素撑满剩余空间的五种方法》在日常开发中，我们经常需要让某个元素占据容器的剩余空间，本文将介绍5种不同的方法来实现这个需求，并分析各种方法的优缺点，感兴趣的朋友一起看看吧... css实现元素撑满剩余空间的5种方法 在日常开发中，我们经常需要让某个元素占据容器的剩余空间。这是一个常见的布局需求
阅读更多...
HTML5 getUserMedia API网页录音实现指南示例小结

HTML5 getUserMedia API网页录音实现指南示例小结

《HTML5getUserMediaAPI网页录音实现指南示例小结》本教程将指导你如何利用这一API，结合WebAudioAPI，实现网页录音功能,从获取音频流到处理和保存录音，整个过程将逐步... 目录1. html5 getUserMedia API简介1.1 API概念与历史1.2 功能与优势1.3
阅读更多...
Java实现删除文件中的指定内容

Java实现删除文件中的指定内容

《Java实现删除文件中的指定内容》在日常开发中,经常需要对文本文件进行批量处理,其中,删除文件中指定内容是最常见的需求之一,下面我们就来看看如何使用java实现删除文件中的指定内容吧... 目录1. 项目背景详细介绍2. 项目需求详细介绍2.1 功能需求2.2 非功能需求3. 相关技术详细介绍3.1 Ja
阅读更多...
使用Python和OpenCV库实现实时颜色识别系统

使用Python和OpenCV库实现实时颜色识别系统

《使用Python和OpenCV库实现实时颜色识别系统》：本文主要介绍使用Python和OpenCV库实现的实时颜色识别系统,这个系统能够通过摄像头捕捉视频流,并在视频中指定区域内识别主要颜色（红... 目录一、引言二、系统概述三、代码解析1. 导入库2. 颜色识别函数3. 主程序循环四、HSV色彩空间详解
阅读更多...
Windows下C++使用SQLitede的操作过程

Windows下C++使用SQLitede的操作过程

《Windows下C++使用SQLitede的操作过程》本文介绍了Windows下C++使用SQLite的安装配置、CppSQLite库封装优势、核心功能（如数据库连接、事务管理）、跨平台支持及性能优... 目录Windows下C++使用SQLite1、安装2、代码示例CppSQLite：C++轻松操作SQ
阅读更多...
PostgreSQL中MVCC 机制的实现

PostgreSQL中MVCC 机制的实现

《PostgreSQL中MVCC机制的实现》本文主要介绍了PostgreSQL中MVCC机制的实现,通过多版本数据存储、快照隔离和事务ID管理实现高并发读写,具有一定的参考价值,感兴趣的可以了解一下... 目录一 MVCC 基本原理python1.1 MVCC 核心概念1.2 与传统锁机制对比二 Postg
阅读更多...

Copyright China编程 23002807@qq.com

沪ICP备2023033072号-2