CN104035059B - 平面回波成像序列图像的重建方法 - Google Patents

Abstract

本发明提供了一种平面回波成像序列的图像重建方法，包括获取平面回波成像数据S_i，并同时采集三条没有经过相位编码的参考回波信号R₁、R₂、R₃；通过所述参考回波信号计算出需要对所述平面回波成像数据进行校正的参数；将所述平面回波成像数据沿读出方向进行一维傅里叶变换得变换结果FS_i，并用所述校正参数校正FS_i，计算出校正后的平面回波成像数据；对校正后的平面回波成像数据沿相位编码方向做一维傅里叶变换得到图像。本发明提供的平面回波成像序列的图像重建方法，在有偏离场存在的情况下，能够保持平面回波序列成像快速成像优点且有效地去除N/2伪影，同时还能够校正由于偏离场存在而导致的图像变形。

Description

平面回波成像序列图像的重建方法

【技术领域】

本发明涉及磁共振成像领域，尤其是涉及一种平面回波成像序列图像的重建方法。

【背景技术】

目前在磁共振成像过程中，平面回波成像序列是一种快速的、只需要单个激发脉冲就可获得完整图像的磁共振成像序列。磁共振梯度系统周期性变换，产生一系列梯度回波；受激平面的图像通过对产生的回波序列使用傅里叶变换而得到。

在传统的平面回波成像中，序列的相位编码梯度和频率编码梯度分别如图1中表示的a和b，由于硬件系统的不完美，比如涡电流，数据采集延迟，磁场漂移等原因，在正和负的频率编码梯度下读出的数据存在差异，如回波中心不一致，k空间每条线的间隔不等等情况。由于回波序列中奇数和偶数回波时间的不匹配或者相位差导致在重建之后得到的正常图像在偏离实际的图像视野一半的位置出现低强度的附加图像，此即为N/2伪影，图像出现了整体漂移以及变形如图5所示。

为了尽可能地消除N/2伪影，现有技术中采用了一些校正方法，主要为两种：相位校正和单边读出梯度法。

相位校正法包括一维相位校正方法以及二维相位校正法。一维相位校正法，在读取数据之前，采集几个不加相位编码的回波信号，如图1中的回波信号c，根据第一个和第二个回波信号，算出它们相位差异，把这些相位差作为校正量来校正采集到的图像数据，即校正所有偶数项信号和奇数项信号之间的相位差。由于不是所有偶数项信号与奇数项信号之间的相位差都等于求得的校正量，所以该方法虽然简单，但是效果不好，并不能有效地消除N/2伪影，同时目前的一维相位校正法并不能同时校正平面回波图像的变形，这在低磁场的磁共振系统尤为明显。

二维相位校正法，同样是在读取图像数据前，采集一序列加相位编码的回波信号，利用这些采集到的数据纠正真正的图像数据。该方法消除N/2伪影的效果比较好，但是序列采集时间延长，平面回波成像序列失去了快速成像的优点。

单边读出梯度法只在一个极性的读出梯度时采集数据，该方法可以完全避免N/2伪影问题，但是不能纠正图像变形，浪费扫描时间。而且由于回波时间间隔比较长，图像在相位编码方向上的位移，以及因磁场不均匀导致的图像变形更为严重。

综上所述，现有技术中能够有效消除N/2伪影的校正方法采集时间较长或者不能够校正由于偏离场而导致的图像变形，采集时间较短的校正方法去除N/2伪影的效果比较差且同时也不能够校正图像的变形。因此，需要提供一种更为精准的平面回波成像序列的图像重建方法。

【发明内容】

为了解决现有技术中平面回波成像序列中消除N/2伪影效果差，不能校正图像变形，以及采集时间较长的问题，本发明提供了一种更为精准的平面回波成像序列的图像重建方法。

一种平面回波成像序列的图像重建方法，包括如下步骤：

获取平面回波成像数据S_i，并同时采集三条没有经过相位编码的参考回波信号R₁、R₂、R₃，所述三条参考回波信号分别为偶信号、奇信号以及偶信号；

通过所述参考回波信号计算出需要对所述平面回波成像数据进行校正的参数；

