VTK 的可视化方法:流线、流管、流面、流带

embedded/2024/9/25 15:32:59/

VTK 的可视化方法:流线、流管、流面、流带

  • VTK 的可视化方法:流线、流管、流面、流带
    • 三种相似的可视化方法
    • 流线生成使用的类
    • 实例:单一流线
    • 实例:流管
    • 实例:多条流线
    • 实例:流面
    • 实例:流带
    • 完整代码

VTK__1">VTK 的可视化方法:流线、流管、流面、流带

2维流线示例:

在这里插入图片描述

本文章主要讲解3维的流线、流管、流面、流带的构造方法。

三种相似的可视化方法

  1. 流线(Streamlines):每个点速度切线方向连成的线。
  2. 迹线(Pathlines):粒子实际的轨迹线。
  3. 脉线(Streaklines):连续时刻出发的粒子在某一时刻的连线。

流线生成使用的类

  • vtkRungeKutta4:四阶龙格库塔 (Runge-Kutta) 求解微分。
  • vtkStreamTracer:通过整合矢量场生成流线。

实例:单一流线

完整代码:

#include "VTKStreamline.h"#include <vtkConeSource.h>
#include <vtkMultiBlockPLOT3DReader.h>
#include <vtkDataSet.h>
#include <vtkMultiBlockDataSet.h>
#include <vtkRungeKutta4.h>
#include <vtkStreamTracer.h>
#include <vtkDataArray.h>
#include <vtkPointData.h>
#include <vtkShrinkPolyData.h>
#include <vtkStructuredGridGeometryFilter.h>
#include <vtkStructuredGridOutlineFilter.h>
#include <vtkContourFilter.h>
#include <vtkPolyDataMapper.h>
#include <vtkActor.h>
#include <vtkRenderer.h>
#include <vtkRenderWindow.h>VTKStreamline::VTKStreamline(QWidget* parent): QMainWindow(parent)
{ui.setupUi(this);_pVTKWidget = new QVTKOpenGLNativeWidget();this->setCentralWidget(_pVTKWidget);// this->showMaximized();// 1. generate data// vtkSmartPointer<vtkConeSource> cone = vtkSmartPointer<vtkConeSource>::New();// or, read data// vtkMultiBlockPLOT3DReader 是一个读取器对象,用于读取 PLOT3D 格式的文件并在输出时生成结构化网格vtkSmartPointer<vtkMultiBlockPLOT3DReader> plot3dReader = vtkSmartPointer<vtkMultiBlockPLOT3DReader>::New();plot3dReader->SetXYZFileName("combxyz.bin");plot3dReader->SetQFileName("combq.bin");plot3dReader->SetScalarFunctionNumber(100);plot3dReader->SetVectorFunctionNumber(202);qDebug() << plot3dReader->GetOutput()->GetNumberOfBlocks(); // 0// 反向更新管线plot3dReader->Update();qDebug() << plot3dReader->GetOutput()->GetNumberOfBlocks(); // 1vtkDataSet* plot3dOutput = (vtkDataSet*)(plot3dReader->GetOutput()->GetBlock(0));// 四阶龙格库塔 (Runge-Kutta) 求解微分vtkSmartPointer<vtkRungeKutta4> integ = vtkSmartPointer<vtkRungeKutta4>::New();// 通过整合矢量场生成流线vtkSmartPointer<vtkStreamTracer> streamer = vtkSmartPointer<vtkStreamTracer>::New();streamer->SetIntegrator(integ);streamer->SetInputData(plot3dOutput);streamer->SetStartPosition(15, 5, 32);streamer->SetMaximumPropagation(100);streamer->SetInitialIntegreationStep(0.1);streamer->SetIntegreationDirectionToBackward();// 2. filter// 产生结构化栅格边界的一个线轮廓vtkSmartPointer<vtkStructuredGridOutlineFilter> outline = vtkSmartPointer<vtkStructuredGridOutlineFilter>::New();outline->SetInputData(plot3dOutput);// 3. mappervtkSmartPointer<vtkPolyDataMapper> outlineMapper = vtkSmartPointer<vtkPolyDataMapper>::New();vtkSmartPointer<vtkPolyDataMapper> singleMapper = vtkSmartPointer<vtkPolyDataMapper>::New();singleMapper->SetScalarRange(plot3dOutput->GetPointData()->GetScalars()->GetRange());// 4. actorvtkSmartPointer<vtkActor> outlineActor = vtkSmartPointer<vtkActor>::New();vtkSmartPointer<vtkActor> singleActor = vtkSmartPointer<vtkActor>::New();// 5. renderervtkSmartPointer<vtkRenderer> renderer = vtkSmartPointer<vtkRenderer>::New();renderer->SetBackground(0.3, 0.6, 0.3); // Background Color: Green// 6. connectoutlineMapper->SetInputConnection(outline->GetOutputPort());singleMapper->SetInputConnection(streamer->GetOutputPort());outlineActor->SetMapper(outlineMapper);singleActor->SetMapper(singleMapper);renderer->AddActor(outlineActor);renderer->AddActor(singleActor);this->_pVTKWidget->renderWindow()->AddRenderer(renderer);this->_pVTKWidget->renderWindow()->Render();
}VTKStreamline::~VTKStreamline()
{}

