VTK使用交互器来从三维体数据中提取二维切片

VTK中鼠标消息是在交互类型对象(interactorstyle)中响应,因此通过为交互类型对象(interactorstyle)添加观察者(observer)来监听相应的消息,当消息触发时,由命令模式执行相应的回调函数

vtkImageInteractionCallback继承自vtkCommand类,并覆盖父类函数Execute()。
该类提供了两个接口:SetImageReslice和SetInteractor。
SetImageReslice用以设置vtkImageSlice对象,vtkImageSlice根据设置的变换矩阵提取三维图像切片。SetInteractor用以设置vtkRenderWindowInteractor,vtkRenderWindowInteractor类对象负责每次提取切片后刷新视图

class vtkImageInteractionCallback : public vtkCommand
{
public:static vtkImageInteractionCallback *New() //回调函数初始化函数{return new vtkImageInteractionCallback;}vtkImageInteractionCallback(){this->Slicing = 0;//切片提取标志位,为1时提取切片。用来检验鼠标左键是否已经按下this->ImageReslice = 0;//vtkImageSlice对象的变换矩阵,用来确定切面在三维图像中的位置this->Interactor = 0;}void SetImageReslice(vtkImageReslice *reslice){this->ImageReslice = reslice;}vtkImageReslice *GetImageReslice(){return this->ImageReslice;}void SetInteractor(vtkRenderWindowInteractor *interactor){this->Interactor = interactor;}vtkRenderWindowInteractor *GetInteractor(){return  this->Interactor;}virtual void Execute(vtkObject * ,unsigned long event,void *){//具体内容放下面}
private:int Slicing;vtkImageReslice *ImageReslice;vtkRenderWindowInteractor *Interactor;
};

   需要设置四个参数,窗位(ColorLevel)、窗宽(ColorWindow)、切片(Slice)和切片方向(Orientation)。按下鼠标左键拖动鼠标,可以调节图像的窗宽窗位,从而显示不同灰度范围内容;按下鼠标右键拖动鼠标可以放缩图像。当然这些交互操作可以由用户根据需要自己定义vtkInteractorStyle子类,并响应相应的操作。

  而显示三维图像时,需要确定当前显示切片和方向。vtkImageViewer2提供了SetSlice()函数设置切片号,SetSliceOrientationToXY()则将切片的方向设置为垂直XY平面方向。此外还可以设置为垂直YZ或者XZ平面方向,其对应函数分别为SetSliceOrientationToYZ()和SetSliceOrientationToXZ()。默认情况下切片方向为垂直于XY平面即沿着Z轴方向,根据设置的切片号获取Z轴方向的具体切片进行显示。

下面重点看Execute函数,该函数提供了具体的切片提取功能。在该函数里面,主要监听了三个消息:
vtkCommand::LeftButtonPressEvent,
vtkCommand::LeftButtonReleaseEvent,
vtkCommand::MouseMoveEvent,
前两个消息分别是鼠标左键的按下和弹起消息。当鼠标左键按下时,就设置切片提取标志为1,而当弹起时,将标志置为0。这样在鼠标移动时,只有在确定切片提取标志为1时,执行切片提取功能。
vtkCommand::MouseMoveEvent即为鼠标移动消息。当检测到该消息时,首先检查切片提取标志,当为1时提取切片。提取切片时,需要为vtkImageSlice对象设置变换矩阵。这里在函数开始时,首先获取了鼠标滑动的前后两次点的位置lastPos和currPos。然后根据两点的Y坐标差deltaY,计算新的中心点center并变换至vtkImageSlice当前变换矩阵中,得到变换中心点,将其设置到原来的变换矩阵matrix中,并设置到vtkImageSlice中,最后执行interactor->Render()即可不断的根据鼠标移动刷新图像。