将所述平面回波成像数据沿读出方向进行一维傅里叶变换得变换结果FS_i，并用所述校正参数校正FS_i，计算出校正后的平面回波成像数据；

对校正后的平面回波成像数据沿相位编码方向做一维傅里叶变换得到图像。

可选地，所述校正参数包括零次项相位偏差因子一次项相位偏差因子

可选地，对所述平面回波成像数据的偶数项信号或者奇数项信号进行校正，得到校正后的数据CFS_i(n)，

其中i为平面回波成像信号序数中的偶数项，或者

其中i为平面回波成像信号序数中的奇数项；

其中n为读出方向样本点的序数，Nro为读出方向总样本点的个数，当所述参考回波信号与平面回波成像信号的相应奇偶梯度极性一致时，所述sign＝+1，当所述参考回波信号与平面回波成像信号的相应奇偶梯度极性相反时，所述sign＝-1。

可选地，对所述平面回波成像数据的偶数项信号与奇数项信号进行校正，得到校正后的数据CFS_i(n)，

其中i为平面回波成像信号序数中的偶数项，以及

其中i为平面回波成像信号序数中的奇数项；

其中n为读出方向样本点的序数，Nro为读出方向总样本点的个数，当所述参考回波信号与平面回波成像信号的相应奇偶梯度极性一致时，所述sign＝+1，当所述参考回波信号与平面回波成像信号的相应奇偶梯度极性相反时，所述sign＝-1。

可选地，所述校正参数包括零次项相位偏移因子一次项相位偏移因子

可选地，对所述平面回波成像数据的偶数项信号与奇数项进行校正，得到校正后的数据CFS_i(n)，

其中i为平面回波成像信号序数；

所述n为读出方向样本点的序数，Nro为读出方向总样本点的个数，Npe为相位编码方向样本点总数。

可选地，对所述平面回波成像数据的偶数项信号与奇数项信号进行校正，得到校正后的数据CFS_i(n)，

其中i为平面回波成像信号序数中的偶数项，以及

其中i为平面回波成像信号序数中的奇数项；或者

其中i为平面回波成像信号序数中的偶数项，以及

其中i为平面回波成像信号序数中的奇数项；或者

其中i为平面回波成像信号序数中的偶数项，以及

其中i为平面回波成像信号序数中的奇数项；

其中n为读出方向样本点的序数，Nro为读出方向总样本点的个数，Npe为相位编码方向样本点总数，当所述参考回波信号与平面回波成像信号的相应奇偶梯度极性一致时，所述sign＝+1，当所述参考回波信号与平面回波成像信号的相应奇偶梯度极性相反时，所述sign＝-1。

可选地，所述一次项相位偏差因子所述一次项相位偏移因子与第一一次项相位第二一次项相位关联，所述第一一次项相位为所述参考回波信号R₁与R₂的一次项相位差，所述第二一次项相位为所述参考回波信号R₂与R₃的一次项相位差。

可选地，所述一次项相位偏差因子所述一次项相位偏移因子

可选地，所述零次项相位偏差因子所述零次项相位偏移因子与第一零次项相位第二零次项相位关联，所述第一零次项相位为所述参考回波信号R₁与R₂沿读出方向的零次项相位差，所述第二零次项相位为所述参考回波信号R₂与R₃沿读出方向的零次项相位差。

可选地，所述零次项相位偏差因子所述零次项相位偏移因子

可选地，对所述参考回波信号R₁、R₂、R₃傅里叶变换后复数共轭相乘，用来计算需要对所述平面回波成像数据进行校正的参数。

可选地，将所述参考回波信号R₁、R₂、R₃傅里叶变换得变换结果：FR₁、FR₂与FR₃，

FR₁＝FT_1D{R₁}

FR₂＝FT_1D{R₂}

FR₃＝FT_1D{R₃}

然后将所述FR₁与FR₂、FR₂与FR₃复数共轭相乘得运算结果：P₁(n)、P₂(n)，

P₁(n)＝FR₁(n)×conj[PR₂(n)]

P₂(n)＝FR₃(n)×conj[FR₂(n)]

其中n的取值范围为[1，Nro]；

