CN102488497A - 磁共振温度测量方法及系统 - Google Patents

Abstract

一种磁共振温度测量方法，包括如下步骤：进行稀疏采样获取导航数据集和图像数据集，并重建得到相应的相位图；通过所述相位图获取待测区域的相位变化值；根据所述相位变化值得到待测区域的温度变化值。上述磁共振温度测量方法及系统中，通过视野采集两组数据，并重建出高时空分辨率的图像，同时在重建过程中不影响其中的相位信息，得到相应的相位图，从而根据相位图测量出温度的变化值，实现了磁共振中对温度变化的实时监控，同时也达到了对图像高时间分辨率和高空间分辨率的要求。

Description

磁共振温度测量方法及系统

【技术领域】

本发明涉及磁共振成像技术，特别是涉及一种磁共振温度测量方法及系统。

【背景技术】

肿瘤热消融术具有创伤小、术后复发率低等优点而广泛应用于前列腺、子宫肌瘤的治疗中，肿瘤热消融术是通过加热使得病变组织温度升高并将其杀死，因此，为保证肿瘤热消融术在治疗中的顺利进行，必须快速准确地测量病变组织中的温度变化，进而医生才能了解到当前的治疗效果、预测坏死区域、及时修改治疗方案以及避免损伤病变区域周围的健康组织。

由于基于磁共振所进行的温度测量技术不仅能够对病变组织进行精确定位，还能够对肿瘤热消融治疗过程中的温度变化进行定量测量，因此，在传统的肿瘤热消融过程中主要通过磁共振监控成像来实现病变组织的温度测量。

传统的肿瘤热消融术运用了传统的扫描序列进行成像，这些扫描序列包括平面回波成像(Echo planar Imaging，简称EPI)以及梯度回波序列等，通过这些扫描序列进行图像采集所得到的图像虽然能够较好地保持图像中的相位信息，能够较为准确地通过保持的相位信息测量出温度的变化，但是在动态成像或者是对成像的时间和空间分辨率要求较高的情况下，并不能够达到要求，因此，无法实现对温度变化的实时监控。

【发明内容】

基于此，有必要提供一种能对温度变化进行实时监控的磁共振温度测量方法。

此外，还有必要提供一种能对温度变化进行实时监控的磁共振温度测量系统。

一种磁共振温度测量方法，包括如下步骤：

采集得到导航数据集和图像数据集，并重建得到相应的相位图；

通过所述相位图获取待测区域的相位变化值；

根据所述相位变化值得到待测区域的温度变化值。

优选地，所述采集得到导航数据集和图像数据集，并重建得到相应的相位图的步骤为：

通过部分可分离函数进行采样获取到导航数据集和图像数据集，并通过所述部分可分离函数依据所述导航数据集和图像数据集对图像进行重建得到相应的相位图。

优选地，所述通过所述相位图获取待测区域的相位变化值的步骤为：

对所述相位图中的待测区域进行相位的解卷绕得到所述待测区域的相位；

从所述相位图中获取所述待测区域的邻近区域的相位；

根据所述待测区域的相位与邻近区域的相位之间的差值计算得到相位变化值。

优选地，所述通过所述相位图获取待测区域的相位变化值的步骤为：

对所述相位图中的邻近区域进行相位的解卷绕得到所述邻近区域的相位；

根据所述邻近区域的相位进行多项式拟合得到待测区域的推测相位；

从所述相位图中获取待测区域的相位；

根据所述待测区域的相位与所述推测相位之间的差值计算得到待测区域的相位变化值。

优选地，所述根据所述相位变化值得到待测区域的温度变化值的步骤为：

根据所述待测区域的相位变化值按照图像中相位变化和温度变化之间的线性关系得到待测区域的温度变化值。

一种磁共振温度测量系统，包括：

成像模块，用于采集得到导航数据集和图像数据集，并重建得到相应的相位图；

相位处理模块，用于通过所述相位图获取待测区域的相位变化值；

温度计算模块，用于根据所述相位变化值得到待测区域的温度变化值。

优选地，所述成像模块还用于通过部分可分离函数进行采样获取到导航数据集和图像数据集，并通过所述部分可分离函数依据所述导航数据集和图像数据集对图像进行重建得到相应的相位图。

优选地，所述相位处理模块包括：

相位获取单元，用于对所述相位图中的待测区域进行相位的解卷绕得到所述待测区域的相位，从所述相位图中获取所述待测区域的邻近区域的相位；

