计算机图形学：实验四——Bezier曲线的绘制

1. 实验目的
练习Bezier曲线的绘制和de Casteljau算法。
2. 实验内容和要求
按要求完成如下一个作业，提交纸质实验报告，同时提交实验报告和代码的电子版。
实现Bezier曲线的de Casteljau递推算法，能够对任意介于0和1之间的参数t计算Bezier曲线上的点，然后依次连接这些点生成Bezier曲线。要求：
(1) .对[0,1]参数区间进行100等分。
(2) 控制点的数目至少为5个，即Bezier曲线的次数不低于4次。
(3) de Casteljau算法用一个函数单独实现。
(4) 绘制Bezier曲线的同时还要绘制其控制多边形。
(5) 至少绘制两条Bezier曲线，具有不同的次数，颜色和曲线宽度。
3. 算法描述
生成Bezier曲线：用一个二维数组存放初始的五个点，按照下图，每有一个t，用下面图示方法将生成一个点。将t从0到1进行100等分，所以最后会生成100个点，将这些点存放在一个新的二维数组中，最后使用glBegin(GL_LINE_STRIP)将这些点连接起来，就可已生成Bezier曲线。
第一次递推的结果应保存在一个新的二维数组中，之后用这个新的二维数组递推就可以。
生成折线：直接将最开始给的那几个点用glBegin(GL_LINE_STRIP)连接起来就可以生成折线（多边形）。
de Casteljau递推算法示意图：

4. 源程序代码

#include"pch.h"
#include<GL/glut.h>
#include<stdio.h>

GLfloat Point[5][2] = { {-200,0},{-300,100},{0,200},{200,100},{0,-150} };   //存放最初的五个点
GLfloat Point_x[6][2] = { {-300,0},{-400,200},{0,300},{300,200},{0,-200},{-100,-250} };
GLfloat Point_2[100][2];
GLfloat Point_3[100][2];    //存放最后生成的一百个点


void init(void) {
	glClearColor(1.0, 1.0, 1.0, 1.0);

	glMatrixMode(GL_PROJECTION);
	gluOrtho2D(-400.0, 400.0, -300.0, 300.0);
}
//绘制曲线函数
void setPixel(GLfloat arry[][2]) {
	int i;
	glBegin(GL_LINE_STRIP);
	for (i = 0; i < 100; i++) {
		glVertex2f(arry[i][0], arry[i][1]);
	}
	glEnd();
}
//de Casteljau算法的实现
void deCasteljau(GLfloat Point_1[][2], GLint n) {
	int i, u = 0, j;
	for (double t = 0.0; t < 1; t += 0.01) {
		for (i = 1; i < n; ++i) {
			for (j = 0; j < n - i; ++j) {
				// i==1时,第一次迭代,由已知控制点计算 
				if (i == 1) {
					Point_2[j][0] = Point_1[j][0] * t + Point_1[j + 1][0] * (1 - t);
					Point_2[j][1] = Point_1[j][1] * t + Point_1[j + 1][1] * (1 - t);
					continue;
				}
				// i != 1时,通过上一次迭代的结果计算  
				Point_2[j][0] = Point_2[j][0] * t + Point_2[j + 1][0] * (1 - t);
				Point_2[j][1] = Point_2[j][1] * t + Point_2[j + 1][1] * (1 - t);
			}
		}
		Point_3[u][0] = Point_2[0][0];
		Point_3[u][1] = Point_2[0][1];
		u++;
	}
	setPixel(Point_3);
}


void lineSegment(void) {
	int j;
	glClear(GL_COLOR_BUFFER_BIT);
	glColor3f(1, 0, 0);
	glLineWidth(3.0);
	//将几个点连接起来
	glBegin(GL_LINE_STRIP);
	for (j = 0; j < 5; j++) {
		glVertex2f(Point[j][0], Point[j][1]);
	}
	glEnd();
	deCasteljau(Point, 5);
	glColor3f(0, 0, 1);
	glLineWidth(4.0);
	//将几个点连接起来
	glBegin(GL_LINE_STRIP);
	for (j = 0; j < 6; j++) {
		glVertex2f(Point_x[j][0], Point_x[j][1]);
	}
	glEnd();
	deCasteljau(Point_x, 6);

	glFlush();
}

void main(int argc, char** argv) {
	glutInit(&argc, argv);
	glutInitDisplayMode(GLUT_SINGLE | GLUT_RGB);
	glutInitWindowPosition(100, 100);
	glutInitWindowSize(800, 600);
	glutCreateWindow("de Casteljau算法");

	init();
	glutDisplayFunc(lineSegment);

	glutMainLoop();
}

5.运行结果