分别取P₁(n)、P₂(n)的幅值最大值AMAX₁、AMAX₂；

将所述P₁(n)及P₂(n)邻近两点的复数共轭相乘得运算结果：DP₁(n)、DP₂(n)，

DP₁(n)＝P₁(n+1)×conj[P₁(n)]

DP₂(n)＝P₂(n+1)×conj[P₂(n)]

将所述DP₁(n)对应所述P₁(n)幅值大于τ×AMAX₁的所有点相加得运算结果SDP₁(n)，所述DP₂(n)对应所述P₂(n)幅值大于τ×AMAX₂的所有点相加得运算结果SDP₂(n)，其中τ为预设参数，δ为0或1，当所述P₁(n)幅值大于τ×AMAX₁或所述P₂(n)幅值大于τ×AMAX₂时，δ＝1，否则δ＝0；

分别取所述SDP₁(n)与SDP₂(n)的相位即可得所述第一一次项相位所述第二一次项相位

可选地，对所述P₁(n)及P₂(n)的相位进行校正，以消除和的影响，得到校正结果CP₁(n)、CP₂(n)，

将所述CP₁(n)对应所述P₁(n)幅值大于τ×AMAX₁的所有点相加得运算结果SCP₁(n)，所述CP₂(n)对应所述P₂(n)幅值大于τ×AMAX₂的所有点相加得运算结果SCP₂(n)，其中τ为预设参数，δ为0或1，当所述P₁(n)幅值大于τ×AMAX₁或所述P₂(n)幅值大于τ×AMAX₂时，δ＝1，否则δ＝0；

分别取所述SCP₁(n)与SCP₂(n)的相位即可得所述第一零次项相位所述第二零次项相位

可选地，所述预设参数τ＝0.5。

可选地，若使用多通道采集平面回波数据，则将所述FR₁与FR₂、FR₂与FR₃复数共轭相乘然后进行多通道合并，包括：

$P_1\left(n\right)=\underset{l=1}{\overset{\mathrm{Ncoil}}{\mathrm{\Σ}}}{\mathrm{FR}}_1^l\left(n\right)\×\mathrm{conj}\left[{\mathrm{FR}}_2^l\left(n\right)\right]$

$P_2\left(n\right)=\underset{l=1}{\overset{\mathrm{Ncoil}}{\mathrm{\Σ}}}{\mathrm{FR}}_3^l\left(n\right)\×\mathrm{conj}\left[{\mathrm{FR}}_2^l\left(n\right)\right]$

其中l与Ncoil分别为接收通道序数、接收通道的总数。

本发明提供了平面回波成像序列的图像重建方法，保持了一维相位校正法快速简单的优点，且能够在有偏离场存在的情况下有效地去除N/2伪影，同时能够校正由于偏离场存在而导致的图像变形。

【附图说明】

图1现有技术中的平面回波成像序列示意图；

图2为本发明中的平面回波成像序列图像重建方法的流程示意图；

图3为平面回波数据的相位偏差和相位偏移校正示意图；

图4为本发明中所述参考数据线的示意图；

图5为未经过校正的平面回波图像；

图6为采用本发明中技术方案采集的平面回波信号K空间数据；

图7为采用本发明中技术方案校正的平面回波图像。

【具体实施方式】

为使本发明的上述目的、特征和优点能够更为明显易懂，下面结合附图和实施例对本发明的具体实施方式做详细的说明。

在下面的描述中阐述了很多具体细节以便于充分理解本发明，但是本发明还可以采用其他不同于在此描述的方式来实施，因此本发明不受下面公开的具体实施例的限制。

如背景技术中提到的在平面回波成像序列图像重建过程中，由于读出梯度正负极性的差异，导致读出的数据并不完全一致，如回波中心不一致，K空间上每条线的间隔不等等，在图像上即表现为N/2伪影，同时又由于磁场的不均匀性即偏离场的存在，图像会出现一定的变形，因此需要提供一种平面回波成像序列的图像重建方法，以克服上述问题。

如图2所示，一种平面回波成像序列的图像重建方法，包括如下步骤：