运行效果:

在这里插入图片描述

实例:流管

我们再加上一个过滤器,把流线变成流管:

	vtkSmartPointer<vtkTubeFilter> streamTube = vtkSmartPointer<vtkTubeFilter>::New();streamTube->SetInputConnection(streamer->GetOutputPort());streamTube->SetRadius(0.06);streamTube->SetNumberOfSides(12);

流管实际上是用一个圆柱面包裹住流线,流线依旧存在。这样的显示效果会更好一点:

在这里插入图片描述

实例:多条流线

在之前的代码中,新增:

	vtkSmartPointer<vtkLineSource> seeds = vtkSmartPointer<vtkLineSource>::New();// 设置线段的端点seeds->SetPoint1(15, -5, 32);seeds->SetPoint2(15, 5, 32);seeds->SetResolution(21);	vtkSmartPointer<vtkStreamTracer> streamer2 = vtkSmartPointer<vtkStreamTracer>::New();streamer2->SetIntegrator(integ);streamer2->SetInputData(plot3dOutput);// streamer2->SetStartPosition(15, 5, 32);streamer2->SetMaximumPropagation(100);streamer2->SetInitialIntegreationStep(0.1);streamer2->SetIntegreationDirectionToBackward();streamer2->SetSourceConnection(seeds->GetOutputPort());vtkSmartPointer<vtkPolyDataMapper> multipleMapper = vtkSmartPointer<vtkPolyDataMapper>::New();multipleMapper->SetScalarRange(plot3dOutput->GetPointData()->GetScalars()->GetRange());vtkSmartPointer<vtkActor> multipleActor = vtkSmartPointer<vtkActor>::New();multipleMapper->SetInputConnection(streamer2->GetOutputPort());multipleActor->SetMapper(multipleMapper);renderer->AddActor(multipleActor);

在这里插入图片描述

实例:流面

新增一个过滤器 vtkRuledSurfaceFilter,把多条流线合并成一个流面:

	vtkSmartPointer<vtkRuledSurfaceFilter> scalarSurface = vtkSmartPointer<vtkRuledSurfaceFilter>::New();scalarSurface->SetInputConnection(streamer2->GetOutputPort());// 设置生成方法scalarSurface->SetRuledModeToPointWalk();// multipleMapper->SetInputConnection(streamer2->GetOutputPort());multipleMapper->SetInputConnection(scalarSurface->GetOutputPort());

运行结果:

在这里插入图片描述

实例:流带

流带其实是按每条流线拓展而成的一条条带状的面,比起流面,更能便于展示走势。

我们只需要在前面的代码中新增一行代码:

	vtkSmartPointer<vtkRuledSurfaceFilter> scalarSurface = vtkSmartPointer<vtkRuledSurfaceFilter>::New();scalarSurface->SetInputConnection(streamer2->GetOutputPort());// 设置生成方法scalarSurface->SetRuledModeToPointWalk();scalarSurface->SetOnRatio(2); // 新增代码

运行结果:

在这里插入图片描述

完整代码

