3-4次Bezier曲面_2021秋季《计算机图形学》_基于《计算机图形学（第四版）》D.H.&M.P.B.&W.R.C.

最新推荐文章于 2024-03-15 09:46:33 发布
最新推荐文章于 2024-03-15 09:46:33 发布
阅读量976 收藏 2
点赞数
分类专栏： 计算机图形学 文章标签： c++ 图形学
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/rd142857/article/details/121447786
版权

一、实验目的

  • 利用 d e − C a s t e l j a u {\rm de-Casteljau} de−Casteljau递推算法绘制3-4次 B e z i e r {\rm Bezier} Bezier曲面

二、实验环境

  • Visual Studio 2019
  • Windows 10

三、算法分析与设计

曲线

将 t t t从 0 0 0取至 1 1 1（指与线段总长之比），点 P P P的轨迹就形成了一条由端点加权赋得的曲线。

示意图如下：
在这里插入图片描述

也就是每次由 n + 1 n+1 n+1个点推出 n n n个新点，最终得到目标点 P P P
P 0 , P 1 , P 2 , P 3 ⇒ Q 0 , Q 1 , Q 2 ⇒ R 0 , R 1 ⇒ P P_0,P_1,P_2,P_3\Rightarrow Q_0,Q_1,Q_2\Rightarrow R_0,R_1\Rightarrow P P0​,P1​,P2​,P3​⇒Q0​,Q1​,Q2​⇒R0​,R1​⇒P

曲面

将曲面视为由相互正交的线组成的网格，于是求曲面的过程可以转化为求组成网格的各曲线。

  1. 对每一列的 n + 1 n+1 n+1个点确定一条曲线
  2. 对每一行的 n + 1 n+1 n+1个点确定一条曲线

四、实验结果

三次 B e z i e r {\rm Bezier} Bezier曲面

在这里插入图片描述

四次 B e z i e r {\rm Bezier} Bezier曲面

在这里插入图片描述

五、附录

三次 B e z i e r {\rm Bezier} Bezier曲面

#include<GL/glew.h>
#include <GLFW/glfw3.h>
#include<bits/stdc++.h>

#define PI acos(-1)
#define MAXN 3600
#define M 1000

using namespace std;

GLfloat ctrlPoints[4][4][3] =
// ctrlPoints[i][j]：该点在网格的第i行第j列上
{
    {
        { -0.8f, -0.7f, 0.2f }, { -0.3f, -0.6f, 0.2f }, { 0.2f, -0.65f, 0.3f }, { 0.7f, -0.7f, 0.2f }
    },
    {
        { -0.9f, -0.2f, 0.3f }, { -0.3f, 0.2f, 0.5f }, { 0.3f, -0.2f, 0.4f }, { 0.75f, -0.2f, 0.3f },
    },
    {
        { -0.9f, 0.3f, 0.3f }, { -0.3f, 0.2f, 0.5f }, { 0.25f, 0.25f, 0.6f }, { 0.8f, 0.3f, 0.3f },
    },
    {
        { -0.8f, 0.8f, 0.1f }, { -0.3f, 0.8f, 0.2f }, { 0.2f, 0.85f, 0.1f }, { 0.7f, 0.8f, 0.1f },
    }
};

void plotPoint()
{
    glPointSize(5.0);
    glColor3f(1.0, 1.0, 0.0);
    for (int i = 0; i < 4; i++) 
    {
        glBegin(GL_POINTS);
        for (int j = 0; j < 4; j++)
            glVertex3fv(&ctrlPoints[i][j][0]);
        glEnd();
    }

    glColor3f(0.0, 1.0, 1.0);
    for (int i = 0; i < 4; i++) {
        // 固定行标，画横线
        glBegin(GL_LINE_STRIP);
        for (int j = 0; j < 4; j++)
            glVertex3fv(&ctrlPoints[i][j][0]);
        glEnd();

        // 固定列标，画竖线
        glBegin(GL_LINE_STRIP);
        for (int j = 0; j < 4; j++)
            glVertex3fv(&ctrlPoints[j][i][0]);
        glEnd();
    }
}

