利用VTK对数据进行可视化的实例

来源：互联网发布：unity服务器端源码编辑：程序博客网时间：2024/04/30 11:05

使用Visualization Toolkit 对数据进行可视化十分方便，程序的结构也正如我们所设计的可视化流程一样，十分简洁明了。下面我们将利用Visualization Toolkit，就可视化应用研究中一个十分热门的领域——计算机断层扫描CT，简单剖析一下利用VisualizationToolkit 进行数据可视化的实现。
我们将要进行处理的CT 数据是一个人体头部的切片数据，共有93 个切片。切片的间距是1.5mm ，每个切片由有12 个灰度级的间距为0.8 毫米的256*256 像素构成。我们打算由这些切片数据恢复出皮肤和骨骼的表面。为此我们首先应选取合适的算法，考虑到我们所处理的数据量是非常庞大的，超过了12 兆比特，而且我们只打算对表面进行重建，所以我们选择基于表面重建的经典算法Marching Cubes 算法。而对于Marching Cubes算法，Visualization Toolkit 中已经有封装好的vtkMarchingCubes 给予支持，这样就进一步简化了我们的工作。

       Marching Cube 算法是Lorensen 等人于1987 提出的，是三维数据场等值面生成的经典算法，是体素单元内等值面抽取技术的代表，并一直沿用至今。
       Marching Cube 算法的基本思想是在数据体中将位于两个相临切片上的8 个相临的体素构成一个立方体（cube），逐个处理数据场中的立方体，分类出与等值面相交的立方体，采用插值计算出等值面与立方体边的交点。根据立方体每一顶点与等值面的相对位置，将等值面与立方体边的交点按一定方式连接生成等值面，作为等值面在该立方体内的一个逼近表示。
4.1 读取数据
       首先，我们要做的事情是读取切片数据，并将其转换为我们的开发工具VisualizationToolkit 所支持的一种数据表达形式；然后根据其物理结构建立起相应的模型，我们给CT数据建立的是比较抽象的等值面模型；最后将物理组件与抽象的模型结合在一起来建立对CT 数据的可视化，以帮助用户正确理解数据。我们所要进行处理的是有结构点阵数据，其拓扑和几何都是隐含知道的，所以我们只需要知道数据的维数、数据源和数据空间。利用Visualization Toolkit 中的vtkVolume16Reader 我们可以很方便的读取切片数据，只需要告诉读取数据对象我们的CT 数据的一些参数，如切片之间的间距、切片上像素之间的间距以及所读取切片的起始段（如从第1 个切片到45 个切片），读取数据的代码如下所示：

vtkVolume16Reader *Reader = vtkVolume16Reader::New(); //建立一个读取对象

Reader->SetDataDimensions(256,256); //设置数据的维数

Reader->SetDataByteOrderToLittleEndian ();

Reader->SetFilePrefix (".. /.. /.. /vtkdata/headsq/quarter"); //设置所读取切片数据文件的路径

Reader->SetImageRange(1, 93); //设置读取切片的起始段

Reader->SetDataSpacing (0.8, 0.8, 1.5); //设置切片之间的间距和像素之间的间距

4.2 提取等值面
下面我们就可以用Marching Cubes 算法对所读取的数据进行处理了。首先利用vtkMarchingCubes 类来提取出某一CT 值的等值面，但这时的等值面其实仍只是一些三角面片，还必须由vtkStripper 类将其拼接起来形成连续的等值面。这样就把读取的原始数据经过处理转换为应用数据，也即由原始的点阵数据转换为多边形数据然后由 vtkPolyDataMapper 将其映射为几何数据，并将其属性赋给窗口中代表它的演员，将结果显示出来。

// 从切片数据中提取出皮肤

vtkMarchingCubes *skinExtractor = vtkMarchingCubes::New(); //建立一个Marching Cubes 算法的对象

skinExtractor->SetInput(Reader->GetOutput()); //获得所读取的CT 数据