virtual void Execute(vtkObject * ,unsigned long event,void *)
{vtkRenderWindowInteractor *interactor = GetInteractor();int lastPos[2];interactor->GetLastEventPosition(lastPos);int currPos[2];interactor->GetEventPosition(currPos);if (event == vtkCommand::LeftButtonPressEvent)//如果当前事件是鼠标左键按下{this->Slicing = 1; //标志位。当鼠标左键按下时,就设置切片提取标志为1}else if (event == vtkCommand::LeftButtonReleaseEvent)/如果当前事件是鼠标右键按下{this->Slicing = 0; //标志位 }else if (event == vtkCommand::MouseMoveEvent)/如果当前事件是鼠标移动{if (this->Slicing)//检验鼠标左键已经按下{vtkImageReslice *reslice = this->ImageReslice;//设置当前切片图像位置的矩阵为this->ImageReslice;但在更新之前,不会将这个矩阵数据填充到reslice中去//记下鼠标Y向变化的幅值大小int deltaY = lastPos[1] - currPos[1];//计算鼠标滑动的前后两次点的位置lastPos和currPos。然后计算两点的Y坐标的差deltaYreslice->Update();//使用更新操作,为reslice获取当前切片图像位置的矩阵double sliceSpacing = reslice->GetOutput()->GetSpacing()[2];获取reslice矩阵中的参数vtkMatrix4x4 *matrix = reslice->GetResliceAxes();//定义需要使用的新的中心店center转换至vtkImageSlice当前变换矩阵中的转换矩阵,但这里没有更新,就没有填充//重新定位切片需要经过的中心点double point[4];double center[4];point[0] = 0;point[1] = 0;point[2] = sliceSpacing*deltaY;point[3] = 1.0;matrix->MultiplyPoint(point, center);//填充矩阵matrix->SetElement(0, 3, center[0]);//填充矩阵matrix->SetElement(1, 3, center[1]);//填充矩阵matrix->SetElement(2, 3, center[2]);//填充矩阵interactor->Render();}else//没有触发上述三个事件{vtkInteractorStyle *style = vtkInteractorStyle::SafeDownCast(interactor->GetInteractorStyle());if (style){style->OnMouseMove();}}}
}

Command对象定义完毕后,即可为交互对象InteractorStyle添加观察者,响应鼠标消息。
这里主要是定义了vtkImageInteractionCallback对象,并设置vtkImageSlice对象和vtkRenderWindowInteractor对象。然后为交互对象vtkInteractorStyle添加观察者来监控相应的消息。

vtkSmartPointer<vtkRenderWindowInteractor> renderWindowInteractor =vtkSmartPointer<vtkRenderWindowInteractor>::New();//定义交互对象vtkSmartPointer<vtkInteractorStyleImage> imagestyle =vtkSmartPointer<vtkInteractorStyleImage>::New();//定义交互模式renderWindowInteractor->SetInteractorStyle(imagestyle);//设置交互模式renderWindowInteractor->SetRenderWindow(renderWindow);//连接交互器与显示窗口renderWindowInteractor->Initialize();//交互器初始化//****************建立 观察者-命令 模式****************//vtkSmartPointer<vtkImageInteractionCallback> callback =vtkSmartPointer<vtkImageInteractionCallback>::New();//设置回调函数callback->SetImageReslice(reslice);//为vtkImageSlice对象设置变换矩阵,确定矩阵对应的切片位置callback->SetInteractor(renderWindowInteractor);imagestyle->AddObserver(vtkCommand::MouseMoveEvent, callback);//鼠标移动对应的回调函数imagestyle->AddObserver(vtkCommand::LeftButtonPressEvent, callback);//按下鼠标左键对应的回调函数imagestyle->AddObserver(vtkCommand::LeftButtonReleaseEvent, callback);//按下鼠标右键对应的回调函数renderWindowInteractor->Start();

