发明内容
本发明的一个目的在于提供一种加速MR温度成像的方法,用来使MR温度成像兼具较好的时间分辨率和空间分辨率。
本发明的另一目的在于提供一种与上述加速MR温度成像的方法相应的装置。
为实现上述目的,本发明提出一种加速磁共振MR温度成像的方法,应用于MR监控的高强度聚焦超声HIFU治疗中,其包括:在MR温度成像过程中,确定超声波焦点处的温度变化;根据所述焦点处的温度变化确定数据采集所需的理想加速速率;根据所确定的理想加速速率调整K空间的变密度数据采样;重建采样得到的数据,得到温度图像。
其中,所述根据温度变化确定数据采集所需的理想加速速率包括:在温度变化快的情况下,确定较高速率作为理想加速速率;在温度变化慢的情况下,确定较低速率作为理想加速速率。
其中,在相位编码方向上,K空间被分为若干个区段,所述根据所确定的理想加速速率调整K空间的变密度数据采样包括:根据所确定的理想加速速率调整K空间变密度采样的区段划分和各区段的加速系数。
其中,所述根据所确定的理想加速速率调整K空间变密度采样的区段划分和各区段的加速系数包括:在温度变化快的情况下,减少加速系数较低的区段的大小和/或增加外部区段的加速系数;在温度变化慢的情况下,增加加速系统较低的区段的大小和/或降低外部区段的加速系数。
其中,所述根据所确定的理想加速速率调整K空间的变密度数据采样包括:如果在多个连续帧中应用相同的理想加速速率,对所述多个连续帧中的相位编码线进行交错采样。
所述方法还包括:
在重建采样得到的数据之前,在K空间的频率编码方向应用离散傅立叶逆变换,选择包含感兴趣区的相位编码线用于重建;则所述重建采样得到的数据包括:重建在所述感兴趣区中的相位编码线采样得到的数据。
其中,所述重建采样得到的数据包括:采用回波共享方法和/或迭代k-t方法来重建采样得到的数据。
其中,所述根据成像需求,采用回波共享方法或迭代k-t方法来重建采样得到的数据包括:在需要实时更新温度图像的情况下,采用回波共享方法重建所述采样得到的数据;在需要进行热剂量估计的情况下,采用迭代k-t方法重建所述采样得到的数据;在既需要实时更新温度图像又需要进行热剂量估计的情况下,同时采用回波共享方法和迭代k-t方法重建所述采样得到的数据。
其中,所述确定超声波焦点处的温度变化包括:预先估计超声波焦点处的温度变化。
其中,所述确定焦点的温度变化包括:定时反馈从重建图像中得到的焦点处的温度值;根据一段时间内焦点处的温度值的变化预测后续的温度变化。
其中,所述定时反馈从重建图像中得到的焦点处的温度值包括:在焦点处的温度值变化较大时,减小从重建图像中得到的焦点处的温度值的反馈周期;在焦点处的温度值变化不大时,增大从重建图像中得到的焦点处的温度值的反馈周期。
为实现上述目的,本发明还提出一种加速磁共振MR温度成像的装置,应用于MR监控的高强度聚焦超声HIFU治疗中,其包括:加速速率确定单元,用于在MR温度成像过程中确定超声波焦点处的温度变化,并根据所述焦点处的温度变化确定数据采集所需的理想加速速率;信号处理单元,用于根据所述加速速率确定单元所确定的理想加速速率调整K空间的变密度数据采样,并重建采样得到的数据,得到温度图像。
其中,所述加速速率确定单元包括:温度变化确定模块,用于在MR温度成像中确定超声波焦点处的温度变化;加速速率确定模块,用于根据所述温度变化确定模块确定的所述温度变化确定数据采集所需的理想加速速率。
其中,所述信号处理单元包括:数据采样模块,用于根据所述加速速率确定单元所确定的理想加速速率调整K空间的变密度数据采样;数据重建模块,用于对所述数据采样模块采集的数据进行重建,得到温度图像。
所述装置还包括:温度反馈单元,用于定时将从重建图像中测量的焦点处的温度值反馈至所述加速速率确定单元;所述加速速率确定单元根据一段时间内所述温度反馈单元反馈的温度值确定所述温度变化。
所述装置还包括:温度反馈控制单元,用于控制所述温度反馈单元向所述加速速率确定单元反馈温度值的周期。
从以上的技术方案可以看出,采用了本发明的方案,由于根据超声波焦点的温度变化来确定加速速率,并调整K空间VD数据采样,提高了加速MR温度成像的能力,兼具较好的时间分辨率和空间分辨率,具有较高的灵活性和稳定性。此外,由于可以在两种重建方法中灵活选择重建方法对采样数据进行重建,兼顾了数据重建的实时性和准确性,能够根据实际需求达到最佳的重建效果,使本发明的方案更加可行。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。
本发明提出的加速MR温度成像的方案,即在温度成像过程中,根据焦点处的温度变化确定数据采集所需的理想加速速率,根据所确定的理想加速速率调整K空间的VD采样,然后对采样的数据进行重建。优选地,在对采样数据进行重建时,根据不同的成像需求交替使用两种重建方法。
图4是根据本发明一实施例的加速MR温度成像的方法流程图。从图4可见,在本实施例中,加速MR温度成像的方法大致包括如下步骤:
步骤S401:根据对温度变化的预计情况确定后续数据采集所需的理想加速速率。
当一组帧的理想的加速速率可以事先确定,例如当温度变化大体能够预知或控制时,预先估计温度变化,然后根据该预先估计的温度变化确定一个固定的加速速率。
理想的加速速率也可以实时确定。例如,利用通过计算历史温度变化(例如通过回波共享方法反馈的温度值),然后利用一段时间内的温度值的变化预测后续的温度变化情况。根据该预测结果,适应性地调整后续数据采集所需的加速速率。这尤其适用于捕捉斜率较大的温度变化曲线的目的。
步骤S402:根据所确定的理想加速速率进行K空间VD(VD,Variable-Density,空间变密度)数据采样。
首先介绍K空间VD采样,在相位编码方向上,K空间被分为几个区段。图5是K空间VD采样技术的示意图。如图5所示,假设将K空间分为5个区段:R-2,R-1,R0,R1,R2,R0位于K空间的中心。在R-2和R2中,采集四个相位编码线中的一个,因此加速系数是4。类似地,在R-1和R1,加速系数是2。R0被全部采样。加速系数从K空间的内部向外部增加。
根据温度的变化速度调整区段划分以及相应的加速系数。例如,当温度改变较快时,可通过减少加速系数较低的区段的大小和/或增加外部区段的加速系数来进一步减少相位编码线,从而获得更高的时间分辨率。反之,当温度变化比较平缓时,相反的调整可达到更高的SNR(或重建图像的更高质量)。如果相同的K空间VD采样模式应用在几个连续的帧中,这些帧中的相位编码线可以被交错采样。图6示出的是交错采样的一个示例。如图6所示,横轴为相位编码方向,纵轴为时间轴,在时间轴上相邻几帧的相位编码线上所采样的信号是交错排列的。
步骤S403:根据实际需要交替或同时使用两种重建方法进行重建,得到温度图像。
首先,可以在频率编码方向应用离散傅立叶逆变换(IDFT,Inverse Discrete FourierTransform),然后选择包含感兴趣区的相位编码线用于重建。以避免重建整个图像,这可节省相当一部分不必要的工作。
所述两种重建方法包括回波共享(Echo Sharing)方法和迭代k-t(Iterative k-t)方法。
回波共享方法通过在相邻几帧中共享数据,用于需要实时更新温度图像和/或温度反馈的情况,采用回波共享方法,重建速度快,但获得的重建精度较低。迭代k-t方法可提供(精度?)更高的重建,在需要热剂量估计时,采用迭代k-t方法有助于精确地估计热剂量,但其计算速度非常低。由于这两种方法互相补偿,因此两者的适当的选择或者结合能够提供比较强大的性能。两种重建方法的选择主要取决于计算负载。如果迭代k-t方法无法在给定时间完成工作,则优选回波共享方法。应该理解,如果同时需要实时更新的温度图像和准确重建,可以同时使用两种重建方法进行重建,两者的结合能够提供更强大的性能。
另外,在本发明中,也可以引入温度控制机制,并通过温度控制机制来控制特定帧的理想的加速速率,获得可变的反馈速率。例如,在温度变化斜率较大的上升沿,需要快速的反馈来避免严重的过冲(overshoot),此时温度控制机制可任意调整温度反馈的速率,以通过设置加速速率来优化性能。
图7是根据本发明一实施例的加速MR温度成像的装置结构图。从图7可以看出,在本实施例中,加速MR温度成像的装置主要包括加速速率确定单元701和信号处理单元702。加速速率确定单元701在MR温度成像过程中确定超声波焦点处的温度变化,并根据焦点处的温度变化确定数据采集所需的理想加速速率;信号处理单元702根据加速速率确定单元701所确定的理想加速速率调整K空间的VD数据采样,并重建采样得到的数据,得到温度图像。
加速速率确定单元701包括温度变化确定模块7011和加速速率确定模块7012,其中:温度变化确定模块7011在MR温度成像中确定超声波焦点处的温度变化;加速速率确定模块7012,用于根据温度变化确定模块7021确定的温度变化确定数据采集所需的理想加速速率。
信号处理单元702包括数据采样模块7021和数据重建模块7022,其中:数据采样模块7021根据加速速率确定单元701所确定的理想加速速率调整K空间的VD数据采样;数据重建模块7022对所述数据采样模块7021采集的数据进行重建,得到温度图像。
优选地,该装置还包括温度反馈单元703,用于定时将从重建图像中测量的焦点处的温度值反馈至加速速率确定单元701;加速速率确定单元701根据一段时间内温度反馈单元703反馈的温度值确定温度变化。
优选地,该装置还包括温度反馈控制单元704,用于控制温度反馈单元703向加速速率确定单元701反馈温度值的周期。
本发明的方案具体以下优点:
高加速能力:通过利用k-t空间的冗余并结合迭代k-t方法,该新的加速方案可达到比其他仅在k空间或时域空间的方法更高的加速速率,同时保持相同的空间分辨率;
高灵活性:如果相同的K空间采样模式既应用在快速变化时间和慢速变化时间,存在两个问题中的一种或两种:一方面,时间精度太粗,以至于无法有效捕捉温度曲线;另一方面,由于在慢速变化时间内过分加速,重建图像会出现不必要的SNR损失或准确度损失。本发明的方案中采用的VD采样模式可根据温度变化速度适应性地调整加速速度。与现行的并行成像技术相比,本发明的方案可以获得连续变化的加速速率,从而避免加速速率的调整而引起较大的重建精度的波动。
可行性:灵活地结合两种重建方法,使本发明更可行,一方面,回波共享方法适用于实时更新的情况,但准确度低;另一方面,因为迭代k-t方法的一次重建需要大量帧,计算负载显著,很难用于实时监视。然而,迭代k-t方法比回波共享方法提供更精确的重建图像。VD采样模式和这两种重建方法的结合使得本发明的方案更可行;
高稳定性:在临床医疗应用中,超声加热监视中的稳定性是至关重要的。通常,大多数表示变化信号的能量聚焦在动态成像中K空间的中心,因此中心部分数据在稳定和准确重建中起特殊作用。然而,一些并行成像(Parallel Imaging)技术在均匀的K空间中采样。不合适的校准会导致严重的残留伪影出现在重建图像中。与这些技术不同,本发明的方案可工作在灵活的K空间采样模式下,通过在K空间的中心进行大量的更密的采样,本发明的方案与现有方案相比能够避免在重建图像中产生严重的残余伪影,提供更准确稳定的结果;
k-t相关性高:与时域空间的均匀采样不同,VD采样模式中,时间分辨率可根据温度变化速度改变,因此由于快速温度变化导致的下降的k-t相关性可通过时间分辨率的增加得到补偿,这有助于进一步改善准确性。图8示出了如何在温度变化快时通过在时域空间中的可变采样补偿k-t相关性。图8的左上图为时域均匀采样的曲线,左下图为时域均匀采样所获得的温度序列;右上图为时域VD采样的曲线,右下图为时域VD采样所获得的温度序列。从图8中时域均匀采样和时域VD采样所得到的温度序列的比较来看,采用时域VD采样,由于改善了时间分辨率,其所获得的温度序列的变化速率降低;
温度控制机制得到改善:VD采样为温度控制器提供了可控的反馈速率,以改善其性能,同时得到快速的温度变化,温度控制器可以任意根据特定目的调整反馈速率,虽然需要在精度和速度之间进行一些权衡。
下面介绍本发明方案的仿真测试结果。
温度模型:
从相位图像中提取温度变化:
此处,ΔT(x,y)是温度变化,Δφ(x,y)是MR图像中两帧的相对相位变化,δ是常数,δ=-0.01ppm/℃,γ是磁化率系数,B0是静磁场强度,TE是温度序列的回波时间。
温度的空间分布模型:
此处,f(t)是超声的焦点温度变化的函数,为了简化,设置σx和σy作为加热过程中的常数。
通过光纤温度测量法在HIFU超声的加热过程中测量焦点附近的722个温度值。采样率是每秒1个。图9是用光纤温度测量计测量到的HIFU加热温度变化过程曲线。所采用的fGRE序列的参数是:TR=70ms,TE=10ms。表示一秒中只能应用14个相位编码步骤,在722秒中总共应用722×14个相位编码步骤。应用相位编码时,在线性插值之后,每个数据表示近似温度值,在插值之后使用这些数据作为f(t)。
仿真加热区域:
由fGRE温度序列采集的浸在水中的牛的肝脏磁共振图像作为参考图像。通过公式(2)计算出的变化的空间温度轮廓通过公式(1)转换为相位轮廓并覆盖在磁共振图像上,以仿真加热过程。图10示出了一幅加热图像。通过DFT(Discrete Fourier Transform,离散傅立叶变换)将这些图像转换为原始数据。位于图10中心的矩形组织是牛的肝脏,环绕牛肝脏的水的信号密度较低。模拟的加热区域位于牛肝脏的中心箭头所示的位置。
几点假设:
a)温度在采集单个相位编码线时不变,因此体模在722秒中具有总共722×14种状态。加热区域的相位信息可逆地从其温度中计算出来;
b)图像矩阵为256×256,表示需要256个相位编码线进行全采样。
仿真测试的目的是测试时间分辨率如何影响快速温度变化的捕捉能力以及回波共享和迭代k-t方法交替使用的可行性。
在仿真时,在K空间全采样,验证捕获温度变化的能力(总共采集16帧,4096条相位编码线)。仿真结果显示在图11中。
图11中,粗实线表示真实的温度变化,细实线表示采集K空间中心的时间点,虚线表示测量的温度,该虚线是利用从采集的16帧数据中重建得到的16个温度值的线性插值拟合而成的。
从图11可得到几个结论:a)重建MR图像的测量值近似于采集K空间时的温度值;b)由于低时间分辨率,快速变化曲线无法有效捕获,如粗实线和虚线之间的差别。
采用本发明的方案,将自适应VD采样应用于捕获温度变化(总共16帧,2516相位编码线)。仿真结果显示在图12中。
图12中,粗实线表示真实的温度变化,细实线表示采集K空间中心的时间点,三角和虚线表示测量的温度,其中,三角表示从回波共享方法重建的图像中测量的温度,虚线表示从迭代k-t方法重建的图像中测量的温度。
比较图12和图11,两种采样模式都在低变化时间内达到高图像质量,然而本发明的VD采样可以适应性地调整时间分辨率以捕获快速的温度变化,而全采样模式损失了这些信息。图12也显示出本发明的重建方法的可行性和准确度。
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。