发明内容
本发明要解决的问题是现有磁共振成像方法中成像速度较慢且计算量较大。
为解决上述问题,本发明技术方案提供一种磁共振成像方法,包括:
通过第一变换操作将欠采集的K空间域变换至中间域,所述中间域位于K空间域与图像域之间;
分离出校准数据,并进行线圈合并系数的计算;
在所述中间域中,基于第二应采集数据以及所述线圈合并系数进行数据线性加权加和操作,以实现欠采集数据的填补;所述第二应采集数据由对在K空间域中采集的第一应采集数据进行所述第一变换操作后得到;
通过第二变换操作将完成欠采集数据填补之后的中间域变换至图像域以获得图像。
可选的,所述分离出校准数据,并进行线圈合并系数的计算是在所述中间域进行,所述校准数据为对在K空间域中采集的第一校准数据进行所述第一变换操作后得到的第二校准数据。
可选的,所述分离出校准数据,并进行线圈合并系数的计算是在K空间域进行,所述校准数据为在K空间域中采集的第一校准数据。
可选的,所述数据线性加权加和操作中引入或未引入频率编码方向的数据。
可选的,所述第一变换操作为频率编码方向的一维傅里叶变换,所述第二变换操作为相位编码方向的一维傅里叶变换。
为解决上述问题,本发明技术方案还提供一种磁共振成像装置,包括:
第一变换单元,适于通过第一变换操作将欠采集的K空间域变换至中间域,所述中间域位于K空间域与图像域之间;
计算单元,适于分离出校准数据,并进行线圈合并系数的计算;
填补单元,适于在所述中间域中,基于第二应采集数据以及所述线圈合并系数进行数据线性加权加和操作,以实现欠采集数据的填补;所述第二应采集数据由对在K空间域中采集的第一应采集数据进行所述第一变换操作后得到;
第二变换单元,适于通过第二变换操作将完成欠采集数据填补之后的中间域变换到图像域以获得图像。
与现有技术相比,本发明技术方案至少具有以下优点:
先通过第一变换操作将欠采集的K空间域变换至中间域,在所述中间域中完成填补欠采集数据的操作,再通过第二变换操作将完成欠采集数据填补的中间域变换至图像域以实现图像的并行采集重建,由于第一变换操作中无需对填补的欠采集数据进行变换,因此减少了第一变换操作在计算上的开销,总体上加快了计算速度,从而提高了磁共振成像的速度。
具体实施方式
为使本发明的上述目的、特征和优点能够更为明显易懂,下面结合附图对本发明的具体实施方式做详细的说明。在以下描述中阐述了具体细节以便于充分理解本发明。但是本发明能够以多种不同于在此描述的其它方式来实施,本领域技术人员可以在不违背本发明内涵的情况下做类似推广。因此本发明不受下面公开的具体实施方式的限制。
如背景技术所述,现有的磁共振成像方法,通过某种变换操作将并行采集重建后的K空间域变换至图像域的过程中,由于这种变换操作会消耗很多的时间,导致成像速度较慢,而且计算量较大。发明人经研究后分析如下:
现有的并行采集重建K空间的过程中,计算线圈合并系数和填补欠采集数据的操作全部是在K空间域中进行的,最终要得到图像,还需要对填补好的K空间域进行某种变换操作(如傅里叶变换)使其转换到图像域,该变换操作的对象包含采集的数据(即应采集数据和校准数据)以及计算填补的数据(即填补的欠采集数据)。由于该变换操作过程中需要变换的数据量较大,会消耗很多的时间,所以如何尽量减少空间变换操作,是提高图像产生速度的一个很重要的手段。
从K空间域转换到图像域的过程,一般可以分解为如图3所示的两个子过程:(a)通过变换操作1将K空间域转换至某个中间域;(b)通过变换操作2将该中间域转换至图像域。在实际实施时,一般都会将变换操作1和变换操作2合为一个变换操作过程,若该变换操作过程以二维傅里叶变换为例,则图3所示的中间域,可以认为是对K空间域进行了一维傅里叶变换(对应变换操作1)后的数据空间,将得到的数据空间再次进行一维傅里叶变换(对应变换操作2)便转换至图像域。
因此,发明人考虑,如果在K空间域不进行对欠采集数据进行填补的操作,而是通过变换操作1将欠采集的K空间域先变换至某个非K空间域的数据空间,该数据空间即上述的中间域,然后在所述中间域中实现对欠采集数据进行填补的操作,再通过变换操作2将完成填补操作之后的中间域变换至图像域,这样便可以从总体上减少变换操作过程中的计算量,提高磁共振成像的速度,理由如下:
由于一般的并行重建方法(如GRAPPA)是在K空间域对欠采集数据进行填补,填补的欠采集数据、应采集数据以及校准数据都需要进行某种变换操作(如二维傅里叶变换),才能转换为图像域的数据;而将填补欠采集数据的过程放在非K空间域(即本发明实施方式中所述的中间域)进行,则填补的欠采集数据是所述中间域中的数据,该数据不需要进行从K空间域到中间域的变换操作1,这样将节省浪费在变换操作1上的计算开销,由此能加快计算速度,进而提高磁共振成像的速度。
基于上述分析,本发明实施方式提供一种磁共振成像方法,在非K空间域的某一个数据空间内,进行数据线性加权加和操作以完成对欠采集数据的填补,从而实现磁共振图像的并行采集重建。该方法能够在一方面加速成像采集速度的同时,另一方面还可以提高重建图像的速度。
图4是本发明实施方式提供的磁共振成像方法的流程示意图。如图4所示,所述磁共振成像方法包括:
步骤S101,通过第一变换操作将欠采集的K空间域变换至中间域,所述中间域位于K空间域与图像域之间;
步骤S102,分离出校准数据,并进行线圈合并系数的计算;
步骤S103,在所述中间域中,基于第二应采集数据以及所述线圈合并系数进行数据线性加权加和操作,以实现欠采集数据的填补;所述第二应采集数据由对在K空间域中采集的第一应采集数据进行所述第一变换操作后得到;
步骤S104,通过第二变换操作将完成欠采集数据填补之后的中间域变换至图像域以获得图像。
下面以具体实施例对上述磁共振成像方法作详细说明。
在具体实施时,需要先采集K空间域的数据,所采集的数据包括应采集数据和校准数据,为了区别在中间域中对应的应采集数据和校准数据,因此将在K空间域中采集的应采集数据和校准数据分别称为第一应采集数据和第一校准数据。
在K空间域采集到所述第一应采集数据和第一校准数据之后,执行步骤S101,通过第一变换操作将欠采集的K空间域变换至中间域,所述中间域位于K空间域与图像域之间。
本实施例中,步骤S101中所述的“欠采集的K空间域”是指采集到所述第一应采集数据和第一校准数据之后形成的K空间域,该K空间域既区别于“未采集的K空间域”,也区别于“完成填补的K空间域”或者“并行采集重建后的K空间域”。因此,所述欠采集的K空间域中的数据包括在K空间域中采集的第一应采集数据和第一校准数据。
本实施例中,所述第一变换操作具体为频率编码方向的一维傅里叶变换,所述中间域具体为对所述欠采集的K空间域进行一维傅里叶变换之后得到的数据空间。欠采集的K空间域变换至中间域的过程,具体是对所述欠采集的K空间域中的数据进行所述第一变换操作的过程。本发明实施方式中,将对所述第一应采集数据进行第一变换操作之后得到的数据称为第二应采集数据,将对所述第一校准数据进行第一变换操作之后得到的数据称为第二校准数据。因此,所述中间域中的数据包括所述第二应采集数据和第二校准数据。
执行步骤S102,分离出校准数据,并进行线圈合并系数的计算。
本实施例中,步骤S102中所述的分离出校准数据,并进行线圈合并系数的计算是在所述中间域进行的,此时,步骤S102在步骤S101之后执行,步骤S102中所述的校准数据具体为所述第二校准数据。
在其他实施例中,所述分离出校准数据,并进行线圈合并系数的计算也可以在K空间域进行,此时,步骤S102在步骤S101之前执行,步骤S102中所述的校准数据具体为所述第一校准数据。
若所述分离出校准数据以及计算线圈合并系数的过程在K空间域进行,则可以采用现有的实现方式进行,此处不再详细描述;若所述分离出校准数据以及计算线圈合并系数的过程在所述中间域进行,其具体实现仍然可以参照在K空间域的实现方式进行,因为发明人发现,虽然经过所述第一变换操作,应采集数据和校准数据的数据内容发生了变化,但第二应采集数据和第二校准数据之间的对应关系依然保持不变,据此拟合出的线圈合并系数也与在K空间域中所计算出的线圈合并系数是相同的,这样通过步骤S102计算出线圈合并系数后,在后续步骤(步骤S103)便能够基于该线圈合并系数以及第二应采集数据在中间域进行数据线性加权加和操作,实现欠采集数据的填补,即在中间域对欠采集数据进行填补是可以实施的。
本实施例中,无论所述数据线性加权加和操作中是否引入频率编码方向的数据(即表明所述数据线性加权加和操作是否受到频率编码方向的影响),对于在中间域计算出的线圈合并系数并无影响,均与在K空间域中计算出的线圈合并系数相同,下面将对此进行详细说明。
1)不引入频率编码方向对数据线性加权加和操作的影响,有:
设计算的基矩阵(对应于图1中的黑色数据点)为Kb,校准线矩阵(对应于图1中的灰色数据点)为L,线圈合并系数为N,则有关系:
Kb·N=L (1)
则经过频率编码方向的一维傅里叶变换后,进入另一线性空间(中间域)其基矩阵变换为Fb,即:
Fb=F·Kb (2)
校准线矩阵变换为Fl,即:
Fl=F·L (3)
其中,F为傅里叶变换矩阵。
在式子(1)两边左乘F,可得:
Fb·N=Fl (4)
所以在中间域里线圈合并系数N不变,如果式子(4)对应的方程为超定的,则可以求得:
N=(FbH·Fb)-1FbH·Fl (5)
利用式子(5)得到的线圈合并系数N,可以在中间域对欠采集数据进行填补,即填补如图1所示的白色空缺点。
2)引入频率编码方向对数据线性加权加和操作的影响,有:
图5是频率编码方向影响下K空间域的数据分布示意图。引入频率编码方向对数据线性加权加和操作的影响,以图5为例说明。
引入频率编码方向的影响,在如图5所示的K空间域中,将会额外多出两个基矩阵,为Kb1和Kb2,分别为Kb向上向下平移一个ΔRo,ΔRo表示频率编码方向的步长。若记:
则有:
Kb1=T1·Kb (6)
Kb2=T2·Kb (7)
则在K空间域存在关系:
Kb1·N1+Kb2·N2+Kb·N3=L (8)
可以将式子(8)转为:
因为:
再将式子(6)、(7)代入式子(9),可以得到:
(T1 T2 E)·Kb·N=L (10)
经过频率编码方向的一维傅里叶变换之后进入中间域,(2)、(3)式仍然成立,对于(10)式两边左乘F,有:
(F·T1 F·T2 F)·Kb·N=F·L (11)
由于F为满秩矩阵,所以存在F-1,则记:
M1=F·T1·F-1(12)
M2=F·T2·F-1(13)
可知:
F·T1=M1·F (14)
F·T2=M2·F (15)
将式子(14)、(15)代入式子(11),可得:
(M1·F M2·F F)·Kb·N=F·L (16)
可以将式子(16)转为:
(M1 M2 E)·F·Kb·N=F·L (17)
将式子(2)、(3)代入式子(17),可得:
(M1 M2 E)·Fb·N=Fl (18)
式子(18)中,M1,M2,E,Fb,Fl都为已知量,由此可知,线圈合并系数N相对于K空间域不变。
记:
Fβ=(M1 M2 E)·Fb (19)
则有:
N=(FβH·Fβ)-1·FβH·Fl (20)
利用式子(20)得到的线圈合并系数N,可以在中间域对欠采集数据进行填补,即填补如图5所示的白色空缺点。
由于傅里叶变换以及T1、T2矩阵的一些特性,M1和M2都为对角阵,
若记:
则有:
在实际操作中,M1·Fb的操作,可以简化为Fb第一行至最后一行分别乘一个常数(对应于M1对角线上的数值,如第j行则乘以Ej)。
将K空间域转换至中间域以及计算出线圈合并系数之后,执行步骤S103,在所述中间域中,基于第二应采集数据以及所述线圈合并系数进行数据线性加权加和操作,以实现欠采集数据的填补;所述第二应采集数据由对在K空间域中采集的第一应采集数据进行所述第一变换操作后得到。
如前所述,无论是否受到频率编码方向的影响,均不影响步骤S103中所述的数据线性加权加和操作的结果,所述数据线性加权加和操作的具体实施可以参照公式(4)及公式(18)以及现有的在K空间域中GRAPPA算法填补操作的实现方式,此处不再详细描述。
需要说明的是,现有技术中,对于欠采集数据的填补在K空间域进行,线圈合并系数N可以通过下式得到:
N=(KbH·Kb)-1KbH·Fl (21)
由此可知,与在中间域计算出的线圈合并系数是相同的。
此外,通过步骤S103对欠采集数据进行填补后数据,与现有技术在K空间域对欠采集数据进行填补后数据有所区别,前者为中间域中的数据,而后者为K空间域中的数据,对前者进行所述第一变换操作可以得到后者。
完成数据线性加权加和操作,以实现对欠采集数据的填补之后,执行步骤S104,通过第二变换操作将完成欠采集数据填补之后的中间域变换至图像域以获得图像。
本实施例中,步骤S104中所述的第二变换操作具体为相位编码方向的一维傅里叶变换。完成欠采集数据填补之后的中间域中的数据,再经过一维傅里叶变换,即可以变换至图像域,实现磁共振成像。
对应于上述磁共振成像方法,本发明实施方式还提供一种磁共振成像装置。图6是本发明实施方式提供的磁共振成像装置的结构示意图。如图6所示,所述磁共振成像装置包括:第一变换单元101,适于通过第一变换操作将欠采集的K空间域变换至中间域,所述中间域位于K空间域与图像域之间;计算单元102,适于分离出校准数据,并进行线圈合并系数的计算;填补单元103,适于在所述中间域中,基于第二应采集数据以及所述线圈合并系数进行数据线性加权加和操作,以实现欠采集数据的填补;所述第二应采集数据由对在K空间域中采集的第一应采集数据进行所述第一变换操作后得到;第二变换单元104,适于通过第二变换操作将完成欠采集数据填补之后的中间域变换到图像域以获得图像。