一种前瞻性采集校正序列触发的磁共振成像方法和装置与流程

1.本发明涉及医学器械技术领域，特别是涉及一种前瞻性采集校正(prospective acquisition correction，pace)序列触发的磁共振成像方法和装置。


背景技术：

2.磁共振成像(magnetic resonance imaging，mri)是随着计算机技术、电子电路技术、超导体技术的发展而迅速发展起来的一种生物磁学核自旋成像技术。它利用磁场与射频脉冲使人体组织内进动的氢核(即h )发生振动产生射频信号，经计算机处理而成像。当把物体放置在磁场中，用适当的电磁波照射它，使之共振，然后分析它释放的电磁波，就可以得知构成这一物体的原子核的位置和种类，据此可以绘制成物体内部的精确立体影像。比如，可以通过磁共振成像扫描人类大脑获得的一个连续切片的动画，由头顶开始，一直到脚部。
3.单激发快速自旋回波序列是一种高速磁共振成像序列，适用于许多临床应用，尤其是运动敏感成像应用。仅激发一次单激发快速自旋回波序列，就可以获得成像对象的一层完整图像。单激发快速自旋回波序列通常包括半傅里叶采集单次激发快速自旋回波序列(half-fourier acquisition single-shot turbo spin-echo，haste)、单次激发tse(turbo spin echo-single-shot，tse-ssh)序列、单次激发fse序列(single shot fast spin echo，ss-fse)序列，等等。在对胸腹部进行磁共振成像时，由于生理运动的影响(例如患者的呼吸运动)，图像中经常出现伪影。因此，需要对检查对象的生理运动进行监测，从而消除由此产生的伪影。通过结合pace序列的触发机制，单激发快速自旋回波序列可以获取几乎没有运动伪影的成像图像。
4.在目前的单激发快速自旋回波序列中，当从pace序列接收到触发信号后，单激发快速自旋回波序列在每个生理运动周期内分别获取单层的成像数据，具有成像时间长的缺点。


技术实现要素：