变化计算单元，用于根据所述待测区域的相位与邻近区域的相位之间的差值计算得到相位变化值。

优选地，所述相位处理模块包括：

邻近相位计算单元，用于对所述相位图中的邻近区域进行相位的解卷绕得到所述邻近区域的相位；

拟合单元，用于根据所述邻近区域的相位进行多项式拟合得到待测区域的推测相位；

待测相位获取单元，用于从所述相位图中获取待测区域的相位；

相位差值处理单元，用于根据所述待测区域的相位与所述推测相位之间的差值计算得到待测区域的相位变化值。

优选地，所述温度计算模块还用于根据所述待测区域的相位变化值按照图像中相位变化和温度变化之间的线性关系得到待测区域的温度变化值。

上述磁共振温度测量方法及系统中，通过视野采集两组数据，并重建出高时空分辨率的图像，同时在重建过程中不影响其中的相位信息，得到相应的相位图，从而根据相位图测量出温度的变化值，实现了磁共振中对温度变化的实时监控，同时也达到了对图像高时间分辨率和高空间分辨率的要求。

【附图说明】

图1为一个实施例中磁共振温度测量方法的流程图；

图2为一个实施例中部分可分离函数的采集示意图；

图3为传统的全采样方式的示意图；

图4为一个实施例中通过部分可分离函数得到的重建图像；

图5为与图4相对应的相位图；

图6为通过传统的方式重建得到的图像；

图7为与图6相对应的相位图；

图8为图5中的相位变化图；

图9为图7中的相位变化图；

图10为一个实施例中通过相位图获取待测区域的相位变化值的方法流程图；

图11为一个实施例中待测区域和邻近区域的示意图；

图12为另一个实施例中通过相位图获取待测区域的相位变化值的方法流程图；

图13为一个实施例中磁共振温度测量系统的结构示意图；

图14为一个实施例中相位处理模块的结构示意图；

图15为另一个实施例中相位处理模块的结构示意图。

【具体实施方式】

在一个实施例中，如图1所示，一种磁共振温度测量方法，包括如下步：

步骤S110，进行采样获取导航数据集和图像数据集，并重建得到相应的相位图。

本实施例中，为了获得高时空分辨率的图像来对温度进行实时监控，，对视野进行采样以获取具有高时间分辨率的导航数据集和高空间分辨率的图像数据集，进而通过这两个数据集进行图像重建。在完成了数据的采集和重建之后，不仅可以得到幅值图，还能够得到对应的相位图。

在一个实施例中，上述步骤S110的具体过程为：通过部分可分离函数进行采样获取到导航数据集和图像数据集，并通过部分可分离函数依据导航数据集和图像数据集对图像进行重建得到相应的相位图。

本实施例中，在应用部分可分离函数(Partially Separable Functions，简称PSF)的数据采集和图像重建中，通常认为图像函数的空间变化和时间变化是L阶分离的，将图像数据表示为和时间上两个独立的变量函数

和

进而利用K空间和时间域上的信号相关来重建丢失数据。

具体地，通过K空间进行采样以得到导航数据集和图像数据集，根据所采集的导航数据集得到时间基函数和频率成分参数，由图像数据集得到空间基函数，然后通过时间基函数和空间基函数重建出完整的K空间数据。

在磁共振成像中，磁共振信号与自旋质子密度实际上是空间

和时间t的函数，即如以下公式所示：

$S(\stackrel{\→}{k},t)={\∫}_{+\∞}^{-\∞}\mathrm{\ρ}(\stackrel{\→}{r},t)e^{-i2\mathrm{\πk}\stackrel{\→}{r}}d\stackrel{\→}{r}+\mathrm{\η}(\stackrel{\→}{k},t)$

其中，

是K空间数据，在采集过程中常受到噪声η的影响，

是理想的图像域数据，为了以较高的时空分辨率恢复出部分可分离函数是认为图像函数

的空间变化和时间变化是L阶分离的，那么利用部分可分离函数的性质，

的L阶部分可分离函数模型可定义为：

为时间基函数，可以通过导航数据进行奇异值分解后得到。为空间基函数，可以由图像数据集

来获得。

如图2所示，部分可分离函数的图像序列主要由两个数据集组成，即高时间、低空间分辨率的导航数据集

以及高空间、低时间分辨率的图像数据集

