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

2024-05-03 09:52

本文主要是介绍VTK 的可视化方法:流线、流管、流面、流带,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

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

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

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()
{}

这篇关于VTK 的可视化方法:流线、流管、流面、流带的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/956462

相关文章

SpringBoot中使用Flux实现流式返回的方法小结

《SpringBoot中使用Flux实现流式返回的方法小结》文章介绍流式返回(StreamingResponse)在SpringBoot中通过Flux实现,优势包括提升用户体验、降低内存消耗、支持长连... 目录背景流式返回的核心概念与优势1. 提升用户体验2. 降低内存消耗3. 支持长连接与实时通信在Sp

Conda虚拟环境的复制和迁移的四种方法实现

《Conda虚拟环境的复制和迁移的四种方法实现》本文主要介绍了Conda虚拟环境的复制和迁移的四种方法实现,包括requirements.txt,environment.yml,conda-pack,... 目录在本机复制Conda虚拟环境相同操作系统之间复制环境方法一:requirements.txt方法

Nginx 重写与重定向配置方法

《Nginx重写与重定向配置方法》Nginx重写与重定向区别:重写修改路径(客户端无感知),重定向跳转新URL(客户端感知),try_files检查文件/目录存在性,return301直接返回永久重... 目录一.try_files指令二.return指令三.rewrite指令区分重写与重定向重写: 请求

MySQL 打开binlog日志的方法及注意事项

《MySQL打开binlog日志的方法及注意事项》本文给大家介绍MySQL打开binlog日志的方法及注意事项,本文通过实例代码给大家介绍的非常详细,对大家的学习或工作具有一定的参考借鉴价值,需要... 目录一、默认状态二、如何检查 binlog 状态三、如何开启 binlog3.1 临时开启(重启后失效)

Python中提取文件名扩展名的多种方法实现

《Python中提取文件名扩展名的多种方法实现》在Python编程中,经常会遇到需要从文件名中提取扩展名的场景,Python提供了多种方法来实现这一功能,不同方法适用于不同的场景和需求,包括os.pa... 目录技术背景实现步骤方法一:使用os.path.splitext方法二:使用pathlib模块方法三

Python打印对象所有属性和值的方法小结

《Python打印对象所有属性和值的方法小结》在Python开发过程中,调试代码时经常需要查看对象的当前状态,也就是对象的所有属性和对应的值,然而,Python并没有像PHP的print_r那样直接提... 目录python中打印对象所有属性和值的方法实现步骤1. 使用vars()和pprint()2. 使

CSS实现元素撑满剩余空间的五种方法

《CSS实现元素撑满剩余空间的五种方法》在日常开发中,我们经常需要让某个元素占据容器的剩余空间,本文将介绍5种不同的方法来实现这个需求,并分析各种方法的优缺点,感兴趣的朋友一起看看吧... css实现元素撑满剩余空间的5种方法 在日常开发中,我们经常需要让某个元素占据容器的剩余空间。这是一个常见的布局需求

Python常用命令提示符使用方法详解

《Python常用命令提示符使用方法详解》在学习python的过程中,我们需要用到命令提示符(CMD)进行环境的配置,:本文主要介绍Python常用命令提示符使用方法的相关资料,文中通过代码介绍的... 目录一、python环境基础命令【Windows】1、检查Python是否安装2、 查看Python的安

Maven 配置中的 <mirror>绕过 HTTP 阻断机制的方法

《Maven配置中的<mirror>绕过HTTP阻断机制的方法》:本文主要介绍Maven配置中的<mirror>绕过HTTP阻断机制的方法,本文给大家分享问题原因及解决方案,感兴趣的朋友一... 目录一、问题场景:升级 Maven 后构建失败二、解决方案:通过 <mirror> 配置覆盖默认行为1. 配置示

SpringBoot排查和解决JSON解析错误(400 Bad Request)的方法

《SpringBoot排查和解决JSON解析错误(400BadRequest)的方法》在开发SpringBootRESTfulAPI时,客户端与服务端的数据交互通常使用JSON格式,然而,JSON... 目录问题背景1. 问题描述2. 错误分析解决方案1. 手动重新输入jsON2. 使用工具清理JSON3.