CN101173978B - 用于使磁共振成像数据同步于身体运动的设备 - Google Patents

用于使磁共振成像数据同步于身体运动的设备 Download PDF

Info

Publication number
CN101173978B
CN101173978B CN2007101938423A CN200710193842A CN101173978B CN 101173978 B CN101173978 B CN 101173978B CN 2007101938423 A CN2007101938423 A CN 2007101938423A CN 200710193842 A CN200710193842 A CN 200710193842A CN 101173978 B CN101173978 B CN 101173978B
Authority
CN
China
Prior art keywords
sequence information
time sequence
equipment
imaging data
data
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.)
Expired - Fee Related
Application number
CN2007101938423A
Other languages
English (en)
Other versions
CN101173978A (zh
Inventor
A·C·拉尔森
R·D·怀特
O·P·西莫内蒂
G·劳比
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.)
Cleveland Clinic Foundation
Siemens Medical Solutions USA Inc
Original Assignee
Cleveland Clinic Foundation
Siemens Medical Solutions USA Inc
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 Cleveland Clinic Foundation, Siemens Medical Solutions USA Inc filed Critical Cleveland Clinic Foundation
Publication of CN101173978A publication Critical patent/CN101173978A/zh
Application granted granted Critical
Publication of CN101173978B publication Critical patent/CN101173978B/zh
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Images

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
    • 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
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/72Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7271Specific aspects of physiological measurement analysis
    • A61B5/7285Specific aspects of physiological measurement analysis for synchronising or triggering a physiological measurement or image acquisition with a physiological event or waveform, e.g. an ECG signal
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/72Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7271Specific aspects of physiological measurement analysis
    • A61B5/7285Specific aspects of physiological measurement analysis for synchronising or triggering a physiological measurement or image acquisition with a physiological event or waveform, e.g. an ECG signal
    • A61B5/7289Retrospective gating, i.e. associating measured signals or images with a physiological event after the actual measurement or image acquisition, e.g. by simultaneously recording an additional physiological signal during the measurement or image acquisition

Landscapes

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

Abstract

用于使磁共振成像数据同步于身体运动的设备。通过从MR成像数据本身抽取时序信号,而不是仅仅依赖仅为了时序所获取的附加数据,使MR图像同步于患者的运动,例如同步于心脏的跳动、肺的呼吸或肢的运动。

Description