在进行导航数据集和图像数据集的采集过程中应当满足如下三个条件：(1)脉冲重复时间T_R必须满足导航数据的时间奈奎斯特速率；(2)相位编码方向的采样间隔Δk_y必须满足图像数据集的空间奈奎斯特速率；(3)从图像数据集获取的采样帧数N必须大于或者等于阶数L。与图3所示的传统的全采样方式相比较，通过部分可分离函数的数据采集和图像重建中，数据量大大减少。

在完成了数据采集之后，可通过导航数据集和图像数据集得到时间基函数数

和空间基函数

导航数据集

可由如下矩阵C表示(t∈T₁)：

对矩阵C进行奇异值分解得到：

$C=\underset{l=1}{\overset{\min \{m,n\}}{\text{\Σ}}}{\mathrm{\λ}}_l{\mathrm{\μ}}_l{v}_l^H$

其中，{λ_l}为矩阵C的按照降序排列的奇异值，{μ_l}和{v_l}是矩阵C的左奇异值向量和右奇异值向量。

对矩阵C进行奇异值分解后得到左奇异值向量，即时间基函数

然后根据时间基函数

和图像数据集

通过最小二乘法得到空间基函数

即：

此时，在计算得到时间基函数

和空间基函数

之后，通过L阶部分可分离函数模型进行计算得到重建的K空间数据。

在完成了K空间数据的重建之后得到了相应的幅值图和相位图，为验证通过部分可分离函数的数据采集和重建得到的相位图是否影响到图像所包含的相位信息，从图4中可以查看到通过部分可分离函数重建得到的图像，图5为与之对应的相位图，图6为通过传统的方式重建得到的图像，图7为与其相对应的相位图，分别从图5和图7中抽取出相位信息，形成如图8和图9所示的相位变化图，从图中可以看出，通过部分可分离函数重建得到的图像相位变化和传统的方式得到的图像相位变化的趋势保持一致，由此说明通过部分可分离函数进行采样和重建并不会对图像相位变化产生影响。

步骤S130，通过相位图获取待测区域的相位变化值。

本实施例中，从相位图中获取相关的相位，以通过获取的相位进行计算得到待测区域的相位变化值。在肿瘤热消融术中，为实现病变组织的温度变化，加热区域即为待测区域，在相位图中获取该加热区域的相位变化值，以通过加热区域的相位变化值实现加热区域的温度变化监测。

在一个实施例中，如图10所示，上述步骤S130的具体过程为：

步骤S131，对相位图中的待测区域进行相位的解卷绕得到待测区域的相位。

本实施例中，在计算相频特性，即相位角时，需要用户到反正切函数ATAN2，反正切函数ATAN2的规定，第一和第二象限的角度为0～π，第三和第四象限的角度为0～-π，由此，若一个角度从0变到2π，但实际得到的结果是0～π，再由-π～0，在ω＝π处出现了跳变，跳变幅度为2π，这种现象称为相位的卷绕(wrapping)，为了得到连续的相频曲线，可在发生2π跳变的以后各处都加上(或减去)2π，这种做法称为相位的解卷绕(unwrapping)。解卷绕后相位在主值区间端点附近连续变化，从而得到连续的“自然态”的相位值。

步骤S132，从相位图中获取待测区域的邻近区域的相位。

本实施例中，邻近区域是相位图中与待测区域相接的区域，如图11所示，在肿瘤热消融术中，待测区域111和邻近区域113可以看作是隔离的。

步骤S133，根据待测区域的相位与邻近区域的相位之间的差值计算得到相位变化值。

本实施例中，基于质子共振频率的测温法，计算待测区域的相位和邻近区域的相位之间的差值即可得到该待测区域的相位变化值。

在另一个实施例中，如图12所示，上述步骤S130的具体过程为：

步骤S135，对相位图中的邻近区域进行相位的解卷绕得到邻近区域的相位。

本实施例中，在发生相同的相位变化的情况下，磁场B₀越稳定所测量的温度变化精确度越高，但是保证磁场稳定是较为困难的，尤其在低磁场的情况下，磁场漂移是很严重的，因此为进一步提高精确度，首先需要获取到邻近区域的相位。

步骤S136，根据邻近区域的相位进行多项式拟合得到待测区域的推测相位。

本实施例中，根据邻近区域的相位，通过拟合得到多项式Φ(x，y)，进而求出待测区域的原始相位，即推测相位Φ_extrapolated，多项式Φ(x，y)如以下公式所示：

