Qt OpenGL（07）递归细分四面体法绘制球体

news2024/11/25 22:37:00

文章目录

  • Qt OpenGL通过递归细分逼近球面
  • 思路
  • 下面就是绘制的代码：
    • Widget.cpp
    • 顶点着色器
    • 片段着色器

Qt OpenGL通过递归细分逼近球面

在这里插入图片描述

在OpenGL中绘制球面，不是太简单的事情。因为球面和圆都不是OpenGL所支持的图元，因此我们将通过一种称为递归细分(recursive subdivision)的方法使用三角形来逼近球面, 我们曾经在绘制Sierpinski镂垫时使用过这种技术。

递归细分方法是一种强有力的技术, 利用它可以以任意精度逼近曲线和曲面。

我们将讨论如何利用正四面体来近似球面。

我们的出发点是一个正四面体，正四面体是由4个等边三角形围成的，由4个顶点确定。这4个顶点都位于以原点为中心的单位球面上。

以一个三角形为例，将其分为4个小三角形

在这里插入图片描述

这一过程和 Siepinski 镂垫极其类似，但为啥最终会显示为球面？最重要的步骤就是顶点的归一化。归一化后的坐标就是贴近球面的坐标。

在这里插入图片描述 在这里插入图片描述

      没有归一化               归一化的效果

思路

引自《交互式计算机图形学 基于OpenGL着色器的自顶向下方法 第6版_中文版》一书，略作修改

  1. 正四面体坐标

    QVector3D v[4] = {
    QVector3D(  0.0f,       0.0f,       1.0f      ),
    QVector3D(  0.0f,       0.942809f, -0.333333f ),
    QVector3D( -0.816497f, -0.471405f, -0.333333f ),
    QVector3D(  0.816497f, -0.471405f, -0.333333f )
};

  2. 定义三角形的顶点数据的容器

    std::vector<QVector3D> vdata;

  3. 保存三角形的函数

    void triangle(QVector3D a, QVector3D b, QVector3D c){
    vdata.push_back (a);
    vdata.push_back (b);
    vdata.push_back (c);
}

  4. 细分方法

    void divide_triangle(QVector3D a, QVector3D b, QVector3D c, int n ){
    QVector3D v1, v2, v3;
    if( n > 0 ) {
        // 如果使用归一化，顶点和的结果是否除以 2 都不影响结果  
        v1 = ( a + b )/2;
        v2 = ( a + c )/2;
        v3 = ( b + c )/2;

        // 归一化 的坐标贴近球面。但如果不进行归一化，显示的就是一个细分的四面体
        v1.normalize ();
        v2.normalize ();
        v3.normalize ();

        divide_triangle(a, v2, v1, n-1);
        divide_triangle(c, v3, v2, n-1);
        divide_triangle(b, v1, v3, n-1);
        divide_triangle(v1,v2, v3, n-1);
    }
    else
        triangle (a,b,c);
}

  5. 定义一个递归函数

    void tetrahedron( int n ){
//    divide_triangle( v[0], v[2], v[1], n);
//    divide_triangle( v[0], v[3], v[2], n);
//    divide_triangle( v[0], v[1], v[3], n);
//    divide_triangle( v[1], v[2], v[3], n);
    divide_triangle( v[0], v[1], v[2], n);
    divide_triangle( v[3], v[2], v[1], n);
    divide_triangle( v[0], v[3], v[1], n);
    divide_triangle( v[0], v[2], v[3], n);
}

下面就是绘制的代码：

其它文件略过, 可参考前面的文章。

Widget.cpp

#include "Widget.h"
#include "Camera.hpp"
#include <QApplication>
#include <QMouseEvent>
#include <QThread>
#include <vector>
int rotateAngle = 0;
#define qRandom   QRandomGenerator::global ()
#define qout if( 1 ) qDebug() << __FILE__ << __LINE__ << ": "

QVector3D v[4] = {
    QVector3D(  0.0f,       0.0f,       1.0f      ),
    QVector3D(  0.0f,       0.942809f, -0.333333f ),
    QVector3D( -0.816497f, -0.471405f, -0.333333f ),
    QVector3D(  0.816497f, -0.471405f, -0.333333f )
};

