棋盘格用途
在进行图像配准时,棋盘格供用户浏览和查看,观察配准前后两组数据的变化、视觉上可以检验配准的效果,一般情况下,CBCT图像引导系统软件均包含有这类功能(如瓦里安)。
基本流程
假设两组数据分别为image1和image2:
- 先对两组数据进行叠加,有两种方法,一种是利用vtkImageBlend 添加两组数据,一种是将两组数据放在actor里,叠加actor;
- 实现vtkCommand的派生类,用于处理交互事件 ;
- 在派生类交互事件中,根据鼠标点击区域和截取范围,利用vtkExtractVOI 分别提取两组数据的ROI;
- 如果ROI1与ROI2的尺寸大小不一致时,后面会发生报错,因此需要先检测尺寸是否一致,不一致则需要重采样;
- 这里笔者之前采用vtkImageChangeInformation,但是发现其尺寸不好改变,即使发生改变,也不能在XYZ三个方向同时满足指定尺寸,因此采用vtkImageReslice ,将ROI2的尺寸与ROI1一致;
- 最终将两组大小一致的ROI放入vtkImageCheckerboard 生成棋盘格;
- 棋盘格的显示和隐藏是通过vtkImageActor来实现的。
部分代码
vtkImageReslice改变尺寸大小:
double origin[3]; double spacing[3]; int ext[6];
extractVOI1->GetOutput()->GetOrigin(origin);
extractVOI1->GetOutput()->GetSpacing(spacing);
extractVOI1->GetOutput()->GetExtent(ext);
//Reslice current overlay roi to the same size of reference image
vtkSmartPointer<vtkImageReslice> overlayExtentReslice = vtkSmartPointer<vtkImageReslice>::New();
overlayExtentReslice->SetOutputExtent(ext[0], ext[1], ext[2], ext[3], ext[4], ext[5]);
overlayExtentReslice->SetOutputSpacing(spacing[0], spacing[1], spacing[2]);
overlayExtentReslice->SetOutputOrigin(origin[0], origin[1], origin[2]);
overlayExtentReslice->SetInputConnection(extractVOI2->GetOutputPort());
overlayExtentReslice->SetOutputDimensionality(3);
overlayExtentReslice->SetInterpolationModeToLinear();
overlayExtentReslice->AutoCropOutputOff();
overlayExtentReslice->SetResliceAxesOrigin(0, 0, 0);
overlayExtentReslice->Update();
checker->SetInputConnection(0, extractVOI1->GetOutputPort());
checker->SetInputConnection(1, overlayExtentReslice->GetOutputPort());
三视图方向变化:
// Reset checkboard orientation
switch (orientation)
{
case 0:
checker->SetNumberOfDivisions(2, 2, 1);
break;
case 1:
checker->SetNumberOfDivisions(2, 1, 2);
break;
case 2:
checker->SetNumberOfDivisions(1, 2, 2);
break;
}
图例