本发明涉及成像技术领域,特别涉及一种基于平面回波成像的磁共振水脂分离和定量方法及装置。
背景技术:
相关技术中,磁共振成像相对高亮的脂肪信号会影响组织结构的观察以及一些重要指标的计算,例如,adc(apparentdiffusioncoefficient,表观扩散系数)等;同时,脂肪含量本身也可以作为临床诊断指标,因此脂肪信号的分离或定量在临床中具有很重要的应用。为了消除脂肪信号的影响,压脂技术得到了广泛应用,临床中常用的压脂技术包括spir(spectralpre-saturationbyinversionrecovery,选择性部分反转恢复法),spair(spectralattenuatedinversionrecovery,频谱衰减反转恢复),stir(shorttiinversionrecovery,短ti反转恢复),这些方法在某些情况下,尤其是在场不均匀性十分严重的区域,例如头颈区域,前额叶等,可能无法彻底的压制脂肪信号,特别地stir还会造成信噪比低的问题。
技术实现要素:
本申请是基于发明人对以下问题的认识和发现作出的:
epi(echoplanarimaging,平面回波采集技术)的快速采集特点使其具有成像速度快,对运动不敏感等优势,在临床中得到了广泛应用,尤其是单次激发epi,在一次rf(radiofrequency,射频)激发后完成整个k空间的采集,在对成像速度要求高的应用中具有非常重要的价值,例如扩散成像、灌注成像、心脏成像以及实时成像等。然而,epi采集也有它本身的不足,较长的读出时间会引入
单次激发epi与并行采集技术的结合可以减少读出窗的长度以及esp(effectiveechospacing,有效回波间隔),减少
与一般压脂技术不同,水脂分离或定量技术在消除脂肪信号对水相的影响的同时,还可以用于脂肪定量分析,常用如dixon方法,该方法利用水信号与脂肪信号共振频率的不同所导致的相位差异,采集不同回波时刻的图像,进而计算得出水相以及脂肪相。在此基础上发展出了多种水脂分离或定量方法,例如多点dixon方法,迭代重建的ideal(iterativedecompositionofwaterandfatwithechoasymmetryandleast-squareestimation)方法,hernando-gc(hernnado-graphcut)方法,qpbo(quadraticpseudobooleanoptimization)方法等。通过多次回波采集(如ffe(fastfieldecho,多回波的快速梯度回波))或者多次偏移读出窗口采集(如fse(fastspinecho,多采集的快速自旋回波)),计算得出水脂信号成分。
本发明旨在至少在一定程度上解决相关技术中的技术问题之一。
为此,本发明的一个目的在于提出一种基于平面回波成像的磁共振水脂分离和定量方法,该方法可以有效提高分离和定量的适用性,采集效率高,无变形、无t2*模糊效应、高信噪比,简单易实现。
本发明的另一个目的在于提出一种基于平面回波成像的磁共振水脂分离和定量装置。
为达到上述目的,本发明一方面实施例提出了一种基于平面回波成像的磁共振水脂分离和定量方法,包括以下步骤:使用无变形的点扩散函数编码的平面回波采集成像序列psf-epi进行数据采集,以得到多通道降采的k空间数据,根据所述k空间数据通过tilted-caipi恢复得到完整的三维k空间数据;使用所述三维k空间数据进行三维逆傅里叶变换,并通过通道合并得到三维幅值图像;从所述三维k空间数据抽取多个二维k空间数据,通过二维逆傅里叶变换得到复数图像并计算相应回波偏移时间;对所述复数图像和回波偏移时间进行水脂分离和定量计算,以得到每个像素中水信号和每个频率峰值处的脂肪信号比例;根据所述水信号及每个频率峰值处的脂肪信号比例对所述三维幅值图像在相位编码方向进行加权平方和计算,以得到无变形的水相以及各脂肪相。
本发明实施例的基于平面回波成像的磁共振水脂分离和定量方法,通过psf编码的epi采集无变形解剖结构图像或扩散图像,并将水脂分离或定量技术与其相结合,可以不需要额外采集数据或偏移数据读出窗口,便得到水脂分离后的图像或者脂肪百分比,直接同时计算得出水脂成分,从而有效提高分离和定量的适用性,采集效率高,无变形、无t2*模糊效应、高信噪比,简单易实现。
另外,根据本发明上述实施例的基于平面回波成像的磁共振水脂分离和定量方法还可以具有以下附加的技术特征:
进一步地,在本发明的一个实施例中,所述根据tilted-caipi恢复三维k空间数据,进一步包括:将相位编码方向和所述点扩散函数编码方向的所述多通道降采的k空间数据同时进行重建得到所述三维k空间数据,其中,重建过程为:
其中,di(a,b,c)是在(a,b,c)处第i个线圈需要被插值的数据,di′(a′,b′,c′)是在kernelk中采集到的数据,w为估计得到的插值权重系数,nc是线圈数目。
进一步地,在本发明的一个实施例中,所述三维幅值图像的三维幅值数据为:
i(s,y)=ρ(s)h(s,y),
其中,h(s,y)为点扩散函数,图像域内s,y分别对应k空间域的ks,ky方向,ρ(s)为无变形的图像信号强度,δ为冲激函数,δ(s)为图像变形所导致的点扩散函数偏移。
进一步地,在本发明的一个实施例中,所述回波偏移时间为:
其中,
进一步地,在本发明的一个实施例中,所述每个像素中水信号为:
其中,ρw,q以及ρf,q分别代表像素点q处的水信号以及脂肪信号成分,fb,q为该处的场强偏移量,
进一步地,在本发明的一个实施例中,所述水相以及各脂肪相为:
其中,*代表共轭操作,不同峰值处的脂肪相依次计算为
为达到上述目的,本发明另一方面实施例提出了一种基于平面回波成像的磁共振水脂分离和定量装置,包括:采集模块,用于采用无变形的点扩散函数编码的平面回波采集成像序列psf-epi进行数据采集,以得到多通道降采的k空间数据,并根据所述k空间数据通过tilted-caipi恢复得到三维k空间数据;变换合并模块,用于对所述k空间数据进行三维逆傅里叶变换,并通过通道合并得到三维幅值图像;抽取模块,用于从所述三维k空间数据抽取多个二维k空间数据,通过二维逆傅里叶变换得到复数图像并计算相应的回波偏移时间;分离与计算模块,用于使用所述复数图像和回波偏移时间进行水脂分离和定量计算,以得到每个像素中水信号和每个频率峰值处的脂肪信号比例;计算模块,用于根据所述水信号及每个频率峰值处的脂肪信号比例对所述三维幅值图像在相位编码方向进行加权平方和计算,以得到无变形的水相以及各脂肪相。
本发明实施例的基于平面回波成像的磁共振水脂分离和定量装置,通过psf编码的epi采集无变形解剖结构图像或扩散图像,并将水脂分离或定量技术与其相结合,可以不需要额外采集数据或偏移数据读出窗口,便得到水脂分离后的图像或者脂肪百分比,直接同时计算得出水脂成分,从而有效提高分离和定量的适用性,采集效率高,无变形、无t2*模糊效应、高信噪比,简单易实现。
另外,根据本发明上述实施例的基于平面回波成像的磁共振水脂分离和定量装置还可以具有以下附加的技术特征:
进一步地,在本发明的一个实施例中,所述采集模块进一步用于将相位编码方向和所述点扩散函数编码方向的所述多通道降采的k空间数据同时进行重建得到所述三维k空间数据,其中,重建过程为:
其中,di(a,b,c)是在(a,b,c)处第i个线圈需要被插值的数据,di′(a′,b′,c′)是在kernelk中采集到的数据,w为估计得到的插值权重系数,nc是线圈数目。
进一步地,在本发明的一个实施例中,所述三维幅值图像的三维幅值数据为:
i(s,y)=ρ(s)h(s,y),
其中,h(s,y)为点扩散函数,图像域内s,y分别对应k空间域的ks,ky方向,ρ(s)为无变形的图像信号强度,δ为冲激函数,δ(s)为图像变形所导致的点扩散函数偏移。
进一步地,在本发明的一个实施例中,其中,所述回波偏移时间为:
其中,
其中,ρw,q以及ρf,q分别代表像素点q处的水信号以及脂肪信号成分,fb,q为该处的场强偏移量,
其中,*代表共轭操作,不同峰值处的脂肪相依次计算为
本发明附加的方面和优点将在下面的描述中部分给出,部分将从下面的描述中变得明显,或通过本发明的实践了解到。
附图说明
本发明上述的和/或附加的方面和优点从下面结合附图对实施例的描述中将变得明显和容易理解,其中:
图1为根据本发明一个实施例的基于平面回波成像的磁共振水脂分离和定量方法的流程图;
图2为根据本发明一个具体实施例的基于平面回波成像的磁共振水脂分离和定量方法的流程图;
图3为根据本发明一个实施例的点扩散函数编码的平面回波采集磁共振成像扫描序列示意图;
图4为根据本发明一个实施例的幅值图像获取流程图;
图5为根据本发明一个实施例的二维空间数据抽取流程图;
图6为根据本发明一个实施例的基于平面回波成像的磁共振水脂分离和定量装置的结构示意图。
具体实施方式
下面详细描述本发明的实施例,所述实施例的示例在附图中示出,其中自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。下面通过参考附图描述的实施例是示例性的,旨在用于解释本发明,而不能理解为对本发明的限制。
下面参照附图描述根据本发明实施例提出的基于平面回波成像的磁共振水脂分离和定量方法及装置,首先将参照附图描述根据本发明实施例提出的基于平面回波成像的磁共振水脂分离和定量方法。
图1是本发明一个实施例的基于平面回波成像的磁共振水脂分离和定量方法的流程图。
如图1所示,该基于平面回波成像的磁共振水脂分离和定量方法包括以下步骤:
在步骤s101中,使用无变形的点扩散函数编码的平面回波采集成像序列psf-epi进行数据采集,以得到多通道降采的k空间数据,并根据k空间数据通过tilted-caipi恢复得到完整的三维k空间数据。
可以理解的是,如图2所示,本发明实施例使用基于psf编码的epi成像序列进行数据采集,获取多通道降采的k空间数据,并使用tilted-caipi技术对降采样的k空间信号进行恢复,获得完整的psf编码的三维k空间数据。
进一步地,在本发明的一个实施例中,根据tilted-caipi恢复三维k空间数据,进一步包括:将相位编码方向和点扩散函数编码方向的多通道降采的k空间数据同时进行重建得到三维k空间数据,其中,重建过程为:
其中,di(a,b,c)是在(a,b,c)处第i个线圈需要被插值的数据,di′(a′,b′,c′)是在kernelk中采集到的数据,w为估计得到的插值权重系数,nc是线圈数目。
具体而言,点扩散函数编码的磁共振成像扫描序列如图3所示(以自旋回波的扩散磁共振成像举例),在传统的单次激发epi序列的相位编码之前会施加一个额外的点扩散函数相位编码,其k空间编码间隔应与基础的epi相位编码保持一致,编码范围取决于需求的图像分辨率,从而采集得到的k空间数据将存在三个维度,分别为读出方向kx,epi相位编码方向(简称epi-pe)ky,点扩散函数相位编码方向(简称psf-pe)ks。针对扩散磁共振成像,为了进行重建过程中的相位矫正,在成像回波之后,可施加180°回聚脉冲以采集二维导航回波信号。psf-pe编码方向的并行成像加速倍数可以使用如10倍、14倍、18倍等,epi-pe编码方向的并行成像加速倍数可以使用如2倍、4倍等,同时在psf-pe及epi-pe方向均可以实施部分傅里叶采集,本发明实施例对成像加速倍数及部分傅里叶采集比例不加限制。
高倍降采的psf编码的epi成像数据恢复采用tilted-caipi技术。该方法将相位编码和psf编码方向的降采数据同时进行重建,如图4a所示,重建过程可以被表示为:
该式中di(a,b,c)是在(a,b,c)(对应kx,ky,ks)处第i个线圈需要被插值的数据;di′(a′,b′,c′)是在kernelk中采集到的数据,w为估计得到的插值权重系数,nc是线圈数目。在此重建中使用的权重w包含了线圈敏感度信息,以及相位编码和psf编码方向的相关性信息(磁场不均匀性信息)。为了训练得到该权重,需要采集校准数据。校准数据需要在psf和相位编码方向满采样,或者恢复为满采样,并提供足够的数据来进行kernel计算,由于所采集的k空间数据在对角线方向具有最高的相关性,所以插值的kernel沿对角线方向,如图4a所示。插值完成后,可以使用如pocs方法进行三维部分傅里叶重建。tilted-caipi技术能够将psf-epi序列加速20倍以上,实现4-8次激发的高分辨率无几何失真成像。
本发明实施例所采用的基于epi的成像序列不局限于实施例中的自旋回波epi采集(se-epi),亦可以使用梯度回波epi采集(ge-epi),反转恢复回波epi采集(ir-epi),受激回波epi采集(steamepi)等;本发明对成像对比度不加限制,包括各种t1加权、t2加权、t2*加权、质子密度(pd)加权、扩散加权(dwi)等;本发明实施例也可使用的扩散成像的扩散准备序列不局限于实施例中的pgse(pulsedgradientspinecho,脉冲梯度自旋回波序列),亦可以使用ste(simulatedecho,受激回波扩散准备序列),ogse(oscillatinggradientspinecho,振荡梯度自旋回波扩散准备序列),dde(doublediffusionencoding,双重扩散编码扩散准备序列),code(convexoptimizeddiffusionencoding,凸优化扩散编码扩散准备序列)等;本发明实施例中扩散磁共振成像的相位矫正的方式不局限于实施例中的额外采集导航回波数据方法,亦可以采用自导航的重建方式(由降采数据自身重建得到低分辨率图像利用其相位信息用于相位矫正)等,在此不做具体限定。
在步骤s102中,对三维k空间数据进行三维逆傅里叶变换,并通过通道合并得到三维幅值图像。
可以理解的是,如图2所示,本发明实施例对重建得到的k空间数据进行三维逆傅里叶变换并使用通道合并获取包含psf信息的幅值图像。
进一步地,在本发明的一个实施例中,三维幅值图像的三维幅值数据为:
i(s,y)=ρ(s)h(s,y),
其中,h(s,y)为点扩散函数,图像域内s,y分别对应k空间域的ks,ky方向,ρ(s)为无变形的图像信号强度,δ为冲激函数,δ(s)为图像变形所导致的点扩散函数偏移。
具体而言,对降采恢复之后的k空间数据进行三维逆傅里叶变换并使用通道合并,由于读出方向独立于epi-pe及psf-pe方向,所得到的三维幅值数据可以简化表示为:
i(s,y)=ρ(s)h(s,y),[2]
其中:
即为点扩散函数。图像域内s,y分别对应k空间域的ks,ky方向,ρ(s)为无变形的图像信号强度,δ为冲激函数,δ(s)为图像变形所导致的点扩散函数偏移,如图4b所示。对该幅值图像按下式分别沿y方向和s方向进行积分可以得到无变形及有变形的图像,如图4c所示:
i(s)=∫i(s,y)dy,[4]
i(y)=∫i(s,y)ds,[5]
本发明实施例不限制所采用的通道合并方法,本领域技术人员可以根基实际情况选择相应的方法,例如,对多次激发的多通道的图像域数据进行合并的方法主要包括平方和sos方法、最优化信噪比方法、自适应重建acc方法、主成分分析pca方法和奇异值分解svd方法等,在此不做具体限定。本发明实施例可通过但并不限于上述任意一种方法将所有多次激发的多通道的图像域数据合并为一幅图像。
在步骤s103中,从三维k空间数据抽取多个二维k空间数据,并通过二维逆傅里叶变换得到复数图像并计算相应回波偏移时间。
可以理解的是,如图2所示,从恢复的三维k空间数据中使用全部或抽取若干用于水脂分离或定量的二维k空间数据(例如,5或7个),通过二维逆傅里叶变换获得相应的复数图像并计算对应的回波偏移时间。
进一步地,在本发明的一个实施例中,回波偏移时间为:
其中,
具体而言,考虑重建恢复的三维k空间数据,对于每个epi相位编码步骤nδky而言,其中n为编码步数,δky为编码间隔(如图5a所示,以最大编码步数n为6举例,n的取值范围为-6至6),kx-ks平面可以看成一个传统的二维k空间数据,以
其中,
从
其中,ρw(x)为位于x处的水信号成分,ρf(x)为位于x处的脂肪信号成分,fb(x)为该处的场不均匀性,ff为脂肪信号相对水信号的共振频率偏移,tn为不同的回波偏移时间,s(x;tn)为回波偏移为tn时位于x处的图像信号。
在步骤s104中,使用复数图像和回波偏移时间进行水脂分离和定量计算,以得到每个像素中水信号和每个频率峰值处的脂肪信号比例。
可以理解的是,如图2所示,本发明实施例使用步骤s103中获得的复数图像以及对应的回波偏移时间采用水脂分离或定量方法进行水脂分离以及定量计算,计算得到每个像素中水信号和各频率峰值处的脂肪信号的各自比例。
进一步地,在本发明的一个实施例中,每个像素中水信号为:
其中,ρw,q以及ρf,q分别代表像素点q处的水信号以及脂肪信号成分,fb,q为该处的场强偏移量,
具体而言,如上述介绍可以得到l个位于不同回波时刻的图像
其中,ρw,q以及ρf,q分别代表像素点q处的水信号以及脂肪信号成分,fb,q为该处的场强偏移量,
本发明实施例不限制所采用的水脂分离方法,常用的水脂分离及定量计算方法有:ideal(iterativedecompositionofwaterandfatwithechoasymmetryandleast-squareestimation)方法,hierarchicalideal方法,max-ideal方法,qpbo(quadraticpseudobooleanoptimization)方法,qpbo-gc(quadraticpseudobooleanoptimizationgraphcut)方法,hernando-gc(graphcut)方法,asr(asafest-firstregion-growing)方法,mrgs(multi-resolutiongoldensectionsearch)等。本发明实施例可通过但并不限于上述任意一种方法进行水脂分离以及定量计算。
在步骤s105中,根据水信号及每个频率峰值处的脂肪信号比例对三维幅值图像在相位编码方向进行加权平方和计算,以得到无变形的水相以及各脂肪相。
可以理解的是,如图2所示,本发明实施例利用计算得到的水脂信号比例对步骤s102中的三维幅值图像在相位编码方向进行加权平方和计算,从而得到无变形的水相以及各脂肪相。
进一步地,在本发明的一个实施例中,水相以及各脂肪相为:
其中,*代表共轭操作,不同峰值处的脂肪相依次计算为
具体而言,为了进一步提高图像的信噪比,按照步骤s104中计算得出水、脂信号成分后,可以计算出各像素点水信号与各峰值处脂肪信号的定量比重,分别为
其中,*代表共轭操作,不同峰值处的脂肪相可以依次计算为:
综上,本发明实施例将psf-epi采集技术与水脂分离或定量方法相结合,无需额外的回波采集或偏移读出窗口采集,利用所采集的psf-epi数据本身,提取可用于水脂分离或定量的回波信号,使用前述的常用技术进行水脂分离和定量。本发明实施例方法继承了epi的快速采集优势以及psf-epi的无变形、无t2*模糊效应、高信噪比等优点,可以在psf-epi的各种成像场合(如各种解剖结构像、扩散成像等)直接计算得出各水脂成分,用于脂肪定量分析或者去除脂肪信号对各种定量分析的影响,如扩散系数adc和各向异性fa等。
根据本发明实施例提出的基于平面回波成像的磁共振水脂分离和定量方法,通过psf编码的epi采集无变形解剖结构图像或扩散图像,并将水脂分离或定量技术与其相结合,可以不需要额外采集数据或偏移数据读出窗口,便得到水脂分离后的图像或者脂肪百分比,直接同时计算得出水脂成分,从而有效提高分离和定量的适用性,采集效率高,无变形、无t2*模糊效应、高信噪比,简单易实现。
其次参照附图描述根据本发明实施例提出的基于平面回波成像的磁共振水脂分离和定量装置。
图6是本发明一个实施例的基于平面回波成像的磁共振水脂分离和定量装置的结构示意图。
如图6所示,该基于平面回波成像的磁共振水脂分离和定量装置10包括:采集模块100、变换合并模块200、抽取模块300、分离与计算模块400和计算模块500。
其中,采集模块100用于采用无变形的点扩散函数编码的平面回波采集成像序列psf-epi进行数据采集,以得到多通道降采的k空间数据,并根据k空间数据通过tilted-caipi恢复得到三维k空间数据。变换合并模块200用于对三维k空间数据进行三维逆傅里叶变换,并通过通道合并得到三维幅值图像。抽取模块300用于从三维k空间数据抽取多个二维k空间数据,通过二维逆傅里叶变换得到复数图像并计算相应的回波偏移时间。分离与计算模块400用于使用复数图像和回波偏移时间进行水脂分离和定量计算,以得到每个像素中水信号和每个频率峰值处的脂肪信号比例。计算模块500用于根据水信号及每个频率峰值处的脂肪信号比例对三维幅值图像在相位编码方向进行加权平方和计算,以得到无变形的水相以及各脂肪相。本发明实施例的装置10通过psf编码的epi采集无变形解剖结构图像或扩散图像,并将水脂分离或定量技术与其相结合,从而有效提高分离和定量的适用性,采集效率高,无变形、无t2*模糊效应、高信噪比,简单易实现。
进一步地,在本发明的一个实施例中,采集模块100进一步用于将相位编码方向和点扩散函数编码方向的多通道降采的k空间数据同时进行重建得到三维k空间数据,其中,重建过程为:
其中,di(a,b,c)是在(a,b,c)处第i个线圈需要被插值的数据,di′(a′,b′,c′)是在kernelk中采集到的数据,w为估计得到的插值权重系数,nc是线圈数目。
进一步地,在本发明的一个实施例中,三维幅值图像的三维幅值数据为:
i(s,y)=ρ(s)h(s,y),
其中,h(s,y)为点扩散函数,图像域内s,y分别对应k空间域的ks,ky方向,ρ(s)为无变形的图像信号强度,δ为冲激函数,δ(s)为图像变形所导致的点扩散函数偏移。
进一步地,在本发明的一个实施例中,其中,回波偏移时间为:
其中,
其中,ρw,q以及ρf,q分别代表像素点q处的水信号以及脂肪信号成分,fb,q为该处的场强偏移量,
其中,*代表共轭操作,不同峰值处的脂肪相依次计算为
需要说明的是,前述对基于平面回波成像的磁共振水脂分离和定量方法实施例的解释说明也适用于该实施例的基于平面回波成像的磁共振水脂分离和定量装置,此处不再赘述。
根据本发明实施例提出的基于平面回波成像的磁共振水脂分离和定量装置,通过psf编码的epi采集无变形解剖结构图像或扩散图像,并将水脂分离或定量技术与其相结合,可以不需要额外采集数据或偏移数据读出窗口,便得到水脂分离后的图像或者脂肪百分比,直接同时计算得出水脂成分,从而有效提高分离和定量的适用性,采集效率高,无变形、无t2*模糊效应、高信噪比,简单易实现。
此外,术语“第一”、“第二”仅用于描述目的,而不能理解为指示或暗示相对重要性或者隐含指明所指示的技术特征的数量。由此,限定有“第一”、“第二”的特征可以明示或者隐含地包括至少一个该特征。在本发明的描述中,“多个”的含义是至少两个,例如两个,三个等,除非另有明确具体的限定。
在本说明书的描述中,参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本发明的至少一个实施例或示例中。在本说明书中,对上述术语的示意性表述不必须针对的是相同的实施例或示例。而且,描述的具体特征、结构、材料或者特点可以在任一个或多个实施例或示例中以合适的方式结合。此外,在不相互矛盾的情况下,本领域的技术人员可以将本说明书中描述的不同实施例或示例以及不同实施例或示例的特征进行结合和组合。
尽管上面已经示出和描述了本发明的实施例,可以理解的是,上述实施例是示例性的,不能理解为对本发明的限制,本领域的普通技术人员在本发明的范围内可以对上述实施例进行变化、修改、替换和变型。