std::vector<QVector3D> vdata;

void triangle(QVector3D a, QVector3D b, QVector3D c){
    vdata.push_back (a);
    vdata.push_back (b);
    vdata.push_back (c);
}

void divide_triangle(QVector3D a, QVector3D b, QVector3D c, int n ){
    QVector3D v1, v2, v3;
    if( n > 0 ) {
        // 如果使用归一化，顶点和的结果是否除以 2 都不影响结果
        v1 = ( a + b )/2;
        v2 = ( a + c )/2;
        v3 = ( b + c )/2;

        // 归一化 的坐标贴近球迷。但如果不进行归一化，显示的就是一个细分的四面体
        v1.normalize ();
        v2.normalize ();
        v3.normalize ();

        divide_triangle(a, v2, v1, n-1);
        divide_triangle(c, v3, v2, n-1);
        divide_triangle(b, v1, v3, n-1);
        divide_triangle(v1,v2, v3, n-1);
    }
    else
        triangle (a,b,c);
}

void tetrahedron( int n ){
    divide_triangle( v[0], v[1], v[2], n);
    divide_triangle( v[3], v[2], v[1], n);
    divide_triangle( v[0], v[3], v[1], n);
    divide_triangle( v[0], v[2], v[3], n);
}

Widget::Widget(QWidget *parent)
    : QOpenGLWidget(parent)
{
    setWindowTitle ("08_Sphere");
    resize (200,200);
    ratio = qApp->devicePixelRatio ();

    tetrahedron(4);
    qout << sizeof(QVector3D)*vdata.size () ;
}

Widget::~Widget()
{
    makeCurrent ();
    glDeleteBuffers (1,&VBO);
    glDeleteVertexArrays (1,&VAO);
    doneCurrent ();
}

void Widget::initializeGL()
{
    initializeOpenGLFunctions ();
    // ---------------------------------

    glGenBuffers (1,&VBO);
    glBindBuffer (GL_ARRAY_BUFFER,VBO);

    glBufferData (GL_ARRAY_BUFFER,
                  sizeof(QVector3D)*vdata.size (),
                  NULL,
                  GL_STATIC_DRAW);


    glBufferSubData (GL_ARRAY_BUFFER,
                     0,
                     sizeof(QVector3D)*vdata.size (),
                     vdata.data ());

    // ---------------------------------
    glGenVertexArrays (1,&VAO);
    glBindVertexArray(VAO);
    glVertexAttribPointer(0,    // vao 索引，对应顶点着色器中的变量的位置 ： “layout (location = 0) in vec3 aPos;”
                          3,    // 变量中元素的个数， 比如 aPos变量 有 3 个数据【x,y,z】组成
                          GL_FLOAT, // 类型
                          GL_FALSE, // 标准化，是否在 [-1,1] 之间
                          sizeof(QVector3D),  // 步长，表示下个元组的首元素 和 该元组首元素之间的大小，因为顶点数据中 可能还夹杂颜色、法向量等数据
                          (void*)0 );   // 变量的偏移量，在多个变量混合时指定变量的偏移
    glEnableVertexAttribArray(0); // 使用 location = 0 的索引

    // ---------------------------------

    QString filename = ":/shader";
    shaderProgram.addShaderFromSourceFile (QOpenGLShader::Vertex,   filename+".vert");
    shaderProgram.addShaderFromSourceFile (QOpenGLShader::Fragment, filename+".frag");
    shaderProgram.link ();

    glBindBuffer (GL_ARRAY_BUFFER,0);

    glPolygonMode (GL_FRONT_AND_BACK,GL_LINE);

    // 深度测试
    glEnable(GL_DEPTH_TEST);
}

