学习三维图像切面的提取后,我们可以实现一个稍微复杂的程序——通过滑动鼠标来切换三维图像切片,这也是医学图像处理软件中一个很基本的功能。实现该功能难点是怎样在VTK中控制鼠标来实时提取图像切片。我们采用观察者/命令(Observer/Command)模式机制来实现。
VTK中鼠标消息是在交互类型对象(interactorstyle)中响应,因此通过为交互类型对象(interactorstyle)添加观察者(observer)来监听相应的消息,当消息触发时,由命令模式执行相应的回调函数。
1.鼠标滑动提取三维图像切片--重新定义切片中心
//记下鼠标Y向变化的幅值大小 int deltaY = lastPos[1] - currPos[1]; //记下鼠标X向变化的幅值大小 int deltaX = lastPos[0] - currPos[0]; //reslice->Update(); double sliceSpacing = reslice->GetOutput()->GetSpacing()[2]; vtkMatrix4x4 *matrix = reslice->GetResliceAxes(); //重新定位切片需要经过的中心点 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]);鼠标竖直滑动重新定义切片的中心位置。注意:理解MultiplyPoint 是两个矩阵相乘,表示 center4*1 = matrix4*4 乘以 point4*1。
2.鼠标滑动设置窗宽、窗位--重新定义窗位的位置
static double colorWindowLevel = 500; static double colorWindowWith = 2000; colorWindowLevel += deltaX; colorMap->GetLookupTable()->SetRange(colorWindowLevel - colorWindowWith / 2, colorWindowLevel + colorWindowWith / 2);鼠标横着滑动,会改变窗位的位置,而窗宽是固定不变的。
关于窗宽、窗位的理解:
窗宽、窗位不是指可视化的窗口的宽度,而是指设置CT图像的灰度值的范围,用于可视化的灰度值范围,CT灰度值范围是0到65536。窗位定以后,根据窗宽值,则可视化的灰度值范围:[窗位-窗宽/2,窗位+窗宽/2]
3.代码
#include <vtkAutoInit.h> VTK_MODULE_INIT(vtkRenderingOpenGL2); #include <vtkSmartPointer.h> #include <vtkMetaImageReader.h> #include <vtkDicomImageReader.h> #include <vtkMatrix4x4.h> #include <vtkLookupTable.h> #include <vtkImageMapToColors.h> #include <vtkImageActor.h> #include <vtkRenderer.h> #include <vtkRenderWindow.h> #include <vtkRenderWindowInteractor.h> #include <vtkInte

本文介绍如何使用VTK库实现实时三维图像切片提取和鼠标交互控制,包括通过鼠标滑动调整切片位置和窗宽窗位,以及代码实现细节。
最低0.47元/天 解锁文章
527

被折叠的 条评论
为什么被折叠?



