一种磁共振温度成像方法

技术领域

本申请涉及图像处理技术领域，更具体地说，涉及一种磁共振温度成像方法。

背景技术

磁共振温度成像(Magnetic Resonance Temperature Imaging，MRTI)可以实现无创、实时、在体监测被试物体内部温度分布及变化，在微创和无创热疗，例如磁共振间质热疗、聚焦超声治疗等治疗的监测过程中具有重要的用途。

当前的磁共振温度成像方法之一为基于质子共振频率(Proton ResonanceFrequency，PRF)位移的测温法，在实践过程中发现，基于质子共振频率位移的测温法受到磁共振线圈磁场均匀性，组织磁化率分布不均匀，呼吸/血流搏动等引起的组织运动等客观环境因素的影响较大，引起的常见误差包括相位解包裹误差，磁化率急速变化导致的误差，运动引起的误差等，容易造成最终获取的温度图与实际温度相差较大的问题，这使得该温度图失去了参考意义。

发明内容

为解决或缓解上述技术问题中的一项或更多项，本申请提供了一种磁共振温度成像方法及相关装置，以实现降低最终获得的温度图的误差，提高温度图的准确性的目的。

第一个方面，本申请实施例提供了一种磁共振温度成像方法，其特征在于，包括：

使用含有i个不同回波时间的梯度回波序列获取目标部位的磁共振影像数据，i为大于或等于2的正整数，所述磁共振影像数据包括与所述回波时间对应的相位图；

选取至少两组对应不同回波时间的相位图获得对应的多个温度差图；

根据所述温度差图获得温度图。

进一步地，本发明的方法中获得所述温度差图的步骤包括：

使用任一时刻的相位图减去基准时刻的相位图得到该时刻的相位差图，选取最小回波时间对应的相位差图作为参考相位差图，基于所述参考相位差图对其他回波时间对应的待校准相位差图进行校准，得到经校准的相位差图；

使用所述参考相位差图和所述经校准的相位差图计算温度差图。

进一步地，得到经校准的相位差图包括以下步骤：

根据相位差图与回波时间成正比例的关系，基于回波时间和参考相位差图，计算得到待校准相位差图的估计值；

使用所述估计值，根据相位周期性对所述待校准相位差图进行解包裹，得到经校准后相位差图。

可选地，本发明的方法中，参考相位差图对应的回波时间不超过18ms，17ms，16ms，15ms，14ms，13ms，12ms，11ms，10ms，9ms，8ms，7ms，6ms，5ms或4ms。

可选地，本发明的方法还包括消除磁共振系统引起的相位漂移的步骤，所述消除磁共振系统引起的相位漂移的步骤在相位差图或温度差图上进行。

可选地，本发明的方法还包括磁化率校正的步骤，磁化率校正的步骤在相位差图或温度差图上进行；

在温度差图上进行磁化率校正的步骤包括：

根据所述参考相位差图得到第一温度图，根据所述经校准相位差图得到对应的第二温度图；

判断所述第二温度图中每个像素对应的温度值与所述第一温度图中相应像素对应的温度值的差值的绝对值是否超过预设温度阈值，如果是，则对所述第二温度图中相应像素对应的温度值进行校正；

在相位差图上进行磁化率校正的步骤包括：

判断所述经校准相位差图中每个像素对应的相位差值与所述参考相位图中相应像素对应的相位差值的差值绝对值是否超过预设相位差阈值，如果是，则对经校准相位差图中相应像素对应的相位差进行校正。

可选地，本发明的方法还包括在相位差图或温度图上进行的校正运动引起的相位误差的步骤；

在相位差图上进行的校正运动引起的相位误差的步骤包括：

在相位差图在每个像素处的线性最小二乘拟合将运动引起的相位误差去除；

在温度差图上进行的校正运动引起的相位误差的步骤包括：

在温度差图的每个像素处用线性最小二乘拟合将运动引起的相位误差去除。