用于使磁共振成像数据同步于身体运动的设备
技术领域
本发明涉及磁共振(MR)成像,且更具体地涉及MR成像数据与患者运动的同步。
背景技术
使MR图像同步于患者的运动,例如心脏的跳动、肺的呼吸或肢的运动提供给诊断医师具有公知的与所感兴趣的运动,例如与心动周期的阶段对应的图像。这种同步对于电影(cine)和静止MR图像均可以是有用的。为了取得同步,有必要拥有指示身体的阶段或位置的时序信号。例如在心脏成像时时序信号可能指示每个心动周期的开始。
用于提供心脏成像中的时序信号的一个技术是:将适当的电极连接到患者用于监视患者的ECG,同时MR成像数据被收集。然而,用于MR中的脉冲磁场梯度和磁场可以干扰ECG信号的收集。已经研制出特殊算法来努力克服这些困难。Chia等人,“Performance of QRS Detection for CardiacMagnetic Resonance Imaging with a Novel VectorcardiographicTriggering Method,”,Journal of Magnetic Resonance Imaging,12:678-688(2000)。Fischer等人,“Novel Real-Time k-Wave DetectionAlgorithm Based on the Vectorcardiogram for Accurate Gated MagneticResonance Acquisitions,”Magnetic Resonance in Medicine,42:361-370(1999)。将ECG电极施加于患者、并且建立以获取ECG数据相对复杂且费时。并且对于能维持下去的ECG信号,在成像期间重新定位电极可以变得有必要,并且这典型地要求MR图像采集被停止且患者被从MR单元的孔中取出。
使MR图像同步在分段或交错的电影成像中尤其具有价值,其中每个图像的数据源自所感兴趣运动的不同周期。为了使图像具有意义,只有来自不同周期的对应阶段的数据应该被加以组合,因此需要使成像数据与所述运动同步。例如,在分段的心脏电影成像中,每个图像的k-空间(k-space)线可来自在单次呼吸保持期间存在的14个或更多个不同的心动周期。
分段或交错的电影成像的同步既可以预期地也可以回顾地实现。当被预期地实现时,响应于源自ECCI电极或作为选择地源自手指脉冲血氧计或被设计成监视与心动周期同步的生理信号的其它设备的时序信号,成像数据的采集开始。数据采集典型地持续一固定的时间间隔,其足够长以覆盖收缩阶段和舒张阶段的开始。然后,典型地存在一静止期间直至下一个时序信号。作为选择地,通过与ECCI基时序信号异步地连续获取成像数据,可以回顾地实现同步,同时每行数据相对于最后触发信号被获取的时间得到记录。采集之后,基于所记录的时序数据,数据被分配给心动周期的适当阶段。
在现有技术中已经尝试直接地从MR数据得到时序信息,以为了消除ECG或其它附加的时序测量的需要。但是这些努力一直依赖于收集在产生MR图像时并不被使用的附加MR数据。
例如,Spraggins的U.S.专利4,961,426和Spraggins的WirelessRetrospective Gating:Application to Cine Cardiac Imaging,MagneticResonance Imaging,8:675-681(1990)讲授获取从中可以得到时序信号的附加“时序片”。以未进行相位编码的回波形式的时序数据与成像数据采集相交错(每隔采集的是时序数据),并且可以从与正在被成像区域不同的心脏区域(例如其中运动更加可见的区域)获取。Kim等人所著的‘Extraction of Cardiac and Respiratory Motion Cycles by Use ofProjection Data and Its Application to NMR Imaging,’MagneticResonance in Medicine,13:25(1990)使用类似的方案,只是附加数据被变换成空间域(Spragg ins已经直接地使用频率域数据)以提供表示图像片层到沿着时间数据采集方向被定向的线上的投影信号。
另一方案可见于Vasanwala等人的“Prospective MR Signal-BasedCardiac Triggering”,Magnetic Resonance in Medicine,42:82-86(1999),其中特殊的“triggering sequence”被用来获取表示主动脉血速的经速度编码的数据。当发现触发事件时,系统从触发序列转换到成像序列。
在呼吸门控领域,被公知为导航器门控或导航器回波的技术从额外的非成像数据得到时序信号。Ehman & Felnilee,“Adaptive technique forhigh-definition MR1 of moving structures.”Radiology,173:255-263(1988)。典型地,垂直于隔膜的投影被获取,同时边缘探测算法被用来确定呼吸周期的位置。
发明内容
总体上,在第一方面,本发明的特点是:通过从MR成像数据抽取时序信号而不是仅仅依赖仅为了时序所获取的附加数据,使MR成像数据与患者的运动(例如心脏的跳动)同步。通过从MR成像数据得到时序信息,优良的图像质量是可能的。例如,在其中图像是基于在一次呼吸保持期间所收集的数据的心脏成像中,在所述呼吸保持期间更多的可用时间可用于收集成像数据。
优良的图像质量还可由MR数据与影响它的运动的直接同步,而不是以外部生理信号或在图像生成中未被使用的MR信号形式的一些非直接测量所产生。
临床生产力得到增加,因为数据收集时间被缩短,以及因为使患者做好MR研究的准备需要更少的时间(因为没有必要缚上ECG电极)。
本发明解决了在MR单元的敌对环境中获取ECG信号的问题。无能力获取可靠的ECG信号是失败的心脏MR检查的常见原因。不需要昂贵且复杂的设备(例如ECG监视设备)来产生时序信息。
宽范围的身体运动可以被同步,包括随意运动象肌肉收缩或咀嚼、以及非随意运动象呼吸。呼吸运动是胸及腹部MR图象中的假象和差图像质量的主要原因。呼吸运动信息可以直接地从MR数据中抽取并且可以被用来使图像数据与呼吸周期的静止期间同步,从而避免了假象。这允许难以控制其呼吸的患者(例如老年患者、婴儿)在自由呼吸期间使数据得到收集,从而避免为了胸或腹部的MRI对呼吸保持的常见要求。
时序信号可以源自胎儿的MR数据,以避免测量胎儿ECG的复杂性,从而使能高时间和空间分辨率的胎儿心脏MR图像的采集,这是迄今为止还一直不可能的事情。
本发明第一方面的优选实施可结合下述中的一个或更多个:
可沿着径向或螺旋k-空间轨迹获取成像数据,这样时序信息可从处于或接近原点的被频繁收集到的k-空间点中抽取。取决于抽取的方法,时序信息可从原始的k-空间数据或从被变换成空间域的k-空间中获取。时序信息可基于k-空间的中心点,其被公知为原始数据的回波峰值。它还可基于通过将原始MR数据变换成空间域而计算出的1维投影或2维图像。多于一个投影可被用来允许图像质心的计算。投影可到一k-空间线上,其取向被选择成增强对所感兴趣的运动的敏感性。时序信息可基于低分辨率图像的相关,所述低分辨率图像可通过使用交错数据采集(例如使每个交错包括覆盖分散的k-空间区域的一组k-空间轨迹)而被获取。时序信息可从低分辨率图像的被选择区域中被抽取。
被抽取的时序信息可被处理以提供与运动的时间对应(例如表示运动开始的时间或在运动期间另一事件时间的时间值)。处理可包括抽取信号的峰值、相位或时间变化率。
时序信息可被用来回顾地或预期地使MR成像数据与所述运动同步。患者的运动可是周期性的(例如,心脏或肺的周期性运动)。时序信息可包括随着运动期间数值发生变化的时间变化信号。MR成像数据可是分段的电影成像数据。
时序信息可独自地从MR成像数据中被抽取或它可从MR成像数据与附加非成像数据的组合中被抽取。
所述方法可通过使用被定位到正在移动的身体部分上的RF线圈(例如RF线圈可被定位到心脏上)而执行。
总体上,在第二方面中,本发明的特点在于MR成像方法,其包括将产生RF信号的脉冲序列同时施加到至少第一和第二RF线圈上;处理在所述第一线圈上所获取的RF信号,以抽取MR成像数据;处理在所述第二线圈上所获取的RF信号,以抽取除MR成像数据以外的数据;以及使用MR成像数据来产生MR图像。
在本发明这个方面的优选实施中,可结合下述中的一个或更多个:本发明可被应用到使MR图像数据与患者的运动同步,在所述情况下时序数据从第二RF线圈中被获取、指示运动的时序信息被从时序数据中抽取、以及时序信息被用来使MR成像数据与所述运动同步。可存在形成线圈阵列的多个第一RF线圈,并且每个被确定大小且被放置主要地用于获取MR成像数据。第二RF线圈可被定位到正在运动的身体部分上。MR成像数据可包括心脏成像数据并且第二RF线圈可被定位到心脏附近的胸部区域。时序信息专门地自从第二RF线圈所获取的时序数据中被抽取。MR成像数据专门地自从第一RF线圈所获取的RF信号中被抽取。
总体上,在第三方面,本发明的特点在于:通过从MR数据中而不是从ECG其它非信号中抽取时序信息,使MR胎儿心脏成像数据与胎儿心脏的运动同步。
在本发明这个方面所优选的实施中,可结合下述中的一个或更多个:时序信息可从MR成像数据中被抽取(即相同于本发明的第一方面)。作为选择地,时序信息可从并不作为MR成像数据使用的MR数据中被抽取。
从附图及详细说明以及从权利要求中,本发明的其它特点和优点将是显而易见的。
附图说明
图1是描绘由本发明的许多实施所遵循的过程的方框图。
图2是示出径向k-空间采集的k-空间图。
图3是示出螺旋k-空间采集的k-空间图。
图4和5是k-空间图,其每个均示出径向k-空间轨迹的不同组(交错)。
图6是示出回顾地被同步的成像数据采集的一部分的时序图。
图7是用在本发明一些实施中的RF接收线圈的示意图。
图8示出用在一些实施中用于提供径向k-空间轨迹的一个可能的脉冲序列。
具体实施方式
本发明存在可能的很多不同实施,太多以致于不可能在此都加以说明。下面对目前所优选的一些可能的实施加以说明。然而,无论怎样强烈地强调也不过分,这些是对本发明实施的说明,且不是对本发明的说明,其并不被限制于在这个段中所说明的详细实施而是在权利要求中以更广泛的术语被加以说明。
图1示出许多这些实施所共同的步骤。MR成像数据被获取。时序信息被从MR成像数据中抽取。时序信息被加以处理以提供与所感兴趣的运动的时间对应。成像数据被与运动同步。被同步的MR成像数据被用来产生MR图像。
MR成像数据以这样的方式被获取,以便于允许时序信息的抽取反映患者的运动。优选地这是通过沿着频繁地通过k-空间中心(或原点)的k-空间轨迹获取数据而完成的,因为在k-空间原点处的值包含表示空间域图像DC(或平均)值的回波峰值。具有图像的平均值通常对于抽取时序信息是非常有益的,因为身体运动趋向于与平均值的变化相关联。几乎任何MR图像采集将包含通过k-空间中心的一些k-空间轨迹,但是总体上所述轨迹太不频繁地通过所述中心以致于无益于获取时序信号。
存在各种数据采集策略,其可以被实施以确保k-空间原点数据被频繁地得到收集。例如,成像数据可以沿着径向或螺旋k-空间轨迹得到收集(图2和3),因为这样做固有地引起数据在k-空间原点10处被频繁地加以收集。用在所谓的Propeller ME中的同心矩形条也可是有益的。Pipe,JamesU.,“Motion Correction with PROPELLER MRI:Application to Head Motionand Free Breathing Cardiac Imaging,”Magnetic Resonance in Medicine42:963-969(1999)。并且甚至笛卡尔(Cartesian)轨迹可被采用以便于数据可能足够经常地从k-空间的中央区域被获取,以允许有益的时序信号的抽取。
没有必要每个数据采集轨迹通过k-空间原点,但是优选地这相对频繁地发生在数据采集期间,这样所抽取的时序信息提供高到足以表示所感兴趣的身体运动的时间分辨率。同样,并不绝对必要的是数据采集正好通过k-空间原点,因为有可能时序信息被从接近但是相隔于原点的k-空间位置被抽取。
优选地,用于获取成像数据的k-空间轨迹靠其自身足以提供时序信息,但是在一些实施中,时序信息主要从k-空间成像数据中、并且还从并不用于成像的一些附加k-空间数据中被抽取。
存在许多可能的用于抽取时序信息的具体技术。例如,时序可以从下述中被抽取:(1)原始的k-空间中心点;(2)1维空间域投影;(3)多个1维空间域投影(记波图);(4)低分辨率的2维图像。在所述第一实例中,时序信息直接从原始的k-空间数据中被获取。在其它实例中,时序信息从被变换的空间域数据中被抽取。通过使用原始的k-空间数据以及被变换的数据,许多其它技术是可能的。
第一实例简单地从k-空间的中心点抽取时序信息。任何频繁地通过k-空间中心的k-空间轨迹(例如,径向的或螺旋的)将以这个方案工作。
第二和第三实例对沿着径向k-空间轨迹所获取的数据加以处理以产生沿着这些轨迹的1维图像投影。时序信息可以从单1维投影(第二实例)或从来自处于变化角度的k-空间线所取得的多个投影(第三实例)中被抽取。计算沿着k-空间线所采集的k-空间数据的傅立叶变换给出沿着垂直于所述线方向的图像投影。使用单投影可以提供沿着所选择方向的质心变化,其应该考虑到正在被检查的解剖体被加以选择。时序信息可以从一个这样的1维投影中(第二实例)、或从每个对应于沿着不同角度的径向线的几个I维投影中(第三实例)被抽取。将来自处于不同角度径向线的投影进行组合允许人们计算可提供时序信息的图像片层的两维质心。
在第四实例中,数据采集是利用提供一系列低分辨率2维图像的交错(螺旋或径向的)k-空间轨迹而执行的,所述2维图像与参考图像进行比较(例如,被交叉相关),并且时序信息是基于所述比较的结果(例如,所计算的相关系数的时间变化可以提供时序信息)。低分辨率的图像可以是有可能包括所感兴趣的移动目标的整个图像或子区域。图4和5是交错的径向k-空间轨迹的实例。在图4中被标注为A至D的径向轨迹集合是经过至少一个完整的运动周期被重复地获取,紧跟其后的是图5中被标注为E至H的径向轨迹集合,其也是经过至少一个运动的满周期被重复地获取(紧接其后的是许多更多的交错的轨迹集合)。低分辨率的2维图像可以从每个轨迹集合(例如从图4中的A至1集合及从图5中的B至H集合)中得到。作为选择地,如果螺旋k-空间轨迹被使用,则每个螺旋轨迹可以靠其自身提供低分辨的2维图像,从所述图像可以抽取时序信息。
在时序信息被抽取之后,信息被进行处理以提供与所感兴趣的运动的时间对应。许多不同的处理技术可以被使用。在一些实施中时序信息可是对应于运动的时间变化信号,并且可以使用算法来识别在信号中对应于所感兴趣的时间事件(或所感兴趣运动的特定期)的特征。但是在通常意义上时序信息作为时间变化信号存在并不是必要的,因为用于同步所需要的信息可仅是针对每个正在研究的运动周期的单个时间值,且因此所述算法可只是处理时序信息,以产生输出时间值,而无需存在时间变化信号的中间步骤。在一些实施中,在时序信息上所执行的处理的结果还可是被用来启动图像数据采集的触发信号。
在其中时序信息是时间变化信号的情况下,信号的各种特征---例如,量值、相位或变化率对于使相联系的成像数据同步是有益的。在心脏脊成像中,时间变化信号可是表示流进及流出由MR单元正在监视的组织片层的血量。在典型的短轴或长轴方位上,周期性的波动由心室血库量的变化而产生。通过使用利用FALSE或TueFISP类型序列的采集技术,时间变化信号在心动周期的舒张期的结束可具有最大的量值,且在收缩期的结束可具有最小的量值。因此,信号的量值可直接有益于确定心脏处于心动周期的某个阶段的时间。
心脏成像还提供时间变化信号的相位是有益的一实例。取决于成像取向,连同在RF激励之后所施加的梯度,信号的相位可与循环血库的速度变化相关联。这可使相位可能被用于将时间标记分配到心脏成像数据中。
信号的量值及相位已经提到,但是在一些实施中,信号的实或虚的组分可直接有益于同步化。
时间变化信号的变化率还可是有价值的。例如,信号中的最大或最小变化率可有益于使所收集的成像数据同步于心动周期、呼吸周期、或患者的其它不随意运动或随意运动中的特定点。
在一些实施中,可在时间变化信号上执行各种信号调节过程,以增强其用于使成像数据同步的价值。例如,在将信号通过峰值探测算法以确定信号的最大量值之前,可使信号通过一低通或带通滤波器。
在时序信号已经被加以处理以提供所希望的与感兴趣运动的时间对应之后,成像数据被与所述运动同步。如上面所说明,所述处理可产生(例如运动开始的、或在运动期间另一阶段或时间事件的)时间值,或它可产生一些类型的触发信号。使成像数据同步可以被预期地(例如通过使用触发信号输出以启动成像数据的采集)或回顾地(例如,通过比较运动开始的时间值及与成像数据相联系的时间标记)完成。
回顾地被同步的实施被示意性地示例于图6。数据以分段型式被获取,这样在给定的采集窗内,对应电影图像帧的仅一部分k-空间线被获取。在图6所示出的实例中,在被建立成至少与最长的预计心动周期一样长的采集窗期间,每个k-空间段被采集12次。每个段具有11个k线,且数据在发生于同一呼吸保持中的15个采集窗上被采集(在图中仅示出两个周期)。对于每个图像帧的总计165个k-空间线,每个图像帧由每个11个k-空间线的15个段组成。段被连续地获取,并且在数据采集完成之后,从成像数据所抽取的时序信息被用来回顾地使被获取的数据同步。为了将成像数据同步于心动周期,所获取的成像数据的时间标记(例如到每个所获取的线)建立了相对于相联系的心动周期开始的时间获取数据的时间(那个时间已经从自成像数据被抽取的时序信息中得到)。
在预期性地被触发的实施中,时序信息的抽取和处理被进行得足够快,以致于提供实时的或接近实时的触发信号来影响成像数据的采集。例如,在心脏分段电影成像情况下,当触发信号指示已经探测到R-波时成像数据段的采集可开始。随着预期的触发,触发启动k-空间的一个段的采集。在电影成像中,数据的一个或更多个线的这个相同段被重复地需要,直至下一个触发发生。在那个点上,采集切换到k-空间数据的下一段。
传统的MR单元可被用来获取MR成像数据,使数据采集如上所讨论沿着k-空间轨迹。一全身RF发射线圈可被用来发射RF脉冲,并且四至八个RF接收线圈的表面阵列可被用来获取RF信号(图7示出四个这样的表面阵列线圈)。众所周知的梯度线圈可被使用。MR单元可以使用极宽种类的脉冲序列,包括图8中所示的一个(其产生径向的相隔k的轨迹)。能够提供所希望的k-空间轨迹并且能够足够快地获取成像数据以满足时序要求的任何序列将足够。例如,对于心脏的同步,可希望地是至少大约每10毫秒获取调谐信息,并且因此能够以那个速率获取数据的任何序列将足够。
典型地是在胸或腹的MRI中使用多个独立RF(射频)线圈元件的阵列(如图7中示意性所示)以优化图像的信噪比。这些线圈的每个在图像视场的一个区域贡献信号,并且线圈的敏感区域典型地彼此重叠。来自多个线圈的信号被独立的接收器信道获取,并且被独立地重建。来自多个线圈的信息组合典型地发生在图像域而不是在原始数据中。虽然阵列中的所有线圈可向图像贡献有用的信息,但是与心脏运动相关的图像数据的变化将在直接被放置到心脏上的线圈中是最强的。这个线圈的敏感区域将由来自心脏的信号支配,并且因此受其运动影响最大。因此,有利地是选择仅一个成像线圈(使选择基于因运动导致的信号变化的幅度被自动进行)作为成像数据源,从此成像数据抽取时序信息。
作为选择地,还可有利地是将附加的小局部H线圈放置在正在被成像的解剖体(例如在用于心脏成像的患者胸上)近旁,并且从由那个线圈所接收的信号中抽取时序信息。附加的线圈在大小可是比用于获取成像数据的线圈较小或较大,但是可能有利地是使所述附加线圈更小。这将趋向于改善运动与被抽取的时序信息的相关。在这种情况下,来自附加线圈的数据虽然从相同的成像脉冲序列产生且基本上与由线圈阵列所获取的成像数据被同时获取,但其并不用于成像。这样的实施是本发明第二方面的一些可能的实施。
被定位的RF线圈-不管它是还用于成像的线圈之一还是未被用于成像的附加定位线圈-充当将所获取的成像数据限制到仅包含心脏(或其它运动的目标)的小局部区域,由此降低对与所希望的时序信息无关的噪声源的敏感性。
除了上述说明的那些以外的本发明的许多其它实施在本发明的范围内,其由随后的权利要求来定义。当权利要求提到“MR成像数据”时,它们是指被用来产生MR图像的数据,且不是指排他地为其它目的而获取的数据,例如仅为提供时序信息而获取的附加数据。
本发明一个方面所要求的特征并没有必要要求去实践本发明的其它方面。例如,在其最广泛的方面,没有必要抽取时序信息以实践本发明的第二方面,并且在实施中当时序信息被抽取时它从M成像数据中被抽取是没有必要的(例如它可以抽取自从上面说明的附加被定位线圈所获取的时序数据)。进一步的实例是:实践本发明的第三方面并不要求从MR成像数据中获取时序信息。
本发明第一方面的其它实施的几个实例包括如下:
时序信息并不需要被排他地从成像数据中抽取,因为通过主要从成像数据中抽取时序信息但与仅为了时序目的而得到的附加数据相结合,可以拥有来自本发明的一些益处。
多于一个技术可被用于抽取时序信息,并且基于在取得所希望的与运动的时间对应时技术的表现,来选择采用哪个技术,或所抽取的时序信息可是两个技术的结果的加权组合。
利用数据选择性的滤波,心脏和呼吸的时序信息可典型地从相同的MR成像数据中被抽取。
本发明的实施可与用于取得时序信息的其它技术相组合。例如,呼吸运动可利用本发明的实施被同步化,虽然心脏成像数据的同步由传统的基于ECG时序来处理。