#include "VTKStreamline.h"#include <vtkConeSource.h>
#include <vtkLineSource.h>
#include <vtkMultiBlockPLOT3DReader.h>
#include <vtkDataSet.h>
#include <vtkMultiBlockDataSet.h>
#include <vtkRungeKutta4.h>
#include <vtkStreamTracer.h>
#include <vtkTubeFilter.h>
#include <vtkDataArray.h>
#include <vtkPointData.h>
#include <vtkShrinkPolyData.h>
#include <vtkStructuredGridGeometryFilter.h>
#include <vtkStructuredGridOutlineFilter.h>
#include <vtkContourFilter.h>
#include <vtkRuledSurfaceFilter.h>
#include <vtkPolyDataMapper.h>
#include <vtkActor.h>
#include <vtkRenderer.h>
#include <vtkRenderWindow.h>VTKStreamline::VTKStreamline(QWidget* parent): QMainWindow(parent)
{ui.setupUi(this);_pVTKWidget = new QVTKOpenGLNativeWidget();this->setCentralWidget(_pVTKWidget);// this->showMaximized();// 1. generate datavtkSmartPointer<vtkLineSource> seeds = vtkSmartPointer<vtkLineSource>::New();// 设置线段的端点seeds->SetPoint1(15, -5, 32);seeds->SetPoint2(15, 5, 32);seeds->SetResolution(21);// or, read data// vtkMultiBlockPLOT3DReader 是一个读取器对象,用于读取 PLOT3D 格式的文件并在输出时生成结构化网格vtkSmartPointer<vtkMultiBlockPLOT3DReader> plot3dReader = vtkSmartPointer<vtkMultiBlockPLOT3DReader>::New();plot3dReader->SetXYZFileName("combxyz.bin");plot3dReader->SetQFileName("combq.bin");plot3dReader->SetScalarFunctionNumber(100);plot3dReader->SetVectorFunctionNumber(202);qDebug() << plot3dReader->GetOutput()->GetNumberOfBlocks(); // 0// 反向更新管线plot3dReader->Update();qDebug() << plot3dReader->GetOutput()->GetNumberOfBlocks(); // 1vtkDataSet* plot3dOutput = (vtkDataSet*)(plot3dReader->GetOutput()->GetBlock(0));// 四阶龙格库塔 (Runge-Kutta) 求解微分vtkSmartPointer<vtkRungeKutta4> integ = vtkSmartPointer<vtkRungeKutta4>::New();// 通过整合矢量场生成流线vtkSmartPointer<vtkStreamTracer> streamer = vtkSmartPointer<vtkStreamTracer>::New();streamer->SetIntegrator(integ);streamer->SetInputData(plot3dOutput);streamer->SetStartPosition(15, 5, 32);streamer->SetMaximumPropagation(100);streamer->SetInitialIntegreationStep(0.1);streamer->SetIntegreationDirectionToBackward();vtkSmartPointer<vtkStreamTracer> streamer2 = vtkSmartPointer<vtkStreamTracer>::New();streamer2->SetIntegrator(integ);streamer2->SetInputData(plot3dOutput);// streamer2->SetStartPosition(15, 5, 32);streamer2->SetMaximumPropagation(100);streamer2->SetInitialIntegreationStep(0.1);streamer2->SetIntegreationDirectionToBackward();streamer2->SetSourceConnection(seeds->GetOutputPort());// 2. filter// 产生结构化栅格边界的一个线轮廓vtkSmartPointer<vtkStructuredGridOutlineFilter> outline = vtkSmartPointer<vtkStructuredGridOutlineFilter>::New();outline->SetInputData(plot3dOutput);vtkSmartPointer<vtkTubeFilter> streamTube = vtkSmartPointer<vtkTubeFilter>::New();streamTube->SetInputConnection(streamer->GetOutputPort());streamTube->SetRadius(0.06);streamTube->SetNumberOfSides(12);vtkSmartPointer<vtkRuledSurfaceFilter> scalarSurface = vtkSmartPointer<vtkRuledSurfaceFilter>::New();scalarSurface->SetInputConnection(streamer2->GetOutputPort());// 设置生成方法scalarSurface->SetRuledModeToPointWalk();scalarSurface->SetOnRatio(2);// 3. mappervtkSmartPointer<vtkPolyDataMapper> outlineMapper = vtkSmartPointer<vtkPolyDataMapper>::New();vtkSmartPointer<vtkPolyDataMapper> singleMapper = vtkSmartPointer<vtkPolyDataMapper>::New();singleMapper->SetScalarRange(plot3dOutput->GetPointData()->GetScalars()->GetRange());vtkSmartPointer<vtkPolyDataMapper> multipleMapper = vtkSmartPointer<vtkPolyDataMapper>::New();multipleMapper->SetScalarRange(plot3dOutput->GetPointData()->GetScalars()->GetRange());// 4. actorvtkSmartPointer<vtkActor> outlineActor = vtkSmartPointer<vtkActor>::New();vtkSmartPointer<vtkActor> singleActor = vtkSmartPointer<vtkActor>::New();vtkSmartPointer<vtkActor> multipleActor = vtkSmartPointer<vtkActor>::New();// 5. renderervtkSmartPointer<vtkRenderer> renderer = vtkSmartPointer<vtkRenderer>::New();renderer->SetBackground(0.3, 0.6, 0.3); // Background Color: Green// 6. connectoutlineMapper->SetInputConnection(outline->GetOutputPort());singleMapper->SetInputConnection(streamTube->GetOutputPort());// multipleMapper->SetInputConnection(streamer2->GetOutputPort());multipleMapper->SetInputConnection(scalarSurface->GetOutputPort());outlineActor->SetMapper(outlineMapper);singleActor->SetMapper(singleMapper);multipleActor->SetMapper(multipleMapper);renderer->AddActor(outlineActor);renderer->AddActor(singleActor);renderer->AddActor(multipleActor);this->_pVTKWidget->renderWindow()->AddRenderer(renderer);this->_pVTKWidget->renderWindow()->Render();
}VTKStreamline::~VTKStreamline()
{}

