k空间运动伪影矫正方法及装置
技术领域
本发明涉及k空间伪影矫正,特别涉及一种k空间运动伪影矫正方法及装置。
背景技术
磁共振成像设备是一种能够将人体器官组织中氢原子信号进行检测并形成可用于医疗诊断图像的影像设备。磁共振设备在检查区域内部提供了高强度的极化磁场、用于空间编码的梯度场,以及使人体器官组织中氢原子发生共振现象产生共振信号的射频场,其产生的共振信号被射频接收线圈接收,经过后续设备的信号处理最终转化为能够分辨出不同组织的诊断图像。其对于软组织有较好的成像效果,因此在脑部、肌肉、心脏以及癌细胞的诊断上相比于其他医学影像设备如CT、x光等应用更为广泛。
大部分磁共振成像所需时间较长,在对某些器官成像过程中,如当人呼吸时,内脏器官也会随着隔膜而移动,又如心脏会周期性运动,因此在需要进行此类器官组织的动态磁共振成像时,运动伪影较难消除。
现有技术中在进行与时间相关的磁共振成像时,传统的方法需要患者屏气来完成图像的采集。成像时间较长时,对患者带来很大的负担。目前也有一些其他成像方法不需患者屏气,减轻患者的痛苦,如导航技术、PROPELLER(Periodically Rotated OverlappingParallEL Lines with Enhanced Reconstruction,螺旋桨)技术。导航技术的原理为其在磁共振成像的同时,利用导航回波对成像器官组织的运动进行检测,将采集磁共振信号的时间选择确定为满足一定条件的时间内或者由导航回波信息对实际采集的磁共振数据相对于器官组织的运动进行必要的修正。PROPELLER技术的基本思想为,将k空间采样分成几个部分顺序进行,如图1所示,每个部分叫做一个k空间条,在采集每个k空间条时由于采集很快完成可以近似认为受检者是静止的,而运动仅仅发生在采集k空间条之间,运动的方向和幅度可以由k空间条之间重叠采样区域的数据估计得到,对运动加以补偿后,就可以重建出不受运动伪影干扰的优质图像。
然而导航技术需要关注成像以外的其他信息,如膈肌运动。而PROPELLER技术在k空间内重复采集,其每次采集均重复采集中心区域,重叠严重,影响了成像的时间,且只适用于单帧图像成像,没有在动态磁共振成像中运动伪影矫正的应用。
发明内容
本发明要解决的问题是提供一种k空间运动伪影矫正方法及装置,在磁共振和时间相关的成像方法中,矫正成像物体的刚体运动伪影,缓解非预期运动对成像造成的影响(如心脏扫描中的呼吸运动造成的影响),减轻成像所带来的对病人的影响(如心脏电影时所需屏气的痛苦),并且减少采集数据,缩短成像时间。
为解决上述问题,本发明提供了一种k空间运动伪影矫正方法,其特征在于,包括以下步骤:
S101,在成像区域连续采集多帧k空间数据,所述各帧内的k空间数据由至少一次射频激发采集得到,相邻帧数据的采集方向不同;
S102,将所述各帧中由一次射频激发采集得到的数据分为一组,所述各帧内至少有一个数据组;
S103,以前一帧中的一个数据组作为标准组,后一帧中的数据组作为待矫正组,利用所述标准组的数据对待矫正组的数据进行数据矫正,依次完成各帧数据组的数据矫正;
S104,将数据矫正后的各帧数据组填回相应帧的k空间中,得到各帧矫正完毕的k空间数据。
优选的,所述步骤S101中采集多帧k空间数据时,奇数帧数据采集方向为kx方向,偶数帧数据采集方向为ky方向。
优选的,所述步骤S103中,所述数据矫正包括:k空间中心矫正和/或旋转角度矫正和/或平移矫正。
优选的,所述k空间中心矫正包括以下步骤:
将待矫正组k空间数据DataZ1变换到图像域得到数据ImageZ1;
将待矫正组k空间数据DataZ1加窗后变换到图像域得到数据ImageZ1Window;
取所述数据ImageZ1的模值为矫正后图像数据Image的模值;
取所述数据ImageZ1和所述数据ImageZ1Window的相位差为所述矫正后图像数据Image的相位;
将所述图像数据Image变换到k空间得到矫正完成数据DataZ1Corr。
优选的,所述旋转角度矫正包括以下步骤:
将标准组数据和待矫正组数据分别取模值,所述标准组数据和待矫正组数据在ky和kx方向大小分别为K1和K2,所述标准组数据在ky方向步长为ΔK1,所述待矫正组数据在kx方向步长为ΔK2;
将取模值后的所述待矫正组数据在k空间内旋转α角度并插值到公共网格空间得到数据组MCi,所述公共网格空间在ky方向大小为K22,步长为ΔK22,在kx方向大小为K11,步长为ΔK11,所述MCi表示经过第i次旋转后插值到公共网格空间得到的数据组,i为正整数,α为0°到360°之间任一值;
将取模值后的所述标准组数据插值到所述公共网格空间得到数据组MB,计算MB与MCi的相关性;
找到最大相关性值对应的旋转角度α值,以所述旋转角度值为待矫正组数据作旋转角度矫正。
优选的,所述平移矫正包括以下步骤:
将标准组数据和待矫正组数据插值到同一公共网格空间,所述标准组数据和待矫正组数据在ky和kx方向大小分别为K1和K2,所述标准组数据在ky方向步长为ΔK1,所述待矫正组数据在kx方向步长为ΔK2,所述公共网格空间在ky方向大小为K22,步长为ΔK22,在kx方向大小为K11,步长为ΔK11;
将公共网格空间内所述标准组数据与所述待矫正组数据的共轭的乘积变换到图像域得到图像,拟合找到图像最大值位置;
所述图像最大值位置与所述图像的中心位置的偏移量作为平移运动造成的偏移量,以所述偏移量为待矫正组作平移矫正。
优选的,所述K11=K1,K22=K2,ΔK11=ΔK1,ΔK22=ΔK2。
本发明技术方案还提供了一种k空间运动伪影矫正装置,包括:
采集单元101,适于在成像区域连续采集多帧k空间数据,所述各帧内的k空间数据由至少一次射频激发采集得到,相邻帧数据的采集方向不同;
分组单元102,适于将所述各帧中由一次射频激发采集得到的数据分为一组,所述各帧内至少有一个数据组;
矫正单元103,适于将前一帧中的一个数据组作为标准组,后一帧中的数据组作为待矫正组,利用所述标准组的数据对待矫正组的数据进行数据矫正,依次完成各帧数据组的数据矫正;
填补单元104,适于将矫正后的各帧数据组填回相应帧的k空间中,得到各帧矫正完毕的k空间数据。
优选的,所述采集单元采集的多帧k空间数据,其中奇数帧数据采集方向为kx方向,偶数帧数据采集方向为ky方向。
优选的,所述矫正单元103还包括:k空间中心矫正单元1031,适于进行k空间中心矫正;旋转矫正单元1032,与所述k空间中心矫正单元1031相连,适于进行旋转矫正;平移矫正单元1033,与所述旋转矫正单元1032相连,适于进行平移矫正;其中一个或多个矫正单元可省略。
优选的,所述矫正单元103还包括:k空间中心矫正单元1031,适于进行k空间中心矫正;旋转矫正单元1032,与所述k空间中心矫正单元1031相连,适于进行旋转矫正;平移矫正单元1033,与所述旋转矫正单元1032相连,适于进行平移矫正;其中一个或多个矫正单元可省略。
优选的,所述k空间中心矫正单元1031适于执行以下步骤:
将待矫正组k空间数据DataZ1变换到图像域得到数据ImageZ1;
将待矫正组k空间数据DataZ1加窗后变换到图像域得到数据ImageZ1Window;
取所述数据ImageZ1的模值为矫正后图像数据Image的模值;
取所述数据ImageZ1和所述数据ImageZ1Window的相位差为所述矫正后图像数据Image的相位;
将所述图像数据Image变换到k空间得到矫正完成数据DataZ1Corr。
优选的,所述旋转角度矫正单元1032适于执行以下步骤:
将标准组数据和待矫正组数据分别取模值,所述标准组数据和待矫正组数据在ky和kx方向大小分别为K1和K2,所述标准组数据在ky方向步长为ΔK1,所述待矫正组数据在kx方向步长为ΔK2;
将取模值后的所述待矫正组数据在k空间内旋转α角度并插值到公共网格空间得到数据组MCi,所述公共网格空间在ky方向大小为K22,步长为ΔK22,在kx方向大小为K11,步长为ΔK11,所述MCi表示经过第i次旋转后插值到公共网格空间得到的数据组,i为正整数,α为0°到360°之间任一值;
将取模值后的所述标准组数据插值到所述公共网格空间得到数据组MB,计算MB与MCi的相关性;
找到最大相关性值对应的旋转角度α值,以所述旋转角度值为待矫正组数据作旋转角度矫正。
优选的,所述平移矫正单元1033适于执行以下步骤:
将标准组数据和待矫正组数据插值到同一公共网格空间,所述标准组数据和待矫正组数据在ky和kx方向大小分别为K1和K2,所述标准组数据在ky方向步长为ΔK1,所述待矫正组数据在kx方向步长为ΔK2,所述公共网格空间在ky方向大小为K22,步长为ΔK22,在kx方向大小为K11,步长为ΔK11;
将公共网格空间内所述标准组数据与所述待矫正组数据的共轭的乘积变换到图像域得到图像,拟合找到图像最大值位置;
所述图像最大值位置与所述图像的中心位置的偏移量作为平移运动造成的偏移量,以所述偏移量为待矫正组作平移矫正。
优选的,所述K11=K1,K22=K2,ΔK11=ΔK1,ΔK22=ΔK2。
与现有技术相比,本发明提供的k空间运动伪影矫正方法及装置,利用邻近帧数据的近似性,对采集的数据进行运动伪影校正。利用本方法可以缓解非预期运动对成像造成的影响(如心脏扫描中的呼吸运动造成的影响),减轻成像所带来的对病人的影响(如心脏电影时所需屏气的痛苦),采集数据较少,加快成像速度。
附图说明
图1为现有技术PROPELLER采集示意图;
图2为本发明k空间运动伪影矫正流程图;
图3为本发明实施例一多帧k空间数据采集示意图;
图4为本发明实施例一数据组在公共空间插值示意图;
图5为本发明k空间运动伪影矫正装置结构示意图;
图6为本发明实施例二多帧k空间数据采集示意图。
具体实施方式
为使本发明的上述目的、特征和优点能够更为明显易懂,下面结合附图对本发明的具体实施方式做详细的说明。
如背景技术所述,现有的k空间运动伪影矫正方案中,屏气对病人带来极大负担,导航技术需要关注成像以外的信息,如隔膜运动、心脏周期性跳动,PROPELLER技术在采集过程中数据重叠严重,且只关注于单帧图像矫正,没有多帧磁共振成像中的应用。本发明实施过程中,能够进行多帧图像运动伪影矫正,在数据采集时,相邻帧数据采集方向不同,无大量重叠数据,借用邻近帧数据近似性,在保证运动伪影矫正数据准确性的情况下减少采集多余重复数据的时间损耗,提高成像速度。
图2为本发明k空间运动伪影矫正方法的流程示意图,包括以下步骤:
步骤S101,在成像区域连续采集多帧k空间数据,所述各帧内的k空间数据由至少一次射频激发采集得到,相邻帧数据的采集方向不同;
步骤S102,将所述各帧中由一次射频激发采集得到的数据分为一组,所述各帧内至少有一个数据组;
步骤S103,以前一帧中的一个数据组作为标准组,后一帧中的数据组作为待矫正组,利用所述标准组的数据对待矫正组的数据进行数据矫正,依次完成各帧数据组的数据矫正;
步骤S104,将数据矫正后的各帧数据组填回相应帧的k空间中,得到各帧矫正完毕的k空间数据。
下面以具体实施例对上述k空间运动伪影矫正方法作详细说明。
实施例一
首先执行步骤S101,在成像区域连续采集多帧k空间数据,所述各帧内的k空间数据由至少一次射频激发采集得到,相邻帧数据的采集方向不同。这里需要说明的是,与现有技术PROPELLER采集方式相比,PROPELLER技术不是采集多帧k空间数据,其没有帧的概念,其在采集时是对一帧图像采集多组k空间数据,各组数据为连续步进式地在多个采集方向进行旋转采集得到。本发明中的采集方式为采集多帧k空间数据,相邻帧k空间数据采集方向不同即可,各帧采集方向间不需要固定的旋转角度间隔,且这样采集的好处在于,没有过多的重复采集数据。
实施例一中本发明的一种优选实施方式,其k空间采集方式为奇数帧数据采集方向均为kx方向,偶数帧数据采集方向均为ky方向,其k空间数据采集方式如图3所示。请参见图3,这里只示意出了前四帧k空间数据分布情况,之后省略部分与前面部分类似,其中RO为读出方向(也称为采集方向),PE为相位编码方向,kx和ky为k空间内的空间维度,t表示时间维度,其符号标注均为本领域技术人员通用标准。此外,每帧内的粗线和细线为该帧内分别经过一次射频激发采集得到的k空间数据。一次射频激发的数据采集间隔时间短,可以认为一次激发后采集到的数据间没有发生刚体运动。
以上选取奇数帧数据采集方向为kx方向,偶数帧数据采集方向为ky方向为方便计算所选取的,实际实施过程中,奇数帧数据与偶数帧数据的采集方向可互换,其并不影响本实施例的实施。进一步的,各帧数据采集方向可任意选取,由于本发明在于利用邻近帧进行运动伪影矫正,因此相邻帧k空间数据采集方向不同即可。
接着,执行步骤S102,将所述各帧中由一次射频激发采集得到的数据分为一组,所述各帧内至少有一个数据组。如图3所示,以由第一帧k空间数据对第二帧k空间数据进行矫正为例,将每帧数据分为两组,粗线部分数据为一组,细线部分数据为一组。这里分组的目的在于,当采集高分辨率图像的k空间时,其对应k空间中的相位编码数目较多,因此为了采集足够的数据需要利用射频脉冲激发多次,每次激发后采集得到的数据为一组数据,因为一组内数据采集速度足够快,这样可以认为一组内的数据没有刚体运动发生。即图3中第一帧或第二帧的粗线或细线部分数据为一组,粗线组或实线组内数据没有刚体运动发生。在进行刚体运动矫正时,每次将前一帧矫正好的一组数据和后一帧任意一组未矫正的数据取出进行矫正。
然后,执行步骤S103,以前一帧中的一个数据组作为标准组,后一帧中的数据组作为待矫正组,利用所述标准组的数据对待矫正组的数据进行数据矫正,依次完成各帧数据组的数据矫正。
数据矫正的过程依次分为k空间中心矫正、旋转角度矫正以及平移矫正,其矫正方法与PROPELLER的技术思想类似。首先矫正其数据组中心与K空间中心的偏移。以第一帧为例,将其k空间原始数据DataZ1变换至图像域得到数据ImageZ1,再将DataZ1加一个二维的三角窗编程变换至图像域得到数据ImageZ1Window,其图像域数据ImageZ1和ImageZ1Window均为复数。令ImageZ1与ImageZ1Window进行相位相减,并取ImageZ1的模值,即令:
abs(Image)=abs(ImageZ1);
angle(Image)=angle(ImageZ1-ImageZ1Window);可以得到图像域数据Image,然后将Image变换回K空间,得到k空间数据DataZ1Corr,即为矫正好的数据,第二帧以及之后各帧的k空间数据均以同样的方法进行数据矫正,即完成了各帧的k空间中心矫正。
接下来进行旋转角度矫正,请参见图4。如图4所示,左上方图表示取出的第一帧内数据,第一帧数据采集方向为kx方向,其实线代表第一帧中取出的标准组数据,右上方图表示取出的第二帧内数据,第二帧数据采集方向为ky方向,其虚线代表第二帧中取出的任意一组待矫正组数据。其标准组数据和待矫正组数据在ky和kx方向大小分别为K1和K2,标准组数据在ky方向步长为ΔK1,待矫正组数据在kx方向步长为ΔK2。图4下方表示为公共网格空间,其中每个黑色圆点代表一个经过插值后得到的数据点,该公共网格空间在ky方向大小为K22,步长为ΔK22,在kx方向大小为K11,步长为ΔK11。
旋转角度矫正的过程为,首先将两帧的数据组数据取模值,然后将取模值后的第一帧数据插值到图4中黑色圆点表示的公共网格空间,插值后数据记为MB。之后将第二帧数据旋转一系列α的角度(α从0度到360度),并插值到图4中公共网格空间得到数据MCi(这里表示为第二帧待矫正组数据经过第i次旋转插值到公共网格空间的数据,i为正整数)。然后计算数据组MB与MCi的相关性(Correlation),求得各旋转系列中得到的相关性最佳值,即求得最大相关性值所对应的旋转角度α值,认为α即是两帧数据组在采集过程中运动伪影的旋转角度变化。将第二帧待矫正组数据按照旋转角度α进行修正,即完成了旋转角度矫正。
最后进行平移矫正。将标准数据组数据和经旋转修正的待矫正组数据共同插值到公共网格空间,分别得到数据组DB和DC,计算F(DB*conj(DC))的值,这里运算符F表示傅里叶变换,conj表示取共轭值,即将数据组DB与DC的共轭进行相乘,然后变换至图像域得到图像,之后将所得图像拟合找到最大值位置。其最大值所在位置与图像的中心位置的偏移量即为两帧数据采集时平移运动造成的偏移。根据计算得出的偏移量,将第二帧的待矫正数据组数据相对于第一帧标准数据组数据的平移移除,即完成了平移矫正。
这里以第一帧中的一个数据组作为标准组完成了第二帧中一个数据组的数据矫正,第二帧中的其他数据组可按照相同方式进行数据矫正(以第一帧中一个数据组作为标准组数据),之后其他各帧数据组均可按照类似的数据矫正方法进行。
进一步说明,这里为了计算方便,在以上进行插值选取公共网格空间时可以选取K11=K1,K22=K2,ΔK11=ΔK1,ΔK22=ΔK2,但也可以选取不相等的情况。至此,第二帧选出的待矫正组数据的运动偏移量矫正完毕。
本实施例的另一种实施方式中,也可以利用刚矫正完毕的第二帧数据组,来矫正第二帧其余数据组的刚体运动伪影。数据矫正过程中选取的标准组和待矫正组应尽量选取相邻帧(或相近帧)进行,这样可以保证运动的物体形变(非刚体运动)不大。
最后执行步骤S104,将数据矫正后的各帧数据组填回相应帧的k空间中,得到各帧矫正完毕的k空间数据。
这里需要说明的是,上述几个步骤中,如果某一种运动影响不大,则可在步骤中忽略。如心脏电影中的旋转运动影响不大,则可略去旋转矫正的步骤直接进行之后的平移矫正。
本实施例中对应上述k空间运动伪影矫正方法还提供了一种k空间运动伪影矫正装置,如图5所示,其包括:
采集单元101,适于在成像区域连续采集多帧k空间数据,所述各帧内的k空间数据由至少一次射频激发采集得到,相邻帧数据的采集方向不同;
分组单元102,与采集单元101相连,适于将所述各帧中由一次射频激发采集得到的数据分为一组,所述各帧内至少有一个数据组;
矫正单元103,与分组单元102相连,适于将前一帧中的一个数据组作为标准组,后一帧中的数据组作为待矫正组,利用所述标准组的数据对待矫正组的数据进行数据矫正,依次完成各帧数据组的数据矫正;
填充单元104,与矫正单元103相连,适于将矫正后的各帧数据组填回相应帧的k空间中,得到各帧矫正完毕的k空间数据。
在本实施例中,矫正单元103又包括:k空间中心矫正单元1031,适于进行k空间中心矫正;旋转矫正单元1032,与所述k空间中心矫正单元1031相连,适于进行旋转矫正;平移矫正单元1033,与所述旋转矫正单元1032相连,适于进行平移矫正。矫正单元1031-1033其中的一个或多个可省略。
所述k空间运动伪影矫正装置的具体实施可参考所述k空间运动伪影矫正方法的实施,在此不再赘述。
实施例二
在本实施例中,与实施例一的区别在于,本实施例中各帧k空间数据采集方式如图6所示,这里只示意出了第一帧k空间的采集数据,其数据经过三次射频激发采集得到,之后各帧的数据采集方式相同。因此在本实施例中,将每帧数据分为三组,每一组为一次射频激发所采集得到的数据,其每组数据分别对应如图6所示的实线部分、虚线部分和点划线部分。在本实际例的实施中,其具体分组情况根据数据采集情况而定,经过一次激发采集得到的数据即分为一组,其分组的数目并不限定于两组或三组。
其刚体运动矫正过程与实施例一中的方法相同,取出前一帧中的一组数据作为标准组数据,后一帧中的数据组作为待矫正组,利用标准组的数据对待矫正组的数据进行数据矫正,依次完成各帧数据组的数据矫正。以第一帧对第二帧进行矫正为例,从第一帧取出任一组数据作为标准组数据(实线、虚线或点划线任一组),第二帧未矫正数据任取出一组数据作为待矫正组数据(实线、虚线或点划线任一组),然后依次进行k空间中心矫正、旋转角度矫正以及平移矫正。若第一帧取出实线组作为标准组数据,第二帧取出实线组作为待矫正组数据,则待第二帧实线组矫正完毕后,第二帧虚线组和点划线组可按照第一帧实线组或者第二帧实线组进行矫正均可。待第二帧每组数据都矫正完毕后,接下来以同样的方式矫正第三帧、第四帧等数据。
本发明已以较佳实施例揭示如上,然其并非用以限定本发明,任何本领域技术人员,在不脱离本发明的精神和范围内,当可作些许的修改和完善,因此本发明的保护范围当以权利要求书所界定的为准。