CN109752683A - 基于单扫描混合时空编码磁共振成像的图像畸变校正方法 - Google Patents

Abstract

本发明涉及一种基于单扫描混合时空编码磁共振成像的图像畸变校正方法，其特征在于包括：步骤1、首先对成像物体进行感兴趣区域定位；步骤2、在核磁共振成像仪上导入事先编译好的单扫描混合时空编码序列；步骤3、进行数据采样；步骤4、确定每个采样回波有效的峰值中心位置；步骤5、计算不均匀场下每个回波峰值中心位置相对于理想情况下回波中心位置的偏离点数；步骤6、计算不均匀场下每个回波峰值位置偏移引起的读梯度偏移面积；步骤7、计算每个回波信号的有效回波时间；步骤8、计算等效不均匀场磁化率梯度值G_inh；步骤9、利用步骤8计算得到的等效不均匀场梯度值对畸变的图像进行校正。本发明具备对传统相位编码方向不均匀场的鲁棒性。

Description

基于单扫描混合时空编码磁共振成像的图像畸变校正方法

技术领域

本发明涉及一种基于单扫描混合时空编码磁共振成像的图像畸变校正方法。

背景技术

人类大脑计划(Human Brain Project,HBP)是继人类基因组计划之后，又一国际性科研大计划。磁共振成像在进行大脑解密研究中发挥着重要的作用。不仅可以实现对大脑解剖结构的分区，并且还可以进行功能成像对大脑的高级认知功能进行解密。因此磁共振功能成像被公认为当今最有效的解析大脑功能的影像手段。超快速成像技术在功能成像中发挥着重要的作用，如对脑组织的血流灌注进行定量测量的动脉自旋标记(ArterialSpin Labeling,ASL)技术、揭示大脑神经束走向的弥散张量成像(Diffusion TensorImaging,DTI)、研究大脑分区活动的功能成像(functional MRI,fMRI)。在众多的超快速方法中，单扫描的平面回波成像(echo planar imaging,EPI)凭借它高的时间分辨率成为最受欢迎的超快速成像技术。但是EPI成像方法对B₀非均匀场非常敏感，很容易出现图像畸变和伪影。

近来，以色列威兹曼研究所的Frydman教授小组提出了一种单扫描的混合时空编码 (Spatiotemporally-encoded，SPEN)方法，这种方法与平面回波成像的成像原理有着本质的不同。对于EPI成像来说，每一个采样时刻获得的信号均源自整个激发样品没质子自旋的贡献。而SPEN方法，某一时刻的采样信号主要源自对应解码空间内质子自旋的贡献。这种超快速的空间编码序列只是将EPI序列中的相位编码梯度更换成时空编解码组合，且保留了EPI序列中的正负切换读出梯度，因此SPEN序列的时间分辨率与EPI 序列几乎相同。由于SPEN在传统的相位编码方向采用了时空编码的方式对该方向的质子自旋进行编解码，因此可以在采样阶段使用较大的梯度进行信号采集，从而相对的减弱了不均匀场的影响。对于SPEN成像，该种抑制非均匀场的能力主要依赖于所用的编解码梯度，编解码梯度的强度越大，SPEN抑制非均匀场的能力就越强。因此，在传统的相位编码方向，SPEN可以根据实际情况调整编解码参数，进而使非均匀场的影响降低。但是对于频率编码方向，SPEN仍然采用傅里叶编码规则进行采样，因此所受到的非均匀场的影响与EPI一样。为了加速采样效率，正负切换的回波链长度一般较大，从而使每个采样回波所受到不均匀场的影响时间不一致，从而图像出现楔形畸变，从而影响图像质量。

为了提高图像质量，基于EPI和SPEN方法的畸变校正先后被国际上的学者提出。但是通过总结以往学者的研究发现，不管是基于EPI超快速成像的畸变校正方法，还是基于SPEN超快速成像的畸变校正方法都存在一些限制。例如，基于EPI方法的畸变校正大多数需要额外的B₀场图扫描，这样就无法适用于一些MR导航的介入手术和成像部位不断运动的情况；基于SPEN超快速成像的畸变校正主要是在快速采样的连续激发和重聚序列(Rapidacquisition by sequential excitation and refocusing，RASER)和正交时空编码成像(biaxial spatiotemporally-encoded,bi-SPEN)序列的变种序列基础进行的。这两种序列的特点是含有一个90°激发脉冲和两个180°的重聚脉冲，这种序列特点就是较长的回波时间和较高的脉冲特异吸收率(Specific Absorption Rate，SAR)。较长的回波时间会引起图像的T2模糊现象和较低的信噪比，而较高的SAR值会引起成像部位的发热，尤其是在DTI和fMRI的扫描中表现尤为突出。