void Widget::paintGL()
{
    glClearColor(0.2f, 0.3f, 0.3f, 1.0f);   // 设置背景色
    glClear(GL_COLOR_BUFFER_BIT | GL_DEPTH_BUFFER_BIT);

    glViewport ((max-min)/2 ,0,min,min);

    shaderProgram.bind ();
    glBindVertexArray(VAO);

    QMatrix4x4 model;
//    model.rotate ( 15, 1.0f, 0.0f, 0.0f);
//    model.rotate ( rotateAngle-15, 0.0f, 1.0f, 0.0f);
    QMatrix4x4 view;
    QMatrix4x4 projection;
    view = m_Camera.GetViewMatrix ();
    projection.perspective (m_Camera.Zoom, 1.0f ,0.01f, -100.0f);
//    projection.ortho (-1,1,-1,1 ,-1, 100);

    model.rotate ( rotateAngle-15, 0.0f, 1.0f, 0.0f);

    shaderProgram.setUniformValue ("model",model);
    shaderProgram.setUniformValue ("view",view);
    shaderProgram.setUniformValue ("projection",projection);

    shaderProgram.setUniformValue ("u_color",1.0f, 0.0f, 0.0f, 1.0f); // 红色

    glDrawArrays (GL_TRIANGLES, 0, (int) vdata.size () );

//    qout << this->windowState ();
//    if( this->windowState () != Qt::WindowMinimized && 1 ){
    if( 1 ){
        QThread::currentThread ()->msleep (50);
        rotateAngle += 1 ;
        update ();
    }
}

void Widget::resizeGL(int w, int h)
{
    qout << "resizeGL";
    min = std::min (w,h) * ratio;
    max = std::max (w,h) * ratio;
}

顶点着色器

#version 330 core
layout (location = 0) in vec3 aPos;
uniform mat4 model;
uniform mat4 view;
uniform mat4 projection;
out float z;
out vec3 color;
void main()
{
   gl_Position = projection * view * model * vec4(aPos.x, aPos.y, aPos.z, 1.0);
   color = vec3(sqrt(aPos.x*aPos.x),sqrt(aPos.y * aPos.y),sqrt(aPos.z * aPos.z));
   z = aPos.z;
};

这种着色方法参考了此博客 ： OpenGL递归细分四面体法绘制球体_走肖暂时无法接通的博客-CSDN博客

片段着色器

#version 330 core
uniform vec4 u_color;
in float z;
in vec3 color;
void main()
{
    gl_FragColor = vec4( color *( 0.3 + z ),1.0 );  // 使用坐标转换过来的颜色
//    gl_FragColor = u_color;
}

本文来自互联网用户投稿，该文观点仅代表作者本人，不代表本站立场。本站仅提供信息存储空间服务，不拥有所有权，不承担相关法律责任。如若转载，请注明出处：http://www.coloradmin.cn/o/154083.html

如若内容造成侵权/违法违规/事实不符，请联系多彩编程网进行投诉反馈，一经查实，立即删除！

相关文章

飞宇科技在创业板IPO终止：长江证券撤回保荐，吴玉飞为实控人

飞宇科技在创业板IPO终止：长江证券撤回保荐，吴玉飞为实控人

2023年1月10日&#xff0c;深圳证券交易所披露的信息显示&#xff0c;江苏飞宇医药科技股份有限公司&#xff08;下称“飞宇科技”&#xff09;的保荐机构&#xff08;长江证券&#xff09;提交了撤回飞宇科技首次公开发行股票并在创业板上市申请的申请因此&#xff0c;深圳证券…
阅读更多...
BOM（二）： 元素偏移量offset 、元素可视区client 、元素滚动 scroll、动画

BOM（二）： 元素偏移量offset 、元素可视区client 、元素滚动 scroll、动画

PC端网页特效元素偏移量offset元素可视区client元素滚动 scroll 系列动画元素偏移量offset 1.offset概述 offset 系列相关属性可以动态的得到该元素的位置&#xff08;偏移&#xff09;、大小等。 获得元素距离带有定位父元素的位置获得元素自身的大小&#xff08;宽度高度&…
阅读更多...
Python 字典(Dictionary)操作详解

Python 字典(Dictionary)操作详解

这篇文章主要介绍了Python 字典(Dictionary)的详细操作方法&#xff0c;需要的朋友可以参考下。 Python字典是另一种可变容器模型&#xff0c;且可存储任意类型对象&#xff0c;如字符串、数字、元组等其他容器模型。 一、创建字典 字典由键和对应值成对组成。字典也被称作关…
阅读更多...
微信小程序云开发之新闻博客社区项目debug后的项目代码

微信小程序云开发之新闻博客社区项目debug后的项目代码

