MPR在旋转轴线时图像周围出现灰色边界的解决方法

本文是由Markdown语法编辑器编辑完成。

1. 问题提出：

近日公司开发的一款医学影像处理软件，在MPR显示的界面，当旋转MPR上的轴线时，在图像的后面会出现灰色的背景。这让用户觉得很是不方便，会影响对图像的观察和处理。

如下图所示为，当旋转MPR的轴线时，随着MPR轴线的旋转会有一个背景为灰色的平面跟随着轴线一起旋转。这便是本文需要解决的问题。

2. 解决过程：

从界面显示的效果直观上看，应该是把MPR输出图像的背景设置成黑色就可以了。由于MPR的实现底层是基于VTK，而VTK实现数据渲染是基于管道机制实现的。关键的代码如下：

......
this->m_spPlaneSource=vtkSmartPointer<vtkPlaneSource>::New();

this->m_spResliceAxes= vtkSmartPointer<vtkMatrix4x4>::New();

this->m_spReslice=vtkSmartPointer<vtkImageSlabReslice>::New();
this->m_spReslice->SetBlendModeToMax();
this->m_spReslice->SetSlabThickness(0.0);
this->m_spReslice->SetSlabResolution(1.0);
this->m_spReslice->TransformInputSamplingOff();
this->m_spReslice->AutoCropOutputOn();
this->m_spReslice->SetResliceAxes(m_spResliceAxes);

this->m_spTexture=vtkSmartPointer<vtkTexture>::New();
this->m_spTexture->SetInputConnection(this->m_spReslice->GetOutputPort());
this->m_spTexture->SetInterpolate(1);
this->m_spTexture->SetQualityTo32Bit();
this->m_spTexture->MapColorScalarsThroughLookupTableOff();
this->m_spTexture->SetInterpolate(1);
this->m_spTexture->RepeatOff();
this->m_spTexture->SetLookupTable(GetLookupTable());

vtkPolyDataMapper* texturePlaneMapper = vtkPolyDataMapper::New();
texturePlaneMapper->SetInputConnection(m_spPlaneSource->GetOutputPort());
texturePlaneMapper->SetResolveCoincidentTopologyToPolygonOffset();

this->m_spTexturePlaneActor=vtkSmartPointer<vtkActor>::New();
this->m_spTexturePlaneActor->SetMapper(texturePlaneMapper);
this->m_spTexturePlaneActor->SetTexture(this->m_spTexture);
this->m_spTexturePlaneActor->PickableOn();
texturePlaneMapper->Delete();

GetRenderer()->GetActiveCamera()->ParallelProjectionOn();
GetRenderer()->AddActor(this->m_spTexturePlaneActor);
......

从源代码中可以看出，最后显示在界面中的Actor是m_spTexturePlaneActor。而它的两个接口SetMapper()和SetTexture()是为这个Actor提供数据来源。Mapper的数据来源是PlaneSource，而Texture的数据来源便是vtkImageReslice切割体数据后形成的二维图像数据。

解决思路1：

为了设置PlaneSource的背景色。通过查阅文献及VTK的例子，发现vtkPlaneSource是无法为自己设置背景颜色的。如果要修改vtkPlaneSource的背景颜色，有两种解决方法。第一种，是将vtkPlaneSource设置成为一个独立的vtkActor，然后给这个Actor设置背景颜色；第二种，是vtkPlaneSource经过vtkPolyDataMapper，为Mapper映射一个vtkLookpuTable，通过为lookupTable只设置一种颜色对照，实现颜色的改变。
这两种方法的相关vtk例子源码见以下链接：
http://www.vtk.org/Wiki/VTK/Examples/Cxx/GeometricObjects/Plane
http://www.vtk.org/Wiki/VTK/Examples/Cxx/PolyData/ColorCells

但是，经过测试后，这两种思路均无法解决上述的问题。

解决思路2：

既然修改vtkPlaneSource无法解决这个问题。那么，只能从vtkImageReslice来考虑了。既然是修改背景颜色，那么应该寻找vtkImageReslice中与设置背景(background)这个关键词相关的接口函数。

在vtkImageReslice中与background相关的接口有两个：

//Set the background color (for multi-component images).
1. SetBackgroundColor(double,double,double,double);

//Set background grey level (for single-component images).
2. SetBackgroundLevel(double);

当时只明白了第一个接口的作用是设置背景颜色。于是便想当然地把颜色设置成黑色，即(0,0,0,1)。但是设置完成后，依然没有任何效果。

3. 解决方案：

最终的解决方案，其实就是调用了vtkImageReslice的一个之前没有注意过的接口。
vtkImageReslice->SetBackgroundLevel(double);

而这个传入的参数是关键。这个是一个double型的数值，它应该是vtkImageReslice要切割图像的灰度最小值。这样，MPR的背景的灰阶是始终保持和整个图像的灰度最小值相当的。

对于CT图像来说，灰度的最小值会出现负数，因此这里的double值是一个负数；而对于MRI图像，灰度的最小值为0，因此这里的double应该为0。