05-VTK在图像处理中的应用(4)

最新推荐文章于 2023-12-07 09:57:25 发布

DolingStudio

最新推荐文章于 2023-12-07 09:57:25 发布

阅读量1.5w

点赞数 6

分类专栏： VTK系列教程 vtk应用开发文章标签： VTK 命令模式图像处理观察者模式

本文链接：https://blog.csdn.net/www_doling_net/article/details/8551523

版权

本文详细介绍了VTK库在图像处理中如何进行感兴趣区域（VOI）提取和三维图像切片提取。通过vtkExtractVOI类提取图像的子区域，并展示了如何利用vtkImageReslice类提取不同平面的切片，包括平行于XY、YZ、XZ平面的切片以及任意切片的设置方法。此外，文章还提供了一个交互式程序示例，允许用户通过鼠标切换切片。

摘要由CSDN通过智能技术生成

5.7 区域提取

5.7.1 提取感兴趣区域

感兴趣区域（Volum of Interest）是指图像内部的一个子区域。在VTK中vtkExtractVOI类实现由用户指定的区域范围提取图像的子图像。该Filter的输入和输出都是一个vtkImageData，因此其结果可以直接作为图像保存。

1: vtkSmartPointer<vtkBMPReader> reader =

2: vtkSmartPointer<vtkBMPReader>::New();

3: reader->SetFileName ( "lena.bmp" );

4: reader->Update();

6: int dims[3];

7: reader->GetOutput()->GetDimensions(dims);

9: vtkSmartPointer<vtkExtractVOI> extractVOI =

10: vtkSmartPointer<vtkExtractVOI>::New();

11: extractVOI->SetInputConnection(reader->GetOutputPort());

12: extractVOI->SetVOI(dims[0]/4.,3.*dims[0]/4.,dims[1]/4.,3.*dims[1]/4., 0, 0);

13: extractVOI->Update();

上例代码实现了提取一副图像的子区域。首先读取一个图像，并获取图像的维数。然后定义vtkExtractVOI对象，该对象接收两个输入一个是图像数据，第二个是区域大小。设置区域大小的函数原型：

void SetVOI(int _arg1, int _arg2, int _arg3, int _arg4, int _arg5, int _arg6)

void SetVOI(int _arg[])

其参数是提取的区域各个方向的大小，共6个参数，依次表示x方向最小值，x方向最大值，y方向最小值，y方向最大值，z方向最小值和z方向最大值。上例中由于读取的是二维图像，因此z方向的区域为[0,0]，而在x方向范围为[ dims[0]/4 , 3*dims[0]/4 ]，y方向范围为[ dims[1]/4 , 3*dims[1]/4 ]，即提取图像原图中间1/4图像。执行结果如下：

图5.18 提取感兴趣区域

5.7.2 三维图像切片提取

切片是指三维图像中的一个切面对应的图像。切面可以是过图像内部一点且平行于XY、YZ、XZ平面的平面，也可以是任意的过三维图像内部一点任意方向的平面。通过提取切片可以方便的浏览和分析图像内部组织结构，是医学图像浏览软件中的一个重要的功能。在VTK中vtkImageReslice类实现图像切片提取功能。下面首先看一段切片提取的代码。

1: vtkSmartPointer<vtkMetaImageReader> reader =

2: vtkSmartPointer<vtkMetaImageReader>::New();

3: reader->SetFileName ( " brain.mhd" );

4: reader->Update();

6: int extent[6];

7: double spacing[3];

8: double origin[3];

10: reader->GetOutput()->GetExtent(extent);

11: reader->GetOutput()->GetSpacing(spacing);

12: reader->GetOutput()->GetOrigin(origin);

13:

14: double center[3];