void bezierSurface()
{
    int n = 4;  // 控制点的数目
    float xarray[4];
    float yarray[4];
    float zarray[4];
    GLfloat ps[11][4][3];
    // 由4*4=16个点确定4*11=44个点，即对每列（共4列）的四个点确定一条曲线
    for (int v = 0; v < 4; v++)
    {
        int u = 0;
        for (double t = 0.0; t <= 1; t += 0.1)
        {
            for (int i = 1; i < n; ++i)
            {
                for (int j = 0; j < n - i; ++j)
                {
                    if (i == 1)
                    // 第一次，使用控制点
                    {
                        xarray[j] = ctrlPoints[j][v][0] * (1 - t) + ctrlPoints[j + 1][v][0] * t;
                        yarray[j] = ctrlPoints[j][v][1] * (1 - t) + ctrlPoints[j + 1][v][1] * t;
                        zarray[j] = ctrlPoints[j][v][2] * (1 - t) + ctrlPoints[j + 1][v][2] * t;
                        continue;
                    }
                    // 非第一次，更新点数组
                    xarray[j] = xarray[j] * (1 - t) + xarray[j + 1] * t;
                    yarray[j] = yarray[j] * (1 - t) + yarray[j + 1] * t;
                    zarray[j] = zarray[j] * (1 - t) + zarray[j + 1] * t;
                }

            }
            ps[u][v][0] = xarray[0];
            ps[u][v][1] = yarray[0];
            ps[u][v][2] = zarray[0];
            u++;
        }
    }

    //glColor3f(255.0, 0.0, 0.0);
    //for (int i = 0; i < 11; i++)
    //{
    //    glBegin(GL_LINE_STRIP);
    //    for (int j = 0; j < 4; j++)
    //        glVertex3fv(&ps[i][j][0]);
    //    glEnd();
    //}

    GLfloat ps1[11][11][3];
    // 由4*11=44个点确定11*11=121个点，即对每行（共11行）的四个点确定一条曲线
    for (int v = 0; v < 11; v++) {
        int u = 0;
        for (double t = 0.0; t <= 1; t += 0.1)
        {
            for (int i = 1; i < n; ++i)
            {
                for (int j = 0; j < n - i; ++j)
                {
                    if (i == 1)
                    {
                        xarray[j] = ps[v][j][0] * (1 - t) + ps[v][j + 1][0] * t;
                        yarray[j] = ps[v][j][1] * (1 - t) + ps[v][j + 1][1] * t;
                        zarray[j] = ps[v][j][2] * (1 - t) + ps[v][j + 1][2] * t;
                        continue;
                    } 
                    xarray[j] = xarray[j] * (1 - t) + xarray[j + 1] * t;
                    yarray[j] = yarray[j] * (1 - t) + yarray[j + 1] * t;
                    zarray[j] = zarray[j] * (1 - t) + zarray[j + 1] * t;
                }

            }
            ps1[v][u][0] = xarray[0];
            ps1[v][u][1] = yarray[0];
            ps1[v][u][2] = zarray[0];
            u++;
        }
    }

    glColor3f(1.0, 1.0, 1.0);
    for (int i = 0; i < 11; i++) 
    {
        glBegin(GL_LINE_STRIP);
        for (int j = 0; j < 11; j++)
            glVertex3fv(&ps1[i][j][0]);
        glEnd();

        glBegin(GL_LINE_STRIP);
        for (int j = 0; j < 11; j++)
            glVertex3fv(&ps1[j][i][0]);
        glEnd();
    }


    glFlush();
}


int main(void)
{
    GLFWwindow* window;

    /* Initialize the library */
    if (!glfwInit())
        return -1;

    /* Create a windowed mode window and its OpenGL context */
    window = glfwCreateWindow(960, 960, "bezierSurface", NULL, NULL);
    if (!window)
    {
        glfwTerminate();
        return -1;
    }

    /* Make the window's context current */
    glfwMakeContextCurrent(window);

    if (glewInit() != GLEW_OK)
        cout << "Error!" << endl;
    
    glRotatef(300.0, 1.0, 0.0, 0.0);

    /* Loop until the user closes the window */
    while (!glfwWindowShouldClose(window))
    {

        /* Render here */
        glClearColor(0.2f, 0.3f, 0.3f, 1.0f);
        glClear(GL_COLOR_BUFFER_BIT);
        plotPoint();
        bezierSurface();
        /* Swap front and back buffers */
        glfwSwapBuffers(window);

        /* Poll for and process events */
        glfwPollEvents();
    }

    glfwTerminate();
    return 0;
}

四次 B e z i e r {\rm Bezier} Bezier曲面

#include<GL/glew.h>
#include <GLFW/glfw3.h>
#include<bits/stdc++.h>

