计算机图形学、贝塞尔曲线及绘制方法、反走样问题的解决（附完整代码）

贝塞尔曲线

  • 1. 本次作业实现的函数及简单描述（详细代码见后）
  • 2. 与本次作业有关的基础知识整理
  • 3. 代码描述（详细）
  • 4. 完整代码
  • 5. 参考文献

（本篇为作者学习计算机图形学时根据作业所撰写的笔记， 如有同课程请勿Crtl+c/v ）

1. 本次作业实现的函数及简单描述（详细代码见后）

这一块主要介绍本次主要实现的几个函数及其功能

1.bezier(const std::vector<cv::Point2f> &control_points, cv::Mat &window)

可以绘制不止四个点的贝塞尔函曲线，并且处理了反走样问题
升级版： 利用距离比例，进行像素点的判别，解决贝塞尔曲线的反走样问题

2.recursive_bezier(const std::vector<cv::Point2f> &control_points, float t)

根据已有的点vector和给定的t值， 计算出来给定t值对应到贝塞尔曲线上的坐标点

2. 与本次作业有关的基础知识整理

这一部分主要介绍一些跟贝塞尔曲线有关的前置知识，以及贝塞尔曲线的画法

贝塞尔曲线：

Bézier curve(贝塞尔曲线) 是应用于二维图形应用程序的数学曲线。

贝塞尔曲线基础定义：

只要求一定经过起止点，起止点之间的若干控制点用于控制曲线弯曲的方向，最终形成一条经过起止点的光滑曲线被称为贝塞尔曲线。

贝塞尔曲线是线性插值的结果 ： “选出两点之间的一个点”

P(t) = P0 + (P1−P0) = (1−t)P0 + t**P1, t∈[0,1] t=P0P1/P0P1

根据控制点的个数，贝塞尔曲线被分为一阶贝塞尔曲线（0个控制点）、二阶贝塞尔曲线（1个控制点）、三阶贝塞尔曲线（2个控制点）等等。（本次实验用到的是三阶贝塞尔曲线（2个控制点））

二阶贝塞尔曲线：

img

三阶贝塞尔曲线（本次作业需要实现的曲线）

外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传

img

对于贝塞尔曲线，最重要的点是数据点和控制点

数据点： 指一条路径的起始点和终止点。
控制点：控制点决定了一条路径的弯曲轨迹

贝塞尔曲线特点：

特点一： 曲线通过始点和终点，并与特征多边形首末两边相切于始点和终点，中间点将曲线拉向自己。
特点二： 平面离散点控制曲线的形状，改变一个离散点的坐标，曲线的形状将随之改变（点对曲线具有整体控制性）。
特点三： 曲线落在特征多边形的凸包之内，它比特征多边形更趋于光滑。

特点四： 起始点的方向为响应控制点的切线方向

逐段贝塞尔曲线：

当控制点太多会影响控制点的效果的时候， 我们选择每四个控制点定义一条贝塞尔曲线，然后再将他们连接起来

img

常见应用场景：

计算机辅助设计和计算机辅助制造应用（CAD/CAM），Adobe Illustrator, Photoshop, Inkscape, Gimp等等。还可以应用在一些图形技术中，像矢量图形（SVG），所以说在可视化学习的重要基础知识。

基本绘制方法： 德卡斯特里奥算法

  • 引入参数t (范围为 0~1)
  • 取b0到b1, b1到b2上t位置上的点b0’,b1’
  • 将b0’,b1’ 连接，取b0’到b1’上t位置上的点b0’’
  • 将所有t∈（0，1） 的点都便利一遍相连即可得到贝塞尔曲线
  • 若有n个控制点则将上面步骤进行递归操作直到找到最后bn’

解决反走样问题

程序在进行画线时是以点的形式，如果放大有较为明显的不连续，因此可以采用反走样来使其平滑过渡;

对于每一个点，其必然会包含在一个像素之中，可按照比例和着色进行插值，以求得一个较为平滑的过渡;

如下图所示，该点一定会在下图所示的黄色方框内部，而该点离像素点的最远距离为根号2，按照该点到其他四个点的距离进行插值；
在这里插入图片描述

3. 代码描述（详细）

这一部分主要介绍实现贝塞尔曲线的各个函数的功能，以及注意要点

de Casteljau的算法来计算贝塞尔曲线的点（步骤如下）

  • 计算点的总个数
  • 判断点的个数，只有一个点直接返回，多个点的话进行计算
  • 递归计算多个点t值对应的贝塞尔曲线上的点，直到计算到最后一个点
