VTK 的可视化方法:Cutting Probing

2024-05-03 09:44

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

VTK 的可视化方法:Cutting & Probing

  • VTK 的可视化方法:Cutting & Probing
    • Cutting
    • Probing

VTK 的可视化方法:Cutting & Probing

Cutting 和 Probing 都是对三维数据截取二维平面的方法。它们的区别是:

  1. Cutting 使用 vtkCutter 类,采用隐式函数截取数据; Probing使用 vtkProbeFilter 类,利用另一个数据来截取数据,或者说采样。
  2. Cutting 截面的大小依据原数据集;Probing 截取数据的大小等于数据截面的大小。
  3. Cutting 截面的分辨率等于原数据集的分辨率,网格也是原数据的网格;Probing 截面的分辨率和网格与数据截面相同,可以自定义。

Cutting

采用vtkCutter类做截取,设置截取平面隐函数cuttingPlane和范围range,即可等间距截取二维平面。

#include "VTKCuttingAndProbing.h"#include <vtkConeSource.h>
#include <vtkSTLReader.h>
#include <vtkMultiBlockPLOT3DReader.h>
#include <vtkPlane.h>
#include <vtkCutter.h>
#include <vtkDataSet.h>
#include <vtkMultiBlockDataSet.h>
#include <vtkShrinkPolyData.h>
#include <vtkStructuredGridGeometryFilter.h>
#include <vtkStructuredGridOutlineFilter.h>
#include <vtkLookupTable.h>
#include <vtkPolyDataMapper.h>
#include <vtkActor.h>
#include <vtkRenderer.h>
#include <vtkRenderWindow.h>VTKCuttingAndProbing::VTKCuttingAndProbing(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));vtkSmartPointer<vtkPlane> cuttingPlane = vtkSmartPointer<vtkPlane>::New();// 设置中心cuttingPlane->SetOrigin(plot3dOutput->GetCenter());// 设置法向量cuttingPlane->SetNormal(1, 0, 0.3);vtkSmartPointer<vtkCutter> planeCut = vtkSmartPointer<vtkCutter>::New();planeCut->SetInputData(plot3dOutput);// 设置截取平面隐函数planeCut->SetCutFunction(cuttingPlane);double range[2] = { -5, 5 };planeCut->GenerateValues(3, range);// 2. filter// 提取作为多边形几何(点,线,表面)的栅格的一部分vtkSmartPointer<vtkStructuredGridGeometryFilter> plane = vtkSmartPointer<vtkStructuredGridGeometryFilter>::New();plane->SetInputData(plot3dOutput);plane->SetExtent(1, 100, 1, 100, 7, 7);// 产生结构化栅格边界的一个线轮廓vtkSmartPointer<vtkStructuredGridOutlineFilter> outline = vtkSmartPointer<vtkStructuredGridOutlineFilter>::New();outline->SetInputData(plot3dOutput);// 颜色映射表vtkSmartPointer<vtkLookupTable> lut = vtkSmartPointer<vtkLookupTable>::New();lut->SetNumberOfColors(256); // 指定颜色映射表中有多少种颜色lut->Build();// 3. mappervtkSmartPointer<vtkPolyDataMapper> planeMapper = vtkSmartPointer<vtkPolyDataMapper>::New();vtkSmartPointer<vtkPolyDataMapper> outlineMapper = vtkSmartPointer<vtkPolyDataMapper>::New();vtkSmartPointer<vtkPolyDataMapper> cuttingMapper = vtkSmartPointer<vtkPolyDataMapper>::New();// 4. actorvtkSmartPointer<vtkActor> planeActor = vtkSmartPointer<vtkActor>::New();vtkSmartPointer<vtkActor> outlineActor = vtkSmartPointer<vtkActor>::New();vtkSmartPointer<vtkActor> cuttingActor = vtkSmartPointer<vtkActor>::New();// 5. renderervtkSmartPointer<vtkRenderer> renderer = vtkSmartPointer<vtkRenderer>::New();renderer->SetBackground(0.3, 0.6, 0.3); // Background Color: Green// 6. connectplaneMapper->SetLookupTable(lut);planeMapper->SetInputConnection(plane->GetOutputPort());planeMapper->SetScalarRange(plot3dOutput->GetScalarRange()); // 设置标量值的范围outlineMapper->SetInputConnection(outline->GetOutputPort());cuttingMapper->SetInputConnection(planeCut->GetOutputPort());cuttingMapper->SetScalarRange(plot3dOutput->GetScalarRange()); // 设置标量值的范围planeActor->SetMapper(planeMapper);outlineActor->SetMapper(outlineMapper);cuttingActor->SetMapper(cuttingMapper);// renderer->AddActor(planeActor);renderer->AddActor(outlineActor);renderer->AddActor(cuttingActor);this->_pVTKWidget->renderWindow()->AddRenderer(renderer);this->_pVTKWidget->renderWindow()->Render();
}VTKCuttingAndProbing::~VTKCuttingAndProbing()
{}