Claims (37)

1.一种使MR成像数据与患者的运动同步的设备,包括:
用于在所述运动期间从所述患者获取MR成像数据的装置,该装置包括:
用于为多个图像帧获取图像数据的装置;和
用于从通过k-空间中心的多个k-空间轨迹获取用于一帧的数据的装置;
用于从至少一部分所述MR成像数据中抽取时序信息的装置,该抽取时序信息的装置包括用于从多个k-空间轨迹获取用于一帧的时序信息的装置,该时序信息指示所述运动;
用于响应于所述时序信息而使所述MR成像数据与所述运动同步的装置;以及
用于响应于被同步的MR成像数据而产生MR图像的装置。
2.根据权利要求1的设备,其中所述k-空间轨迹包括径向或螺旋轨迹。
3.根据权利要求2的设备,其中所述抽取时序信息的装置包括用于从原始k-空间数据中抽取时序信息的装置。
4.根据权利要求2的设备,其中所述抽取时序信息的装置包括用于从被变换成空间域的k-空间数据中抽取时序信息的装置。
5.根据权利要求3的设备,其中所述获取时序信息的装置包括用于从位于k-空间的中心的位置获取时序信息。
6.根据权利要求4的设备,其中所述抽取时序信息的装置包括:
用于基于一个或多个被计算的1维投影而抽取时序信息的装置;和
用于变换沿着径向k-空间轨迹而收集的原始k-空间数据的装置。
7.根据权利要求6的设备,其中所述抽取时序信息的装置包括用于基于多个被计算的1维投影而抽取时序信息的装置。
8.根据权利要求7的设备,其中所述抽取时序信息的装置包括用于基于被计算的图像的近似质心而抽取时序信息的装置。
9.根据权利要求6的设备,其中所述抽取时序信息的装置包括用于基于计算机1-维投影而抽取时序信息的装置,该投影基于一k-空间线,所述k-空间线的取向被选择成增强所述时序信息对运动的敏感性。
10.根据权利要求4的设备,其中所述抽取信息的装置包括用于基于使用交错数据采集所获取的低分辨率图像的相关性而抽取时序信息的装置,其中该交错包括覆盖分散的k-空间区域的k-空间轨迹组。
11.根据权利要求1的设备,其中所述抽取时序信息的装置包括用于抽取在运动进程中值发生变化的时间变化信号的装置。
12.根据权利要求1的设备,还包括:
用于处理所述时序信息的装置;和
用于响应于所述处理装置以提供与所述运动的时间对应的装置。
13.根据权利要求12的设备,其中所述提供时间对应的装置包括用于提供表示运动开始时间或在运动期间另一事件时间的时间值的装置。
14.根据权利要求12的设备,其中所述时序信息包括时间变化信号,该时间变化信号的值在运动阶段发生变化。
15.根据权利要求14的设备,其中所述提供时间对应的装置包括用于抽取该时间变化信号的峰值、相位或速率的装置。
16.根据权利要求1的设备,其中所述用于同步的装置包括回顾地使MR成像数据与运动同步。
17.根据权利要求1的设备,其中所述MR成像数据包括分段的电影成像数据。
18.根据权利要求1的设备,还包括被定位到患者的正在移动的部分的RF线圈。
19.根据权利要求18的设备,其中所述患者的正在移动的部分是心脏,并且该RF线圈被定位在心脏。
20.根据权利要求18的设备,其中所述RF线圈是RF接收线圈,且一个单独的RF发射线圈包括用于递送RF激励的装置。
21.根据权利要求1的设备,其中所述患者的运动是周期性的。
22.根据权利要求21的设备,其中所述运动包括心脏的周期性运动。
23.根据权利要求21的设备,其中所述运动包括在呼吸期间肺的周期性运动。
24.根据权利要求21的设备,其中所述运动包括:
在呼吸期间的肺的周期性运动;和
心脏的周期性运动;并且
所述用于抽取的装置包括用于为多个运动选择性地对成像数据进行滤波的装置。
25.根据权利要求1的设备,其中所述用于同步的装置包括响应于所述时序信息而单独地同步所述患者的多个运动的装置。
26.根据权利要求25的设备,其中所述用于同步的装置包括:
用于在呼吸期间同步肺的周期性运动的装置;和
用于同步心脏的周期性运动的装置,其中所述用于同步心脏的周期性运动的装置响应于ECG信号。
27.根据权利要求1的设备,其中所述抽取时序信息的装置包括用于抽取包含MR成像数据和非成像数据的时序信息的装置。
28.根据权利要求27的设备,其中所述MR成像数据包括所述时序信息的主源。
29.根据权利要求1的设备,其中所述用于抽取时序信息的装置包括用于仅从MR成像数据中抽取时序信息的装置。
30.一种用于使MR成像数据与患者的运动同步的设备,包括:
用于施加一个脉冲序列的装置;
用于在至少第一和第二RF线圈上同时获取RF信号的装置,其中该第二RF线圈包括用于获取时序数据的装置;
用于响应于在所述第一线圈上所获取的RF信号以获取MR成像数据的装置;
用于响应于在所述第二线圈上所获取的RF信号以获取时序数据的装置;
用于响应于所述时序数据而抽取指示所述运动的时序信息;
用于响应于所述时序信息而使所述MR成像与所述运动同步的装置;以及
用于响应于所述MR成像数据而产生MR图像的装置。
31.根据权利要求30的设备,还包括用于在脉冲序列期间防止所述第一和第二RF线圈发射RF激励的装置。
32.根据权利要求30的设备,其中所述第一RF线圈包括其中包含用于获取MR成像数据的装置的多个第一RF线圈的阵列。
33.根据权利要求30的设备,其中所述第二RF线圈被定位到身体的正在运动的部分上。
34.根据权利要求33的设备,其中所述MR成像数据包括心脏成像数据,并且所述第二RF线圈被定位到心脏附近的区域。
35.根据权利要求30的设备,其中所述抽取时序信息的装置包括用于从自所述第二RF线圈所获取的时序数据抽取时序信息的装置。
36.根据权利要求30的设备,其中所述获取MR成像数据的装置包括从自所述第一RF线圈所获取的RF信号抽取数据的装置。
37.根据权利要求1的设备,其中
所述用于获取MR成像数据的装置包括用于从胎儿获取MR数据的装置;和
用于从对应于胎儿心脏的运动的MR成像数据的至少一部分抽取时序信息的装置。
CN2007101938423A 2003-02-06 2004-02-06 用于使磁共振成像数据同步于身体运动的设备 Expired - Fee Related CN101173978B (zh)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
US10/360,530 US6798199B2 (en) 2003-02-06 2003-02-06 Method for synchronizing magnetic resonance imaging data to body motion
US10/360530 2003-02-06