cv::Point2f recursive_bezier(const std::vector<cv::Point2f> &control_points, float t) 
{// TODO: Implement de Casteljau's algorithmint n = control_points.size();if (n == 1) return control_points[0];std::vector<cv::Point2f> res_control_points;for (int i = 0; i < n - 1; i++) {res_control_points.push_back(cv::Point2f((1 - t) * control_points[i].x + t * control_points[i + 1].x,(1 - t) * control_points[i].y + t * control_points[i + 1].y));}return recursive_bezier(res_control_points, t);
}

基础版贝塞尔： (简单的将所给点画出)

  • 对每个t值循环计算出来对应的点然后进行绘制
void bezier(const std::vector<cv::Point2f> &control_points, cv::Mat &window) 
{// TODO: Iterate through all t = 0 to t = 1 with small steps, and call de Casteljau's // recursive Bezier algorithm.for (double t = 0.0; t <= 1.0; t += 0.001)  // 循环得到点，并记录输出{cv::Point2f point = recursive_bezier(control_points, t); window.at<cv::Vec3b>(point.y, point.x)[1] = 255;}
}

升级版贝塞尔（可以解决曲线的反走样问题）

  • 先找到根据t值计算出来的贝塞尔曲线上的点，其周围四个坐标
  • 计算该点到周围四个像素点的距离，距离越大，颜色应该越浅 {max-distance}
  • 按照distance/sum(distance)重新给四个点的像素赋值，并与原始色素的值进行比较，不能超过255.
  • 根据计算出来的点和对应的像素值进行重新绘制
void bezier(const std::vector<cv::Point2f> &control_points, cv::Mat &window) 
{// TODO: Iterate through all t = 0 to t = 1 with small steps, and call de Casteljau's // recursive Bezier algorithm.for (float t = 0.0; t <= 1.0; t = t + 0.001){cv::Point2f point = recursive_bezier(control_points, t);  //接收每一个t所求得的点//找到该点周围的四个点坐标，先找到一个右上的顶点，然后其他的通过这个点求出cv::Point2f point_1 = cv::Point2f(point.x - std::floor(point.x) < 0.5 ? std::floor(point.x) : std::ceil(point.x),   //右上的点point.y - std::floor(point.y) < 0.5 ? std::floor(point.y) : std::ceil(point.y));cv::Point2f point_2 = cv::Point2f(point_1.x - 1.0, point_1.y);cv::Point2f point_3 = cv::Point2f(point_1.x - 1.0, point_1.y - 1.0);cv::Point2f point_4 = cv::Point2f(point_1.x, point_1.y - 1.0);//用vector容器存储刚刚求得的四个像素点坐标;std::vector<cv::Point2f> distance_dot{ point_1,point_2,point_3,point_4 };float MaxDistance = sqrt(2.0);float SumDistance = 0.0f;float pi_distance = 0.0f;std::vector<float> distance_List = {};window.at<cv::Vec3b>(point.y, point.x)[1] = 255;//Mat.at<存储类型名称>(行，列)[通道]for (int i = 0; i < 4; i++){cv::Point2f point_coordinate(distance_dot[i].x + 0.5, distance_dot[i].y + 0.5);//记录像素中心点;// 距离像素点越近的点 实际上和颜色占比应该越小 因此采用 max-实际距离  的方法进行数据预处理pi_distance = MaxDistance - sqrt(std::pow(point.x - point_coordinate.x, 2) + std::pow(point.y - (point_coordinate.y), 2));distance_List.push_back(pi_distance);SumDistance += pi_distance;       // 计算一个总的距离值}for (int i = 0; i < 4; i++){float d = distance_List[i] / SumDistance;  // d是距离比例系数window.at<cv::Vec3b>(distance_dot[i].y, distance_dot[i].x)[1] = std::min(255.f, window.at<cv::Vec3b>(distance_dot[i].y, distance_dot[i].x)[1] + 255.f * d);//不超过255.f，要进行最小值比较}}
}

4. 完整代码