步骤S10，获取平面回波成像数据S_i，并同时采集三条没有经过相位编码的参考回波信号R₁、R₂、R₃，所述三条回波信号分别为偶信号、奇信号以及偶信号；

步骤S11，通过所述参考回波信号计算出需要对所述平面回波成像数据进行校正的参数；

步骤S12，将所述平面回波成像数据沿读出方向进行一维傅里叶变换得变换结果FS_i，并用所述校正参数校正FS_i，计算出校正后的平面回波成像数据；

步骤S13，对校正后的平面回波成像数据沿相位编码方向做一维傅里叶变换得到图像。

下面结合实例具体说明平面回波成像序列图像的重建过程。

执行步骤S10，获取平面回波成像数据S_i，并同时采集三条没有经过相位编码的参考回波信号R₁、R₂、R₃，所述三条回波信号分别为偶信号、奇信号以及偶信号。

脉冲激发以及磁场梯度周期性变换，产生一系列梯度回波，平面回波成像数据S_i即这些梯度回波信号，同时采集三条没有经过相位编码的回波信号，这里需要说明的是，三个参考回波信号的采集是在成像数据采集的同一次扫描中进行的。

执行步骤S11，通过所述参考回波信号计算出需要对所述平面回波成像数据进行校正的参数；

所述校正参数包括零次项相位偏差因子一次项相位偏差因子零次项相位偏移因子一次项相位偏移因子

如图3所示，图中的水平直线代表了在理想情况下平面回波信号的相位随回波位置的变化。与之相交的有一定角度的直线代表了在有偏离场ΔB存在的情况下回波信号的相位随回波位置的变化。同时由于正负极性差异而散落在直线两侧的数据即表示正负回波信号的相位差异。偏离场ΔB的存在导致平面回波图像存在一定的变形，而散落在直线两侧所表示的回波信号差异表现在图像上即是N/2伪影。

为了得到清晰完整的能够用于诊断的高质量图像，必须要消除伪影以及校正图像的变形。当偏离场ΔB的主要来源为读出(以x轴表示)小梯度的切换时，它与空间位置x呈一定的线性关系。因此而产生的回波信号相位偏差及偏移也成线性变化。

用表述沿读出方向的回波信号相位，则有 其中零次项相位偏移一次项相位偏移零次项相位偏差一次项相位偏差

从上面可以看出，校正零次项相位偏差及一次项相位偏差即可消除N/2伪影带来的影响，而校正零次项相位偏移及一次项相位偏移即可校正ΔB带来的图像变形。

下面详细介绍零次项相位偏差因子一次项相位偏差因子零次项相位偏移因子一次项相位偏移因子这四个校正参数的计算的具体过程，如下：

如图4所示，分别为三条采集的参考回波信号幅值，明显可见奇偶回波在k空间的位置差异，参考回波信号R₁与R₃同为偶数项信号或者奇数项信号，参考回波信号R₂明显向右偏移。需要说明的是图中的参考的回波信号的幅值、频率等信息都是在一次实验中获取的数据，并不表示采用本技术方案的参考回波信号必须与图示中的幅值频率一致。

将所述参考回波信号R₁、R₂、R₃傅里叶变换后复数共轭相乘，用来计算需要对所述平面回波成像数据进行校正的参数。下面以R₁和R₂、R₂和R₃傅里叶变换后复数共轭相乘为例进行说明：

将所述参考回波信号R₁、R₂、R₃傅里叶变换得变化结果FR₁、FR₂以及FR₃，

FR₁＝FT_1D{R₁}

FR₂＝FT_1D{R₂}

FR₃＝FT_1D{R₃}

然后将所述FR₁与FR₂、FR₂与FR₃复数共轭相乘得运算结果P₁(n)、P₂(n)，

P₁(n)＝FR₁(n)×conj[FR₂(n)]

P₂(n)＝FR₃(n)×conj[FR₂(n)]

其中n的取值范围为[1，Nro]，n及所述Nro分别为读出方向的样本点序数、总样本点的个数。分别取P₁(n)、P₂(n)的幅值最大值AMAX₁、AMAX₂；

将所述P₁(n)及P₂(n)邻近两点的复数共轭相乘得运算结果DP₁(n)、DP₂(n)，