可选地，本发明的方法还包括获得加权温度图的步骤，该步骤使用所述参考相位差图和所述经校准的相位差图计算该时刻的温度差图：

使用所述参考相位差图和经校准的相位差图计算温度差图，并进行加权以获得待测目标的加权温度图。

另一方面，本发明还提供了一种磁共振温度成像系统，其包括：存储器和处理器；

所述存储器用于存储程序代码，所述处理器用于调用所述程序代码，所述程序代码用于执行本发明的磁共振温度成像方法。

进一步的，本发明的磁共振温度成像方法并非追溯性算法或迭代算法，运算量较小，可提供几乎实时的温度图，具有较高的参考意义。

附图说明

为了更清楚地说明本申请实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本申请的实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据提供的附图获得其他的附图。

图1为现有技术在离体和在体环境中获得的磁共振数据的幅值、相位图和温度图；

图2为本申请的一个实施例提供的一种磁共振温度成像方法的流程示意图；

图3为本申请的一个实施例提供的获取的相位图、相位差图和温度图；

图4为本申请的另一个实施例提供的一种磁共振温度成像方法的流程示意图；

图5为本申请的一个实施例提供的实验装置示意图；

图6为本申请的一个实施例提供的离体猪肉实验的激光间质热疗温度图的局部放大；

图7为本申请的一个实施例提供的在组织模拟物(图7(a))或离体猪肉(图7(b))实验期间，温度随时间的变化示意图；

图8为本申请的一个实施例提供的体内实验中的狗01的代表性温度图；

具体实施方式

磁共振温度成像可以引导多种能量输送型治疗手段，例如激光间质热疗、聚焦超声治疗、射频消融等，监控目标组织温度和治疗效果。发明人通过研究发现获取的温度图中错误的主要来源是相位解包裹错位导致的相位误差、磁化率误差和运动导致的相位误差。随着输入能量剂量的改变，磁化率畸变会导致图像振幅减小以及图像相位中的相应误差，从而破坏加热中心及其周围的重建温度图。重建温度图的错误可能会导致消融区域的估计错误，从而可能导致治疗效果的变化以及对关键组织的热损伤。因此，准确的温度成像对于治疗的有效性和安全性至关重要，尤其是应用于脑组织的消融区域时。

基于质子共振频率位移的测温法基于以下原理：氢质子的共振频率随水分子中的温度而变化。对于含水组织，局部磁场随温度的变化可描述为：

其中，α是随温度变化的质子共振频率系数，本发明中取0.008-0.015ppm/℃。受温度影响的水质子的相应共振频率变化可以表示为：

Δf＝αγB

其中，ΔT表示温度变化，Δf表示共振频率变化，γ表示旋磁比，B

可以在复杂的磁共振成像的相位中观察到由于温度变化引起的共振频率的变化。对于给定的梯度回波序列的间隔时间TE，可以根据相位差Δφ计算相对温度变化ΔT，该方程可表示为：

梯度召回回波(Gradient recalled echo,GRE)脉冲序列，简称梯度回波序列，是基于质子共振频率位移的测温法中使用的序列，包括对应不同回波时间的序列，例如可以包括i个回波时间对应的回波序列，i为不小于2的正整数。根据公式(3)可知，梯度回波序列中回波时间越长的序列，相同的温度变化可能导致相位差越大，可以获得更高的温度灵敏度。

参考图1，示出了现有技术中，使用单一回波时间的回波序列获得的温度图存在的多种常见问题，随着梯度回波序列的回波时间增加，相位对比和相位包裹都增加，这表明使用较长回波时间的序列，温度灵敏度更高，相位包裹更多。图1中，在(a)(离体，猪脑)和(b)(体内)中通过本申请实施例使用的含有4个不同回波时间的梯度回波序列获得的第一至第四回波的幅值(上排)和相位(第二行)。使用传统的PRF算法根据每个TE(回波时间)设置计算温度图(下排)。在较长的回波时间内会出现更多的相位包裹，同时图像对比度也相应增加。强烈的加热会由于磁化率的变化而导致信号损失，并且还会转化为加热中心周围像素的相位和温度误差。在体内实验中，脑脊液(Cerebrospinal Fluid，CSF)运动会在磁共振温度成像上引起不正确的高温，这在前一个回波中更明显，因为较短的TE对所引入的相似相位误差的耐受性较低。