Related Parent Applications (1)

Application Number Title Priority Date Filing Date
CNA2004100038113A Division CN1557253A (zh) 2003-02-06 2004-02-06 用于使磁共振成像数据同步于身体运动的方法

Publications (2)

Publication Number Publication Date
CN101173978A CN101173978A (zh) 2008-05-07
CN101173978B true CN101173978B (zh) 2012-08-08

Family

ID=32824027

Family Applications (2)

Application Number Title Priority Date Filing Date
CNA2004100038113A Pending CN1557253A (zh) 2003-02-06 2004-02-06 用于使磁共振成像数据同步于身体运动的方法
CN2007101938423A Expired - Fee Related CN101173978B (zh) 2003-02-06 2004-02-06 用于使磁共振成像数据同步于身体运动的设备

Family Applications Before (1)

Application Number Title Priority Date Filing Date
CNA2004100038113A Pending CN1557253A (zh) 2003-02-06 2004-02-06 用于使磁共振成像数据同步于身体运动的方法

Country Status (4)

Country Link
US (1) US6798199B2 (zh)
JP (1) JP2004237109A (zh)
KR (1) KR20040071643A (zh)
CN (2) CN1557253A (zh)

Families Citing this family (51)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6724923B2 (en) * 2001-04-13 2004-04-20 Ge Medical Systems Global Technology Co., Llc Automatic coil selection of multi-receiver MR data using fast prescan data analysis
DE10247299A1 (de) * 2002-10-10 2004-04-22 Philips Intellectual Property & Standards Gmbh Bildverarbeitungseinheit und Verfahren für die Zuordnung von gespeicherten zu aktuellen Aufnahmen
WO2005039253A1 (en) * 2003-10-16 2005-04-28 Philips Intellectual Property & Standards Gmbh Device and method for providing an angiographic image
US7844317B2 (en) * 2003-11-26 2010-11-30 General Electric Company Method and system for estimating three-dimensional respiratory motion
US20050113673A1 (en) * 2003-11-26 2005-05-26 Avinash Gobal B. Method and system to reduce motion-related image artifacts during breath holding
US7756565B2 (en) * 2003-11-26 2010-07-13 General Electric Company Method and system for composite gating using multiple inputs
US8064983B2 (en) * 2003-11-26 2011-11-22 General Electric Company Method and system for prospective gating using multiple inputs
US7367953B2 (en) * 2003-11-26 2008-05-06 Ge Medical Systems Global Technology Company Method and system for determining a period of interest using multiple inputs
US8285359B2 (en) * 2003-11-26 2012-10-09 General Electric Company Method and system for retrospective gating using multiple inputs
US9301704B2 (en) 2004-03-26 2016-04-05 Kabushiki Kaisha Toshiba Magnetic resonance imaging system for non-contrast MRA and magnetic resonance signal acquisition method employed by the same
JP4820567B2 (ja) * 2004-03-26 2011-11-24 株式会社東芝 磁気共鳴イメージング装置及び磁気共鳴信号の収集方法
US7693569B1 (en) * 2004-10-12 2010-04-06 General Electric Company Method and system of determining motion in a region-of-interest directly and independently of k-space trajectory
US7477771B2 (en) * 2004-12-03 2009-01-13 Genral Electric Company Method and system for extracting information about the cardiac cycle from CT projection data
JP4612000B2 (ja) * 2005-01-20 2011-01-12 株式会社日立メディコ 磁気共鳴イメージング装置
US7574249B2 (en) * 2005-02-08 2009-08-11 General Electric Company Device-less gating of physiological movement for improved image detection
JP4950466B2 (ja) * 2005-09-22 2012-06-13 株式会社東芝 磁気共鳴イメージング装置
KR100735253B1 (ko) 2005-10-11 2007-07-03 삼성전자주식회사 디지털 방송 수신 단말기의 디지털 방송 화면 캡쳐 방법
EP1948022A2 (en) * 2005-11-10 2008-07-30 Philips Intellectual Property & Standards GmbH Signal processing unit for producing images
KR100726739B1 (ko) 2006-01-16 2007-06-11 한국과학기술원 방사형 케이-공간 경로법을 이용한 병렬 자기 공명 영상획득 방법, 그 장치 및 이를 실행하기 위한 컴퓨터로 읽을수 있는 기록매체.
JP2008178592A (ja) * 2007-01-25 2008-08-07 Ge Medical Systems Global Technology Co Llc 磁気共鳴イメージング装置、スキャン装置、磁気共鳴イメージング方法、および、そのプログラム
DE102007009182B4 (de) * 2007-02-26 2016-09-22 Siemens Healthcare Gmbh Verfahren und Vorrichtung zur Bilddarstellung von sich zyklisch bewegenden Objekten
US8155389B2 (en) * 2007-10-02 2012-04-10 The University Of Utah Research Foundation Method and system for motion correction in imaging systems
US8170307B2 (en) 2008-09-23 2012-05-01 The Methodist Hospital Automated wall motion quantification in aortic dissections
DE102009027356A1 (de) * 2009-05-28 2010-12-02 Sirona Dental Systems Gmbh Bildgebendes System zur Erzeugung eines 3D-Datensatzes sowie Verfahren zu dessen Betrieb
JP5569714B2 (ja) * 2009-06-16 2014-08-13 株式会社国際電気通信基礎技術研究所 磁気共鳴画像化装置および磁気共鳴画像化方法
JP5388749B2 (ja) * 2009-08-11 2014-01-15 株式会社東芝 磁気共鳴イメージング装置
US9008753B2 (en) 2009-11-10 2015-04-14 Deutsches Herzzentrum Berlin Look-locker IR-SSFP for cardiac MR imaging with simultaneous generation of cardiac T1 maps, cine images and IR-prepared images
JP5858791B2 (ja) * 2010-01-22 2016-02-10 株式会社日立メディコ 磁気共鳴イメージング装置及び血管画像撮像方法
US9241645B2 (en) 2010-10-13 2016-01-26 Kabushiki Kaisha Toshiba Multiple MR fluid flow imaging at predetermined temporal resolution within selected period of cardiac cycle determined by multiple MR imaging at different temporal resolution
US8543187B2 (en) * 2010-10-13 2013-09-24 Kabushiki Kaisha Toshiba Non-contrast MRI signal acquisition
DE102011085399A1 (de) * 2011-10-28 2013-05-02 Siemens Aktiengesellschaft Verfahren zur Aufnahme von Magnetresonanzbilddaten bei Verwendung einer Beatmungsvorrichtung
KR101845222B1 (ko) 2011-11-09 2018-04-05 삼성전자주식회사 고차 확산자기공명 영상의 인공물 저감 장치 및 방법
CN102540116A (zh) * 2011-12-08 2012-07-04 中国科学院深圳先进技术研究院 磁共振成像方法和系统
CN102540115A (zh) * 2011-12-12 2012-07-04 中国科学院深圳先进技术研究院 磁共振成像方法及系统
GB201121406D0 (en) 2011-12-13 2012-01-25 Isis Innovation Systems and methods
GB201121404D0 (en) 2011-12-13 2012-01-25 Isis Innovation Diagnosis
CA2866089C (en) 2012-03-01 2020-12-29 Syracuse University Enhanced electronic external fetal monitoring system
US9002079B2 (en) * 2012-05-21 2015-04-07 General Electric Company Systems and methods for motion detecting for medical imaging
WO2013190451A1 (en) * 2012-06-21 2013-12-27 Koninklijke Philips N.V. Magnetic resonance examination system with motion detection
KR101967240B1 (ko) * 2012-10-10 2019-04-09 삼성전자주식회사 피검체의 움직임을 보정하여 mri 영상을 획득하는 방법 및 mri 장치
GB201304728D0 (en) 2013-03-15 2013-05-01 Isis Innovation Medical imaging
US10261160B2 (en) * 2013-08-07 2019-04-16 Koninklijke Philips N.V. ECG-based triggering for magnetic resonance imaging
US9835705B2 (en) * 2013-11-07 2017-12-05 Beth Israel Deaconess Medical Center, Inc. System and method for free-breathing volumetric imaging of cardiac tissue
CN105793722B (zh) 2013-12-02 2019-06-28 皇家飞利浦有限公司 用于稳态mr序列的实时自适应生理同步和门控
JP6489744B2 (ja) * 2014-02-28 2019-03-27 ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー 磁気共鳴装置およびプログラム
US9702956B2 (en) * 2014-03-25 2017-07-11 Beth Israel Deaconess Medical Center, Inc. (Bidmc, Inc.) MRI methods and apparatus for flexible visualization of any subset of an enlarged temporal window
DE102014225282B4 (de) * 2014-12-09 2016-07-21 Siemens Healthcare Gmbh Deformationsberechnung bei zyklischer Bewegung eines Untersuchungsobjekts
DE102015205937A1 (de) * 2015-04-01 2016-10-06 Siemens Healthcare Gmbh Erfassen von physiologischen Signalen und Bilddaten
EP3470869A1 (en) * 2017-10-10 2019-04-17 Koninklijke Philips N.V. Mr imaging using motion-dependent radial or spiral k-space sampling
US11119176B2 (en) 2017-11-27 2021-09-14 Siemens Healthcare Gmbh Automatic capture of cardiac motion by pre-scan and automated data evaluation for determination of motionless periods within the RR-interval
US20230181041A1 (en) * 2020-04-27 2023-06-15 Children's Medical Center Corporation Fetal cardiac mri using self-gating with a cartesian k-space trajectory

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4961426A (en) * 1988-08-19 1990-10-09 Siemens Medical Systems, Inc. Method for retrospectively gating NMR data
CN1151858A (zh) * 1995-12-15 1997-06-18 通用电器横河医疗系统株式会社 扩散敏化成像方法和磁共振成像装置
JP2646663B2 (ja) * 1988-06-07 1997-08-27 株式会社日立製作所 動体イメージング方法およびその装置
US6073041A (en) * 1996-07-18 2000-06-06 Regents Of The University Of Minnesota Physiological corrections in functional magnetic resonance imaging
US6461315B1 (en) * 1998-03-16 2002-10-08 Siemens-Elema Ab Apparatus for improving the distribution of gas in the lungs of a patient receiving respiratory treatment
CN1380983A (zh) * 2000-03-24 2002-11-20 皇家菲利浦电子有限公司 利用子采样的磁共振成像方法

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4728890A (en) * 1985-08-16 1988-03-01 Picker International, Inc. Motion artifact suppression technique of magnetic resonance imaging
US4710717A (en) * 1986-12-29 1987-12-01 General Electric Company Method for fast scan cine NMR imaging
US6539074B1 (en) * 2000-08-25 2003-03-25 General Electric Company Reconstruction of multislice tomographic images from four-dimensional data
US6644976B2 (en) * 2001-09-10 2003-11-11 Epoch Innovations Ltd Apparatus, method and computer program product to produce or direct movements in synergic timed correlation with physiological activity
US6922580B2 (en) * 2002-06-04 2005-07-26 Koninklijke Philips Electronics N.V. Blood flow gated MRI

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2646663B2 (ja) * 1988-06-07 1997-08-27 株式会社日立製作所 動体イメージング方法およびその装置
US4961426A (en) * 1988-08-19 1990-10-09 Siemens Medical Systems, Inc. Method for retrospectively gating NMR data
CN1151858A (zh) * 1995-12-15 1997-06-18 通用电器横河医疗系统株式会社 扩散敏化成像方法和磁共振成像装置
US6073041A (en) * 1996-07-18 2000-06-06 Regents Of The University Of Minnesota Physiological corrections in functional magnetic resonance imaging
US6461315B1 (en) * 1998-03-16 2002-10-08 Siemens-Elema Ab Apparatus for improving the distribution of gas in the lungs of a patient receiving respiratory treatment
CN1380983A (zh) * 2000-03-24 2002-11-20 皇家菲利浦电子有限公司 利用子采样的磁共振成像方法

