第1关:Bezier曲线和B样条曲线
一.任务描述
根据提示,在右侧修改代码,并自己绘制出图形。平台会对你编写的代码进行测试。
1.本关任务
了解曲线的生成原理; 掌握几种常见的曲线生成算法,利用VC+OpenGL实现Bezier曲线生成算法和B样条曲线算法。
2.预期输出
3.渲染要求
(1).背景色为黑色,用 glclearcolor 来完成; (2).结合公式及代码示范代码了解Bezier曲线生成原理与算法实现; (3).参考代码文件中的Bezier曲线生成代码,根据B样条公式,修改B样条曲线代码来实现B样条的生成。要求B样条的控制顶点坐标为Point pt2[4]={ { 450, 100}, { 540,300},{650, 320}, {690, 120} }。
二.相关知识
Bezier曲线是通过一组多边形折线的顶点来定义的。如果折线的顶点固定不变,则由其定义的Bezier曲线是唯一的。在折线的各顶点中,只有第一点和最后一点在曲线上且作为曲线的起始点和终止点,其他的点用于控制曲线的形状及阶次。曲线的形状趋向于多边形折线的形状,要修改曲线,只要修改折线的各顶点就可以了。因此,多边形折线又称Bezier曲线的控制多边形,其顶点称为控制点。 三次Bezier曲线,有四个控制点,其数学表示如下:
开始你的任务吧,祝你成功!
三.实验代码
// 提示:在合适的地方修改或添加代码
#include <GL/freeglut.h>
#include<stdio.h>
#include <stdlib.h>
#include <vector>
// 评测代码所用头文件-开始
#include<opencv2/core/core.hpp>
#include<opencv2/highgui/highgui.hpp>
#include<opencv2/imgproc/imgproc.hpp>
// 评测代码所用头文件-结束
using namespace std;
struct Point
{
int x, y;
};
Point bz[11], bspt[11]; //bz为Bezier曲线,bspt为B样条曲线
int nInput;
//Bezier曲线控制点
Point pt[4] = { { 50, 100}, { 140,300},{250, 320}, {290, 120} };
//B样条曲线控制点
/**********请在此处添加你的代码*******/
/************** Begin ************/
Point pt2[4]={{ 450, 100}, { 540,300},{650, 320}, {690, 120}};
/*************** end **************/
void CalcBZPoints() //Bezier曲线算法
{
float a0, a1, a2, a3, b0, b1, b2, b3;
a0 = pt[0].x;
a1 = -3 * pt[0].x + 3 * pt[1].x;
a2 = 3 * pt[0].x - 6 * pt[1].x + 3 * pt[2].x;
a3 = -pt[0].x + 3 * pt[1].x - 3 * pt[2].x + pt[3].x;
b0 = pt[0].y;
b1 = -3 * pt[0].y + 3 * pt[1].y;
b2 = 3 * pt[0].y - 6 * pt[1].y + 3 * pt[2].y;
b3 = -pt[0].y + 3 * pt[1].y - 3 * pt[2].y + pt[3].y;
float t = 0;
float dt = 0.01;
for (int i = 0; t < 1.1; t += 0.1, i++)
{
bz[i].x = a0 + a1 * t + a2 * t * t + a3 * t * t * t;
bz[i].y = b0 + b1 * t + b2 * t * t + b3 * t * t * t;
}
}
void CalcBSPoints() //B样条曲线
{
/**********请在此处添加你的代码*******/
/************** Begin ************/
float a0, a1, a2, a3, b0, b1, b2, b3;
a3 = -1.0/6*pt2[0].x+1.0/2*pt2[1].x-1.0/2*pt2[2].x+1.0/6*pt2[3].x;
a2 = 1.0/2 * pt2[0].x - pt2[1].x+1.0/2*pt2[2].x;
a1 = -1.0/2 * pt2[0].x + 1.0/2 * pt2[2].x;
a0 = 1.0/6*pt2[0].x +2.0/3*pt2[1].x+1.0/6*pt2[2].x;
b3= -1.0/6*pt2[0].y+1.0/2*pt2[1].y-1.0/2*pt2[2].y+1.0/6*pt2[3].y;
b2 = 1.0/2 * pt2[0].y - pt2[1].y+1.0/2*pt2[2].y;
b1 = -1.0/2 * pt2[0].y + 1.0/2 * pt2[2].y;
b0=1.0/6* pt2[0].y +2.0/3*pt2[1].y+1.0/6*pt2[2].y;
float t = 0;
float dt = 0.01;
for (int i = 0; t < 1.1; t += 0.1, i++)
{
bspt[i].x =a0 + a1 * t + a2 * t * t + a3 * t * t * t;
bspt[i].y =b0 + b1 * t + b2 * t * t + b3 * t * t * t;
}
}
/************* end ************/
void ControlPoint()
{
glPointSize(2);
for (int i = 0; i < 4; i++)
{
glBegin(GL_POINTS);
glColor3f(1.0f, 0.0f, 0.0f);
glVertex2i(pt[i].x, pt[i].y);
glEnd();
}
}
void PolylineGL(Point* pt, int num)
{
glBegin(GL_LINE_STRIP);
for (int i = 0; i < num; i++)
{
glColor3f(0.0f, 1.0f, 0.0f);
glVertex2i(pt[i].x, pt[i].y);
}
glEnd();
}
void PolylineGL1(Point* pt, int num)
{
glBegin(GL_LINE_STRIP);
for (int i = 0; i < num; i++)
{
glColor3f(1.0f, 1.0f, 1.0f);
glVertex2i(pt[i].x, pt[i].y);
}
glEnd();
}
void display()
{
glClear(GL_COLOR_BUFFER_BIT);
glColor3f(1.0f, 1.0f, 1.0f);
ControlPoint(); //画4个控制点
PolylineGL1(pt, 4); //画4个点之间的线段
CalcBZPoints();
PolylineGL(bz, 11);
PolylineGL1(pt2, 4);
CalcBSPoints();
PolylineGL(bspt, 11);
glFlush();
}
void init()
{
glClearColor(0.0, 0.0, 0.0, 0.0);
glShadeModel(GL_SMOOTH);
}
void reshape(int w, int h)
{
glViewport(0, 0, (GLsizei)w, (GLsizei)h);
glMatrixMode(GL_PROJECTION);
glLoadIdentity();
gluOrtho2D(0.0, (GLdouble)w, 0.0, (GLdouble)h);
}
int main(int argc, char* argv[])
{
glutInit(&argc, argv);
glutInitWindowPosition(100, 100);
glutInitWindowSize( 800 , 400 ); //设置初始化窗口大小
glutCreateWindow("三维观察");
init();
glutDisplayFunc(display);
glutReshapeFunc(reshape);
glutMainLoopEvent();
/*************以下为评测代码,与本次实验内容无关,请勿修改**************/
GLubyte* pPixelData = (GLubyte*)malloc(800 * 400 * 3);//分配内存
GLint viewport[4] = { 0 };
glReadBuffer(GL_FRONT);
glPixelStorei(GL_UNPACK_ALIGNMENT, 4);
glGetIntegerv(GL_VIEWPORT, viewport);
glReadPixels(viewport[0], viewport[1], viewport[2], viewport[3], GL_RGB, GL_UNSIGNED_BYTE, pPixelData);
cv::Mat img;
std::vector<cv::Mat> imgPlanes;
img.create(400, 800, CV_8UC3);
cv::split(img, imgPlanes);
for (int i = 0; i < 400; i++) {
unsigned char* plane0Ptr = imgPlanes[0].ptr<unsigned char>(i);
unsigned char* plane1Ptr = imgPlanes[1].ptr<unsigned char>(i);
unsigned char* plane2Ptr = imgPlanes[2].ptr<unsigned char>(i);
for (int j = 0; j < 800; j++) {
int k = 3 * (i * 800 + j);
plane2Ptr[j] = pPixelData[k];
plane1Ptr[j] = pPixelData[k + 1];
plane0Ptr[j] = pPixelData[k + 2];
}
}
cv::merge(imgPlanes, img);
cv::flip(img, img, 0);
cv::namedWindow("openglGrab");
cv::imshow("openglGrab", img);
//cv::waitKey();
cv::imwrite("../img_step3/test.jpg", img);
return 0;
}