由此可见，开发一种无需额外B₀场图扫描并且序列本身能保证较高的信噪比和较低的SAR值的校正方法是非常有必要的。

发明内容

本发明所要解决的技术问题是针对上述现有技术提供一种基于单扫描混合时空编码序列的畸变校正方法，该方法不需要额外的扫描就能对畸变图像进行校正，并且对于极大的全局不均匀场仍有良好的表现，同时序列本身扫描具备较高的信噪比及较低SAR值。

本发明解决上述技术问题所采用的技术方案为：一种基于单扫描混合时空编码磁共振成像的图像畸变校正方法，其特征在于包括以下步骤：

步骤1、首先对成像物体进行感兴趣区域定位，然后进行自动匀场、功率和频率校正；

步骤2、在核磁共振成像仪上，导入事先编译好的单扫描混合时空编码序列；

所述事先编译好的单扫描混合时空编码序列包括合时空编码模块、第一回波延时模块、180度层选重聚模块、第二回波延时模块、采样回波链；

所述混合时空编码模块是由90度的线性扫描脉冲与作用在y方向的梯度G_exc组成，用于完成对y方向的空间编码；

其中所述第一回波延时模块的作用范围是从线性扫描脉冲的中间点至180度层选重聚模块的中间点；

所述的180度层选重聚模块是由180度的sinc脉冲与作用在z方向的层选梯度G_ss组成，用于完成空间层选及信号重聚，且在180度的sinc脉冲前后都有x,y,z三个方向的破坏梯度作用；

所述第二回波延时模块的作用范围是从所述180度的sinc脉冲中间点到所述采样回波链的中间点；

第一回波延时模块的作用时间与第二回波延时模块的作用时间之和等于整个所述采样回波链的采样时间；

所述采样回波链由分别作用在x方向的重聚梯度G_ro和作用在y方向的梯度回波链组成；

步骤3、执行步骤2设置好的单扫描混合时空编码序列，进行数据采样，数据采样完成后，执行下一步骤，否则继续采样；从而得到混合时空编码磁共振数据；

步骤4、当数据全部采样完成后，首先采用能量谱分析方法确定每个采样回波有效的峰值中心位置，将第n个采样回波的峰值中心位置标记为

所述能量谱分析方法是指对每个回波信号按照以下公式逐点进行能量计算：

其中M_ro表示每个回波的采样点数，st_ro(n,p)为第n个回波信号的离散表达式，power(n,k)表示对第n个回波信号st_ro(n,p)逐点删除后剩下信号的能量和；当采样点数减少到时，回波信号的能量急速下降，此时就认为是该回波信号的峰值中心位置点；

步骤5、计算不均匀场下每个回波峰值中心位置相对于理想情况下回波中心位置的偏离点数：

其中为均匀场下每个回波的峰值中心位置；由于在均匀场下，每个采样回波没有受到不均匀场的影响，因此每个回波的峰值中心点位置都是相同的，都为M_ro/2；

步骤6、计算不均匀场下每个回波峰值位置偏移引起的读梯度偏移面积：

Area(n)＝△k_n×G_ro

其中G_ro为读梯度的强度，Area(n)表示第n个回波中心位置偏移引起的读梯度偏移面积；该梯度偏差是由于等效的不均匀场磁化率梯度引起的，因此每个读梯度的偏移面积Area(n)进一步表示为：

其中G_inh为等效不均匀场磁化率梯度的强度，为第n个回波采集到的磁化矢量所经历的有效的驰豫时间；

步骤7、计算每个回波信号的有效回波时间：

其中T_pre为180度sinc重聚选层脉冲中心至第一个采样回波梯度之间的延迟时间，T_esp为每个采样回波的时间间隔，dt为读梯度采样时间间隔，N_pe为时空编码解码步数， TE为理想回波时间，为均匀场下线性扫频脉冲中心至第N_pe/2个采样回波峰值中心点的持续时间，T_exc为线性扫频脉冲的持续时间；

步骤8、计算等效不均匀场磁化率梯度值G_inh：

步骤9、利用步骤8计算得到的等效不均匀场梯度值对畸变的图像进行校正：

先对x方向的每个回波信号进行畸变校正，然后再沿着y方向进行时空编码重建；

