Dicom图像方位解析

问题

任意给定一张Dicom图像或者MPR重建图像，怎么解析这个图像的具体方位呢？

前置

1. Dicom图像中通常使用LPS来表示X，Y，Z轴正方向。
2. Dicom图像中的tag(0020,0037)表示图像X轴与Y轴在标准LPS坐标系里的方位向量。

解析

对于任意给定的方向向量（x,y,z）

1. x>0，表示该方向与L同向；x<0，表示该方向与R同向；
2. y>0，表示该方向与P同向；x<0，表示该方向与A同向；
3. z>0，表示该方向与S/H同向；x<0，表示该方向与I/F同向；

另外单位向量有x^2 + y^2 + z^2 = 1，

1. 解析1个方向，考虑分量的绝对值大于sqrt(1/2.0)；
2. 解析2个方向，考虑分量的绝对值大于sqrt(1/3.0)；

char* ComputeOrientation(std::vector<double> vecs)
{char *orientation = new char[4];char *optr = orientation;*optr = '\0';// RAI（F）char x_desc = vecs[0] < 0 ? 'R' : 'L';char y_desc = vecs[1] < 0 ? 'A' : 'P';char z_desc = vecs[2] < 0 ? 'I' : 'S';double abs_x = fabs(vecs[0]);double abs_y = fabs(vecs[1]);double abs_z = fabs(vecs[2]);// double eps = 0.001	// 3个方向// double eps = 0.577;	// 2个方向double eps = 0.707;		// 1个方向for (int i = 0; i < 3; ++i){if (abs_x > eps && abs_x > abs_y && abs_x > abs_z){*optr++ = x_desc;abs_x = 0;}else if (abs_y > eps && abs_y > abs_x && abs_y > abs_z){*optr++ = y_desc;abs_y = 0;}else if (abs_z > eps && abs_z > abs_x && abs_z > abs_y){*optr++ = z_desc;abs_z = 0;}else break;*optr = '\0';}return orientation;
}