Φ(x，y)≈a₀+a₁x+a₂y+a₃x²+a₄xy+a₅y²+...

其中，多项式中的系数a_i可通过最小二乘法得到。

步骤S137，从相位图中获取待测区域的相位。

步骤S138，根据待测区域的相位与推测相位之间的差值计算得到待测区域的相位变化值。

本实施例中，待测区域的相位Φ和推测相位Φ_extrapolated之差Φ-Φ_extrapolated即可作为待测区域的相位变化值，将待测区域的相位和推测相位的差值作为待测区域的相位变化值能够克服磁场不稳定所造成的影响，且测量的精度也比较高。

步骤S150，根据相位变化值得到待测区域的温度变化值。

本实施例中，基于质子共振频率的测温方法是通过相位变化来测量出温度的变化的，详细计算过程如以下公式所示：

$\mathrm{\ΔT}=\frac{\mathrm{\Δ\Φ}}{\mathrm{\γ\α}{B}_0{T}_E}$

其中，α＝-0.01ppm/℃，γ为旋磁比，B₀为磁场，T_E为回波时间。

在一个实施例中，上述步骤S150的具体过程为：根据待测区域的相位变化值按照图像中相位变化和温度变化之间的线性关系得到待测区域的温度变化值。

本实施例中，在质子共振频率的测温方法中，相位变化与温度变化之间是存在着线性关系的，因此测量得到准确的待测区域的温度变化在于计算得到准确的相位变化。

在一个实施例中，如图13所示，一种磁共振温度测量系统，包括成像模块10、相位处理模块30以及温度计算模块50。

成像模块10，用于进行采样获取导航数据集和图像数据集，并重建得到相应的相位图。

本实施例中，为了获得高时空分辨率的图像来对温度进行实时监控，成像模块10对视野进行采样以获取具有高时间分辨率的导航数据集和高空间分辨率的图像数据集，进而通过这两个数据集进行图像重建。在完成了数据的采集和重建之后，不仅可以得到幅值图，还能够得到对应的相位图。

在一个实施例中，成像模块10还用于通过部分可分离函数进行采样获取到导航数据集和图像数据集，并通过部分可分离函数依据导航数据集和图像数据集对图像进行重建得到相应的相位图。

本实施例中，在应用部分可分离函数(Partially Separable Functions，简称PSF)的数据采集和图像重建中，通常认为图像函数的空间变化和时间变化是L阶分离的，将图像数据表示为和时间上两个独立的变量函数

和

进而利用K空间和时间域上的信号相关来重建丢失数据。

具体地，成像模块10通过K空间进行采样以得到导航数据集和图像数据集，根据所采集的导航数据集得到时间基函数和频率成分参数，根据图像数据集得到空间基函数，然后根据时间基函数和空间基函数重建出完整的K空间数据。

在磁共振成像中，磁共振信号与自旋质子密度实际上是空间

和时间t的函数，即如以下公式所示：

$S(\stackrel{\→}{k},t)={\∫}_{+\∞}^{-\∞}\mathrm{\ρ}(\stackrel{\→}{r},t)e^{-i2\mathrm{\πk}\stackrel{\→}{r}}d\stackrel{\→}{r}+\mathrm{\η}(\stackrel{\→}{k},t)$

其中，

是K空间数据，在采集过程中常受到噪声η的影响，

是理想的图像域数据，为了以较高的时空分辨率恢复出

部分可分离函数是认为图像函数

的空间变化和时间变化是L阶分离的，那么利用部分可分离函数的性质，

的L阶部分可分离函数模型可定义为

为时间基函数，可以通过导航数据

进行奇异值分解后得到。

为空间基函数，可以由图像数据集

来获得。

在完成了数据采集之后，可通过导航数据集和图像数据集得到时间基函数数

和空间基函数

导航数据集

可由如下矩阵C表示(t∈T₁)：

对矩阵C进行奇异值分解得到：

$C=\underset{l=1}{\overset{\min \{m,n\}}{\text{\Σ}}}{\mathrm{\λ}}_l{\mathrm{\μ}}_l{v}_l^H$

其中，{λ_l}为矩阵C的按照降序排列的奇异值，{μ_l}和{v_l}是矩阵C的左奇异值向量和右奇异值向量。

对矩阵C进行奇异值分解后得到左奇异值向量，即时间基函数然后根据时间基函数和图像数据集