#define PI acos(-1)
#define MAXN 3600
#define M 1000

using namespace std;

GLfloat ctrlPoints[5][5][3] =
// ctrlPoints[i][j]：该点在网格的第i行第j列上
{
    {
        { -0.8f, -0.7f, 0.2f }, { -0.3f, -0.6f, 0.2f }, {0.0f,-0.6f,0.2f}, { 0.2f, -0.65f, 0.3f }, { 0.7f, -0.7f, 0.2f }
    },
    {
        { -0.9f, -0.2f, 0.3f }, { -0.3f, 0.0f, 0.5f }, {0.1f,-0.2f,0.3f},{ 0.3f, -0.2f, 0.4f }, { 0.75f, -0.2f, 0.3f },
    },
    {
        { -0.8f, 0.0f, 0.4f }, { -0.3f, 0.0f, 0.2f }, {0.0f,-0.1f,0.4f},{ 0.2f, 0.0f, 0.5f }, { 0.7f, 0.0f, 0.1f },
    },
    {
        { -0.9f, 0.3f, 0.3f }, { -0.3f, 0.2f, 0.5f }, {0.0f,0.2f,0.6f},{ 0.25f, 0.25f, 0.6f }, { 0.8f, 0.3f, 0.3f },
    },
    {
        { -0.8f, 0.8f, 0.1f }, { -0.3f, 0.8f, 0.2f }, {0.1f,0.8f,0.1f},{ 0.2f, 0.85f, 0.1f }, { 0.7f, 0.8f, 0.1f },
    },
};

void plotPoint()
{
    glPointSize(5.0);
    glColor3f(1.0, 1.0, 0.0);
    for (int i = 0; i < 5; i++) 
    {
        glBegin(GL_POINTS);
        for (int j = 0; j < 5; j++)
            glVertex3fv(&ctrlPoints[i][j][0]);
        glEnd();
    }

    glColor3f(0.0, 1.0, 1.0);
    for (int i = 0; i < 5; i++) {
        // 固定行标，画横线
        glBegin(GL_LINE_STRIP);
        for (int j = 0; j < 5; j++)
            glVertex3fv(&ctrlPoints[i][j][0]);
        glEnd();

        // 固定列标，画竖线
        glBegin(GL_LINE_STRIP);
        for (int j = 0; j < 5; j++)
            glVertex3fv(&ctrlPoints[j][i][0]);
        glEnd();
    }
}

void bezierSurface()
{
    int n = 5;  // 控制点的数目
    float xarray[5];
    float yarray[5];
    float zarray[5];
    GLfloat ps[11][5][3];
    for (int v = 0; v < 5; v++)
    {
        int u = 0;
        for (double t = 0.0; t <= 1; t += 0.1)
        {
            for (int i = 1; i < n; ++i)
            {
                for (int j = 0; j < n - i; ++j)
                {
                    if (i == 1)
                    // 第一次，使用控制点
                    {
                        xarray[j] = ctrlPoints[j][v][0] * (1 - t) + ctrlPoints[j + 1][v][0] * t;
                        yarray[j] = ctrlPoints[j][v][1] * (1 - t) + ctrlPoints[j + 1][v][1] * t;
                        zarray[j] = ctrlPoints[j][v][2] * (1 - t) + ctrlPoints[j + 1][v][2] * t;
                        continue;
                    }
                    // 非第一次，更新点数组
                    xarray[j] = xarray[j] * (1 - t) + xarray[j + 1] * t;
                    yarray[j] = yarray[j] * (1 - t) + yarray[j + 1] * t;
                    zarray[j] = zarray[j] * (1 - t) + zarray[j + 1] * t;
                }

            }
            ps[u][v][0] = xarray[0];
            ps[u][v][1] = yarray[0];
            ps[u][v][2] = zarray[0];
            u++;
        }
    }

    //glColor3f(255.0, 0.0, 0.0);
    //for (int i = 0; i < 11; i++)
    //{
    //    glBegin(GL_LINE_STRIP);
    //    for (int j = 0; j < 5; j++)
    //        glVertex3fv(&ps[i][j][0]);
    //    glEnd();
    //}

    GLfloat ps1[11][11][3];
    for (int v = 0; v < 11; v++) {
        int u = 0;
        for (double t = 0.0; t <= 1; t += 0.1)
        {
            for (int i = 1; i < n; ++i)
            {
                for (int j = 0; j < n - i; ++j)
                {
                    if (i == 1)
                    {
                        xarray[j] = ps[v][j][0] * (1 - t) + ps[v][j + 1][0] * t;
                        yarray[j] = ps[v][j][1] * (1 - t) + ps[v][j + 1][1] * t;
                        zarray[j] = ps[v][j][2] * (1 - t) + ps[v][j + 1][2] * t;
                        continue;
                    } 
                    xarray[j] = xarray[j] * (1 - t) + xarray[j + 1] * t;
                    yarray[j] = yarray[j] * (1 - t) + yarray[j + 1] * t;
                    zarray[j] = zarray[j] * (1 - t) + zarray[j + 1] * t;
                }

            }
            ps1[v][u][0] = xarray[0];
            ps1[v][u][1] = yarray[0];
            ps1[v][u][2] = zarray[0];
            u++;
        }
    }

    glColor3f(1.0, 1.0, 1.0);
    for (int i = 0; i < 11; i++) 
    {
        glBegin(GL_LINE_STRIP);
        for (int j = 0; j < 11; j++)
            glVertex3fv(&ps1[i][j][0]);
        glEnd();

        glBegin(GL_LINE_STRIP);
        for (int j = 0; j < 11; j++)
            glVertex3fv(&ps1[j][i][0]);
        glEnd();
    }


    glFlush();
}