DP₁(n)＝P₁(n+1)×conj[P₁(n)]

DP₂(n)＝P₂(n+1)×conj[P₂(n)]

将所述DP₁(n)对应所述P₁(n)幅值大于τ×AMAX₁的所有点相加得运算结果SDP₁(n)，

${\mathrm{SDP}}_1\left(n\right)={\∑}_{n=1}^{\mathrm{Nro}-1}\mathrm{\δ}\×{\mathrm{DP}}_1\left(n\right)$

所述DFP₂(n)对应所述P₂(n)幅值大于τ×AMAX₂的所有点相加得运算结果SDP₂(n)，

${\mathrm{SDP}}_2\left(n\right)={\∑}_{n=1}^{\mathrm{Nro}-1}\mathrm{\δ}\×{\mathrm{DP}}_2\left(n\right)$

其中τ为预设参数，比如τ＝0.5，δ为0或1，当所述P₁(n)幅值大于τ×AMAX₁或所述P₂(n)幅值大于τ×AMAX₂时，δ＝1，否则δ＝0。

分别取所述SDP₁(n)与SDP₂(n)的相位即可得第一一次项相位第二一次项相位

对所述P₁(n)及P₂(n)的相位进行校正，以消除和的影响，得到校正结果CP₁(n)、CP₂(N)，

将所述CP₁(n)对应所述P₁(n)幅值大于τ×AMAX₁的所有点相加得运算结果SCP₁(n)，

${\mathrm{SCP}}_1\left(n\right)={\∑}_{n=1}^{\mathrm{Nro}-1}\mathrm{\δ}\×{\mathrm{CP}}_1\left(n\right)$

所述CP₂(n)对应所述P₂(n)幅值大于τ×AMAX₂的所有点相加得运算结果SCP₂(n)，

${\mathrm{SCP}}_2\left(n\right)={\∑}_{n=1}^{\mathrm{Nro}-1}\mathrm{\δ}\×{\mathrm{CP}}_2\left(n\right)$

其中τ为预设参数，比如τ＝0.5，δ为0或1，当所述FP₁(n)幅值大于τ×AMAX₁或所述FP₂(n)幅值大于τ×AMAX₂时，δ＝1，否则δ＝0。

分别取所述SCP₁(n)与SCP₂(n)的相位即可得第一零次项相位第二零次项相位

所述一次项相位偏差因子所述一次项相位偏移因子与所述第一一次项相位所述第二一次项相位关联，所述第一一次项相位为第一与第二参考信号沿读出方向的一次项相位差，所述第二一次项相位为第三与第二参考信号沿读出方向的一次项相位差。

则所述一次项相位偏差因子所述一次项相位偏移因子

所述零次项相位偏差因子所述零次项相位偏移因子与第一零次项相位第二零次项相位关联，所述第一零次项相位为第一(R₁)与第二(R₂)参考信号沿读出方向的零次项相位差，所述第二零次项相位为第三(R₃)与第二(R₂)参考信号沿读出方向的零次项相位差。

则所述零次项相位偏差因子所述零次项相位偏移因子

执行步骤S12，对所述平面回波成像数据沿读出方向进行一维傅里叶变换得FS_i，并用所述校正参数校正FS_i，计算出校正后的平面回波成像数据。

如果只对平面回波数据进行N/2伪影校正，只需要使用零次项相位偏差因子一次项相位偏差因子通常有两种方法：

方法一：对所述平面回波成像数据的偶数项信号或者奇数项信号进行校正：

其中i为平面回波成像信号序数中的偶数项；或者

其中i为平面回波成像信号序数中的奇数项。

方法二：对所述平面回波成像数据的偶数项信号与奇数项信号进行校正：

其中i为平面回波成像信号序数中的偶数项；以及

其中i为平面回波成像信号序数中的奇数项。

如果只对平面回波数据进行变形校正，只需要使用零次项相位偏移因子一次项相位偏移因子包括对所述平面回波成像数据的偶数项信号与奇数项进行校正：

其中i为平面回波成像信号的序数。

如果需要同时对平面回波数据进行N/2伪影以及变形校正，包括：