在对沿x方向的每个回波信号进行畸变校正时，参照以下公式：

其中s(t_x,n)为第n个回波采集到的磁共振信号；t_x为第n个回波的采样时刻；L_x为读方向的视野大小；ρ(x)为x方向畸变校正后的图像数据；γ为氢核的旋磁比；T_ro为读梯度的持续时间，T_ro＝M_ro×dt，G_inh为等效的不均匀场梯度强度，G_ro为读梯度的强度；为第n个回波采集到的磁化矢量所经历的有效的驰豫时间。

将以上操作步骤依次实施后就可以获得高分辨率的、能正确反映空间信息的混合时空编码磁共振图像。

作为改进，所述步骤4中，首先将采集到的混合时空编码磁共振数据依照采样参数转换成二维矩阵显示，再沿着x方向采用能量谱分析方法获取每个回波实际的峰值中心位置。

与现有技术相比，本发明的优点在于：本专利首先提出一种单扫描的混合时空编码方法，该方法不仅能通过调整时空编解码参数使序列本身具备对传统相位编码方向不均匀场的鲁棒性，并且通过能量谱分析法获取每个频率编码方向采样回波的实际峰值中心位置，进而获得到该方向等效的不均匀场磁化率梯度，然后将该不均匀场磁化率梯度代入到混合时空编码的重建中，最后得到高分辨的、无楔形畸变的高品质磁共振图像；该方法能承受更大的全局和局部不均匀场的影响，同时不需要额外的场图扫描，对所采数据回波进行能量谱分析就可以获得每个回波的偏移量，进而得到等效的不均匀场磁化率梯度，相比平面回波的畸变校正方法，本发明的方法效率更高、稳定性更好，而相比基于RASER和bi-SPEN变种序列的畸变校正方法，本发明的方法序列实施简单、信噪比较高、SAR值较低，能更好的适用于高维的DTI以及长时间累计扫描的fMRI中。

附图说明

图1为本发明实施例中事先编译好的单扫描混合时空编码序列图。

图2为本发明实施例中均匀场及不均匀场下回波中心位置示意图；

图3为本发明实施例中不同序列、不同场条件下的大鼠脑部磁共振像:

A为经过图像畸变校正后的单扫描混合时空编码图像；

B为均匀场下获得的用来参考的单扫描混合时空编码图像；

C畸变的单扫描混合时空编码图像；

D畸变的平面回波图像；

具体实施方式

以下结合附图实施例对本发明作进一步详细描述。

本实施例提供的基于单扫描混合时空编码磁共振成像的图像畸变校正方法，具体包括如下步骤：

步骤1、首先对成像物体进行感兴趣区域定位，然后进行自动匀场、功率和频率校正，自动匀场、功率和频率校正的操作方法为核磁共振成像仪的常规操作方法；

步骤2、在核磁共振成像仪上，导入事先编译好的单扫描混合时空编码序列；并设置所述单扫描混合时空编码序列的各个参数；

参考图1，所述导所述事先编译好的单扫描混合时空编码序列包括时空编码模块、第一回波延时模块、180度层选重聚模块、第二回波延时模块、采样回波链；

所述混合时空编码模块是由90度的线性扫描脉冲，简称90ochirp脉冲，与作用在y方向的空间编码梯度G_exc组成，用于完成对y方向的空间编码；混合时空编码模块的作用时间为T_exc；

其中所述第一回波延时模块TE/2的作用时间范围是从线性扫描脉冲的中间点至180度层选重聚模块的中间点；

所述的180度层选重聚模块是由180度的sinc脉冲，简称180o sinc与作用在z方向的层选梯度G_ss组成，用于完成空间层选及信号重聚，且在180度的sinc脉冲前后都有x,y,z三个方向的破坏梯度作用；

所述第二回波延时模块TE/2的作用时间范围是从所述180度的sinc脉冲中间点到所述采样回波链的中间点；

第一回波延时模块的作用时间TE/2与第二回波延时模块的作用时间TE/2之和等于整个所述采样回波链的采样时间；

所述采样回波链由分别作用在x方向的重聚梯度G_ro和作用在y方向的梯度回波链组成；作用在x方向的重聚梯度G_ro包括一系列正负切换的梯度构成，且每个梯度的面积是由频率编码点数、频率编码梯度大小以及采样间隔决定的；y方向的梯度回波链是由一系列大小相等的方向空间解码梯度G_acq构成，且所述方向空间解码梯度G_acq梯度的总面积和等于所述空间编码梯度G_exc的面积，正负和所述空间编码梯度G_exc一致；在所述采样回波链之前，x方向施加了预先重聚梯度G_ror，所述预先重聚梯度G_ror的面积是x 方向第一个梯度面积的一半，方向与之相反；

