当前位置：首页 > news >正文

医学图像（DICOM数据）读取及显示（横断面、冠状面、矢状面、3D显示）为什么用ITK+VTK，单独用ITK或者VTK能实一样功能吗？

news 来源：原创 2025/4/25 11:59:14

在医学图像处理中，结合使用 ITK 和 VTK 是常见的做法，因为它们各自专注于不同的核心功能。以下是逐步解释为何代码中同时使用两者，以及单独使用是否可行的分析：

1. 为什么用ITK处理DICOM数据？

1.1 ITK的DICOM处理优势

复杂DICOM系列支持
ITK的 GDCMSeriesFileNames 类能自动识别目录中的多个DICOM系列（例如不同扫描序列），并通过 GetSeriesUIDs() 获取所有唯一标识符。这在包含多个患者或多次扫描的场景中非常关键。
```
std::vector<std::string> vecSeries = nameGenarator->GetSeriesUIDs();
```
代码中通过遍历所有系列，选择文件最多的一个（假设为最完整的扫描），这需要ITK的灵活API支持。
递归目录搜索
ITK支持递归遍历子目录（SetRecursive(true)），确保分散存储的DICOM文件也能被正确读取。
```
nameGenarator->SetRecursive(true);
```

元数据解析
ITK的 GDCMImageIO 能高效提取DICOM元数据（如患者姓名、扫描参数）：

// 使用ITK读取DICOM系列
itk::SmartPointer<itk::ImageSeriesReader<itk::Image<short, 3>>> seriesReader = itk::ImageSeriesReader<itk::Image<short, 3>>::New();
seriesReader->SetFileNames(vecFiles);  // 设置文件列表
// 配置GDCM图像IO以处理DICOM元数据
itk::SmartPointer<itk::GDCMImageIO> dicomIO = itk::GDCMImageIO::New();
seriesReader->SetImageIO(dicomIO);  // 关联IO对象
seriesReader->Update();  // 执行读取操作，加载图像数据
char name[512];
dicomIO->GetPatientName(name); // 直接获取患者姓名

而直接通过元数据字典访问特定标签（如"0010|0010"）也是ITK的强项：

// 通过元数据字典查找特定标签（如患者姓名）
const itk::MetaDataDictionary& dic = dicomIO->GetMetaDataDictionary();
auto tagItr = dic.Find("0010|0010"); // 查找患者姓名标签

1.2 单独使用VTK读取DICOM的局限性

vtkDICOMImageReader的限制
VTK的DICOM阅读器虽然能读取单一系列，但面对多系列、嵌套目录时，需要手动处理文件分组，缺乏ITK的自动化能力。代码中注释的以下部分暗示了这一点：
```
// vtkDICOMImageReader仅支持单个序列，不支持深度读取
// dicomReader->SetDirectoryName(strFolder.toLocal8Bit().data());
```

2. 为什么需要转换到VTK？

2.1 VTK的可视化优势

三维渲染与交互
VTK专为科学可视化设计，支持体绘制、表面渲染和交互式操作（如调整切片）。
坐标系校正
DICOM图像的坐标系（如患者方向）可能与VTK的显示坐标系不一致。可能需要翻转图像，确保显示方向正确：

2.2 ITK的可视化局限性

ITK的 ImageViewer 类仅提供基础的二维图像显示，缺乏三维交互和高级渲染功能。若需构建包含多平面重建（MPR）和三维模型的医学影像浏览器，必须依赖VTK。

3. 能否单独使用ITK或VTK？

3.1 仅用ITK的场景

可行性
如果仅需图像处理（如分割、滤波）和简单二维显示，ITK可以独立完成任务。例如：
```
// ITK的简单图像显示（无交互）
itk::ImageViewer<ImageType>::Show(seriesReader->GetOutput());
```
局限性
无法实现三维可视化、交互式切片浏览或复杂的UI集成。

3.2 仅用VTK的场景

可行性
对于结构简单的DICOM系列（单层扫描、无嵌套目录），可使用 vtkDICOMImageReader 直接读取：
```
vtkSmartPointer<vtkDICOMImageReader> reader = ...;
reader->SetDirectoryName("path/to/dicom");
reader->Update();
```
局限性
1. 无法自动选择最优系列（需手动处理多系列情况）。
2. 元数据提取功能较弱（如获取患者姓名需手动解析标签）。
3. 缺乏ITK的高级图像处理算法（如配准、分割）。

4. 总结：为什么混合使用ITK和VTK？

库	核心优势	在代码中的作用
ITK	DICOM读取、元数据解析、图像处理	递归搜索DICOM文件，选择最佳系列，提取元数据
VTK	三维可视化、交互式渲染	显示二维切片、三维模型，校正图像方向