大家好&#xff0c;我是csdn的博主&#xff1a;lqj_本人 这是我的个人博客主页&#xff1a;lqj_本人的博客_CSDN博客-微信小程序,html特效,vue2基础领域博主 本次文章主要时为我最近在哔哩哔哩上的新发布的视频做一个映射&#xff01; 哔哩哔哩欢迎关注&#xff1a;小淼前端 哔…
阅读更多...
【Python从入门到进阶】3、运行python代码

【Python从入门到进阶】3、运行python代码

接上篇《2、Python环境的安装》 上一篇我们学习了如何下载安装Python编程环境&#xff0c;以及如何使用pip管理Python包。本篇我们来学习如何使用终端和执行文件运行python代码。 一、终端运行 我们可以在命令行终端进入python解释程序&#xff0c;直接输入python代码&#x…
阅读更多...
Spring Cloud OpenFeign @SpringQueryMap注解

Spring Cloud OpenFeign @SpringQueryMap注解

概述 我们在使用GetMapping方法是&#xff0c;对于多个参数都是通过RequestParam参数一个一个接取的&#xff0c;多个参数时会比较麻烦&#xff0c;能像RequestBody那样接取PostMapping的参数么&#xff1f;答案是可以的。 SpringQueryMap介绍 SpringQueryMap 注解是 spring…
阅读更多...
关于数据管理、数据治理和数据资产，你搞得清楚吗？

关于数据管理、数据治理和数据资产，你搞得清楚吗？

全球数字化的趋势愈演愈烈&#xff0c;我国信息化规模日渐增大&#xff0c;我们也每天都在接收大量的数据&#xff0c;要想深度挖掘数据的价值&#xff0c;首先从搞清楚数据的概念入手。今天我们就来聊一聊数据管理、数据治理和数据资产的概念和它们之间的相互关系。 数据管理…
阅读更多...
从低代码来看，何谓客户体验自动化

从低代码来看，何谓客户体验自动化

在当今竞争激烈的大环境下&#xff0c;拥有出色的产品或服务并不意味着100%成功&#xff0c;反而是客户在购买流程中的体验变得比以往任何时候都更加重要。 随着企业业务的发展&#xff0c;提供个性化的客户体验变得愈加困难&#xff0c;许多企业因无法提供令人满意的客户体验而…
阅读更多...
JavaScript 与 C++ 的差异

JavaScript 与 C++ 的差异

在正式开始向各位前端开发者介绍 C 语言之前&#xff0c;我们先讨论一下 C 语言与 JavaScript 语言的差异&#xff08;为了约束讨论的范围&#xff0c;这里就不提 HTML 和 CSS 了&#xff09;。 C 语言于 1979 年由贝尔实验室的 Bjarne Stroustrup&#xff08;本贾尼斯特劳斯特…
阅读更多...
【web安全】——HTTP请求头注入

【web安全】——HTTP请求头注入

作者名&#xff1a;Demo不是emo主页面链接&#xff1a; 主页传送门创作初心&#xff1a; 舞台再大&#xff0c;你不上台&#xff0c;永远是观众&#xff0c;没人会关心你努不努力&#xff0c;摔的痛不痛&#xff0c;他们只会看你最后站在什么位置&#xff0c;然后羡慕或鄙夷座右…
阅读更多...
uni-app开发常用点

uni-app开发常用点

前言 应该会不定时更新&#xff0c;主要记录自己在uni-app开发过程中经常用到的东西。 组件库推荐&#xff1a; 1、https://v1.uviewui.com/ 我们老项目用到是这个的1.x版本&#xff0c;新版本没用过 2、https://vant-contrib.gitee.io/vant/v2/#/zh-CN/ 3、https://uniapp.d…
阅读更多...
智能在线客服系统源码 国际版多语言多商户智能机器人源码

智能在线客服系统源码 国际版多语言多商户智能机器人源码

一套智能在线客服系统源码 多商户网页客服系统源码 支持二十种国际语言 带机器人自动回复。 框架&#xff1a;Thinkphp5workerman&#xff0c; 环境&#xff1a;nginxphp7.3mysql5.6 支持H5公众号APP小程序 私信了解更多源码内容&#xff01; 系统功能特点&#xff1a; 1、…
阅读更多...
k8s创建数据库