具体地，由于大脑中脑脊液(CSF)的运动，扫描间运动是基于质子共振频率位移的测温法测量的温度图中的一个问题。脑脊液的幅值和相位信号经常通过脑脊液的正常动态运动在脉冲梯度回波序列上改变，这可能会混淆温度估计。脑脊液运动还可能导致心室内和周围心室的像素移动，从而导致相差图出现错误。如图1b所示，体内温度图显示由于脑脊液运动，第三脑室内的虚假高温。温度误差在较短回波时间的脉冲梯度回波序列上更为明显，因为根据公式(3)，它们对脑脊液流动引入的相移强度的容忍度更小。

实际上，水质子的局部磁场也应考虑磁化率x

其中，

激光加热会在激光尖端周围的GRE成像中引起明显的磁化伪影。仍然参考图1，温度急剧变化的加热中心(如图1(a)中箭头所示)在较长的回波时间幅值上显示严重的信号损失。体素内自旋相移(intra-voxel spin dephasing)是由温度和磁化率变化引起的局部磁场不均匀引起的。

加热引起的磁化伪影，尤其是较长回波时间的梯度回波序列对应的图像中的磁化伪影，是造成误差的重要原因。仍然参考图1，在离体或体内实验中，加热中心周围的相位误差转化为磁共振热成像上的人为低温。通常情况下，在梯度回波磁共振成像中，建议使用具有尽可能短的回波时间延迟时间以最大程度地减小磁化率伪影。但是，更长的回波时间的梯度脉冲序列可以提供更好的温度灵敏度和信噪比，这是当前面临的两难选择。

为了兼顾温度灵敏度、信噪比和低误差，本申请实施例提供了一种磁共振温度成像方法，包括：

使用含有i个不同回波时间的梯度回波序列获取目标部位的磁共振影像数据，i为大于或等于2的正整数，所述磁共振影像数据包括与所述回波时间对应的相位图；

选取至少两组对应不同回波时间的相位图获得对应的温度差图；

根据所述温度差图获得温度图。

其中，获得所述温度差图的步骤包括：

使用任一时刻的相位图减去基准时刻的相位图得到该时刻的相位差图，选取最小回波时间对应的相位差图作为参考相位差图，基于参考相位差图对其他回波时间对应的待校准相位差图进行校准，得到经校准的相位差图；

使用参考相位差图和经校准的相位差图计算温度差图。

梯度回波序列的回波时间与磁化率伪影的大小成正比关系，因此最小回波时间对应的梯度回波序列获得的相位图受到的由于加热而导致的磁化率变化的影响最小，其图像数据仍保持正确的相位，因此可以选取最小回波时间对应的相位差图作为参考相位差图对其他相位图进行校准，以降低最终获得的温度图的误差，提高获得的温度图的准确性的目的。

进一步的，所述磁共振温度成像方法并非追溯性算法或迭代算法，运算量较小，可提供几乎实时的温度图，具有较高的参考意义。

本发明的磁共振温度成像方法的一个具体实例，如图2所示，包括：

S101：使用含有i个不同回波时间的梯度回波序列对待测目标进行扫描，得到不同时刻与所述回波时间对应的i组相位图，所述i为大于或等于2的正整数；

在步骤S101中，所述梯度回波序列中回波时间的最小取值和最大取值均可根据实际需求而定，一般情况下，为了尽量降低磁化率变化导致的磁化率伪影，所述梯度回波序列中回波时间的最小取值可以取磁共振温度成像设备能够取到的最小值，所述梯度回波序列中回波时间的最大取值一般不超过对待测目标成像的回波时间的取值范围上限。例如，在对头部成像来说，其可选的梯度时间序列的回波时间的取值范围为3～30ms，所述梯度回波序列中包含的回波时间的具体取值均要在该取值范围内。

