CN114076913B - 一种前瞻性采集校正序列触发的磁共振成像方法和装置 - Google Patents

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

Info

Publication number
CN114076913B
CN114076913B CN202010836538.1A CN202010836538A CN114076913B CN 114076913 B CN114076913 B CN 114076913B CN 202010836538 A CN202010836538 A CN 202010836538A CN 114076913 B CN114076913 B CN 114076913B
Authority
CN
China
Prior art keywords
physiological motion
physiological
sampling window
imaging
sequence
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202010836538.1A
Other languages
English (en)
Other versions
CN114076913A (zh
Inventor
王玉宇
董芳
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Siemens Shenzhen Magnetic Resonance Ltd
Original Assignee
Siemens Shenzhen Magnetic Resonance Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Siemens Shenzhen Magnetic Resonance Ltd filed Critical Siemens Shenzhen Magnetic Resonance Ltd
Priority to CN202010836538.1A priority Critical patent/CN114076913B/zh
Priority to US17/405,475 priority patent/US20220057466A1/en
Publication of CN114076913A publication Critical patent/CN114076913A/zh
Application granted granted Critical
Publication of CN114076913B publication Critical patent/CN114076913B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00Arrangements or instruments for measuring magnetic variables
    • G01R33/20Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48NMR imaging systems
    • G01R33/54Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
    • G01R33/56Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
    • G01R33/567Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution gated by physiological signals, i.e. synchronization of acquired MR data with periodical motion of an object of interest, e.g. monitoring or triggering system for cardiac or respiratory gating
    • G01R33/5676Gating or triggering based on an MR signal, e.g. involving one or more navigator echoes for motion monitoring and correction
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00Arrangements or instruments for measuring magnetic variables
    • G01R33/20Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48NMR imaging systems
    • G01R33/54Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
    • G01R33/56Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
    • G01R33/565Correction of image distortions, e.g. due to magnetic field inhomogeneities
    • G01R33/56509Correction of image distortions, e.g. due to magnetic field inhomogeneities due to motion, displacement or flow, e.g. gradient moment nulling
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/05Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves 
    • A61B5/055Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves  involving electronic [EMR] or nuclear [NMR] magnetic resonance, e.g. magnetic resonance imaging
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00Arrangements or instruments for measuring magnetic variables
    • G01R33/20Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48NMR imaging systems
    • G01R33/54Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
    • G01R33/56Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
    • G01R33/561Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution by reduction of the scanning time, i.e. fast acquiring systems, e.g. using echo-planar pulse sequences
    • G01R33/5615Echo train techniques involving acquiring plural, differently encoded, echo signals after one RF excitation, e.g. using gradient refocusing in echo planar imaging [EPI], RF refocusing in rapid acquisition with relaxation enhancement [RARE] or using both RF and gradient refocusing in gradient and spin echo imaging [GRASE]
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00Arrangements or instruments for measuring magnetic variables
    • G01R33/20Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48NMR imaging systems
    • G01R33/54Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
    • G01R33/56Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
    • G01R33/561Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution by reduction of the scanning time, i.e. fast acquiring systems, e.g. using echo-planar pulse sequences
    • G01R33/5615Echo train techniques involving acquiring plural, differently encoded, echo signals after one RF excitation, e.g. using gradient refocusing in echo planar imaging [EPI], RF refocusing in rapid acquisition with relaxation enhancement [RARE] or using both RF and gradient refocusing in gradient and spin echo imaging [GRASE]
    • G01R33/5617Echo train techniques involving acquiring plural, differently encoded, echo signals after one RF excitation, e.g. using gradient refocusing in echo planar imaging [EPI], RF refocusing in rapid acquisition with relaxation enhancement [RARE] or using both RF and gradient refocusing in gradient and spin echo imaging [GRASE] using RF refocusing, e.g. RARE
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00Arrangements or instruments for measuring magnetic variables
    • G01R33/20Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48NMR imaging systems
    • G01R33/54Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
    • G01R33/56Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
    • G01R33/567Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution gated by physiological signals, i.e. synchronization of acquired MR data with periodical motion of an object of interest, e.g. monitoring or triggering system for cardiac or respiratory gating
    • G01R33/5673Gating or triggering based on a physiological signal other than an MR signal, e.g. ECG gating or motion monitoring using optical systems for monitoring the motion of a fiducial marker

Landscapes

  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • General Health & Medical Sciences (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Radiology & Medical Imaging (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Condensed Matter Physics & Semiconductors (AREA)
  • Signal Processing (AREA)
  • Biophysics (AREA)
  • Physiology (AREA)
  • Pulmonology (AREA)
  • Power Engineering (AREA)
  • Cardiology (AREA)
  • Medical Informatics (AREA)
  • Veterinary Medicine (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Biomedical Technology (AREA)
  • Pathology (AREA)
  • Public Health (AREA)
  • Molecular Biology (AREA)
  • Animal Behavior & Ethology (AREA)
  • Surgery (AREA)
  • Magnetic Resonance Imaging Apparatus (AREA)

Abstract

本发明实施方式公开了一种前瞻性采集校正序列触发的磁共振成像方法和装置。方法包括:确定单激发快速自旋回波序列的重复时间和采样窗时间;基于所述采样窗时间和所述重复时间确定每个生理运动周期内的最大成像层数N,其中N为大于等于2的正整数;当前瞻性采集校正序列生成触发信号时,使能单激发快速自旋回波序列在至少一个采样窗时间内获取M层成像数据,其中M为大于等于2且小于等于N的正整数。本发明实施方式在单个采样窗时间内获取多层成像数据,可以提高扫描速度,并减少成像时间。

Description

一种前瞻性采集校正序列触发的磁共振成像方法和装置
技术领域
本发明涉及医学器械技术领域,特别是涉及一种前瞻性采集校正(ProspectiveAcquisition CorrEction,PACE)序列触发的磁共振成像方法和装置。
背景技术
磁共振成像(Magnetic Resonance Imaging,MRI)是随着计算机技术、电子电路技术、超导体技术的发展而迅速发展起来的一种生物磁学核自旋成像技术。它利用磁场与射频脉冲使人体组织内进动的氢核(即H+)发生振动产生射频信号,经计算机处理而成像。当把物体放置在磁场中,用适当的电磁波照射它,使之共振,然后分析它释放的电磁波,就可以得知构成这一物体的原子核的位置和种类,据此可以绘制成物体内部的精确立体影像。比如,可以通过磁共振成像扫描人类大脑获得的一个连续切片的动画,由头顶开始,一直到脚部。
单激发快速自旋回波序列是一种高速磁共振成像序列,适用于许多临床应用,尤其是运动敏感成像应用。仅激发一次单激发快速自旋回波序列,就可以获得成像对象的一层完整图像。单激发快速自旋回波序列通常包括半傅里叶采集单次激发快速自旋回波序列(half-fourier acquisition single-shot turbo spin-echo,HASTE)、单次激发TSE(Turbo spin echo-single-shot,TSE-SSH)序列、单次激发FSE序列(single shot fastspin echo,SS-FSE)序列,等等。在对胸腹部进行磁共振成像时,由于生理运动的影响(例如患者的呼吸运动),图像中经常出现伪影。因此,需要对检查对象的生理运动进行监测,从而消除由此产生的伪影。通过结合PACE序列的触发机制,单激发快速自旋回波序列可以获取几乎没有运动伪影的成像图像。
在目前的单激发快速自旋回波序列中,当从PACE序列接收到触发信号后,单激发快速自旋回波序列在每个生理运动周期内分别获取单层的成像数据,具有成像时间长的缺点。
发明内容
本发明实施方式提出一种前瞻性采集校正序列触发的磁共振成像方法和装置。
本发明实施方式的技术方案如下:
一种PACE序列触发的磁共振成像方法,包括:
确定单激发快速自旋回波序列的重复时间和采样窗时间;
基于所述采样窗时间和所述重复时间确定每个生理运动周期内的最大成像层数N,其中N为大于等于2的正整数;
当PACE序列生成触发信号时,使能所述单激发快速自旋回波序列在至少一个采样窗时间内获取M层成像数据,其中M为大于等于2且小于等于N的正整数。
可见,不同于现有技术仅在采样窗时间内获取单层成像图像,在本发明实施方式在采样窗时间内可以获取多层成像数据,从而提高扫描速度,并减少成像时间。
在一个实施方式中,所述确定采样窗时间包括:
施加所述PACE序列以探测生理信号曲线;
基于所述生理信号曲线确定生理运动周期;
基于所述生理运动周期确定采样窗时间。
因此,本发明实施方式可以基于PACE序列探测出的生理信号曲线确定采样窗时间,提高了采样窗时间的准确度。
在一个实施方式中,所述基于采样窗时间和重复时间确定每个生理运动周期内的最大成像层数N包括:
确定所述最大成像层数N,其中N=int(TA/TR);
其中TR为所述重复时间,TA为所述采样窗时间,int(TA/TR)为(TA/TR)的取整值。
可见,本发明实施方式可以基于采样窗时间和重复时间快速确定最大成像层数。
在一个实施方式中,所述使能单激发快速自旋回波序列在至少一个采样窗时间内获取M层成像数据包括:
当预定的总成像层数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层成像数据。
因此,本发明实施方式考虑到总成像层数能否被M整除的不同状况,实现了在采样窗时间内采集多层成像数据的不同处理方式。
一种PACE序列触发的磁共振成像装置,包括:
时间确定模块,用于确定单激发快速自旋回波序列的重复时间和采样窗时间;
最大成像层数确定模块,用于基于所述采样窗时间和所述重复时间确定每个生理运动周期内的最大成像层数N,其中N为大于等于2的正整数;
成像数据获取模块,用于当PACE序列生成触发信号时,使能所述单激发快速自旋回波序列在至少一个采样窗时间内获取M层成像数据,其中M为大于等于2且小于等于N的正整数。
可见,不同于现有技术仅在采样窗时间内获取单层成像图像,在本发明实施方式在采样窗时间内可以获取多层成像数据,从而提高扫描速度,并减少成像时间。
在一个实施方式中,时间确定模块,用于施加所述PACE序列以探测生理信号曲线;基于所述生理信号曲线确定生理运动周期;基于所述生理运动周期确定采样窗时间。
因此,本发明实施方式可以基于PACE序列探测出的生理信号曲线确定采样窗时间,提高了采样窗时间的准确度。
在一个实施方式中,最大成像层数确定模块,用于确定所述最大成像层数N,其中N=int(TA/TR);其中TR为所述重复时间,TA为所述采样窗时间,int(TA/TR)为(TA/TR)的取整值。
可见,本发明实施方式可以基于采样窗时间和重复时间快速确定最大成像层数。
在一个实施方式中,成像数据获取模块,用于当预定的总成像层数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层成像数据。
因此,本发明实施方式考虑到总成像层数能否被M整除的不同状况,实现了在采样窗时间内采集多层成像数据的不同处理方式。
一种磁共振成像系统的控制主机,包括:存储器;处理器;其中所述存储器中存储有可被所述处理器执行的应用程序,用于使得所述处理器执行如上任一项所述的PACE序列触发的磁共振成像方法。
可见,本发明实施方式的控制主机,可以控制在采样窗时间内获取多层成像数据,从而提高扫描速度,并减少成像时间。
一种计算机可读存储介质,所述计算机可读存储介质上存储计算机程序,所述计算机程序被处理器执行时实现如上任一项所述的PACE序列触发的磁共振成像方法。
因此,本发明实施方式还提出了一种计算机可读存储介质,可以控制在采样窗时间内获取多层成像数据,从而提高扫描速度,并减少成像时间。
附图说明
图1为现有技术中PACE序列触发的HASTE序列获取成像数据的示意图。
图2为根据本发明实施方式的PACE序列触发的磁共振成像方法的流程图。
图3为根据本发明实施方式的HASTE序列获取单层成像数据的示意图。
图4为根据本发明实施方式的PACE序列触发的HASTE序列获取成像数据的示意图。
图5为根据本发明实施方式的PACE序列触发的磁共振成像装置的结构图。
图6为根据本发明实施方式的磁共振成像系统的控制主机的结构图。
其中,附图标记如下:
101 生理信号曲线
102 导航序列
103 单层成像数据
104 接收窗
200 PACE序列触发的磁共振成像方法
201~203 步骤
300 HASTE序列获取单层成像数据的执行过程
301 射频激发脉冲
302 射频回聚脉冲
401 第一个单层成像数据
402 第二个单层成像数据
403 第三个单层成像数据
404 第四个单层成像数据
500 PACE序列触发的磁共振成像装置
501 时间确定模块
502 最大成像层数确定模块
503 成像数据获取模块
600 磁共振成像系统的控制主机
601 存储器
602 处理器
具体实施方式
为了使本发明的技术方案及优点更加清楚明白,以下结合附图及实施方式,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施方式仅仅用以阐述性说明本发明,并不用于限定本发明的保护范围。
为了描述上的简洁和直观,下文通过描述若干代表性的实施方式来对本发明的方案进行阐述。实施方式中大量的细节仅用于帮助理解本发明的方案。但是很明显,本发明的技术方案实现时可以不局限于这些细节。为了避免不必要地模糊了本发明的方案,一些实施方式没有进行细致地描述,而是仅给出了框架。下文中,“包括”是指“包括但不限于”,“根据……”是指“至少根据……,但不限于仅根据……”。由于汉语的语言习惯,下文中没有特别指出一个成分的数量时,意味着该成分可以是一个也可以是多个,或可理解为至少一个。
申请人发现:PACE技术作为一种呼吸门控技术,可以有效去除呼吸伪影。当PACE技术应用于呼吸门控时,不断扫描小角度激发梯度回波序列,采集低分辨率的二维图像,获取触发点以触发诸如单激发快速自旋回波序列之类的成像序列的扫描。在目前的单激发快速自旋回波序列中,当从PACE序列接收到触发信号后,单激发快速自旋回波序列在单个生理运动周期内获取单层的成像数据。然而,单个生理运动周期只能获取单层的成像数据,因此扫描速度较慢,成像时间长。
以单激发快速自旋回波序列实施为HASTE序列为例进行说明。
图1为现有技术中PACE序列触发的HASTE序列获取成像数据的示意图。由图1可见,在时间轴T上,包含多个生理运动周期(比如呼吸周期)的生理信号曲线101周期性延伸。在每个生理运动周期中,在PACE序列发出触发信号之前,导航序列102被执行。当PACE序列的接收窗104与生理信号曲线101接触时,PACE序列向HASTE序列发出如箭头S所示的触发信号。在生理信号曲线101的两个相邻波谷(比如,点A与点B)所形成的一个生理运动周期中,当从PACE序列接收到触发信号后,HASTE序列执行成像扫描以获取单层成像数据103。即,在每个呼吸周期内,HASTE序列只获取一层的成像数据。可见,生理运动周期内的其它时间段并没有被利用,因此扫描速度较慢,成像时间长。
申请人发现了现有技术中的上述缺陷。经过进一步研究,申请人还发现:如果不在单个生理运动周期内获取单层的成像数据,而是在单个生理运动周期内获取多层的成像数据,则可以显著提高扫描速度,并由此降低成像时间。
图2为根据本发明实施方式的PACE序列触发的磁共振成像方法的流程图。
如图2所示,PACE序列触发的磁共振成像方法200包括:
步骤201:确定单激发快速自旋回波序列的重复时间(repetition time)和采样窗时间(acquisition window time)。
在一个实施方式中,单激发快速自旋回波序列可以实施为HASTE序列、TSE-SSH序列或SS-FSE序列,等等。
单激发快速自旋回波序列的重复时间为一个激励脉冲的施加时刻与下一个激励脉冲的施加时刻之间的时间长度,它决定了这两个激励脉冲之间纵向磁化矢量的恢复程度。也就是,单激发快速自旋回波序列的重复时间为单激发快速自旋回波序列中各次重复射频脉冲之间的时间间隔。
可以由用户自行设置重复时间。比如,基于用户对成像图像的分辨率的设置,实现间接设置重复时间。通常情况下,设置的分辨率越高,则重复时间越长。
采样窗时间为单激发快速自旋回波序列在生理运动周期内的扫描(即数据获取)时间。每个生理运动周期都具有相同的采样窗时间。采样窗时间的起始点为单激发快速自旋回波序列从PACE序列接收到触发信号的时刻,采样窗时间的结束点可以为生理运动周期的终止时刻之前的任意时刻。
在一个实施方式中,可以基于用户设置确定出适用于所有用户的采样窗时间。比如,考虑到大多数人群的生理运动周期通常为2秒~4秒之间,因此可以将采样窗时间设置为1~2秒之间。此时,该采样窗时间适用于所有用户。
在一个实施方式中,可以基于用户的个性化呼吸特征确定出适用于该用户的采样窗时间。比如,施加PACE序列以探测该用户的生理信号曲线;基于该用户的生理信号曲线确定该用户的生理运动周期;基于该用户的生理运动周期确定该用户的采样窗时间。优选地,将生理运动周期的一半时间设置为该用户的采样窗时间。比如,假定PACE序列检测出用户的生理运动周期为2.4秒,可以将该用户的采样窗时间设置为1.2秒。
以上示范性描述了单激发快速自旋回波序列、重复时间和采样窗时间的典型实例,本领域技术人员可以意识到,这种描述仅是示范性的,并不用于限定本发明实施方式的保护范围。
步骤202:基于采样窗时间和重复时间确定每个生理运动周期内的最大成像层数N,其中N为大于等于2的正整数。
在这里,最大成像层数N为在每个生理运动周期内所能采集到的成像层数的最大值。
在一个实施方式中,确定最大成像层数N,其中N=int(TA/TR);其中TR为重复时间,TA为采样窗时间,int(TA/TR)为(TA/TR)的取整值。
举例,假定重复时间TR为400毫秒(ms),TA为1.3秒(s),则TA/TR=1.3(s)/400(ms)=3.25,int(TA/TR)=3,最大成像层数N等于3。
再举例,假定重复时间TR为500毫秒,TA为1.3秒,则TA/TR=1.3(s)/500(ms)=2.6,int(TA/TR)=2,则最大成像层数N等于2。
步骤203:当PACE序列生成触发信号时,使能单激发快速自旋回波序列在至少一个生理运动周期的采样窗时间内获取M层成像数据,其中M为大于等于2且小于等于N的正整数。M值可以由用户自行设置,也可以为磁共振成像系统的控制主机自行随机确定。
在一个实施方式中,当预定的总成像层数L不能被M整除时,确定生理运动周期的数目K,其中K=(int(L/M)+1),int(L/M)为(L/M)的取整值;使能单激发快速自旋回波序列在第一个生理运动周期到第(K-1)个生理运动周期中的每个生理运动周期的采样窗时间内分别获取M层成像数据;使能单激发快速自旋回波序列在第K个生理运动周期内获取(L-(K-1)*M)层成像数据。
当预定的总成像层数L不能被M整除时,M值可以为小于等于N的任意正整数,包括下列情形:
情形(一):M小于N。
举例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次扫描。
举例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次扫描。
情形(二):M等于N。
举例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次扫描。
举例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次扫描。
在一个实施方式中,当预定的总成像层数L能被M整除时,确定生理运动周期的数目K,其中K=(L/M),使能单激发快速自旋回波序列在第一个生理运动周期到第K个生理运动周期中的每个生理运动周期的采样窗时间内分别获取M层成像数据。
当预定的总成像层数L能被M整除时,M值可以为小于等于N的任意正整数,包括下列情形:
情形(一):M小于N。
举例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次扫描。
举例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次扫描。
情形(二):M等于N。
举例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次扫描。
举例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次扫描。
以上示范性描述了在采样窗时间内执行M次扫描的具体方式,本领域技术人员可以意识到,这种描述仅是示范性的,并不用于限定本发明实施方式的保护范围。
优选地,同一生理运动周期内所获取的M层成像数据的层号相互错开(优选为等间隔错开),从而防止同一生理运动周期内的数据干扰。比如,假定生理运动周期共为4个,且在每个生理运动周期内的采样窗时间内执行3次扫描。那么,第1个生理运动周期内,可以分别采集第1层成像数据、第5层成像数据和第9层成像数据;第2个生理运动周期内,可以分别采集第2层成像数据、第6层成像数据和第10层成像数据;第3个生理运动周期内,可以分别采集第3层成像数据、第7层成像数据和第11层成像数据;第4个生理运动周期内,可以分别采集第4层成像数据、第8层成像数据和第12层成像数据。
以上示范性描述了M层成像数据的层号相互错开的典型实例,本领域技术人员可以意识到,这种描述仅是示范性的,并不用于限定本发明实施方式的保护范围。
下面以单激发快速自旋回波序列具体实施为HASTE序列为实例,对本发明实施方式进行说明。
图3为根据本发明实施方式的HASTE序列获取单层成像数据的示意图。
在图3中,从上到下分别标出了射频脉冲(RF),选层梯度(Gs),相位编码梯度(Gp),读出梯度(Gro)和模数转换(ADC)在时间轴上的执行顺序。其中,ADC即为数据采集步骤。
第一步:在RF上施加90度的射频激发脉冲301,同时在Gs方向施加选层梯度,以便将所选层内的磁化强度矢量翻转到横向平面。
第二步:紧接着在Gs方向施加补偿梯度,以补偿从激发脉冲中心到脉冲结束这段时间内选层梯度引起的磁化强度矢量的散相,同时在Gro方向施加预散相梯度。
第三步:在RF上施加180度的射频回聚脉冲302,同时在Gs方向施加选层梯度,以便使所选层内的磁化强度矢量在随后的两个180度射频回聚脉冲的中间回聚形成回波。紧接着在Gp方向施加相位编码梯度。
第四步:在Gro方向施加读出梯度同时在ADC上进行数据采集,紧接着再在Gp方向施加与前述相位编码梯度方向相反的梯度,以便将前一个相位编码梯度的效应抵消。
第五步:根据预先设置的在相位编码方向的分辨率,重复执行第三步和第四步直到一层完整的数据采集完毕。随后重复第一步到第五步进行下一层数据的采集。
可见,HASTE序列是一种仅激发一次(一个90度激发脉冲),就可以获得一层完整图像的单激发快速成像序列。
其中,当成像图像的分辨率为256*256时,射频回聚脉冲302的数目为136个;当成像图像的分辨率为128*128时,射频回聚脉冲302的数目为72个;当成像图像的分辨率为320*320时,射频回聚脉冲302的数目为184个。可见,通过设置不同的分辨率,可以改变射频回聚脉冲302的数目,从而改变HASTE序列的重复时间。
图4为根据本发明实施方式的PACE序列触发的HASTE序列获取成像数据的示意图。
由图4可见,在时间轴T上,包含多个生理运动周期(比如呼吸周期)的生理信号曲线101周期性延伸。在每个生理运动周期中,在PACE序列发出触发信号之前,导航序列102被执行。当PACE序列的接收窗104与生理信号曲线101接触时,PACE序列向HASTE序列发出如箭头S所示的触发信号。假定基于重复时间TR和采样窗时间TA确定的最大成像层数N为4,且M等于N。
在生理信号曲线101的两个相邻波谷(比如,点A与点B)所形成的一个生理运动周期(比如呼吸周期)中,当从PACE序列接收到触发信号后,HASTE序列执行成像扫描以获取第一个单层成像数据401、第二个单层成像数据402、第三个单层成像数据403和第四个单层成像数据404。也就是说,当接收到触发信号之后,HASTE序列在一个生理运动周期之内执行4次如图3所示的HASTE序列获取单层成像数据的执行过程300,以获取4层成像数据。
另外,由图4可见,TA的起始点可以为HASTE序列从PACE序列接收到触发信号的时刻,TA的结束点可以为生理运动周期的终止时刻之前的任意时刻。
在图3-图4中以单激发快速成像序列实施为HASTE序列为例进行说明,本领域技术人员可以意识到,单激发快速成像序列还可以实施为TSE-SSH序列或SS-FSE序列,等等。
基于上述描述,本发明实施方式还提出了PACE序列触发的磁共振成像装置。
图5为根据本发明实施方式的PACE序列触发的磁共振成像装置的结构图。
如图5所示,PACE序列触发的磁共振成像装置500包括:
时间确定模块501,用于确定单激发快速自旋回波序列的重复时间和采样窗时间;
最大成像层数确定模块502,用于基于所述采样窗时间和所述重复时间确定每个生理运动周期内的最大成像层数N,其中N为大于等于2的正整数;
成像数据获取模块503,用于当从PACE序列接收到触发信号时,使能所述单激发快速自旋回波序列在至少一个采样窗时间内获取M层成像数据,其中M为大于等于2且小于等于N的正整数。
在一个实施方式中,时间确定模块501,用于施加PACE序列以探测生理信号曲线;基于所述生理信号曲线确定生理运动周期;基于所述生理运动周期确定采样窗时间。
在一个实施方式中,最大成像层数确定模块502,用于确定所述最大成像层数N,其中N=int(TA/TR);其中TR为所述重复时间,TA为所述采样窗时间,int(TA/TR)为(TA/TR)的取整值。
在一个实施方式中,成像数据获取模块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层成像数据。
基于上述描述,本发明实施方式还提出了磁共振成像系统的控制主机。
图6为根据本发明实施方式的磁共振成像系统的控制主机的结构图。
如图6所示,控制主机600包括处理器601、存储器602及存储在存储器602上并可在处理器601上运行的计算机程序,计算机程序被处理器601执行时实现如上任一项的PACE序列触发的磁共振成像方法。
其中,存储器602具体可以实施为电可擦可编程只读存储器(EEPROM)、快闪存储器(Flash memory)、可编程程序只读存储器(PROM)等多种存储介质。处理器601可以实施为包括一或多个中央处理器或一或多个现场可编程门阵列,其中现场可编程门阵列集成一或多个中央处理器核。具体地,中央处理器或中央处理器核可以实施为CPU或MCU或DSP等等。
综上所述,在本发明实施方式中,确定单激发快速自旋回波序列的重复时间和采样窗时间;基于采样窗时间和重复时间确定每个生理运动周期内的最大成像层数N,其中N为大于等于2的正整数;当PACE序列生成触发信号时,使能单激发快速自旋回波序列在至少一个采样窗时间内获取M层成像数据,其中M为大于等于2且小于等于N的正整数。由此可见,本发明实施方式在单个采样窗时间内获取多层成像数据,可以提高扫描速度,并减少成像时间。
而且,本发明实施方式可以基于PACE序列探测出的生理信号曲线确定采样窗时间,提高了采样窗时间的准确度。
另外,本发明实施方式考虑到总成像层数能否被M整除的不同状况,实现了在采样窗时间内采集多层成像数据的不同处理方式。
可以遵循一定规范的应用程序接口,将本发明实施方式所提出的前瞻性心电触发方法编写为安装到磁共振成像控制主机、个人电脑、移动终端等中的插件程序,也可以将其封装为应用程序以供用户自行下载使用。
可以通过指令或指令集存储的储存方式将本发明实施方式所提出的前瞻性心电触发方法存储在各种存储介质上。这些存储介质包括但是不局限于:软盘、光盘、DVD、硬盘、闪存等。另外,还可以将本发明实施方式所提出的前瞻性心电触发方法应用到基于闪存(Nand flash)的存储介质中,比如U盘、CF卡、SD卡、SDHC卡、MMC卡、SM卡、记忆棒、xD卡等。
以上所述,仅为本发明的较佳实施例而已,并非用于限定本发明的保护范围。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。

Claims (10)

1.一种前瞻性采集校正序列触发的磁共振成像方法(200),其特征在于,包括:
确定单激发快速自旋回波序列的重复时间和采样窗时间(201);
基于所述采样窗时间和所述重复时间确定每个生理运动周期内的最大成像层数N,其中N为大于等于2的正整数(202);
当前瞻性采集校正序列生成触发信号时,使能所述单激发快速自旋回波序列在至少一个采样窗时间内获取M层成像数据,其中M为大于等于2且小于等于N的正整数(203)。
2.根据权利要求1所述的前瞻性采集校正序列触发的磁共振成像方法(200),其特征在于,所述确定采样窗时间包括:
施加所述前瞻性采集校正序列以探测生理信号曲线;
基于所述生理信号曲线确定生理运动周期;
基于所述生理运动周期确定采样窗时间。
3.根据权利要求1所述的前瞻性采集校正序列触发的磁共振成像方法(200),其特征在于,所述基于采样窗时间和重复时间确定每个生理运动周期内的最大成像层数N包括:
确定所述最大成像层数N,其中N=int(TA/TR);
其中TR为所述重复时间,TA为所述采样窗时间,int(TA/TR)为(TA/TR)的取整值。
4.根据权利要求1所述的前瞻性采集校正序列触发的磁共振成像方法(200),其特征在于,
所述使能单激发快速自旋回波序列在至少一个采样窗时间内获取M层成像数据包括:
当预定的总成像层数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层成像数据。
5.一种前瞻性采集校正序列触发的磁共振成像装置(500),其特征在于,包括:
时间确定模块(501),用于确定单激发快速自旋回波序列的重复时间和采样窗时间;
最大成像层数确定模块(502),用于基于所述采样窗时间和所述重复时间确定每个生理运动周期内的最大成像层数N,其中N为大于等于2的正整数;
成像数据获取模块(503),用于当前瞻性采集校正序列生成触发信号时,使能所述单激发快速自旋回波序列在至少一个采样窗时间内获取M层成像数据,其中M为大于等于2且小于等于N的正整数。
6.根据权利要求5所述的前瞻性采集校正序列触发的磁共振成像装置(500),其特征在于,
时间确定模块(501),用于施加所述前瞻性采集校正序列以探测生理信号曲线;基于所述生理信号曲线确定生理运动周期;基于所述生理运动周期确定采样窗时间。
7.根据权利要求5所述的前瞻性采集校正序列触发的磁共振成像装置(500),其特征在于,
最大成像层数确定模块(502),用于确定所述最大成像层数N,其中N=int(TA/TR);其中TR为所述重复时间,TA为所述采样窗时间,int(TA/TR)为(TA/TR)的取整值。
8.根据权利要求5所述的前瞻性采集校正序列触发的磁共振成像装置(500),其特征在于,
成像数据获取模块(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层成像数据。
9.一种磁共振成像系统的控制主机(600),其特征在于,包括:存储器(601);处理器(602);其中所述存储器(601)中存储有可被所述处理器(602)执行的应用程序,用于使得所述处理器(602)执行如权利要求1至4中任一项所述的前瞻性采集校正序列触发的磁共振成像方法(200)。
10.一种计算机可读存储介质,其特征在于,所述计算机可读存储介质上存储计算机程序,所述计算机程序被处理器执行时实现如权利要求1至4中任一项所述的前瞻性采集校正序列触发的磁共振成像方法(200)。
CN202010836538.1A 2020-08-19 2020-08-19 一种前瞻性采集校正序列触发的磁共振成像方法和装置 Active CN114076913B (zh)

Priority Applications (2)

Application Number Priority Date Filing Date Title
CN202010836538.1A CN114076913B (zh) 2020-08-19 2020-08-19 一种前瞻性采集校正序列触发的磁共振成像方法和装置
US17/405,475 US20220057466A1 (en) 2020-08-19 2021-08-18 Accelerated haste for pace triggering

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010836538.1A CN114076913B (zh) 2020-08-19 2020-08-19 一种前瞻性采集校正序列触发的磁共振成像方法和装置

Publications (2)

Publication Number Publication Date
CN114076913A CN114076913A (zh) 2022-02-22
CN114076913B true CN114076913B (zh) 2024-05-17

Family

ID=80269501

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010836538.1A Active CN114076913B (zh) 2020-08-19 2020-08-19 一种前瞻性采集校正序列触发的磁共振成像方法和装置

Country Status (2)

Country Link
US (1) US20220057466A1 (zh)
CN (1) CN114076913B (zh)

Citations (20)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4717879A (en) * 1986-05-05 1988-01-05 Duke University Pulse sequence for NMR image acquisition
US5420509A (en) * 1992-02-26 1995-05-30 Kabushiki Kaisha Toshiba Magnetic resonance imaging apparatus
JPH0871060A (ja) * 1994-08-31 1996-03-19 Shimadzu Corp 核磁気共鳴イメージング装置
DE19616388A1 (de) * 1996-04-24 1997-11-06 Siemens Ag Pulssequenz zur Selektion mehrerer Schichten in einem Kernspintomographiegerät
JP2008136851A (ja) * 2006-11-10 2008-06-19 Toshiba Corp 磁気共鳴イメージング装置および磁気共鳴イメージング方法
CN101281242A (zh) * 2006-11-27 2008-10-08 西门子公司 在记录磁共振测量数据时确定运动的方法及相关设备
CN102078192A (zh) * 2009-11-27 2011-06-01 西门子公司 借助磁共振技术采集呼吸着的检查对象的测量数据的方法
CN102680930A (zh) * 2011-03-17 2012-09-19 西门子公司 用于设置匀场电流和射频中心频率的方法以及磁共振设备
CN104042216A (zh) * 2014-07-01 2014-09-17 中国科学院武汉物理与数学研究所 一种基于预扫描和非均匀采样的薄层快速磁共振成像方法
WO2014154544A1 (en) * 2013-03-25 2014-10-02 Fatnav Ekonomisk Förening Real-time motion correction for mri using fat navigators
CN105785297A (zh) * 2014-12-18 2016-07-20 西门子(深圳)磁共振有限公司 多片层数据采集方法及其磁共振成像方法
CN106597337A (zh) * 2016-12-09 2017-04-26 深圳先进技术研究院 一种磁共振t2*加权快速成像方法及装置
CN107212887A (zh) * 2016-03-22 2017-09-29 西门子(深圳)磁共振有限公司 一种心脏弥散成像的受激回波序列扫描方法和装置
CN107367703A (zh) * 2017-08-29 2017-11-21 上海联影医疗科技有限公司 磁共振扫描方法、系统、装置及存储介质
CN108113702A (zh) * 2017-12-21 2018-06-05 飞依诺科技(苏州)有限公司 一种超声波扫查方法和装置
CN108294753A (zh) * 2018-01-03 2018-07-20 上海东软医疗科技有限公司 磁共振定量信息图的获取方法和装置
CN108652624A (zh) * 2017-03-31 2018-10-16 西门子(深圳)磁共振有限公司 一种前瞻性心电触发方法和装置
CN109061533A (zh) * 2018-06-29 2018-12-21 深圳安科高技术股份有限公司 一种磁共振成像数据采集方法、装置及存储介质
GB201909212D0 (en) * 2019-06-26 2019-08-07 Cerebriu As An improved medical scan protocol for in-scanner patient acquisition analysis
CN111352055A (zh) * 2020-03-17 2020-06-30 无锡鸣石峻致医疗科技有限公司 一种前瞻性相位校正平面回波成像技术

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7209777B2 (en) * 2000-11-30 2007-04-24 General Electric Company Method and apparatus for automated tracking of non-linear vessel movement using MR imaging
DE10318428A1 (de) * 2003-04-23 2004-11-25 Siemens Ag Echtzeit-gesteuerte optimierte Magnet-Resonanz-Bildgebung unter Berücksichtigung von Geräte- und Patienten-spezifischen Grenzwerten
US8604787B2 (en) * 2006-04-27 2013-12-10 Stefan Posse Magnetic resonance spectroscopy with real-time correction of motion and frequency drift, and real-time shimming
DE102011005445B4 (de) * 2011-03-11 2014-10-09 Siemens Aktiengesellschaft Normalisierung von Magnetresonanzbilddaten bei bewegtem Tisch
DE102011007850B4 (de) * 2011-04-21 2015-05-13 Siemens Aktiengesellschaft Verfahren zur getriggerten Akquisition von Messdaten mittels eines Magnetresonanzgeräts unter Unterdrückung von unerwünschten Signalen, Magnetresonanzgerät, Computerprogramm und elektronisch lesbarer Datenträger
CN106572817B (zh) * 2014-08-11 2020-01-31 皇家飞利浦有限公司 利用针对4d磁共振成像的回顾性验证的前瞻性呼吸触发
DE102014218653A1 (de) * 2014-09-17 2016-03-17 Siemens Aktiengesellschaft Prospektive Bewegungskorrektur

Patent Citations (20)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4717879A (en) * 1986-05-05 1988-01-05 Duke University Pulse sequence for NMR image acquisition
US5420509A (en) * 1992-02-26 1995-05-30 Kabushiki Kaisha Toshiba Magnetic resonance imaging apparatus
JPH0871060A (ja) * 1994-08-31 1996-03-19 Shimadzu Corp 核磁気共鳴イメージング装置
DE19616388A1 (de) * 1996-04-24 1997-11-06 Siemens Ag Pulssequenz zur Selektion mehrerer Schichten in einem Kernspintomographiegerät
JP2008136851A (ja) * 2006-11-10 2008-06-19 Toshiba Corp 磁気共鳴イメージング装置および磁気共鳴イメージング方法
CN101281242A (zh) * 2006-11-27 2008-10-08 西门子公司 在记录磁共振测量数据时确定运动的方法及相关设备
CN102078192A (zh) * 2009-11-27 2011-06-01 西门子公司 借助磁共振技术采集呼吸着的检查对象的测量数据的方法
CN102680930A (zh) * 2011-03-17 2012-09-19 西门子公司 用于设置匀场电流和射频中心频率的方法以及磁共振设备
WO2014154544A1 (en) * 2013-03-25 2014-10-02 Fatnav Ekonomisk Förening Real-time motion correction for mri using fat navigators
CN104042216A (zh) * 2014-07-01 2014-09-17 中国科学院武汉物理与数学研究所 一种基于预扫描和非均匀采样的薄层快速磁共振成像方法
CN105785297A (zh) * 2014-12-18 2016-07-20 西门子(深圳)磁共振有限公司 多片层数据采集方法及其磁共振成像方法
CN107212887A (zh) * 2016-03-22 2017-09-29 西门子(深圳)磁共振有限公司 一种心脏弥散成像的受激回波序列扫描方法和装置
CN106597337A (zh) * 2016-12-09 2017-04-26 深圳先进技术研究院 一种磁共振t2*加权快速成像方法及装置
CN108652624A (zh) * 2017-03-31 2018-10-16 西门子(深圳)磁共振有限公司 一种前瞻性心电触发方法和装置
CN107367703A (zh) * 2017-08-29 2017-11-21 上海联影医疗科技有限公司 磁共振扫描方法、系统、装置及存储介质
CN108113702A (zh) * 2017-12-21 2018-06-05 飞依诺科技(苏州)有限公司 一种超声波扫查方法和装置
CN108294753A (zh) * 2018-01-03 2018-07-20 上海东软医疗科技有限公司 磁共振定量信息图的获取方法和装置
CN109061533A (zh) * 2018-06-29 2018-12-21 深圳安科高技术股份有限公司 一种磁共振成像数据采集方法、装置及存储介质
GB201909212D0 (en) * 2019-06-26 2019-08-07 Cerebriu As An improved medical scan protocol for in-scanner patient acquisition analysis
CN111352055A (zh) * 2020-03-17 2020-06-30 无锡鸣石峻致医疗科技有限公司 一种前瞻性相位校正平面回波成像技术

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
Prospective acquisition correction for head motion with image-based tracking for real-time fMRI;Stefan Thesen 等;《Magnetic Resonance in Medicine》;457-465 *
结合导航回波技术进行伪影矫正研究;董芳;《中国优秀硕士学位论文全文数据库 (医药卫生科技辑)》;E060-17 *

Also Published As

Publication number Publication date
CN114076913A (zh) 2022-02-22
US20220057466A1 (en) 2022-02-24

Similar Documents

Publication Publication Date Title
US7649354B2 (en) Method and apparatus for acquiring magnetic resonance imaging data
US5565777A (en) Method/apparatus for NMR imaging using an imaging scheme sensitive to inhomogeneity and a scheme insensitive to inhomogeneity in a single imaging step
US7689263B1 (en) Method and apparatus for acquiring free-breathing MR images using navigator echo with saturation RF pulse
US8831704B2 (en) Magnetic resonance method and apparatus for triggered acquisition of magnetic resonance measurement data
JP5854575B2 (ja) 磁気共鳴イメージング装置
US20190033410A1 (en) System and Method for Optimized Diffusion - Weighted Imaging
KR101809213B1 (ko) T2-이완 시간의 다중 슬라이스 영상화를 위한 방법 및 장치
JP2003528667A (ja) 適応的に選択されたフリップ角を用いる磁気共鳴撮像方法及びシステム
EP1227332B1 (en) Acquisition of high-temporal free-breathing MR images
US11033199B2 (en) Echo-planar imaging magnetic resonance elastography pulse sequence
US9791532B2 (en) Method and apparatus for magnetic resonance imaging
US10254373B2 (en) Method and apparatus for acquiring magnetic resonance data with navigator scans on opposite sides of the examination region
CN109242847B (zh) 一种磁共振弥散加权成像方法和装置
EP3168636A2 (en) Gradient echo mri using velocity-selective excitation
CN108261199B (zh) 用于利用分组数据采集来执行对象的时间分辨mr成像的mri系统和方法
CN112394311B (zh) 具有改进的导航器的mri系统
US20180100906A1 (en) Magnetic resonance imaging apparatus and method of operating the same
JP4230875B2 (ja) 磁気共鳴イメージング装置
CN114076913B (zh) 一种前瞻性采集校正序列触发的磁共振成像方法和装置
US10667760B2 (en) Stimulated echo sequence scanning magnetic resonance method and apparatus for heart diffusion imaging
US20150094562A1 (en) Magnetic resonance imaging with dynamic inversion preparation
US11300641B2 (en) Preventing interruptions in magnetic resonance measurements
CN112285622A (zh) 核磁共振成像控制方法、装置及计算机可读存储介质
CN114814685B (zh) 一种射频场强度测量方法、装置、电子设备及存储介质
US20230408616A1 (en) Single-voxel spectroscopy for quantitation of myocardial metabolites

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant