VTK-医学三维图像四视图显示以及鼠标滑轮控制切片交互

一、四视图介绍

vs直接创建Qt Widgets Application项目

四视图界面，采用vtk和QVTKOpenGLNativeWidget实现

示例：pandas 是基于NumPy 的一种工具，该工具是为了解决数据分析任务而创建的。

二、ui界面

1、如何在QDesigner中添加QVTKOpenGLNativeWidget或者其他跟vtk相关的widget

将vtk编译好的QVTKWidgetPlugin.dll（…\bin\Release）

放入QT的designer的路径下（…l\QT\5.15.0\msvc_64\plugins\designer）

上述完成后，Qt deginger中出现QVTK组件

2、如果没有QVTKWidgetPlugin.dll怎么办？

手动配置vtk的cmake中的QT相关编译项

3、在QTdesigner中拖动widget，修改为QVTKOpenGLNativeWidget

三、读入数据

直接在QMainWindow中使用vtk读取三维医学图像数据

代码如下（示例）：

vtkSmartPointer<vtkMetaImageReader> reader = vtkSmartPointer<vtkMetaImageReader>::New();reader->SetFileName("xx.mhd");reader->Update();

该处使用的url网络请求的数据。

四、3D显示

针对原始图像采用体绘制方法（如果是分割mask，可以采用面绘制方法）

代码如下（示例）：

vtkSmartPointer<vtkFixedPointVolumeRayCastMapper> volumeMapper = vtkSmartPointer<vtkFixedPointVolumeRayCastMapper>::New();volumeMapper->SetInputData(image);vtkSmartPointer<vtkVolumeProperty> volumeProperty = vtkSmartPointer<vtkVolumeProperty>::New();volumeProperty->SetInterpolationTypeToLinear();volumeProperty->ShadeOn(); //打开或者关闭阴影测试volumeProperty->SetAmbient(0.4);volumeProperty->SetDiffuse(1.6); //漫反射volumeProperty->SetSpecular(0.2); //镜面反射//设置不透明度vtkSmartPointer<vtkPiecewiseFunction> compositeOpacity = vtkSmartPointer<vtkPiecewiseFunction>::New();compositeOpacity->AddPoint(70, 0.00);compositeOpacity->AddPoint(90, 0.40);compositeOpacity->AddPoint(180, 0.60);volumeProperty->SetScalarOpacity(compositeOpacity); //设置不透明度传输函数vtkSmartPointer<vtkPiecewiseFunction> volumeGradientOpacity = vtkSmartPointer<vtkPiecewiseFunction>::New();volumeGradientOpacity->AddPoint(10, 0.0);volumeGradientOpacity->AddPoint(90, 0.5);volumeGradientOpacity->AddPoint(100, 1.0);vtkSmartPointer<vtkColorTransferFunction> color = vtkSmartPointer<vtkColorTransferFunction>::New();color->AddRGBPoint(0.000, 0.00, 0.00, 0.00);color->AddRGBPoint(64.00, 1.00, 0.52, 0.30);color->AddRGBPoint(190.0, 1.00, 1.00, 1.00);color->AddRGBPoint(220.0, 0.20, 0.20, 0.20);volumeProperty->SetColor(color);vtkSmartPointer<vtkVolume> volume = vtkSmartPointer<vtkVolume>::New();volume->SetMapper(volumeMapper);volume->SetProperty(volumeProperty);this->m_3DRenderder = vtkSmartPointer<vtkRenderer >::New();this->m_3DRenderder->SetBackground(0.3, 0.4, 0.3);this->m_3DRenderder->AddVolume(volume);this->m_3DRenderder->ResetCamera();ui.m_3DWidget->GetRenderWindow()->AddRenderer(this->m_3DRenderder);ui.m_3DWidget->GetRenderWindow()->SetSize(640, 480);ui.m_3DWidget->GetRenderWindow()->Render();ui.m_3DWidget->GetRenderWindow()->SetWindowName("VolumeRendering");vtkRenderWindowInteractor *iren = ui.m_3DWidget->GetInteractor();iren->Render();ui.m_3DWidget->show();

备注：ui.m_3DWidget为界面中的QVTKOpenGLNativeWidget组件。

this->m_3DRenderder 为头文件中定义的vtkRenderer。

五、2D切片显示

1、交互函数

继承自vtkInteractorStyleImage，重写滑轮滚动事件

代码如下（示例）：

class myVtkInteractorStyleImage : public vtkInteractorStyleImage{public:static myVtkInteractorStyleImage* New();vtkTypeMacro(myVtkInteractorStyleImage, vtkInteractorStyleImage);protected:vtkImageViewer2* ImageViewer;int Slice;int MinSlice;int MaxSlice;public:void SetImageViewer(vtkImageViewer2* imageViewer){this->ImageViewer = imageViewer;this->MinSlice = imageViewer->GetSliceMin();this->MaxSlice = imageViewer->GetSliceMax();this->Slice = (this->MinSlice + this->MaxSlice) / 2;this->ImageViewer->SetSlice(this->Slice);this->ImageViewer->Render();}protected:virtual void OnMouseWheelForward(){if (this->Slice < this->MaxSlice){this->Slice += 1;this->ImageViewer->SetSlice(this->Slice);this->ImageViewer->Render();}}virtual void OnMouseWheelBackward(){if (this->Slice > this->MinSlice){this->Slice -= 1;this->ImageViewer->SetSlice(this->Slice);this->ImageViewer->Render();}}};vtkStandardNewMacro(myVtkInteractorStyleImage);

2、显示

m_Axialviewer = vtkSmartPointer<vtkImageViewer2>::New();m_Axialviewer->SetInputConnection(this->imageC);m_Axialviewer->SetColorLevel(500);m_Axialviewer->SetColorWindow(2000);m_Axialviewer->SetSlice(40);m_Axialviewer->SetSliceOrientationToXY();m_Axialviewer->GetRenderer()->SetBackground(1, 1, 1);vtkSmartPointer<vtkGenericOpenGLRenderWindow> rwopengl = vtkSmartPointer<vtkGenericOpenGLRenderWindow>::New();m_Axialviewer->SetRenderWindow(rwopengl);ui.m_AxialWidget->SetRenderWindow(m_Axialviewer->GetRenderWindow());vtkSmartPointer<vtkRenderWindowInteractor> rwi = vtkSmartPointer<vtkRenderWindowInteractor>::New(); vtkSmartPointer<myVtkInteractorStyleImage> imagestyle = vtkSmartPointer<myVtkInteractorStyleImage>::New();imagestyle->SetImageViewer(m_Axialviewer);m_Axialviewer->GetRenderWindow()->GetInteractor()->SetInteractorStyle(imagestyle);ui.m_AxialWidget->show();

备注：ui.m_AxialWidget为界面中的QVTKOpenGLNativeWidget组件。

this->m_Axialviewer 为头文件中定义的vtkImageViewer2。