CN109717869B - 磁共振成像过程中的运动监测方法、计算机程序、存储设备 - Google Patents
磁共振成像过程中的运动监测方法、计算机程序、存储设备 Download PDFInfo
- Publication number
- CN109717869B CN109717869B CN201711047981.5A CN201711047981A CN109717869B CN 109717869 B CN109717869 B CN 109717869B CN 201711047981 A CN201711047981 A CN 201711047981A CN 109717869 B CN109717869 B CN 109717869B
- Authority
- CN
- China
- Prior art keywords
- noise
- imaging
- coil
- repetition time
- navigation
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 36
- 238000012544 monitoring process Methods 0.000 title claims abstract description 19
- 238000002595 magnetic resonance imaging Methods 0.000 title claims abstract description 15
- 238000004590 computer program Methods 0.000 title claims description 10
- 238000003860 storage Methods 0.000 title claims description 5
- 230000008569 process Effects 0.000 title abstract description 7
- 238000003384 imaging method Methods 0.000 claims abstract description 49
- 230000005284 excitation Effects 0.000 claims description 8
- 230000001960 triggered effect Effects 0.000 claims description 3
- 238000001914 filtration Methods 0.000 abstract description 6
- 230000000241 respiratory effect Effects 0.000 description 25
- 230000008859 change Effects 0.000 description 9
- 238000001208 nuclear magnetic resonance pulse sequence Methods 0.000 description 9
- 238000010586 diagram Methods 0.000 description 8
- 238000002474 experimental method Methods 0.000 description 6
- 238000004364 calculation method Methods 0.000 description 4
- 230000029058 respiratory gaseous exchange Effects 0.000 description 4
- 238000013459 approach Methods 0.000 description 3
- 230000036387 respiratory rate Effects 0.000 description 3
- 238000012360 testing method Methods 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000005684 electric field Effects 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 238000007781 pre-processing Methods 0.000 description 2
- 230000002411 adverse Effects 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012806 monitoring device Methods 0.000 description 1
- 230000000737 periodic effect Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
- G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
- G01R33/48—NMR imaging systems
- G01R33/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
- G01R33/56—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
- G01R33/565—Correction of image distortions, e.g. due to magnetic field inhomogeneities
- G01R33/56509—Correction of image distortions, e.g. due to magnetic field inhomogeneities due to motion, displacement or flow, e.g. gradient moment nulling
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
- G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
- G01R33/48—NMR imaging systems
- G01R33/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
- G01R33/56—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
- G01R33/5608—Data processing and visualization specially adapted for MR, e.g. for feature analysis and pattern recognition on the basis of measured MR data, segmentation of measured MR data, edge contour detection on the basis of measured MR data, for enhancing measured MR data in terms of signal-to-noise ratio by means of noise filtering or apodization, for enhancing measured MR data in terms of resolution by means for deblurring, windowing, zero filling, or generation of gray-scaled images, colour-coded images or images displaying vectors instead of pixels
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
- G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
- G01R33/48—NMR imaging systems
- G01R33/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
- G01R33/56—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
- G01R33/567—Image 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/5673—Gating 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)
- General Physics & Mathematics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Condensed Matter Physics & Semiconductors (AREA)
- High Energy & Nuclear Physics (AREA)
- Signal Processing (AREA)
- General Health & Medical Sciences (AREA)
- Radiology & Medical Imaging (AREA)
- Biophysics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Physiology (AREA)
- Pulmonology (AREA)
- Power Engineering (AREA)
- Cardiology (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
本发明的实施例提供了一种磁共振成像过程中的运动监测方法,包括:在成像扫描序列的每个成像重复时间之前或之后采集接收线圈的噪声;确定接收线圈中与扫描对象的运动相关联的主要线圈通道;确定各主要线圈通道的噪声的幅度的平方和;以及,对主要线圈通道的噪声的幅度的平方和进行滤波以获取扫描对象的运动轨迹。
Description
技术领域
本发明涉及医疗检测领域,尤其涉及一种磁共振成像过程中的运动监测方法。
背景技术
在利用磁共振(Magnetic Resonance,MR)扫描技术对人体进行成像时,常常会出现由于患者的呼吸运动产生的图像伪影,为了去除或降低呼吸运动伪影,通常需要对呼吸运动进行实时监测。
现有技术中,一种方法是通过外部的监测设备来监测呼吸运动,并根据实时监测结果来确定何时进行图像数据的采集,这种方法具有以下几个缺点:1)增加了设备设置的时间,2)需要花费较多的时间对患者进行摆位,3)增加了额外的设备成本。还有一种导航回波技术,这种技术需要在图像采集序列之前施加额外的射频脉冲和梯度脉冲来进行呼吸门控或呼吸补偿,这种方法会延长扫描时间,降低扫描效率,并且容易对图像采集过程造成不利的影响。
另外,本领域技术人员也曾提出,通过测量射频线圈的阻抗来监视呼吸运动,但是这种方法需要为磁共振成像***提供额外的拾取线圈(pick-up coil)来与用于图像采集的射频线圈产生耦合,需要提供专门的射频监视脉冲来代替成像序列中的射频激发脉冲,并且,还需要对射频放大器进行漂移补偿。还有一种类似的方法是使用额外的定向耦合器来代替该拾取线圈,在实施这些方法时,需要充分考虑成本、可行性方面的影响。
发明内容
本发明的一个目的在于提供一种准确性较高、成本较低的呼吸运动监测方法。
本发明的实施例提供了一种磁共振成像过程中的运动监测方法,包括:在成像扫描序列的每个成像重复时间之前或之后采集接收线圈的噪声;确定接收线圈中与扫描对象的运动相关联的主要线圈通道;确定各主要线圈通道的噪声的幅度的平方和;以及,对主要线圈通道的噪声的幅度的平方和进行滤波以获取扫描对象的运动轨迹。
本发明的示例性实施例还提供了一种计算机程序,当该计算机程序运行于一磁共振成像***中时,使该磁共振成像***执行上述的运动监测方法。
本发明的示例性实施例还提供了一种存储设备,用于存储上述计算机程序。
通过下面的详细描述、附图以及权利要求,其他特征和方面会变得清楚。
附图说明
通过结合附图对于本发明的示例性实施例进行描述,可以更好地理解本发明,在附图中:
图1为本发明一个实施例的运动监测方法的流程图;
图2示例性地示出了本发明实施例中在采集接收线圈的噪声时执行的导航序列;
图3为在每个成像重复时间之前施加导航序列的示例图;
图4为在每个成像重复时间之后施加导航序列的示例图;
图5为本发明一个示例中提供的呼吸率计算脉冲序列的示意图;
图6-图9分别示出了利用本发明的实施例获取呼吸运动轨迹的各流程阶段的数据;
图10-图13分别示出了利用本发明的实施例监测志愿者的呼吸运动轨迹的各流程阶段的数据;
图14为利用外部监测仪得到的志愿者的呼吸运动轨迹。。
具体实施方式
以下将描述本发明的具体实施方式,需要指出的是,在这些实施方式的具体描述过程中,为了进行简明扼要的描述,本说明书不可能对实际的实施方式的所有特征均作详尽的描述。应当可以理解的是,在任意一种实施方式的实际实施过程中,正如在任意一个工程项目或者设计项目的过程中,为了实现开发者的具体目标,为了满足***相关的或者商业相关的限制,常常会做出各种各样的具体决策,而这也会从一种实施方式到另一种实施方式之间发生改变。此外,还可以理解的是,虽然这种开发过程中所作出的努力可能是复杂并且冗长的,然而对于与本发明公开的内容相关的本领域的普通技术人员而言,在本公开揭露的技术内容的基础上进行的一些设计,制造或者生产等变更只是常规的技术手段,不应当理解为本公开的内容不充分。
除非另作定义,权利要求书和说明书中使用的技术术语或者科学术语应当为本发明所属技术领域内具有一般技能的人士所理解的通常意义。本发明专利申请说明书以及权利要求书中使用的“第一”、“第二”以及类似的词语并不表示任何顺序、数量或者重要性,而只是用来区分不同的组成部分。“一个”或者“一”等类似词语并不表示数量限制,而是表示存在至少一个。“包括”或者“包含”等类似的词语意指出现在“包括”或者“包含”前面的元件或者物件涵盖出现在“包括”或者“包含”后面列举的元件或者物件及其等同元件,并不排除其他元件或者物件。“连接”或者“相连”等类似的词语并非限定于物理的或者机械的连接,也不限于是直接的还是间接的连接。
图1为本发明一个实施例的运动监测方法的流程图,图2、图3分别示例性地示出了利用该运动监测方法进行呼吸门控的示意图。如图1至图3所示,该方法包括步骤S110、S120、S130和S140。
在步骤S110中,在成像扫描序列的每个成像重复时间之前或之后采集接收线圈的噪声。本领域技术人员应当理解,上述成像扫描序列通常包括用于激发射频线圈的激发脉冲,用于控制梯度场的梯度脉冲等,其成像重复时间(TR)为一个成像脉冲序列开始到结束的时间。
研究表明,接收线圈的噪声变化与其有效电阻有关,具体如下式(1)所示:
其中,Vnoise代表接收线圈的噪声变化,k代表玻尔兹曼(Boltzman)常数,T代表开氏(Kelvin)温度,BW代表检测***的带宽,Reff代表接收线圈的有效电阻。
而接收线圈的有效电阻与靠近其的人体的电阻有关,具体如下式(2)所示:
Reff=Rbody+Rcoil+Relectronics (2)
其中,Rbody代表人体的电阻,Rcoil代表接收线圈自身的电阻,Relectronics代表其它电子器件的电阻,实际中,Rcoil和Relectronics的值相对Rbody小很多,因此接收线圈的有效电阻主要来自人体的电阻。
人体的电阻可以通过下式(3)来表示:
其中,σ代表人体组织的导电率,E代表由接收线圈中单位幅度的电流引起的人体中的电场,*表示共轭转置,I表示接收线圈的电流,V表示接收线圈的电压,r代表空间位置,t代表时间。
研究表明,人体的运动能够改变其自身的导电率、电场等,这也导致人体电阻的改变,通过上式(2)可知,接收线圈的电阻的大部分是人体贡献的。当将人体置于磁共振成像***的扫描腔中时,能够引起靠近人体的接收线圈的阻抗的改变。
因此,本发明通过在步骤S110中采集接收线圈的噪声,并在后续的步骤中对采集的噪声进行适当的分析,以获取与之相关的扫描对象的运动轨迹,该运动轨迹可以为,例如,人体的呼吸运动。
在步骤S120中,确定接收线圈中与扫描对象的运动相关联的主要线圈通道。本领域技术人员应当理解,接收线圈一般具有由多个线圈通道,该多个线圈通道可以以阵列形式分布于不同的位置以分别靠近人体的不同部位,一般来说,越靠近人体的运动部位的线圈通道,其阻抗变化越明显,因此,本步骤通过确定与扫描对象的运动相关联的线圈通道作为主要线圈通道,可以提升后续分析的准确性和简便性。
步骤S130中,求该主要线圈通道的噪声的幅度的平方和。
步骤S140中,对该主要线圈通道的噪声的幅度的平方和进行滤波以获取扫描对象的运动轨迹。
图2示例性地示出了本发明实施例中在采集接收线圈的噪声时执行的导航序列。如图2所示,在步骤S110中,“采集接收线圈的噪声”是在停止施加射频激发脉冲期间进行的。例如,可以在成像扫描序列的每个成像重复时间之前或之后执行导航序列,其中,在执行该导航序列期间,停止施加射频激发脉冲RF,并采集接收线圈的噪声。
图2中示例的导航序列为二维笛卡尔梯度回波脉冲序列,在其它实施例中,也可以采用其它形式的脉冲序列,只要在该脉冲序列中不施加射频激发脉冲即可,例如,也可以停止施加梯度脉冲Gx、Gy、Gz。通过设置这样的导航序列,使得在切换导航序列和成像扫描序列时,仅需要切换射频激发脉冲的触发和关闭,无需对脉冲序列进行进行其它的改变。可以通过采用较大的采样带宽(例如±500kHz)来使噪声采集时间达到最小,因此,能够达到足够的噪声实时采集点,例如一个导航重复时间内具有500个噪声采集点,通过在这些采集点采集的噪声进一步地可以用于分析呼吸运动轨迹。
图3为在每个成像重复时间之前施加导航序列的示例图,其中曲线为接收线圈的噪声变化曲线,两条横线min和max分别表示噪声幅度的最小门限值和最大门限值,该第一预设范围可以大于或等于该最小门限值并小于或等于该最大门限值。如图3所示,在每个成像重复时间之前施加的导航序列可以被设置为:1)与相邻的成像重复时间相间隔;2)包括多个导航重复时间,且相邻导航重复时间互相间隔。
在图3所示的示例中,当在执行该导航序列期间采集的噪声处于第一预设范围内时,触发该成像扫描序列的一个成像重复时间,以采集图像数据。该第一预设范围可以对应于能够较好地采集到感兴趣区域的图像数据的位置,例如,人体在呼吸过程中,其横膈膜靠近扫描切片的中心位置并且图像数据的读出方向垂直于该横膈膜。可以通过这种方式进行呼吸门控,使得成像扫描序列可以在恰当的时机被触发,以减少由于人体呼吸运动带来的图像伪影。并且,本发明实施例可以在图像采集期间持续地监视接收线圈的噪声变化,相较传统的触发方式,在成像重复时间之后无需额外的时间来等待由于射频激发脉冲引起的纵向弛豫的恢复。
图4为在每个成像重复时间之后施加导航序列的示例图,其中曲线为接收线圈的噪声变化曲线,两条横线min和max分别表示噪声幅度的最小门限值和最大门限值,该第二预设范围可以大于或等于该最小门限值并小于或等于该最大门限值。如图4所示,在每个成像重复时间之后施加的导航序列可以包括一个导航重复时间,并且该导航重复时间紧随其前面的成像重复时间,即与其前面的成像重复时间之间没有时间间隔。
在图4所示的示例中,当在任一导航重复时间内采集的接收线圈的噪声的幅度在第二预设范围外时,将在该导航重复时间之前的成像重复时间内采集的图像数据确定为无效图像数据;当在任一导航重复时间内采集的接收线圈的噪声的幅度在第二预设范围内时,将在该导航重复时间之前的成像重复时间内采集的图像数据确定为有效图像数据。
该第二预设范围可以与第一预设范围相同。通过这种方式进行呼吸门控,使得能够通过噪声导航来确定在相应的成像扫描序列中采集的图像数据是否被接受,在后续的图像处理过程中,可以选择丢弃该无效图像数据,而仅基于该有效图像数据进行图像重建。以减少由于呼吸运动带来的图像伪影。
图5为本发明一个示例中提供的呼吸率计算脉冲序列的示意图,如图5所示,可以在步骤S110之前通过预先执行该呼吸率计算脉冲序列来获取合适的第一预设范围和第二预设范围,该呼吸率计算脉冲序列包括依次设置的多个导航重复时间,通过在每个导航重复时间中实时记录感兴趣区域(例如横膈膜)的位置随接收线圈的噪声的变化信息,通过对记录的该位置变化信息进行分析,可以找出感兴趣区域处于适于进行成像扫描的位置时的噪声值,并根据找出的噪声值确定第一预设范围和第二预设范围,例如,该第一预设范围可以为找出的噪声值的最大值和最小值之间的范围,也可以进一步根据需要对该范围进行合适的调整。上述感兴趣区域处于适于进行成像扫描的位置可以为,例如,横膈膜靠近切片中心且图像读出方向垂直于横膈膜的位置。
需要说明的是,在图3-图5的示意图中,其中标有“imaging”的任意线框用于示意成像扫描序列的一个成像重复时间,其中标有“nNAV”的任意线框用于示意导航序列的一个导航重复时间。
在步骤S110中,可以通过以上任一实施方式采集接收线圈的每一个线圈通道的噪声,并在步骤S120中,可以根据采集的噪声确定出其中与呼吸运动相关联的主要线圈通道。
可选地,在步骤S110中,“采集接收线圈的噪声”还可以为包括对采集的原始的噪声数据进行预处理的步骤,则在步骤S120中可以基于预处理后的该噪声确定上述主要线圈。对采集的原始的噪声数据进行预处理的过程可以包括:
步骤一:利用卡尔曼滤波器对原始的噪声数据进行滤波;和/或
步骤二:对原始的噪声数据进行去耦合,例如,去除接收线圈阵列的非对角线上的噪声耦合。
步骤S120可以包括以下步骤:
步骤一:分析从各线圈通道采集的噪声的频域信息。
步骤二:将频率波动较大的频域信息所对应的线圈通道确定该主要线圈通道。
图6-图9分别示出了利用本发明的实施例获取呼吸运动轨迹的各流程阶段的数据,其中图6为分别从多个线圈通道channel 1~channel 12采集的多组噪声幅度信息,图7为该多组噪声幅度信息的频域信息,图8为将主要线圈通道的噪声幅度信息合并后的噪声幅度信息,图9为通过对图8的噪声幅度信息进行滤波得到的呼吸运动轨迹。例如图6、图7所示,可以通过对该多组噪声幅度信息进行傅立叶变换来获取对应的多组频域信息,并通过分析该多组频域信息可知,其中左下角两组频域信息中,具有较大的频率波动。因此可以将对应的该两组频域信息分别进行反傅立叶变换来重新获取对应的噪声幅度信息,如图6中被圆圈标注的两组噪声幅度信息,通过将该两组噪声幅度信息与从各线圈通道采集的噪声幅度信息进行匹配,可以确定与呼吸运动相关联的主要线圈(channel 10,channel 11)。
因此,在步骤S130之前,还可以包括以下步骤:将从各线圈通道采集的噪声的频域信息变换为对应的噪声幅度信息。例如,通过上述的反傅立叶变换,得到如图6中被圆圈标注的两组噪声幅度信息。
在步骤S130中,通过求该两组噪声幅度信息中的噪声幅度的平方和,以得到如图8所示的噪声幅度信息。
在步骤S140中,对该噪声幅度的平方和进行滤波,可以获取图9所示的较为平滑的呼吸运动轨迹,滤波的方式可以为低通滤波。
本发明的实施例还可以提供一种计算机程序,当该计算机程序运行于一磁共振成像***中时,使该磁共振成像***执行上述任一实施例的运动监测方法。例如,该计算机程序包括多条指令以用于执行上述运动监测方法的各步骤。上述计算机程序可以安装于,例如磁共振成像***的计算机控制台中,其可以通过计算机控制台的处理器加载并执行。
本发明的实施例还可以提供一种存储设备,其用于存储该计算机程序,
以下将描述为了解决本发明的技术问题进行的两组试验:
该两组试验基于磁共振扫描成像产品Kizuna 1.5T进行,其中具有12个表面线圈通道组成的表面线圈阵列,将体模放置在扫描台上,该表面线圈阵列相对于主磁场固定。
在第一组试验中,利用软件“table rocker”周期性地前后移动该扫描台,并将该体模作为扫描对象来实施本发明的实施例。
在第二组试验中,关闭软件“table rocker”,扫描台不再移动。
在第一组试验的条件下将该体模作为扫描对象来实施本发明的实施例时,各阶段的试验结果如图6-图9所示,不再赘述。在第二组试验的条件下将该体模作为扫描对象来实施本发明的实施例时,未采集到噪声的周期性信号。两组试验的对比显示,接收线圈的噪声动态变化主要是由体模的运动引起的。并且,即使将扫描台的位移、速度等调整到最大值,该运动引起的噪声变化曲线都基本呈正弦曲线。另外,接收线圈阵列中仅部分通道对有效噪声变化具有贡献。
图10-图13分别示出了利用本发明的实施例监测志愿者的呼吸运动轨迹的各流程阶段的数据,图14为利用外部监测仪得到的志愿者的呼吸运动轨迹。其中利用本发明的实施例监测得到的呼吸运动轨迹如图13所示,其与图14所示的利用外部监测仪监测的运动轨迹相匹配。
上述试验和临床结果表明,通过本发明的实施例,能够有效地监测扫描对象的运动轨迹,例如呼吸运动轨迹。
本发明的实施例通过在成像扫描序列的成像重复时间之前或之后执行套行序列的重复时间,来采集接收线圈的噪声,由于该噪声的变化主要由扫描对象的运动引起,因此,通过对采集的噪声进行进一步处理,无需额外的硬件,便可以得到扫描对象的运动曲线,成本较低、准确性较高,能够帮助减少对扫描对象的成像的运动伪影。
另外,由于可以在停止施加脉冲扫描的情况下执行导航序列,避免了额外的等待时间,因此避免了扫描时间的延长。
上面已经描述了一些示例性实施例,然而,应该理解的是,可以做出各种修改。例如,如果所描述的技术以不同的顺序执行和/或如果所描述的***、架构、设备或电路中的组件以不同方式被组合和/或被另外的组件或其等同物替代或补充,则可以实现合适的结果。相应地,其他实施方式也落入权利要求的保护范围内。
Claims (10)
1.一种磁共振成像过程中的运动监测方法,包括:
在成像扫描序列的每个成像重复时间之前或之后采集接收线圈的噪声;
确定接收线圈中与扫描对象的运动相关联的主要线圈通道;
确定各主要线圈通道的噪声的幅度的平方和;
对所述主要线圈通道的噪声的幅度的平方和进行滤波以获取扫描对象的运动轨迹。
2.如权利要求1所述的方法,其特征在于,所述“采集接收线圈的噪声”包括:
在成像扫描序列的每个成像重复时间之前或之后执行导航序列,所述导航序列执行期间,停止施加射频激发脉冲并采集所述接收线圈的噪声。
3.根据权利要求2所述的方法,其特征在于,在每个成像重复时间之前施加的导航序列被设置为:
与相邻的成像重复时间相间隔;
包括多个导航重复时间,且相邻导航重复时间互相间隔。
4.根据权利要求3所述的方法,其特征在于,当采集的接收线圈的噪声在第一预设范围内时,触发所述成像扫描序列的一个成像重复时间,以采集图像数据。
5.根据权利要求2所述的方法,其特征在于,在每个成像重复时间之后施加的导航序列包括一个导航重复时间,所述一个导航重复时间紧随其前面的成像重复时间。
6.根据权利要求5所述的方法,其特征在于,当在任一导航重复时间内采集的接收线圈的噪声的幅度在第二预设范围外时,将所述导航重复时间之前的成像重复时间内采集的图像数据确定为无效图像数据;当在任一导航重复时间内采集的接收线圈的噪声的幅度在第二预设范围内时,将在所述导航重复时间之前的成像重复时间内采集的图像数据确定为有效图像数据。
7.如权利要求1所述的方法,其特征在于,所述“确定接收线圈中与运动相关联的主要线圈通道”包括:
分析从各线圈通道采集的噪声的频域信息;
将频率波动较大的频域信息所对应的线圈通道确定为所述主要线圈通道。
8.如权利要求7所述的方法,其特征在于,所述“确定各主要线圈通道的噪声的幅度的平方和”之前还包括:
将从各线圈通道采集的噪声的频域信息变换为对应的噪声幅度信息。
9.一种计算机程序,当所述计算机程序运行于一磁共振成像***中时,使所述磁共振成像***执行权利要求1-8任一项所述的运动监测方法。
10.一种存储设备,用于存储权利要求9所述的计算机程序。
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711047981.5A CN109717869B (zh) | 2017-10-31 | 2017-10-31 | 磁共振成像过程中的运动监测方法、计算机程序、存储设备 |
US16/174,366 US10859659B2 (en) | 2017-10-31 | 2018-10-30 | Motion monitoring method during MR imaging, computer process, and storage device |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711047981.5A CN109717869B (zh) | 2017-10-31 | 2017-10-31 | 磁共振成像过程中的运动监测方法、计算机程序、存储设备 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109717869A CN109717869A (zh) | 2019-05-07 |
CN109717869B true CN109717869B (zh) | 2024-05-14 |
Family
ID=66242841
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711047981.5A Active CN109717869B (zh) | 2017-10-31 | 2017-10-31 | 磁共振成像过程中的运动监测方法、计算机程序、存储设备 |
Country Status (2)
Country | Link |
---|---|
US (1) | US10859659B2 (zh) |
CN (1) | CN109717869B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113100741B (zh) * | 2020-01-13 | 2024-02-27 | 上海联影医疗科技股份有限公司 | 一种磁共振扫描方法、设备及存储介质 |
DE102020111776A1 (de) * | 2020-04-30 | 2021-11-04 | Schölly Fiberoptic GmbH | Verfahren zur Berechnung eines Ausgabebildes und/oder einer Folge von Ausgabebildern aus Bildrohdaten, Bildaufnahmevorrichtung, Bildaufnahmesystem und Gesamtsystem |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2006314491A (ja) * | 2005-05-12 | 2006-11-24 | Hitachi Medical Corp | 磁気共鳴撮影装置 |
CN1951323A (zh) * | 2005-10-19 | 2007-04-25 | Ge医疗***环球技术有限公司 | 磁共振成像设备 |
CN102078196A (zh) * | 2009-11-27 | 2011-06-01 | 西门子公司 | 借助磁共振采集相应于对象的呼吸运动的信号数据的方法 |
WO2012137563A1 (ja) * | 2011-04-01 | 2012-10-11 | 株式会社日立メディコ | 磁気共鳴イメージング装置および磁気共鳴イメージング方法 |
CN103271740A (zh) * | 2013-06-14 | 2013-09-04 | 刘再毅 | 核磁共振成像方法和*** |
WO2013165571A1 (en) * | 2012-04-30 | 2013-11-07 | The General Hospital Corporation | System and method for quiet magnetic resonance imaging |
CN105078459A (zh) * | 2014-05-20 | 2015-11-25 | 西门子公司 | 用于测量患者的呼吸过程的方法、测量设备和磁共振装置 |
CN106233154A (zh) * | 2014-04-22 | 2016-12-14 | 皇家飞利浦有限公司 | 使用预脉冲和导航器的具有运动校正的磁共振成像 |
CN106291417A (zh) * | 2016-09-13 | 2017-01-04 | 深圳市贝斯达医疗股份有限公司 | 磁共振成像方法和装置 |
Family Cites Families (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2005525183A (ja) * | 2002-05-13 | 2005-08-25 | コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ | 磁気共鳴イメージング |
US6788063B1 (en) * | 2003-02-26 | 2004-09-07 | Ge Medical Systems Technology Company, Llc | Method and system for improving transient noise detection |
DE102004019590B4 (de) * | 2004-04-22 | 2014-01-02 | Siemens Aktiengesellschaft | Verfahren zur kernspintomographischen Messung dynamischer Signale |
EP2069815A1 (en) * | 2006-09-13 | 2009-06-17 | Eidgenössische Technische Hochschule (ETH) | Method and apparatus for mri coil array compression |
US20080315879A1 (en) * | 2007-06-19 | 2008-12-25 | General Electric Company | System and apparatus for electromagnetic noise detection in an mr imaging scanner environment |
US8427153B2 (en) | 2010-01-15 | 2013-04-23 | Beth Israel Deaconess Medical Center, Inc. | Method for motion correction in magnetic resonance imaging using radio frequency coil arrays |
AU2015311825B2 (en) * | 2014-09-05 | 2019-01-17 | Hyperfine Operations, Inc. | Low field magnetic resonance imaging methods and apparatus |
-
2017
- 2017-10-31 CN CN201711047981.5A patent/CN109717869B/zh active Active
-
2018
- 2018-10-30 US US16/174,366 patent/US10859659B2/en active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2006314491A (ja) * | 2005-05-12 | 2006-11-24 | Hitachi Medical Corp | 磁気共鳴撮影装置 |
CN1951323A (zh) * | 2005-10-19 | 2007-04-25 | Ge医疗***环球技术有限公司 | 磁共振成像设备 |
CN102078196A (zh) * | 2009-11-27 | 2011-06-01 | 西门子公司 | 借助磁共振采集相应于对象的呼吸运动的信号数据的方法 |
WO2012137563A1 (ja) * | 2011-04-01 | 2012-10-11 | 株式会社日立メディコ | 磁気共鳴イメージング装置および磁気共鳴イメージング方法 |
WO2013165571A1 (en) * | 2012-04-30 | 2013-11-07 | The General Hospital Corporation | System and method for quiet magnetic resonance imaging |
CN103271740A (zh) * | 2013-06-14 | 2013-09-04 | 刘再毅 | 核磁共振成像方法和*** |
CN106233154A (zh) * | 2014-04-22 | 2016-12-14 | 皇家飞利浦有限公司 | 使用预脉冲和导航器的具有运动校正的磁共振成像 |
CN105078459A (zh) * | 2014-05-20 | 2015-11-25 | 西门子公司 | 用于测量患者的呼吸过程的方法、测量设备和磁共振装置 |
CN106291417A (zh) * | 2016-09-13 | 2017-01-04 | 深圳市贝斯达医疗股份有限公司 | 磁共振成像方法和装置 |
Also Published As
Publication number | Publication date |
---|---|
CN109717869A (zh) | 2019-05-07 |
US20190128990A1 (en) | 2019-05-02 |
US10859659B2 (en) | 2020-12-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US8744154B2 (en) | System and method for acquiring magnetic resonance imaging (MRI) data | |
US6546275B2 (en) | Determination of the arterial input function in dynamic contrast-enhanced MRI | |
US9897677B2 (en) | Method for correcting errors associated with asynchronous timing offsets between transmit and receive clocks in MRI wireless radiofrequency coils | |
DE60318438T2 (de) | Blutströmungs-gate-mri | |
US9804246B2 (en) | Magnetic resonance imaging apparatus | |
US20140343399A1 (en) | Seed-Based Connectivity Analysis in Functional MRI | |
EP2730935A1 (en) | Motional error correction in functional magnetic resonance imaging | |
EP2010933A1 (en) | Mri of a continuously moving object involving motion compensation | |
US20060244445A1 (en) | Motion compensation for magnetic resonance imaging | |
JP6559431B2 (ja) | Ecgシグナル中に現れるmriシークエンシングノイズの動的解除 | |
JPH05154130A (ja) | 体動アーチファクト除去方法 | |
JP6534532B2 (ja) | Ecgシグナル中に現れるmriシークエンシングノイズのアナログ解除 | |
EP1227332B1 (en) | Acquisition of high-temporal free-breathing MR images | |
US10739420B2 (en) | Motion-robust transmit radio frequency field mapping in magnetic resonance imaging using interleaved bloch-siegert shifting | |
US20150346307A1 (en) | Method and magnetic resonance apparatus for acquiring magnetic resonance data with a prospective motion correction | |
US10330764B2 (en) | Magnetic resonance measuring apparatus and image processing method | |
CN109717869B (zh) | 磁共振成像过程中的运动监测方法、计算机程序、存储设备 | |
US10330760B2 (en) | System and method for assessing T2-relaxation times with improved accuracy | |
US11163029B2 (en) | MRI system with improved navigator | |
JP4349647B2 (ja) | 磁気共鳴イメージング装置 | |
JP3699963B2 (ja) | 磁気共鳴映像装置 | |
US11079448B2 (en) | System and method for proton density mapping and receiver bias correction using magnetic resonance fingerprinting (MRF) | |
WO2002086530A1 (en) | Determination of the arterial input function in dynamic contrast-enhanced mri | |
US11497412B2 (en) | Combined oxygen utilization, strain, and anatomic imaging with magnetic resonance imaging | |
US10139466B2 (en) | System and method for simultaneous magnetic resonance elastography and water-fat imaging |
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 |