这里主要是三个消息:
vtkCommand::LeftButtonPressEvent,
vtkCommand::LeftButtonReleaseEvent,
vtkCommand::MouseMoveEvent,
当响应到这三个消息时,立即执行vtkImageInteractionCallback的Execute函数,以便实现切片的实时提取和更新。完成以后,运行程序,当鼠标在图像上移动时,会发现图像会跟着鼠标的移动而变化。

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/794449.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

设计模式总结-组合模式

组合设计模式 模式动机模式定义模式结构组合模式实例与解析实例一&#xff1a;水果盘实例二&#xff1a;文件浏览 更复杂的组合总结 模式动机 对于树形结构&#xff0c;当容器对象&#xff08;如文件夹&#xff09;的某一个方法被调用时&#xff0c;将遍历整个树形结构&#x…

市面上的常见传感器价格

自动窗帘传感器的价格因品牌、型号、功能以及市场供需关系等因素而有所不同。一般来说&#xff0c;普通自动窗帘传感器的价格可能在几十到几百元人民币不等。然而&#xff0c;一些高端、智能型或具有特殊功能的传感器可能会更贵&#xff0c;价格可能达到几百元甚至更高。 购买时…

Spark面试整理-解释Spark中的内存管理和持久化机制

在Apache Spark中,内存管理和持久化机制是核心特性,它们对于提高大规模数据处理的效率和性能至关重要。 内存管理 统一的内存管理:Spark使用统一的内存管理模型,将执行内存(用于计算如shuffle、join等)和存储内存(用于缓存数据如RDDs)合并在一起。这种模型提供了更高的…

刷题之Leetcode27题(超级详细)

27. 移除元素 力扣题目链接(opens new window)https://leetcode.cn/problems/remove-element/ 给你一个数组 nums 和一个值 val&#xff0c;你需要 原地 移除所有数值等于 val 的元素&#xff0c;并返回移除后数组的新长度。 不要使用额外的数组空间&#xff0c;你必须仅使用…

面试算法-132-和至少为 K 的最短子数组

题目 给你一个整数数组 nums 和一个整数 k &#xff0c;找出 nums 中和至少为 k 的 最短非空子数组 &#xff0c;并返回该子数组的长度。如果不存在这样的 子数组 &#xff0c;返回 -1 。 子数组 是数组中 连续 的一部分。 示例 1&#xff1a; 输入&#xff1a;nums [1], …

一篇文章带你掌握二叉树(附带二叉树基本操作完整代码演示,和两种思路)

【本长内容】 1. 掌握树的基本概念 2. 掌握二叉树概念及特性 3. 掌握二叉树的基本操作 4. 完成二叉树相关的面试题练习 1. 树形结构 1.1 概念 树是一种非线性的数据结构&#xff0c;它是由n&#xff08;n>0&#xff09;个有限结点组成一个具有层次关系的集合。把它叫做树是…

Vue - 2( 10000 字 Vue 入门级教程)

一&#xff1a;Vue 1.1 绑定样式 1.1.1 绑定 class 样式 <!DOCTYPE html> <html><head><meta charset"UTF-8" /><title>绑定样式</title><style>......</style><script type"text/javascript" src&…

从零开始:成功进入IT行业的方法与技巧

如今&#xff0c;信息技术&#xff08;IT&#xff09;行业成为了就业市场上的热门领域。由于其快速发展和广阔的职业机会&#xff0c;许多人希望能够进入这个行业。然而&#xff0c;对于没有任何相关背景知识的人来说&#xff0c;要成功进入IT行业可能会面临一些挑战。本文将分…

基于java web的超市管理系统

摘要 随着社会经济的不断发展&#xff0c;人们的生活水平不断提高。越来越多的零售行业得到了快速的发展&#xff0c;以最常见的超市最为明显。零售行业繁荣的背后也随之带来了许多行业隐患&#xff0c;越来越激烈的行业竞争不断的要求经营者更加高要求的管理超市内部的整个供…

代码随想录第31天 | 455.分发饼干 、376. 摆动序列、53. 最大子序和

