找到在同一扫描会话中获取的两个MR数据集的相同切片位置的坐标(mm)

时间:2016-05-03 05:35:21

标签: dicom image-registration mri

我有两次MR采集,其中第一次是3D采集(1x1x1 mm3),第二次是2D采集(2.24 x 2.24 x 5.00 mm,轴向切片)。高分辨率数据集是一个全头3D采集,如果沿轴向重新定向,则会产生176个切片(在采集中最初是定向的)。 2D采集仅包含3个轴向切片,这些切片被选择用于靶向大脑中的特定区域并且被连续采集。

是否有可能知道3D高分辨率数据集中的哪些切片准确对应于2D数据集中的3个切片,假设主体在扫描之间没有移动并且这些数据集是在同一扫描会话中获取的? 我正在研究这两个数据集的dicoms并试图使用ImageOrientationPatient和ImagePositionPatient标签来尝试准确找出2D扫描中切片相对于磁体等中心的坐标。这样我可以告诉2D采集中第一个切片相对于磁体等中心的精确坐标,假设两个扫描共享相同的原点,那么我可以准确地知道坐标对应的3D扫描中的哪个轴向切片? 我面临的问题是ImageOrientationPatient向量对于两次采集都是不同的,因为高位是通过矢状方向特定来获得的(虽然3d并且可以在任何方向上获得切片),而2D数据被特别地获取为轴向切片。

有经验使用dicom处理的人能否对我如何链接两次扫描有所了解?由于它是相同的扫描会话,我假设这两次采集的参考位置应该是相同的。这是对的吗?

1 个答案:

答案 0 :(得分:3)

只有在一次扫描中拍摄所有图像时,ImagePositionPatient才能可靠地引用相同的坐标系。不幸的是,没有对患者进行校准"这可确保特定坐标始终引用同一患者的相同位置。

假设您的扫描就是这种情况,那么任务非常简单。每次扫描的完整几何形状由以下定义:

  • ImagePositionPatient(0020,0032) - 每个切片左上角像素的坐标
  • ImageOrientationPatient(0020,0037) - 每个切片的方向矢量,即像素行和列对齐的轴

  • PixelSpacing(0028,0030) - 每个像素的高度和宽度(是的,y维度首先在这里!)

现在您已准确知道两次扫描中每个切片的每个像素的位置,剩下的任务是在3D扫描的尺寸中表示2D扫描的左上角像素。我会使用坐标转换来做到这一点,但其他方法也可以。