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。

相关文章

Python中注释使用方法举例详解

Python中注释使用方法举例详解

《Python中注释使用方法举例详解》在Python编程语言中注释是必不可少的一部分,它有助于提高代码的可读性和维护性,：本文主要介绍Python中注释使用方法的相关资料,需要的朋友可以参考下... 目录一、前言二、什么是注释？示例：三、单行注释语法：以 China编程# 开头，后面的内容为注释内容示例：示例：四
阅读更多...
一文详解Git中分支本地和远程删除的方法

一文详解Git中分支本地和远程删除的方法

《一文详解Git中分支本地和远程删除的方法》在使用Git进行版本控制的过程中,我们会创建多个分支来进行不同功能的开发,这就容易涉及到如何正确地删除本地分支和远程分支,下面我们就来看看相关的实现方法吧... 目录技术背景实现步骤删除本地分支删除远程www.chinasem.cn分支同步删除信息到其他机器示例步骤
阅读更多...
在Golang中实现定时任务的几种高效方法

在Golang中实现定时任务的几种高效方法

《在Golang中实现定时任务的几种高效方法》本文将详细介绍在Golang中实现定时任务的几种高效方法,包括time包中的Ticker和Timer、第三方库cron的使用,以及基于channel和go... 目录背景介绍目的和范围预期读者文档结构概述术语表核心概念与联系故事引入核心概念解释核心概念之间的关系
阅读更多...
在Linux终端中统计非二进制文件行数的实现方法

在Linux终端中统计非二进制文件行数的实现方法

《在Linux终端中统计非二进制文件行数的实现方法》在Linux系统中,有时需要统计非二进制文件（如CSV、TXT文件）的行数,而不希望手动打开文件进行查看,例如,在处理大型日志文件、数据文件时,了解... 目录在linux终端中统计非二进制文件的行数技术背景实现步骤1. 使用wc命令2. 使用grep命令
阅读更多...
Python中Tensorflow无法调用GPU问题的解决方法

Python中Tensorflow无法调用GPU问题的解决方法

《Python中Tensorflow无法调用GPU问题的解决方法》文章详解如何解决TensorFlow在Windows无法识别GPU的问题,需降级至2.10版本,安装匹配CUDA11.2和cuDNN... 当用以下代码查看GPU数量时，gpuspython返回的是一个空列表，说明tensorflow没有找到
阅读更多...
XML重复查询一条Sql语句的解决方法

XML重复查询一条Sql语句的解决方法

《XML重复查询一条Sql语句的解决方法》文章分析了XML重复查询与日志失效问题,指出因DTO缺少@Data注解导致日志无法格式化、空指针风险及参数穿透,进而引发性能灾难,解决方案为在Controll... 目录一、核心问题：从SQL重复执行到日志失效二、根因剖析：DTO断裂引发的级联故障三、解决方案：修复
阅读更多...
C++ 检测文件大小和文件传输的方法示例详解

C++ 检测文件大小和文件传输的方法示例详解

《C++检测文件大小和文件传输的方法示例详解》文章介绍了在C/C++中获取文件大小的三种方法,推荐使用stat()函数,并详细说明了如何设计一次性发送压缩包的结构体及传输流程,包含CRC校验和自动解... 目录检测文件的大小✅ 方法一：使用 stat() 函数（推荐）✅ 用法示例：✅ 方法二：使用 fsee
阅读更多...
Java继承映射的三种使用方法示例

Java继承映射的三种使用方法示例

《Java继承映射的三种使用方法示例》继承在Java中扮演着重要的角色,它允许我们创建一个类(子类),该类继承另一个类(父类)的所有属性和方法,：本文主要介绍Java继承映射的三种使用方法示例,需... 目录前言一、单表继承（Single Table Inheritance）1-1、原理1-2、使用方法1-
阅读更多...
SpringBoot中使用Flux实现流式返回的方法小结

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

《SpringBoot中使用Flux实现流式返回的方法小结》文章介绍流式返回（StreamingResponse）在SpringBoot中通过Flux实现,优势包括提升用户体验、降低内存消耗、支持长连... 目录背景流式返回的核心概念与优势1. 提升用户体验2. 降低内存消耗3. 支持长连接与实时通信在Sp
阅读更多...
Conda虚拟环境的复制和迁移的四种方法实现

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

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

Copyright China编程 23002807@qq.com

沪ICP备2023033072号-2