设置所述单扫描混合时空编码序列的各个参数的步骤为：首先打开所述单扫描混合时空编码序列的各个相关模块，包括括时空编码模块、第一回波延时模块、180度层选重聚模块、第二回波延时模块、采样回波链；接着根据感兴趣区域大小设置进行相应的实验参数，包括90度的线性扫描脉冲的激发频率宽度△Ο₉₀，激发时间T_exc，功率tpwr₉₀， 180度的sinc脉冲的激发频率宽度△Ο₁₈₀，激发时间T₁₈₀，功率tpwr₁₈₀，x方向采样点数 M_ro，y方向采样点数N_pe，x方向的采样带宽sw，y方向空间解码梯度G_acq的方向sign，第一回波延时模块的延时时间，第一回波延时模块的延时时间，x方向的成像视野L_x， y方向的成像视野L_y，成像层厚度thk，破坏梯度的强度和时间；

步骤3、执行步骤2设置好的单扫描混合时空编码序列，进行数据采样，数据采样完成后，执行下一步骤，否则继续采样；从而得到混合时空编码磁共振数据；

有别于常规多扫描成像序列，本方法可以一次激发获得整幅磁共振像，从而可以避免运动伪影；区别于常规的单扫描序列的是本方法可以进行二维的小视野成像，而不会发生折叠伪影，可以提高图像的分辨率或者进一步加快成像速度。并且相比常规的单扫描成像方法而言，本方法在不均匀场上有更好的表现。如果要进行多层成像或者功能成像，首先序列要延迟一段时间，目的是为了让磁化率矢量弛豫恢复过来；接着，脉冲序列的各个部分依次对成像物体进行作用演化，即时空编码、层选、延时作用直至数据采样结束，最后得到二维的混合时空编码磁共振像数据。获得一幅磁共振像的时间大约几十毫秒到几百毫秒；

步骤4、当数据全部采样完成后，首先采用能量谱分析方法确定每个采样回波有效的峰值中心位置，将第n个采样回波的峰值中心位置标记为

所述能量谱分析方法是指对每个回波信号按照以下公式逐点进行能量计算：

其中M_ro表示每个回波的采样点数，st_ro(n,p)为第n个回波信号的离散表达式，power(n,k)表示对第n个回波信号st_ro(n,p)逐点删除后剩下信号的能量和；当采样点数减少到时，回波信号的能量急速下降，此时就认为是该回波信号的峰值中心位置点；

步骤5、计算不均匀场下每个回波峰值中心位置相对于理想情况下回波中心位置的偏离点数，参见图2所示；

其中为均匀场下每个回波的峰值中心位置；由于在均匀场下，每个采样回波没有受到不均匀场的影响，因此每个回波的峰值中心点位置都是相同的，都为M_ro/2；

步骤6、计算不均匀场下每个回波峰值位置偏移引起的读梯度偏移面积：

Area(n)＝△k_n×G_ro

其中G_ro为读梯度的强度，Area(n)表示第n个回波中心位置偏移引起的读梯度偏移面积；该梯度偏差是由于等效的不均匀场磁化率梯度引起的，因此每个读梯度的偏移面积Area(n)进一步表示为：

其中G_inh为等效不均匀场磁化率梯度的强度，为第n个回波采集到的磁化矢量所经历的有效的驰豫时间；

步骤7、计算每个回波信号的有效回波时间：

其中T_pre为180度sinc重聚选层脉冲中心至第一个采样回波梯度之间的延迟时间，T_esp为每个采样回波的时间间隔，dt为读梯度采样时间间隔，N_pe为时空编码解码步数， TE为理想回波时间，为均匀场下线性扫频脉冲中心至第N_pe/2个采样回波峰值中心点的持续时间，T_exc为线性扫频脉冲的持续时间；

步骤8、计算等效不均匀场磁化率梯度值G_inh：

步骤9、利用步骤8计算得到的等效不均匀场梯度值对畸变的图像进行校正：

先对x方向的每个回波信号进行畸变校正，然后再沿着y方向进行时空编码重建，在对沿x方向的每个回波信号进行畸变校正时，参照以下公式：

