【作业】计算机图形学——Casteljau算法绘制Bézier曲线（cpp）

程序功能：给定平面上15个不同的点(数据文件输入)，用Casteljau算法编程生成平面Bézier曲线。

程序语言：c++

程序说明：

1. 头文件引入：程序使用了iostream、fstream、vector以及graphics.h和conio.h头文件。其中iostream、fstream和vector用于处理文件读取和存储控制点，而graphics.h和conio.h则用于图形的显示和用户输入的处理。
2. drawLine函数：用于在图形窗口中绘制直线，参数包括起始点的坐标、终点的坐标和颜色。
3. casteljau函数：实现了Casteljau算法，用于计算Bézier曲线上的点。该函数接受控制点集合、参数t以及输出参数x和y，通过递推的方式计算出曲线上的点的坐标。
4. main函数：程序的入口函数。首先从文件中读取控制点的坐标，并存储在control_points向量中。然后初始化图形显示窗口，绘制控制多边形和Bézier曲线，最后等待用户按下任意键来关闭图形窗口。
5. 文件操作：程序通过ifstream从指定路径的文件中读取控制点的坐标，文件格式为每行两个浮点数表示一个点的坐标。如果无法打开文件，则输出错误信息并退出程序。
6. 图形绘制：利用initgraph初始化图形窗口，利用drawLine函数绘制控制多边形和Bézier曲线，最后通过_getch等待用户按键并关闭图形窗口。

程序代码：

#include <iostream>
#include <fstream>
#include <vector>
#include <graphics.h>
#include <conio.h>

using namespace std;

// 直线绘制
void drawLine(int x1, int y1, int x2, int y2, int color)
{
    setcolor(color);
    line(x1, y1, x2, y2);
}

// Casteljau算法：计算Bézier曲线上的点
void casteljau(const vector<pair<float, float>>& control_points, float t, float& x, float& y)
{
    vector<pair<float, float>> points = control_points;
    int n = points.size();

    for (int r = 1; r < n; ++r)
    {
        for (int i = 0; i < n - r; ++i)
        {
            points[i].first = (1 - t) * points[i].first + t * points[i + 1].first;
            points[i].second = (1 - t) * points[i].second + t * points[i + 1].second;
        }
    }

    x = points[0].first;
    y = points[0].second;
}

int main()
{
    ifstream infile("F:/Program/CG_HM5/HM5/points.txt");
    if (!infile)
    {
        cout << "无法打开文件！" << endl;
        return 1;
    }

    vector<pair<float, float>> control_points; // 存储控制点

    float x, y;
    while (infile >> x >> y)
    {
        control_points.push_back(make_pair(x, y));
    }
    infile.close();

    // 初始化图形显示窗口
    initgraph(800, 800);

    // 绘制控制多边形
    for (size_t i = 0; i < control_points.size() - 1; ++i)
    {
        drawLine(control_points[i].first, control_points[i].second, control_points[i + 1].first, control_points[i + 1].second, GREEN);
    }

    // 计算并绘制Bézier曲线
    const int num_points = 100; // 曲线上的点数
    for (int i = 0; i < num_points - 1; ++i)
    {
        float t1 = (float)i / (num_points - 1);
        float t2 = (float)(i + 1) / (num_points - 1);
        float x1, y1, x2, y2;
        casteljau(control_points, t1, x1, y1);
        casteljau(control_points, t2, x2, y2);
        drawLine(static_cast<int>(x1), static_cast<int>(y1), static_cast<int>(x2), static_cast<int>(y2), WHITE);
    }

    cout << "按下任意键关闭图形窗口…";
    _getch();
    closegraph();

    return 0;
}

仅供学习交流参考