DICOM:校准

Question

我对预先使用dicom files.Thanks相对较新。

我有两份同一病人的迪克姆文件，以不同的时间间隔拍摄。它们不是完全相同的维度。

第一个是 cube1 104X163X140，第二个是cube2 107X164X140的级。我想在原点对齐两个立方体，并比较它们。

第一个文件的ImagePositionPatient是：-207.4748，-151.3715 -198.7500

第二个文件的ImagePositionPatient是：-207.4500，-156.3500 -198.7500

两个文件具有相同的ImageOrientationPatient -1 0 0 0 1 0 0

有人能给我举个例子吗？我不知道如何将物理平面映射回图像平面？

提前谢谢你，

灰分

===============================================================

加: 23/2/17

在我的例子中，我使用了以下基于链接的矩阵公式：

立方体的IPP (Sxyz) = -207.4748，-151.3715-198.7500

Xxyz (IOP) = 1,0,0

Yxyz (IOP) = 1,0,0

delta_i = 2.5

delta_j = 2.5

​

​

因此，对于cube1的i= 0: 103和j= 0:162的值，我应该计算Pxyz的值。

下一步是什么？对不起，我看不出这将如何帮助我将这两个立方体与不同的IPP对齐到图像平面？。

对新手的问题很抱歉..。

Answer 1

我没有验证你建立的矩阵。但是，如果计算正确，则可以在卷坐标系(x1、y1、z1)之间进行转换，其中x1 = column、y1 = row和z1 = slice number转换为病人坐标系(x2、y2、z2) --这些坐标定义了患者在毫米波中的位置。

通过转换矩阵，您可以从PCS转换回VCS。

比方说，第1卷:= M1的转换矩阵和第2卷:= M2的转换矩阵PCS -> VCS。然后，您可以将第1卷中的点p1转换为第2卷中相应的点p2，方法是使用M1将其转换为PCS，并使用M2‘(倒置M2)将其转换为第2卷。

通过将M1和M2‘相乘，可以直接计算从volume1到volume2的矩阵。

所以：

p2 = (M1 * M2') * p1