运行结果:

在这里插入图片描述

按W键切换成网格模式:

在这里插入图片描述

Probing

#include "VTKCuttingAndProbing.h"#include <vtkConeSource.h>
#include <vtkPlaneSource.h>
#include <vtkSTLReader.h>
#include <vtkMultiBlockPLOT3DReader.h>
#include <vtkPlane.h>
#include <vtkCutter.h>
#include <vtkTransform.h>
#include <vtkDataSet.h>
#include <vtkMultiBlockDataSet.h>
#include <vtkTransformPolyDataFilter.h>
#include <vtkStructuredGridGeometryFilter.h>
#include <vtkOutlineFilter.h>
#include <vtkStructuredGridOutlineFilter.h>
#include <vtkProbeFilter.h>
#include <vtkLookupTable.h>
#include <vtkPolyDataMapper.h>
#include <vtkAppendPolyData.h>
#include <vtkActor.h>
#include <vtkRenderer.h>
#include <vtkRenderWindow.h>VTKCuttingAndProbing::VTKCuttingAndProbing(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));// CuttingvtkSmartPointer<vtkPlane> cuttingPlane = vtkSmartPointer<vtkPlane>::New();// 设置中心cuttingPlane->SetOrigin(plot3dOutput->GetCenter());// 设置法向量cuttingPlane->SetNormal(1, 0, 0.3);vtkSmartPointer<vtkCutter> planeCut = vtkSmartPointer<vtkCutter>::New();planeCut->SetInputData(plot3dOutput);// 设置截取平面隐函数planeCut->SetCutFunction(cuttingPlane);double range[2] = { -5, 5 };planeCut->GenerateValues(3, range);// 截面数据vtkSmartPointer<vtkPlaneSource> probePlane = vtkSmartPointer<vtkPlaneSource>::New();probePlane->SetResolution(50, 50);// 对于一个数据,我们需要坐标变换才能移动它的位置vtkSmartPointer<vtkTransform> transform1 = vtkSmartPointer<vtkTransform>::New();transform1->Translate(3.7, 0, 28.37);transform1->Scale(5, 5, 5);transform1->RotateY(90);vtkSmartPointer<vtkTransform> transform2 = vtkSmartPointer<vtkTransform>::New();transform2->Translate(9.2, 0, 31.2);transform2->Scale(5, 5, 5);transform2->RotateY(90);vtkSmartPointer<vtkTransform> transform3 = vtkSmartPointer<vtkTransform>::New();transform3->Translate(13.27, 0, 33.3);transform3->Scale(5, 5, 5);transform3->RotateY(90);// 2. filter// 提取作为多边形几何(点,线,表面)的栅格的一部分vtkSmartPointer<vtkStructuredGridGeometryFilter> plane = vtkSmartPointer<vtkStructuredGridGeometryFilter>::New();plane->SetInputData(plot3dOutput);plane->SetExtent(1, 100, 1, 100, 7, 7);// 产生结构化栅格边界的一个线轮廓vtkSmartPointer<vtkStructuredGridOutlineFilter> outline = vtkSmartPointer<vtkStructuredGridOutlineFilter>::New();outline->SetInputData(plot3dOutput);// 颜色映射表vtkSmartPointer<vtkLookupTable> lut = vtkSmartPointer<vtkLookupTable>::New();lut->SetNumberOfColors(256); // 指定颜色映射表中有多少种颜色lut->Build();// 截面 1vtkSmartPointer<vtkTransformPolyDataFilter> tpd1Filter = vtkSmartPointer<vtkTransformPolyDataFilter>::New();tpd1Filter->SetInputConnection(probePlane->GetOutputPort());tpd1Filter->SetTransform(transform1); // 坐标变换// 截面 1 的 outlinevtkSmartPointer<vtkOutlineFilter> outlineTpd1 = vtkSmartPointer<vtkOutlineFilter>::New();outlineTpd1->SetInputConnection(tpd1Filter->GetOutputPort());// 截面 2vtkSmartPointer<vtkTransformPolyDataFilter> tpd2Filter = vtkSmartPointer<vtkTransformPolyDataFilter>::New();tpd2Filter->SetInputConnection(probePlane->GetOutputPort());tpd2Filter->SetTransform(transform2); // 坐标变换// 截面 2 的 outlinevtkSmartPointer<vtkOutlineFilter> outlineTpd2 = vtkSmartPointer<vtkOutlineFilter>::New();outlineTpd2->SetInputConnection(tpd2Filter->GetOutputPort());// 截面 3vtkSmartPointer<vtkTransformPolyDataFilter> tpd3Filter = vtkSmartPointer<vtkTransformPolyDataFilter>::New();tpd3Filter->SetInputConnection(probePlane->GetOutputPort());tpd3Filter->SetTransform(transform3); // 坐标变换// 截面 3 的 outlinevtkSmartPointer<vtkOutlineFilter> outlineTpd3 = vtkSmartPointer<vtkOutlineFilter>::New();outlineTpd3->SetInputConnection(tpd3Filter->GetOutputPort());// 截面集合vtkSmartPointer<vtkAppendPolyData> append = vtkSmartPointer<vtkAppendPolyData>::New();append->AddInputConnection(tpd1Filter->GetOutputPort());append->AddInputConnection(tpd2Filter->GetOutputPort());append->AddInputConnection(tpd3Filter->GetOutputPort());// ProbevtkSmartPointer<vtkProbeFilter> probeFilter = vtkSmartPointer<vtkProbeFilter>::New();probeFilter->SetInputConnection(append->GetOutputPort());probeFilter->SetSourceData(plot3dOutput);// 3. mappervtkSmartPointer<vtkPolyDataMapper> planeMapper = vtkSmartPointer<vtkPolyDataMapper>::New();vtkSmartPointer<vtkPolyDataMapper> outlineMapper = vtkSmartPointer<vtkPolyDataMapper>::New();vtkSmartPointer<vtkPolyDataMapper> cuttingMapper = vtkSmartPointer<vtkPolyDataMapper>::New();vtkSmartPointer<vtkPolyDataMapper> outlineTpd1Mapper = vtkSmartPointer<vtkPolyDataMapper>::New();vtkSmartPointer<vtkPolyDataMapper> outlineTpd2Mapper = vtkSmartPointer<vtkPolyDataMapper>::New();vtkSmartPointer<vtkPolyDataMapper> outlineTpd3Mapper = vtkSmartPointer<vtkPolyDataMapper>::New();vtkSmartPointer<vtkPolyDataMapper> probeMapper = vtkSmartPointer<vtkPolyDataMapper>::New();// 4. actorvtkSmartPointer<vtkActor> planeActor = vtkSmartPointer<vtkActor>::New();vtkSmartPointer<vtkActor> outlineActor = vtkSmartPointer<vtkActor>::New();vtkSmartPointer<vtkActor> cuttingActor = vtkSmartPointer<vtkActor>::New();vtkSmartPointer<vtkActor> outlineTpd1Actor = vtkSmartPointer<vtkActor>::New();vtkSmartPointer<vtkActor> outlineTpd2Actor = vtkSmartPointer<vtkActor>::New();vtkSmartPointer<vtkActor> outlineTpd3Actor = vtkSmartPointer<vtkActor>::New();vtkSmartPointer<vtkActor> probeActor = vtkSmartPointer<vtkActor>::New();// 5. renderervtkSmartPointer<vtkRenderer> renderer = vtkSmartPointer<vtkRenderer>::New();renderer->SetBackground(0.3, 0.6, 0.3); // Background Color: Green// 6. connectplaneMapper->SetLookupTable(lut);planeMapper->SetInputConnection(plane->GetOutputPort());planeMapper->SetScalarRange(plot3dOutput->GetScalarRange()); // 设置标量值的范围outlineMapper->SetInputConnection(outline->GetOutputPort());cuttingMapper->SetInputConnection(planeCut->GetOutputPort());cuttingMapper->SetScalarRange(plot3dOutput->GetScalarRange()); // 设置标量值的范围outlineTpd1Mapper->SetInputConnection(outlineTpd1->GetOutputPort());outlineTpd2Mapper->SetInputConnection(outlineTpd2->GetOutputPort());outlineTpd3Mapper->SetInputConnection(outlineTpd3->GetOutputPort());probeMapper->SetInputConnection(probeFilter->GetOutputPort());probeMapper->SetScalarRange(plot3dOutput->GetScalarRange()); // 设置标量值的范围planeActor->SetMapper(planeMapper);outlineActor->SetMapper(outlineMapper);cuttingActor->SetMapper(cuttingMapper);outlineTpd1Actor->SetMapper(outlineTpd1Mapper);outlineTpd2Actor->SetMapper(outlineTpd2Mapper);outlineTpd3Actor->SetMapper(outlineTpd3Mapper);probeActor->SetMapper(probeMapper);// renderer->AddActor(planeActor);renderer->AddActor(outlineActor);// renderer->AddActor(cuttingActor);renderer->AddActor(outlineTpd1Actor);renderer->AddActor(outlineTpd2Actor);renderer->AddActor(outlineTpd3Actor);renderer->AddActor(probeActor);this->_pVTKWidget->renderWindow()->AddRenderer(renderer);this->_pVTKWidget->renderWindow()->Render();
}VTKCuttingAndProbing::~VTKCuttingAndProbing()
{}

