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
我对提前使用 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