[发明专利]一种背景相位提取方法及系统

摘要

本发明适用于通信领域，提供了一种背景相位提取方法，所述方法包括如下步骤选择采集梯度回波序列的多回波时间TE中任一TE；计算任一TE的滤除磁化率干扰的相位图；对该相位图进行多项式拟合得到系数，将系数代入相位图表达式得到相位图的表达式，依据该表达式计算得到被拟合区域的参考相位；计算任一TE相位图中被拟合区域的原始相位和参考相位的复数差通过计算得到温度变化值ΔT，该ΔT即为最终的温度图像。本发明具体实施方式提供的技术方案具有消除运动伪影和磁化率伪影等的影响，测量精度准确的优点。

主权项

一种背景相位提取方法，其特征在于，所述方法包括如下步骤：S11、选择采集梯度回波序列的多回波时间TE中任一TE；S12、计算任一TE的滤除磁化率干扰的相位图其中，x,y为任一TE相位图的拟合区域像素的坐标值；S13、对该相位图进行多项式拟合得到系数{Wn(m)}，将{Wn(m)}代入相位图表达式得到相位图的表达式，依据该表达式计算得到被拟合区域的参考相位其中，n为x的阶次，m为y的阶次；S14、计算任一TE相位图中被拟合区域的原始相位和参考相位的复数差通过计算得到温度变化值ΔT，该ΔT即为最终的温度图像；所述S12的实现方法具体可以为：δB*＝AχB*；A＝F‑1FdF；χB*=argminχB||w(δ-d⊗χB)||;]]>其中，为滤除磁化率干扰的相位图，为多回波序列采集的任一TE的相位图；δB*为过滤掉磁化率干扰的背景局部磁场，F为Fourier变换矩阵，F‑1为F的逆矩阵，Fd为d(r)的Fourier矩阵，d(r)为沿着主磁场B0方向的单位磁偶极矩磁场分布，χB为非组织区域磁环率分布，δ为由磁化率分布χ(r)引起的磁场不均匀性，r为向量，w为幅值图像；通过多回波相位图像得到局部场图：其中，为初始相位，γ为磁旋比，δB为磁场的不均匀性，δB由公式2计算得到。