1 #include "itkImage.h"
2
3 // Function to simulate getting mouse click fr an image
4 static itk::Image< unsigned short, 3 >::IndexType GetIndexFrMouseClick()
5
12 int main(int, char *[])
13 }; //Size along
19 const ImageType::IndexType start = }; // First index on
20
21 ImageType::RegionType region;
22 region.SetSize( size );
23 region.SetIndex( start );
24
25 image>SetRegions( region );
26 image>Allocate(true); // initialize buffer to zero
27
28 //图像类中处理原点和间距的方法
29 //创建一个和图像数据类型相一致的数列spacing
30 ImageType::SpacingType spacing;
31 //设定X、Y、Z方向间距
32 spacing[0] = 0.33; // spacing along X X方向上相邻两像素中心的间距
33 spacing[1] = 0.33; // spacing along Y Y方向上相邻两像素中心的间距
34 spacing[2] = 1.20; // spacing along Z Z方向上相邻两像素中心的间距
35 //使用 SetSpacing( ) 方法指向数列spacing
36 image>SetSpacing( spacing );
37 //使用 GetSpacing( ) 方法可以从图像中得到间距信息, const 表示数列是不可修改的
38 const ImageType::SpacingType& sp = image>GetSpacing();
39 //输出读取到的图像X、Y、Z方向的间距信息
40 std::cout << "Spacing = ";
41 std::cout << sp[0] << ", " << sp[1] << ", " << sp[2] << std::endl;
42
43 //初始化图像原点的变量newOrigin的创建和分配
44 ImageType::PointType newOrigin;
45 newOrigin.Fill(0.0);
46 image>SetOrigin( newOrigin );
47 // GetOrigin( ) 方法可以从图像中读取原点
48 const ImageType::PointType & origin = image>GetOrigin();
49 //输出读取到图像的原点坐标
50 std::cout << "Origin = ";
51 std::cout << origin[0] << ", "
52 << origin[1] << ", "
53 << origin[2] << std::endl;
54
55
56 ImageType::DirectionType direction;
57 direction.SetIdentity();
58 image>SetDirection( direction );
59
60 const ImageType::DirectionType& direct = image>GetDirection();
61
62 std::cout << "Direction = " << std::endl;
std::cout << direct << std::endl;
64 //将物理空间映射到读取最近像素内容的图像index中
65 //声明一个 itk::Point 类型。这个 Point 类型在用来表示坐标的类型和空间大小之上模块化
66 typedef itk::Point< double, ImageType::ImageDimension > PointType;
67
68 PointType point;
69 point[0] = 1.45; // x coordinate
70 point[1] = 7.21; // y coordinate
71 point[2] = 9.28; // z coordinate
72 //图像类型中定义的 IndexType 来对 index 对象进行实例化
73 ImageType::IndexType pixelIndex;
74
75 // Point 到 index 的映射和访问图像像素数据的像素 index 的用法
76 const bool isInside =
77 image>TransformPhysicalPointToIndex( point, pixelIndex );
78 if ( isInside )
79
84
85 const ImageType::IndexType LeftEyeIndex = GetIndexFrMouseClick();
86 ImageType::PointType LeftEyePoint;
87 image>TransformIndexToPhysicalPoint(LeftEyeIndex,LeftEyePoint);
88
89 std::cout << "===========================================" << std::endl;
90 std::cout << "The Left Eye Location is " << LeftEyePoint << std::endl;
91
92
93 typedef itk::Matrix<double, Dimension, Dimension> MatrixType;
94 MatrixType SpacingMatrix;
95 SpacingMatrix.Fill( 0.0F );
96
97 const ImageType::SpacingType & ImageSpacing = image>GetSpacing();
98 SpacingMatrix( 0,0 ) = ImageSpacing[0];
99 SpacingMatrix( 1,1 ) = ImageSpacing[1];
100 SpacingMatrix( 2,2 ) = ImageSpacing[2];
101
102 const ImageType::DirectionType & ImageDirectionCosines =
103 image>GetDirection();
104 const ImageType::PointType &ImageOrigin = image>GetOrigin();
105
106 typedef itk::Vector< double, Dimension > VectorType;
107 VectorType LeftEyeIndexVector;
108 LeftEyeIndexVector[0]= LeftEyeIndex[0];
109 LeftEyeIndexVector[1]= LeftEyeIndex[1];
110 LeftEyeIndexVector[2]= LeftEyeIndex[2];
111
112 ImageType::PointType LeftEyePointByHand =
113 ImageOrigin + ImageDirectionCosines * SpacingMatrix * LeftEyeIndexVector;
114 // Software Guide : EndCodeSnippet
115
116 std::cout << "===========================================" << std::endl;
117 std::cout << "Spacing:: " << std::endl << SpacingMatrix << std::endl;
118 std::cout << "===========================================" << std::endl;
119 std::cout << "DirectionCosines:: " << std::endl << ImageDirectionCosines << std::endl;
120 std::cout << "===========================================" << std::endl;
121 std::cout << "Origin:: " << std::endl << ImageOrigin << std::endl;
122 std::cout << "===========================================" << std::endl;
123 std::cout << "The Left Eye Location is " << LeftEyePointByHand << std::endl;
124
125 if ( (LeftEyePointByHand LeftEyePoint).GetNorm() < 0.01F )
126
132
133 return EXIT_SUCCESS;
134 }
上一篇：VTK 实例55：连通区域分析
下一篇：VTK 实例64：体绘制管线（光线投影法体绘制）
ITK