int main(void)
{
    GLFWwindow* window;

    /* Initialize the library */
    if (!glfwInit())
        return -1;

    /* Create a windowed mode window and its OpenGL context */
    window = glfwCreateWindow(960, 960, "bezierSurface", NULL, NULL);
    if (!window)
    {
        glfwTerminate();
        return -1;
    }

    /* Make the window's context current */
    glfwMakeContextCurrent(window);

    if (glewInit() != GLEW_OK)
        cout << "Error!" << endl;
    
    glRotatef(300.0, 1.0, 0.0, 0.0);

    /* Loop until the user closes the window */
    while (!glfwWindowShouldClose(window))
    {

        /* Render here */
        glClearColor(0.2f, 0.3f, 0.3f, 1.0f);
        glClear(GL_COLOR_BUFFER_BIT);
        plotPoint();
        bezierSurface();
        /* Swap front and back buffers */
        glfwSwapBuffers(window);

        /* Poll for and process events */
        glfwPollEvents();
    }

    glfwTerminate();
    return 0;
}
rd142857
关注
  • 0
    点赞
  • 踩
  • 2
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
四次bezier曲线（直接法、分裂法、几何法）
12-24
四次bezier曲线（直接法、分裂法、几何法）
计算机图形学 第4版.pdf.zip
06-22
计算机图形学，搞图形开发基础，大家有需要自行下载哈
[网盘]计算机图形学 第四版 中文版 完整版.2018_03_17
03-17
计算机图形学 第4版 中文版 [（美）赫恩，（美）巴克，（美）卡里瑟斯著][电子工业出版社][2014.11][678页][13684300].pdf “免责声明：源代码仅仅用于内部学习，不得用于任何商业行为和破坏行为，如有侵权请及时告知，以便我及时下架侵权内容！”
[网盘]计算机图形学 第四版 中文版 完整版.2018_03_19
03-19
计算机图形学 第4版 中文版 [（美）赫恩，（美）巴克，（美）卡里瑟斯著][电子工业出版社][2014.11][678页][13684300].pdf “免责声明：源代码仅仅用于内部学习，不得用于任何商业行为和破坏行为，如有侵权请及时告知，以便我及时下架侵权内容！”
计算机图形学Bezier曲线的生成
03-25
本资源是用动画实现计算机图形学里的一个算法Bezier曲线的生成。里面是详细的代码。
[网盘]计算机图形学 第四版 中文版 完整版
03-21
计算机图形学 第4版 中文版 [（美）赫恩，（美）巴克，（美）卡里瑟斯著][电子工业出版社][2014.11][678页][13684300].pdf “免责声明：源代码仅仅用于内部学习，不得用于任何商业行为和破坏行为，如有侵权请及时告知，以便我及时下架侵权内容！”
glo_纹理和光照模型_计算机图形学_bezier曲面OpenGL_bezier曲面光照_
09-29
在计算机图形学领域，纹理和光照模型是两个关键概念，它们极大地影响了图像的真实感和视觉效果。Bezier曲面作为非线性参数曲线和曲面的一种表示方式，广泛应用于3D建模和渲染中。本篇文章将深入探讨如何使用OpenGL库...
bezier-surface-master.zip_bezier matlab_bezier曲面_bezier曲面matl
09-24
Bezier曲面是一种在计算机图形学中广泛使用的参数化曲面表示方法，特别是在三维建模、渲染和动画领域。它的核心思想是通过一组控制点来定义一个连续光滑的曲面，这些控制点决定了曲面的形状和特征。Matlab作为一种...
计算机图形学(MFC)-Bezier曲面
03-24
在本资源中，我们关注的是Bezier曲面，这是一种在计算机图形学中广泛使用的数学工具，用于创建平滑的曲线和曲面。MFC（Microsoft Foundation Classes）是微软提供的一个C++库，用于构建Windows应用程序，它将Windows...
wyc.rar_WYC_bezier_bezier曲面_图形_曲面
09-22
标题“wyc.rar_WYC_bezier_bezier曲面_图形_曲面”以及描述中的“双三次bezìér曲面（计算机图形学大作业）”暗示了这是一个关于使用bezìér技术来创建和理解曲面的项目或教程。 bezìér曲线由一系列控制点定义...
2次三次4次Bezier曲线演示程序(无源码)
06-17
2次三次4次Bezier曲线演示程序(无源码) QT4.6 VS2008编译
四次Bezier曲线动态生成 VC++代码实现
04-16
Bezier曲线为法国Citeroen汽车公司的 De Cateljau于1959年发明，但当时作为公司 技术机密，直到1975年之后才引起人们注意，1962年法国Renault公司的工程师Bezier独立提出Bezier曲线曲面，于1968年成功用于UNISURF汽车造型系统中，该代码很好模拟显示四次Bezier曲线动态显示非常直观，该代码VC++6实现
数值分析Matlab作二次三次四次曲线拟合
12-08
一种商品的需求量与价格有一定关系。现对一定时期内的商品价格（x）与需求量（y）进行观察，取得一下样本数据，分别作出上述数据点的二次、三次、四次多项式拟合曲线f(x)，绘出图形
计算机图形学—bezier曲面
06-05
计算机图形学实验，用TURBO C 2.0编写，用bezier的定义编写一个4X4阶的曲面。
3次曲面拟合matlab程序
10-22
三次曲面拟合程序，可以改编代码为n次曲面拟合，可以对图像进行低频拟合
上海大学计算机图形学课程项目-Bezier曲面（Opengl版）
03-12
上海大学计算机图形学课程项目 基于Opengl开发 Bezier曲面的离散生成算法(de Casteljau算法) 功能： 数据输入：曲面的控制网格数据通过文件输入 数据存储：设计程序内用于存储控制网格数据的数据结构 数据输出：在...
opengl 四次Bezier曲线（学习笔记-仅供参考）
qq_45228525的博客
06-16 838
#include <GL/glut.h> #include<math.h> void init(void) { glClearColor(0.0, 0.0, 0.0, 0.0); glMatrixMode(GL_PROJECTION); gluOrtho2D(0.0, 500.0, 0.0, 500.0); glMatrixMode(GL_MODELVIEW); } void myDisplay() { float x0 = 50, y0.
贝塞尔曲线二（一次到五次）
最新发布
qq_36246937的博客
03-15 1083
在linux中使用sdl库实现贝塞尔曲线一次到五次
TIT 计算机图形学 实验二 黑白颜色填充双三次bezier曲面片
张时贰的博客
12-22 744
TIT 计算机图形学 实验二 黑白颜色填充双三次bezier曲面片 前言 参考视频计算机图形学全套算法讲解和C++编码实现（共23讲配套源码)，计算机图形学案例视频讲解以及主页相关算法。孔老师是我的代课老师，孔教授有十多年教学经验，视频课很不错，所有的源程序都基于他写的函数，并非小张写的。所有源程序都基于C+编译 参考教材《计算机图形学-理论与实践项目化教程》 孔令德著，大家多多支持哇 一、实验项目要求 1.使用均布网格细分曲面 2.使用黑白棋盘颜色填充细分后的小面 二．理论分析或算法分析 首先绘制双三
计算机图形学：规则曲线曲面生成详解及常用方法
计算机图形学第五章深入探讨了曲线曲面的生成及其在工程中的应用，共36页内容详实丰富。这一章节首先介绍了曲线曲面的表示方法，区分了规则曲线和不规则曲线，规则曲线包括可以通过参数方程直接表示的如二维平面的x=...

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值