VTK 的可视化方法：Cutting Probing

2024-05-03 09:44
文章标签 方法 可视化 vtk cutting probing

本文主要是介绍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。

相关文章

PHP轻松处理千万行数据的方法详解

PHP轻松处理千万行数据的方法详解

《PHP轻松处理千万行数据的方法详解》说到处理大数据集,PHP通常不是第一个想到的语言,但如果你曾经需要处理数百万行数据而不让服务器崩溃或内存耗尽,你就会知道PHP用对了工具有多强大,下面小编就... 目录问题的本质php 中的数据流处理：为什么必不可少生成器：内存高效的迭代方式流量控制：避免系统过载一次性
阅读更多...
python获取指定名字的程序的文件路径的两种方法

python获取指定名字的程序的文件路径的两种方法

《python获取指定名字的程序的文件路径的两种方法》本文主要介绍了python获取指定名字的程序的文件路径的两种方法,文中通过示例代码介绍的非常详细,对大家的学习或者工作具有一定的参考学习价值,需要... 最近在做项目，需要用到给定一个程序名字就可以自动获取到这个程序在Windows系统下的绝对路径，以下
阅读更多...
JavaScript中的高级调试方法全攻略指南

JavaScript中的高级调试方法全攻略指南

《JavaScript中的高级调试方法全攻略指南》什么是高级JavaScript调试技巧,它比console.log有何优势,如何使用断点调试定位问题,通过本文,我们将深入解答这些问题,带您从理论到实... 目录观点与案例结合观点1观点2观点3观点4观点5高级调试技巧详解实战案例断点调试：定位变量错误性能分
阅读更多...
Python中 try / except / else / finally 异常处理方法详解

Python中 try / except / else / finally 异常处理方法详解

《Python中try/except/else/finally异常处理方法详解》：本文主要介绍Python中try/except/else/finally异常处理方法的相关资料,涵... 目录1. 基本结构2. 各部分的作用tryexceptelsefinally3. 执行流程总结4. 常见用法（1）多个e
阅读更多...
JavaScript中比较两个数组是否有相同元素(交集)的三种常用方法

JavaScript中比较两个数组是否有相同元素(交集)的三种常用方法

《JavaScript中比较两个数组是否有相同元素(交集)的三种常用方法》：本文主要介绍JavaScript中比较两个数组是否有相同元素(交集)的三种常用方法,每种方法结合实例代码给大家介绍的非常... 目录引言：为什么"相等"判断如此重要？方法1：使用some()+includes()（适合小数组）方法2
阅读更多...
504 Gateway Timeout网关超时的根源及完美解决方法

504 Gateway Timeout网关超时的根源及完美解决方法

《504GatewayTimeout网关超时的根源及完美解决方法》在日常开发和运维过程中,504GatewayTimeout错误是常见的网络问题之一,尤其是在使用反向代理（如Nginx）或... 目录引言为什么会出现 504 错误？1. 探索 504 Gateway Timeout 错误的根源 1.1 后端
阅读更多...
MySQL 表空却 ibd 文件过大的问题及解决方法

MySQL 表空却 ibd 文件过大的问题及解决方法

《MySQL表空却ibd文件过大的问题及解决方法》本文给大家介绍MySQL表空却ibd文件过大的问题及解决方法,本文给大家介绍的非常详细,对大家的学习或工作具有一定的参考借鉴价值,需要的朋友参考... 目录一、问题背景：表空却 “吃满” 磁盘的怪事二、问题复现：一步步编程还原异常场景1. 准备测试源表与数据
阅读更多...
python 线程池顺序执行的方法实现

python 线程池顺序执行的方法实现

《python线程池顺序执行的方法实现》在Python中,线程池默认是并发执行任务的,但若需要实现任务的顺序执行,文中通过示例代码介绍的非常详细,对大家的学习或者工作具有一定的参考学习价值,需要的朋... 目录方案一：强制单线程（伪顺序执行）方案二：按提交顺序获取结果方案三：任务间依赖控制方案四：队列顺序消
阅读更多...
SpringBoot通过main方法启动web项目实践

SpringBoot通过main方法启动web项目实践

《SpringBoot通过main方法启动web项目实践》SpringBoot通过SpringApplication.run()启动Web项目,自动推断应用类型,加载初始化器与监听器,配置Spring... 目录1. 启动入口：SpringApplication.run()2. SpringApplicat
阅读更多...
使用Java读取本地文件并转换为MultipartFile对象的方法

使用Java读取本地文件并转换为MultipartFile对象的方法

《使用Java读取本地文件并转换为MultipartFile对象的方法》在许多JavaWeb应用中,我们经常会遇到将本地文件上传至服务器或其他系统的需求,在这种场景下,MultipartFile对象非... 目录1. 基本需求2. 自定义 MultipartFile 类3. 实现代码4. 代码解析5. 自定
阅读更多...

Copyright China编程 23002807@qq.com

沪ICP备2023033072号-2