运行结果:

在这里插入图片描述

按W键切换成网格模式,可以看出分辨率为50*50:

在这里插入图片描述

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



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

相关文章

Java中流式并行操作parallelStream的原理和使用方法

《Java中流式并行操作parallelStream的原理和使用方法》本文详细介绍了Java中的并行流(parallelStream)的原理、正确使用方法以及在实际业务中的应用案例,并指出在使用并行流... 目录Java中流式并行操作parallelStream0. 问题的产生1. 什么是parallelS

MySQL数据库双机热备的配置方法详解

《MySQL数据库双机热备的配置方法详解》在企业级应用中,数据库的高可用性和数据的安全性是至关重要的,MySQL作为最流行的开源关系型数据库管理系统之一,提供了多种方式来实现高可用性,其中双机热备(M... 目录1. 环境准备1.1 安装mysql1.2 配置MySQL1.2.1 主服务器配置1.2.2 从

Python版本信息获取方法详解与实战

《Python版本信息获取方法详解与实战》在Python开发中,获取Python版本号是调试、兼容性检查和版本控制的重要基础操作,本文详细介绍了如何使用sys和platform模块获取Python的主... 目录1. python版本号获取基础2. 使用sys模块获取版本信息2.1 sys模块概述2.1.1

Python实现字典转字符串的五种方法

《Python实现字典转字符串的五种方法》本文介绍了在Python中如何将字典数据结构转换为字符串格式的多种方法,首先可以通过内置的str()函数进行简单转换;其次利用ison.dumps()函数能够... 目录1、使用json模块的dumps方法:2、使用str方法:3、使用循环和字符串拼接:4、使用字符

Python版本与package版本兼容性检查方法总结

《Python版本与package版本兼容性检查方法总结》:本文主要介绍Python版本与package版本兼容性检查方法的相关资料,文中提供四种检查方法,分别是pip查询、conda管理、PyP... 目录引言为什么会出现兼容性问题方法一:用 pip 官方命令查询可用版本方法二:conda 管理包环境方法

Linux云服务器手动配置DNS的方法步骤

《Linux云服务器手动配置DNS的方法步骤》在Linux云服务器上手动配置DNS(域名系统)是确保服务器能够正常解析域名的重要步骤,以下是详细的配置方法,包括系统文件的修改和常见问题的解决方案,需要... 目录1. 为什么需要手动配置 DNS?2. 手动配置 DNS 的方法方法 1:修改 /etc/res

JavaScript对象转数组的三种方法实现

《JavaScript对象转数组的三种方法实现》本文介绍了在JavaScript中将对象转换为数组的三种实用方法,文中通过示例代码介绍的非常详细,对大家的学习或者工作具有一定的参考学习价值,需要的朋友... 目录方法1:使用Object.keys()和Array.map()方法2:使用Object.entr

SpringBoot中ResponseEntity的使用方法举例详解

《SpringBoot中ResponseEntity的使用方法举例详解》ResponseEntity是Spring的一个用于表示HTTP响应的全功能对象,它可以包含响应的状态码、头信息及响应体内容,下... 目录一、ResponseEntity概述基本特点:二、ResponseEntity的基本用法1. 创

java中判断json key是否存在的几种方法

《java中判断jsonkey是否存在的几种方法》在使用Java处理JSON数据时,如何判断某一个key是否存在?本文就来介绍三种方法,文中通过示例代码介绍的非常详细,对大家的学习或者工作具有一定的... 目http://www.chinasem.cn录第一种方法是使用 jsONObject 的 has 方法

java中ssh2执行多条命令的四种方法

《java中ssh2执行多条命令的四种方法》本文主要介绍了java中ssh2执行多条命令的四种方法,包括分号分隔、管道分隔、EOF块、脚本调用,可确保环境配置生效,提升操作效率,具有一定的参考价值,感... 目录1 使用分号隔开2 使用管道符号隔开3 使用写EOF的方式4 使用脚本的方式大家平时有没有遇到自