一、前言 参考文献&#xff1a;代码随想录 今天的内容是贪心算法&#xff0c;这个算法分为两个极端&#xff0c;一个极端是很简单&#xff0c;靠常识就可以解出来&#xff0c;另外一个是&#xff0c;你怎么想也想不出来&#xff0c;只能看题解的那种。 and 对第一天和第二天…

微信小程序上传到gitee

共三步 1、新建gitee仓库 点号&#xff0c;新建仓库&#xff0c;填入仓库信息新建即可 2、修改版本管理参数 微信开发者工具中点开版本管理&#xff0c;未初始化&#xff0c;需要先点初始化 接下来将设置中的通用、网络认证、远程3个部分的参数填写好 通用&#xff1a;核对…

DashOJ-8.奇偶统计

题目链接&#xff1a; 题目详情 - 奇偶统计 - DashOJ 思路&#xff1a; &#xff08;while循环加if分支语句&#xff09; 巧用死循环 while(1) 然后在里面第一句就判断输入的数字是否等于0 if(x0) &#xff0c;如果 等于0就直接break跳出循环 或者用 while(cin>>x) 代…

Win11解决Fiddler everyone卡在启动界面

问题描述 fiddler在window11卡在启动界面,重新安装也不行 解决方案 右击以管理员运行

echarts 地图 自己圈地图 乡镇街道

这个是方式是我实在不愿意做的&#xff01; 如果有现成的最好&#xff0c;没有办法的情况下再用这个东西。 今天公司有一个项目&#xff0c;地方划分了一块区域&#xff0c;但是国家没有审核&#xff0c;但是项目里面用到了一个地图展示数据&#xff01;然后就需要我们自己把…

Java学习笔记02:面向对象

3.1.3 面向对象特征 唯一性&#xff08;每个对象都是唯一的&#xff09;、分类性、继承性&#xff08;Java是单继承&#xff0c;C是多继承&#xff09;、多态性&#xff08;子类重写父类方法&#xff09; 3.2.1 类的声明 访问符&#xff1a;public、protected、private。修饰…

非conda环境(powershell与cmd)使用opencv-python库

1.安装python 3.11.7 Python Release Python 3.11.7 | Python.org 安装成功后,查看版本为3.11.7 2.安装numpy pip install numpy 3.安装matplotlib pip install matplotlib matplotlib库安装成功

论文阅读-Directed Greybox Fuzzing 定向灰盒模糊(2017)

摘要&#xff1a; 本文介绍了一种新的Directed Greybox Fuzzing (DGF)方法&#xff0c;可以有效地指导模糊测试&#xff0c;以便更快地到达目标程序位置。作者开发了一种基于模拟退火的能量调度算法&#xff0c;逐渐为更接近目标位置的种子分配更多能量&#xff0c;同时减少远离…

k8s CNI Calico 网络模式总结

目录 calico架构图 IPIP模式下的架构图 calico 核心组件 Overlay 网络模式&#xff1a; Pod IP对外暴露 不对外暴露&#xff1a; 实现对外暴露的方法&#xff1a; overlay模式下的网络MTU Iptables & ipvs overlay的主要缺点&#xff1a; Full-mesh Unoverla…

Python使用executemany实现数据批量插入

Python的pymysql中提供有cursor.execute()实现数据的逐行插入&#xff0c;但是当数据量较大时&#xff0c;这种插入方法十分耗时。这时我们可以使用cursor.executemany()函数来实现数据的批量插入。 executemany()用法 在数据库连接后&#xff0c;使用cursor.excutemany(sql,…

【Linux】基础IO—1

提示&#xff1a;文章写完后&#xff0c;目录可以自动生成&#xff0c;如何生成可参考右边的帮助文档 目录 文章目录 前言 先来段代码回顾C文件接口 "w"写文件 "a"追加文件 "r"读文件 输出信息到显示器&#xff0c;你有哪些方法 stdin & std…