在 DICOM 中旋转图像方向（患者）答案

【问题标题】：Rotate Image Orientation(Patient) in DICOM在 DICOM 中旋转图像方向（患者）
【发布时间】：2017-08-15 13:32:25
【问题描述】：

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

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

任何形式的帮助都将不胜感激，谢谢。

西蒙

【问题讨论】：

您的问题到底是什么？寻找变换矩阵来计算两个堆栈中像素的世界（=患者）坐标系坐标？或者计算将位置向量从一个堆栈转换到另一个堆栈的矩阵？还是两者兼而有之？
我想计算将位置向量从一个堆栈转换到另一个堆栈的矩阵。

标签： dicom

【解决方案1】：

This page 解释了如何从 DICOM 标头中的几何信息形成变换矩阵。这些变换矩阵用于从体积坐标系（pixel-x、pixel-y、切片编号）变换到患者/世界坐标系（以毫米为单位）。

从第 1 卷转换到第 2 卷的基本思想是从第 1 卷转换到患者坐标系，从患者坐标转换到第 2 卷坐标系。将这两个矩阵相乘得到矩阵，可以直接从第 1 卷转换到第 2 卷。

注意：显然，不能保证 v1 中的每个坐标都匹配 v2 中的坐标，即堆栈可能具有不同的大小和/或位置。

所以你有：

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

然后

M1 * (M2^(-1)) 是将位置向量从体积 1 转换到体积 2 的矩阵（输入和输出是像素-x，像素-y，切片数）

和

M2 * (M1^(-1)) 是将位置向量从体积 1 转换为体积 2 的矩阵（输入和输出为像素-x，像素-y，切片数）

【讨论】：

感谢您的回答！也许我不够清楚，但我所拥有的只是两个堆栈的图像方向，即行和列方向相对于世界坐标系的方向余弦。我如何理解为了将位置向量从第一组转换为第二组要应用哪些旋转？
那是我在问你已经拥有什么以及你需要什么。我已经相应地更新了我的答案。