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。

相关文章

Nginx设置连接超时并进行测试的方法步骤

Nginx设置连接超时并进行测试的方法步骤

《Nginx设置连接超时并进行测试的方法步骤》在高并发场景下,如果客户端与服务器的连接长时间未响应,会占用大量的系统资源,影响其他正常请求的处理效率,为了解决这个问题,可以通过设置Nginx的连接... 目录设置连接超时目的操作步骤测试连接超时测试方法：总结：设置连接超时目的设置客户端与服务器之间的连接
阅读更多...
Java判断多个时间段是否重合的方法小结

Java判断多个时间段是否重合的方法小结

《Java判断多个时间段是否重合的方法小结》这篇文章主要为大家详细介绍了Java中判断多个时间段是否重合的方法,文中的示例代码讲解详细,感兴趣的小伙伴可以跟随小编一起学习一下... 目录判断多个时间段是否有间隔判断时间段集合是否与某时间段重合判断多个时间段是否有间隔实体类内容public class D
阅读更多...
Python使用国内镜像加速pip安装的方法讲解

Python使用国内镜像加速pip安装的方法讲解

《Python使用国内镜像加速pip安装的方法讲解》在Python开发中,pip是一个非常重要的工具,用于安装和管理Python的第三方库,然而,在国内使用pip安装依赖时,往往会因为网络问题而导致速... 目录一、pip 工具简介1. 什么是 pip？2. 什么是 -i 参数？二、国内镜像源的选择三、如何
阅读更多...
IDEA编译报错“java: 常量字符串过长”的原因及解决方法

IDEA编译报错“java: 常量字符串过长”的原因及解决方法

《IDEA编译报错“java:常量字符串过长”的原因及解决方法》今天在开发过程中,由于尝试将一个文件的Base64字符串设置为常量,结果导致IDEA编译的时候出现了如下报错java:常量字符串过长,... 目录一、问题描述二、问题原因2.1 理论角度2.2 源码角度三、解决方案解决方案①：StringBui
阅读更多...
Linux使用nload监控网络流量的方法

Linux使用nload监控网络流量的方法

《Linux使用nload监控网络流量的方法》Linux中的nload命令是一个用于实时监控网络流量的工具,它提供了传入和传出流量的可视化表示,帮助用户一目了然地了解网络活动,本文给大家介绍了Linu... 目录简介安装示例用法基础用法指定网络接口限制显示特定流量类型指定刷新率设置流量速率的显示单位监控多个
阅读更多...
Java覆盖第三方jar包中的某一个类的实现方法

Java覆盖第三方jar包中的某一个类的实现方法

《Java覆盖第三方jar包中的某一个类的实现方法》在我们日常的开发中,经常需要使用第三方的jar包,有时候我们会发现第三方的jar包中的某一个类有问题,或者我们需要定制化修改其中的逻辑,那么应该如何... 目录一、需求描述二、示例描述三、操作步骤四、验证结果五、实现原理一、需求描述需求描述如下：需要在
阅读更多...
JavaScript中的reduce方法执行过程、使用场景及进阶用法

JavaScript中的reduce方法执行过程、使用场景及进阶用法

《JavaScript中的reduce方法执行过程、使用场景及进阶用法》：本文主要介绍JavaScript中的reduce方法执行过程、使用场景及进阶用法的相关资料,reduce是JavaScri... 目录1. 什么是reduce2. reduce语法2.1 语法2.2 参数说明3. reduce执行过程
阅读更多...
C#中读取XML文件的四种常用方法

C#中读取XML文件的四种常用方法

《C#中读取XML文件的四种常用方法》Xml是Internet环境中跨平台的,依赖于内容的技术,是当前处理结构化文档信息的有力工具,下面我们就来看看C#中读取XML文件的方法都有哪些吧... 目录XML简介格式C#读取XML文件方法使用XmlDocument使用XmlTextReader/XmlTextWr
阅读更多...
C++初始化数组的几种常见方法(简单易懂)

C++初始化数组的几种常见方法(简单易懂)

《C++初始化数组的几种常见方法(简单易懂)》本文介绍了C++中数组的初始化方法,包括一维数组和二维数组的初始化,以及用new动态初始化数组,在C++11及以上版本中,还提供了使用std::array... 目录1、初始化一维数组1.1、使用列表初始化（推荐方式）1.2、初始化部分列表1.3、使用std::
阅读更多...
oracle DBMS_SQL.PARSE的使用方法和示例

oracle DBMS_SQL.PARSE的使用方法和示例

《oracleDBMS_SQL.PARSE的使用方法和示例》DBMS_SQL是Oracle数据库中的一个强大包,用于动态构建和执行SQL语句,DBMS_SQL.PARSE过程解析SQL语句或PL/S... 目录语法示例注意事项DBMS_SQL 是 oracle 数据库中的一个强大包，它允许动态地构建和执行
阅读更多...

Copyright China编程 23002807@qq.com

沪ICP备2023033072号-2