在 DICOM 中旋转图像方向(患者)

Rotate Image Orientation(Patient) in DICOM

我从 MRI 采集中提取了一个 3D 表面,描述该表面的点的坐标(我相信)是相对于该系列第一张图像的参考系统(我的意思是原点对应于图像位置(患者)和轴方向到图像方向(患者))。 我有另一组具有不同图像位置(患者)和不同图像方向(患者)的图像;我想旋转和平移从第一组图像中提取的表面,以使其与第二组图像完全匹配。

我很难找到可以完成这项工作的正确 4x4 矩阵,一旦找到它,我就知道如何将它应用到我的表面上。

如有任何帮助,我们将不胜感激,谢谢。

西蒙

This page 解释了如何从 DICOM headers 中的几何信息形成转换矩阵。这些变换矩阵用于从体积坐标系(pixel-x、pixel-y、切片数)变换到patient/world坐标系(单位为毫米)。

从体积 1 到体积 2 的基本思想是从体积 1 到患者坐标系,从患者坐标到体积 2 坐标系。将两个矩阵相乘得到的矩阵直接从体积 1 转换到体积 2。

注意:显然,不能保证 v1 中的每个坐标都与 v2 中的坐标匹配,即堆栈可能具有不同的大小 and/or 位置。

所以你有:

M1 - 从体积 1 转换到世界坐标系的矩阵和 M2 - 从体积 2 转换到世界坐标系的矩阵

然后

M1 * (M2^(-1)) 是将位置向量从体积 1 转换为体积 2 的矩阵(输入和输出为 pixel-x、pixel-y、切片编号)

M2 * (M1^(-1)) 是将位置向量从体积 1 转换为体积 2 的矩阵(输入和输出为 pixel-x、pixel-y、切片编号)