方法一：对所述平面回波成像数据的偶数项信号进行偏移校正，对偶数项以及奇数项信号进行偏移和偏差校正：

其中i为平面回波成像信号序数中的偶数项，以及

其中i为平面回波成像信号序数中的奇数项；

方法二：对所述平面回波成像数据的偶数项信号进行偏移和偏差校正，对奇数项信号进行偏移校正：

其中i为平面回波成像信号序数中的偶数项，以及

其中i为平面回波成像信号序数中的奇数项。

方法三：对所述平面回波成像数据的偶数项信号和奇数项信号同时进行偏移和1/2偏差校正：

其中i为平面回波成像信号序数中的偶数项，以及

其中i为平面回波成像信号序数中的奇数项；

上述所有公式中，所述n为读出方向样本点的序数，Nro为读出方向总样本点的个数，Npe为相位编码方向总样本点的个数，当所述参考回波信号与平面回波成像信号的相应奇偶梯度极性一致时，所述sign＝+1，当所述参考回波信号与平面回波成像信号的相应奇偶梯度极性相反时，所述sign＝-1。

执行步骤S13，对校正后的平面回波成像数据沿相位编码方向做一维傅里叶变换得到图像。即I＝FT_1D{CFS}，图6～7分别为采用本发明的技术方案所采集的成像K空间数据以及经过校正的图像。

在磁共振系统中，目前普遍采用多通道方式采集成像数据，因此在本技术方案中，若通过多通道方式采集平面回波数据，则需要将经过傅里叶变换的参考回波信号复数共轭相乘后进行多通道合并，具体计算方式如下：

$P_1\left(n\right)=\underset{l=1}{\overset{\mathrm{Ncoil}}{\mathrm{\Σ}}}{\mathrm{FR}}_1^l\left(n\right)\×\mathrm{conj}\left[{\mathrm{FR}}_2^l\left(n\right)\right]$

$P_2\left(n\right)=\underset{l=1}{\overset{\mathrm{Ncoil}}{\mathrm{\Σ}}}{\mathrm{FR}}_3^l\left(n\right)\×\mathrm{conj}\left[{\mathrm{FR}}_2^l\left(n\right)\right]$

其中l与Ncoil分别为接收通道序数、接收通道的总数。

需要说明的是，通过以上的实施方式的描述，本领域的技术人员可以清楚地了解到本发明的部分或全部可借助软件并结合必需的通用硬件平台来实现。基于这样的理解，本发明的技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来，该计算机软件产品可包括其上存储有机器可执行指令的一个或多个机器可读介质，这些指令在由诸如计算机、计算机网络或其他电子设备等一个或多个机器执行时可使得该一个或多个机器根据本发明的实施例来执行操作。机器可读介质可包括，但不限于，软盘、光盘、CD-ROM(紧致盘-只读存储器)、磁光盘、ROM(只读存储器)、RAM(随机存取存储器)、EPROM(可擦除可编程只读存储器)、EEPROM(电可擦除可编程只读存储器)、磁卡或光卡、闪存、或适于存储机器可执行指令的其他类型的介质/机器可读介质。

本发明可用于众多通用或专用的计算系统环境或配置中。例如：个人计算机、服务器计算机、手持设备或便携式设备、平板型设备、多处理器系统、基于微处理器的系统、置顶盒、可编程的消费电子设备、网络PC、小型计算机、大型计算机、包括以上任何系统或设备的分布式计算环境等。

本发明可以在由计算机执行的计算机可执行指令的一般上下文中描述，例如程序模块。一般地，程序模块包括执行特定任务或实现特定抽象数据类型的例程、程序、对象、组件、数据结构等等。也可以在分布式计算环境中实践本申请，在这些分布式计算环境中，由通过通信网络而被连接的远程处理设备来执行任务。在分布式计算环境中，程序模块可以位于包括存储设备在内的本地和远程计算机存储介质中。

综上所述，本发明提供的平面回波成像序列的图像重建方法，能够在保持平面回波序列成像快速成像特点的情况下，有效地去除N/2伪影，同时还能够校正由于偏离场存在而导致的图像变形。

