Bezier曲线的绘制

原文地址：http://blog.csdn.net/mylovestart/article/details/8434310

Bezier曲线是参数多项式曲线,它由一组控制多边形折线(控制多边形)的顶点唯一定义,在控制多边形的各顶点中,只有第一个和最后一个顶点在曲线上,其他的顶点则用以定义曲线的导数,阶次和形状

Bezier曲线的数学基础是能够在第一个和最后一个顶点之间进行插值的一个多项式混合函数,对于有n+1个控制点的Bezier曲线段用参数方程表示如下：

式中pk(xk,yk,zk),k=0,1,2…..n是控制多边形的n+1个顶点,BENk,n(t)是Bernstein基函数

下面输入n个点并绘制Bezier曲线

[cpp] view plain copy print ?

1. #include<stdio.h> 
2. #include<string.h> 
3. #include<conio.h> 
4. #include<windows.h> 
5. #include<graphics.h>  
6. #include<math.h> 
7. #define N 1000 
8. int n;//控制点的个数 
9. struct point  //控制点的坐标 
10. { 
11.     double x; 
12.     double y; 
13. }point[N]; 
14. void init()  //输入控制点的坐标 
15. { 
16.     int i; 
17.     printf("please input the number of the points: "); 
18.     scanf("%d",&n); 
19.     printf("please input the location of the points\n"); 
20.     for(i=0;i<n;i++) 
21.     scanf("%lf %lf",&point[i].x,&point[i].y); 
22. } 
23. void sol1()  //绘制控制多边形的轮廓 
24. { 
25.     int i; 
26.     setcolor(RED); 
27.     for(i=0;i<n-1;i++) 
28.     line((int)point[i].x,(int)point[i].y,(int)point[i+1].x,(int)point[i+1].y); 
29. } 
30. double sol2(int nn,int k)  //计算多项式的系数C(nn,k) 
31. { 
32.     int i; 
33.     double sum=1; 
34.     for(i=1;i<=nn;i++) 
35.     sum*=i; 
36.     for(i=1;i<=k;i++) 
37.     sum/=i; 
38.     for(i=1;i<=nn-k;i++) 
39.     sum/=i; 
40.     return sum; 
41. } 
42. void sol3(double t)  //计算Bezier曲线上点的坐标 
43. { 
44.     double x=0,y=0,Ber; 
45.     int k; 
46.     for(k=0;k<n;k++) 
47.     { 
48.         Ber=sol2(n-1,k)*pow(t,k)*pow(1-t,n-1-k); 
49.         x+=point[k].x*Ber; 
50.         y+=point[k].y*Ber; 
51.     } 
52.     putpixel((int)x,(int)y,GREEN); 
53. } 
54. void sol4()  //根据控制点，求曲线上的m个点 
55. { 
56.     int m=500,i; 
57.     for(i=0;i<=m;i++) 
58.     sol3((double)i/(double)m); 
59. } 
60. int main() 
61. { 
62.     init(); 
63.     initgraph(640, 480);  
64.     sol1(); 
65.     sol4(); 
66.     getch(); 
67.     closegraph(); 
68.     return 0; 
69. } 

#include<stdio.h>
#include<string.h>
#include<conio.h>
#include<windows.h>
#include<graphics.h> 
#include<math.h>
#define N 1000
int n;//控制点的个数
struct point  //控制点的坐标
{
	double x;
	double y;
}point[N];
void init()  //输入控制点的坐标
{
	int i;
    printf("please input the number of the points: ");
	scanf("%d",&n);
	printf("please input the location of the points\n");
	for(i=0;i<n;i++)
	scanf("%lf %lf",&point[i].x,&point[i].y);
}
void sol1()  //绘制控制多边形的轮廓
{
    int i;
	setcolor(RED);
	for(i=0;i<n-1;i++)
	line((int)point[i].x,(int)point[i].y,(int)point[i+1].x,(int)point[i+1].y);
}
double sol2(int nn,int k)  //计算多项式的系数C(nn,k)
{
	int i;
	double sum=1;
    for(i=1;i<=nn;i++)
	sum*=i;
	for(i=1;i<=k;i++)
	sum/=i;
	for(i=1;i<=nn-k;i++)
	sum/=i;
	return sum;
}
void sol3(double t)  //计算Bezier曲线上点的坐标
{
    double x=0,y=0,Ber;
	int k;
	for(k=0;k<n;k++)
	{
		Ber=sol2(n-1,k)*pow(t,k)*pow(1-t,n-1-k);
		x+=point[k].x*Ber;
		y+=point[k].y*Ber;
	}
	putpixel((int)x,(int)y,GREEN);
}
void sol4()  //根据控制点，求曲线上的m个点
{
	int m=500,i;
	for(i=0;i<=m;i++)
	sol3((double)i/(double)m);
}
int main()
{
    init();
	initgraph(640, 480); 
	sol1();
    sol4();
	getch();
	closegraph();
	return 0;
}

给出一组测试用例

(红色为控制多边形,绿色为Bezier曲线)