含有i个不同回波时间的梯度回波序列信息的获取可以是从服务器或其他存储设备中读取或接收，也可以是根据工作人员的设定实时获取的，本申请对获取多个梯度回波序列的具体方法并不做限定，具体视实际情况而定。

参考图3，具体参考图3(a)，获取的相位图如图3(a)所示。

S102：选取任一时刻的至少两组对应不同回波时间的相位图获得对应该时刻的相位差图，参见图3(b)，其中示出了相位包裹；接着，选取最小回波时间对应的相位差图作为参考相位差图，基于参考相位差图对其他回波时间对应的待校准相位差图进行校准，得到经校准的相位差图，结果参见图3(c)，消除了相位包裹；

S103，然后根据相位差图获得温度图，参见图3(d)，其中总示出了磁化率引起的误差。

下面对本申请实施例提供的磁共振温度成像方法的各个步骤的可行执行方式进行描述。

在上述实施例的基础上，在本申请的一个实施例中，获得温度差图的具体步骤包括：

使用任一时刻的相位图减去基准时刻的相位图得到该时刻的相位差图，基准时刻为对目标组织传输能量(例如热能、光能、射频消融、冷冻消融)之前的任意时刻，优选地为进行能量传输之前不久的时刻，例如即将进行能量传输的时刻；

可选的，参考相位差图相位解包裹的回波时间的取值小于或等于12ms，优选地，不超过10ms，更优选地，不超过6ms。

所述选取在该时刻最小回波时间对应的相位差图作为参考相位差图，对其他回波时间对应的待校准相位差图进行校准包括以下步骤：

使用所述参考相位差图和待校准的回波时间对应的相位差图，根据相位差与回波时间成正比例的关系，基于回波时间和参考相位差图的相位差，计算得到所述待校准的回波时间对应的相位差图的相位差的估计值；然后使用所述估计值，根据相位周期性对所述待校准的相位差进行解包裹，得到经校准后的相位差。

在上述实施例的基础上，在本申请的另一个实施例中，所述磁共振温度成像方法还包括：

S104消除磁共振系统引起的相位漂移的步骤，该步骤在相位差图或温度图上进行。

在相位差图上进行消除磁共振系统引起的相位漂移的步骤包括：

选取多个热参考点(Region of Interest，ROI)，通过从每个相位差图减去所述热参考点的平均相位差，物理温度稳定无变化且组织均匀的区域均可作为热参考点；

在温度差图上消除磁共振系统引起的相位漂移的步骤包括：

在温度差图中减去任一所述热参考点的平均温度差，进行校正。

在上述实施例的基础上，在本申请的另一个实施例中，所述磁共振温度成像方法还包括：

S105：磁化率校正的步骤，磁化率校正的步骤在在相位差图或温度图上进行；

在温度图上进行磁化率校正的步骤包括：

根据所述参考相位差图得到第一温度图，根据所述经校准相位差图得到对应的第二温度图；

可选地，在所述第一温度图和各个所述第二温度图中确定预设区域；

判断所述第二温度图中每个像素对应的温度值与所述第一温度图中所述预设区域中相应像素对应的温度值的差值的绝对值是否超过预设温度阈值，如果是，则对所述第二温度图中相应像素对应的温度值进行校正；校正可以有多种方法，例如可以使用第一温度图的温度值替换第二温度图的温度值，或者使用第二温度图中相邻像素的温度值替换第二温度图中该像素的温度值，或者基于相邻像素的温度值和第一温度图的温度值拟合一个近似温度替代第二温度图的温度值；

在相位差图上进行磁化率校正的步骤包括：

在所述参考相位差图和经校准的相位差图中确定预设区域；

