涡流校正方法、装置、移动终端及可读存储介质制造方法及图纸

本发明专利技术适用于属于磁共振序列设计领域，提供了一种涡流校正方法、装置、移动终端及可读存储介质，包括：步骤S1，利用双极性读出梯度回波序列采集多回波图像；步骤S2，估计采集的所述多回波图像中的涡流引入的额外相位项的一阶项系数；步骤S3，去除估计的所述一阶项系数，并估计采集的所述多回波图像中的涡流引入的额外相位项的零阶项系数；步骤S4，根据估计的所述一阶项系数和所述零阶项系数去除涡流引入的所述额外相位项的误差；本发明专利技术提供的涡流校正方法可以去除采集的图像中的涡流引起的相位误差，从而保证后续水脂分离算法结果的正确性。

【技术实现步骤摘要】
涡流校正方法、装置、移动终端及可读存储介质
本专利技术属于磁共振序列设计领域，尤其涉及一种涡流校正方法、装置、移动终端及可读存储介质。
技术介绍
基于磁共振成像的水脂分离或脂肪定量具有重要的临床意义。多回波梯度回波(GradientRecalledEcho，GRE)序列是最常见的用于水脂分离的磁共振成像序列。该序列中，为了减小回波间距和重复周期，通常采用双极性读出梯度采集多回波图像，但双极性读出方式往往受到系统涡流的影响，使得奇偶回波相位不一致，多回波数据不满足理论水脂混合模型，导致后续基于复数模型的水脂分离算法失败。具体地，含有两种成分(常见为水和脂肪)的理论水脂混合模型为：其中，Sn是在回波时间TEn下的信号强度，N为回波个数；ρw和ρf是水和脂肪的信号强度值；脂肪含有P个波峰分量，每个分量对应的相对幅度是αp，满足fF,p为相应的化学位移，其值已知；fB＝γ△B0局部磁场强度，其中γ为氢质子旋磁比(42.576MHz/T)，R2*为表观横向弛豫率。通常都利用多回波梯度回波序列采集多回波图像，如图1所示。为了提高采集速度，在正负读出梯度都开启模数转换器采集图像，然而在这种双极性读本文档来自技高网...

【技术保护点】
一种涡流校正方法，其特征在于，包括：步骤S1，利用双极性读出梯度回波序列采集多回波图像；步骤S2，估计采集的所述多回波图像中的涡流引入的额外相位项的一阶项系数；步骤S3，去除估计的所述一阶项系数，并估计采集的所述多回波图像中的涡流引入的额外相位项的零阶项系数；步骤S4，根据估计的所述一阶项系数和所述零阶项系数去除涡流引入的所述额外相位项的误差。

【技术特征摘要】
1.一种涡流校正方法，其特征在于，包括：步骤S1，利用双极性读出梯度回波序列采集多回波图像；步骤S2，估计采集的所述多回波图像中的涡流引入的额外相位项的一阶项系数；步骤S3，去除估计的所述一阶项系数，并估计采集的所述多回波图像中的涡流引入的额外相位项的零阶项系数；步骤S4，根据估计的所述一阶项系数和所述零阶项系数去除涡流引入的所述额外相位项的误差。2.如权利要求1所述的涡流校正方法，其特征在于，所述步骤S1中，采集的所述多回波图像至少包括6个回波图像，其中，1、3、5为奇数回波，2、4、6为偶数回波，且采集的任2个连续回波之间的时间间距相等。3.如权利要求1所述的涡流校正方法，其特征在于，所述步骤S2包括：步骤S21，在采集的所述多回波图像中寻找空间分布上组织分布均匀的像素点，所述像素点在n＝1,2,3,….下均满足以下条件：abs{Sn(x,y,z)}>t1,abs{Sn(x-1,y,z)}>t1其中，abs(.)表示取复数信号的幅值，t1，δ1，δ2为预设的阈值，(x,y,z)表示图像坐标，Sn(x,y,z)表示第n个回波的坐标位置为(x,y,z)处的信号强度；步骤S22，对寻找的所述像素点进行空间相位差分，得到其中，定义：步骤S23，结合空间相位差分得到的结果求取涡流引入的额外相位项的一阶项系数α，其中，m表示步骤S21寻找出的所有像素点，angle(.)表示对复数取相位角；其中，定义：其中，M表示上式中参与计算的DXn的个数。4.如权利要求1所述的涡流校正方法，其特征在于，所述步骤S3包括：步骤S31，根据如下公式去除涡流引入的额外相位项的所述一阶项系数α；其中，Sn是在回波时间TEn下的信号强度，N为回波个数；ρw和ρf是水和脂肪的信号强度值；脂肪含有P个波峰分量，每个分量对应的相对幅度是αp，满足fF,p为相应的化学位移，其值已知；fB＝γ△B0为局部磁场强度，其中，γ为氢质子旋磁比(42.576MHz/T)，R2*为表观横向弛豫率，β为涡流引入的额外相位项的零阶项系数。步骤S32，从去除所述一阶项系数的所有回波图像中寻找满足如下条件的信噪比较高的像素点；abs{SCn(x,y,z)}>t1,abs{SCn(x-1,y,z)}>t1,n＝1,2,3,...步骤S33，从前三个偶数回波图像寻找出的信噪比较高的像素点中筛选出满足预设水脂比例的像素点；首先，选取前三个偶数回波图像：I1＝SC2,I2＝SC4,I3＝SC6定义：其中，δ为脂肪中波谱主峰与水的化学位移差，B0表示主磁场强度，△TE表示回波时间间距；然后，仅对于寻找出的所述信噪比较高的像素点计算：最后，筛选出满足如下不等式的像素点；Rwf>c1orRwf<c2其中，Rwf表示水脂比例，c1、c2表示预设的水脂比例常数；步骤S34，对于筛选出的满足预设水脂比例的像素点，计算：并估计涡流引入的额外相位项的零阶项系数：其中，k表示步骤S33筛选出的所有像素点，angle(.)表示对复数取相位角。5.如权利要求1所述的涡流校正方法，其特征在于，所述步骤S4包括：结合估计的所述一阶项系数和所述零阶项系数并根据如下公式去除涡流引入的所述额外相位项的误差，得到校正后的数据SEn；其中，Sn是在回波时间TEn下的信号强度，N为回波个数；ρw和ρf是水和脂肪的信号强度值；脂肪含有P个波峰分量，每个分量对应的相对幅度是αp，满足fF,p为相应的化学位移，其值已知；fB＝γ△B0为局部磁场强度，其中γ为氢质子旋磁比(42.576MHz/T)，R2*为表观横向弛豫率。6.一种涡流校正装置，其特征在于，包括：采集模块，用于利用双极性读出梯度回波序列采集多回波图像；第一估计模块，用于估计采集的所述多回波图像中的涡流引入的额外相位项的一阶项系数；第二估计模块，用于去除估计的所述一阶项系数，并估计采集的所述多回波图像...

【专利技术属性】
技术研发人员：邹超，程传力，刘新，郑海荣，
申请(专利权)人：深圳先进技术研究院，
类型：发明
国别省市：广东,44