5.本发明实施方式提出一种前瞻性采集校正序列触发的磁共振成像方法和装置。
6.本发明实施方式的技术方案如下：
7.一种pace序列触发的磁共振成像方法，包括：
8.确定单激发快速自旋回波序列的重复时间和采样窗时间；
9.基于所述采样窗时间和所述重复时间确定每个生理运动周期内的最大成像层数n，其中n为大于等于2的正整数；
10.当pace序列生成触发信号时，使能所述单激发快速自旋回波序列在至少一个采样窗时间内获取m层成像数据，其中m为大于等于2且小于等于n的正整数。
11.可见，不同于现有技术仅在采样窗时间内获取单层成像图像，在本发明实施方式在采样窗时间内可以获取多层成像数据，从而提高扫描速度，并减少成像时间。
12.在一个实施方式中，所述确定采样窗时间包括：
13.施加所述pace序列以探测生理信号曲线；
14.基于所述生理信号曲线确定生理运动周期；
15.基于所述生理运动周期确定采样窗时间。
16.因此，本发明实施方式可以基于pace序列探测出的生理信号曲线确定采样窗时间，提高了采样窗时间的准确度。
17.在一个实施方式中，所述基于采样窗时间和重复时间确定每个生理运动周期内的最大成像层数n包括：
18.确定所述最大成像层数n，其中n＝int(ta/tr)；
19.其中tr为所述重复时间，ta为所述采样窗时间，int(ta/tr)为(ta/tr)的取整值。
20.可见，本发明实施方式可以基于采样窗时间和重复时间快速确定最大成像层数。
21.在一个实施方式中，所述使能单激发快速自旋回波序列在至少一个采样窗时间内获取m层成像数据包括：
22.当预定的总成像层数l不能被m整除时，确定生理运动周期的数目k，其中k＝(int(l/m) 1)，int(l/m)为(l/m)的取整值；使能单激发快速自旋回波序列在第一个生理运动周期到第(k-1)个生理运动周期中的每个生理运动周期的采样窗时间内分别获取m层成像数据；使能单激发快速自旋回波序列在第k个生理运动周期内获取(l-(k-1)*m)层成像数据；或
23.当预定的总成像层数l能被m整除时，确定生理运动周期的数目k，其中k＝(l/m)；使能单激发快速自旋回波序列在第一个生理运动周期到第k个生理运动周期中的每个生理运动周期的采样窗时间内分别获取m层成像数据。
24.因此，本发明实施方式考虑到总成像层数能否被m整除的不同状况，实现了在采样窗时间内采集多层成像数据的不同处理方式。
25.一种pace序列触发的磁共振成像装置，包括：
26.时间确定模块，用于确定单激发快速自旋回波序列的重复时间和采样窗时间；
27.最大成像层数确定模块，用于基于所述采样窗时间和所述重复时间确定每个生理运动周期内的最大成像层数n，其中n为大于等于2的正整数；
28.成像数据获取模块，用于当pace序列生成触发信号时，使能所述单激发快速自旋回波序列在至少一个采样窗时间内获取m层成像数据，其中m为大于等于2且小于等于n的正整数。
29.可见，不同于现有技术仅在采样窗时间内获取单层成像图像，在本发明实施方式在采样窗时间内可以获取多层成像数据，从而提高扫描速度，并减少成像时间。
30.在一个实施方式中，时间确定模块，用于施加所述pace序列以探测生理信号曲线；基于所述生理信号曲线确定生理运动周期；基于所述生理运动周期确定采样窗时间。
31.因此，本发明实施方式可以基于pace序列探测出的生理信号曲线确定采样窗时间，提高了采样窗时间的准确度。
32.在一个实施方式中，最大成像层数确定模块，用于确定所述最大成像层数n，其中n＝int(ta/tr)；其中tr为所述重复时间，ta为所述采样窗时间，int(ta/tr)为(ta/tr)的取整值。
33.可见，本发明实施方式可以基于采样窗时间和重复时间快速确定最大成像层数。
34.在一个实施方式中，成像数据获取模块，用于当预定的总成像层数l不能被m整除时，确定生理运动周期的数目k，其中k＝(int(l/m) 1)，int(l/m)为(l/m)的取整值；使能单激发快速自旋回波序列在第一个生理运动周期到第(k-1)个生理运动周期中的每个生理运动周期的采样窗时间内分别获取m层成像数据；使能单激发快速自旋回波序列在第k个生理运动周期内获取(l-(k-1)*m)层成像数据；或，当预定的总成像层数l能被m整除时，确定生理运动周期的数目k，其中k＝(l/m)；使能单激发快速自旋回波序列在第一个生理运动周期到第k个生理运动周期中的每个生理运动周期的采样窗时间内分别获取m层成像数据。
35.因此，本发明实施方式考虑到总成像层数能否被m整除的不同状况，实现了在采样窗时间内采集多层成像数据的不同处理方式。
36.一种磁共振成像系统的控制主机，包括：存储器；处理器；其中所述存储器中存储有可被所述处理器执行的应用程序，用于使得所述处理器执行如上任一项所述的pace序列触发的磁共振成像方法。
37.可见，本发明实施方式的控制主机，可以控制在采样窗时间内获取多层成像数据，从而提高扫描速度，并减少成像时间。
38.一种计算机可读存储介质，所述计算机可读存储介质上存储计算机程序，所述计算机程序被处理器执行时实现如上任一项所述的pace序列触发的磁共振成像方法。
39.因此，本发明实施方式还提出了一种计算机可读存储介质，可以控制在采样窗时间内获取多层成像数据，从而提高扫描速度，并减少成像时间。
附图说明
40.图1为现有技术中pace序列触发的haste序列获取成像数据的示意图。
41.图2为根据本发明实施方式的pace序列触发的磁共振成像方法的流程图。
42.图3为根据本发明实施方式的haste序列获取单层成像数据的示意图。
43.图4为根据本发明实施方式的pace序列触发的haste序列获取成像数据的示意图。
44.图5为根据本发明实施方式的pace序列触发的磁共振成像装置的结构图。
45.图6为根据本发明实施方式的磁共振成像系统的控制主机的结构图。
46.其中，附图标记如下：
47.101生理信号曲线102导航序列103单层成像数据104接收窗200pace序列触发的磁共振成像方法201～203步骤300haste序列获取单层成像数据的执行过程301射频激发脉冲302射频回聚脉冲401第一个单层成像数据402第二个单层成像数据
403第三个单层成像数据404第四个单层成像数据500pace序列触发的磁共振成像装置501时间确定模块502最大成像层数确定模块503成像数据获取模块600磁共振成像系统的控制主机601存储器602处理器
具体实施方式
48.为了使本发明的技术方案及优点更加清楚明白，以下结合附图及实施方式，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施方式仅仅用以阐述性说明本发明，并不用于限定本发明的保护范围。
49.为了描述上的简洁和直观，下文通过描述若干代表性的实施方式来对本发明的方案进行阐述。实施方式中大量的细节仅用于帮助理解本发明的方案。但是很明显，本发明的技术方案实现时可以不局限于这些细节。为了避免不必要地模糊了本发明的方案，一些实施方式没有进行细致地描述，而是仅给出了框架。下文中，“包括”是指“包括但不限于”，“根据
……”
是指“至少根据
……
，但不限于仅根据
……”
。由于汉语的语言习惯，下文中没有特别指出一个成分的数量时，意味着该成分可以是一个也可以是多个，或可理解为至少一个。
50.申请人发现：pace技术作为一种呼吸门控技术，可以有效去除呼吸伪影。当pace技术应用于呼吸门控时，不断扫描小角度激发梯度回波序列，采集低分辨率的二维图像，获取触发点以触发诸如单激发快速自旋回波序列之类的成像序列的扫描。在目前的单激发快速自旋回波序列中，当从pace序列接收到触发信号后，单激发快速自旋回波序列在单个生理运动周期内获取单层的成像数据。然而，单个生理运动周期只能获取单层的成像数据，因此扫描速度较慢，成像时间长。
51.以单激发快速自旋回波序列实施为haste序列为例进行说明。
52.图1为现有技术中pace序列触发的haste序列获取成像数据的示意图。由图1可见，在时间轴t上，包含多个生理运动周期(比如呼吸周期)的生理信号曲线101周期性延伸。在每个生理运动周期中，在pace序列发出触发信号之前，导航序列102被执行。当pace序列的接收窗104与生理信号曲线101接触时，pace序列向haste序列发出如箭头s所示的触发信号。在生理信号曲线101的两个相邻波谷(比如，点a与点b)所形成的一个生理运动周期中，当从pace序列接收到触发信号后，haste序列执行成像扫描以获取单层成像数据103。即，在每个呼吸周期内，haste序列只获取一层的成像数据。可见，生理运动周期内的其它时间段并没有被利用，因此扫描速度较慢，成像时间长。
53.申请人发现了现有技术中的上述缺陷。经过进一步研究，申请人还发现：如果不在单个生理运动周期内获取单层的成像数据，而是在单个生理运动周期内获取多层的成像数据，则可以显著提高扫描速度，并由此降低成像时间。
54.图2为根据本发明实施方式的pace序列触发的磁共振成像方法的流程图。
55.如图2所示，pace序列触发的磁共振成像方法200包括：
56.步骤201：确定单激发快速自旋回波序列的重复时间(repetition time)和采样窗时间(acquisition window time)。
57.在一个实施方式中，单激发快速自旋回波序列可以实施为haste序列、tse-ssh序列或ss-fse序列，等等。
58.单激发快速自旋回波序列的重复时间为一个激励脉冲的施加时刻与下一个激励脉冲的施加时刻之间的时间长度，它决定了这两个激励脉冲之间纵向磁化矢量的恢复程度。也就是，单激发快速自旋回波序列的重复时间为单激发快速自旋回波序列中各次重复射频脉冲之间的时间间隔。
59.可以由用户自行设置重复时间。比如，基于用户对成像图像的分辨率的设置，实现间接设置重复时间。通常情况下，设置的分辨率越高，则重复时间越长。
60.采样窗时间为单激发快速自旋回波序列在生理运动周期内的扫描(即数据获取)时间。每个生理运动周期都具有相同的采样窗时间。采样窗时间的起始点为单激发快速自旋回波序列从pace序列接收到触发信号的时刻，采样窗时间的结束点可以为生理运动周期的终止时刻之前的任意时刻。
61.在一个实施方式中，可以基于用户设置确定出适用于所有用户的采样窗时间。比如，考虑到大多数人群的生理运动周期通常为2秒～4秒之间，因此可以将采样窗时间设置为1～2秒之间。此时，该采样窗时间适用于所有用户。
62.在一个实施方式中，可以基于用户的个性化呼吸特征确定出适用于该用户的采样窗时间。比如，施加pace序列以探测该用户的生理信号曲线；基于该用户的生理信号曲线确定该用户的生理运动周期；基于该用户的生理运动周期确定该用户的采样窗时间。优选地，将生理运动周期的一半时间设置为该用户的采样窗时间。比如，假定pace序列检测出用户的生理运动周期为2.4秒，可以将该用户的采样窗时间设置为1.2秒。
63.以上示范性描述了单激发快速自旋回波序列、重复时间和采样窗时间的典型实例，本领域技术人员可以意识到，这种描述仅是示范性的，并不用于限定本发明实施方式的保护范围。
64.步骤202：基于采样窗时间和重复时间确定每个生理运动周期内的最大成像层数n，其中n为大于等于2的正整数。
65.在这里，最大成像层数n为在每个生理运动周期内所能采集到的成像层数的最大值。
66.在一个实施方式中，确定最大成像层数n，其中n＝int(ta/tr)；其中tr为重复时间，ta为采样窗时间，int(ta/tr)为(ta/tr)的取整值。
67.举例，假定重复时间tr为400毫秒(ms)，ta为1.3秒(s)，则ta/tr＝1.3(s)/400(ms)＝3.25，int(ta/tr)＝3，最大成像层数n等于3。
68.再举例，假定重复时间tr为500毫秒，ta为1.3秒，则ta/tr＝1.3(s)/500(ms)＝2.6，int(ta/tr)＝2，则最大成像层数n等于2。
69.步骤203：当pace序列生成触发信号时，使能单激发快速自旋回波序列在至少一个生理运动周期的采样窗时间内获取m层成像数据，其中m为大于等于2且小于等于n的正整数。m值可以由用户自行设置，也可以为磁共振成像系统的控制主机自行随机确定。
70.在一个实施方式中，当预定的总成像层数l不能被m整除时，确定生理运动周期的数目k，其中k＝(int(l/m) 1)，int(l/m)为(l/m)的取整值；使能单激发快速自旋回波序列在第一个生理运动周期到第(k-1)个生理运动周期中的每个生理运动周期的采样窗时间内分别获取m层成像数据；使能单激发快速自旋回波序列在第k个生理运动周期内获取(l-(k-1)*m)层成像数据。
71.当预定的总成像层数l不能被m整除时，m值可以为小于等于n的任意正整数，包括下列情形：
72.情形(一)：m小于n。
73.举例1：假定预定的总成像层数l为25，在步骤202中确定的最大成像层数n为3，则m值为小于3且大于等于2的正整数，具体为2。此时，确定生理运动周期的数目k＝(int(l/m) 1)＝int(25/2) 1＝13。因此，全部生理运动周期包括在时间轴上递增的第1个生理运动周期、第2个生理运动周期、第3个生理运动周期
…
第13个生理运动周期。其中：在第1个生理运动周期内，执行2次扫描；在第2个生理运动周期内，执行2次扫描；在第3个生理运动周期内，执行2次扫描；
…
在第12个生理运动周期内，执行2次扫描。而且，在第13个生理运动周期内，执行((l-(k-1)*m))＝1次扫描。
74.举例2：假定预定的总成像层数l为25，在步骤202中确定的最大成像层数n为4，则m值可以为小于4且大于等于2的正整数，比如为2或3，假定m为3。此时，确定生理运动周期的数目k＝int(l/m) 1＝int(25/3) 1＝9，即生理运动周期的数目为9个。因此，全部生理运动周期包括在时间轴上递增的第1个生理运动周期、第2个生理运动周期、第3个生理运动周期
…
第9个生理运动周期。其中：在第1个生理运动周期内，执行3次扫描；在第2个生理运动周期内，执行3次扫描；在第3个生理运动周期内，执行3次扫描；
…
在第8个生理运动周期内，执行3次扫描。而且，在第9个生理运动周期内，执行((l-(k-1)*m))＝1次扫描。
75.情形(二)：m等于n。
76.举例1：假定预定的总成像层数l为25，在步骤202中确定的最大成像层数n为3。m值等于n，具体为3。此时，确定生理运动周期的数目k＝int(l/m) 1＝int(25/3) 1＝9，即生理运动周期的数目为9个。因此，全部生理运动周期包括在时间轴上递增的第1个生理运动周期、第2个生理运动周期、第3个生理运动周期
…
第9个生理运动周期。其中：在第1个生理运动周期内，执行3次扫描；在第2个生理运动周期内，执行3次扫描；在第3个生理运动周期内，执行3次扫描；
…
在第8个生理运动周期内，执行3扫描。而且，在第9个生理运动周期内，执行((l-(k-1)*m))＝1次扫描。
77.举例2：假定预定的总成像层数l为30，在步骤202中确定的最大成像层数n为4。m值等于n，具体为4。此时，确定生理运动周期的数目k＝int(l/m) 1＝int(30/4) 1＝8，即生理运动周期的数目为8个。因此，全部生理运动周期包括在时间轴上递增的第1个生理运动周期、第2个生理运动周期、第3个生理运动周期
…
第8个生理运动周期。其中：在第1个生理运动周期内，执行4次扫描；在第2个生理运动周期内，执行4次扫描；在第3个生理运动周期内，执行4次扫描；
…
在第7个生理运动周期内，执行4次扫描。而且，在第8个生理运动周期内，执行((l-(k-1)*m))＝2次扫描。
78.在一个实施方式中，当预定的总成像层数l能被m整除时，确定生理运动周期的数目k，其中k＝(l/m)，使能单激发快速自旋回波序列在第一个生理运动周期到第k个生理运
动周期中的每个生理运动周期的采样窗时间内分别获取m层成像数据。
79.当预定的总成像层数l能被m整除时，m值可以为小于等于n的任意正整数，包括下列情形：
80.情形(一)：m小于n。
81.举例1：假定预定的总成像层数l为24，在步骤202中确定的最大成像层数n为3，则m值为小于3且大于等于2的正整数，具体为2。此时，确定生理运动周期的数目k＝24/2＝12，即生理运动周期的数目为12个。因此，全部生理运动周期包括在时间轴上递增的第1个生理运动周期、第2个生理运动周期、第3个生理运动周期
…
第12个生理运动周期。其中：在第1个生理运动周期内，执行2次扫描；在第2个生理运动周期内，执行2次扫描；在第3个生理运动周期内，执行2次扫描；
…
在第12个生理运动周期内，执行2次扫描。
82.举例2：假定预定的总成像层数l为24，在步骤202中确定的最大成像层数n为4，则m值可以为小于4且大于等于2的正整数，比如为2或3，假定m为3。此时，确定生理运动周期的数目k＝24/3＝8，即生理运动周期的数目为8个。因此，全部生理运动周期包括在时间轴上递增的第1个生理运动周期、第2个生理运动周期、第3个生理运动周期
…
第8个生理运动周期。其中：在第1个生理运动周期内，执行3次扫描；在第2个生理运动周期内，执行3次扫描；在第3个生理运动周期内，执行3次扫描；
…
在第8个生理运动周期内，执行3次扫描。
83.情形(二)：m等于n。
84.举例1：假定预定的总成像层数l为30，在步骤202中确定的最大成像层数n为3。m值等于n，具体为3。此时，确定生理运动周期的数目k＝30/3＝10，即生理运动周期的数目为10个。因此，全部生理运动周期包括在时间轴上递增的第1个生理运动周期、第2个生理运动周期、第3个生理运动周期
…
第10个生理运动周期。其中：在第1个生理运动周期内，执行3次扫描；在第2个生理运动周期内，执行3次扫描；在第3个生理运动周期内，执行3次扫描；
…
在第10个生理运动周期内，执行3次扫描。
85.举例2：假定预定的总成像层数l为28，在步骤202中确定的最大成像层数n为4。m值等于n，具体为4。此时，确定生理运动周期的数目k＝28/4＝7，即生理运动周期的数目为7个。因此，全部生理运动周期包括在时间轴上递增的第1个生理运动周期、第2个生理运动周期、第3个生理运动周期
…
第7个生理运动周期。其中：在第1个生理运动周期内，执行4次扫描；在第2个生理运动周期内，执行4次扫描；在第3个生理运动周期内，执行4次扫描；
…
在第7个生理运动周期内，执行4次扫描。
86.以上示范性描述了在采样窗时间内执行m次扫描的具体方式，本领域技术人员可以意识到，这种描述仅是示范性的，并不用于限定本发明实施方式的保护范围。
87.优选地，同一生理运动周期内所获取的m层成像数据的层号相互错开(优选为等间隔错开)，从而防止同一生理运动周期内的数据干扰。比如，假定生理运动周期共为4个，且在每个生理运动周期内的采样窗时间内执行3次扫描。那么，第1个生理运动周期内，可以分别采集第1层成像数据、第5层成像数据和第9层成像数据；第2个生理运动周期内，可以分别采集第2层成像数据、第6层成像数据和第10层成像数据；第3个生理运动周期内，可以分别采集第3层成像数据、第7层成像数据和第11层成像数据；第4个生理运动周期内，可以分别采集第4层成像数据、第8层成像数据和第12层成像数据。
88.以上示范性描述了m层成像数据的层号相互错开的典型实例，本领域技术人员可
以意识到，这种描述仅是示范性的，并不用于限定本发明实施方式的保护范围。
89.下面以单激发快速自旋回波序列具体实施为haste序列为实例，对本发明实施方式进行说明。
90.图3为根据本发明实施方式的haste序列获取单层成像数据的示意图。
91.在图3中，从上到下分别标出了射频脉冲(rf)，选层梯度(gs)，相位编码梯度(gp)，读出梯度(gro)和模数转换(adc)在时间轴上的执行顺序。其中，adc即为数据采集步骤。
92.第一步：在rf上施加90度的射频激发脉冲301，同时在gs方向施加选层梯度，以便将所选层内的磁化强度矢量翻转到横向平面。
93.第二步：紧接着在gs方向施加补偿梯度，以补偿从激发脉冲中心到脉冲结束这段时间内选层梯度引起的磁化强度矢量的散相，同时在gro方向施加预散相梯度。
94.第三步：在rf上施加180度的射频回聚脉冲302，同时在gs方向施加选层梯度，以便使所选层内的磁化强度矢量在随后的两个180度射频回聚脉冲的中间回聚形成回波。紧接着在gp方向施加相位编码梯度。
95.第四步：在gro方向施加读出梯度同时在adc上进行数据采集，紧接着再在gs方向施加与前述相位编码梯度方向相反的梯度，以便将前一个相位编码梯度的效应抵消。
96.第五步：根据预先设置的在相位编码方向的分辨率，重复执行第三步和第四步直到一层完整的数据采集完毕。随后重复第一步到第五步进行下一层数据的采集。
97.可见，haste序列是一种仅激发一次(一个90度激发脉冲)，就可以获得一层完整图像的单激发快速成像序列。
98.其中，当成像图像的分辨率为256*256时，射频回聚脉冲302的数目为136个；当成像图像的分辨率为128*128时，射频回聚脉冲302的数目为72个；当成像图像的分辨率为320*320时，射频回聚脉冲302的数目为184个。可见，通过设置不同的分辨率，可以改变射频回聚脉冲302的数目，从而改变haste序列的重复时间。
99.图4为根据本发明实施方式的pace序列触发的haste序列获取成像数据的示意图。
100.由图4可见，在时间轴t上，包含多个生理运动周期(比如呼吸周期)的生理信号曲线101周期性延伸。在每个生理运动周期中，在pace序列发出触发信号之前，导航序列102被执行。当pace序列的接收窗104与生理信号曲线101接触时，pace序列向haste序列发出如箭头s所示的触发信号。假定基于重复时间tr和采样窗时间ta确定的最大成像层数n为4，且m等于n。
101.在生理信号曲线101的两个相邻波谷(比如，点a与点b)所形成的一个生理运动周期(比如呼吸周期)中，当从pace序列接收到触发信号后，haste序列执行成像扫描以获取第一个单层成像数据401、第二个单层成像数据402、第三个单层成像数据403和第四个单层成像数据404。也就是说，当接收到触发信号之后，haste序列在一个生理运动周期之内执行4次如图3所示的haste序列获取单层成像数据的执行过程300，以获取4层成像数据。
102.另外，由图4可见，ta的起始点可以为haste序列从pace序列接收到触发信号的时刻，ta的结束点可以为生理运动周期的终止时刻之前的任意时刻。
103.在图3-图4中以单激发快速成像序列实施为haste序列为例进行说明，本领域技术人员可以意识到，单激发快速成像序列还可以实施为tse-ssh序列或ss-fse序列，等等。
104.基于上述描述，本发明实施方式还提出了pace序列触发的磁共振成像装置。
105.图5为根据本发明实施方式的pace序列触发的磁共振成像装置的结构图。
106.如图5所示，pace序列触发的磁共振成像装置500包括：
107.时间确定模块501，用于确定单激发快速自旋回波序列的重复时间和采样窗时间；
108.最大成像层数确定模块502，用于基于所述采样窗时间和所述重复时间确定每个生理运动周期内的最大成像层数n，其中n为大于等于2的正整数；
109.成像数据获取模块503，用于当从pace序列接收到触发信号时，使能所述单激发快速自旋回波序列在至少一个采样窗时间内获取m层成像数据，其中m为大于等于2且小于等于n的正整数。
110.在一个实施方式中，时间确定模块501，用于施加pace序列以探测生理信号曲线；基于所述生理信号曲线确定生理运动周期；基于所述生理运动周期确定采样窗时间。
111.在一个实施方式中，最大成像层数确定模块502，用于确定所述最大成像层数n，其中n＝int(ta/tr)；其中tr为所述重复时间，ta为所述采样窗时间，int(ta/tr)为(ta/tr)的取整值。
112.在一个实施方式中，成像数据获取模块503，用于当预定的总成像层数l不能被m整除时，确定生理运动周期的数目k，其中k＝(int(l/m) 1)，int(l/m)为(l/m)的取整值；使能单激发快速自旋回波序列在第一个生理运动周期到第(k-1)个生理运动周期中的每个生理运动周期的采样窗时间内分别获取m层成像数据；使能单激发快速自旋回波序列在第k个生理运动周期内获取(l-(k-1)*m)层成像数据；或，当预定的总成像层数l能被m整除时，确定生理运动周期的数目k，其中k＝(l/m)；使能单激发快速自旋回波序列在第一个生理运动周期到第k个生理运动周期中的每个生理运动周期的采样窗时间内分别获取m层成像数据。
113.基于上述描述，本发明实施方式还提出了磁共振成像系统的控制主机。
114.图6为根据本发明实施方式的磁共振成像系统的控制主机的结构图。
115.如图6所示，控制主机600包括处理器601、存储器602及存储在存储器602上并可在处理器601上运行的计算机程序，计算机程序被处理器601执行时实现如上任一项的pace序列触发的磁共振成像方法。
116.其中，存储器602具体可以实施为电可擦可编程只读存储器(eeprom)、快闪存储器(flash memory)、可编程程序只读存储器(prom)等多种存储介质。处理器601可以实施为包括一或多个中央处理器或一或多个现场可编程门阵列，其中现场可编程门阵列集成一或多个中央处理器核。具体地，中央处理器或中央处理器核可以实施为cpu或mcu或dsp等等。
117.综上所述，在本发明实施方式中，确定单激发快速自旋回波序列的重复时间和采样窗时间；基于采样窗时间和重复时间确定每个生理运动周期内的最大成像层数n，其中n为大于等于2的正整数；当pace序列生成触发信号时，使能单激发快速自旋回波序列在至少一个采样窗时间内获取m层成像数据，其中m为大于等于2且小于等于n的正整数。由此可见，本发明实施方式在单个采样窗时间内获取多层成像数据，可以提高扫描速度，并减少成像时间。
118.而且，本发明实施方式可以基于pace序列探测出的生理信号曲线确定采样窗时间，提高了采样窗时间的准确度。
119.另外，本发明实施方式考虑到总成像层数能否被m整除的不同状况，实现了在采样窗时间内采集多层成像数据的不同处理方式。
120.可以遵循一定规范的应用程序接口，将本发明实施方式所提出的前瞻性心电触发方法编写为安装到磁共振成像控制主机、个人电脑、移动终端等中的插件程序，也可以将其封装为应用程序以供用户自行下载使用。
121.可以通过指令或指令集存储的储存方式将本发明实施方式所提出的前瞻性心电触发方法存储在各种存储介质上。这些存储介质包括但是不局限于：软盘、光盘、dvd、硬盘、闪存等。另外，还可以将本发明实施方式所提出的前瞻性心电触发方法应用到基于闪存(nand flash)的存储介质中，比如u盘、cf卡、sd卡、sdhc卡、mmc卡、sm卡、记忆棒、xd卡等。
122.以上所述，仅为本发明的较佳实施例而已，并非用于限定本发明的保护范围。凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。