虽然本发明已以较佳实施例揭示如上，然其并非用以限定本发明，任何本领域技术人员，在不脱离本发明的精神和范围内，当可作些许的修改和完善，因此本发明的保护范围当以权利要求书所界定的为准。

Claims (16)

1.一种平面回波成像序列的图像重建方法，其特征在于，包括如下步骤：

获取平面回波成像数据S_i，并同时采集三条没有经过相位编码的参考回波信号R₁、R₂、R₃，所述三条参考回波信号分别为偶信号、奇信号以及偶信号；

通过所述参考回波信号计算出需要对所述平面回波成像数据进行校正的参数；

将所述平面回波成像数据沿读出方向进行一维傅里叶变换得变换结果FS_i，并用所述校正参数校正FS_i，计算出校正后的平面回波成像数据；

对校正后的平面回波成像数据沿相位编码方向做一维傅里叶变换得到图像。

2.如权利要求1所述平面回波成像序列图像重建方法，其特征在于，所述校正参数包括零次项相位偏差因子一次项相位偏差因子

3.如权利要求2所述的平面回波成像序列图像重建方法，其特征在于，对所述平面回波成像数据的偶数项信号或者奇数项信号进行校正，得到校正后的数据CFS_i(n)，

其中i为平面回波成像信号序数中的偶数项，或者

其中i为平面回波成像信号序数中的奇数项；

其中n为读出方向样本点的序数，Nro为读出方向总样本点的个数，当所述参考回波信号与平面回波成像信号的相应奇偶梯度极性一致时，所述sign＝+1，当所述参考回波信号与平面回波成像信号的相应奇偶梯度极性相反时，所述sign＝-1。

4.如权利要求2所述的平面回波成像序列图像重建方法，其特征在于，对所述平面回波成像数据的偶数项信号与奇数项信号进行校正，得到校正后的数据CFS_i(n)，

其中i为平面回波成像信号序数中的偶数项，以及

其中i为平面回波成像信号序数中的奇数项；

其中n为读出方向样本点的序数，Nro为读出方向总样本点的个数，当所述参考回波信号与平面回波成像信号的相应奇偶梯度极性一致时，所述sign＝+1，当所述参考回波信号与平面回波成像信号的相应奇偶梯度极性相反时，所述sign＝-1。

5.如权利要求1或2所述的平面回波成像序列图像重建方法，其特征在于，所述校正参数包括零次项相位偏移因子一次项相位偏移因子

6.如权利要求5所述的平面回波序列图像重建方法，其特征在于，对所述平面回波成像数据的偶数项信号与奇数项信号进行校正，得到校正后的数据CFS_i(n)，

其中i为平面回波成像信号序数；

所述n为读出方向样本点的序数，Nro为读出方向总样本点的个数，Npe为相位编码方向样本点总数。

7.如权利要求5所述的平面回波成像序列图像重建方法，其特征在于，对所述平面回波成像数据的偶数项信号与奇数项信号进行校正，得到校正后的数据CFS_i(n)，

其中i为平面回波成像信号序数中的偶数项，以及

其中i为平面回波成像信号序数中的奇数项；或者

其中i为平面回波成像信号序数中的偶数项，以及

其中i为平面回波成像信号序数中的奇数项；或者

其中i为平面回波成像信号序数中的偶数项，以及

其中i为平面回波成像信号序数中的奇数项；

其中n为读出方向样本点的序数，Nro为读出方向总样本点的个数，Npe为相位编码方向样本点总数，当所述参考回波信号与平面回波成像信号的相应奇偶梯度极性一致时，所述sign＝+1，当所述参考回波信号与平面回波成像信号的相应奇偶梯度极性相反时，所述sign＝-1。

8.如权利要求5所述的平面回波成像序列图像重建方法，其特征在于，所述一次项相位偏差因子所述一次项相位偏移因子与第一一次项相位第二一次项相位关联，所述第一一次项相位为所述参考回波信号R₁与R₂的一次项相位差，所述第二一次项相位为所述参考回波信号R₂与R₃的一次项相位差。

9.如权利要求8所述的平面回波成像序列图像重建方法，其特征在于，所述一次项相位偏差因子所述一次项相位偏移因子