判断所述经校准相位差图中每个像素对应的相位差值与所述参考相位图中所述预设区域中相应像素对应的相位差值的差值绝对值是否超过预设相位差阈值，如果是，则对经校准相位差图中的相位差进行校正。

在上述实施例的基础上，在本申请的又一个实施例中，所述磁共振温度成像方法还包括：

S106：在相位差图或温度图上进行的校正运动引起的相位误差的步骤；

在相位差图上进行的校正运动引起的相位误差的步骤包括：

通过使用所述参考相位差图和所述经校准相位差图在每个像素处的线性最小二乘拟合将运动引起的相位误差去除；

在温度图上进行的校正运动引起的相位误差的步骤包括：

根据所述参考相位差图得到第一温度图，根据所述经相位差图得到对应的第二温度图；

通过使用所述第一温度图和所述第二温度图在每个像素处的线性最小二乘拟合将运动引起的相位误差去除。

仍然参考图4(c)，图4(c)示出了在没有(左图)和有(右图)运动误差(MotionError)校正的情况下作为时间的函数的相差(第一行)和相对温度变化(第二行)。对于较短的回波时间，相位误差Δφ(x,y)

如前文所述，步骤S104、S105和S106均即可在相位差图层面上进行，也可在温度图上进行。即所述消除磁共振系统引起的相位漂移的步骤在相位差图和/或温度图上进行，所述磁化率校正的步骤在相位差图和/或温度图上进行。所述校正运动引起的相位误差的步骤在相位差图和/或温度图上进行。

仍然参考图3，图3(c)表示经校正的相差图，图3(d)表示温度图，图3(e)表示磁化率校正后的图像，图3(f)表示运动误差校正后的图像。

在上述实施例的基础上，在本申请的又一个实施例中，所述使用所述参考相位差图和所述经校准的相位差图计算该时刻的温度图包括：

S1031：使用所述参考相位差图和经校准的相位差图计算温度，并对计算得到的温度进行加权以获得待测目标的温度图，加权可以在温度差图或温度图二者上进行，加权可以是各种加权方法，例如平均加权，或者可以为单独的一个回波时间对应的温度图，即该温度图的加权系数为1，其他相位温度图的加权系数为0。

在步骤S1031之后还可包括：

S107：对所述待测目标的温度图进行多次插值处理，并利用插值处理后的所述待测目标的温度图计算消融区域边界。

对所述待测目标的温度图进行多次插值处理的目的是为了获得更平滑的消融区域边界，差值处理的具体次数可以是2或3次。

在计算消融区域边界的过程中，具体利用如下公式：

其中，E

下面结合具体实例对本申请实施例提供的磁共振温度成像方法进行验证。

使用包括10W，980nm二极管激光器和冷却的激光施加器系统的激光消融系统对组织模拟物(凝胶体模)进行加热。相位图像是在3T MR扫描仪(Ingenia，PhilipsHealthcare，Best，荷兰)使用16条接收线圈使用多回波时间梯度回波序列获取的：翻转角＝30°，TE＝6/12/18/24ms，TR＝22ms，矩阵＝176×176，FOV＝200x200mm

如图5所示，还将两个MR兼容的光纤温度探头插入到组织模拟物中，探头尖端位于靠近消融光纤的位置，以获取各点的凝胶温度。由于加热过程中光纤探头受到消融光纤的影响，因此温度计仅监视冷却阶段。图5中。具体地，图5中显示在组织模拟物中插入了消融光纤和两个光纤温度探头，填充凝胶的参考管固定在周围作为绝缘参考。

使用与组织模拟物实验相同的扫描参数进行猪肉和猪脑的离体实验。对每种类型的组织(凝胶，猪肉，猪脑)进行了两次实验，其中一种进行了若干次激光循环加热，另一种进行了持续的加热和冷却。计算出MR测量温度和光纤测量温度之间的均方根误差，作为温度精度的测量值。