其中s(t_x,n)为第n个回波采集到的磁共振信号；t_x为第n个回波的采样时刻；L_x为读方向的视野大小；ρ(x)为x方向畸变校正后的图像数据；γ为氢核的旋磁比；T_ro为读梯度的持续时间，T_ro＝M_ro×dt，G_inh为等效的不均匀场梯度强度，G_ro为读梯度的强度；为第n个回波采集到的磁化矢量所经历的有效的驰豫时间；

对y方向进行时空编码重建的方法采用现有技术中的常规方法，此处不再详述。

将以上操作步骤依次实施后就可以获得高分辨率的、能正确反映空间信息的混合时空编码磁共振图像。

实施例：

我们将单扫描混合时空编码方法和平面回波成像方法在活体Wistar大鼠上进行实施例展示，用来验证本发明的可行性。实验测试是在一台小动物7T成像仪(AgilentTechnologies,Santa Clara,CA,USA)下进行。实验所采用的成像样本是约200克的活体Wistar大鼠，实验时采用异氟烷与氧气混合气体对大鼠进行全程气体麻醉，等老鼠进入休眠状态后进行相应的实验。进行实验之前，先把大鼠固定在实验床上，然后导入磁共振成像仪中。在磁共振成像仪操作台上，打开成像仪相应的操作软件，对感兴趣的大鼠部位进行定位，此处我们选择对大鼠的脑部进行coronal面成像。等成像区域定位好后，进行自动匀场、频率和功率校正。虽然进行了自动匀场，但是此时的不均匀场仍然很强，为了能体现本方法固有抵抗不均匀场的能力和校正方法的有效性，在此场环境下进行单扫描的混合时空编码成像实验，同时也进行平面回波成像实验作为对比；按照上述单扫描混合时空编码磁共振成像方法和畸变校正方法的操作流程，在进行单扫描的混合时空编码之前，首先用一维的时空编码序列分别测量所用到的90度线性扫频脉冲和180度 sinc脉冲的功率。然后导入编译好的单扫描混合时空编码序列，打开脉冲序列的各个相关模块，包括时空编码模块、回波延时模块TE、采样回波链模块，设置实验参数，具体对于本实施例所采用的样本，其试验参数设置如下：90度线性扫频脉冲的激发频率宽度△Ο90为64kHz，激发时间T_exc为3毫秒，x方向采样点数M_ro为64，y方向采样点数N_pe为64，x方向的采样带宽sw为250kHz，回波延时时间TE为31.97毫秒，x方向的成像视野L_x为4.5厘米，y方向的成像视野L_y为4.5厘米，成像层厚度thk为2毫米，破坏梯度的强度为3.0高斯每厘米，时间为1毫秒。将以上实验参数设置好后，直接运行整个序列的采样时间约为60毫秒。在此不均匀场环境下，获得单扫描混合时空编码数据经过高分辨重建后，得到高分辨图像如图3C所示；在相同场环境下，获得的同层的平面回波图像如图3D所示；在不均匀场下，获得单扫描混合时空编码和平面回波数据后，进行复杂繁琐的自动和手动匀场，使场能更加均匀，然后在均匀场下进行单扫描混合时空编码实验，获得良好磁共振像作为校正参考图像如图3B所示；

数据采样完后，按照上述步骤5至9的处理过程对在不均匀场下获得的单扫描混合时空编码数据进行畸变校正，经过处理后可以获得本实施例的高分辨无畸变的高品质磁共振像，最终结果如图3A所示。从图3可以看出，在相同的不均匀场下，平面回波图像，如图3D所示，畸变更为严重，并且出现折叠现象。经过畸变校正的混合时空编码图像，如图3A所示，与在均匀场下获得的图像，如图3B所示，有着相似的图像质量。由图3可以证明，单扫描的混合时空编码成像方法相比平面回波成像方法，在不均匀场下有着更好的表现，经过基于混合时空编码成像方法的畸变校正后，图像质量有着很大的提升。

综上所述，本发明所提出的基于单扫描混合空编码磁共振成像的图像畸变校正方法能够承受更大的全局和局部不均匀场的影响，并且采用能量谱分析方法获得每个回波的中心偏移量，进而得到等效的不均匀场磁化率梯度，再将该不均匀场磁化率梯度代入到混合时空编码的高分辨重建中，最后就可以获得高分辨的、无畸变的磁共振像。

Claims (2)

1.一种基于单扫描混合时空编码磁共振成像的图像畸变校正方法，其特征在于包括以下步骤：

