Bezier曲线原理及实现代码(c++)
2017-05-24 16:36
399 查看
转载于:http://devres.zoomquiet.io/data/20110728232822/index.html
一、原理:
贝塞尔曲线于1962年,由法国工程师皮埃尔·贝塞尔(Pierre
Bézier)所广泛发表,他运用贝塞尔曲线来为汽车的主体进行设计。贝塞尔曲线最初由 Paul
de Casteljau 于1959年运用 de
Casteljau 算法开发,以稳定数值的方法求出贝塞尔曲线。
线性贝塞尔曲线
给定点 P0、P1,线性贝塞尔曲线只是一条两点之间的直线。这条线由下式给出:
![](http://devres.zoomquiet.io/data/20110728232822/05c4210c69ffb1358ceb8eb83a1a06fe.png)
且其等同于线性插值。
二次方贝塞尔曲线的路径由给定点 P0、P1、P2 的函数 B(t)
追踪:
![](http://devres.zoomquiet.io/data/20110728232822/8adc5cc34ea9649d6e546043fce9c407.png)
。TrueType 字型就运用了以贝塞尔样条组成的二次贝塞尔曲线。
P0、P1、P2、P3 四个点在平面或在三维空间中定义了三次方贝塞尔曲线。曲线起始于 P0 走向 P1,并从 P2 的方向来到 P3。一般不会经过 P1 或 P2;这两个点只是在那里提供方向资讯。 P0 和 P1 之间的间距,决定了曲线在转而趋进P3 之前,走向 P2 方向的“长度有多长”。
曲线的参数形式为:
![](http://devres.zoomquiet.io/data/20110728232822/597ecc5022fa7ab65509d5edfa9c148c.png)
。现代的成象系统,如 PostScript、Asymptote 和 Metafont,运用了以贝塞尔样条组成的三次贝塞尔曲线,用来描绘曲线轮廓。
P0、P1、…、Pn,其贝塞尔曲线即
![](http://devres.zoomquiet.io/data/20110728232822/5cdbbad5d4698c786a3f02a1b9fb5f2d.png)
。例如 :
![](http://devres.zoomquiet.io/data/20110728232822/db41ef611ab1eedb0606d9bce9012c30.png)
。如上公式可如下递归表达: 用
![](http://devres.zoomquiet.io/data/20110728232822/12482104eb5115c14219436749bcbbc0.png)
表示由点 P0、P1、…、Pn 所决定的贝塞尔曲线。则
![](http://devres.zoomquiet.io/data/20110728232822/99a7fa8da2d1abbcb3fb8a5723eb75bb.png)
用平常话来说, 阶贝塞尔曲线之间的插值。
一些关于参数曲线的术语,有
![](http://devres.zoomquiet.io/data/20110728232822/3d8f42d9608ecbd44e2f57758d159eb2.png)
即多项式
![](http://devres.zoomquiet.io/data/20110728232822/9dd8bd0a2b2706e270e89cbf22626e0e.png)
又称作 n 阶的伯恩斯坦基底多项式,定义
00 = 1。
点 Pi 称作贝塞尔曲线的控制点。多边形以带有线的贝塞尔点连接而成,起始于 P0 并以 Pn 终止,称作贝塞尔多边形(或控制多边形)。贝塞尔多边形的凸包(convex
hull)包含有贝塞尔曲线。
线性贝塞尔曲线函数中的 t 会经过由 P0 至P1 的 B(t) 所描述的曲线。例如当 t=0.25 时,B(t)
即一条由点 P0 至 P1 路径的四分之一处。就像由 0 至 1 的连续 t,B(t) 描述一条由 P0 至 P1 的直线。
![](http://devres.zoomquiet.io/data/20110728232822/bezier_1_big.gif)
为建构二次贝塞尔曲线,可以中介点 Q0 和 Q1 作为由 0 至 1 的 t:
由 P0 至 P1 的连续点 Q0,描述一条线性贝塞尔曲线。
由 P1 至 P2 的连续点 Q1,描述一条线性贝塞尔曲线。
由 Q0 至 Q1 的连续点 B(t),描述一条二次贝塞尔曲线。
为建构高阶曲线,便需要相应更多的中介点。对于三次曲线,可由线性贝塞尔曲线描述的中介点 Q0、Q1、Q2,和由二次曲线描述的点 R0、R1 所建构:
对于四次曲线,可由线性贝塞尔曲线描述的中介点 Q0、Q1、Q2、Q3,由二次贝塞尔曲线描述的点 R0、R1、R2,和由三次贝塞尔曲线描述的点 S0、S1 所建构:
P(t)=(1-t)P0+tP1 ,
![](http://devres.zoomquiet.io/data/20110728232822/mo6_44.gif)
。
矩阵表示为:
![](http://devres.zoomquiet.io/data/20110728232822/mo6_5f40.gif)
,
![](http://devres.zoomquiet.io/data/20110728232822/mo6_41.gif)
。
P(t)=(1-t)2P0+2t(1-t)P1+t2P2,
![](http://devres.zoomquiet.io/data/20110728232822/mo6_44.gif)
。
矩阵表示为:
![](http://devres.zoomquiet.io/data/20110728232822/mo6_42.gif)
,
![](http://devres.zoomquiet.io/data/20110728232822/mo6_41.gif)
。
P(t)=(1-t)3P0+3t(1-t)2P1+3t2(1-t)P2+t3P3
矩阵表示为:
![](http://devres.zoomquiet.io/data/20110728232822/mo6_43.gif)
,
![](http://devres.zoomquiet.io/data/20110728232822/mo6_41.gif)
。
(6-3-2)
![](http://devres.zoomquiet.io/data/20110728232822/mo6_45.gif)
,
![](http://devres.zoomquiet.io/data/20110728232822/mo6_41.gif)
。
在(6-3-2)式中,Mn+1是一个n+1阶矩阵,称为n次Bezier矩阵。
![](http://devres.zoomquiet.io/data/20110728232822/dingyi.gif)
(6-3-3)
![](http://devres.zoomquiet.io/data/20110728232822/mo6_46.gif)
。
其中,
![](http://devres.zoomquiet.io/data/20110728232822/mo6_47.gif)
利用(6-3-3)式,我们可以得到任意次Bezier矩阵的显式表示,例如4次和5次Bezier矩阵为:
![](http://devres.zoomquiet.io/data/20110728232822/mo6_48.gif)
,
![](http://devres.zoomquiet.io/data/20110728232822/mo6_49.gif)
可以证明,n次Bezier矩阵还可以表示为递推的形式:
![](http://devres.zoomquiet.io/data/20110728232822/dingyi.gif)
(6-3-4)
![](http://devres.zoomquiet.io/data/20110728232822/mo6_50.gif)
二、算法(c++)
工程目录是:Win32App
vc6.0
#include<windows.h>
#include<stdlib.h>
#include<time.h>
#define NUM 10
LRESULT CALLBACK Winproc(HWND,UINT,WPARAM,LPARAM);
int WINAPI WinMain(HINSTANCE hInstance,HINSTANCE hPrevInstanc,LPSTR lpCmdLine,int nShowCmd)
{
MSG msg;
static TCHAR szClassName[] = TEXT("::Bezier样条计算公式由法国雷诺汽车公司的工程师Pierm Bezier于六十年代提出");
HWND hwnd;
WNDCLASS wc;
wc.cbClsExtra =0;
wc.cbWndExtra =0;
wc.hbrBackground = (HBRUSH)GetStockObject(WHITE_BRUSH);
wc.hCursor = LoadCursor(NULL,IDC_ARROW);
wc.hIcon = LoadIcon(NULL,IDI_APPLICATION);
wc.hInstance = hInstance;
wc.lpfnWndProc = Winproc;
wc.lpszClassName = szClassName;
wc.lpszMenuName = NULL;
wc.style = CS_HREDRAW|CS_VREDRAW;
if(!RegisterClass(&wc))
{
MessageBox(NULL,TEXT("注册失败"),TEXT("警告框"),MB_ICONERROR);
return 0;
}
hwnd = CreateWindow(szClassName,szClassName,
WS_OVERLAPPEDWINDOW,
CW_USEDEFAULT,CW_USEDEFAULT,
CW_USEDEFAULT,CW_USEDEFAULT,
NULL,NULL,hInstance,NULL);
ShowWindow(hwnd,SW_SHOWMAXIMIZED);
UpdateWindow(hwnd);
while(GetMessage(&msg,NULL,0,0))
{
TranslateMessage(&msg);
DispatchMessage(&msg);
}
return msg.wParam;
}
LRESULT CALLBACK Winproc(HWND hwnd,UINT message, WPARAM wparam,LPARAM lparam)
{
HDC hdc;
static POINT pt[NUM];
TEXTMETRIC tm;
static int cxClient,cyClient;
HPEN hpen;
int i,j,k,n,t;
switch(message)
{
case WM_CREATE:
static int cxchar;
hdc = GetDC(hwnd);
GetTextMetrics(hdc,&tm);
cxchar = tm.tmAveCharWidth;
ReleaseDC(hwnd,hdc);
case WM_SIZE:
cxClient = LOWORD(lparam);
cyClient = HIWORD(lparam);
return 0;
case WM_PAINT:
hdc = GetDC(hwnd);
srand(time(0));
Rectangle(hdc,0,0,cxClient,cyClient);
for(i=0; i<500; i++)
{
SelectObject(hdc,GetStockObject(WHITE_PEN));
PolyBezier(hdc,pt,NUM);
for(j=0; j<NUM; j++)
{
pt[j].x = rand()%cxClient;
pt[j].y = rand()%cyClient;
}
hpen = CreatePen(PS_INSIDEFRAME,3,RGB(rand()%256,rand()%256,rand()%256));
DeleteObject(SelectObject(hdc,hpen));
PolyBezier(hdc,pt,NUM);
for(k=0; k<50000000;k++);
}
for(i=0; i<100;i++)
{
Ellipse(hdc,rand()%cxClient,rand()%cyClient,rand()%cxClient,rand()%cyClient);
Pie(hdc,j=rand()%cxClient,k=rand()%cyClient,n=rand()%cxClient,t=rand()%cyClient,rand()%cxClient,rand()%cyClient,rand()%cxClient,rand()%cyClient) ;
}
if((n=(n+j)/2)>cxchar*20) n=cxchar*20;
SetTextColor(hdc,RGB(rand()%256,rand()%256,rand()%256));
TextOut(hdc,n/2,(t+k)/2,TEXT("瑾以此向Pierm Bezier致敬!"),lstrlen(TEXT("瑾以此向Pierm Bezier致敬!")));
ReleaseDC(hwnd,hdc);
DeleteObject(hpen);
ValidateRect(hwnd,NULL);
return 0;
case WM_DESTROY:
PostQuitMessage(0);
return 0;
}
return DefWindowProc(hwnd,message,wparam,lparam);
}
Bezier曲线原理及实现代码(c++)
一、原理:贝塞尔曲线于1962年,由法国工程师皮埃尔·贝塞尔(Pierre
Bézier)所广泛发表,他运用贝塞尔曲线来为汽车的主体进行设计。贝塞尔曲线最初由 Paul
de Casteljau 于1959年运用 de
Casteljau 算法开发,以稳定数值的方法求出贝塞尔曲线。
线性贝塞尔曲线
给定点 P0、P1,线性贝塞尔曲线只是一条两点之间的直线。这条线由下式给出:
![](http://devres.zoomquiet.io/data/20110728232822/05c4210c69ffb1358ceb8eb83a1a06fe.png)
且其等同于线性插值。
追踪:
![](http://devres.zoomquiet.io/data/20110728232822/8adc5cc34ea9649d6e546043fce9c407.png)
。TrueType 字型就运用了以贝塞尔样条组成的二次贝塞尔曲线。
曲线的参数形式为:
![](http://devres.zoomquiet.io/data/20110728232822/597ecc5022fa7ab65509d5edfa9c148c.png)
。现代的成象系统,如 PostScript、Asymptote 和 Metafont,运用了以贝塞尔样条组成的三次贝塞尔曲线,用来描绘曲线轮廓。
一般化
P0、P1、…、Pn,其贝塞尔曲线即
![](http://devres.zoomquiet.io/data/20110728232822/5cdbbad5d4698c786a3f02a1b9fb5f2d.png)
。例如 :
![](http://devres.zoomquiet.io/data/20110728232822/db41ef611ab1eedb0606d9bce9012c30.png)
。如上公式可如下递归表达: 用
![](http://devres.zoomquiet.io/data/20110728232822/12482104eb5115c14219436749bcbbc0.png)
表示由点 P0、P1、…、Pn 所决定的贝塞尔曲线。则
![](http://devres.zoomquiet.io/data/20110728232822/99a7fa8da2d1abbcb3fb8a5723eb75bb.png)
用平常话来说, 阶贝塞尔曲线之间的插值。
![](http://devres.zoomquiet.io/data/20110728232822/3d8f42d9608ecbd44e2f57758d159eb2.png)
即多项式
![](http://devres.zoomquiet.io/data/20110728232822/9dd8bd0a2b2706e270e89cbf22626e0e.png)
又称作 n 阶的伯恩斯坦基底多项式,定义
00 = 1。
点 Pi 称作贝塞尔曲线的控制点。多边形以带有线的贝塞尔点连接而成,起始于 P0 并以 Pn 终止,称作贝塞尔多边形(或控制多边形)。贝塞尔多边形的凸包(convex
hull)包含有贝塞尔曲线。
线性贝塞尔曲线函数中的 t 会经过由 P0 至P1 的 B(t) 所描述的曲线。例如当 t=0.25 时,B(t)
即一条由点 P0 至 P1 路径的四分之一处。就像由 0 至 1 的连续 t,B(t) 描述一条由 P0 至 P1 的直线。
![](http://devres.zoomquiet.io/data/20110728232822/bezier_1_big.gif)
为建构二次贝塞尔曲线,可以中介点 Q0 和 Q1 作为由 0 至 1 的 t:
由 P0 至 P1 的连续点 Q0,描述一条线性贝塞尔曲线。
由 P1 至 P2 的连续点 Q1,描述一条线性贝塞尔曲线。
由 Q0 至 Q1 的连续点 B(t),描述一条二次贝塞尔曲线。
![]() | ![]() | |
![]() | ![]() | |
![]() | ![]() | |
![](http://devres.zoomquiet.io/data/20110728232822/mo6_44.gif)
。
矩阵表示为:
![](http://devres.zoomquiet.io/data/20110728232822/mo6_5f40.gif)
,
![](http://devres.zoomquiet.io/data/20110728232822/mo6_41.gif)
。
P(t)=(1-t)2P0+2t(1-t)P1+t2P2,
![](http://devres.zoomquiet.io/data/20110728232822/mo6_44.gif)
。
矩阵表示为:
![](http://devres.zoomquiet.io/data/20110728232822/mo6_42.gif)
,
![](http://devres.zoomquiet.io/data/20110728232822/mo6_41.gif)
。
P(t)=(1-t)3P0+3t(1-t)2P1+3t2(1-t)P2+t3P3
矩阵表示为:
![](http://devres.zoomquiet.io/data/20110728232822/mo6_43.gif)
,
![](http://devres.zoomquiet.io/data/20110728232822/mo6_41.gif)
。
(6-3-2)
![](http://devres.zoomquiet.io/data/20110728232822/mo6_45.gif)
,
![](http://devres.zoomquiet.io/data/20110728232822/mo6_41.gif)
。
在(6-3-2)式中,Mn+1是一个n+1阶矩阵,称为n次Bezier矩阵。
![](http://devres.zoomquiet.io/data/20110728232822/dingyi.gif)
(6-3-3)
![](http://devres.zoomquiet.io/data/20110728232822/mo6_46.gif)
。
其中,
![](http://devres.zoomquiet.io/data/20110728232822/mo6_47.gif)
利用(6-3-3)式,我们可以得到任意次Bezier矩阵的显式表示,例如4次和5次Bezier矩阵为:
![](http://devres.zoomquiet.io/data/20110728232822/mo6_48.gif)
,
![](http://devres.zoomquiet.io/data/20110728232822/mo6_49.gif)
可以证明,n次Bezier矩阵还可以表示为递推的形式:
![](http://devres.zoomquiet.io/data/20110728232822/dingyi.gif)
(6-3-4)
![](http://devres.zoomquiet.io/data/20110728232822/mo6_50.gif)
二、算法(c++)
工程目录是:Win32App
vc6.0
#include<windows.h>
#include<stdlib.h>
#include<time.h>
#define NUM 10
LRESULT CALLBACK Winproc(HWND,UINT,WPARAM,LPARAM);
int WINAPI WinMain(HINSTANCE hInstance,HINSTANCE hPrevInstanc,LPSTR lpCmdLine,int nShowCmd)
{
MSG msg;
static TCHAR szClassName[] = TEXT("::Bezier样条计算公式由法国雷诺汽车公司的工程师Pierm Bezier于六十年代提出");
HWND hwnd;
WNDCLASS wc;
wc.cbClsExtra =0;
wc.cbWndExtra =0;
wc.hbrBackground = (HBRUSH)GetStockObject(WHITE_BRUSH);
wc.hCursor = LoadCursor(NULL,IDC_ARROW);
wc.hIcon = LoadIcon(NULL,IDI_APPLICATION);
wc.hInstance = hInstance;
wc.lpfnWndProc = Winproc;
wc.lpszClassName = szClassName;
wc.lpszMenuName = NULL;
wc.style = CS_HREDRAW|CS_VREDRAW;
if(!RegisterClass(&wc))
{
MessageBox(NULL,TEXT("注册失败"),TEXT("警告框"),MB_ICONERROR);
return 0;
}
hwnd = CreateWindow(szClassName,szClassName,
WS_OVERLAPPEDWINDOW,
CW_USEDEFAULT,CW_USEDEFAULT,
CW_USEDEFAULT,CW_USEDEFAULT,
NULL,NULL,hInstance,NULL);
ShowWindow(hwnd,SW_SHOWMAXIMIZED);
UpdateWindow(hwnd);
while(GetMessage(&msg,NULL,0,0))
{
TranslateMessage(&msg);
DispatchMessage(&msg);
}
return msg.wParam;
}
LRESULT CALLBACK Winproc(HWND hwnd,UINT message, WPARAM wparam,LPARAM lparam)
{
HDC hdc;
static POINT pt[NUM];
TEXTMETRIC tm;
static int cxClient,cyClient;
HPEN hpen;
int i,j,k,n,t;
switch(message)
{
case WM_CREATE:
static int cxchar;
hdc = GetDC(hwnd);
GetTextMetrics(hdc,&tm);
cxchar = tm.tmAveCharWidth;
ReleaseDC(hwnd,hdc);
case WM_SIZE:
cxClient = LOWORD(lparam);
cyClient = HIWORD(lparam);
return 0;
case WM_PAINT:
hdc = GetDC(hwnd);
srand(time(0));
Rectangle(hdc,0,0,cxClient,cyClient);
for(i=0; i<500; i++)
{
SelectObject(hdc,GetStockObject(WHITE_PEN));
PolyBezier(hdc,pt,NUM);
for(j=0; j<NUM; j++)
{
pt[j].x = rand()%cxClient;
pt[j].y = rand()%cyClient;
}
hpen = CreatePen(PS_INSIDEFRAME,3,RGB(rand()%256,rand()%256,rand()%256));
DeleteObject(SelectObject(hdc,hpen));
PolyBezier(hdc,pt,NUM);
for(k=0; k<50000000;k++);
}
for(i=0; i<100;i++)
{
Ellipse(hdc,rand()%cxClient,rand()%cyClient,rand()%cxClient,rand()%cyClient);
Pie(hdc,j=rand()%cxClient,k=rand()%cyClient,n=rand()%cxClient,t=rand()%cyClient,rand()%cxClient,rand()%cyClient,rand()%cxClient,rand()%cyClient) ;
}
if((n=(n+j)/2)>cxchar*20) n=cxchar*20;
SetTextColor(hdc,RGB(rand()%256,rand()%256,rand()%256));
TextOut(hdc,n/2,(t+k)/2,TEXT("瑾以此向Pierm Bezier致敬!"),lstrlen(TEXT("瑾以此向Pierm Bezier致敬!")));
ReleaseDC(hwnd,hdc);
DeleteObject(hpen);
ValidateRect(hwnd,NULL);
return 0;
case WM_DESTROY:
PostQuitMessage(0);
return 0;
}
return DefWindowProc(hwnd,message,wparam,lparam);
}
相关文章推荐
- Bezier曲线原理及实现代码(c++)
- Bezier曲线原理及实现代码(c++)
- Bezier曲线原理及实现代码(c++)
- Bezier曲线原理及实现代码(c++)
- Bezier曲线原理及实现代码(c++)
- Bezier曲线原理及实现代码(c++)
- Bezier曲线原理及实现代码(c++)
- Bezier曲线原理及实现代码(c++)
- Bezier曲线原理及实现代码(c++)
- Bezier曲线原理及实现代码(c++)
- 归并排序的原理,c++代码实现
- UTF-8编码的字符串拆分成单字、获取UTF-8字符串的字符个数的代码及原理(c++实现)
- BP神经网络原理及C++代码实现
- 平衡二叉搜索树(AVL树)的原理及实现源代码(有图文详解和C++、Java实现代码)
- 伸展树的原理及实现源代码(有图文详解和C++实现代码)
- C/C++实现如下 根据Merge Sort原理, 自己实现的归并排序算法+详细注释+代码(C#,C/C++)
- B+树和R树的原理和C++代码实现
- 图像去雾之何凯明暗通道先验去雾算法原理及c++代码实现
- BP神经网络原理分析及c++代码实现(上)
- ID3决策树算法原理及C++实现(其中代码转自别人的博客)