通过最小二乘法得到空间基函数

即：

此时，在计算得到时间基函数

和空间基函数

之后，通过L阶部分可分离函数模型进行计算得到重建的K空间数据。

相位处理模块30，用于通过相位图获取待的相位变化值。

本实施例中，相位处理模块30从相位图中获取相关的相位，以通过获取的相位进行计算得到待测区域的相位变化值。在肿瘤热消融术中，为实现病变组织的温度测量，加热区域即为待测区域，在相位图中获取该加热区域的相位变化值，以通过加热区域的相位变化值实现加热区域的温度变化监测。

在一个实施例中，如图14所示，上述相位处理模块30包括相位获取单元310以及变化计算单元320。

相位获取单元310，用于对相位图中的待测区域进行相位的解卷绕得到待测区域的相位，从相位图中获取待测区域的邻近区域的相位。

本实施例中，邻近区域是相位图中与待测区域相接的区域。在计算相频特性，即相位角时，需要用户到反正切函数ATAN2，反正切函数ATAN2的规定，第一和第二象限的角度为0～π，第三和第四象限的角度为0～-π，由此，若一个角度从0变到2π，但实际得到的结果是0～π，再由-π～0，在ω＝π处出现了跳变，跳变幅度为2π，这种现象称为相位的卷绕(wrapping)，为了得到连续的相频曲线，可在发生2π跳变的以后各处都加上(或减去)2π，这种做法称为相位的解卷绕(unwrapping)。解卷绕后相位在主值区间端点附近连续变化，从而得到连续的“自然态”的相位值。

变化计算单元330，用于根据待测区域的相位与邻近区域的相位之间的差值计算得到相位变化值。

本实施例中，基于质子共振频率的测温法，变化计算单元330计算待测区域的相位和邻近区域的相位之间的差值即可得到该待测区域的相位变化值。

在另一个实施例中，如图15所示，上述相位处理模块30包括邻近相位计算单元350、拟合单元360、待测相位获取单元370以及相位差值处理单元380。

邻近相位计算单元350，用于对相位图中的邻近区域进行相位的解卷绕得到邻近区域的相位。

本实施例中，在发生相同的相位变化的情况下，磁场B₀越稳定所测量的温度变化精确度越高，但是保证磁场稳定是非常困难的，尤其在低磁场的情况下，磁场漂移是很严重的，因此为进一步提高精确度，首先需要邻近相位计算单元350获取到邻近区域的相位。

拟合单元360，用于根据邻近区域的相位进行多项式拟合得到待测区域的推测相位。

本实施例中，拟合单元360根据邻近区域的相位，通过拟合得到多项式Φ(x，y)，进而求出待测区域的原始相位，即推测相位Φ_extrapolated，多项式Φ(x，y)如以下公式所示：

Φ(x，y)≈a₀+a₁x+a₂y+a₃x²+a₄xy+a₅y²+...

其中，多项式中的系数a_i可通过最小二乘法得到。

待测相位获取单元370，用于从相位图中获取待测区域的相位。

相位差值处理单元380，用于根据待测区域的相位与推测相位之间的差值计算得到待测区域的相位变化值。

本实施例中，待测区域的相位Φ和推测相位Φ_extrapolated之差Φ-Φ_extrapolated即可作为待测区域的相位变化值，相位差值处理单元380将待测区域的相位和推测相位的差值作为待测区域的相位变化值能够克服磁场不稳定所造成的影响，且测量的精度也比较高。

温度计算模块50，用于根据相位变化值得到待测区域的温度变化值。

本实施例中，基于质子共振频率的测温方法是通过相位变化来测量出温度的变化的，详细计算过程如以下公式所示：

$\mathrm{\ΔT}=\frac{\mathrm{\Δ\Φ}}{\mathrm{\γ\α}{B}_0{T}_E}$

其中，α＝-0.01ppm/℃，γ为旋磁比，B₀为磁场，T_E为回波时间。

在一个实施例中，温度计算模块50还用于根据待测区域的相位变化值按照图像中相位变化和温度变化之间的线性关系得到待测区域的温度变化值。

本实施例中，在质子共振频率的测温方法中，相位变化与温度变化之间是存在着线性关系的，因此测量得到准确的待测区域的温度变化在于计算得到准确的相位变化。

