QVTKWidget+itk+vtk 交互问题

时间:2021-08-18 04:30:49

现在做到单张dicom图片读取头信息已经显示的部分,按照东灵工作室写的教程例子先写了一下,他们使用ImageToVTKImageFilter类来对itk和vtk进行数据连接,在用vtk的vtkImageViewer2类连接qt的qvtkwidget进行显示。原理是(来源 http://codeplanet.sinaapp.com/?p=361):

ITK的数据流采用如图1所示流水线(pipeline)作业方式。遵循此方式,所处理的数据来龙去脉清晰明了。在图1中,数据流左边的ImageFile代表待处理的图像文件,而右边ImageFile则代表处理后的图像文件。ImageFileReade:将存贮图像文件读人到内存,形成Image,再调用所希望处理的Fileter算法(比如滤波、分割、配准等具体算法),将对读人的图像数据进行处理。最后,ImageFileWriter将处理的图像数据结果保存,或输人给图形显示系统(如VTK),进行数据可视化。

QVTKWidget+itk+vtk 交互问题

VTK是C++语言开发的、公开源码的、面向对象的数据可视化软件开发包,并不单纯用于医学图像数据的可视化。VTK支持跨平台的编译,可以应用于Windows,Linux等系统。VTK支持包括数量、向量、张量、结构和测定体积等方法一系列可视化算法,VTK还支持包括建模。VTK的数据流同样采用流水线作业方式,如图2所示。

QVTKWidget+itk+vtk 交互问题

vtkSource为整个可视化流的开始,比如读取、生成数据等。vtkFilte对输入源数据进行各种处理,与ITK的Fileter类似。原始数据经过各种Filte:的处理后,再由vtkMappe:转换(或映射)成可供其它算法模块直接使用的数据形式,比如为几何数据。vtkActor类用来描述绘制场景中一个实体,即绘制场景的某个角色(Actor)。通过SetMapper()方法,可以将几何数据的属性传递给Actor,再由vtkRenderer类将处理后的数据结果显示出来。
由于ITK和VTK的数据流都是采用流水线作业形式,比较方便地利用类库本身提供的类型转换类来实现,如图3所示。通过ITKtoVTKImageFilter类操作将ITK和VTK的图像数据连接了起来,VTK负责图像数据的显示,而ITK负责图像的处理,其中VTK的Image Viewe:是一个很方便实现的图像显示类,它在内部管理着vtkImageWindow, vtkRenderer, vtkActor2D和vtkImageMapper等图象类,实现图像查看、显示和可视化功能。
QVTKWidget+itk+vtk 交互问题


因此,东灵他们例子的代码是这么写的:

 // 将reader的输出作为m_pImageViewer的输入,并设置m_pImageViewer与渲染器m_pRenderer的关联  
       m_pImageViewer->SetInput(reader->GetOutput());  
       m_pImageViewer->UpdateDisplayExtent();  
       m_pImageViewer->SetRenderWindow(m_QVTKWidget->GetRenderWindow());  
       m_pImageViewer->SetRenderer(m_pRenderder);  
       m_pImageViewer->SetupInteractor(m_QVTKWidget->GetRenderWindow()->GetInteractor());  
       m_pImageViewer->SetSliceOrientationToXY(); //默认就是这个方向的  
       m_pImageViewer->GetImageActor()->InterpolateOff();  
       m_pRenderder->ResetCamera();  
       m_pRenderder->DrawOn();  
       m_QVTKWidget->GetRenderWindow()->Render();  

我用来实现dicom图像显示及交互,然后发现了只要一点击图片,程序就发生崩溃的错误。网上找了好多原因也没检查出来,后来选择不用vtkimageview2这个类而改用传统的写法,实例化vtk的render、actor、renderwindow,renderwindowinteractor后进行,成功实现。

之后再传到qvtkwidget,也成功了。

后来才发现东灵他们的代码没有new一个vtkRenderWindowInteractor,我在程序初始化的时候增加一个

vtkSmartPointer< vtkRenderWindowInteractor >  m_pInteractor = vtkSmartPointer< vtkRenderWindowInteractor >::New();
然后设置

qvtkWidget->GetRenderWindow()->SetInteractor(m_pInteractor);
测试发现可以成功实现。


ps:在网上搜索解决问题的时候也看到了一个重点就是

重载的vtkInteractorStyleTrackballCamera 也可以是子类 vtkInteractorImage. 同时要注意qvtk 与vtkImageViewer的 RenderWindow与Interactor的共享, 多以qvtk为准. SetInteractor 一定要在 SetRenderWindow之前, 以避免error出现.

-----------------------------------------------------------------------------------------------------

参考博客:

http://blog.csdn.net/www_doling_net/article/details/8668870

http://codeplanet.sinaapp.com/?p=361

http://blog.csdn.net/colddie/article/details/8966521

感谢作者提供的经验。

————————————————