基于装置成像坐标的MR图像配准方法、分析方法及其装置制造方法及图纸

本发明专利技术公开了一种基于装置成像坐标的MR图像配准方法及装置，属于核磁共振领域。本发明专利技术用于通过坐标系变化，实现MR结构像与CEST图像之间的跨模态配准，以便于后续在高分辨率的MR结构像上提取感兴趣区域，提高感兴趣区域勾画的精确度，并计算CEST图像上对应区域的CEST量化指标。本发明专利技术能自动快速实现3D

【技术实现步骤摘要】
基于装置成像坐标的MR图像配准方法、分析方法及其装置


[0001]本专利技术属于磁共振成像领域，具体涉及一种MR图像的配准、分析方法及装置。

技术介绍

[0002]磁共振化学交换饱和转移(Chemical exchange saturation transfer,CEST)成像是在磁化转移技术及化学交换理论基础上发展起来的一种新的磁共振成像方法，可产生与常规磁共振成像序列不一样的对比度。其成像原理是利用特定的偏共振饱和脉冲，对特定物质中的可交换质子或分子进行充分的预饱和，这种饱和通过化学交换，进一步影响自由水的信号强度，因此通过检测水的信号，可间接反映这种物质的信息。CEST MRI可以检测内源性和外源性物质，例如蛋白质和多肽中的酰胺质子(
‑
NH)，谷氨酸和肌酸中的胺质子(
‑
NH2)，糖胺聚糖、糖原和葡萄糖中的羟基质子(
‑
OH)等，具有广泛的研究和应用价值。酰胺质子转移(Amideproton transfer,APT)成像是CEST MRI的一种细分技术，通过检测游离蛋白中的酰胺质子本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于装置成像坐标的MR图像配准方法，用于实现MR结构像与CEST图像之间的跨模态配准，其中MR结构像和CEST图像中一个为参考图像，另一个为待配准图像，其特征在于，MR图像配准方法包括：S1、基于分别从参考图像和待配准图像的DICOM信息中提取的旋转变换矩阵、平移变换矩阵和缩放变换矩阵，计算将参考图像由图像坐标系变换到解剖学空间坐标系所需的第一仿射变换矩阵A1，以及将待配准图像由图像坐标系变换到解剖学空间坐标系所需的第二仿射变换矩阵A2；S2、根据第一仿射变换矩阵A1，将参考图像中每个体素在图像坐标系下的坐标[i1，j1，k1]转换为解剖学空间坐标系下的坐标[x，y，z]：S3、基于待配准图像和参考图像的解剖学空间坐标系的统一性，利用第二仿射变换矩阵A2的逆变换，将参考图像中每个体素在解剖学空间坐标系下的坐标[x，y，z]转换为该体素在待配准图像的图像坐标系下的坐标[i2，j2，k2]：S4、针对参考图像中每个体素经过转换后的坐标[i2，j2，k2]，通过插值方法利用待配准图像中该坐标周边的体素灰度值计算该坐标处的灰度值，从而得到与参考图像中每个体素一一对应的待配准图像，实现待配准图像在参考图像上的跨模态配准。2.如权利要求1所述的基于装置成像坐标的MR图像配准方法，其特征在于，所述S1中，待配准图像与参考图像由同一患者在同一次MR扫描时采集得到，且两次序列采集之间无明显头动。3.如权利要求1所述的基于装置成像坐标的MR图像配准方法，其特征在于，所述S1中，参考图像和待配准图像的仿射变换矩阵计算方法相同，步骤如下：S11、从目标图像的DICOM信息中提取六元double数组[R
ix
，R
iy
，R
iz
，R
jx
，R
jy
，R
jz
]形式的(0020,0037)Image Orientation(Patient)字段，再由[R
ix
,R
iy
,R
iz
]和[R
jx
,R
jy
,R
jz
]的叉积计算得到[R
kx
,R
ky
,R
kz
]，从而得到旋转变换矩阵Rotation：S12、从目标图像的DICOM信息中提取三元double数组[T
x
,T
y
,T
z
]形式的(0020,0032)Image Position(Patient)字段，从而得到平移变换矩阵Translation：S13、从目标图像的DICO...

【专利技术属性】
技术研发人员：张祎，汪文祺，吴丹，
申请(专利权)人：浙江大学，
类型：发明
国别省市：