skinExtractor->SetValue(0, 500); //提取出CT 值为500 的皮肤

vtkStripper *skinStripper = vtkStripper::New(); //建立三角带对象

skinStripper->SetInput(skinExtractor->GetOutput()); //将生成的三角片连接成三角带

vtkPolyDataMapper *skinMapper = vtkPolyDataMapper::New(); //建立一个数据映射对象

skinMapper->SetInput(skinStripper->GetOutput()); //将三角带映射为几何数据

vtkActor *skin = vtkActor::New(); //建立一个代表皮肤的演员

skin->SetMapper(skinMapper); //获得皮肤几何数据的属性

skin->GetProperty()->SetDiffuseColor(1, .49, .25); //设置皮肤颜色的属性

skin->GetProperty()->SetSpecular(.3); //设置反射率

skin->GetProperty()->SetSpecularPower(20); //设置反射光强度

利用同样的方法，我们也可以提取出骨骼的等值面。只是骨骼的CT 值是1150 左右而已。所以只要在SetValue()方法中将参数设置为1150 就可以了。而且Visualization Toolkit支持多表面重建，所以在实际应用中我们可以设置多个参数值，提取出多个等值面并同时显示出来。在这个应用实例中我们只对皮肤和骨骼地的等值面进行了重建。
4.3 显示结果
通过前面这些工作，我们基本上已经完成了对数据的读取处理映射等步骤，下面我们就要对数据进行显示了。

vtkRenderer *ren = vtkRenderer::New(); //建立绘制者

vtkRenderWindow *renWindow = vtkRenderWindow::New(); //建立绘制窗口

renWindow->AddRenderer(ren); //将绘制者加入绘制窗口

vtkRenderWindowInteractor *iren = vtkRenderWindowInteractor::New(); //对绘制结果进行交互操作

iren->SetRenderWindow (renWindow); // 告诉绘制者将要在绘制窗口中进行显示的演员

ren->AddActor(skin); //皮肤

ren->AddActor(bone); //骨骼

//在我的VTK 5.0.2中，相应的测试程序代码如下：

#include "vtkVolume16Reader.h"
       #include "vtkRenderWindowInteractor.h"
     #include "vtkRenderer.h"
      #include "vtkRenderWindow.h"
     #include "vtkMarchingCubes.h"
      #include "vtkStripper.h"
      #include "vtkActor.h"
      #include "vtkPolyDataMapper.h"
      #include "vtkProperty.h"
//读取RAW文件，提取等值面。

int main()
{
      vtkVolume16Reader *reader=vtkVolume16Reader ::New();
      reader->SetDataDimensions(64,64);
      reader->SetDataByteOrderToLittleEndian();
      reader->SetFilePrefix("D://headsq//quarter");
      reader->SetImageRange(1,93);
      reader->SetDataSpacing(3.2,3.2,1.5);

vtkMarchingCubes *boneExtractor=vtkMarchingCubes::New();
      boneExtractor->SetInput((vtkDataObject *)reader->GetOutput());
      boneExtractor->SetValue(0,500);
      vtkStripper *boneStripper=vtkStripper::New();
      boneStripper->SetInput(boneExtractor->GetOutput());

       vtkPolyDataMapper *boneMapper=vtkPolyDataMapper::New();
       boneMapper->SetInput(boneStripper->GetOutput());

vtkActor *bone=vtkActor::New();
bone->SetMapper(boneMapper);

bone->GetProperty()->SetDiffuseColor(.1,.94,.52);
       bone->GetProperty()->SetSpecular(.3);
       bone->GetProperty()->SetSpecularPower(20);
       vtkRenderer *ren=vtkRenderer::New();
       vtkRenderWindow *renWindow=vtkRenderWindow::New();
       renWindow->AddRenderer(ren);
       vtkRenderWindowInteractor *iren=vtkRenderWindowInteractor::New();
       ren->AddActor(bone);
       iren->Initialize();
       iren->Start();
       reader->Delete();
       iren->Delete();

return 0;
}