上述磁共振温度测量方法及系统中，与传统的序列成像监控温度变化的方法相比较，可以看出基于部分可分离函数能够在获得高时空分辨率图像的前提下，不损失图像相位信息，重建效果具有很强的鲁棒性，因此，根据质子共振频率的测温方法，根据重建图像后获得的相位信息，对温度进行准确的实时监控。

上述磁共振温度测量方法及系统中，通过视野采集两组数据，并重建出高时空分辨率的图像，同时在重建过程中不影响其中的相位信息，得到相应的相位图，从而根据相位图测量出温度的变化值，实现了磁共振中对温度变化的实时监控，同时也达到了对图像高时间分辨率和高空间分辨率的要求。

以上所述实施例仅表达了本发明的几种实施方式，其描述较为具体和详细，但并不能因此而理解为对本发明专利范围的限制。应当指出的是，对于本领域的普通技术人员来说，在不脱离本发明构思的前提下，还可以做出若干变形和改进，这些都属于本发明的保护范围。因此，本发明专利的保护范围应以所附权利要求为准。

Claims (10)

1.一种磁共振温度测量方法，包括如下步骤：

采集得到导航数据集和图像数据集，并重建得到相应的相位图；

通过所述相位图获取待测区域的相位变化值；

根据所述相位变化值得到待测区域的温度变化值。

2.根据权利要求1所述的磁共振温度测量方法，其特征在于，所述采集得到导航数据集和图像数据集，并重建得到相应的相位图的步骤为：

通过部分可分离函数进行采样获取到导航数据集和图像数据集，并通过所述部分可分离函数依据所述导航数据集和图像数据集对图像进行重建得到相应的相位图。

3.根据权利要求1所述的磁共振温度测量方法，其特征在于，所述通过所述相位图获取待测区域的相位变化值的步骤为：

对所述相位图中的待测区域进行相位的解卷绕得到所述待测区域的相位；

从所述相位图中获取所述待测区域的邻近区域的相位；

根据所述待测区域的相位与邻近区域的相位之间的差值计算得到相位变化值。

4.根据权利要求1所述的磁共振温度测量方法，其特征在于，所述通过所述相位图获取待测区域的相位变化值的步骤为：

对所述相位图中的邻近区域进行相位的解卷绕得到所述邻近区域的相位；

根据所述邻近区域的相位进行多项式拟合得到待测区域的推测相位；

从所述相位图中获取待测区域的相位；

根据所述待测区域的相位与所述推测相位之间的差值计算得到待测区域的相位变化值。

5.根据权利要求4所述的磁共振温度测量方法，其特征在于，所述根据所述相位变化值得到待测区域的温度变化值的步骤为：

根据所述待测区域的相位变化值按照图像中相位变化和温度变化之间的线性关系得到待测区域的温度变化值。

6.一种磁共振温度测量系统，其特征在于，包括：

成像模块，用于采集得到导航数据集和图像数据集，并重建得到相应的相位图；

相位处理模块，用于通过所述相位图获取待测区域的相位变化值；

温度计算模块，用于根据所述相位变化值得到待测区域的温度变化值。

7.根据权利要求6所述的磁共振温度测量系统，其特征在于，所述成像模块还用于通过部分可分离函数进行采样获取到导航数据集和图像数据集，并通过所述部分可分离函数依据所述导航数据集和图像数据集对图像进行重建得到相应的相位图。

8.根据权利要求6所述的磁共振温度测量系统，其特征在于，所述相位处理模块包括：

相位获取单元，用于对所述相位图中的待测区域进行相位的解卷绕得到所述待测区域的相位，从所述相位图中获取所述待测区域的邻近区域的相位；

变化计算单元，用于根据所述待测区域的相位与邻近区域的相位之间的差值计算得到相位变化值。

9.根据权利要求6所述的磁共振温度测量系统，其特征在于，所述相位处理模块包括：

邻近相位计算单元，用于对所述相位图邻近区域进行相位的解卷绕得到所述邻近区域的相位；

拟合单元，用于根据所述邻近区域的相位进行多项式拟合得到待测区域的推测相位；

待测相位获取单元，用于从所述相位图中获取待测区域的相位；

相位差值处理单元，用于根据所述待测区域的相位与所述推测相位之间的差值计算得到待测区域的相位变化值。

10.根据权利要求9所述的磁共振温度测量系统，其特征在于，所述温度计算模块还用于根据所述待测区域的相位变化值按照图像中相位变化和温度变化之间的线性关系得到待测区域的温度变化值。

Images