VTK 实例23:DCM医学图像信息的访问(图像维数、原点坐标和像素间隔)

发布时间 2023-08-15 11:29:17作者: 一杯清酒邀明月
 1 #include "vtkAutoInit.h" 
 2 VTK_MODULE_INIT(vtkRenderingOpenGL2);
 3 VTK_MODULE_INIT(vtkInteractionStyle);
 4 
 5 #include <vtkSmartPointer.h>
 6 #include <vtkImageViewer2.h>
 7 #include <vtkRenderWindow.h>
 8 #include <vtkRenderWindowInteractor.h>
 9 #include <vtkRenderer.h>
10 #include <vtkImageReader2Factory.h>
11 #include <vtkImageReader2.h>
12 //#include <vtkBMPReader.h>
13 #include <vtkImageData.h>
14 #include <vtkDICOMImageReader.h>//DCM医学文件读取类
15 
16 //测试图像:../data/lena.bmp
17 int main(int argc, char* argv[])
18 {
19     /*vtkSmartPointer<vtkBMPReader> reader =
20         vtkSmartPointer<vtkBMPReader>::New();*/
21 
22     //实例化一个对象dcmReader
23     vtkDICOMImageReader* reader = vtkDICOMImageReader::New();
24     reader->SetFileName("C:\\Users\\Administrator\\Desktop\\VTK2\\hellovtk\\CT\\123.dcm");//读单张切片
25 
26     //reader->SetFileName("lena.bmp");//读取图片
27     reader->Update();
28 
29     int dims[3];
30     reader->GetOutput()->GetDimensions(dims);
31     std::cout<<"图像维数:" <<dims[0]<<" "<<dims[1]<<" "<<dims[2]<<std::endl;
32 
33     double origin[3];
34     reader->GetOutput()->GetOrigin(origin);
35     std::cout<<"图像原点:" <<origin[0]<<" "<<origin[1]<<" "<<origin[2]<<std::endl;
36 
37     double spaceing[3];
38     reader->GetOutput()->GetSpacing(spaceing);
39     std::cout<<"像素间隔:" <<spaceing[0]<<" "<<spaceing[1]<<" "<<spaceing[2]<<std::endl;
40 
41     vtkSmartPointer<vtkImageViewer2> imageViewer =
42         vtkSmartPointer<vtkImageViewer2>::New();
43     imageViewer->SetInputConnection(reader->GetOutputPort());
44 
45     vtkSmartPointer<vtkRenderWindowInteractor> renderWindowInteractor =
46         vtkSmartPointer<vtkRenderWindowInteractor>::New();
47     imageViewer->SetupInteractor(renderWindowInteractor);
48     imageViewer->Render();
49     imageViewer->GetRenderer()->ResetCamera();
50     imageViewer->Render();
51 
52     imageViewer->GetRenderer()->SetBackground(1.0, 1.0, 1.0);
53     imageViewer->SetSize(640, 480);
54     imageViewer->GetRenderWindow()->SetWindowName("GetImageInformationExample");
55 
56     renderWindowInteractor->Start();
57 
58     return EXIT_SUCCESS;
59 }