Also Published As

Publication number Publication date
JP2004237109A (ja) 2004-08-26
CN1557253A (zh) 2004-12-29
KR20040071643A (ko) 2004-08-12
US20040155653A1 (en) 2004-08-12
CN101173978A (zh) 2008-05-07
US6798199B2 (en) 2004-09-28

Similar Documents

Publication Publication Date Title
CN101173978B (zh) 用于使磁共振成像数据同步于身体运动的设备
US11026636B2 (en) Method for generating a medical data set of a moving body part
Bartula et al. Camera-based system for contactless monitoring of respiration
US9572519B2 (en) Method and apparatus for invasive device tracking using organ timing signal generated from MPS sensors
EP3413076B1 (en) Mri involving the generation of a physiological motion signal using a pilot tone navigator
US6683454B2 (en) Shifting of artifacts by reordering of k-space
US7389136B2 (en) Method and system using a non-electrical sensor for gating
US11096630B2 (en) Generating a movement signal of a part of the human or animal body
JP2003525679A (ja) Ecgとppuの組合わせ制御による磁気共鳴方法及び装置
JP2006198407A (ja) 磁気共鳴(mr)イメージングにおける動き補償方法および動き補償システム
EP3105710A1 (en) Delay coordinate analysis of periodic data
JP6466912B2 (ja) 画像データから胸郭及び腹部呼吸信号を決定する装置及び方法
JP7318083B2 (ja) 磁気共鳴容積イメージング
JP2011147561A (ja) 磁気共鳴映像装置
CN109222928B (zh) 呼吸信号的提取方法、装置、处理设备和系统
JP3993799B2 (ja) 磁気共鳴イメージング装置
CN105266813B (zh) 一种自门控三维心脏成像的运动信号提取方法及装置
RU2346653C2 (ru) Способ и система для анализа сердечно-сосудистых звуков
CN103083030B (zh) D4d‑ct成像方法、设备和系统
US20100030572A1 (en) Temporal registration of medical data
US20160345926A1 (en) Acquisition of projection data for motion-corrected computed tomography images
EP4141465A1 (en) Assignment of mr images to cardiac phases
EP3554341A1 (en) Automated computation of trigger delay for triggered magnetic resonance imaging sequences
JP6510193B2 (ja) 磁気共鳴イメージング装置及び画像処理装置
Unay et al. Single syllable tongue motion analysis using tagged cine MRI

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20120808

Termination date: 20170206