杜宾狗的体内实验已获得清华大学伦理审查委员会的批准。九只成年杜宾狗接受了激光间质热疗的热疗。加热过程在3T MR扫描仪(Ingenia，Philips Healthcare，Best，荷兰)上通过32条接收头线圈使用多回波时间梯度回波序列进行监控。

仍然参考图3和图4，图3和图4说明了本申请实施例提供的磁共振温度成像方法的一个实例。图3示出了磁共振温度成像的一个具体实例步骤。图3(a)，在激光热疗期间获得的一个时间，首先通过多TE回波序列获得线圈组合相位图像；图3(b)，然后获得相位差图，白色箭头指示在加热中心周围的相位图上发生的相位包裹；图3(c)显示经过相位解包裹和B0漂移校正的相位差图。图3(d)根据PRF偏移方法从图3(c)计算得到的温度图。白色箭头突出显示了磁化率引起的误差。图3(e)磁化率校正后的温度图。白色箭头显示了残留的脑脊液运动(motion)引起的误差。图3(f)运动校正的温度图。

图4中，关于代表性像素的一个示例性方法流程包括：步骤1，获取相位差图(左上)以及参考相位差图估算的解包裹获得的相位图(右上)，如黑色箭头所示，在温度快速变化的情况下，18ms和24ms回波时间对应的相位差图发生了包裹。步骤2，获取相位解包裹图，步骤3，静态磁场强度(B

模体和离体实验结果：

图6示出了激光间质热疗期间离体猪肉实验的代表性温度图。从热循环期间获取的300帧(3s/帧)中选择六个代表性图像(#50表示第50帧，#146表示第146帧，依次类推)。第一行和第二行分别是使用传统的相位解包裹方法和本申请实施例中提出的基于多回波时间的相位解包裹方法。使用现有技术的相位展开方法，由于激光热量引起的磁化率变化，温度图上的像素会严重损坏，即使不再使用激光也无法恢复。技术原理如下：将现有技术的相位展开方法应用于时间维度进行相位跳变检测，如果当前帧的相差图被错误地解包裹，则随后的所有帧都会受到影响。另一方面，本发明提出的磁共振温度成像方法是在多回波维度的基础上进行的，因此避免了来自先前帧的干扰。第三行是经过相位解包裹和磁化率校正的单个回波时间温度图，加热中心周围的损坏像素已正确恢复。最后一行是使用本申请实施例提供的磁共振温度成像方法进行的多回波时间数据联合结果。显示出最终的磁共振热成像在热点的温度更均匀。

图7示出了在组织模拟物(图7(a))或离体猪肉(图7(b))实验期间随时间变化的温度，分别由两个测温光纤测量(红色线条)和本申请实施例提供的方法计算得出的(虚线黑色线条)。在多次加热(图7(a))或单次加热(图7(b))的情况下，质子共振频率(PRF)计算的温度-时间行为在降温阶段与测温光纤测得的温度-时间非常匹配。表2列出了MR计算值和测温光纤测量值之间的均方根误差(RMSE)值，这些值代表了所提出算法的温度精度。实验1进行了几次激光循环加热，而实验2则是连续的加热和冷却阶段。结果表明，在大多数情况下，凝胶，猪肉或猪脑组织的均方根误差均小于0.5℃。

表2.本申请实施例提供的方法在光纤测量的温度和MR计算的温度之间的比较。

缩写：RMSE，均方根误差；实验，实验L(R)，左侧(右侧)的光纤温度探头。

图8显示了体内实验中的狗01的代表性温度图。需要注意的是，消融区域位于靠近第三脑室和侧脑室的位置。选择在激光消融过程中获取的100帧(3s/帧)图像叠加在消融后T2w磁共振热成像上。从上至下是通过现有技术算法分别根据单个回波时间(TE)数据(TE＝6ms和TE＝24ms)计算出的温度图，并使用本发明所提出的算法根据多(联合)TE回波序列计算出的温度图。第一行(TE＝6ms)显示了第三侧脑室和侧脑室内的人为高温，表明短TE计算温度严重受CSF流动伪影的影响。CSF诱发的第三脑室内伪影(白色箭头指示)仍存在于第二行(TE＝24ms)，但被所提出的多TE回波序列算法很好地抑制了。第二行显示，与较短的TE(TE＝6ms)相比，更长的TE(TE＝24ms)可以提供更平滑的边界和更好的温度SNR，但是如上所述，由于磁化率的变化，加热中心周围的像素会损坏。另一方面，我们提出的方法整合了多个回波的信息，因此获得的温度图同时消除了CSF引起的误差和磁化率引起的误差，显示出更均匀和对称的加热区域。

上面的实验结果表明，可以使用本申请实施例提供的磁共振温度成像方法校正由于加热激光器本身引起的质子共振频率温度图中磁化率引起的误差。我们首先提出通过质子共振频率移位方法将多回波时间梯度回波脉冲序列应用于磁共振热成像。代替单回波序列，多梯度回波序列可提供更多信息，而无需额外的扫描时间，并且为相位解包裹和伪影消除提供了新方法。

较短的回波时间可以忍受磁化率伪像，但对噪声敏感，而较长的回波时间则具有较好的温度敏感性和信噪比，但受磁化率伪像的影响很大。本发明提出的磁共振温度成像方法融合了不同回波的优势，获得了更好的温度图测量结果。而且，本发明的磁共振温度成像方法可以显着提高磁共振热成像的鲁棒性和信噪比，从而避免了由于误估了低温而对健康组织造成的损害。

本发明的方法还具有出色的脑脊液流动模体抑制能力，并且可以在心室内或周围提供准确的温度测量。补偿脑脊液运动引起的错误对于激光间质热疗治疗脑室周围脑部病变在临床上很重要。此外，提出的算法是在线兼容的，不需要迭代计算，因此非常适合于磁共振热成像，因为需要非常接近实时的温度图。

下面对本申请实施例提供的磁共振温度成像系统进行描述，下文描述的磁共振温度成像系统可与上文描述的磁共振温度成像方法相互对应参照。

相应的，本申请实施例还提供了一种磁共振温度成像系统，其特征在于，包括：存储器和处理器；

所述存储器用于存储程序代码，所述处理器用于调用所述程序代码，所述程序代码用于执行本发明的磁共振温度成像方法。

综上所述，本申请实施例提供了一种磁共振温度成像方法及相关装置，其中，所述磁共振温度成像方法基于含有i个不同回波时间的梯度回波序列获得i组相位图，选取至少两组对应不同回波时间的相位图获得对应的多个相位差图，并根据所述温度差图获得温度图。发明人经过研究发现，梯度回波序列的回波时间与磁化率伪影的大小成正比关系，因此与较小回波时间对应的梯度回波序列获得的相位图受到的由于加热而导致的磁化率变化的影响最小，其图像数据仍保持正确的相位，因此可以基于含有i个不同回波时间的梯度回波序列获得的i组相位图以及相位差图进行温度图的获取，以降低最终获得的温度图的误差，提高获得的温度图的准确性的目的。

进一步的，所述磁共振温度成像方法并非追溯性算法或迭代算法，运算量较小，可提供几乎实时的温度图，具有较高的参考意义。

下面将结合本申请实施例中的附图，对本申请实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本申请一部分实施例，而不是全部的实施例。基于本申请中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本申请保护的范围。

本说明书中各实施例中记载的特征可以相互替换或者组合，每个实施例重点说明的都是与其他实施例的不同之处，各个实施例之间相同相似部分互相参见即可。

对所公开的实施例的上述说明，使本领域专业技术人员能够实现或使用本申请。对这些实施例的多种修改对本领域的专业技术人员来说将是显而易见的，本文中所定义的一般原理可以在不脱离本申请的精神或范围的情况下，在其它实施例中实现。因此，本申请将不会被限制于本文所示的这些实施例，而是要符合与本文所公开的原理和新颖特点相一致的最宽的范围。