http://www.ppmy.cn/embedded/29929.html

相关文章

C语言 循环语句 (2) do-while 循环语句

上文 C语言 循环语句 (1) 讲述循环概念演示while语句 我们讲述了循环的基本概念 且 演示了 while 循环语句 那么 本文 我们继续 我们来看 第二个循环语句 do-while 语句 语句结构如下图 先是 关键字 do 然后 执行循环结构 然后 while 表达式判断 执行流程图如下图 它的结…

速盾:cdn真的抗打吗?

CDN&#xff08;Content Delivery Network&#xff09;是一种分布式网络架构&#xff0c;旨在提供高速、可靠的内容传输服务。它通过将内容存储在位于全球各地的服务器上&#xff0c;使用户可以从最近的服务器获取内容&#xff0c;从而提高访问速度和性能。然而&#xff0c;尽管…

Node.js -- MongoDB

文章目录 1. 相关介绍2. 核心概念3. 命令行交互3.1数据库命令3.2 集合命令3.3 文档命令 4. 数据库应用场景4.1 新增4.2 删除4.3 更新4.4 查询 5. 图形化工具Robo 3T 1. 相关介绍 一、简介 Mongodb是什么 MongoDB是一个基于分布式文件存储的数据库&#xff0c;官方地址https://…

Nuxt3路由跳转

学习这篇文章之前&#xff0c;需要具备使用Nuxt3框架搭建基础项目 ,如果有这个基础可以跳过.直接看下面的声明式路由跳转 声明式路由跳转 创建两个基础文件 ~/pages/index.vue 和 ~/pages/About.vue ⚠️一定要是index.vue不能够大写Index.vue&#xff0c; 否则报错. <!--…

SQL注入漏洞--报错/union/布尔盲注/时间盲注

之前介绍了数据库的基本操作&#xff0c;今天这篇文章就来实操SQL注入。 阅读本文前可以先看一下基本操作&#xff0c;有助于更好理解本文。。。 https://blog.csdn.net/weixin_60885144/article/details/138356410?spm1001.2014.3001.5502 what SQL---结构化查询语言---S…

RISC-V异常处理相关内容

异常处理相关内容 异常处理相关内容异常处理准备工作异常处理函数Opensbi系统调用的注册异常处理相关内容 异常处理准备工作 这里需要特殊强调的是异常处理构建的相关内容: 这里会将a4寄存器中的值存储到CSR_MTVEC这个状态寄存器,也就是异常处理程序的的入口;如果遇到异常…

【StreetGen】城市级交通路网生成算法

街道面积广阔、多样化&#xff0c;用于多种&#xff08;可能相互冲突的&#xff09;交通方式以及社会和文化活动。 适当的规划至关重要&#xff0c;并且需要数据。 手动制作代表街道的数据&#xff08;街道重建&#xff09;容易出错且耗时。 由于细节的多样性、大小和规模&…

在远程服务器上安装anaconda以及配置pytorch虚拟环境

目录 第一步&#xff1a;官网或者清华源下载Anaconda。 第二步&#xff1a;创建虚拟环境。 第三步&#xff1a;在服务器终端输入nvidia-smi查看服务器信息。 第四步&#xff1a;在pytorch官网找到对应版本cuda的命令。 第一步&#xff1a;官网或者清华源下载Anaconda。 官网…