DICOM:对齐

DICOM: alignment

我对提前使用 dicom files.Thanks 比较陌生。

我有 2 个不同时间间隔拍摄的同一患者的 dicom 文件。 它们的尺寸不完全相同。

第一个是:立方体1的尺寸 104X163X140,第二个是立方体2的尺寸[=64] =] 107X164X140。我想在原点对齐两个立方体并进行比较。

第一个文件的ImagePositionPatient为:[-207.4748, -151.3715 -198.7500]

第二个文件的ImagePositionPatient为:[-207.4500, -156.3500 -198.7500]

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

有人可以给我举个例子吗?我不确定如何将物理平面映射回图像平面?

非常感谢,

小灰

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

添加:23/2/17

我根据 link 使用了下面的矩阵公式,在我的例子中:

立方体 1 的 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 的两个立方体与图像平面对齐?

抱歉新手问题...

我没有验证你构建的矩阵。但如果计算正确,您可以在体积坐标系 (VCS) (x1,y1,z1) 之间转换,其中 x1 = 列,y1 = 行,z1 = 切片编号到患者坐标系 (PCS) (x2, y2, z2) - 这些坐标以毫米为单位定义了患者体内的点。

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

假设,体积 1 的变换矩阵:= M1 和体积 2 的变换矩阵 PCS -> VCS := M2。然后,您可以将点 p1 从体积 1 转换为体积 2 中的对应点 p2,方法是使用 M1 将其转换为 PCS,并使用 M2'(倒置的 M2)从 PCS 转换为体积 2。

将M1和M2'相乘,可以计算出一个从volume1直接转换到volume2的矩阵。

所以:

p2 = (M1 * M2') * p1