C++实现itk读取dicom文件并转成vtk数据

以下是一个简单的示例，展示如何使用ITK和VTK库读取DICOM文件并将其转换为VTK数据：

#include <vtkSmartPointer.h>
#include <vtkDICOMImageReader.h>
#include <vtkImageData.h>
#include <vtkImageViewer2.h>

#include <itkImage.h>
#include <itkImageFileReader.h>

int main(int argc, char* argv[])
{
    // 读取DICOM文件并将其转换为ITK图像
    typedef itk::Image<short, 3> ImageType;
    typedef itk::ImageFileReader<ImageType> ReaderType;
    ReaderType::Pointer reader = ReaderType::New();
    reader->SetFileName("example.dcm");
    reader->Update();
    ImageType::Pointer itkImage = reader->GetOutput();

    // 将ITK图像转换为VTK图像
    vtkSmartPointer<vtkImageData> vtkImage = vtkSmartPointer<vtkImageData>::New();
    vtkImage->SetDimensions(itkImage->GetLargestPossibleRegion().GetSize()[0],
                            itkImage->GetLargestPossibleRegion().GetSize()[1],
                            itkImage->GetLargestPossibleRegion().GetSize()[2]);
    vtkImage->SetSpacing(itkImage->GetSpacing()[0], itkImage->GetSpacing()[1], itkImage->GetSpacing()[2]);
    vtkImage->SetOrigin(itkImage->GetOrigin()[0], itkImage->GetOrigin()[1], itkImage->GetOrigin()[2]);
    vtkImage->AllocateScalars(VTK_SHORT, 1);
    vtkIdType count = 0;
    for (int z = 0; z < itkImage->GetLargestPossibleRegion().GetSize()[2]; z++) {
        for (int y = 0; y < itkImage->GetLargestPossibleRegion().GetSize()[1]; y++) {
            for (int x = 0; x < itkImage->GetLargestPossibleRegion().GetSize()[0]; x++) {
                vtkImage->SetScalarComponentFromDouble(x, y, z, 0, itkImage->GetPixel({ x, y, z }));
                count++;
            }
        }
    }

    // 显示VTK图像
    vtkSmartPointer<vtkImageViewer2> viewer = vtkSmartPointer<vtkImageViewer2>::New();
    viewer->SetInputData(vtkImage);
    viewer->Render();
    viewer->GetRenderWindow()->SetWindowName("DICOM Viewer");
    viewer->GetRenderWindow()->Render();
    viewer->Start();

    return EXIT_SUCCESS;
}

需要注意的是，这只是一个简单的示例，仅适用于读取和显示单张DICOM文件。在实际应用中，可能需要处理多个DICOM文件、根据DICOM标签进行图像处理等等。