#include <chrono>
#include <iostream>
#include <opencv2/opencv.hpp>std::vector<cv::Point2f> control_points;void mouse_handler(int event, int x, int y, int flags, void *userdata) 
{if (event == cv::EVENT_LBUTTONDOWN) {std::cout << "Left button of the mouse is clicked - position (" << x << ", "<< y << ")" << '\n';control_points.emplace_back(x, y);}     
}void naive_bezier(const std::vector<cv::Point2f> &points, cv::Mat &window) 
{auto &p_0 = points[0];auto &p_1 = points[1];auto &p_2 = points[2];auto &p_3 = points[3];for (double t = 0.0; t <= 1.0; t += 0.001) {auto point = std::pow(1 - t, 3) * p_0 + 3 * t * std::pow(1 - t, 2) * p_1 +3 * std::pow(t, 2) * (1 - t) * p_2 + std::pow(t, 3) * p_3;window.at<cv::Vec3b>(point.y, point.x)[2] = 255;}
}cv::Point2f recursive_bezier(const std::vector<cv::Point2f> &control_points, float t) 
{// TODO: Implement de Casteljau's algorithmint n = control_points.size();if (n == 1) return control_points[0];std::vector<cv::Point2f> res_control_points;for (int i = 0; i < n - 1; i++) {res_control_points.push_back(cv::Point2f((1 - t) * control_points[i].x + t * control_points[i + 1].x,(1 - t) * control_points[i].y + t * control_points[i + 1].y));}return recursive_bezier(res_control_points, t);
}void bezier(const std::vector<cv::Point2f> &control_points, cv::Mat &window) 
{// TODO: Iterate through all t = 0 to t = 1 with small steps, and call de Casteljau's // recursive Bezier algorithm.for (float t = 0.0; t <= 1.0; t = t + 0.001){cv::Point2f point = recursive_bezier(control_points, t);  //接收每一个t所求得的点//找到该点周围的四个点坐标，先找到一个右上的顶点，然后其他的通过这个点求出cv::Point2f point_1 = cv::Point2f(point.x - std::floor(point.x) < 0.5 ? std::floor(point.x) : std::ceil(point.x),   //右上的点point.y - std::floor(point.y) < 0.5 ? std::floor(point.y) : std::ceil(point.y));cv::Point2f point_2 = cv::Point2f(point_1.x - 1.0, point_1.y);cv::Point2f point_3 = cv::Point2f(point_1.x - 1.0, point_1.y - 1.0);cv::Point2f point_4 = cv::Point2f(point_1.x, point_1.y - 1.0);//用vector容器存储刚刚求得的四个像素点坐标;std::vector<cv::Point2f> distance_dot{ point_1,point_2,point_3,point_4 };float MaxDistance = sqrt(2.0);float SumDistance = 0.0f;float pi_distance = 0.0f;std::vector<float> distance_List = {};window.at<cv::Vec3b>(point.y, point.x)[1] = 255;//Mat.at<存储类型名称>(行，列)[通道]for (int i = 0; i < 4; i++){cv::Point2f point_coordinate(distance_dot[i].x + 0.5, distance_dot[i].y + 0.5);//记录像素中心点;// 距离像素点越近的点 实际上和颜色占比应该越小 因此采用 max-实际距离  的方法进行数据预处理pi_distance = MaxDistance - sqrt(std::pow(point.x - point_coordinate.x, 2) + std::pow(point.y - (point_coordinate.y), 2));distance_List.push_back(pi_distance);SumDistance += pi_distance;       // 计算一个总的距离值}for (int i = 0; i < 4; i++){float d = distance_List[i] / SumDistance;  // d是距离比例系数window.at<cv::Vec3b>(distance_dot[i].y, distance_dot[i].x)[1] = std::min(255.f, window.at<cv::Vec3b>(distance_dot[i].y, distance_dot[i].x)[1] + 255.f * d);//不超过255.f，要进行最小值比较}}}int main() 
{cv::Mat window = cv::Mat(700, 700, CV_8UC3, cv::Scalar(0));cv::cvtColor(window, window, cv::COLOR_BGR2RGB);cv::namedWindow("Bezier Curve", cv::WINDOW_AUTOSIZE);cv::setMouseCallback("Bezier Curve", mouse_handler, nullptr);int key = -1;while (key != 27) {window.setTo(0);for (auto &point : control_points) {cv::circle(window, point, 3, {255, 255, 255}, 3);}if (control_points.size() >= 4) {naive_bezier(control_points, window);   // 所有点的绘制的贝塞尔曲线用绿色表示bezier(control_points, window);   // 前四个点绘制的贝塞尔曲线用红色显示cv::imshow("Bezier Curve", window);cv::imwrite("my_bezier_curve.png", window);key = cv::waitKey(1);}cv::imshow("Bezier Curve", window);key = cv::waitKey(20);}return 0;
}

5. 参考文献

1. 贝塞尔曲线

2. opencv 处理鼠标事件的方法

3. 贝塞尔曲线简单介绍

4. 解决反走样问题

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

如若内容造成侵权/违法违规/事实不符，请联系一条长河网进行投诉反馈，一经查实，立即删除！

相关文章

数字时代古文的传承———云南文化瑰宝“爨文化“(我为家乡发声)

数字时代古文的传承———云南文化瑰宝“爨文化“(我为家乡发声)

文章目录 前言⭐ "爨"意味着什么&#xff0c;究竟何为"爨文化"&#xff1f;⭐ 爨文化鲜明的特点1.经济生活2.政治生活3.文化艺术 ⭐ 数字时代古文的传承与传播1.藏品数字化2.建立数据库3.传播大众化 前言 爨文化是继古滇文化之后崛起于珠江正源南盘江流域…
阅读更多...
根据GWAS数据估算样本量N和使用千人基因组填充maf的参考文献

根据GWAS数据估算样本量N和使用千人基因组填充maf的参考文献

https://github.com/GenomicSEM/GenomicSEM/wiki/2.1-Calculating-Sum-of-Effective-Sample-Size-and-Preparing-GWAS-Summary-Statistics
阅读更多...
【SQL】mysql创建定时任务执行存储过程--20230928

【SQL】mysql创建定时任务执行存储过程--20230928

1.先设定时区 https://blog.csdn.net/m0_46629123/article/details/133382375 输入命令show variables like “%time_zone%”;&#xff08;注意分号结尾&#xff09;设置时区&#xff0c;输入 set global time_zone “8:00”; 回车,然后退出重启&#xff08;一定记得重启&am…
阅读更多...
【C++】手撕vector（vector的模拟实现）

【C++】手撕vector（vector的模拟实现）

手撕vector目录&#xff1a; 一、基本实现思路方针 二、vector的构造函数剖析&#xff08;构造歧义拷贝构造&#xff09; 2.1构造函数使用的歧义问题 2.2 vector的拷贝构造和赋值重载&#xff08;赋值重载不是构造哦&#xff0c;为了方便写在一起&#xff09; 三、vector的…
阅读更多...
Apollo自动驾驶系统概述（文末参与活动赠送百度周边）

Apollo自动驾驶系统概述（文末参与活动赠送百度周边）

前言 「作者主页」&#xff1a;雪碧有白泡泡 「个人网站」&#xff1a;雪碧的个人网站 「推荐专栏」&#xff1a; ★java一站式服务 ★ ★ React从入门到精通★ ★前端炫酷代码分享 ★ ★ 从0到英雄&#xff0c;vue成神之路★ ★ uniapp-从构建到提升★ ★ 从0到英雄&#xff…
阅读更多...
算法基础课第二部分

算法基础课第二部分

算法基础课 第四讲 数学知识AcWing1381. 阶乘(同余&#xff0c;因式分解) 质数AcWing 866. 质数的判定---试除法AcWing 868. 质数的判定---埃氏筛AcWing867. 分解质因数---试除法AcWing 197. 阶乘---分解质因数---埃式筛 约数AcWing 869. 求约数---试除法AcWing 870. 约数个数-…
阅读更多...
云可观测性安全平台——掌动智能

云可观测性安全平台——掌动智能

云可观测性安全平台是一个跨架构、跨平台的可观测性方案&#xff0c;实现对云环境下的细粒度数据可视化&#xff0c;满足安全部门对云内部安全领域的多场景诉求&#xff0c;包括敏感数据动态监管、云网攻击回溯分析、攻击横移风险监控、云异常流量分析。本文将介绍掌动智能云可…
阅读更多...
腾讯云中使用ubuntu安装属于自己的overleaf

腾讯云中使用ubuntu安装属于自己的overleaf

在自己的云服务器上安装overleaf的需求是从写论文开始的&#xff0c;总担心自己的论文放在一个网站上被泄露&#xff0c;所以想要在自己的服务器上安装自己的overleaf&#xff0c;正好手边有一个云服务器&#xff0c;现在开始。 配置腾讯云 因为使用overleaf的优势就是在不同…
阅读更多...
SQLAlchemy常用数据类型

SQLAlchemy常用数据类型

目录 SQLAlchemy常用数据类型 代码演示 代码分析 SQLAlchemy常用数据类型 SQLAlchemy 是一个Python的SQL工具库和对象关系映射(ORM)工具&#xff0c;它提供了一种在Python中操作数据库的高效方式。下面是SQLAlchemy中常用的一些数据类型&#xff1a; Integer&#xff1a;整形&…
阅读更多...
【数据结构】选择排序 堆排序（二）

【数据结构】选择排序 堆排序（二）

目录 一&#xff0c;选择排序 1&#xff0c;基本思想 2&#xff0c; 基本思路 3&#xff0c;思路实现 二&#xff0c;堆排序 1&#xff0c;直接选择排序的特性总结&#xff1a; 2&#xff0c;思路实现 3&#xff0c;源代码 最后祝大家国庆快乐&#xff01; 一&#xf…
阅读更多...
数组和切⽚ - Go语言从入门到实战

数组和切⽚ - Go语言从入门到实战

数组和切⽚ - Go语言从入门到实战 数组的声明 package main import "fmt" func main() { var a [3]int //声明并初始化为默认零值 a[0] 1 fmt.Println("a:", a) // 输出: a: [1 0 0] b : [3]int{1, 2, 3} //声明同时初始化 fmt.Println("b:…
阅读更多...
Python3数据科学包系列(二):数据分析实战

Python3数据科学包系列(二):数据分析实战

Python3中类的高级语法及实战 Python3(基础|高级)语法实战(|多线程|多进程|线程池|进程池技术)|多线程安全问题解决方案 Python3数据科学包系列(一):数据分析实战 Python3数据科学包系列(二):数据分析实战 一&#xff1a;通过read_table函数读取数据创建(DataFrame)数据框 #…
阅读更多...
C++指针的使用

C++指针的使用

文章目录 1.C指针1.1 定义指针1.2 使用指针 2.空指针和野指针2.1 空指针2.2 野指针 3.指针所占空间4.使用const修饰指针4.1 const修饰指针4.2 const修饰常量4.3 const 既修饰指针也修饰常量 5.指针操作数组6.指针做函数参数7.使用指针知识实现冒泡排序 1.C指针 指针其实就是一…
阅读更多...
mysql主从同步

mysql主从同步

原理 概述 将主库的数据变更同步到从库&#xff0c;从而保证主库和从库数据一致 数据备份&#xff0c;失败迁移&#xff0c;读写分离&#xff0c;降低单库读写压力 原理 主数据库设置 docker run --restartalways --name mysql-master -p 3306:3306 -v /home/apps/mysql-…
阅读更多...
偏微分方程的人工智能

偏微分方程的人工智能

9 偏微分方程的人工智能 在本节中&#xff0c;我们详细介绍了用于解决偏微分方程&#xff08;Partial Differential Equations&#xff0c;PDEs&#xff09;的人工智能领域的进展。我们在第9.1节中概述了PDE建模的一般形式&#xff0c;并阐述了在这个背景下使用机器学习方法的…
阅读更多...
竞赛 大数据房价预测分析与可视

竞赛 大数据房价预测分析与可视

0 前言 &#x1f525; 优质竞赛项目系列&#xff0c;今天要分享的是 &#x1f6a9; 大数据房价预测分析与可视 &#x1f947;学长这里给一个题目综合评分(每项满分5分) 难度系数&#xff1a;3分工作量&#xff1a;3分创新点&#xff1a;4分 该项目较为新颖&#xff0c;适合…
阅读更多...
JavaScript中如何确定this的值？如何指定this的值？

JavaScript中如何确定this的值？如何指定this的值？

&#x1f380;JavaScript中的this 在绝大多数情况下&#xff0c;函数的调用方法决定了this的值&#xff08;运行时绑定&#xff09;。this不能在执行期间被赋值&#xff0c;并且在每次函数呗调用时this的值也可能会不同。 &#x1f37f;如何确定this的值&#xff1a; 在非严格…
阅读更多...
百度交易中台之内容分润结算系统架构浅析

百度交易中台之内容分润结算系统架构浅析

作者 | 交易中台团队 导读 随着公司内容生态的蓬勃发展&#xff0c;内容产出方和流量提供方最关注的“收益结算”的工作&#xff0c;也就成为重中之重。本文基于内容分润结算业务为入口&#xff0c;介绍了实现过程中的重难点&#xff0c;比如千万级和百万级数据量下的技术选型和…
阅读更多...
C++八股

C++八股

1、简述一下C中的多态 在面向对象中&#xff0c;多态是指通过基类的指针或引用&#xff0c;在运行时动态调用实际绑定对象函数的行为&#xff0c;与之相对应的编译时绑定函数称为静态绑定。 静态多态 静态多态是编译器在编译期间完成的&#xff0c;编译器会根据实参类型来选择…
阅读更多...
UART相关参数和Modbus协议

UART相关参数和Modbus协议

温湿度数据和风速风向数据的读取和计算方法 文章目录 温湿度数据和风速风向数据的读取和计算方法1 串行通信数据格式1.1 协议介绍1.2 UART相关参数1.3 UART通信过程 2 USB转串口模块的使用3 串口调试助手的使用3.1 串口控制区3.2 发送控制区3.3 接收控制区 4 GY-39气象信息模块…
阅读更多...
最新文章

Copyright @ 2022~2023