实现Bezier样条曲线
1.给出n+1 个控制点pk=(xk,yk,zk),这里k可取值0-n,多项式函数公式如下
获取的单个点的代码
void zmBezier::getPoint(float u, double p[3])
{int n = m_count - 1;double x = 0, y = 0, z = 0;for(int k = 0; k <= n; k++){x += m_ctrlPoints[k][0] * BEZ_k_n(n, k, u);y += m_ctrlPoints[k][1] * BEZ_k_n(n, k, u);z += m_ctrlPoints[k][2] * BEZ_k_n(n, k, u);}p[0] = x;p[1] = y;p[2] = z;
}
2.混合函数是如下的多项式
double zmBezier::BEZ_k_n(int n, int k, double u)
{return C_n_k(n, k) * pow(u, k) * pow(1 - u, n - k);
}
3.二项式系数
double zmBezier::C_n_k(int n, int k)
{n = m_count - 1;return factorial(n) / (factorial(k) * factorial(n - k));
}
4.Bezier样条完整代码,全部用指针表示点
/**
Bezier曲线
给定n+1个控制点 Pk=(Xk,Yk,Zk),k取值0-n
多项式函数
-----------------------------------nP(u)= Σ Pk × BEZ(u) 0≤u≤1k=0 k,n
-----------------------------------
混合函数
-----------------------------------k n-kBEZ(u)=C(n,k) × u × (1-u) 0≤u≤1k,n-----------------------------------
二项式系数
-----------------------------------n!C(n,k)=——————————————————k! × (n-k)!-----------------------------------不想使用 点 结构,全部用指针数组表示点集
*/
#ifndef ZMBEZIER_H
#define ZMBEZIER_Hclass zmBezier
{
public:zmBezier();~zmBezier();zmBezier(int n, double (*points)[3]);void getPoint(float u, double p[3]); //获取参数u时的某一点void getCurve(int n, double (*curve)[3]); //获取n个插值点,代表曲线void setCtrlPoints(int n, double (*points)[3]); //设置控制点void getCtrlPoints(int &n, double (*points)[3]); //获取控制点private:inline double factorial(double n); //阶乘inline double C_n_k(int n, int k); //二项式系数,参数n为了形式上更接近二项式inline double BEZ_k_n(int n, int k, double u); //混合函数private:int m_count; //控制点数量double (*m_ctrlPoints)[3]; //控制点坐标
};#endif // ZMBEZIER_H
#include "zmBezier.h"#include<cmath>
#include<string>zmBezier::zmBezier(){m_count = 0;m_ctrlPoints = nullptr;
}zmBezier::zmBezier(int n, double(*points)[3])
{m_count = n;m_ctrlPoints = new double[n][3];memcpy_s(m_ctrlPoints, sizeof (double)*n * 3, points, sizeof (double)*n * 3);
}zmBezier::~zmBezier()
{m_count = 0;delete [] m_ctrlPoints;
}double zmBezier::C_n_k(int n, int k)
{n = m_count - 1;return factorial(n) / (factorial(k) * factorial(n - k));
}double zmBezier::factorial(double n)
{return tgamma(n + 1);
}void zmBezier::getPoint(float u, double p[3])
{int n = m_count - 1;double x = 0, y = 0, z = 0;for(int k = 0; k <= n; k++){x += m_ctrlPoints[k][0] * BEZ_k_n(n, k, u);y += m_ctrlPoints[k][1] * BEZ_k_n(n, k, u);z += m_ctrlPoints[k][2] * BEZ_k_n(n, k, u);}p[0] = x;p[1] = y;p[2] = z;
}double zmBezier::BEZ_k_n(int n, int k, double u)
{return C_n_k(n, k) * pow(u, k) * pow(1 - u, n - k);
}void zmBezier::getCurve(int count, double (*curve)[3])
{double point[3] = {0};for(int k = 0; k < count; k++) {getPoint(1.0 * k / (count - 1), point);curve[k][0] = point[0];curve[k][1] = point[1];curve[k][2] = point[2];}
}void zmBezier::setCtrlPoints(int n, double(*points)[3])
{delete [] m_ctrlPoints;m_count = n;m_ctrlPoints = new double[n][3];int size = sizeof (double) * n * 3;memcpy_s(m_ctrlPoints, size, points, size);
}void zmBezier::getCtrlPoints(int &n, double (*points)[3])
{n = m_count;if(m_count){int size = sizeof (double) * n * 3;memcpy_s(points, size, m_ctrlPoints, size);}
}
5. 继承QWidget,定义可显示的控制点
#ifndef MYCTRLPOINT_H
#define MYCTRLPOINT_H#include <QWidget>class myCtrlPoint : public QWidget
{Q_OBJECT
public:myCtrlPoint(QWidget *parent);QPoint getPosition();void setPostion(const QPoint &point);
protected:void paintEvent(QPaintEvent *event) override;void mousePressEvent(QMouseEvent *event) override;void mouseMoveEvent(QMouseEvent *event) override;void mouseReleaseEvent(QMouseEvent *event) override;private:QPoint m_clicked;};#endif // MYCTRLPOINT_H
主要是实现鼠标事件:
5.1 鼠标左键单击,点变成绿色
5.2 鼠标左键拖动,点在父窗口中移动
5.3 鼠标右键,从父类中删除自己
#include"myCanvas.h"
#include"myCtrlPoint.h"#include<QKeyEvent>
#include<QPainter>
#include<QMouseEvent>myCtrlPoint::myCtrlPoint(QWidget *parent): QWidget(parent)
{setFixedSize(20, 20);
}void myCtrlPoint::paintEvent(QPaintEvent *event)
{QWidget::paintEvent(event);QPainter painter(this);if(m_clicked != QPoint()) {painter.setBrush(Qt::green);}else {painter.setBrush(Qt::lightGray);}painter.drawRect(rect());
}void myCtrlPoint::mousePressEvent(QMouseEvent *event)
{if(event->button() == Qt::LeftButton){m_clicked = event->globalPos();update();}else if(event->button() == Qt::RightButton){myCanvas *canvase = (myCanvas *)parent();canvase->m_ctrlWidgets.removeOne(this);this->deleteLater();canvase->update();}}void myCtrlPoint::mouseMoveEvent(QMouseEvent *event)
{if(m_clicked == QPoint()){QWidget::mouseMoveEvent(event);}else{QPoint cur = event->globalPos();QPoint dis = cur - m_clicked;m_clicked = cur;move(mapToParent(QPoint(0, 0)) + dis);((QWidget *)parent())->update();}
}void myCtrlPoint::mouseReleaseEvent(QMouseEvent *event)
{m_clicked = QPoint();update();
}QPoint myCtrlPoint::getPosition()
{return mapToParent(rect().center());
}void myCtrlPoint::setPostion(const QPoint &point)
{QPoint target = point - rect().topLeft();move(target);}
6. 继承QWidget,实现一块画布
#ifndef MYCANVAS_H
#define MYCANVAS_H#include <QWidget>#include"zmBezier.h"class myCtrlPoint;
class myCanvas : public QWidget
{friend class myCtrlPoint;Q_OBJECT
public:explicit myCanvas(QWidget *parent = nullptr);~myCanvas();protected:void paintEvent(QPaintEvent *event) override;void mouseDoubleClickEvent(QMouseEvent *event) override;private:zmBezier m_curve;double m_points[1024][3]; //不想paintEvent中动态分配内存QVector<myCtrlPoint *>m_ctrlWidgets;
};#endif // MYCANVAS_H
6.1 构造时随机生成4个控制点
6.2 绘制事件中绘制控制点之间的连线、绘制Bezier曲线
6.3 鼠标左键双击空白处会添加一个控制点
6.4 因为不想再绘制事件中动态分配内存,所以用了一个比较大的数组
6.5 控制点是画布的友元类,方便控制点删除自己
#include"myCanvas.h"
#include"myCtrlPoint.h"#include<QTime>
#include<QDebug>
#include<QPainter>
#include<QMouseEvent>
#include<QRandomGenerator>myCanvas::myCanvas(QWidget *parent): QWidget(parent)
{QRandomGenerator random(QTime::currentTime().second());for(int i = 0; i < 4; i++){myCtrlPoint *ctrl = new myCtrlPoint(this);m_ctrlWidgets.append(ctrl);ctrl->setPostion(QPoint(random.generateDouble() * 400, random.generateDouble() * 400));}resize(500, 500);
}myCanvas::~myCanvas()
{}void myCanvas::paintEvent(QPaintEvent *event)
{QWidget::paintEvent(event);QPainter painter(this);painter.drawText(20, 20, "1.左键拖动控制点");painter.drawText(20, 40, "2.右键删除控制点");painter.drawText(20, 60, "3.左键双击空白处添加控制点");int n = m_ctrlWidgets.count();if(n){painter.setPen(QPen(Qt::blue, 1, Qt::DotLine));for(int i = 0; i < n - 1; i++){painter.drawLine(m_ctrlWidgets[i]->getPosition(), m_ctrlWidgets[i + 1]->getPosition());}// double (*ctrls)[3] = new double[n][3]; 尽量别动态分配了,下面限制下点数if(n > 1024) {n = 1024;}for(int i = 0; i < n; i++){
// ctrls[i][0] = m_ctrlWidgets[i]->getPosition().x();
// ctrls[i][1] = m_ctrlWidgets[i]->getPosition().y();
// ctrls[i][2] = 0;m_points[i][0] = m_ctrlWidgets[i]->getPosition().x();m_points[i][1] = m_ctrlWidgets[i]->getPosition().y();m_points[i][2] = 0;}m_curve.setCtrlPoints(n, m_points);
// m_curve.setCtrlPoints(n, ctrls);
// delete [] ctrls;int request = 100;
// double (*points)[3] = new double[request][3];// m_curve.getCurve(request, points);m_curve.getCurve(request, m_points);painter.setPen(QPen(Qt::green, 1));for(int i = 0; i < request - 1; i++) {painter.drawLine(m_points[i][0], m_points[i][1],m_points[i + 1][0], m_points[i + 1][1]);}// delete [] points;}
}void myCanvas::mouseDoubleClickEvent(QMouseEvent *event)
{if(event->button() == Qt::LeftButton){QPoint point = event->pos();myCtrlPoint *ctrl = new myCtrlPoint(this);m_ctrlWidgets.append(ctrl);ctrl->setPostion(point);ctrl->show();update();}}
7.直接显示画布
#include<QApplication>#include"myCanvas.h"int main(int argc, char *argv[])
{QApplication a(argc, argv);myCanvas camvas;camvas.show();return a.exec();
}