k8s创建数据库

参考配置1 [rootk8s-master1 ~]# cat deployment.yaml apiVersion: apps/v1 kind: Deployment metadata:labels:app: mysql #为该Deployment设置key为app&#xff0c;value为mysql的标签name: mysqlnamespace: test spec:replicas: 1 #副本数量selector: #标签选择器&#xf…
阅读更多...
深挖你拥有的东西的背后的价值到底是什么，价值转换模型分析。

深挖你拥有的东西的背后的价值到底是什么，价值转换模型分析。

深挖你拥有的东西的背后的价值到底是什么。 我们先把它转换成为观众观看这条视频的动机&#xff0c;甚至是观众关注账号的动机。 例如你能帮助别人赚钱&#xff0c;你能帮助别人变美。你提供的价值是成长的价值&#xff0c;你提供的价值是精神追求的价值&#xff0c;是生活理念…
阅读更多...
【ROS】—— ROS常用组件—rosbag与rqt工具箱(十二)

【ROS】—— ROS常用组件—rosbag与rqt工具箱(十二)

文章目录前言1. rosbag1.1 rosbag使用_命令行2. rosbag使用_编码2.1 C实现2.1.1 写bag2.1.2 读bag2.2 python实现2.2.1 写bag2.2.2 读bag3. rqt工具箱3.1 rqt安装启动与基本使用3.2 rqt常用插件:rqt_graph3.3 rqt常用插件:rqt_console3.4 rqt常用插件:rqt_plot3.5 rqt常用插件:…
阅读更多...
python算法与数据结构2-栈、队列、排序算法

python算法与数据结构2-栈、队列、排序算法

目录1、栈1.1 栈的介绍1.2 栈的代码实现2、队列2.1 队列的介绍2.2 队列的代码实现3、双端队列4、数据结构与算法_排序算法4.1 排序算法的稳定性4.2 冒泡排序4.3 选择排序4.4 插入排序4.5 快速排序1、栈 1.1 栈的介绍 栈&#xff1a;运算受限的线性表&#xff0c;其限制是仅允…
阅读更多...
黑盒测试用例设计 - 功能图法

黑盒测试用例设计 - 功能图法

目录功能图法原理介绍功能图法步骤案例&#xff1a;以QQ登录界面为例功能图法原理介绍 功能图方法其实是一种灰盒测试&#xff08;因其兼有黑盒和白盒测试&#xff09;用例设计方法&#xff1b;通常情况一个程序的功能说明通常由动态说明和静态说明组成&#xff0c;动态说明描…
阅读更多...
UG/NX 二次开发(C#)自动出2D零件图思路

UG/NX 二次开发(C#)自动出2D零件图思路

一、前言 项目需要自动出2D零件图&#xff0c;可是我之前没做过这方面的内容&#xff0c;没有一点思路。然后我就做了下面几件事&#xff1a; 1、百度、google翻了一遍&#xff0c;搜索关键字“UG二次开发自动出图” 2、csdn.com 和 cnblogs.com翻了一遍 3、平时逛的qq技术…
阅读更多...
数字图像处理（第四版)-冈萨雷斯-学习过程的笔记

数字图像处理（第四版)-冈萨雷斯-学习过程的笔记

作者介绍1.绪论自己读书存在一个问题&#xff0c;书太厚&#xff0c;重点难以把握&#xff0c;对如此经典的书籍&#xff0c;希望自己的学习历程和重点记录下来。 目前在持续更新和学习中&#xff0c;觉得有帮助的话可以先收藏和关注我博客的内容--更新于2023/1/102.数字图像基…
阅读更多...
C语言进阶——自定义类型——结构体

C语言进阶——自定义类型——结构体

目录 一. 结构体类型的声明 结构体的声明 特殊的声明 二. 结构的自引用 三. 结构体变量的定义和初始化 定义 初始化 四. 结构体内存对齐 未完待续 一. 结构体类型的声明 结构体的声明 结构是一些值的集合&#xff0c;这些值称为成员变量。结构的每个成员可以…
阅读更多...
推荐文章
最新文章

Copyright @ 2022~2023