CN104873204A - 磁共振装置及程序 - Google Patents

磁共振装置及程序 Download PDF

Info

Publication number
CN104873204A
CN104873204A CN201510089388.1A CN201510089388A CN104873204A CN 104873204 A CN104873204 A CN 104873204A CN 201510089388 A CN201510089388 A CN 201510089388A CN 104873204 A CN104873204 A CN 104873204A
Authority
CN
China
Prior art keywords
signal
magnetic resonance
passage
channel
coil
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.)
Granted
Application number
CN201510089388.1A
Other languages
English (en)
Other versions
CN104873204B (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.)
GE Medical Systems Global Technology Co LLC
Original Assignee
GE Medical Systems Global Technology Co LLC
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 GE Medical Systems Global Technology Co LLC filed Critical GE Medical Systems Global Technology Co LLC
Publication of CN104873204A publication Critical patent/CN104873204A/zh
Application granted granted Critical
Publication of CN104873204B publication Critical patent/CN104873204B/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/28Details of apparatus provided for in groups G01R33/44 - G01R33/64
    • G01R33/32Excitation or detection systems, e.g. using radio frequency signals
    • G01R33/36Electrical details, e.g. matching or coupling of the coil to the receiver
    • G01R33/3664Switching for purposes other than coil coupling or decoupling, e.g. switching between a phased array mode and a quadrature mode, switching between surface coil modes of different geometrical shapes, switching from a whole body reception coil to a local reception coil or switching for automatic coil selection in moving table MR or for changing the field-of-view
    • 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/28Details of apparatus provided for in groups G01R33/44 - G01R33/64
    • G01R33/32Excitation or detection systems, e.g. using radio frequency signals
    • G01R33/34Constructional details, e.g. resonators, specially adapted to MR
    • G01R33/341Constructional details, e.g. resonators, specially adapted to MR comprising surface coils
    • G01R33/3415Constructional details, e.g. resonators, specially adapted to MR comprising surface coils comprising arrays of sub-coils, i.e. phased-array coils with flexible receiver channels
    • 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/483NMR imaging systems with selection of signals or spectra from particular regions of the volume, e.g. in vivo spectroscopy
    • G01R33/4833NMR imaging systems with selection of signals or spectra from particular regions of the volume, e.g. in vivo spectroscopy using spatially selective excitation of the volume of interest, e.g. selecting non-orthogonal or inclined slices
    • G01R33/4835NMR imaging systems with selection of signals or spectra from particular regions of the volume, e.g. in vivo spectroscopy using spatially selective excitation of the volume of interest, e.g. selecting non-orthogonal or inclined slices of multiple slices

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Condensed Matter Physics & Semiconductors (AREA)
  • Engineering & Computer Science (AREA)
  • Health & Medical Sciences (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Radiology & Medical Imaging (AREA)
  • Signal Processing (AREA)
  • General Health & Medical Sciences (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Physiology (AREA)
  • Pulmonology (AREA)
  • Cardiology (AREA)
  • Power Engineering (AREA)
  • Biophysics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Magnetic Resonance Imaging Apparatus (AREA)
  • Vascular Medicine (AREA)
  • Artificial Intelligence (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Electromagnetism (AREA)

Abstract

本发明题为磁共振装置及程序。为了尽可能正确地检测被检体的活动,本发明提供执行用于从包含肝脏的摄影部位产生磁共振信号的扫描的磁共振装置(100),包括:具有多个通道(CH1~CH8)的线圈(4);从通道(CH1~CH8)中选择配置在肝脏的端部(E1)附近的通道(CH1、mCH2、CH5及CH6)的通道选择单元(82);以及基于通道(CH1、mCH2、CH5及CH6)接收的磁共振信号生成呼吸信号的呼吸信号生成单元(83)。

Description

磁共振装置及程序
技术领域
本发明涉及从包含活动部位的摄影部位收集磁共振信号的磁共振装置以及适用于该磁共振装置的程序。
背景技术
作为进行体动校正的方法,已知DC自导航法(参照非专利文献1)。
现有技术文献
非专利文献
非专利文献1:Brau et al.,Magnetic Resonance in Medicine 55:263-270(2006)。
发明内容
发明要解决的课题
在DC自导航法中,收集表示k空间中心的数据的DC信号,通过利用该DC信号,能够校正体动。另外,在DC自导航法中,利用与用于收集成像信号的RF脉冲相同的RF脉冲,能够收集DC信号。因此,无需考虑各自以RF脉冲分别收集成像信号和导航信号的情况下产生的自旋饱和效应,适合使用翻转角(flip angle)大的RF脉冲(例如,90°脉冲)的2D摄影。
另外,被检体的磁共振信号的接收,一般使用具有多个通道的线圈。特别是,近年来,出于适合宽范围部位的摄影等的理由,使用多通道的线圈。
然而,在DC自导航法的情况下,根据摄影部位和通道的位置关系,在线圈所具有的多个通道之中,会包含不适合检测被检体的活动的通道。因此,有难以检测被检体的活动,无法减少活动伪影的情况。基于这样的理由,希望有例如在利用DC自导航法对被检体进行摄影的情况下,尽可能地正确检测被检体的活动的方法。
用于解决课题的方案
本发明的第1观点是执行用于从包含活动部位的摄影部位产生第1磁共振信号的扫描的磁共振装置,其中包括:
线圈,具有接收所述第1磁共振信号的多个通道;
通道选择单元,从所述多个通道之中,选择配置在所述活动部位的端部附近的第1通道;以及
生成单元,基于所述第1通道接收的所述第1磁共振信号,生成包含所述扫描中的表示所述摄影部位的活动的活动信息的生物信号。
本发明的第2观点是一种程序,适用于具备执行用于从包含活动部位的摄影部位产生第1磁共振信号的扫描的扫描部、和具有接收所述第1磁共振信号的多个通道的线圈的磁共振装置,使计算机执行:
通道选择处理,从所述多个通道之中选择配置在所述活动部位的端部附近的第1通道;以及
生成处理,基于所述第1通道接收的所述第1磁共振信号,生成包含所述扫描中的表示所述摄影部位的活动的活动信息的生物信号。
发明效果
由于能够从多个通道之中选择配置在活动部位的端部附近的通道,所以能够取得更加正确的活动信息。
附图说明
图1是本发明的第1方式的磁共振装置的概略图;
图2是线圈4的说明图;
图3是概略示出线圈部AC的通道CH1~CH4与摄影部位的位置关系的图;
图4是概略示出线圈部PC的通道CH5~CH8与摄影部位的位置关系的图;
图5是示出处理器8执行的处理的图;
图6是在第1方式中执行的扫描的说明图;
图7是概略示出通过定位器扫描LS而取得的图像D的一个例子的图。
图8是概略示出由操作员设定的n块的切片L1~Ln的图;
图9是正式扫描MS的说明图;
图10是示出执行定位器扫描LS及正式扫描MS时的MR装置的动作流程的图;
图11是在期间P1中执行序列C1~Cn时的说明图;
图12是在期间P2中执行序列C1~Cn时的说明图;
图13是在期间Pm中执行序列C1~Cn时的说明图;
图14是示出在期间P1中执行序列C1时的情况的图;
图15是示出在期间P1中执行序列C2时的情况的图;
图16是示出在期间P1中执行序列Cn时的情况的图;
图17是示出在期间P1中合成从通道CH1~CH8输出的信号时的情况的图;
图18是示出合成信号A1的积分值S1的图;
图19是在期间P2算出呼吸信号时的说明图;
图20是概略示出各期间中的呼吸信号的图;
图21是示出呼吸信号的图;
图22是概略示出通道CH1及CH3与肝脏的位置关系的图;
图23是概略示出登记在数据库的数据的图;
图24是在期间P1中算出呼吸信号时的说明图;
图25是在期间P2中算出呼吸信号时的说明图;
图26是概略示出以第1方式的方法得到的呼吸信号的图;
图27是判断成像信号的接受/拒绝的方法的说明图;
图28是再收集成像信号B11~B1n时的说明图;
图29是在期间Pm+2中执行序列C1~Cn时的说明图;
图30是在第2方式中处理器执行的处理的说明图;
图31是在第2方式执行的扫描的说明图;
图32是预扫描PS的说明图;
图33是示出第2方式中的MR装置的动作流程的图;
图34是示出通道CH1~CH8接收DC信号A0时的情况的图;
图35是示出将通道CH1~CH8的输出信号A01~A08分别沿z方向进行傅里叶变换时的情况的图;
图36是概略示出曲线F1~F8的z方向的范围与切片Lc的z方向的范围的对应关系的图;
图37是示出中心位置zc的图;
图38是示出按每个曲线算出的积分值Sa及Sb的图;
图39是示出按每个曲线算出的积分值之比的图;
图40是示出第3方式中存储在存储器的数据库的图;
图41是第3方式中处理器执行的处理的说明图;
图42是示出第3方式中的MR装置的动作流程的图;
图43是第3方式中的正式扫描MS的说明图;
图44是第4方式中处理器执行的处理的说明图;
图45是示出第4方式中的MR装置的动作流程的图;
图46是预扫描PS的说明图;
图47是示出将通道CH1~CH4的输出信号A01~A04分别沿z方向进行傅里叶变换时的情况的图;
图48是示出曲线F1~F4的比J1~J4的图。
具体实施方式
以下,说明用于实施发明的方式,但是本发明并不限定于以下的方式。··
(1)第1方式
图1是本发明的第1方式的磁共振装置的概略图。
磁共振装置(以下,称为“MR装置”。MR:Magnetic Resonance)100具有磁体2、工作台3、接收RF线圈(以下,仅称为“线圈”)4等。
磁体2具有容纳被检体12的膛21。另外在磁体2内置有超导线圈、梯度线圈、RF线圈等(未图示)。超导线圈施加静磁场,梯度线圈施加梯度磁场,RF线圈施加RF脉冲。
工作台3具有托架3a。托架3a构成为能够在膛21内移动。利用托架3a可将被检体12输送至膛21。
线圈4安装在被检体12的躯干部。
图2是线圈4的说明图。
线圈4具有线圈部4a和线圈部4b。线圈部4a是配置在被检体12的前侧(腹部侧)的线圈,具有4个通道CH1、CH2、CH3及CH4。4个通道CH1~CH4以2行2列并排。
线圈部4b是配置在被检体12的后侧(脊背侧)的线圈,具有4个通道CH5、CH6、CH7及CH8。4个通道CH5~CH8以2行2列并排。
第1方式中,摄影对象的器官为肝脏,因此线圈部4a及4b安装在肝脏附近。
图3是概略示出线圈部4a的通道CH1~CH4与摄影部位的位置关系的图。图3(a)示出zx面内的通道的位置,图3(b)示出图3(a)的d-d截面中的通道的位置。
通道CH1及CH2沿x方向并排,通道CH3及CH4也沿x方向并排。与通道CH1相比,通道CH3在x方向的位置相同,但z方向的位置不同。另外,与通道CH2相比,通道CH4在x方向的位置相同,但z方向的位置不同。通道CH1及CH2配置在肝脏的端部E1附近,但是通道CH3及CH4配置在从肝脏的端部E1沿-z方向远离的位置。例如通道CH3配置在肝脏的与肺侧相反侧的端部E2附近。
图4是概略示出线圈部4b的通道CH5~CH8与摄影部位的位置关系的图。图4(a)示出zx面内的通道的位置,图4(b)示出在图4(a)的d-d截面中的通道的位置。
通道CH5及CH6沿x方向并排,通道CH7及CH8也沿x方向并排。与通道CH5相比,通道CH7在x方向的位置相同,但是z方向的位置不同。另外,与通道CH6相比,通道CH8在x方向的位置相同,但是z方向的位置不同。通道CH5及CH6配置在肝脏的端部E1附近,通道CH7及CH8配置在从肝脏的端部E1沿-z方向远离的位置。
根据图1继续进行说明。
MR装置100还包括发送器5、梯度磁场电源6、计算机7、操作部10及显示部11等。
发送器5向RF线圈供给电流,梯度磁场电源6向梯度线圈供给电流。此外,将磁体2、发送器5、梯度磁场电源6合在一起,相当于扫描部。
计算机7控制MR装置100的各部的动作,以实现向显示部11传输需要的信息或重构图像等、MR装置100的各种的动作。计算机7具有处理器8及存储器9等。
在存储器9中存储有通过处理器8执行的程序、后述的数据库(参照图23)等。处理器8读出存储在存储器9中的程序,执行描述在程序中的处理。图5示出处理器8执行的处理。处理器8通过读出存储在存储器9中的程序,构成切片设定单元81~判断单元84等。
切片设定单元81基于从操作部10输入的信息,设定切片。
通道选择单元82基于后述的数据库,从线圈4所具有的通道CH1~CH8之中选择配置在肝脏的端部E1(参照图3)附近的通道。
呼吸信号生成单元83基于经通道选择单元82选择的通道的接收信号,生成呼吸信号。
判断单元84判断是否将成像信号作为图像重构的信号接受。
处理器8是构成切片设定单元81~判断单元84的一个例子,通过执行既定程序,作为这些单元发挥功能。
操作部10由操作员操作,向计算机7输入各种信息。显示部11显示各种信息。
MR装置100如上所述地构成。
图6是第1方式中执行的扫描的说明图。
第1方式中,执行定位器扫描LS及正式扫描MS等。
定位器扫描LS是用于取得在设定切片时使用的图像D的扫描。定位器扫描LS中能取得轴向图像、矢状图像、冠状图像。图7中作为通过定位器扫描LS取得的图像D,仅示出冠状图像。
操作员基于图像D设定切片。图8中概略示出由操作员设定的n块的切片L1~Ln。在设定切片L1~Ln后,执行正式扫描MS。
图9是正式扫描MS的说明图。
正式扫描MS是用于利用多切片法取得n块的切片L1~Ln的图像的扫描。在正式扫描MS中,首先,在期间P1执行用于取得切片L1~Ln的图像的序列C1~Cn。图9中概略示出序列C1的例子。序列C1构成为利用DC自导航法能够收集表示k空间的中心的数据的MR信号(以下,称为“DC信号”)A和用于作成图像的MR信号(以下,称为“成像信号”)B。
序列C1具有用于激励切片L1的RF脉冲α。从通过RF脉冲α来激励的切片L1收集成像信号B。另外,RF脉冲α不仅用于收集成像信号B,而且也用于收集DC信号A。DC信号A是在设定在即将施加梯度磁场Gy及Gz之前的待机时间Twait期间收集的。待机时间Twait例如为20μs。
在执行序列C1之后,依次执行用于取得切片L2~Ln的图像的序列C2~Cn。除了RF脉冲α的激励频率之外,以与序列C1相同的序列图表示序列C2~Cn。因此,在期间P1,每次执行序列C1~Cn的各个序列时,收集DC信号A及成像信号B。
期间P1中执行序列C1~Cn后,在下个期间P2也执行序列C1~Cn。以下同样如此,反复执行序列C1~Cn。图9中示出期间P1~Pm中执行的序列C1~Cn。此外,序列Cl~Cn的相位编码量按每个期间变化。
以下,对执行定位器扫描LS及正式扫描MS时的MR装置的动作流程进行详细说明。
图10是示出执行定位器扫描LS及正式扫描MS时的MR装置的动作流程的图。
步骤ST1中,执行定位器扫描LS。通过执行定位器扫描LS,取得图像D(参照图7)。在执行定位器扫描LS后,进入步骤ST2。
步骤ST2中,操作员操作操作部10(参照图1),一边参考图像D一边输入用于设定切片L1~Ln(参照图8)的信息。切片设定单元81(参照图5)基于从操作部10输入的信息,设定切片L1~Ln。在设定切片L1~Ln后,进入步骤ST3。
步骤ST3中,执行正式扫描MS。正式扫描MS中,首先,期间P1中,执行序列C1~Cn(参照图11)。
图11是期间P1中执行序列C1~Cn时的说明图。
图11概略地示出通过在期间P1中执行序列C1~Cn来收集的DC信号A和成像信号B。此外,在图11中,为了区分期间P1中得到的多个DC信号A,对标号A附加有添加字“11”、“12”、...“1n”。同样地,为了区分成像信号B,对标号B也附加有添加字“11”、“12”、...“1n”。
期间P1中,首先,执行序列C1。通过执行序列C1,收集DC信号A11和成像信号B11。成像信号B11作为切片L1的ky=32的直线的数据使用。在执行序列C1后,执行序列C2。
通过执行序列C2,收集DC信号A12和成像信号B12。成像信号B12作为切片L2的ky=32的直线的数据使用。
以下同样如此,依次执行用于从切片L3~Ln的每一个收集DC信号及成像信号的序列。在期间P1的最后,执行用于收集切片Ln的数据的序列Cn。通过执行序列Cn,收集DC信号A1n和成像信号B1n。成像信号B1n作为切片Ln的ky=32的直线的数据使用。
因此,在期间P1的期间,能够收集切片L1~Ln的ky=32的数据。接着,转移到期间P2
图12是在期间P2中执行序列C1~Cn时的说明图。
在图12中,为了区分在期间P2得到的多个DC信号A,对标号A附加添加字“21”、“22”、...“2n”。同样地,为了区分成像信号B,对标号B也附加添加字“21”、“22”、...“2n”。
期间P2中,首先执行序列Cl。通过执行序列Cl,收集DC信号A21和成像信号B21。成像信号B21表示切片L1的ky=31的直线的数据。在执行序列C1后,依次执行序列C2~Cn。因此,在期间P2的期间,能够收集切片L1~Ln的ky=31的数据。
在期间P2收集ky=31的数据后,也反复执行用于收集剩下的ky视图的数据的序列C1~Cn(参照图13)。
图13是在期间Pm中执行序列C1~Cn时的说明图。
在图13中,为了区分在期间Pm得到的多个DC信号A,对标号A附加添加字“m1”、“m2”、...“mn”。同样地,为了区分成像信号B,对标号B也附加添加字“m1”、“m2”、...“mn”。
期间Pm中,首先,执行序列C1。通过执行序列C1,收集DC信号Am1和成像信号Bm1。成像信号Bm1表示切片L1的ky=-32的直线的数据。在执行序列C1后,依次执行序列C2~Cn。因此,在期间Pm的期间,能够收集切片L1~Ln的ky=-32的数据。
另外,通过执行序列C1~Cn,除了成像信号B之外,还收集DC信号A。在第1方式中,利用DC信号A,生成被检体的呼吸信号。以下,对第1方式中生成呼吸信号的方法进行说明。此外,在说明第1方式的呼吸信号的生成方法时,为了明确第1方式的呼吸信号的生成方法的效果,首先,对以不同于第1方式的方式生成呼吸信号的方法的一个例子进行说明,随后,对第1方式的呼吸信号的生成方法进行说明。
图14~图19是以不同于第1方式的方式生成呼吸信号的方法的一个例子的说明图。
首先,如图14所示,在期间P1中,执行序列C1。通过执行序列C1,从切片L1收集到DC信号A11和成像信号B11
由于线圈4具有通道CH1~CH8,所以能由通道CH1~CH8分别接收DC信号A11。在图14的下侧,当通道CH1~CH8接收DC信号A11的情况下,以标号“A11,1”~“A11,8”表示从通道CH1~CH8分别输出的信号。
在执行序列C1后,执行序列C2。图15示出执行序列C2时的情况。通过执行序列C2,从切片L2收集到DC信号A12和成像信号B12
由于线圈4具有通道CH1~CH8,所以DC信号A12与DC信号A11同样,由通道CH1~CH8分别接收。在图15的下侧,当通道CH1~CH8接收DC信号A12的情况下,以标号“A12,1”~“A12,8”表示从通道CH1~CH8分别输出的信号。
以下同样如此,执行用于从切片L3~Ln分别收集DC信号及成像信号的序列。而且,在期间P1的最后,执行用于收集切片Ln的数据的序列Cn。图16示出执行序列Cn时的情况。通过执行序列Cn,从切片Ln收集到DC信号A1n和成像信号B1n
由于线圈4具有通道CH1~CH8,所以DC信号A1n与DC信号A11同样,由通道CH1~CH8分别接收。在图16的下侧,当通道CH1~CH8接收DC信号A1n的情况下,以标号“A1n,1”~“A1n,8”表示从通道CH1~CH8分别输出的信号。
这样,在每次执行序列时,从各通道输出DC信号。
接着,在期间P1中合成从通道CH1~CH8输出的信号(参照图17)。
图17是示出在期间P1中合成从通道CH1~CH8输出的信号时的情况的图。
在图17中,示出在期间P1中通过相加从通道CH1~CH8输出的信号,合成通道CH1~CH8的信号的例子。通过相加通道CH1~CH8的全部信号,能得到合成信号A1
在取得合成信号A1后,算出合成信号A1的积分值。在图18中以标号“S1”表示算出的合成信号A1的积分值。积分值S1作为期间P1中的被检体的呼吸信号的信号值使用。
在期间P1中执行序列后,转移到期间P2
图19是算出期间P2中的呼吸信号的信号值时的说明图。
在期间P2,也与期间P1同样执行地序列,合成通道的信号。而且,求出合成信号A2的积分值S2。积分值S2作为期间P2中的被检体的呼吸信号的信号值使用。
以下同样如此,在各期间执行序列C1~Cn,算出合成信号的积分值。因此,能够求出各期间中的呼吸信号的信号值(参照图20)。
图20概略示出各期间中的呼吸信号的信号值的图。
图20中示出以图14~图19所示的方法求出的呼吸信号Q1和理想的呼吸信号Q2。
为了识别被检体的呼吸状态(呼气、吸气等),如理想的呼吸信号Q2所示,需要呼吸信号根据被检体的呼吸运动与时间一起尽可能大地变化。然而,以图14~图19所示的方法生成呼吸信号的情况下,由于呼吸信号的振幅变小,所以存在难以得到合适的呼吸信号的问题。
因此,本申请发明人,为了使呼吸信号的振幅变小的理由清楚,利用图9所示的序列实际扫描被检体,对从全部通道的合成信号求出的呼吸信号和从1个通道的接收信号求出的呼吸信号之差异进行了研究。以下,对研究结果进行说明。
图21是示出呼吸信号的图。
图21(a)是示出从全部通道的合成信号求出的呼吸信号V0的图。参照图21(a),可知呼吸信号V0的振幅不怎么变大。
图21(b1)~(b8)是示出仅从1个通道的接收信号求出的呼吸信号的图。以下,对各图进行说明。
图21(b1)是示出仅从通道CH1的接收信号求出的呼吸信号VI的图。可知通道CH1的呼吸信号VI(周期T)反映肝脏的活动而较大地变化。
图21(b2)是示出仅从通道CH2的接收信号求出的呼吸信号V2的图。通道CH2的呼吸信号V2也与通道CH1的呼吸信号V1同样,反映肝脏的活动而较大地变化。
图21(b3)是示出仅从通道CH3的接收信号求出的呼吸信号V3的图。通道CH3的呼吸信号V3与通道CH1的呼吸信号V1同样反映肝脏的活动而较大地变化。然而,通道CH3的呼吸信号V3相对于通道CH1的呼吸信号V1,波形在时间方向仅偏移△T。
图21(b4)是示出仅从通道CH4的接收信号求出的呼吸信号V4的图。通道CH4的呼吸信号V4与通道CH1的呼吸信号V1同样反映肝脏的活动而较大地变化。然而,通道CH4的呼吸信号V4相对于通道CH1的呼吸信号V1,波形在时间方向仅偏移△T。
图21(b5)是示出仅从通道CH5的接收信号求出的呼吸信号V5的图。通道CH5的呼吸信号V5与通道CH1的呼吸信号V1同样,反映肝脏的活动而较大地变化。通道CH5的呼吸信号V5相对于通道CH1的呼吸信号V1,几乎看不到波形在时间方向的偏移。
图21(b6)是示出仅从通道CH6的接收信号求出的呼吸信号V6的图。通道CH6的呼吸信号V6与通道CH1的呼吸信号V1同样,反映肝脏的活动而较大地变化。通道CH6的呼吸信号V6相对于通道CH1的呼吸信号V1,几乎看不到波形在时间方向的偏移。
图21(b7)是示出仅从通道CH7的接收信号求出的呼吸信号V7的图。可知通道CH7的呼吸信号V7的振幅不怎么变化,不充分地反映肝脏的活动。
图21(b8)是示出仅从通道CH8的接收信号求出的呼吸信号V8的图。可知通道CH8的呼吸信号V8的振幅不怎么变化,不充分地反映肝脏的活动。
因此,可知相对于通道CH1的呼吸信号V1,包含波形在时间方向仅偏移△T的呼吸信号。例如,通道CH3的呼吸信号V3相对于通道CH1的呼吸信号V1,波形在时间方向仅偏移△T。以下,考察呼吸信号的波形在时间方向偏移的理由。
图22是概略示出通道CH1及CH3与肝脏的位置关系的图。在图22中,以实线示出被检体吐气时的肝脏,以虚线示出被检体吸气时的肝脏的位置。
在被检体吐气的情况下,肝脏的端部E1沿z方向活动,因此肝脏接近通道CH1。因此,通道CH1的接收信号的信号值受到肝脏的影响而增加。另一方面,从通道CH3来看,肝脏离去,因此通道CH3的接收信号的信号值减少。
相对于此,在被检体吸气的情况下,肝脏的端部E1沿-z方向活动,肝脏从通道CH1离开。因此,通道CH1的接收信号的信号值减少。另一方面,从通道CH3来看,肝脏靠近来,因此通道CH3的接收信号的信号值增加。因此,通道CH1的接收信号的增减和通道CH3的接收信号的增减相反。
这样接收信号的增减相反,因此认为从通道CH3的接收信号求出的呼吸信号V3相对于从通道CH1的接收信号求出的呼吸信号V1,波形在时间方向仅偏移△T。因此,若相加呼吸信号V1和V3,信号就互相抵消。
另外,参照图21(bl)~(b8),可知通道CH7及CH8的呼吸信号V7及V8的振幅不怎么变化。作为该理由,认为是由于通道CH7及CH8比其他通道更远离肝脏,所以肝脏的活动造成的信号值的变化不怎么大。
因而,在通道CH1~CH8之中,包含将彼此的信号抵消的通道、不怎么反映肝脏的活动的通道,因此认为在合成全部通道的接收信号的情况下,呼吸信号的振幅不怎么变大。
此外,参照图21(bl)~(b8),可知在通道CH1~CH8之中,包含几乎看不到呼吸信号的波形在时间方向的偏移的通道。具体而言,通道CH1、CH2、CH5及CH6的呼吸信号V1、V2、V5及V6,几乎看不到时间方向的偏移。这些通道CH1、CH2、CH5及CH6,配置在肝脏的端部E1附近,对应于肝脏的活动的信号值的增减以相同定时进行,因此认为几乎看不到呼吸信号的波形在时间方向的偏移。因此,认为如果仅合成通道CH1、CH2、CH5及CH6的接收信号,就能得到跟随被检体的呼吸较大地变化的呼吸信号。因此,第1方式中,仅利用通道CH1~CH8中配置在肝脏的端部E1附近的通道CH1、CH2、CH5及CH6的接收信号,生成呼吸信号。以下,对第1方式中的呼吸信号的生成方法进行说明。
第1方式中,对存储器9(参照图1)存储有包含线圈的通道的信息的数据库。图23中概略示出登记在数据库的数据。
数据库中登记有表示线圈的项目a、表示线圈的通道的项目b、和表示通道是否配置在肝脏的肺侧的端部E1附近的项目c。项目c的记号“○”表示通道配置在肝脏的端部E1附近。在此,通道CH1、CH2、CH5及CH6作为配置在肝脏的端部E1附近的通道而登记。
第1方式中,基于图23的数据库生成呼吸信号。以下,参照图24及图25,对利用数据库生成呼吸信号的顺序进行说明。
首先,如图24所示,期间P1中,执行序列C1。通过执行序列C1,从切片L1收集到DC信号A11和成像信号B11
由于线圈4具有通道CH1~CH8,由通道CH1~CH8分别接收DC信号A11。通道CH1~CH8接收DC信号A11,从而分别输出信号A11,1~A11,8
在执行序列C1后,执行序列C2。通过执行序列C2,从切片L2收集到DC信号A12和成像信号B12
DC信号A12与DC信号A11同样,由通道CH1~CH8分别接收。通道CH1~CH8接收DC信号A12,从而分别输出信号A12,1~A12,8
以下同样如此,执行用于从切片L3~Ln的各个切片收集DC信号及成像信号的序列。而且,在期间P1的最后,执行用于收集切片Ln的数据的序列Cn。通过执行序列Cn,从切片Ln收集到DC信号A1n和成像信号B1n
由通道CH1~CH8分别接收的DC信号A1n。通道CH1~CH8接收DC信号A1n,分别输出信号A1n,1~A1n,8
在期间P1中执行序列C1~Cn后,如以下那样生成呼吸信号。
首先,通道选择单元82(参照图5)参照数据库(参照图23)。而且,通道选择单元82基于数据库的项目c的信息,选择作为配置在肝脏的端部E1附近的通道而登记的通道CH1、CH2、CH5及CH6。
接着,呼吸信号生成单元83(参照图5)丢弃通道CH1~CH8之中未被选择的通道CH3、CH4、CH7及CH8的输出信号,仅合成(加法运算)选择的通道CH1、CH2、CH5及CH6的输出信号。由此,得到合成信号A1
在取得合成信号A1后,呼吸信号生成单元83算出合成信号A1的积分值S1。积分值S1用作期间P1中的被检体的呼吸信号的信号值。
在期间P1中执行序列后,转移到期间P2
图25是期间P2中的算出呼吸信号时的说明图。
在期间P2也与期间1同样执行序列C1~Cn。呼吸信号生成单元83丢弃通道CH3、CH4、CH7及CH8的输出信号,仅合成(加法运算)选择的通道CH1、CH2、CH5及CH6的输出信号,生成合成信号A2。而且,呼吸信号生成单元83求出合成信号A2的积分值S2。积分值S2用作期间P2中的被检体的呼吸信号的信号值。
以下同样如此,在各期间中执行序列C1~Cn。呼吸信号生成单元83丢弃通道CH3、CH4、CH7及CH8的输出信号,仅合成(加法运算)选择的通道CH、CH2、CH5及CH6的输出信号。而且,算出合成信号的积分值。因此,能够求出各期间中的呼吸信号(参照图26)。
图26是概略示出以第1方式的方法得到的呼吸信号的图。
第1方式中,仅合成(加法运算)位于肝脏的端部E1附近的通道CH1、CH2、CH5及CH6的输出信号。由于通道CH1、CH2、CH5及CH6的输出信号以相同定时增减,所以通过仅合成这些通道的输出信号,能够取得对应于被检体的呼吸而较大地变化的呼吸信号Vsyn。
此外,由于肝脏随着呼吸而活动,如果仅用期间P1~Pm收集的成像信号重构图像,就会在图像中显现体动人为因素。因此,第1方式中,为了减轻体动人为因素,基于呼吸信号Vsyn,判断将成像信号作为在图像重构中使用的信号接受,还是拒绝接受成像信号。以下,对其判断方法进行说明。
图27是判断成像信号的接受、拒绝的方法的说明图。
判断单元84(参照图5)首先求出相当于被检体的结束吐气的位置的信号值x0。结束吐气的信号值x0例如能够参考呼吸信号的峰值而求出。接着,求出呼吸信号的最大值与最小值之差△D。而且,以结束吐气的信号值x0为中心,设定差△D的x%(例如,x=20)的范围AW。将这样设定的范围AW定为接受成像信号B的容许范围AW。而且,判断单元84在呼吸信号包含于容许范围AW的情况下,判断为接受成像信号作为图像重构中使用的信号。另一方面,判断单元84在呼吸信号不包括在容许范围AW的情况下,判断为拒绝将成像信号接受为使用于图像重构的信号。观看图27,期间P1的信号值(积分值)S1未包括在容许范围AW,因此拒绝在期间P1收集的成像信号B11~B1n(参照图24)。然而,期间P2的信号值(积分值)S2包括在容许范围AW,因此判断为接受期间P2收集的成像信号B21~B2n(参照图25)。以下同样如此,根据各期间的呼吸信号是否包括在容许范围AW,判断是接受成像信号还是拒绝接受成像信号。
而且,在期间P1~Pm收集的成像信号之中,被拒绝作为使用于图像重构用途的信号而接受的成像信号,在期间Pm之后再收集。例如,在期间P1收集的成像信号B11~B1n(参照图24)被拒绝接受为使用于图像重构用途的信号,因此再收集成像信号B11~B1n(参照图28)。
图28是再收集成像信号B11~B1n时的说明图。
在期间Pm+1中,执行用于收集成像信号B11~B1n的序列C1~Cn。通过执行序列C1~Cn,再收集DC信号A11~A1n和成像信号B11~B1n
DC信号A11~A1n及成像信号B11~B1n,由通道CH1~CH8各自接收。此外,在图28中,说明的方便起见,仅示出DC信号A11~A1n被通道CH1~CH8各自接收的情况。通道CH1~CH4分别输出信号A11,1~A11,8
呼吸信号生成单元83生成通道CH1、CH2、CH5及CH6的输出信号的合成信号,算出合成信号Am+1的积分值Sm+1。由此,能够得到期间Pm+1中的呼吸信号Sm+1
接着,判断单元84判断呼吸信号Sm+1是否包括在容许范围AW。在图28中,呼吸信号Sm+1不包括在容许范围AW,因此在期间Pm+1收集的成像信号B11~B1n不能接受为图像重构用的数据。因此,拒绝成像信号B11~B1n。在该情况下,接着,下个期间Pm+2中,执行用于再收集成像信号B11~B1n的序列C1~Cn(参照图29)。
图29是在期间Pm+2中执行序列C1~Cn时的说明图。
在期间Pm+2,也与期间Pm+1同样,执行用于再收集成像信号B11~B1n的序列C1~Cn。通过执行序列C1~Cn,再收集DC信号A11~A1n和成像信号B11~B1n
DC信号A11~A1n及成像信号B11~B1n,由通道CH1~CH8各自接收。此外,在图29中,说明的方便起见,仅示出通道CH1~CH8各自接收DC信号A11~A1n的情况。通道CH1~CH4分别输出信号A11,1~A11,8
呼吸信号生成单元83生成通道CH1、CH2、CH5及CH6的输出信号的合成信号,算出合成信号Am+2的积分值Sm+2。由此,能够得到期间Pm+2中的呼吸信号Sm+2。
接着,判断单元84判断呼吸信号Sm+2是否包括在容许范围AW。在图29中,由于呼吸信号Sm+2不包括在容许范围AW,所以判断为期间Pm+2中收集的成像信号B11~B1n被接受为图像重构用的数据。
以下同样如此,在再收集其他拒绝的成像信号的情况下,也反复执行序列,直至呼吸信号包括在容许范围AW。因此,能够将呼吸信号包括在容许范围AW时收集的ky=-32~32的成像信号作为图像重构用的数据取得。
这样,再收集拒绝的成像信号后,进行图像重构。
第1方式中,合成配置在肝脏的端部E1附近的通道接收的DC信号,因此能够取得根据被检体的呼吸而较大地变化的呼吸信号Vsyn。因此,能够确定被检体吐完气时的呼吸信号的大体范围AW。另外,第1方式中,在呼吸信号不包括在范围AW的情况下,再收集成像信号,直至呼吸信号包括在范围AW为止,因此能够得到体动人为因素减少的图像。
此外,第1方式中,作为配置在肝脏的端部E1附近的通道,登记有4个通道CH1、CH2、CH5及CH6。然而,也可以不全部登记4个通道CH1、CH2、CH5及CH6,而仅登记4个通道CH1、CH2、CH5及CH6中的1个通道、2个通道或3个通道。如参照图21进行说明的那样,在通道CH1、CH2、CH5及CH6的哪个通道,都能得到充分地反映肝脏的活动的呼吸信号。因此,如果登记有4个通道CH1、CH2、CH5及CH6中的至少一个以上的通道,就能得到充分地反映肝脏的活动的呼吸信号。
另外,也可以取代位于肝脏的端部E1的通道CH1、CH2、CH5及CH6,而预先登记例如位于肝脏的端部E2(参照图3)附近的通道CH3。如参照图21进行说明的那样,从通道CH3得到的呼吸信号,相对于从通道CH1得到的呼吸信号,波形在时间方向仅偏移△T,但充分地反映肝脏的活动。因此,取代通道CH1、CH2、CH5及CH6而登记通道CH3,也能得到充分地反映肝脏的活动的呼吸信号。
第1方式中,基于操作员从操作部10输入的信息,切片设定单元81设定切片。然而,基于图像D,切片设定单元81也可以自动设定切片。
(2)第2方式
第1方式中,预先在数据库登记配置在肝脏的端部E1附近的通道CH1、CH2、CH5及CH6,通过参照数据库的信息,从通道CH1~CH8之中选择通道CH1、CH2、CH5及CH6。第2方式中,对数据库中不登记通道,而从通道CH1~CH8之中选择配置在肝脏的端部E1附近的通道CH1、CH2、CH5及CH6的例子进行说明。此外,MR装置的硬件构成与第1方式相同。
图30是第2方式中执行处理器的处理的说明图。
处理器8通过读出存储在存储器9中的程序,构成切片设定单元81~判断单元84等。
切片设定单元81基于从操作部10输入的信息设定切片。
曲线作成单元811基于通过后述的预扫描PS(参照图32)收集的MR信号,作成表示摄影部位的z方向上的各位置和信号值的关系的曲线。
通道选择单元82基于由曲线作成单元811作成的曲线,从线圈4所具有的通道CH1~CH8之中,选择配置在肝脏的端部E1(参照图3)附近的通道。
呼吸信号生成单元83基于由通道选择单元82选择的通道的输出信号,生成呼吸信号。
判断单元84判断是否将成像信号接受为图像重构的信号。
处理器8是构成切片设定单元81~判断单元84的一个例子,通过执行既定程序,作为这些单元发挥功能。
图31是第2方式中执行的扫描的说明图。
第2方式中,执行定位器扫描LS、预扫描PS及正式扫描MS。第2方式与第1方式相比,执行定位器扫描LS及正式扫描MS这一点相同,但在定位器扫描LS与正式扫描MS之间执行预扫描PS这一点不同。
图32是预扫描PS的说明图。
图32中示出在预扫描PS中执行的序列H。序列H除了不沿Gy方向施加相位编码梯度脉冲这一点外,与图9所示的脉冲序列相同。
预扫描PS是为了从通道CH1~CH8之中选择配置在肝脏的端部E1附近的通道CH1、CH2、CH5及CH6而执行的扫描。对于预扫描PS,将在后面进行详细说明。
以下,对第2方式中执行定位器扫描LS、预扫描PS及正式扫描MS时的MR装置的动作流程进行说明。
图33是示出第2方式中的MR装置的动作流程的图。
步骤ST1及ST2与第1方式相同,因此省略详细的说明。步骤ST2中设定切片L1~Ln(参照图8)后,进入步骤ST21。
步骤ST21中执行预扫描PS。预扫描PS是为了从通道CH1~CH8之中选择配置在肝脏的端部E1附近的通道CH1、CH2、CH5及CH6而执行的扫描。以下,对预扫描PS进行说明(参照图34)。
图34是预扫描PS的说明图。
预扫描PS中,只有切片L1~Ln中的任意一块的切片被激励,从被激励的一块切片收集DC信号A0及成像信号B0。第2方式中,激励切片L1~Ln中的中央的切片Lc。因此,从切片Lc收集到DC信号A0和成像信号B0
通过执行序列H,从切片Lc收集到DC信号A0及成像信号B0。DC信号A0及成像信号B0分别由通道CH1~CH8接收。在图34中,说明的方便起见,仅示出DC信号A0被通道CH1~CH8分别接收的情况。此外,在DC信号A0及成像信号B0之中,通道的选择所使用的信号为成像信号B0,而DC信号A0不用在通道的选择。
由于通道CH1~CH8接收成像信号B0,所以通道CH1~CH8分别输出信号B01~B08
在执行预扫描PS后,进入步骤ST22。
在步骤ST22中,曲线作成单元811(参照图30)将通道CH1~CH8的输出信号B01~B08分别沿z方向进行傅里叶变换(FT:Fourier Transformation)。由此,如图35所示,能够按每个通道,作成表示z方向的各位置和信号值的关系的曲线(F1~F8)。图36中概略地示出曲线F1~F8的z方向的范围。在图36的左侧示出曲线F1~F4的z方向的范围,在图36的右侧示出曲线F5~F8的z方向的范围。
曲线F1~F8的z方向的范围以标号“za”及“zb”来表示。za是肝脏的端部E2附近的位置,zb是横切肺的位置。
在作成曲线F1~F8后,进入步骤ST23。
在步骤ST23中,通道选择单元82(参照图30)求出表示曲线F1~F8各自的特性的特性值,基于特性值,从通道CH1~CH8之中选择配置在肝脏的端部E1附近的通道。以下,首先对求出曲线CH1~CH8各自的特性值的方法进行说明,接着,对基于特性值选择通道的方法进行说明。
图37~图39是求出曲线CH1~CH8各自的特性值的方法的说明图。
通道选择单元82首先规定将曲线F1~F8各自的z方向的范围za~zb二分的中心位置zc。在图37示出中心位置zc。在规定中心位置zc后,通道选择单元82针对各曲线算出区间za~zc中的积分值Sa和区间zc~zb中的积分值Sb。图38中示出按每个曲线算出的积分值Sa及Sb。
在算出积分值Sa及Sb后,通道选择单元82按每个曲线算出积分值Sb与Sa的比。图39中示出按每个曲线算出的积分值的比。在图39中,曲线F1~F8的比由标号“J1”~“J8”表示。第2方式中,积分值的比作为曲线的特性值而求出。
将比J1~J8进行比较时,比J1~J8能够根据通道的配置位置分为具有大的值的比和具有小的值的比。以下,对该理由进行说明。
首先,考虑比J1~J8中的4个比J1~J4(参照图39的左侧)。
通道CH1及CH2相对于中心位置zc沿z方向配置,但是通道CH3及CH4相对于中心位置zc配置在-z方向。因此,在范围zc~zb中,通道CH1及CH2具有比通道CH3及CH4高的灵敏度。因此,通道CH1及CH2的曲线F1及F2的积分值Sb成为比通道CH3及CH4的曲线F3及F4的积分值Sb大的值。另一方面,在范围za~zc中,通道CH1及CH2具有比通道CH3及CH4低的灵敏度。因此,通道CH1及CH2的曲线F1及F2的积分值Sa成为比通道CH3及CH4的曲线F3及F4的积分值Sa小的值。
因此,可知通道CH1及CH2的比J1及J2变得比通道CH3及CH4的比J大。
上述说明中,虽然对通道CH1~CH4的比J1~J4进行了说明,但也能同样地说明通道CH5~CH8的比J5~J8。通道CH5及CH6的比J5及J6变得比通道CH7及CH8的比J7及J8大。
因此,可知通过从比J1~J8之中确定值大的比,能够选择配置在肝脏的端部E1附近的通道CH1、CH2、CH5及CH6。
因此,通道选择单元82将比J1~J8按照从大到小的值顺序排列,按值的从大到小的顺序确定4个比。在此,确定J1、J2、J5及J6作为值大的4个比。因此,能够从通道CH1~CH8之中选择配置在肝脏的端部E1附近的通道CH1、CH2、CH5及CH6。
在选择通道后,进入步骤ST3。
在步骤ST3中,执行正式扫描MS。正式扫描MS中与第1方式同样,仅合成通道CH1、CH2、CH5及CH6的输出信号,生成呼吸信号。
而且,与第1方式同样,基于期间P1~Pm的呼吸信号,设定接受成像信号B的容许范围AW(参照图27),在呼吸信号不包括在容许范围AW的情况下再收集数据,结束流程。
第2方式中执行预扫描PS,基于通过预扫描PS收集的MR信号,作为通道CH1~CH8的曲线F1~F8。而且,算出曲线F1~F8的比J1~J8。比J1~J8的值能够根据通道的配置位置分为大的值和小的值,因此基于比J1~J8,能够选择配置在肝脏的端部E1附近的通道。另外,即使是与线圈4不同的线圈,也通过执行预扫描PS,能够从该不同的线圈所具有的通道之中选择配置在肝脏的端部E1附近的通道。因此,无需在数据库按每个摄影所使用的线圈先登记通道,因此也能够减轻花费在数据库的维护的负担。
在第2方式中,作为曲线的特性值,算出曲线的积分值的比(J1~J8)。然而,只要能够区别通道CH1、CH2、CH5及CH6和通道CH3、CH4、CH7及CH8,也可以取代积分值的比而求出其他特性值。例如,算出范围za~zc的信号值的最大值和范围zc~zb的信号值的最大值,将最大值的比作为曲线的特性值也可。
此外,第2方式中,作为配置在肝脏的端部E1附近的通道,通道选择单元82选择4个通道CH1、CH2、CH5及CH6。然而,不选择4个通道CH1、CH2、CH5及CH6的全部,而仅选择4个通道CH1、CH2、CH5及CH6中的1个通道、2个通道、或3个通道也可。如参照图21进行说明的那样,在通道CH1、CH2、CH5及CH6的哪个通道都能得到充分地反映肝脏的活动的呼吸信号。因此,如果选择4个通道CH1、CH2、CH5及CH6中的至少一个以上的通道,就能得到充分地反映肝脏的活动的呼吸信号。
第2方式中,在预扫描PS中,从切片Lc收集磁共振信号,作为各通道的曲线。然而,从与切片Lc不同的切片收集磁共振信号,作成各通道的曲线也可。另外,从多个切片收集磁共振信号,作成各通道的曲线也可。进而,第2方式中,预扫描PS是2D扫描,但也可为3D扫描。
(3)第3方式
第3方式中,对线圈4具有多个线圈模式的情况进行说明。此外,线圈4以外的MR装置的硬件构成与第1方式(参照图1)相同。
第3方式中,线圈4根据摄影条件,构成为能够按以下的线圈模式接收MR信号。
(1)线圈模式Ml(通道CH1+CH2+CH3+CH4)
(2)线圈模式M2(通道CH5+CH6+CH7+CH8)
(3)线圈模式M3(通道CH1+CH2+CH3+CH4+CH5+CH6+CH7+CH8)。
线圈模式M1表示以4个通道CH1~CH4接收MR信号的模式。线圈模式M2表示以4个通道CH5~CH8接收MR信号的模式。线圈模式M3表示以8个通道CH1~CH8接收MR信号的模式。
图40是示出第3方式中存储在存储器9的数据库的图。
数据库中登记有:表示线圈4的项目a;表示线圈4所具有的线圈模式的项目a1;表示线圈模式所具有的通道的项目b;以及表示通道是否配置在肝脏的肺侧的端部E1附近的项目c。项目c的记号“○”表示通道配置在肝脏的端部E1附近。在此,作为配置在肝脏的端部E1附近的通道,登记有通道CH1、CH2、CH5及CH6。
图41是第3方式中处理器执行的处理的说明图。
处理器8通过读出存储在存储器9的程序,构成线圈模式选择单元80~判断单元84等。
线圈模式选择单元80基于从操作部10输入的信息,从线圈模式M1~M3中选择摄影所使用的线圈模式。
切片设定单元81基于从操作部10输入的信息,设定切片。
通道选择单元82基于数据库(参照图40),从所选择的线圈模式所具有的通道中,选择配置在肝脏的端部E1(参照图3)附近的通道。
呼吸信号生成单元83基于由通道选择单元82选择的通道的输出信号,生成呼吸信号。
判断单元84判断是否将成像信号接受为图像重构的信号。
处理器8是构成线圈模式选择单元80~判断单元84的一个例子,通过执行既定程序,作为这些单元起作用。
以下对第3方式中的MR装置的动作流程进行说明。
图42是示出第3方式中的MR装置的动作流程的图。
步骤ST0中,在执行定位器扫描LS之前,操作员对操作部10进行操作,输入用于从线圈模式M1~M3中选择摄影被检体时使用的线圈模式的信息。当该信息输入时,线圈模式选择单元80(参照图41)基于从操作部10输入的信息,从线圈模式M1~M3中选择摄影被检体时使用的线圈模式。在此,设为选择了线圈模式Ml。在选择线圈模式Ml之后,进入步骤ST1。
步骤ST1中,利用线圈模式Ml执行定位器扫描LS。通过执行定位器扫描LS,取得图像D(参照图7)。在执行定位器扫描LS之后,进入步骤ST2。
步骤ST2中,操作员基于图像D,设定切片L1~Ln(参照图8)。在设定切片L1~Ln之后,进入步骤ST3。
步骤ST3中,执行正式扫描MS。
图43是第3方式中的正式扫描MS的说明图。
首先,在期间P1中,执行序列C1。通过执行序列C1,从切片L1收集DC信号A11和成像信号B11
第3方式中,选择线圈模式Ml,因此DC信号A11及成像信号B11被线圈模式M1的通道CH1~CH4分别接收。此外,在图43中,说明的方便起见,仅示出DC信号A11被线圈模式Ml的通道CH1~CH4分别接收的情况。通道CH1~CH4分别输出信号A11,1~A11,4
在执行序列C1之后,执行序列C2。通过执行序列C2,从切片L2收集到DC信号A12和成像信号B12
DC信号A12及成像信号B12被线圈模式Ml的通道CH1~CH4分别接收。此外,在图43中,说明的方便起见,仅示出DC信号A12被线圈模式M1的通道CH1~CH4分别接收的情况。通道CH1~CH4分别输出信号A12,1~A12,4
以下同样如此,执行用于从切片L3~Ln各自收集DC信号及成像信号的序列。而且,在期间P1的最后,执行用于收集切片Ln的数据的序列Cn。通过执行序列Cn,从切片Ln收集到DC信号A1n和成像信号B1n
DC信号Aln及成像信号B1n被线圈模式Ml的通道CH1~CH4分别接收。此外,在图43中,说明的方便起见,仅示出DC信号Aln被线圈模式M1的通道CH1~CH4分别接收的情况。通道CH1~CH4分别输出信号A1n,1~A1n,4
期间P1中执行序列C1~Cn后,如以下那样生成呼吸信号。
首先,通道选择单元82(参照图41)参照数据库(参照图40)。而且,通道选择单元82基于数据库的项目c的信息,从线圈模式Ml的通道CH1~CH4中选择作为配置在肝脏的端部E1附近的通道而登记的通道CH1及CH2。
接着,呼吸信号生成单元83(参照图41)丢弃线圈模式Ml的通道CH1~CH4中未被选择的通道CH3及CH4的输出信号,仅合成(加法运算)所选择的通道CH1及CH2的输出信号。由此,得到合成信号A1
在取得合成信号A1后,呼吸信号生成单元83算出合成信号A1的积分值S1。积分值S1用作为期间P1中的被检体的呼吸信号的信号值。
以下同样如此,在期间P2~Pm中也执行序列C1~Cn。呼吸信号生成单元83丢弃通道CH3及CH4的输出信号,仅合成(加法运算)所选择的通道CH1及CH2的输出信号。而且,算出合成信号的积分值。因此,能够求出期间P1~Pm中的呼吸信号。
而且,与第1方式同样,基于期间P1~Pm的呼吸信号,设定接受成像信号B的容许范围AW(参照图27),在呼吸信号不包括在容许范围AW的情况下,再收集数据,并结束流程。
第3方式中,按每个线圈模式,使配置在肝脏的端部E1附近的通道对应(参照图40)。因此,无论使用哪个线圈模式,都能生成反映肝脏的活动的良好的呼吸信号。
(4)第4方式
第4方式中,与第3方式同样,对线圈4具有线圈模式Ml~M3的情况进行说明。但是,第4方式中,对不在数据库登记通道,而利用第2方式的预扫描PS(参照图32)选择通道的例子进行说明。此外,线圈4以外的MR装置的硬件构成与第1方式(参照图1)相同。
图44是第4方式中处理器执行的处理的说明图。
处理器8通过读出存储在存储器9的程序,构成线圈模式选择单元80~判断单元84等。
线圈模式选择单元80基于从操作部10输入的信息,从线圈模式M1~M3中选择摄影所使用的线圈模式。
切片设定单元81基于从操作部10输入的信息,设定切片。
曲线作成单元811基于通过预扫描PS而收集的MR信号,作成表示摄影部位的z方向上的各位置和信号强度的关系的曲线。
通道选择单元82基于曲线作成单元811作成的曲线,从所选择的线圈模式所具有的通道中,选择配置在肝脏的端部E1(参照图3)附近的通道。
呼吸信号生成单元83基于通道选择单元82选择的通道的输出信号,生成呼吸信号。
判断单元84判断是否将成像信号接受为图像重构的信号。
处理器8是构成线圈模式选择单元80~判断单元84的一个例子,通过执行既定程序,作为这些单元起作用。
以下,对第4方式中的MR装置的动作流程进行说明。
图45是示出第4方式中的MR装置的动作流程的图。
步骤ST0中,选择线圈模式。在第4方式中也与第3方式同样,设为选择线圈模式Ml。在选择线圈模式Ml后,进入步骤ST1。
步骤ST1及步骤ST2与第3方式相同,因此省略详细的说明。在步骤ST2中设定切片L1~Ln(参照图8)后,进入步骤ST21。
步骤ST21中,利用线圈模式Ml执行预扫描PS。
图46是预扫描PS的说明图。
预扫描PS中,仅激励切片L1~Ln中的任意一块的切片,从被激励的一块切片收集DC信号A0及成像信号B0。第2方式中,设为激励切片L1~Ln中的中央的切片Lc。因此,从切片Lc收集到DC信号A0和成像信号B0
第4方式中,选择线圈模式Ml,因此DC信号A0及成像信号B0被线圈模式M1的通道CH1~CH4分别接收。在图46中,说明的方便起见,仅示出成像信号B0被线圈模式Ml的通道CH1~CH4分别接收的情况。此外,DC信号A0及成像信号B0中,通道的选择所使用的信号为成像信号B0,而DC信号A0不在通道的选择中使用。
通道CH1~CH4接收成像信号B0,从而通道CH1~CH4分别输出信号B01~B04
在执行预扫描PS后,进入步骤ST22。
步骤ST22中,曲线作成单元811(参照图44)将通道CH1~CH8的输出信号B01~B08分别沿z方向进行傅里叶变换(FT:Fourier Transformation)。由此,如图47所示,作成曲线F1~F4。
如果作成了曲线F1~F4,就会按每个曲线算出积分值Sb和Sa的比,并算出积分值的比。图48中以标号“J1”~“J4”表示曲线F1~F4。
通道选择单元82(参照图44)以值的从大到小的顺序排列比J1~J4,按值的从大到小的顺序确定两个比。由此,能够从通道CH1~CH4中选择配置在肝脏的端部E1附近的通道CH1及CH2。
在选择通道后,进入步骤ST3。
步骤ST3中,执行正式扫描MS。第4方式中的正式扫描MS以与第3方式中的正式扫描MS(参照图43)同样的顺序执行。
第4方式中,与第3方式同样,无论使用哪个线圈模式,都能生成反映肝脏的活动的良好的呼吸信号。
另外,在第4方式中,执行预扫描PS,基于通过预扫描PS而收集的MR信号,选择配置在肝脏的端部E1附近的通道。因此,无需在数据库按摄影所使用的每个线圈模式登记通道,因此也能减轻数据库的维护花费的负担。
此外,第4方式中,在预扫描PS中,从切片Lc收集磁共振信号,作成各通道的曲线。然而,也可以从与切片Lc不同的切片收集磁共振信号,从而作成各通道的曲线。另外,也可以从多个切片收集磁共振信号,从而作成各通道的曲线。进而,第4方式中,预扫描PS为2D扫描,但也可为3D扫描。
另外,第3及第4方式中,基于操作员从操作部10输入的信息,线圈模式选择单元80选择线圈模式。然而,也可以使用自动线圈选择(Auto Coil Selection)的方法,使线圈模式选择单元自动选择线圈模式。
第1~第4方式中,通过相加通道所接收的信号,得到合成信号。然而,信号的合成并不限定于相加,例如,既可以加权相加信号而得到合成信号,也可以通过将信号相乘而得到合成信号。进而,第1~第4方式中,采用合成信号的积分值作为呼吸信号的信号值。然而,也可以采用与合成信号的积分值不同的值(例如,合成信号的最大值)作为呼吸信号的信号值。
在第1~第4方式中,基于表示k空间的中心的数据的DC信号作成呼吸信号。然而,也可以基于与DC信号不同的MR信号生成呼吸信号。
在第1~第4方式中,正式扫描MS为2D扫描,但也可为3D扫描。
另外,在第1~第4方式中,示出了取得呼吸信号的例子。然而,本发明不限定于呼吸信号的取得。例如,在摄影心脏的情况下,能够取得包含心脏的跳动的信息的生物信号。
标号说明
2 磁体;3 工作台;3a 托架;4 线圈;5 发送器;6梯度磁场电源;7 控制部;8 处理器;9 存储器;10 操作部;11 显示部;12 被检体;21 膛;80 线圈模式选择单元;81 切片设定单元;82 通道选择单元;83 呼吸信号生成单元;84 判断单元;100 MR装置。

Claims (24)

1. 一种磁共振装置,执行用于从包含活动部位的摄影部位产生第1磁共振信号的扫描,其中包括:
线圈,具有接收所述第1磁共振信号的多个通道;
通道选择单元,从所述多个通道之中,选择配置在所述活动部位的端部附近的第1通道;以及
生成单元,基于所述第1通道接收的所述第1磁共振信号,生成包含所述扫描中的表示所述摄影部位的活动的活动信息的生物信号。
2. 如权利要求1所述的磁共振装置,其中,
所述通道选择单元,基于包含用于从所述多个通道中确定配置在所述活动部位的端部附近的通道的信息的数据库,从所述多个通道中选择所述第1通道。
3. 如权利要求2所述的磁共振装置,其中,
所述通道选择单元基于所述数据库,从所述多个通道中选择至少两个第1通道,
所述生成单元合成所述至少两个第1通道接收的所述第1磁共振信号,并基于经合成而得到的合成信号,生成所述生物信号。
4. 如权利要求1~3中的任一项所述的磁共振装置,其中,
所述多个通道包含配置在相比所述第1通道远离所述活动部位的端部的位置的第2通道,
所述第2通道接收的所述第1磁共振信号不用于生成所述生物信号。
5. 如权利要求1~4中的任一项所述的磁共振装置,其中,
所述线圈具有多个由所述多个通道中的至少一个以上的通道构成的线圈模式。
6. 如权利要求5所述的磁共振装置,其中包括:
线圈模式选择单元,从所述线圈所具有的多个线圈模式中,选择用于接收所述第1磁共振信号的第1线圈模式,
所述第1线圈模式具有所述第1通道。
7. 如权利要求6所述的磁共振装置,其中,
所述第1线圈模式具有配置在相比所述第1通道远离所述活动部位的端部的位置的第2通道,
所述第2通道接收的所述第1磁共振信号不用于生成所述生物信号。
8. 如权利要求1所述的磁共振装置,其中包括:
扫描部,在所述扫描之前,执行用于从所述摄影部位产生第2磁共振信号的其他扫描;以及
曲线作成单元,基于所述其他扫描中按所述线圈的每个通道收集的所述第2磁共振信号,按每个所述通道作成表示与所述摄影部位的既定方向相关的各位置的信号值的曲线,
所述通道选择单元基于所述曲线,从所述多个通道中选择所述第1通道。
9. 如权利要求8所述的磁共振装置,其中,
所述通道选择单元基于所述曲线,从所述多个通道中选择至少两个第1通道,
所述生成单元合成所述至少两个第1通道接收的所述第1磁共振信号,基于经合成而得到的合成信号,生成所述生物信号。
10. 如权利要求8或9所述的磁共振装置,其中,
所述多个通道包含配置在相比所述第1通道远离所述活动部位的端部的位置的第2通道,
所述第1通道和所述第2通道的所述既定方向的位置不同。
11. 如权利要求8~10中的任一项所述的磁共振装置,其中,
所述通道选择单元算出表示所述曲线的特性的特性值,基于所述特性值,选择所述第1通道。
12. 如权利要求11所述的磁共振装置,其中,
所述通道选择单元将所述曲线的所述既定方向的范围分为第1范围和第2范围,算出所述第1范围上的第1积分值和所述第2范围上的第2积分值,
基于所述第1积分值和所述第2积分值算出所述特性值。
13. 如权利要求8~12中的任一项所述的磁共振装置,其中,
所述其他扫描从横切所述摄影部位的切片产生所述第2磁共振信号,
所述切片与所述既定方向平行。
14. 如权利要求13所述的磁共振装置,其中,所述切片为矢状切片。
15. 如权利要求8~14中的任一项所述的磁共振装置,其中,
所述线圈具有多个由所述多个通道中的至少一个以上的通道构成的线圈模式。
16. 如权利要求15所述的磁共振装置,其中包括:
线圈模式选择单元,从所述线圈所具有的多个线圈模式中,选择用于接收所述第1磁共振信号及所述第2磁共振信号的第1线圈模式,
所述第1线圈模式具有所述第1通道。
17. 如权利要求16所述的磁共振装置,其中,
所述第1线圈模式具有配置在相比所述第1通道远离所述活动部位的端部的位置的第2通道,
所述第2通道接收的所述第1磁共振信号不用于生成所述生物信号。
18. 如权利要求1~17中的任一项所述的磁共振装置,其中,
所述扫描从所述摄影部位产生用于重构所述摄影部位的图像的第3磁共振信号。
19. 如权利要求18所述的磁共振装置,其中包括:
判断单元,基于所述生物信号,判断是否将所述第3磁共振信号接受为图像重构的信号。
20. 如权利要求1~19中的任一项所述的磁共振装置,其中,
所述第1磁共振信号是表示k空间的中心的数据的信号。
21. 如权利要求1~20中的任一项所述的磁共振装置,其中,
所述活动部位为肝脏,
所述活动部位的端部为肝脏的肺侧的端部。
22. 如权利要求1~20中的任一项所述的磁共振装置,其中,
所述活动部位为肝脏,
所述活动部位的端部是与肝脏的肺侧相反侧的端部。
23. 如权利要求1~22中的任一项所述的磁共振装置,其中,
所述生物信号为呼吸信号。
24. 一种适用于磁共振装置的程序,其中所述磁共振装置具备:扫描部,执行用于从包含活动部位的摄影部位产生第1磁共振信号的扫描;以及线圈,具有接收所述第1磁共振信号的多个通道,所述程序使计算机执行以下处理:
通道选择处理,从所述多个通道之中,选择配置在所述活动部位的端部附近的第1通道;以及
生成处理,基于所述第1通道接收的所述第1磁共振信号,生成包含表示所述扫描中的所述摄影部位的活动的活动信息的生物信号。
CN201510089388.1A 2014-02-28 2015-02-27 磁共振装置及程序 Active CN104873204B (zh)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
JP2014-039390 2014-02-28
JP2014039390A JP6489744B2 (ja) 2014-02-28 2014-02-28 磁気共鳴装置およびプログラム

Publications (2)

Publication Number Publication Date
CN104873204A true CN104873204A (zh) 2015-09-02
CN104873204B CN104873204B (zh) 2019-09-17

Family

ID=53941081

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510089388.1A Active CN104873204B (zh) 2014-02-28 2015-02-27 磁共振装置及程序

Country Status (3)

Country Link
US (1) US20150247911A1 (zh)
JP (1) JP6489744B2 (zh)
CN (1) CN104873204B (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106725508A (zh) * 2016-12-15 2017-05-31 上海联影医疗科技有限公司 生理运动数据采集方法、磁共振成像方法以及装置
CN108882889A (zh) * 2016-03-31 2018-11-23 通用电气公司 磁共振成像装置和程序
CN109009113A (zh) * 2018-08-03 2018-12-18 上海联影医疗科技有限公司 运动检测装置、方法和磁共振成像系统、方法

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP5755285B2 (ja) * 2013-05-31 2015-07-29 ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー 磁気共鳴装置およびプログラム
JP6593688B2 (ja) 2015-08-20 2019-10-23 日立オートモティブシステムズ株式会社 ブレーキ装置およびブレーキシステム
JP6742758B2 (ja) * 2016-02-29 2020-08-19 ゼネラル・エレクトリック・カンパニイ 磁気共鳴装置およびプログラム
US10859646B2 (en) 2018-05-31 2020-12-08 General Electric Company Method and systems for coil selection in magnetic resonance imaging to reduce annefact artifact
US10866292B2 (en) 2018-05-31 2020-12-15 General Electric Company Methods and systems for coil selection in magnetic resonance imaging
US10859645B2 (en) 2018-05-31 2020-12-08 General Electric Company Method and systems for coil selection in magnetic resonance imaging
US10802101B2 (en) 2018-05-31 2020-10-13 General Electric Company Method and systems for coil selection in magnetic resonance imaging to reduce phase wrap artifact

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1821804A (zh) * 2005-02-15 2006-08-23 株式会社东芝 磁共振成像装置及其控制方法
CN101268943A (zh) * 2007-03-19 2008-09-24 株式会社东芝 磁共振摄像装置及方法、磁共振摄像维护装置及系统
JP2008283994A (ja) * 2007-05-15 2008-11-27 Hitachi Medical Corp 磁気共鳴イメージング装置
US20090066329A1 (en) * 2007-08-28 2009-03-12 Kenichi Kanda Magnetic resonance imaging apparatus and navigator data analyzing method
JP2010075573A (ja) * 2008-09-29 2010-04-08 Ge Medical Systems Global Technology Co Llc 磁気共鳴イメージング装置
CN101843484A (zh) * 2009-01-05 2010-09-29 Ge医疗系统环球技术有限公司 磁共振成像装置
CN102349830A (zh) * 2008-06-26 2012-02-15 株式会社东芝 磁共振成像装置及磁共振成像方法

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6798199B2 (en) * 2003-02-06 2004-09-28 Siemens Medical Solutions Usa, Inc. Method for synchronizing magnetic resonance imaging data to body motion
JP5405732B2 (ja) * 2007-10-30 2014-02-05 株式会社日立メディコ 磁気共鳴イメージング装置
WO2009134820A2 (en) * 2008-04-28 2009-11-05 Cornell University Tool for accurate quantification in molecular mri
JP5492033B2 (ja) * 2010-09-07 2014-05-14 株式会社日立メディコ 核磁気共鳴イメージング装置および表面コイルの選択方法
JP5908868B2 (ja) * 2013-06-28 2016-04-26 ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー 磁気共鳴装置およびプログラム

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1821804A (zh) * 2005-02-15 2006-08-23 株式会社东芝 磁共振成像装置及其控制方法
CN101268943A (zh) * 2007-03-19 2008-09-24 株式会社东芝 磁共振摄像装置及方法、磁共振摄像维护装置及系统
JP2008283994A (ja) * 2007-05-15 2008-11-27 Hitachi Medical Corp 磁気共鳴イメージング装置
US20090066329A1 (en) * 2007-08-28 2009-03-12 Kenichi Kanda Magnetic resonance imaging apparatus and navigator data analyzing method
CN102349830A (zh) * 2008-06-26 2012-02-15 株式会社东芝 磁共振成像装置及磁共振成像方法
JP2010075573A (ja) * 2008-09-29 2010-04-08 Ge Medical Systems Global Technology Co Llc 磁気共鳴イメージング装置
CN101843484A (zh) * 2009-01-05 2010-09-29 Ge医疗系统环球技术有限公司 磁共振成像装置

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108882889A (zh) * 2016-03-31 2018-11-23 通用电气公司 磁共振成像装置和程序
CN106725508A (zh) * 2016-12-15 2017-05-31 上海联影医疗科技有限公司 生理运动数据采集方法、磁共振成像方法以及装置
CN106725508B (zh) * 2016-12-15 2021-09-28 上海联影医疗科技股份有限公司 生理运动数据采集方法、磁共振成像方法以及装置
CN109009113A (zh) * 2018-08-03 2018-12-18 上海联影医疗科技有限公司 运动检测装置、方法和磁共振成像系统、方法

Also Published As

Publication number Publication date
US20150247911A1 (en) 2015-09-03
CN104873204B (zh) 2019-09-17
JP6489744B2 (ja) 2019-03-27
JP2015163109A (ja) 2015-09-10

Similar Documents

Publication Publication Date Title
CN104873204A (zh) 磁共振装置及程序
US10222443B2 (en) Generating motion information for an at least partially moving examination region
US10444315B2 (en) MRI with motion correction using navigators acquired using a dixon technique
US9476955B2 (en) Method and control device to operate a magnetic resonance system
US9575153B2 (en) MR imaging using a multi-point dixon technique
EP1015903B1 (en) Simultaneous acquisition of spatial basis functions: ultra-fast imaging with radiofrequency coil arrays
US8116541B2 (en) Method and apparatus for multi-coil magnetic resonance imaging
CN103238082B (zh) 使用多点Dixon技术和低分辨率校准的MR成像
EP2496954B1 (en) Mr imaging using navigators
CN105103001A (zh) 使用高snr同相图像和较低snr至少部分地异相图像的dixon型水/脂肪分离mri
CN102870000A (zh) 使用单个和多个通道接收器线圈的用于同时多切片磁共振成象的方法
CN103211596A (zh) 用于运动修正mr扩散成像的系统
CN107076818B (zh) 零回波时间mr成像
CN106456046A (zh) 经改进的多时相动态对比增强磁共振成像的方法
WO2012038886A1 (en) Magnetic resonance imaging of chemical species
CN105301537A (zh) 同时磁共振成像方法和用于同时多核磁共振成像的装置
JP2019535434A (ja) ディクソンタイプの水/脂肪分離を用いたmrイメージング
JP6363126B2 (ja) 磁気共鳴イメージング装置およびプログラム
JP6742758B2 (ja) 磁気共鳴装置およびプログラム
JP6356412B2 (ja) 磁気共鳴装置およびプログラム
Sodickson Spatial encoding using multiple RF coils: SMASH imaging and parallel MRI
JP6685097B2 (ja) 磁気共鳴装置
JP6554723B2 (ja) 磁気共鳴装置およびプログラム

Legal Events

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