背景技术
在磁共振成像中,如果信号源沿梯度场方向流动,那么所测到的信号相位就会与流体的流速以及梯度场脉冲形状相关联。基于这种原理,通过改变梯度场脉冲形状从而引起的磁共振信号相位改变,可以测量流体的流速和流量。有关这种测量流速流量的磁共振技术详情,可参见文献1:[Pelc NJ,“Flow quantification andanalysis methods”,Magn Reson Imaging Clin N Am3:413-424(1995)]。
简单地说,采用二种用以测量流速的不同梯度场脉冲形状,所采集的磁共振信号可作如下描述:
式中M代表流体的磁共振信号幅度,ν代表流速,β代表主磁场或射频场不均匀等导致的系统相位,m1(1)和m1(2)分别为测量中用到的二种梯度场脉冲形状G(t)的一阶矩:
式中γ代表原子核的磁旋比常数,TE代表回波时间。
流速即可以简单地从二次测量的信号相位差求得:
ν=Δφ/[m1(2)-m1(1)] [4]
Δφ=arg {S1×conj(S2)} [5]
式中arg{…}代表求复数相位;conj(…)代表复数共轭。
当有脂肪信号存在的情况下,以上式[1]和式[2]必须考虑脂肪信号相应改作:
式中W代表流体水信号幅度,F代表脂肪信号的幅度,α代表脂肪信号化学位移引起的相位进动角。从式[6]和[7]可以看出,由于脂肪信号的存在,流速ν再也不能简单地从S1和S2的相位差精确求得,否则所测得的流速存在误差。
美国加州洛杉矶大学的Matthew J.Middione和Daniel B.Ennis对脂肪信号在采用磁共振相位差测量流速中所带来的误差作了详细的分析,参见文献2:[Middione,M.J.and Ennis,D.B.(2012),Chemical shift-induced phase errorsin phase-contrast MRI.Magn Reson Med.doi:10.1002/mrm.24262]。虽然他们也提出了通过对回波时间(TE)和读出带宽的适当选择来减小误差的方法,但还是不能完全消除这种误差。因而有必要开发不受脂肪信号干扰的磁共振相位反差技术,以用于医疗诊断中能精确测量血液流速和流量。
发明内容
本发明所要解决的技术问题是提供一种相位反差磁共振流速、流量测量方法及其测量装置,不受脂肪信号干扰,能够精确测量血液流速和流量。
本发明为解决上述技术问题而采用的技术方案是提供一种相位反差磁共振流速测量方法,包括如下步骤:a)对成像序列分别进行三种不同的流速编码;b)由所述成像序列获取两组混合k空间数据k12和k13,第一组混合k空间数据k12的奇数相位编码线数据对应k空间数据k1,第一组混合k空间数据k12的偶数相位编码线数据对应k空间数据k2,第二组混合k空间数据k13的奇数相位编码线数据对应k空间数据k1,第二组混合k空间数据k13的偶数相位编码线数据对应k空间数据k3,所述混合k空间数据k12和k13的相位编码线步长为正常k空间数据相位编码线步长的一半,所述k空间数据k1、k2和k3分别对应一种步骤a)中所述流速编码;c)将所述第一组混合k空间数据k12和所述第二组混合k空间数据k13分别变换到图像域得到两幅图像,每幅图像均包含AC图和DC图;d)由所述两幅图像中的AC图利用相位反差计算流速。
上述的相位反差磁共振流速测量方法,其中,所述步骤b)由所述成像序列获取两组混合k空间数据k12和k13包括如下步骤:由所述成像序列采集得到k空间数据k1、k2和k3;将所述k空间数据k1和所述k空间数据k2混合得到第一组混合k空间数据k12,将所述k空间数据k1和所述k空间数据k3混合得到第二组混合k空间数据k13。
上述的相位反差磁共振流速测量方法,其中,所述步骤b)中由所述成像序列获取两组混合k空间数据k12和k13包括如下步骤:由所述成像序列交替采集得到第一组混合k空间数据k12;由所述成像序列交替采集得到第二组混合k空间数据k13。
上述的相位反差磁共振流速测量方法,其中,所述采集为:首先进行并行采集,采集完毕后进行k空间并行加速数据重建。
上述的相位反差磁共振流速测量方法,其中,所述采集为:首先进行并行采集,采集完毕后进行图像域并行加速数据重建,并将重建后的图像数据变换回k空间数据。
上述的相位反差磁共振流速测量方法,其中,所述采集为采集完整数据。
上述的相位反差磁共振流速测量方法,其中,所述步骤a)中三种不同的流速编码为通过梯度场脉冲形状的改变而实现,所述三种流速编码分别为m1(1)、 m1(2)和m1(3),其对应为梯度场一阶矩,所述m1(1)、 m1(2)和m1(3)分布在2π/νmax范围内,νmax为所要测量的最大流速。
上述的相位反差磁共振流速测量方法,其中,所述步骤d)中由所述两幅图像中的AC图利用相位反差计算流速的公式如下:
ν=arg{IAC(2)×conj[IAC(1)]}/[m1(3)-m1(2)];
其中,W代表流体水信号幅度,ν代表流速,β代表主磁场或射频场不均匀导致的系统相位,arg{…}代表求复数相位,conj[…]代表复数共轭,IAC(1)代表第一幅混合k空间对应的AC图,IAC(2)代表第二幅混合k空间对应的AC图。
上述的相位反差磁共振流速测量方法,其中,所述步骤a)中成像序列为梯度回波序列、自旋回波序列、平面回波序列或快速自旋回波序列。
本发明为解决上述技术问题还提供一种相位反差磁共振流量测量方法,先采用上述的相位反差磁共振流速测量方法获得流速,然后在预设的流量测定范围内对所述流速积分得到流量。
本发明为解决上述技术问题还提供一种相位反差磁共振流速测量装置,包括:对成像序列分别进行三种不同的流速编码的序列产生装置;由所述成像序列获取两组混合k空间数据k12和k13的数据处理装置,第一组混合k空间数据k12的奇数相位编码线数据对应k空间数据k1,第一组混合k空间数据k12的偶数相位编码线数据对应k空间数据k2,第二组混合k空间数据k13的奇数相位编码线数据对应k空间数据k1,第二组混合k空间数据k13的偶数相位编码线数据对应k空间数据k3,所述混合k空间数据k12和k13的相位编码线步长为正常k空间数据相位编码线步长的一半,所述k空间数据k1、k2和k3分别对应一种序列产生装置中所述流速编码;将所述第一组混合k空间数据k12和所述第二组混合k空间数据k13分别变换到图像域形成两幅图像的第一图像变换装置,每幅图像均包含AC图和DC图;以及,根据所述两幅图像中的AC图利用相位反差计算流速的计算装置。
上述的相位反差磁共振流速测量装置,其中,所述数据处理装置包括采集装置和数据处理子装置,所述采集装置为由所述成像序列采集得到k空间数据k1、k2和k3的第一采集装置,所述数据处理子装置适于将所述k空间数据k1和所述k空间数据k2混合得到第一组混合k空间数据k12,将所述k空间数据k1和所述k空间数据k3混合得到第二组混合k空间数据k13。
上述的相位反差磁共振流速测量装置,其中,所述数据处理装置包括采集装置,所述采集装置为由所述成像序列交替采集得到第一组混合k空间数据k12,由所述成像序列交替采集得到第二组混合k空间数据k13的第二采集装置。
上述的相位反差磁共振流速测量装置,其中,所述采集装置包括并行采集装置和k空间数据重建装置,所述并行采集装置适于进行并行采集,所述k空间数据重建装置适于在采集完毕后进行k空间并行加速数据重建得到三组k空间数据k1、k2和k3。
上述的相位反差磁共振流速测量装置,其中,所述采集装置包括并行采集装置、图像域数据重建装置和第二图像变换装置,所述并行采集装置适于进行并行采集,所述图像域数据重建装置适于在采集完毕后进行图像域并行加速数据重建,所述第二图像变换装置适于将重建后图像变换回k空间。
上述的相位反差磁共振流速测量装置,其中,所述采集装置为全采集装置,其数据采集方式为采集完整数据。
本发明为解决上述技术问题还提供一种相位反差磁共振流量测量装置,包括上述的相位反差磁共振流速测量装置;以及根据预设的流量测定范围对所述流速进行积分的计算装置。
本发明对比现有技术有如下的有益效果:本发明提供的相位反差磁共振流速、流量测量方法及其测量装置,通过在两个混合k空间交替采集对应于三种不同流速编码的磁共振信号,利用源于k空间线与线之间交变的动态信号获取相位差图,从而不受脂肪产生的稳态信号的干扰,能够精确测量血液流速和流量。
具体实施方式
下面结合附图和实施例对本发明作进一步的描述。
图1为本发明相位反差磁共振流速测量流程示意图。
请参见图1,本发明提供的相位反差磁共振流速测量方法包括如下步骤:
步骤S101:对成像序列分别进行三种不同的流速编码。
所述成像序列可以是任何其他与流速编码相结合的成像序列,如梯度回波序列、自旋回波序列、平面回波序列或快速自旋回波序列。以梯度回波序列为例,除选择用于流速编码的梯度场脉冲外,其余部分与传统的梯度回波成像序列相同。三种不同的流速编码为通过梯度场脉冲形状的改变而实现,分别为m1(1)、 m1(2)和m1(3),其对应为梯度场一阶矩,所述m1(1)、m1(2)和m1(3)分布在2π/νmax范围内,即流速编码可以在任意[-π/νmax,π/νmax]范围内取值,当取值在其范围外时,由于其函数周期性同样可将该值转换为范围内的数值,νmax为所要测量的最大流速。所述m1(2)和m1(3)最好对称分布在m1(1)的两侧,以保证采集到的磁共振信号之间的相位差达到最大;
如:
m1(1)=0, max为所要测量的最大流速。
关于流速编码详情可参阅文献3:[Bittoun J,Jolivet O,Herment A,Itti E,Durand E,Mousseaux E and Tasu JP,“Multidimensional MR mapping of multiplecomponents of velocity and acceleration by Fourier phase encoding with a smallnumber of encoding steps”,Magn Reson Med44:723-730(2000)]。
步骤S102,由所述成像序列获取两组混合k空间数据k12和k13,第一组混合k空间数据k12的奇数相位编码线数据对应k空间数据k1,第一组混合k空间数据k12的偶数相位编码线数据对应k空间数据k2,第二组混合k空间数据k13的奇数相位编码线数据对应k空间数据k1,第二组混合k空间数据k13的偶数相位编码线数据对应k空间数据k3,所述混合k空间数据k12和k13的相位编码线步长为正常k空间数据相位编码线步长的一半,所述k空间数据k1、k2和k3分别对应一种步骤S101中所述流速编码。
获取所述两组混合k空间数据k12和k13的过程,在数据的采集顺序上优选的可以分为两种。第一种采集顺序为首先在三种不同的流速编码m1(1)、m1(2)和m1(3)下分别采集得到三组k空间数据k1、k2和k3,所述的流速编码各对应一组k空间数据,之后将k空间数据k1和k2混合得到第一组混合k空间数据k12,将k空间数据k2和k3混合得到第二组混合k空间数据k13。混合k空间示意图可参见图2,如图2所示,图中三种不同的k空间数据线分别为三种不同的流速编码:m1(1)、m1(2)和m1(3)所对应的k空间数据(即k1、k2和k3)。其中k1、k2和k3所对应的相位编码线分别用细线、粗线和虚线表示。图2a为第一组混合k空间数据k12,其为k空间数据k1和k2的混合,其混合过程为:将k空间数据k1保持不变,将k2沿相位编码方向下移半个编码步长,之后插入k空间数据k1的k空间中;图2b为第二组混合k空间数据k13,其为k空间数据k1和k3的混合,其混合过程为:将k空间数据k1保持不变,将k3沿相位编码方向下移半个编码步长,之后插入k空间数据k1的k空间中。对于混合k空间数据k12,其奇数相位编码线数据(即细线数据)为原k1的数据,偶数相位编码线数据(即粗线数据)为原k2的数据;对于混合k空间数据k13同理,其奇数相位编码线数据(即细线数据)为原k1的数据,偶数相位编码线数据(即虚线数据)为原k3的数据。
优选的第二种采集顺序为交替采集得到混合k空间数据k12和混合k空间数据k13。这里的交替采集方式指的是在成像过程中采集一条k1相位编码线数据,再采集一条k2相位编码线数据,再采集一条k1相位编码线数据,再采集一条k2相位编码线数据,如此交替反复的方式。在流速编码m1(1)和m1(2)所对应的成像序列下,交替采集得到第一组混合k空间数据k12,在流速编码m1(1)和m1(3)所对应的成像序列下,交替采集得到第二组混合k空间数据k13。其中,第一组混合k空间数据k12的奇数相位编码线数据(即细线数据)对应所述k空间数据k1,所述第一组混合k空间数据k12的偶数相位编码线数据(即粗线数据)对应所述k空间数据k2;同理,第二组混合k空间数据k13的奇数相位编码线数据(即细线数据)对应所述k空间数据k1,所述第二组混合k空间数据k13的偶数相位编码线数据(即虚线数据)对应所述k空间数据k3。
上述混合k空间相位编码数据线的数目为正常k空间相位编码数据线的两倍,相位编码方向k空间的步长是正常k空间相位编码线步长的一半,若用Δk(pe)表示正常k空间相位编码方向k空间的编码步长,则在混合k空间中其相位编码线步长为Δk(pe)/2,所以经过傅里叶变换后所得图像在相位编码方向的视野为正常k空间所得图像的两倍。这里,正常k空间指的是常规的k空间,即磁共振成像中利用传统方式进行数据采集得到的k空间,其概念属于本领域技术人员已知的范畴。Δk(pe)是从视野计算得到:
Δk(pe)=2π/(γFOVPE) [8]
式中,γ为旋磁比常数。对氢质子而言,γ=2π×42.57兆赫/(特斯拉);FOVPE为相位编码方向的视野。关于混合k空间技术,详情可参阅文献4:[Zhang W,“AQuantitative Analysis of Alternated Line Scanning in k Space and ItsApplication in MRI of Regional Tissue Perfusion by Arterial Spin Labeling”,Journal of Magnetic Resonance Series B 107:165-171(1995)]。
在具体实施过程中还有其他改变采集顺序得到混合k空间数据的采集方式,但只需保证最终得到两组混合k空间数据k12和k13,混合k空间数据的特征为混合k空间相位编码数据线的数目为正常k空间相位编码数据线数目的两倍,相位编码方向k空间的步长是正常k空间相位编码线步长的一半。混合k空间数据的数据形式为奇数相位编码线数据和偶数相位编码线数据各对应一种不同的流速编码,两组混合k空间数据中每组混合k空间数据所包含的两组k空间数据不完全相同即可,如第一组混合k空间数据对应k空间数据k1和k2,则第二组混合k空间数据可以对应k空间数据k1和k3,或k2和k3;第一组混合k空间数据对应k空间数据k1和k3,则第二组混合k空间数据可以对应k空间数据k1和k2,或k2和k3;第一组混合k空间数据对应k空间数据k2和k3,则第二组混合k空间数据可以对应k空间数据k1和k2,或k1和k3。
在数据采集方式上,本发明的技术方案可分为全采集方式和并行采集方式。以上无论是哪种采集顺序,均可由全采集方式或并行采集方式进而获取两组混合k空间数据。全采集方式在上述两种优选采集顺序方式下按行列逐点采集即得到完整的k空间数据k1、k2和k3或混合k空间数据k12和k13。
为了加快信号采集,最新的磁共振成像采用并行成像数据采集和图像重建,并行成像技术可以是图像域技术,如SENSE,参阅文献5:[Pruessmann KP,Weiger M,Scheidegger MB,Boesiger P,“SENSE:Sensitivity encoding for fast MRI”,Magn Reson Med42:952-962(1999)],或k空间技术,如GRAPPA,参阅文献6:[Griswold MA,Jakob PM,Heidemann RM,Nittka M,Jellus V,Wang J,KieferB,Haase A,Generalized autocalibrating partially parallel acquisitions(GRAPPA).Magn Reson Med47:1202–1210(2002)]。
对于并行采集方式下本发明技术方案的实施,当成像序列用并行成像数据采集方式采集进而获取三组k空间数据k1、k2和k3或两组混合k空间数据k12和k13时,由于其为并行采集,其由成像序列采集得到的各组k空间数据并非完整的k空间数据,之后需要经过并行加速成像数据重建得到完整的各组k空间数据。对于在交替采集方式下得到混合k空间数据的情况,则是将混合k空间数据中两组不同流速编码对应的k空间数据分离出来,分别进行并行加速数据重建。
并行加速成像数据重建的方法分为k空间技术和图像域技术,k空间技术(例如GRAPPA技术)是在k空间内进行欠采集数据的计算和填补,得到完整的k空间数据。对于分别采集三组k空间数据之后混合得到混合k空间数据的情况,则是首先由成像序列并行采集三组k空间数据k1’、k2’和k3’,k空间数据k1’、k2’和k3’为不完整的k空间数据,其分别对应流速编码m1(1)、m1(2)和m1(3)下采集得到,之后对k1’、k2’和k3’分别进行k空间数据重建,从而得到三组完整的k空间数据k1、k2和k3,最后混合得到混合k空间数据k12和k13。对于交替采集两组混合k空间数据的情况,则是首先由成像序列在流速编码m1(1)、m1(2)以及m1(1)、m1(3)下并行采集得到两组混合k空间数据k12’和k13’,之后将第一组不完整的混合k空间数据k12’根据不同流速编码所对应的数据分离为k1’和k2’,然后分别对k1’和k2’进行k空间并行加速数据重建得到完整的两组k空间数据k1和k2,最后将两组k空间数据k1和k2放回混合k空间中得到完整的混合k空间数据k12。第二组混合k空间数据的k空间并行加速重建方式相同。
图像域技术(例如SENSE技术)则是在图像域进行并行加速重建得到包含完整数据的图像,之后变换回k空间。对于分别采集三组k空间数据之后混合得到混合k空间数据的情况,则是首先由成像序列并行采集得到三组不完整的k空间数据k1’、k2’和k3’,然后分别将这三组k空间数据变换到图像域,在图像域进行并行加速重建,得到包含完整k空间数据的三组图像,之后将三组图像变换回k空间得到三组完整的k空间数据k1、k2和k3,最后混合得到混合k空间数据k12和k13。对于交替采集两组混合k空间数据的情况,则是首先由成像序列采集得到两组不完整的混合k空间数据k12’和k13’,将第一组不完整的混合k空间数据k12’根据不同流速编码所对应的k空间数据分离为k1’和k2’,然后将k1’和k2’分别变换到图像域,进行图像域的并行加速数据重建,并将经过并行加速数据重建后得到包含完整k空间数据的图像数据变换回k空间数据,得到两组完整的k空间数据k1和k2,将这两组数据放回混合k空间中得到完整的混合k空间数据k12,另一组混合k空间数据k13的图像域并行加速重建过程相同。
步骤S103:将所述第一组混合k空间数据k12和所述第二组混合k空间数据k13分别变换到图像域得到两幅图像,每幅图像均包含AC图和DC图。
第一组混合k空间图像信号:
第二组混合k空间图像信号:
式中j代表相位编码线序号(j=1,2,…N),其他符号与在先的定义相同。
由以上二组混合k空间数据k12、k13,通过傅里叶方法(或者其他数学变换方法,将混合k空间数据k12和k13变换为图像域数据)变换到图像域,可以得到两个混合k空间数据k12、k13所对应的二幅图像。所得每幅图像可以分为AC图和DC图,AC图分布在所得图像顶端到上四分之一处的位置以及图像下四分之一处的位置到图像底端,DC图分布在图像上四分之一到图像下四分之一的位置。之后,将经过傅里叶变换后所得图像沿相位编码方向滚动下移视野的四分之一,每幅图像含一个来自于衡值信号的DC图和一个来自于交变信号的AC图,如图3所示,每幅图像的上半幅图像为AC图,圆形示意为血管;下半幅图像为DC图,圆圈示意为血管周围的脂肪。DC图源于k空间线与线之间不变的静态信号,AC图源于k空间线与线之间交变的动态信号。
步骤S104:由所述两幅图像中的AC图利用相位反差计算流速。因为脂肪产生的是稳态信号,而血流信号为动态,所以从AC图使用相位差来测量流速流量不会受脂肪信号的干扰。
图像组一:
图像组二:
从IAC(1)和IAC(2)可测得流速ν,且不受脂肪信号的干扰:
ν=arg{IAC(2)×conj[IAC(1)]}/[m1(3)-m1(2)] [15]
式中arg{…}代表求复数相位;conj[…]代表复数共轭,IAC(1)代表第一幅混合k空间对应的AC图,IAC(2)代表第二幅混合k空间对应的AC图。
最后,根据需要,如需测量血液流量时,可在预设的流量测定范围内对所述流速积分得到流量。
本发明还提供一种相位反差磁共振流速测量装置,其中包括:
对成像序列分别进行三种不同的流速编码的序列产生装置;
由所述成像序列获取两组混合k空间数据k12和k13的数据处理装置,第一组混合k空间数据k12的奇数相位编码线数据对应k空间数据k1,第一组混合k空间数据k12的偶数相位编码线数据对应k空间数据k2,第二组混合k空间数据k13的奇数相位编码线数据对应k空间数据k1,第二组混合k空间数据k13的偶数相位编码线数据对应k空间数据k3,所述混合k空间数据k12和k13的相位编码线步长为正常k空间数据相位编码线步长的一半,所述k空间数据k1、k2和k3分别对应一种序列产生装置中所述流速编码;
将所述第一组混合k空间数据k12和所述第二组混合k空间数据k13分别变换到图像域形成两幅图像的第一图像变换装置,每幅图像均包含AC图和DC图;以及
根据所述两幅图像中的AC图利用相位反差计算流速的计算装置。
所述数据处理装置包括采集装置和数据处理子装置,所述采集装置为由所述成像序列采集得到k空间数据k1、k2和k3的第一采集装置,所述数据处理子装置适于将所述k空间数据k1和所述k空间数据k2混合得到第一组混合k空间数据k12,将所述k空间数据k1和所述k空间数据k3混合得到第二组混合k空间数据k13。或者所述数据处理装置包括采集装置,所述采集装置为由所述成像序列交替采集得到第一组混合k空间数据k12,由所述成像序列交替采集得到第二组混合k空间数据k13的第二采集装置。
所述采集装置包括并行采集装置和k空间数据重建装置,所述并行采集装置适于进行并行采集,所述k空间数据重建装置适于在采集完毕后进行k空间并行加速数据重建得到三组k空间数据k1、k2和k3。或者所述采集装置包括并行采集装置、图像域数据重建装置和第二图像变换装置,所述并行采集装置适于进行并行采集,所述图像域数据重建装置适于在采集完毕后进行图像域并行加速数据重建,所述第二图像变换装置适于将重建后图像变换回k空间。或者所述采集装置为全采集装置,其数据采集方式为采集完整数据。
所述采集装置、数据重建装置、图像变换装置以及计算装置可以选用现有的嵌入式集成芯片MCU或专用图像处理芯片实现。
本发明还提供一种相位反差磁共振流量测量装置,包括上述的相位反差磁共振流速测量装置;以及根据预设的流量测定范围对所述流速进行积分的计算装置。同理,所述计算装置可以采用现有的各种集成处理芯片,在此不再一一赘述。
虽然本发明已以较佳实施例揭示如上,然其并非用以限定本发明,任何本领域技术人员,在不脱离本发明的精神和范围内,当可作些许的修改和完善,因此本发明的保护范围当以权利要求书所界定的为准。