步骤1、首先对成像物体进行感兴趣区域定位，然后进行自动匀场、功率和频率校正；

步骤2、在核磁共振成像仪上，导入事先编译好的单扫描混合时空编码序列；

所述事先编译好的单扫描混合时空编码序列包括时空编码模块、第一回波延时模块、180度层选重聚模块、第二回波延时模块、采样回波链；

所述混合时空编码模块是由90度的线性扫描脉冲与作用在y方向的梯度G_exc组成，用于完成对y方向的空间编码；

其中所述第一回波延时模块的作用范围是从线性扫描脉冲的中间点至180度层选重聚模块的中间点；

所述的180度层选重聚模块是由180度的sinc脉冲与作用在z方向的层选梯度G_ss组成，用于完成空间层选及信号重聚，且在180度的sinc脉冲前后都有x,y,z三个方向的破坏梯度作用；

所述第二回波延时模块的作用范围是从所述180度的sinc脉冲中间点到所述采样回波链的中间点；

第一回波延时模块的作用时间与第二回波延时模块的作用时间之和等于整个所述采样回波链的采样时间；

所述采样回波链由分别作用在x方向的重聚梯度G_ro和作用在y方向的梯度回波链组成；

步骤3、执行步骤2设置好的单扫描混合时空编码序列，进行数据采样，数据采样完成后，执行下一步骤，否则继续采样；从而得到混合时空编码磁共振数据；

步骤4、当数据全部采样完成后，首先采用能量谱分析方法确定每个采样回波有效的峰值中心位置，将第n个采样回波的峰值中心位置标记为

所述能量谱分析方法是指对每个回波信号按照以下公式逐点进行能量计算：

其中M_ro表示每个回波的采样点数，st_ro(n,p)为第n个回波信号的离散表达式，power(n,k)表示对第n个回波信号st_ro(n,p)逐点删除后剩下信号的能量和；当采样点数减少到时，回波信号的能量急速下降，此时就认为是该回波信号的峰值中心位置点；

步骤5、计算不均匀场下每个回波峰值中心位置相对于理想情况下回波中心位置的偏离点数：

其中为均匀场下每个回波的峰值中心位置；由于在均匀场下，每个采样回波没有受到不均匀场的影响，因此每个回波的峰值中心点位置都是相同的，都为M_ro/2；

步骤6、计算不均匀场下每个回波峰值位置偏移引起的读梯度偏移面积：

Area(n)＝△k_n×G_ro

其中G_ro为读梯度的强度，Area(n)表示第n个回波中心位置偏移引起的读梯度偏移面积；该梯度偏差是由于等效的不均匀场磁化率梯度引起的，因此每个读梯度的偏移面积Area(n)进一步表示为：

其中G_inh为等效不均匀场磁化率梯度的强度，为第n个回波采集到的磁化矢量所经历的有效的驰豫时间；

步骤7、计算每个回波信号的有效回波时间：

其中T_pre为180度sinc重聚选层脉冲中心至第一个采样回波梯度之间的延迟时间，T_esp为每个采样回波的时间间隔，dt为读梯度采样时间间隔，N_pe为时空编码解码步数，TE为理想回波时间，为均匀场下线性扫频脉冲中心至第N_pe/2个采样回波峰值中心点的持续时间，T_exc为线性扫频脉冲的持续时间；

步骤8、计算等效不均匀场磁化率梯度值G_inh：

步骤9、利用步骤8计算得到的等效不均匀场梯度值对畸变的图像进行校正：

先对x方向的每个回波信号进行畸变校正，然后再沿着y方向进行时空编码重建，在对沿x方向的每个回波信号进行畸变校正时，参照以下公式：

其中s(t_x,n)为第n个回波采集到的磁共振信号；L_x为读方向的视野大小；ρ(x)为x方向畸变校正后的图像数据；γ为氢核的旋磁比；T_ro为读梯度的持续时间，T_ro＝M_ro×dt，G_inh为等效不均匀场磁化率梯度值，G_ro为读梯度的强度；为第n个回波采集到的磁化矢量所经历的有效的驰豫时间；t_x为第n个回波的采样时刻，i表示相位。

2.根据权利要求1所述的基于单扫描混合时空编码磁共振成像的图像畸变校正方法，其特征在于：所述步骤4中，首先将采集到的混合时空编码磁共振数据依照采样参数转换成二维矩阵显示，再沿着x方向采用能量谱分析方法获取每个回波实际的峰值中心位置。