10.如权利要求5所述的平面回波成像序列图像重建方法，其特征在于，所述零次项相位偏差因子所述零次项相位偏移因子与第一零次项相位第二零次项相位关联，所述第一零次项相位为所述参考回波信号R₁与R₂的零次项相位差，所述第二零次项相位为所述参考回波信号R₂与R₃的零次项相位差。

11.如权利要求10所述的平面回波成像序列图像重建方法，其特征在于，所述零次项相位偏差因子所述零次项相位偏移因子

12.如权利要求8所述的平面回波成像序列图像重建方法，其特征在于，

对所述参考回波信号R₁、R₂、R₃傅里叶变换后复数共轭相乘，用来计算需要对所述平面回波成像数据进行校正的参数。

13.如权利要求12所述的平面回波成像序列图像重建方法，其特征在于，

将所述参考回波信号R₁、R₂、R₃傅里叶变换得变换结果FR₁、FR₂与FR₃，

FR₁＝FT_1D{R₁}

FR₂＝FT_1D{R₂}

FR₃＝FT_1D{R₃}

然后将所述FR₁与FR₂、FR₂与FR₃复数共轭相乘得运算结果P₁(n)、P₂(n)，

P₁(n)＝FR₁(n)×conj[FR₂(n)]

P₂(n)＝FR₃(n)×conj[FR₂(n)]

其中n的取值范围为[1，Nro]；

分别取P₁(n)、P₂(n)的幅值最大值AMAX₁、AMAX₂；

将所述P₁(n)及P₂(n)邻近两点的复数共轭相乘得运算结果DP₁(n)、DP₂(n)，

DP₁(n)＝P₁(n+1)×conj[P₁(n)]

DP₂(n)＝P₂(n+1)×conj[P₂(n)]

将所述DP₁(n)对应所述P₁(n)幅值大于τ×AMAX₁的所有点相加得运算结果SDP₁(n)，所述DP₂(n)对应所述P₂(n)幅值大于τ×AMAX₂的所有点相加得运算结果SDP₂(n)，其中τ为预设参数，δ为0或1，当所述P₁(n)幅值大于τ×AMAX₁或所述P₂(n)幅值大于τ×AMAX₂时，δ＝1，否则δ＝0；

分别取所述SDP₁(n)与SDP₂(n)的相位即可得所述第一一次项相位所述第二一次项相位

14.如权利要求13所述的平面回波成像序列图像重建方法，其特征在于，

对所述P₁(n)及P₂(n)的相位进行校正，以消除和的影响，得到校正结果CP₁(n)、CP₂(n)，

将所述CP₁(n)对应所述P₁(n)幅值大于τ×AMAX₁的所有点相加得运算结果SCP₁(n)，所述CP₂(n)对应所述P₂(n)幅值大于τ×AMAX₂的所有点相加得运算结果SCP₂(n)，其中τ为预设参数，δ为0或1，当所述P₁(n)幅值大于τ×AMAX₁或所述P₂(n)幅值大于τ×AMAX₂时，δ＝1，否则δ＝0；

分别取所述SCP₁(n)与SCP₂(n)的相位即可得所述第一零次项相位所述第二零次项相位

15.如权要求13或14中任一项所述的平面回波成像序列图像重建方法，其特征在于，所述预设参数τ＝0.5。

16.如权利要求13所述的平面回波成像序列图像重建方法，其特征在于，

若使用多通道采集平面回波数据，则将所述FR₁与FR₂、FR₂与FR₃复数共轭相乘然后进行多通道合并，包括：

$P_1\left(n\right)=\underset{l=1}{\overset{\mathrm{Ncoil}}{\mathrm{\Σ}}}{\mathrm{FR}}_1^l\left(n\right)\×\mathrm{conj}\left[{\mathrm{FR}}_2^l\left(n\right)\right]$

$P_2\left(n\right)=\underset{l=1}{\overset{\mathrm{Ncoil}}{\mathrm{\Σ}}}{\mathrm{FR}}_3^l\left(n\right)\×\mathrm{conj}\left[{\mathrm{FR}}_2^l\left(n\right)\right]$

其中l与Ncoil分别为接收通道序数、接收通道的总数。