1.一种背景相位提取方法，其特征在于，所述方法包括如下步骤：S11、选择采集梯度回波序列的多回波时间TE中任一TE；

S12、计算任一TE的滤除磁化率干扰的相位图 ；

S13、对该相位图 进行多项式拟合得到系数{Wn(m)}，将{Wn(m)}代入相位图表达式得到相位图 的表达式，依据该表达式计算得到被拟合区域的参考相位；

S14、计算任一TE相位图中被拟合区域的原始相位和参考相位的复数差 通过计算得到温度变化值ΔT，该ΔT即为最终的温度图像。

2.根据权利要求1所述的方法，其特征在于，所述S12的实现方法具体可以为：公式1；

其中， 为滤除磁化率干扰的相位图， 为多回波序列采集的任一TE的相位图； 为过滤掉磁化率干扰的背景局部磁场。

通过多回波相位图像得到局部场图， 公式2；其中， 为初始相位，γ为磁旋比，δB为磁场的不均匀性，δB由公式2计算得到；

公式3；

公式4；

θrr'表示向量r、r'之间的夹角，d(r)为沿着主磁场B0方向的单位磁偶极矩磁场分布；

公式5

公式6

δ＝F-1FdFX＝AX公式7其中，F为Fourier变换矩阵,F-1为F的逆矩阵，Fd为d(r)的Fourier矩阵；

χB非组织区域磁环率分布；A=F-1FdF。

3.根据权利要求1所述的方法，其特征在于，所述S13具体为：将 代入相位图表达式公式8进行多项式拟合得到系数{Wn(m)}；

公式8

其中，x,y任一TE相位图的拟合区域像素的坐标值，n为x的阶次，m为y的阶次；

将任意TE相位图的拟合区域像素的坐标值代入公式8中计算得到被拟合区域的参考相位 。

4.根据权利要求1所述的方法，其特征在于，所述S14的方法具体为：公式9；

其中 为任一TE相位图中被拟合区域的原始相位；

公式10

是加热时的相位； 是拟合得到的参考背景相位；α是质子共振频率偏移系数，对于大部分组织α＝-0.01ppm/℃；γ是磁旋比；B0是主磁场强度。