涡流校正方法、磁共振图像校正方法、装置、设备和介质与流程

1.本技术涉及数据处理技术领域，特别是涉及一种涡流校正方法、磁共振图像校正方法、装置、计算机设备和存储介质。


背景技术：

2.在磁共振脉冲成像序列运行过程中，会使用到一系列随时间变化的梯度场，这种随时间变化的梯度场会引起梯度线圈周边空间中的磁通量发生时变,与此同时，梯度周边的部件中的金属导电材料中会产生时变的电流，这种感生的电流被称作涡流。涡流对后期磁共振成像影响较为严重，因此必须尽可能消除涡流对后期图像的影响，而对涡流进行校正就是消除对图像影响的一种手段。
3.相关技术中，在对涡流进行校正时，通常是通过预加重的方式来实现。也就是说，在前期对磁共振系统的硬件(例如梯度线圈等)加上一个反向的梯度，然后通过该反向的梯度实现对涡流的校正。
4.然而，上述技术在对涡流进行校正后，仍然会存在涡流残留，即该方式对涡流校正的准确性较低。


技术实现要素：

5.基于此，有必要针对上述技术问题，提供一种能够提高对磁共振涡流校正的准确性的涡流校正方法、磁共振图像校正方法、装置、计算机设备和存储介质。
6.一种涡流校正方法，该方法包括：
7.获取第一梯度施加后采集的多个回波时间的第一磁共振信号和第二梯度施加后采集的多个回波时间的第二磁共振信号；上述多个回波时间的第一磁共振信号的采集时刻与上述多个回波时间的第二磁共振信号的采集时刻相对应；
8.根据各个回波时间的第一磁共振信号和第二磁共振信号，确定各个回波时间的两个磁共振图像，并确定各个回波时间的两个磁共振图像之间的相位差；
9.根据各相位差，确定至少一个成像回波时间下的涡流场分布；
10.利用成像回波时间下的涡流场分布对所述成像回波时间对应的磁共振信号进行涡流校正。
11.在其中一个实施例中，上述根据各相位差，确定至少一个成像回波时间下的涡流场分布，包括：
12.对各相位差进行处理，确定每个回波时间下的涡流场分布；
13.对每个回波时间下的涡流场分布进行求解，确定至少一个成像回波时间下的涡流场分布。
14.在其中一个实施例中，上述对每个回波时间下的涡流场分布进行求解，确定至少一个成像回波时间下的涡流场分布，包括：
15.根据预设的球谐函数对每个回波时间下的涡流场分布进行系数展开，确定各个成
像回波时间下的涡流球谐系数；
16.根据上述各个成像回波时间下的涡流球谐系数，得到至少一个成像回波时间下的涡流场分布。
17.在其中一个实施例中，上述获取第一梯度施加后采集的多个回波时间的第一磁共振信号和第二梯度施加后采集的多个回波时间的第二磁共振信号，包括：
18.在上述第一梯度施加后，采用预扫描序列进行多层数据交叉采集，获得上述第一梯度施加后多个回波时间的第一磁共振信号；
19.在上述第二梯度施加后，采用预扫描序列进行多层数据交叉采集，获得上述第二梯度施加后多个回波时间的第二磁共振信号。
20.在其中一个实施例中，上述第一梯度和上述第二梯度为梯度强度大小相等，且方向相反的两个梯度；上述根据各相位差，确定各个成像回波时间下的涡流场分布，包括：
21.根据各相位差，确定上述第一梯度施加后或第二梯度施加后各个成像回波时间下的涡流场分布。
22.在其中一个实施例中，上述方法还包括：
23.获取预设的第三梯度的强度和第一梯度的强度之间的第一强度比例，或第三梯度的强度和上述第二梯度的强度之间的第二强度比例；
24.根据上述第一强度比例或第二强度比例，对上述第一梯度施加后和第二梯度施加后各个成像回波时间下的涡流场分布进行线性叠加处理，确定上述第三梯度下各个成像回波时间下的涡流场分布。
25.在其中一个实施例中，上述根据各个回波时间的第一磁共振信号和第二磁共振信号，确定各个回波时间的两个磁共振图像，并确定各个回波时间的两个磁共振图像之间的相位差，包括：
26.对各个回波时间的第一磁共振信号和第二磁共振信号进行图像重建，获得各个回波时间的第一磁共振图像和第二磁共振图像；
27.根据各个回波时间的第一磁共振图像获得各个回波时间的第一相位，以及根据各个回波时间的第二磁共振图像获得各个回波时间的第二相位；
28.对各个回波时间的第一相位和第二相位进行求差运算，得到各个回波时间的第一磁共振图像和第二磁共振图像之间的相位差。
29.一种磁共振图像校正方法，该方法包括：
30.获取扫描对象的多组磁共振图像，每组磁共振图像包括对应同一回波时间的第一磁共振图像和第二磁共振图像，上述第一磁共振图像为第一梯度施加后采集的第一磁共振信号重建获得，上述第二磁共振图像为第二梯度施加后采集的第二磁共振信号重建获得；
31.确定每组磁共振图像中第一磁共振图像和第二磁共振图像之间的相位差；
32.根据上述相位差，确定至少一个成像回波时间的涡流场分布；
33.根据上述涡流场分布对所述成像回波时间对应的待处理磁共振图像进行校正，获取校正的磁共振图像。
34.一种涡流校正装置，该装置包括：
35.信号获取模块，用于获取第一梯度施加后采集的多个回波时间的第一磁共振信号和第二梯度施加后采集的多个回波时间的第二磁共振信号；上述多个回波时间的第一磁共
振信号的采集时刻与上述多个回波时间的第二磁共振信号的采集时刻相对应；
36.第一相位差确定模块，用于根据各个回波时间的第一磁共振信号和第二磁共振信号，确定各个回波时间的两个磁共振图像，并确定各个回波时间的两个磁共振图像之间的相位差；
37.第一涡流场确定模块，用于根据各相位差，确定至少一个成像回波时间下的涡流场分布；
38.涡流校正模块，用于利用成像回波时间下的涡流场分布进行涡流校正。
39.一种磁共振图像校正装置，该装置包括：
40.图像获取模块，用于获取扫描对象的多组磁共振图像，每组磁共振图像包括对应同一回波时间的第一磁共振图像和第二磁共振图像，上述第一磁共振图像为第一梯度施加后采集的第一磁共振信号重建获得，上述第二磁共振图像为第二梯度施加后采集的第二磁共振信号重建获得；
41.第二相位差确定模块，用于确定每组磁共振图像中第一磁共振图像和第二磁共振图像之间的相位差；
42.第二涡流场确定模块，用于根据上述相位差，确定至少一个成像回波时间的涡流场分布；
43.图像校正模块，用于根据上述涡流场分布对所述成像回波时间对应的待处理磁共振图像进行校正，获取校正的磁共振图像。
44.一种计算机设备，包括存储器和处理器，所述存储器存储有计算机程序，所述处理器执行所述计算机程序时实现以下步骤：
45.获取第一梯度施加后采集的多个回波时间的第一磁共振信号和第二梯度施加后采集的多个回波时间的第二磁共振信号；上述多个回波时间的第一磁共振信号的采集时刻与上述多个回波时间的第二磁共振信号的采集时刻相对应；
46.根据各个回波时间的第一磁共振信号和第二磁共振信号，确定各个回波时间的两个磁共振图像，并确定各个回波时间的两个磁共振图像之间的相位差；
47.根据各相位差，确定至少一个成像回波时间下的涡流场分布；
48.利用至少一个成像回波时间下的涡流场分布对成像回波时间对应的磁共振信号进行涡流校正。
49.一种计算机设备，包括存储器和处理器，所述存储器存储有计算机程序，所述处理器执行所述计算机程序时实现以下步骤：
50.获取扫描对象的多组磁共振图像，每组磁共振图像包括对应同一回波时间的第一磁共振图像和第二磁共振图像，上述第一磁共振图像为第一梯度施加后采集的第一磁共振信号重建获得，上述第二磁共振图像为第二梯度施加后采集的第二磁共振信号重建获得；
51.确定每组磁共振图像中第一磁共振图像和第二磁共振图像之间的相位差；
52.根据上述相位差，确定至少一个成像回波时间的涡流场分布；
53.根据上述涡流场分布对成像回波时间对应的待处理磁共振图像进行校正，获取校正的磁共振图像。
54.一种计算机可读存储介质，其上存储有计算机程序，所述计算机程序被处理器执行时实现以下步骤：
55.获取第一梯度施加后采集的多个回波时间的第一磁共振信号和第二梯度施加后采集的多个回波时间的第二磁共振信号；上述多个回波时间的第一磁共振信号的采集时刻与上述多个回波时间的第二磁共振信号的采集时刻相对应；
56.根据各个回波时间的第一磁共振信号和第二磁共振信号，确定各个回波时间的两个磁共振图像，并确定各个回波时间的两个磁共振图像之间的相位差；
57.根据各相位差，确定至少一个成像回波时间下的涡流场分布；
58.利用至少一个成像回波时间下的涡流场分布对成像回波时间对应的磁共振信号进行涡流校正。
59.一种计算机可读存储介质，其上存储有计算机程序，所述计算机程序被处理器执行时实现以下步骤：
60.获取扫描对象的多组磁共振图像，每组磁共振图像包括对应同一回波时间的第一磁共振图像和第二磁共振图像，上述第一磁共振图像为第一梯度施加后采集的第一磁共振信号重建获得，上述第二磁共振图像为第二梯度施加后采集的第二磁共振信号重建获得；
61.确定每组磁共振图像中第一磁共振图像和第二磁共振图像之间的相位差；
62.根据上述相位差，确定至少一个成像回波时间的涡流场分布；
63.根据上述涡流场分布对成像回波时间对应的待处理磁共振图像进行校正，获取校正的磁共振图像。
64.上述涡流校正方法、磁共振图像校正方法、装置、计算机设备和存储介质，通过采集时刻相对应的第一梯度施加后多个回波时间的第一磁共振信号和第二梯度施加后多个回波时间的第二磁共振信号，通过各回波时间的两个磁共振信号确定各回波时间的两个磁共振图像以及两个磁共振图像之间的相位差，根据各相位差确定成像回波时间下的涡流场分布，并利用成像回波时间下的涡流场分布对成像回波时间对应的磁共振信号进行涡流校正。在该方法中，由于可以通过各个回波时间不同梯度施加后的磁共振图像之间的相位差确定任意成像回波时间下的涡流场分布，这样在对涡流进行校正时，就可以准确地对任意成像回波时间的涡流进行涡流校正，这样就不会存在很大涡流残留，从而可以实现对涡流的准确校正，即可以提高对涡流进行校正的准确性。
附图说明
65.图1为一个实施例中计算机设备的内部结构图；
66.图2为一个实施例中涡流校正方法的流程示意图；
67.图3为另一个实施例中涡流校正步骤的流程示意图；
68.图3a为另一个实施例中对应回波时间下两个磁共振图像的相位差的示例图；
69.图3b为另一个实施例中球谐系数的幅度随时间变化的示例图；
70.图4为另一个实施例中涡流校正方法的流程示意图；
71.图4a为另一个实施例中磁共振信号的采集时序示例图；
72.图5为另一个实施例中涡流校正方法的流程示意图；
73.图6为另一个实施例中磁共振图像校正方法的流程示意图；
74.图6a为另一个实施例中所获得的待校正的多个弥散方向的弥散图像与合成图像；
75.图6b为另一个实施例中所获得的校正后的多个弥散方向的弥散图像与校正的合
成图像；
76.图7为一个实施例中涡流校正装置的结构框图。
具体实施方式
77.为了使本技术的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本技术进行进一步详细说明。应当理解，此处描述的具体实施例仅仅用以解释本技术，并不用于限定本技术。
78.本技术实施例提供的涡流校正方法，可以应用于磁共振系统或计算机设备，该磁共振系统可以包括相互连接的磁共振扫描设备和计算机设备，在对检测对象进行扫描时，可以通过磁共振扫描设备对检测对象进行扫描，并将扫描的数据传输至计算机设备进行处理，实现涡流校正。这里以涡流校正方法应用于计算机设备为例，该计算机设备可以是终端或服务器，以计算机设备是终端为例，其内部结构图可以如图1所示。该计算机设备包括通过系统总线连接的处理器、存储器、通信接口、显示屏和输入装置。其中，该计算机设备的处理器用于提供计算和控制能力。该计算机设备的存储器包括非易失性存储介质、内存储器。该非易失性存储介质存储有操作系统和计算机程序。该内存储器为非易失性存储介质中的操作系统和计算机程序的运行提供环境。该计算机设备的通信接口用于与外部的终端进行有线或无线方式的通信，无线方式可通过wifi、运营商网络、nfc(近场通信)或其他技术实现。该计算机程序被处理器执行时以实现一种涡流校正方法。该计算机设备的显示屏可以是液晶显示屏或者电子墨水显示屏，该计算机设备的输入装置可以是显示屏上覆盖的触摸层，也可以是计算机设备外壳上设置的按键、轨迹球或触控板，还可以是外接的键盘、触控板或鼠标等。
79.本领域技术人员可以理解，图1中示出的结构，仅仅是与本技术方案相关的部分结构的框图，并不构成对本技术方案所应用于其上的计算机设备的限定，具体的计算机设备可以包括比图中所示更多或更少的部件，或者组合某些部件，或者具有不同的部件布置。
80.需要说明的是，本技术实施例的执行主体可以是磁共振系统，也可以磁共振系统中的计算机设备，也可以是计算机设备内的涡流校正装置，当然还可以是其他装置，以下就以计算机设备为执行主体来对本技术的技术方案进行说明。
81.在一个实施例中，如图2所示，提供了一种涡流校正方法，该方法可以包括以下步骤：
82.s202，获取第一梯度施加后采集的多个回波时间的第一磁共振信号和第二梯度施加后采集的多个回波时间的第二磁共振信号。
83.其中，上述多个回波时间的第一磁共振信号的采集时刻与上述多个回波时间的第二磁共振信号的采集时刻相对应，这里的对应可以指的是第一梯度施加后的多个回波时间的采集时刻的数量和第二梯度施加后的多个回波时间的采集时刻的数量相等，例如，在第一梯度施加后的多个回波时间的采集时刻可以是t1-tn，在第二梯度施加后的多个回波时间的采集时刻也可以是t1-tn。这里多个回波时间的采集时刻的数量可以预先设定好，例如需要采集50个点，即50个采集时刻，这50个时刻之间可以是等间隔或不等间隔的时刻，例如可以是等间隔5ms、10ms等等的采集时刻。
84.这里第一梯度的强度大小和第二梯度的强度大小可以相等，也可以不相等。第一
梯度的方向和第二梯度的方向可以是相同的方向，也可以是相反的方向，当然也可以是其他类型的方向等等。
85.具体的，在实际对检测对象进行磁共振扫描时，可以通过磁共振设备中的梯度系统产生不同的梯度，并将产生的梯度施加给扫描对象，这里可以先施加第一梯度，并在各个回波时间的采集时刻下，采用预扫描序列采集该第一梯度施加后的磁共振信号，获得各个回波时间的第一磁共振信号。之后，可以给磁共振线圈施加第二梯度，并在与第一梯度施加后的采集时刻数量相同的各个回波时间的采集时刻下，采用相同的预扫描序列采集该第二梯度施加后的磁共振信号，获得各个回波时间的第二磁共振信号。这里获得的各第二磁共振信号与各第一磁共振信号的回波时间的采集时刻相对应，即对应的采集时刻下会获得一个第一磁共振信号和一个第二次磁共振信号。另外，这里的预扫描序列可以是gre序列，当然也可以是其他序列，例如se序列等等。
86.s204，根据各个回波时间的第一磁共振信号和第二磁共振信号，确定各个回波时间的两个磁共振图像，并确定各个回波时间的两个磁共振图像之间的相位差。
87.在本步骤中，在获得各个回波时间下的第一磁共振信号以及第二磁共振信号之后，可以通过图像重建等方式获得各个回波时间下的两个磁共振信号分别对应的磁共振图像，并通过磁共振图像上的参数获得各个磁共振图像的相位。之后，可以通过各个回波时间下的两个磁共振图像作差等方式获得每两个磁共振图像之间的相位差。
88.s206，根据各相位差，确定至少一个成像回波时间下的涡流场分布。
89.其中，上述回波时间指的是执行预扫描序列所设定的回波时间，例如可以是在1ms、2ms、3ms等等；这里的成像回波时间指的是执行成像扫描序列时的回波时间，其可以是任意的回波时间，例如可以是1.2ms、1.5ms、2.3ms等等。
90.在本步骤中，在获得各回波时间下每两个磁共振图像之间的相位差之后，这里每两个磁共振图像是对应回波时间的两个不同梯度施加后的磁共振图像，那么可以采用与涡流的时间特性和空间特性相关的数学算法对每两个磁共振图像的相位差进行数学运算处理，获得各个回波时间下任意空间位置上的涡流场，这里的回波时间即任意回波时间，即各个成像回波时间，那么也就是可以获得各个成像回波时间下任意空间位置上的涡流场。
91.s208，利用成像回波时间下的涡流场分布对成像回波时间对应的磁共振信号进行涡流校正。
92.在本步骤中，在获得各个回波时间下任意空间位置上的涡流场之后，也即是可以获得任意成像回波时间下任意空间位置上的涡流场，那么在实际进行涡流场校正时，就可以在执行成像扫描序列采集的磁共振信号后，采用这里获得的涡流场对任意成像回波时间下以及任意空间位置上的涡流进行校正。在此实施例中，磁共振信号可以指磁共振系统的接收线圈采集原始数据，也可只填充至k空间中的笛卡尔数据线、非笛卡尔数据线。对磁共振信号经过傅里叶变换可以得到待处理磁共振图像。可选的，利用成像回波时间下的涡流场分布对成像回波时间对应的磁共振信号进行涡流校正可以是：直接利用涡流场分布对磁共振信号进行校正，也可以是利用涡流场分布对待处理磁共振图像进行校正。
93.由上述描述可知，这里可以采用获得的涡流场分布对任意成像回波时间及任意空间位置的图像进行涡流校正，那么就不会存在对某一个或多个回波时间下，或某一个或多个空间位置上的图像没有进行涡流校正的情况，从而可以实现精确地对每个回波时间以及
每个空间位置上的图像进行涡流校正，因此最终可以实现准确地涡流校正，获得更准确地磁共振图像。
94.上述涡流校正方法中，通过采集时刻相对应的第一梯度施加后多个回波时间的第一磁共振信号和第二梯度施加后多个回波时间的第二磁共振信号，通过各回波时间的两个磁共振信号确定各回波时间的两个磁共振图像以及两个磁共振图像之间的相位差，根据各相位差确定一个或多个成像回波时间下的涡流场分布，并利用成像回波时间下的涡流场分布进行涡流校正。在该方法中，由于可以通过各个回波时间不同梯度施加后的磁共振图像之间的相位差确定任意成像回波时间下的涡流场分布，这样在对涡流进行校正时，就可以准确地对任意成像回波时间的涡流进行涡流校正，这样就不会存在很大涡流残留，从而可以实现对涡流的准确校正，即可以提高对涡流进行校正的准确性。
95.在另一个实施例中，提供了另一种涡流校正方法，在上述实施例的基础上，如图3所示，上述s206可以包括以下步骤：
96.s302，对各相位差进行处理，确定每个回波时间下的涡流场分布。
97.其中，涡流场分布与相位(或相位差)相关，可以采用如下公式(1)或公式(1)的变形表示，公式(1)如下所示：
98.bz＝φ/(2
·
pi
·
γ
·
te
·
2)
ꢀꢀꢀꢀꢀꢀꢀ
(1)
99.其中，pi指的是圆周率π；γ指的是旋磁比，在磁共振线圈设定好之后可以获知，为已知量；te指的是回波时间，指射频脉冲与相应的回波之间的时间间隔，也为已知量；φ为对应时刻下两个磁共振图像之间的相位差；bz为涡流场分布。
100.需要说明的是，这里对应回波时间下两个磁共振图像之间的相位差为一个可以表征空间的相位差，即该相位差并不简单的是一个数值，可以是一个相位空间分布。示例地，参见图3a所示，在回波时间0-200ms内，每间隔10ms采集一次磁共振信号(即间隔10ms获得两个磁共振图像)，采集20次，即20个时刻，那么可以获得20个采集时刻下两个磁共振图像之间的相位差，所获得的20个相位差分布图如图3a所示。另外，这里图3a只是示例，并不影响本技术实施例的实质内容。
101.在上述获得各个对应回波时间下两个磁共振图像之间的相位差以及相关计算参数之后，由于相位差是带有空间分布的，即可以表征空间位置，那么通过公式(1)或公式(1)的变形进行计算，就可以获得各个成像回波时间下，任意空间位置上的涡流场分布bz。
102.s304，对每个回波时间下的涡流场分布进行求解，确定一个或多个成像回波时间下的涡流场。
103.在本步骤中，可选的，可以采用以下步骤a1和a2来求解多个回波时间下任意空间位置上的涡流场：
104.步骤a1，根据预设的球谐函数对每个回波时间下的涡流场分布进行系数展开，确定各个成像回波时间下的涡流球谐系数。
105.其中，球谐函数与各个点的空间位置相关，可以采用如下公式(2)或公式(2)的变形表示，公式(2)如下所示：
[0106][0107]
其中，这里的bz(x,y,z)与上述的bz相同，只不过这里的bz(x,y,z)中表示了每个点的空间位置x、y、z。公式(2)中的第一项一般被称为涡流b0项，第二到第四项一般称为涡流线性项或一阶项，通过化简之后这三项正好对应于梯度线圈(磁共振线圈)产生了x、y、z三个方向上的梯度。阶数高于以上四项的一般被称为涡流高阶项。
[0108]
具体的，在获得各个对应回波时间下，任意空间位置上的涡流场分布bz(x,y,z)之后，那么也就可以获知各个点的空间位置，之后，可以将各个点的bz(x,y,z)均采用公式(2)或公式(2)的变形表示，那么这些表示的公式中就只有公式(2)中等这些系数为未知量，然后就可以联立这多个表示的公式对这些系数进行求解，就可以获得这些系数的值，由于这里的系数是各个成像回波时间下的系数，因此可以记为各个成像回波时间下的涡流球谐系数。
[0109]
需要说明的是，这里获得的每个球谐系数等均是与时间相关的，即每个球谐系数在每个成像回波时间的采集时刻均会有一个值。参见图3b所示，在0-200ms不同的回波时间内，通过上述求解过程，可以获得各个球谐系数的幅度amplitude(幅值)随着时间time(ms)的变化趋势。另外，这里图3b只是示例，并不影响本技术实施例的实质内容。
[0110]
步骤a2，根据上述各个成像回波时间下的涡流球谐系数，得到各个成像回波时间下的涡流场分布。
[0111]
在本步骤中，在获得各个成像回波时间下的涡流球谐系数之后，通常涡流在时间上一般采用多指数衰减函数模型，一般采用如下公式(3)或公式(3)的变形表示：
[0112][0113]
其中，g
eddy
(t)为与时间相关的涡流大小；e(t)为脉冲冲击响应；n为多指数涡流成分，为已知量；α、τ为每一种涡流成分的幅度和时间常数，也可以是已知量。
[0114]
那么可以采用如上公式(3)或公式(3)的变形获得任意成像回波时间下的涡流，同时结合上面的各个成像回波时间下的涡流球谐系数代入至公式(2)中，就可以获得各个成像回波时间下任意空间位置上的涡流场。在需要获得任意成像回波时间下任意空间位置上的涡流场时，也可以从上述涡流球谐系数中选取相应相邻成像回波时间相应空间位置上的涡流场，进而采用插值算法等拟合算法对上述已知成像回波时间已知空间位置上的涡流场进行拟合运算，获得所需成像回波时间所需空间位置上的涡流场，即获得各个成像回波时间下任意空间位置上的涡流场。
[0115]
本实施例中，通过采用于相位相关的球谐函数对各回波时间下的两个磁共振图像的相位差进行数学运算处理，获得每个回波时间下任意空间位置上的涡流场球谐系数，并
通过多指数衰减模型或者插值算法对每个对应时刻下任意空间位置上的涡流场进行求解获得各个成像回波时间下任意空间位置上的涡流场。这里采用球谐函数对相位差进行数学运算处理确定任意成像回波时间任意空间位置上的涡流场，计算过程直观准确，从而可以进一步提升获得的任意成像回波时间任意空间位置上的涡流场的准确性。进一步地，采用与空间位置相关的球谐函数对各成像回波时间下任意空间位置上的涡流场进行系数求解，可以准确获得每一个球谐系数，从而可以便于后续快速获得各个成像回波时间下任意空间位置上的涡流场。
[0116]
在另一个实施例中，提供了另一种涡流校正方法，在上述实施例的基础上，如图4所示，上述s202可以包括以下步骤：
[0117]
s402，在第一梯度施加后，采用预扫描序列进行多层数据交叉采集，获得上述第一梯度施加后多个回波时间的第一磁共振信号。
[0118]
在本步骤中，在进行数据采集时，采用的是预扫描序列是gre序列，这里采用gre序列进行数据采集，可以实现数据的每一层之间交互激发，在一个时间段内可以同时采集多个时刻的磁共振信号，这样可以提高数据采集效率。
[0119]
另外，这里的多层数据交叉采集指的是在各个回波时间的采集时刻下，将各层的数据进行交叉采集，即各层的数据是融合在一起进行采集的，并不是一层数据采集完才采集另外一层数据，即并不是单一层数据激发采集。
[0120]
示例地，以三层数据采集为例，三层数据分别为slice1、slice2、slice3，的回波时间从t1-tn，每一层数据的相位编码方向的数据分别填充在pe1-peend(pe1-peend分别表述k空间中沿着相位编码方向的填充区域)，那么在采用本实施例的多层数据交叉采集时，图4a所对应的整个序列的采集时序如下所示：
[0121][0122]
具体的，如上所示，整个预扫描包含三次不同的激发，每次激发检测对象的多个层面，且两个或多个层面的回波时间不同：在第一次激发过程中同时采集slice1-slice3三个层面的不同回波时间的磁共振信号，在此实施例中分别采集slice1-t1、slice2-t12、slice3-t3以及slice1-t4的磁共振信号，且属于同次激发的同一片层的不同回波填充在k空间中的不同位置(如上表中的pe1、pe2等)，属于不同片层和/或不同回波时间的磁共振信号分别填充在不同k空间中；第二次激发分别采集slice2-t1、slice3-t12、slice1-t3以及slice2-t4的磁共振信号；第三次激发分别采集slice3-t1、slice1-t12、slice2-t3以及slice3-t4的磁共振信号。采用gre序列通过上述的方式对第一梯度下各个回波时间的磁共振信号进行多层数据交叉采集，就可以得到在第一梯度施加后各个回波时间的磁共振信
号，均记为第一磁共振信号。
[0123]
由上述描述可知，上述在第一梯度施加后下可以同时采集多个回波时间的磁共振信号，即可以同时获得多个回波时间下的磁共振图像，例如从t1-tn多个回波时间的采集时刻内采集了n个被磁共振信号填充的k空间，相对现有技术一个回波时间采集一个被磁共振信号填充的k空间，这里可以将磁共振信号提高n倍，即可以提高磁共振信号或图像的采集效率。
[0124]
另外，采用现有的数据采集方式，采集每一层数据的时间间隔假设是tr，那么由上述整个序列的采集时序可知，本实施例在采集每一层数据时的时间间隔为n*tr，这里的n是采集的数据的层数(例如在上述的示例中，n为3)，由此可见，采用本实施例的方案在采集同一层的磁共振信号时的时间间隔扩大了，每一层的纵向磁化矢量有足够的时间恢复，可避免不同次激发之间的相互影响，每一采集的磁共振信号可以获得信噪比高的图像，减少多次激发产生的信号失真。
[0125]
s404，在第二梯度施加后，采用预扫描序列进行多层数据交叉采集，获得上述第二梯度施加后多个回波时间的第二磁共振信号。
[0126]
在本步骤中，可以采用与上述s402相同的方式，在同样的多个回波时间的采集时刻内，采用gre序列通过上述的方式对第二梯度施加后各个回波时间的磁共振信号进行多层数据交叉采集，就可以得到在第二梯度施加后各个回波时间的磁共振信号，均记为第二磁共振信号。
[0127]
示例地，参见图4a所示，数据采集层数为n，每一层于不同回波时间(te)中分别执行一次激发采集得到的磁共振信号分别获得多个k空间的数据，图中tr(repetition time)表示相邻的两次执行的时间间隔。在第一梯度施加后(例如gtest梯度下)对磁共振信号进行多层数据交叉采集，可以获得不同回波时间的采集时刻下的第一磁共振信号，之后，在第二梯度施加后(例如将gtest梯度反向的梯度)对磁共振信号进行多层数据交叉采集，可以获得不同回波时间的采集时刻下的第二磁共振信号。另外，需要说明的是，这里图4a只是示例，并不影响本技术实施例的实质内容。
[0128]
本实施例中，通过采用预扫描序列对第一梯度以及第二梯度施加后多个回波时间的磁共振信号进行多层数据交叉采集，就可以获得第一梯度施加后多个回波时间的第一磁共振信号以及第二梯度施加后多个回波时间的第二磁共振信号。在本实施例中，一方面，这里在第一梯度以及第二梯度施加后可以同时采集多个回波时间的磁共振信号，即可以同时采集多个回波时间的磁共振图像，从而可以提高信号或图像的采集效率。另一方面，预扫描序列采用gre序列在第一梯度以及第二梯度施加后进行多层数据交叉采集，磁化矢量的有效恢复时间从tr变为n*tr，从而可以有效提高图像质量的信噪比。
[0129]
在另一个实施例中，提供了另一种涡流校正方法，在上述实施例的基础上，上述第一梯度和上述第二梯度为大小相等，且方向相反的两个梯度；上述s206具体可以包括以下步骤b：
[0130]
步骤b，根据各相位差，确定上述第一梯度施加后或第二梯度施加后各个成像回波时间下的涡流场分布。
[0131]
在本步骤中，第一梯度的强度与第二梯度的强度大小相等，且方向相反。那么本步骤中可以在上述s302-s304的步骤下，获得第一梯度下各个成像回波时间下任意空间位置
上的涡流场分布，由于第一梯度和第二梯度的强度相等，那么也可以记为是获得第二梯度施加后各个成像回波时间下任意空间位置上的涡流场分布。也就是说，这里可以获得某一个梯度下的各个成像回波时间下任意空间位置上的涡流场分布。
[0132]
通常，梯度系统是一个线性系统，那么上述在获得某一个梯度施加后的各个成像回波时间下任意空间位置上的涡流场分布之后，可以通过对该梯度施加后的各个成像回波时间下任意空间位置上的涡流场分布进行线性计算，就可以获得任意梯度下的各个成像回波时间下任意空间位置上的涡流场分布。具体的计算方式可以包括以下计算步骤c1和c2：
[0133]
步骤c1，获取预设的第三梯度的强度和第一梯度的强度之间的第一强度比例，或第三梯度的强度和上述第二梯度的强度之间的第二强度比例。
[0134]
在本步骤中，在需要获得任意梯度施加后的涡流场分布时，可以预先获得该任意梯度及其强度，这里记为第三梯度。之后，可以将第三梯度的强度与第一梯度的强度做除法，获得第三梯度的强度在第一梯度的强度上所占的比例，记为第一强度比例。同样的，也可以获得第三梯度的强度在第二梯度的强度上所占的比例，记为第二强度比例。这里由于第一梯度的强度和第二梯度的强度大小相等，所以第一强度比例和第二强度比例通常也是相等的。
[0135]
步骤c2，根据上述第一强度比例或第二强度比例，对上述第一梯度施加后和第二梯度施加后各个成像回波时间下的涡流场分布进行线性叠加处理，确定上述第三梯度下各个成像回波时间下的涡流场分布。
[0136]
在本步骤中，在获得第一强度比例或第二强度比例之后，同时也可以获得第一梯度或第二梯度施加后的各个成像回波时间下任意空间位置上的涡流场，这样就可以通过第一强度比例或第二强度比例与第一梯度或第二梯度施加后各个成像回波时间下任意空间位置上的涡流场做乘积，获得的乘积即为第三梯度施加后各个成像回波时间下任意空间位置上的涡流场。
[0137]
本实施例中，第一梯度的强度和第二梯度的强度大小相等，方向相反，那么可以通过上述各个成像回波时间下的两个磁共振图像的相位差，获得该第一梯度或第二梯度施加后各个成像回波时间下任意空间位置上的涡流场分布，通过此方式可以获得某一特定梯度下的各个成像回波时间下任意空间位置上的涡流场分布，便于对该梯度下的涡流进行准确校正。进一步地，可以通过任意梯度与已知第一梯度或第二梯度之间的强度比例，获得任意梯度下的各个成像回波时间下任意空间位置上的涡流场，这样可以不受涡流的梯度的强度限制，提高涡流校正的适用范围。
[0138]
在另一个实施例中，提供了另一种涡流校正方法，在上述实施例的基础上，如图5所示，上述s204可以包括以下步骤：
[0139]
s502，对各个回波时间的第一磁共振信号和第二磁共振信号进行图像重建，获得各个回波时间的第一磁共振图像和第二磁共振图像。
[0140]
在本步骤中，可以在获得每个第一磁共振信号后分别对当前获得的第一磁共振信号进行图像重建，获得当前回波时间下的第一磁共振图像，这样最终可以获得每一个片层所对应的各个回波时间的第一磁共振图像；当然，也可以是在获得各个回波时间全部的第一磁共振信号之后，分别对每一个片层各个回波时间的第一磁共振信号进行图像重建，获得每一个片层各个回波时间的第一磁共振图像。
[0141]
同样的，可以采用上述方式对每一个片层各个回波时间的第二磁共振信号进行图像重建，获得每一个片层各个回波时间的第二磁共振图像。
[0142]
这里获得的各个回波时间的第一磁共振图像与各个回波时间的第二磁共振图像一一对应，两者对应同一片层。例如第一个回波时间获得的第一磁共振图像与第一回波时间获得的第二磁共振图像为对应回波时间获得的两个磁共振图像，第m个回波时间获得的第一磁共振图像与第m回波时间获得的第二磁共振图像为对应时刻获得的两个磁共振图像，m为任意正整数。
[0143]
另外，在对磁共振信号进行图像重建时，可以选取任意的图像重建算法进行图像重建，以获得第一磁共振图像和第二磁共振图像。
[0144]
s504，根据各个回波时间的第一磁共振图像获得各个回波时间的第一相位，以及根据各个回波时间的第二磁共振图像获得各个回波时间的第二相位。
[0145]
在本步骤中，在上述对各个回波时间的磁共振信号进行图像重建获得磁共振图像之后，也可以获得磁共振图像对应的相位，这里每二个磁共振图像的相位一般是可以表征空间位置的相位分布，一般每个回波时间获得的相位可能是不同的。
[0146]
其中，各个回波时间的第一磁共振图像的相位均可以记为第一相位，各个回波时间的第二磁共振图像的相位均可以记为第二相位，这样可以获得每个对应回波时间的第一相位及第二相位。
[0147]
s506，对各个回波时间的第一相位和第二相位进行求差运算，得到各个回波时间的第一磁共振图像和第二磁共振图像之间的相位差。
[0148]
在本步骤中，在获得每个对应回波时间的第一相位及第二相位之后，可以对每个对应回波时间的第一相位及第二相位进行求差运算，求差运算可以是第一相位减去第二相位，也可以是第二相位减去第一相位，总之可以得到每个对应回波时间下的两个相位之间的差值，即获得各个对应回波时间的第一磁共振图像和第二磁共振图像之间的相位差。
[0149]
本实施例中，通过对不同对应回波时间的第一磁共振信号以及第二磁共振信号进行图像重建，获得各个对应回波时间下的两个磁共振图像及各自的相位，并对回波时间相对应的两个磁共振图像的两个相位进行求差运算，获得回波时间相对应的两个磁共振图像之间的相位差，该计算过程简单直观，因此可以提高获得两个磁共振图像之间的相位差的效率，进一步可以提高涡流校正的效率。
[0150]
以上实施例中介绍了对涡流进行校正，在此基础上还可以对磁共振图像进行校正，以下就对该过程进行说明。
[0151]
在另一个实施例中，提供了一种磁共振图像校正方法，在上述实施例的基础上，参见图6所示，该方法可以包括以下步骤：
[0152]
s602，获取扫描对象的多组磁共振图像，每组磁共振图像包括对应同一回波时间的第一磁共振图像和第二磁共振图像，上述第一磁共振图像为第一梯度施加后采集的第一磁共振信号重建获得，上述第二磁共振图像为第二梯度施加后采集的第二磁共振信号重建获得。扫描对象包括多个片层，多组磁共振图像则包含每个片层的不同回波时间下一系列的第一磁共振图像和第二磁共振图像。
[0153]
其中，扫描对象可以是人体或动物体的身体的任意一个或多个部位。可以通过对扫描对象施加第一梯度，并在该第一梯度下获取各层的第一磁共振信号，即获取各回波时
间下的第一磁共振信号，并对每层的第一磁共振信号进行图像重建获得各第一磁共振图像。
[0154]
对于各第二磁共振图像的获取方式可以与第一磁共振图像获取方式相同，区别仅在于施加了宇第一梯度不同的第二梯度，这里就不再赘述。
[0155]
s604，确定每组磁共振图像中第一磁共振图像和第二磁共振图像之间的相位差。
[0156]
本步骤的解释说明可以参见上述s204的解释说明，这里不再赘述。
[0157]
s606，根据上述相位差，确定至少一个成像回波时间的涡流场分布。
[0158]
本步骤的解释说明可以参见上述s206的解释说明，这里不再赘述。
[0159]
s608，根据上述涡流场分布对成像回波时间对应的待处理磁共振图像进行校正，获取校正的磁共振图像。
[0160]
在本步骤中，以待处理磁共振图像为弥散图像为例，在获得至少一个成像回波时间的涡流场分布之后，就可以采用该涡流场分布对待处理的弥散图像进行图像校正，获得校正后的弥散图像。
[0161]
在一个实施例中，成像扫描序列采用平面回波弥散加权成像(epi-dwi)序列，其激发检测对象后采集得到的不同弥散方向的磁共振信号，将不同弥散方向的磁共振信号经过重建得到多个初始弥散图像(该多个初始弥散图像为待处理的磁共振图像)，可采用涡流场分布分别对多个初始弥散图像进行校正，得到多个校正的弥散图像；以及对多个校正的弥散图像进行合并处理，获得合成的弥散图像。
[0162]
在一个实施例中，可将不同弥散方向的磁共振信号经过重建与合并获得的弥散加权(dwi)图像为待处理的磁共振图像，该dwi图像由多个方向(以三方向为例)弥散图像的像素逐点相乘再开立方得到，如下：
[0163][0164]
其中，vox表示dwi图像中的任一像素点的值；d1、d2和d3分别表示三个不同的弥散方向对应的像素点的值。由于不同方向时，弥散梯度的大小不同，会在采集图像时刻产生不同的涡流，dwi序列采集时，沿着相位编码(pe)方向带宽很低。因而，会造成图像的pe方向有位置偏移。又由于不同方向弥散梯度涡流不同，因而，合成的dwi图像就会模糊。示例性的，待处理的磁共振图像的偏移可以表示为：
[0165][0166]
其中，δy表示沿着相位编码方向位移；t
esp
为epi_dwi的相邻回波时间间隔；gyτy表示pe方向尖峰脉冲的零阶矩。根据前述涡流场分布bz即可得到dwi图像中沿着相位编码方向的位移，根据该相位编码方向的位移即可对dwi图像进行校正。
[0167]
如图6a所示，从左至右分别为本技术一实施所得到的方向1的初始弥散图像、方向2的初始弥散图像、方向3的初始弥散图像以及合成dwi图像。其中，初始弥散图像由于受到不同大小涡流的影响，产生了不同的形变，合成dwi图像存在较为严重的模糊。采用本技术实施例附图6的流程，对合成dwi图像校正涡流影响。如图6b所示，从左至右分别为校正后的dwi图像所对应的方向1的初始弥散图像、方向2的初始弥散图像、方向3的初始弥散图像以及校正后的dwi图像。经过校正处理的弥散图像的单个方向的图像位置趋于一致，校正后的dwi图像边界变得清晰。
[0168]
本实施例中，通过获取扫描对象的多个片层的多组磁共振图像，每组磁共振图像包括对应同一回波时间的第一磁共振图像和第二磁共振图像，确定每组磁共振图像中第一磁共振图像和第二磁共振图像之间的相位差，以及通过相位差确定各成像回波时间的涡流场分布，并通过涡流场分布对待处理弥散图像进行校正，获得校正后的弥散图像。通过这样的方式可以实现精确地对每个回波时间以及每个空间位置上的图像进行涡流校正，并在此基础上对弥散图像进行校正，获得更准确地校正图像。
[0169]
应该理解的是，虽然图2、3、4、5、6的流程图中的各个步骤按照箭头的指示依次显示，但是这些步骤并不是必然按照箭头指示的顺序依次执行。除非本文中有明确的说明，这些步骤的执行并没有严格的顺序限制，这些步骤可以以其它的顺序执行。而且，图2、3、4、5、6中的至少一部分步骤可以包括多个步骤或者多个阶段，这些步骤或者阶段并不必然是在同一时刻执行完成，而是可以在不同的时刻执行，这些步骤或者阶段的执行顺序也不必然是依次进行，而是可以与其它步骤或者其它步骤中的步骤或者阶段的至少一部分轮流或者交替地执行。
[0170]
在一个实施例中，如图7所示，提供了一种涡流校正装置，包括：信号获取模块10、第一相位差确定模块11、第一涡流场确定模块12和涡流校正模块13，其中：
[0171]
信号获取模块10，用于获取第一梯度施加后采集的多个回波时间的第一磁共振信号和第二梯度施加后采集的多个回波时间的第二磁共振信号；上述多个回波时间的第一磁共振信号的采集时刻与上述多个回波时间的第二磁共振信号的采集时刻相对应；
[0172]
第一相位差确定模块11，用于根据各个回波时间的第一磁共振信号和第二磁共振信号，确定各个回波时间的两个磁共振图像，并确定各个回波时间的两个磁共振图像之间的相位差；
[0173]
第一涡流场确定模块12，用于根据各相位差，确定一个或多个成像回波时间下的涡流场分布；
[0174]
涡流校正模块13，用于利用成像回波时间下的涡流场分布进行涡流校正。
[0175]
关于涡流校正装置的具体限定可以参见上文中对于涡流校正方法的限定，在此不再赘述。
[0176]
在另一个实施例中，提供了另一种涡流校正装置，在上述实施例的基础上，上述第一涡流场确定模块12可以包括数学运算处理单元和涡流场确定单元，其中：
[0177]
数学运算处理单元，用于对各相位差进行数学运算处理，确定每个回波时间下的涡流场分布；
[0178]
涡流场确定单元，用于对每个回波时间下的涡流场分布进行求解，确定各个成像回波时间下的涡流场分布。
[0179]
可选的，上述涡流场确定单元可以包括系数展开子单元和涡流场确定子单元，其中：
[0180]
系数展开子单元，用于根据预设的球谐函数对每个回波时间下的涡流场分布进行系数展开，确定各个成像回波时间下的涡流球谐系数；
[0181]
涡流场确定子单元，用于根据上述各个成像回波时间下的涡流球谐系数，得到各个成像回波时间下的涡流场分布。
[0182]
在另一个实施例中，提供了另一种涡流校正装置，在上述实施例的基础上，上述获
取模块10可以包括第一信号采集单元和第二信号采集单元，其中：
[0183]
第一信号采集单元，用于在第一梯度施加后，采用预扫描序列进行多层数据交叉采集，获得上述第一梯度施加后多个回波时间的第一磁共振信号；
[0184]
第二信号采集单元，用于在第二梯度施加后，采用预扫描序列进行多层数据交叉采集，获得上述第二梯度施加后多个回波时间的第二磁共振信号。
[0185]
在另一个实施例中，提供了另一种涡流校正装置，在上述实施例的基础上，上述第一梯度和上述第二梯度为梯度强度大小相等，且方向相反的两个梯度；上述涡流场确定模块12，具体用于根据各相位差，确定上述第一梯度施加后或第二梯度施加后各个成像回波时间下的涡流场分布。
[0186]
在另一个实施例中，提供了另一种涡流校正装置，在上述实施例的基础上，上述相位差确定模块11可以包括图像重建单元、相位获取单元和相位差确定单元，其中：
[0187]
图像重建单元，用于对各个回波时间的第一磁共振信号和第二磁共振信号进行图像重建，获得各个回波时间的第一磁共振图像和第二磁共振图像；
[0188]
相位获取单元，用于根据各个回波时间的第一磁共振图像获得各个回波时间的第一相位，以及根据各个回波时间的第二磁共振图像获得各个回波时间的第二相位；
[0189]
相位差确定单元，用于对各个回波时间的第一相位和第二相位进行求差运算，得到各个回波时间的第一磁共振图像和第二磁共振图像之间的相位差。
[0190]
关于涡流校正装置的具体限定可以参见上文中对于涡流校正方法的限定，在此不再赘述。
[0191]
在另一个实施例中，提供了一种磁共振图像校正装置，在上述实施例的基础上，该装置包括：
[0192]
图像获取模块，用于获取扫描对象的多个片层的多组磁共振图像，每组磁共振图像包括对应同一回波时间的第一磁共振图像和第二磁共振图像，上述第一磁共振图像为第一梯度施加后采集的第一磁共振信号重建获得，上述第二磁共振图像为第二梯度施加后采集的第二磁共振信号重建获得；
[0193]
第二相位差确定模块，用于确定每组磁共振图像中第一磁共振图像和第二磁共振图像之间的相位差；
[0194]
第二涡流场确定模块，用于根据上述相位差，确定至少一个成像回波时间的涡流场分布；
[0195]
图像校正模块，用于根据上述涡流场分布对成像回波时间对应的待处理磁共振图像进行校正，获取校正的磁共振图像。
[0196]
上述涡流校正装置以及磁共振图像校正装置中的各个模块可全部或部分通过软件、硬件及其组合来实现。上述各模块可以硬件形式内嵌于或独立于计算机设备中的处理器中，也可以以软件形式存储于计算机设备中的存储器中，以便于处理器调用执行以上各个模块对应的操作。
[0197]
在一个实施例中，提供了一种计算机设备，包括存储器和处理器，存储器中存储有计算机程序，该处理器执行计算机程序时实现以下步骤：
[0198]
获取第一梯度施加后采集的多个回波时间的第一磁共振信号和第二梯度施加后采集的多个回波时间的第二磁共振信号；上述多个回波时间的第一磁共振信号的采集时刻
与上述多个回波时间的第二磁共振信号的采集时刻相对应；根据各个回波时间的第一磁共振信号和第二磁共振信号，确定各个回波时间的两个磁共振图像，并确定各个回波时间的两个磁共振图像之间的相位差；根据各相位差，确定一个或多个成像回波时间下的涡流场分布；利用成像回波时间下的涡流场分布对成像回波时间对应的磁共振信号进行涡流校正。
[0199]
在一个实施例中，处理器执行计算机程序时还实现以下步骤：
[0200]
对各相位差进行数学运算处理，确定每个回波时间下的涡流场分布；对每个回波时间下的涡流场分布进行求解，确定各个成像回波时间下的涡流场分布。
[0201]
在一个实施例中，处理器执行计算机程序时还实现以下步骤：
[0202]
根据预设的球谐函数对每个回波时间下的涡流场分布进行系数展开，确定各个成像回波时间下的涡流球谐系数；根据上述各个成像回波时间下的涡流球谐系数，得到各个成像回波时间下的涡流场分布。
[0203]
在一个实施例中，处理器执行计算机程序时还实现以下步骤：
[0204]
在第一梯度施加后，采用预扫描序列进行多层数据交叉采集，获得上述第一梯度施加后多个回波时间的第一磁共振信号；在第二梯度施加后，采用预扫描序列进行多层数据交叉采集，获得上述第二梯度施加后多个回波时间的第二磁共振信号。
[0205]
在一个实施例中，处理器执行计算机程序时还实现以下步骤：
[0206]
根据各相位差，确定上述第一梯度施加后或第二梯度施加后各个成像回波时间下的涡流场分布。
[0207]
在一个实施例中，处理器执行计算机程序时还实现以下步骤：
[0208]
对各个回波时间的第一磁共振信号和第二磁共振信号进行图像重建，获得各个回波时间的第一磁共振图像和第二磁共振图像；根据各个回波时间的第一磁共振图像获得各个回波时间的第一相位，以及根据各个回波时间的第二磁共振图像获得各个回波时间的第二相位；对各个回波时间的第一相位和第二相位进行求差运算，得到各个回波时间的第一磁共振图像和第二磁共振图像之间的相位差。
[0209]
在一个实施例中，处理器执行计算机程序时还实现以下步骤：
[0210]
获取扫描对象的多个片层的多组磁共振图像，每组磁共振图像包括对应同一回波时间的第一磁共振图像和第二磁共振图像，上述第一磁共振图像为第一梯度施加后采集的第一磁共振信号重建获得，上述第二磁共振图像为第二梯度施加后采集的第二磁共振信号重建获得；确定每组磁共振图像中第一磁共振图像和第二磁共振图像之间的相位差；根据上述相位差，确定至少一个成像回波时间的涡流场分布；根据上述涡流场分布对成像回波时间对应的待处理磁共振图像进行校正，获取校正的磁共振图像。
[0211]
在一个实施例中，提供了一种计算机可读存储介质，其上存储有计算机程序，计算机程序被处理器执行时实现以下步骤：
[0212]
获取第一梯度施加后采集的多个回波时间的第一磁共振信号和第二梯度施加后采集的多个回波时间的第二磁共振信号；上述多个回波时间的第一磁共振信号的采集时刻与上述多个回波时间的第二磁共振信号的采集时刻相对应；根据各个回波时间的第一磁共振信号和第二磁共振信号，确定各个回波时间的两个磁共振图像，并确定各个回波时间的两个磁共振图像之间的相位差；根据各相位差，确定一个或多个成像回波时间下的涡流场
分布；利用成像回波时间下的涡流场分布对成像回波时间对应的待处理磁共振图像进行校正，获取校正的磁共振图像进行涡流校正。
[0213]
在一个实施例中，计算机程序被处理器执行时还实现以下步骤：
[0214]
对各相位差进行数学运算处理，确定每个回波时间下的涡流场分布；对每个回波时间下的涡流场分布进行求解，确定各个成像回波时间下的涡流场分布。
[0215]
在一个实施例中，计算机程序被处理器执行时还实现以下步骤：
[0216]
根据预设的球谐函数对每个回波时间下的涡流场分布进行系数展开，确定各个成像回波时间下的涡流球谐系数；根据上述各个成像回波时间下的涡流球谐系数，得到各个成像回波时间下的涡流场分布。
[0217]
在一个实施例中，计算机程序被处理器执行时还实现以下步骤：
[0218]
在第一梯度施加后，采用预扫描序列进行多层数据交叉采集，获得上述第一梯度施加后多个回波时间的第一磁共振信号；在第二梯度施加后，采用预扫描序列进行多层数据交叉采集，获得上述第二梯度施加后多个回波时间的第二磁共振信号。
[0219]
在一个实施例中，计算机程序被处理器执行时还实现以下步骤：
[0220]
根据各相位差，确定上述第一梯度施加后或第二梯度施加后各个成像回波时间下的涡流场分布。
[0221]
在一个实施例中，计算机程序被处理器执行时还实现以下步骤：
[0222]
对各个回波时间的第一磁共振信号和第二磁共振信号进行图像重建，获得各个回波时间的第一磁共振图像和第二磁共振图像；根据各个回波时间的第一磁共振图像获得各个回波时间的第一相位，以及根据各个回波时间的第二磁共振图像获得各个回波时间的第二相位；对各个回波时间的第一相位和第二相位进行求差运算，得到各个回波时间的第一磁共振图像和第二磁共振图像之间的相位差。
[0223]
在一个实施例中，计算机程序被处理器执行时还实现以下步骤：
[0224]
获取扫描对象的多个片层的多组磁共振图像，每组磁共振图像包括对应同一回波时间的第一磁共振图像和第二磁共振图像，上述第一磁共振图像为第一梯度施加后采集的第一磁共振信号重建获得，上述第二磁共振图像为第二梯度施加后采集的第二磁共振信号重建获得；确定每组磁共振图像中第一磁共振图像和第二磁共振图像之间的相位差；根据上述相位差，确定至少一个成像回波时间的涡流场分布；根据上述涡流场分布对成像回波时间对应的待处理磁共振图像进行校正，获取校正的磁共振图像。
[0225]
本领域普通技术人员可以理解实现上述实施例方法中的全部或部分流程，是可以通过计算机程序来指令相关的硬件来完成，所述的计算机程序可存储于一非易失性计算机可读取存储介质中，该计算机程序在执行时，可包括如上述各方法的实施例的流程。其中，本技术所提供的各实施例中所使用的对存储器、存储、数据库或其它介质的任何引用，均可包括非易失性和易失性存储器中的至少一种。非易失性存储器可包括只读存储器(read-only memory，rom)、磁带、软盘、闪存或光存储器等。易失性存储器可包括随机存取存储器(random access memory，ram)或外部高速缓冲存储器。作为说明而非局限，ram可以是多种形式，比如静态随机存取存储器(static random access memory，sram)或动态随机存取存储器(dynamic random access memory，dram)等。
[0226]
以上实施例的各技术特征可以进行任意的组合，为使描述简洁，未对上述实施例
中的各个技术特征所有可能的组合都进行描述，然而，只要这些技术特征的组合不存在矛盾，都应当认为是本说明书记载的范围。
[0227]
以上所述实施例仅表达了本技术的几种实施方式，其描述较为具体和详细，但并不能因此而理解为对发明专利范围的限制。应当指出的是，对于本领域的普通技术人员来说，在不脱离本技术构思的前提下，还可以做出若干变形和改进，这些都属于本技术的保护范围。因此，本技术专利的保护范围应以所附权利要求为准。