CN107369153B - 磁共振脂肪抑制图像获取方法、存储介质及扫描*** - Google Patents

磁共振脂肪抑制图像获取方法、存储介质及扫描*** Download PDF

Info

Publication number
CN107369153B
CN107369153B CN201710588221.9A CN201710588221A CN107369153B CN 107369153 B CN107369153 B CN 107369153B CN 201710588221 A CN201710588221 A CN 201710588221A CN 107369153 B CN107369153 B CN 107369153B
Authority
CN
China
Prior art keywords
radio frequency
pulse
fat
angle alpha
pressing
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
CN201710588221.9A
Other languages
English (en)
Other versions
CN107369153A (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.)
Shanghai United Imaging Healthcare Co Ltd
Original Assignee
Shanghai United Imaging Healthcare Co 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 Shanghai United Imaging Healthcare Co Ltd filed Critical Shanghai United Imaging Healthcare Co Ltd
Priority to CN201710588221.9A priority Critical patent/CN107369153B/zh
Priority to US15/798,147 priority patent/US10578699B2/en
Publication of CN107369153A publication Critical patent/CN107369153A/zh
Priority to PCT/CN2017/114102 priority patent/WO2019015212A1/en
Priority to EP17204652.6A priority patent/EP3432020A1/en
Application granted granted Critical
Publication of CN107369153B publication Critical patent/CN107369153B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/0002Inspection of images, e.g. flaw detection
    • G06T7/0012Biomedical image inspection
    • 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/5607Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution by reducing the NMR signal of a particular spin species, e.g. of a chemical species for fat suppression, or of a moving spin species for black-blood imaging
    • 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/565Correction of image distortions, e.g. due to magnetic field inhomogeneities
    • G01R33/5659Correction of image distortions, e.g. due to magnetic field inhomogeneities caused by a distortion of the RF magnetic field, e.g. spatial inhomogeneities of the RF magnetic field
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10072Tomographic images
    • G06T2207/10088Magnetic resonance imaging [MRI]

Landscapes

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

Abstract

本发明涉及一种磁共振脂肪抑制图像获取方法,包括以下步骤:获取单射频压脂脉冲的射频脉冲角度α;通过单射频压脂脉冲的射频脉冲角度α确定双射频压脂脉冲的第一射频脉冲角度α1与第二射频脉冲角度α2;以及基于所述双射频压脂脉冲的第一射频脉冲角度α1与第二射频脉冲角度α2获取脂肪抑制图像。本发明的磁共振脂肪抑制图像获取方法可以显著降低压脂均匀性对射频发射场的敏感性。本发明还涉及一种存储介质及扫描***。

Description

磁共振脂肪抑制图像获取方法、存储介质及扫描***
技术领域
本发明涉及医疗设备领域,尤其涉及一种磁共振脂肪抑制图像获取方法、存储介质及扫描***。
背景技术
传统的频率选择性脂肪抑制方法,通过在选层激发成像前,施加一个适当带宽的射频软脉冲,激发脂肪信号并梯度散相,使得水信号得以保留而脂肪信号受到抑制。压脂均匀性受射频发射场(即B1)不均匀的影响很大,亦即传统方法压脂均匀性对B1很敏感。在T2FSE(fast spin echo,快速自旋回波)成像中,B1不均匀对脂肪抑制均匀性的影响往往更显著。3T***的B1不均匀明显,压脂均匀性成为影响图像质量的一个严峻问题。
发明内容
基于此,提供一种磁共振脂肪抑制图像获取方法、存储介质及扫描***,以显著降低压脂均匀性对射频发射场的敏感性。
一种磁共振脂肪抑制图像获取方法,包括以下步骤:
获取单射频压脂脉冲的射频脉冲角度α;
通过单射频压脂脉冲的射频脉冲角度α确定双射频压脂脉冲的第一射频脉冲角度α1与第二射频脉冲角度α2;以及
基于所述双射频压脂脉冲的第一射频脉冲角度α1与第二射频脉冲角度α2获取脂肪抑制图像。
进一步地,所述获取单射频压脂脉冲的射频脉冲角度α的步骤包括:
根据单射频压脂脉冲的相邻压脂模块的时间间隔Tf,脂肪纵向恢复特征时间T1,脂肪残余量r以及压脂脉冲等效中心距离成像序列激发脉冲等效中心的时间间隔τ,确定出单射频压脂脉冲的射频脉冲角度α:
Figure BDA0001353571760000021
其中
Figure BDA0001353571760000022
进一步地,所述通过单射频压脂脉冲的射频脉冲角度α确定双射频压脂脉冲的第一射频脉冲角度α1与第二射频脉冲角度α2的步骤包括:
获取单射频压脂脉冲的射频脉冲角度α与双射频压脂脉冲的第一射频脉冲角度α1、第二射频脉冲角度α2之间的映射关系,以确定第一射频脉冲角度α1与第二射频脉冲角度α2的预选角度值。
进一步地,所述获取单射频压脂脉冲的射频脉冲角度α与双射频压脂脉冲的第一射频脉冲角度α1、第二射频脉冲角度α2之间的映射关系的步骤包括:
通过以脂肪纵向磁矩的作用与单射频压脂脉冲作用相当,确定出单射频压脂脉冲的射频脉冲角度α与双射频压脂脉冲的第一射频脉冲角度α1、第二射频脉冲角度α2之间的映射关系:cos(α1)cos(α2)=cos(α)。
进一步地,所述通过单射频压脂脉冲的射频脉冲角度α确定双射频压脂脉冲的第一射频脉冲角度α1与第二射频脉冲角度α2的步骤还包括:
限定双射频压脂脉冲的第一射频脉冲角度α1、第二射频脉冲角度α2之间的关系,以从预选角度值中筛选出第一射频脉冲角度α1、第二射频脉冲角度α2对射频发射场最不敏感的角度值。
进一步地,基于以下公式确定所述第一射频脉冲角度α1、第二射频脉冲角度α2对射频发射场最不敏感的角度值:
Figure BDA0001353571760000023
即α1sin(α1)cos(α2)+α2cos(α1)sin(α2)=0,
其中,λ代表射频发射场的均匀性,λ<1表示射频发射场的均匀偏弱,λ>1射频发射场的均匀偏强设。
进一步地,在确定出第一射频脉冲角度α1、第二射频脉冲角度α2对射频发射场最不敏感的角度值之后还包括以下步骤:
存储每个单射频压脂脉冲的射频脉冲角度α对应的双单射频压脂脉冲的第一射频脉冲角度α1、第二射频脉冲角度α2对射频发射场最不敏感的角度值,形成映射表。
进一步地,还包括:
再次获取单射频压脂脉冲的射频脉冲角度α;
通过映射表确定所述第一射频脉冲角度α1、第二射频脉冲角度α2对射频发射场最不敏感的角度值。
一种扫描***,包括磁共振扫描设备和计算机,其中计算机包括处理器、存储介质及存储在存储介质上并可在处理器上运行的计算机程序,所述处理器执行所述程序时可用于执行一种磁共振脂肪抑制图像获取方法,所述方法包括以下步骤:
获取单射频压脂脉冲的射频脉冲角度α;
通过单射频压脂脉冲的射频脉冲角度α确定双射频压脂脉冲的第一射频脉冲角度α1与第二射频脉冲角度α2;以及
基于所述双射频压脂脉冲的第一射频脉冲角度α1与第二射频脉冲角度α2获取脂肪抑制图像。
一种存储介质,其上存储有计算机程序,计算机可读取该存储介质内存储的计算机程序,该程序被处理器执行可用于执行上述磁共振脂肪抑制图像获取方法。
本发明的磁共振脂肪抑制图像获取方法、存储介质及扫描***相对传统单射频压脂脉冲压脂具有通过双射频压脂脉冲压脂,能够在基本不增加扫描时间的基础上显著降低压脂均匀性对射频发射场的敏感性。
附图说明
图1是本发明实施例的扫描***的方框图。
图2是本发明实施例的扫描***的结构示意图。
图3本发明实施例的磁共振脂肪抑制图像获取方法。
图4是单射频压脂脉冲与相应的梯度信号以及主扫面序列的波形图,其中,主扫面序列位于虚线方框内,虚线方框以外的部分为单射频压脂脉冲与相应的梯度信号。
图5是本发明的双射频压脂脉冲与相应的梯度信号以及主扫面序列的波形图,其中虚线左侧为双射频压脂脉冲与相应的梯度信号,虚线右侧为主扫面序列。
图6是利用射频压脂脉冲压脂和利用传统单射频压脂脉冲腰肢,脂肪残余对B1的敏感性图。
图7与图8在一扫描***中利用双射频压脂脉冲与单射频压脂脉冲分别扫描腹部横断与盆腔横断位的效果对比图。
图9是在一扫描***中利用双射频压脂脉冲与单射频压脂脉冲分别扫描髋关节横断位的效果对比图。
图10与图11是在一扫描***中利用双射频压脂脉冲与单射频压脂脉冲分别扫描腹部横断与盆腔横断位。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅用以解释本发明,并不用于限定本发明。
请参阅图1与图2,本发明实施例提供一种扫描***100,包括磁共振扫描设备10及与磁共振扫描设备10连接的计算机12,该计算机12中存储有可以运行的计算机程序,该计算机程序在运行时可以运行本发明所提供的一种磁共振脂肪抑制图像获取方法。
在本实施例中,磁共振扫描设备10包括磁共振机架,机架内有主磁体101,主磁体101可以是由超导线圈构成,用来产生主磁场,在一些情况下也可以采用永磁体。主磁体101可以用来产生0.2特斯拉、0.5特斯拉、1.0特斯拉、1.5特斯拉、3.0特斯拉或者更高的主磁场强度。在磁共振成像时,成像对象150例如患者由患者床106进行承载,随着床板的移动,将成像对象150移入主磁场磁场分布较为均匀的区域105内。通常对于磁共振扫描设备,如图2所示,空间坐标系(即设备的坐标系)的z方向设置为与磁共振***机架的轴向相同,通常将患者的身长方向与z方向保持一致进行成像,磁共振***的水平平面设置为xz平面,x方向与z方向垂直,y方向与x和z方向均垂直。
磁共振扫描设备10还包括射频线圈、梯度线圈102、脉冲控制单元111、梯度信号产生单元112、射频脉冲产生单元116、开关控制单元117、射频接收单元118及图像重建单元121。
在磁共振成像时,脉冲控制单元111控制射频脉冲产生单元116产生射频脉冲,射频脉冲由放大器放大后,经过开关控制单元117,最终由射频线圈(图未标号)发出,对成像对象150进行射频激发。射频线圈包括体线圈103与局部线圈104射频线圈在发出射频脉冲通常通过体线圈103或者局部线圈104发出。成像对象150根据射频激发,会由共振产生相应的射频信号。在接收成像对象150根据激发产生的射频信号时,可以是由体线圈103或者局部线圈104进行接收,射频接收链路可以有很多条,射频信号发送到射频接收单元118后,进一步发送到图像重建单元121进行图像重建,形成磁共振图像。
梯度线圈102可以用来在磁共振成像时对射频信号进行空间编码。脉冲控制单元111控制梯度信号产生单元112产生梯度信号,梯度信号通常会分为三个相互正交方向的信号:x方向、y方向和z方向,不同方向的梯度信号经过梯度放大器(113、114、115)放大后,由梯度线圈102发出,在区域105内产生梯度磁场。
计算机12还包括处理器122、显示单元123、输入/输出设备124、存储介质125及通信端口126。脉冲控制单元111、图像重建单元121与处理器122、显示单元123、输入/输出设备124、存储介质125及通信端口126之间可以通过通信总线127进行数据传输,从而实现对磁共振成像过程的控制。其中,处理器122可以由一个或多个处理器组成。显示单元123可以是提供给用户用来显示图像的显示器。输入/输出设备124可以是键盘、鼠标、控制盒等相关设备,支持输入/输出相应数据流。存储介质125可以是只读存储器(ROM)、随机存取存储器(RAM)、硬盘等,存储介质125可以用来存储需要处理和/或通信使用的各种数据文件,以及处理器122所执行的可能的程序指令(该可能的程序指令包括在运行时可以运行本发明所提供的一种磁共振脂肪抑制图像获取方法的指令)。通信端口105可以实现与其他部件例如:外接设备、图像采集设备、数据库、外部存储以及图像处理工作站等之间进行数据通信。
如图3所示,本发明实施例所提供的一种磁共振脂肪抑制图像获取方法,该方法包括以下步骤:
步骤S1:获取单射频压脂脉冲的射频脉冲角度α。
该过程具体包括:根据单射频压脂脉冲的相邻压脂模块的时间间隔Tf,脂肪纵向恢复特征时间T1,脂肪残余量r以及压脂脉冲等效中心距离成像序列激发脉冲等效中心的时间间隔τ,计算出射频脉冲角度α,射频脉冲角度α的具体计算公式如下:
Figure BDA0001353571760000071
其中
Figure BDA0001353571760000072
在本实施例中,脂肪残余量r的范围为-0.1~0。
需要说明的是,具体如图4所示,图4示出了单射频压脂脉冲的相邻压脂模块的时间间隔Tf、压脂脉冲等效中心距离成像序列激发脉冲等效中心的时间间隔τ及与单射频压脂脉冲相匹配的主脉冲序列,由图4可知,一般压脂序列的时序是这样的:一个压脂模块,紧接着主序列扫第1层;再一个压脂模块,紧接着主序列扫第2层;......再一个压脂模块,紧接着主序列扫第n层。相邻压脂模块时间间隔Tf,即是指对应时间上相继扫描的两层所对应的压脂模块的时间间隔。另外,单脉冲压脂模块的压脂脉冲有一定时间,比如10ms,压脂脉冲后还会有散相梯度,约2ms,再紧接是主序列的选层梯度与射频激发脉冲,约3,4ms,那么从压脂脉冲中心到主序列激发脉冲中心的时间间隔就称为τ。
需要进一步地说明的是,将脂肪残余量r作为参数的目的是:假设不压脂时,图像上脂肪的信号强度约定为1。脂肪残余量,对应着压脂图像上脂肪的亮度。因为医生并不需要把脂肪抑制得漆黑(对应脂肪残余为0),而是希望图像上的脂肪能保持一定的灰度(比如,脂肪残余为0.04)。
步骤S2:通过单射频压脂脉冲的射频脉冲角度α确定双射频压脂脉冲的第一射频脉冲角度α1与第二射频脉冲角度α2
作为其中一个实施例,该过程具体包括:首先确定单射频压脂脉冲的射频脉冲角度α与双射频压脂脉冲的第一射频脉冲角度α1、第二射频脉冲角度α2之间的映射关系,以确定第一射频脉冲角度α1、第二射频脉冲角度α2的角度预选值。
单射频压脂脉冲的射频脉冲角度α与双射频压脂脉冲的第一射频脉冲角度α1、第二射频脉冲角度α2之间的对应关系,具体可通过忽略射频脉冲的时间及两者间的散相梯度时间,以及双脉冲对脂肪纵向磁矩的作用要与单脉冲作用相当来确定,该方法确定出单射频压脂脉冲的射频脉冲角度α与双射频压脂脉冲的第一射频脉冲角度α1、第二射频脉冲角度α2之间存在如下关系:
cos(α1)cos(α2)=cos(α) (e1)。
需要说明的是,如果不忽略两射频脉冲各自的时长,以及两脉冲间的梯度时长,就要考虑脂肪的T1恢复,那么等式左边的表达式将变得很复杂。这是没有必要的,临床实验已经表明,这种简化的表达式所能达到的均匀性改进效果已经非常显著了,而且理论上也可以论证(主要是因为脉冲时间,及梯度时间,相对于脂肪的T1,是非常短的。脂肪的T1为[email protected],380ms@3T),这种简化的表达式确实是足够准确的。另外,“双脉冲的作用要与单脉冲的作用相当”,就对应公式中等号的成立。公式左边是简化处理的双脉冲对脂肪纵向磁矩的作用,公式右边是单脉冲对脂肪纵向磁矩的作用,两者作用相当,意味着等号成立。
可以理解是,在单射频压脂脉冲的射频脉冲角度α值确定的情况下,满足角度预选值的双射频压脂脉冲的第一射频脉冲角度α1与第二射频脉冲角度α2的数值有很多,所以还可以通过进一步限定双射频压脂脉冲中的第一射频脉冲角度α1与第二射频脉冲角度α2之间的关系,以从预选角度值中筛选出更优的一组第一射频脉冲角度α1、第二射频脉冲角度α2
双射频压脂脉冲中的第一射频脉冲角度α1与第二射频脉冲角度α2之间存在如下关系:
Figure BDA0001353571760000091
即α1sin(α1)cos(α2)+α2cos(α1)sin(α2)=0 (e2)
其中,λ表示射频发射场的均匀性分布,λ<1表示射频发射场的均匀偏弱,λ>1表示射频发射场的均匀偏强。
可以理解的是,通过联立e1公式与e2确定出双射频压脂脉冲的第一射频脉冲角度α1与第二射频脉冲角度α2对射频发射场(B1)最不敏感的一组值。
根据本步骤的方法,通常情况下可以计算得出第一射频脉冲角度α1的数值大于90度,第二射频脉冲角度α2的数值小于90度;或者第一射频脉冲角度α1的数值小于90度,第二射频脉冲角度α2的数值大于90度。
步骤S3:基于所述双射频压脂脉冲的第一射频脉冲角度α1、第二射频脉冲角度α2获取脂肪抑制图像。
该过程具体包括:处理器122在获取双射频压脂脉冲的第一射频脉冲角度α1与第二射频脉冲角度α2后(例如可以是对射频发射场最不敏感的第一射频脉冲角度α1与第二射频脉冲角度α2),处理器122将第一射频脉冲角度α1与第二射频脉冲角度α2发送至脉冲控制单元111控制单元,脉冲控制单元111相应的控制射频脉冲产生单元116连续无间隔地输出两个射频压脂脉冲(中间有散相梯度,散相梯度和射频脉冲紧密相连,见序列图,具体参见图5),其中一个射频压脂脉冲具有第一射频脉冲角度α1,另一个射频压脂脉冲具有第二射频脉冲角度α2,可以理解的是,具有第一射频脉冲角度α1的射频压脂脉冲需要伴随相应地散相梯度信号,具有第二射频脉冲角度α2的射频压脂脉冲也需要伴随相应地散相梯度信号,处理器122同时控制梯度信号产生单元112连续无间隔地输出两个散相梯度信号,其中一个散相梯度信号与具有第一射频脉冲角度α1的射频压脂脉冲相对应,另一个与具有第二射频脉冲角度α2的射频压脂脉冲相对应,需要说明的是,需配合主扫描序列调整与双射频压脂脉冲对应的两个散相梯度信号,以确保将横向脂肪信号散相,不产生脂肪伪影。双射频压脂脉冲、与双射频压脂脉冲相对应的散相梯度信号、主扫描序列如图5所示,竖线左侧的代表双射频压脂脉冲、与双射频压脂脉冲相对应的散相梯度信号,竖线右侧的代表主扫描序列。双射频压脂脉冲及相应地散相梯度信号经过射频线圈及梯度线圈发出,主扫描序列的射频信号和梯度信号也经射频线圈及梯度线圈发出。基于以上磁共振序列对成像对象进行扫描并获取磁共振射频信号,再经过磁共振扫描系100处理后形成脂肪抑制图像,形成脂肪抑制图像的过程为现有技术,在此不再赘述。
需要说明的是,在其他实施例中,可以将与每个单射频压脂脉冲的射频脉冲角度α所对应的双射频压脂脉冲的第一射频脉冲角度α1与第二射频脉冲角度α2存储在存储介质125中,在存储介质125中形成映射表,这样在再次获取单射频压脂脉冲的射频脉冲角度α后,处理器122可以从存储介质125中直接调取与其对应的双射频压脂脉冲的第一射频脉冲角度α1与第二射频脉冲角度α2,然后处理器122将第一射频脉冲角度α1与第二射频脉冲角度α2发送至脉冲控制单元111控制单元,脉冲控制单元111相应的控制射频脉冲产生单元116连续无间隔地输出两个射频压脂脉冲,其中一个射频压脂脉冲具有第一射频脉冲角度α1,另一个射频压脂脉冲具有第二射频脉冲角度α2
在本实施例中,以1.5T磁共振***,FSE序列的脂肪抑制为例,脂肪T1取260ms,设TR/Slice=2700ms/12,TR(repetition time)表示重复时间,即射频压脂脉冲每Tf=225ms施加一次。分析得到采用双射频压脂脉冲压脂和传统单射频压脂脉冲压脂,脂肪残余对B1的敏感性如下图6所示,图6中横轴表示射频发射场的不均匀性,数值>1表示射频发射场偏大,1.5表示偏大50%;数值<1表示发射场偏小,0.6表示偏低40%。纵轴表示脂肪纵向磁矩的残余,曲线越平缓,表示在相同的B1不均匀范围内,脂肪纵向残余范围越小,对应脂肪抑制越均匀;曲线越陡峭,表示在相同的B1不均匀范围内,脂肪纵向残余范围越大,对应脂肪抑制越不均匀。其中,图6中虚线表示的曲线代表的是利用单射频压脂脉冲压脂,脂肪抑制对射频发射场的敏感度;图6中实线表的曲线示代表的是利用双射频压脂脉冲压脂,脂肪抑制对射频发射场的敏感度。由图6可知,本发明具有以下有益效果:
(1)、利用双射频压脂脉冲压脂,压脂均匀性对B1的敏感性显著降低,而利用原单射频压脂脉冲压脂,脂肪残余随B1变化迅速,亦即轻微的B1不均匀,就会导致脂肪残余的显著变化,从而图像很容易表现为压脂不均匀。通过双射频压脂脉冲压脂,脂肪残余随B1的变化要缓慢,亦即脂肪残余对B1不敏感,从而图像的压脂均匀性有所提高。定量来比较,以[0,+5%]脂肪残余为衡量标准,利用传统单射频压脂脉冲压脂,压脂允许的B1不均匀性为[-7%,0%];而本发明允许B1的程度增加为[-21%,+24%]。在同样的脂肪残余范围内,采用双射频压脂脉冲压脂,B1不均匀容忍度有显著提高。
(2)、利用双射频压脂脉冲压脂,脂肪残余恒为正,从而使脂肪抑制表现得更均匀。利用传统单射频压脂脉冲压脂,在B1偏强时,脂肪残余为负,B1偏弱时,脂肪残余为正。B1场强弱交替会带来脂肪残余正负交替,使得压脂不均匀表现特别明显。本发明使得无论B1偏强还是偏弱,脂肪残余始终为正,不存在脂肪残余随B1场变化而正负交替的现象,压脂均匀性从而显著改善。
另请参见本发明的具体应用,以更直接对比出本发明通过双射频压脂脉冲压脂相对于传统的采用单射频压脂脉冲压脂的优势。
具体应用一
具体应用一是利用双射频压脂脉冲与单射频压脂脉冲分别扫描腹部横断与盆腔横断位。具体请参阅图7与图8,图7是利用双射频压脂脉冲与单射频压脂脉冲分别扫描腹部横断位的效果对比图,图8是利用双射频压脂脉冲与单射频压脂脉冲分别扫描盆腔横断位的效果对比图。
图7与图8中的图像均使用上海联影医疗科技有限公司的1.5T磁共振***,T2FSE压脂扫描获得。传统单射频压脂脉冲压脂与本发明双射频压脂脉冲压脂的等效带宽和频率偏置基本相同,以确保压脂均匀性的差别都来自对B1不均匀的敏感程度的不同。图7与图8中左图均代表传统单射频压脂脉冲压脂效果,右图均代表本发明中的双射频压脂脉冲压脂效果。可以看到,对于B1不均匀较为明显的大范围体部扫描,皮下脂肪抑制的均匀性改善明显。
具体应用二
具体应用二是利用双射频压脂脉冲与单射频压脂脉冲分别扫描髋关节横断位。具体请参阅图9,图9是利用双射频压脂脉冲与单射频压脂脉冲分别扫描髋关节横断位的效果对比图。
请参阅图9,图9中的图像均使用上海联影医疗科技有限公司的3.0T磁共振***,T2FSE压脂扫描获得。传统单射频压脂脉冲压脂与本发明中双射频压脂脉冲压脂的等效带宽和频率偏置基本相同,以确保压脂均匀性的差别都来自对B1不均匀的敏感程度的不同。在图9中,左图均代表传统单射频压脂脉冲效果,右图均代表本发明射频压脂脉冲效果。可以看到,对于B1不均匀较为明显的大范围体部扫描,皮下脂肪抑制的均匀性改善明显。
具体应用三
具体应用三是利用双射频压脂脉冲与单射频压脂脉冲分别扫描腹部横断与盆腔横断位。具体请参阅图10与图11,图10是利用双射频压脂脉冲与单射频压脂脉冲分别扫描腹部横断位的效果对比图,图11是利用双射频压脂脉冲与单射频压脂脉冲分别扫描盆腔横断位的效果对比图。
图10与图11中的图像均使用上海联影医疗科技有限公司的1.5T磁共振***,T2FSE压脂扫描获得。传统单射频压脂脉冲压脂与本发明双射频压脂脉冲压脂的等效带宽和频率偏置基本相同,以确保压脂均匀性的差别都来自对B1不均匀的敏感程度的不同。图10与图11中左图均代表传统单射频压脂脉冲压脂效果,右图均代表本发明中的双射频压脂脉冲压脂效果。可以看到,对于B1不均匀较为明显的大范围体部扫描,皮下脂肪抑制的均匀性改善明显。
以上所述实施例仅表达了本发明的几种实施方式,其描述较为具体和详细,但并不能因此而理解为对发明专利范围的限制。应当指出的是,对于本领域的普通技术人员来说,在不脱离本发明构思的前提下,还可以做出若干变形和改进,这些都属于本发明的保护范围。因此,本发明专利的保护范围应以所附权利要求为准。

Claims (8)

1.一种基于双射频压脂脉冲的磁共振脂肪抑制图像获取方法,其特征在于,包括以下步骤:
获取单射频压脂脉冲的射频脉冲角度α;
通过单射频压脂脉冲的射频脉冲角度α确定双射频压脂脉冲的第一射频脉冲角度α1与第二射频脉冲角度α2;以及
基于所述双射频压脂脉冲的第一射频脉冲角度α1与第二射频脉冲角度α2获取脂肪抑制图像;
所述通过单射频压脂脉冲的射频脉冲角度α确定双射频压脂脉冲的第一射频脉冲角度α1与第二射频脉冲角度α2的步骤包括:
获取单射频压脂脉冲的射频脉冲角度α与双射频压脂脉冲的第一射频脉冲角度α1、第二射频脉冲角度α2之间的映射关系,以确定第一射频脉冲角度α1与第二射频脉冲角度α2的预选角度值;
所述获取单射频压脂脉冲的射频脉冲角度α与双射频压脂脉冲的第一射频脉冲角度α1、第二射频脉冲角度α2之间的映射关系的步骤包括:
通过以脂肪纵向磁矩的作用与单射频压脂脉冲作用相当,确定出单射频压脂脉冲的射频脉冲角度α与双射频压脂脉冲的第一射频脉冲角度α1、第二射频脉冲角度α2之间的映射关系:cos(α1)cos(α2)=cos(α)。
2.如权利要求1所述的磁共振脂肪抑制图像获取方法,其特征在于,所述获取单射频压脂脉冲的射频脉冲角度α的步骤包括:
根据单射频压脂脉冲的相邻压脂模块的时间间隔Tf,脂肪纵向恢复特征时间T1,脂肪残余量r以及压脂脉冲等效中心距离成像序列激发脉冲等效中心的时间间隔τ,确定出单射频压脂脉冲的射频脉冲角度α:
Figure FDA0002697063160000021
其中
Figure FDA0002697063160000022
3.如权利要求1所述的磁共振脂肪抑制图像获取方法,其特征在于,所述通过单射频压脂脉冲的射频脉冲角度α确定双射频压脂脉冲的第一射频脉冲角度α1与第二射频脉冲角度α2的步骤还包括:
限定双射频压脂脉冲的第一射频脉冲角度α1、第二射频脉冲角度α2之间的关系,以从预选角度值中筛选出第一射频脉冲角度α1、第二射频脉冲角度α2对射频发射场最不敏感的角度值。
4.如权利要求3所述的磁共振脂肪抑制图像获取方法,其特征在于,基于以下公式确定所述第一射频脉冲角度α1、第二射频脉冲角度α2对射频发射场最不敏感的角度值:
Figure FDA0002697063160000023
即α1sin(α1)cos(α2)+α2cos(α1)sin(α2)=0,
其中,λ代表射频发射场的均匀性,λ<1表示射频发射场偏弱,λ>1表示射频发射场偏强。
5.如权利要求3所述的磁共振脂肪抑制图像获取方法,其特征在于,在确定出第一射频脉冲角度α1、第二射频脉冲角度α2对射频发射场最不敏感的角度值之后还包括以下步骤:
存储每个单射频压脂脉冲的射频脉冲角度α对应的双单射频压脂脉冲的第一射频脉冲角度α1、第二射频脉冲角度α2对射频发射场最不敏感的角度值,形成映射表。
6.如权利要求5所述的磁共振脂肪抑制图像获取方法,其特征在于,还包括:
再次获取单射频压脂脉冲的射频脉冲角度α;
通过映射表确定所述第一射频脉冲角度α1、第二射频脉冲角度α2对射频发射场最不敏感的角度值。
7.一种扫描***,包括磁共振扫描设备和计算机,其中计算机包括处理器、存储介质及存储在存储介质上并可在处理器上运行的计算机程序,其特征在于,所述处理器执行所述程序时可用于执行一种磁共振脂肪抑制图像获取方法,所述方法包括以下步骤:获取单射频压脂脉冲的射频脉冲角度α;
通过单射频压脂脉冲的射频脉冲角度α确定双射频压脂脉冲的第一射频脉冲角度α1与第二射频脉冲角度α2;以及
基于所述双射频压脂脉冲的第一射频脉冲角度α1与第二射频脉冲角度α2获取脂肪抑制图像;
所述通过单射频压脂脉冲的射频脉冲角度α确定双射频压脂脉冲的第一射频脉冲角度α1与第二射频脉冲角度α2的步骤包括:
获取单射频压脂脉冲的射频脉冲角度α与双射频压脂脉冲的第一射频脉冲角度α1、第二射频脉冲角度α2之间的映射关系,以确定第一射频脉冲角度α1与第二射频脉冲角度α2的预选角度值;
所述获取单射频压脂脉冲的射频脉冲角度α与双射频压脂脉冲的第一射频脉冲角度α1、第二射频脉冲角度α2之间的映射关系的步骤包括:
通过以脂肪纵向磁矩的作用与单射频压脂脉冲作用相当,确定出单射频压脂脉冲的射频脉冲角度α与双射频压脂脉冲的第一射频脉冲角度α1、第二射频脉冲角度α2之间的映射关系:cos(α1)cos(α2)=cos(α)。
8.一种存储介质,其上存储有计算机程序,计算机可读取该存储介质内存储的计算机程序,其特征在于,该程序被处理器执行可用于执行权利要求1-6任意一项所述的方法。
CN201710588221.9A 2017-07-18 2017-07-18 磁共振脂肪抑制图像获取方法、存储介质及扫描*** Active CN107369153B (zh)

Priority Applications (4)

Application Number Priority Date Filing Date Title
CN201710588221.9A CN107369153B (zh) 2017-07-18 2017-07-18 磁共振脂肪抑制图像获取方法、存储介质及扫描***
US15/798,147 US10578699B2 (en) 2017-07-18 2017-10-30 Systems and methods for fat suppression in magnetic resonance imaging
PCT/CN2017/114102 WO2019015212A1 (en) 2017-07-18 2017-11-30 SYSTEMS AND METHODS FOR FAT SUPPRESSION IN MAGNETIC RESONANCE IMAGING
EP17204652.6A EP3432020A1 (en) 2017-07-18 2017-11-30 Systems and methods for fat suppression in magnetic resonance imaging

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710588221.9A CN107369153B (zh) 2017-07-18 2017-07-18 磁共振脂肪抑制图像获取方法、存储介质及扫描***

Publications (2)

Publication Number Publication Date
CN107369153A CN107369153A (zh) 2017-11-21
CN107369153B true CN107369153B (zh) 2020-11-27

Family

ID=60306700

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710588221.9A Active CN107369153B (zh) 2017-07-18 2017-07-18 磁共振脂肪抑制图像获取方法、存储介质及扫描***

Country Status (2)

Country Link
US (1) US10578699B2 (zh)
CN (1) CN107369153B (zh)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107369153B (zh) 2017-07-18 2020-11-27 上海联影医疗科技股份有限公司 磁共振脂肪抑制图像获取方法、存储介质及扫描***
WO2019015212A1 (en) * 2017-07-18 2019-01-24 Shenzhen United Imaging Healthcare Co., Ltd. SYSTEMS AND METHODS FOR FAT SUPPRESSION IN MAGNETIC RESONANCE IMAGING
US11287504B2 (en) * 2020-08-11 2022-03-29 Uih America, Inc. Systems and methods for magnetic resonance imaging
CN112315449B (zh) * 2020-11-24 2023-08-08 上海联影医疗科技股份有限公司 磁共振成像方法、***、电子装置和存储介质
CN116942134B (zh) * 2023-09-19 2023-11-28 山东奥新医疗科技有限公司 一种磁共振成像方法、装置、设备及存储介质

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1842720A (zh) * 2003-08-27 2006-10-04 皇家飞利浦电子股份有限公司 在磁共振成像中射频线圈的空间灵敏度分布的确定
CN101327123A (zh) * 2007-06-20 2008-12-24 通用电器横河医疗***株式会社 磁共振成像装置和磁共振图像产生方法
CN101396265A (zh) * 2007-09-28 2009-04-01 株式会社东芝 磁共振成像装置及其脂肪抑制效果的分析方法
CN101401723A (zh) * 2007-10-04 2009-04-08 株式会社东芝 Mri装置
CN102772206A (zh) * 2011-05-12 2012-11-14 上海联影医疗科技有限公司 磁共振成像方法
CN103202693A (zh) * 2012-01-13 2013-07-17 株式会社东芝 磁共振摄像装置
CN103561645A (zh) * 2011-06-09 2014-02-05 株式会社日立医疗器械 磁共振成像装置以及翻转角决定方法

Family Cites Families (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8054075B2 (en) * 2005-02-03 2011-11-08 The Johns Hopkins University Method for magnetic resonance imaging using inversion recovery with on-resonant water suppression including MRI systems and software embodying same
US8040134B2 (en) * 2006-10-06 2011-10-18 Hitachi Medical Corporation Magnetic resonance imaging device configured to suppress signals from fat by excluding effect of non-uniformity of irradiated magnetic field
JP5398149B2 (ja) * 2007-03-27 2014-01-29 株式会社東芝 磁気共鳴イメージング装置
US7990140B2 (en) * 2007-10-04 2011-08-02 Kabushiki Kaisha Toshiba MRI apparatus
US8723516B2 (en) * 2010-02-26 2014-05-13 Kabushiki Kaisha Toshiba B1-robust and T1-robust species suppression in MRI
US9008274B2 (en) * 2012-12-24 2015-04-14 General Electric Company Systems and methods for selecting image display parameters
WO2014121065A1 (en) * 2013-01-31 2014-08-07 The Board Of Trustees Of The University Of Illinois Two-dimensional j-resolved laser and semi-laser spectroscopy
DE102013221938B4 (de) 2013-10-29 2018-11-08 Siemens Healthcare Gmbh Verfahren zur Aufnahme von Magnetresonanzdaten mit einer diffusionsgewichteten Magnetresonanzsequenz und Magnetresonanzeinrichtung
US10052033B2 (en) * 2015-03-12 2018-08-21 Siemens Healthcare Gmbh System and method for motion-robust 3D magnetic resonance imaging of vessel walls
CN106466182B (zh) 2015-08-18 2020-03-31 上海联影医疗科技有限公司 磁共振成像方法及***
CN106952261B (zh) 2017-04-06 2019-11-22 深圳安科高技术股份有限公司 一种抑制脂肪信号的图像处理方法及装置
CN107369153B (zh) 2017-07-18 2020-11-27 上海联影医疗科技股份有限公司 磁共振脂肪抑制图像获取方法、存储介质及扫描***

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1842720A (zh) * 2003-08-27 2006-10-04 皇家飞利浦电子股份有限公司 在磁共振成像中射频线圈的空间灵敏度分布的确定
CN101327123A (zh) * 2007-06-20 2008-12-24 通用电器横河医疗***株式会社 磁共振成像装置和磁共振图像产生方法
CN101396265A (zh) * 2007-09-28 2009-04-01 株式会社东芝 磁共振成像装置及其脂肪抑制效果的分析方法
CN101401723A (zh) * 2007-10-04 2009-04-08 株式会社东芝 Mri装置
CN102772206A (zh) * 2011-05-12 2012-11-14 上海联影医疗科技有限公司 磁共振成像方法
CN103561645A (zh) * 2011-06-09 2014-02-05 株式会社日立医疗器械 磁共振成像装置以及翻转角决定方法
CN103202693A (zh) * 2012-01-13 2013-07-17 株式会社东芝 磁共振摄像装置

Also Published As

Publication number Publication date
CN107369153A (zh) 2017-11-21
US20190025395A1 (en) 2019-01-24
US10578699B2 (en) 2020-03-03

Similar Documents

Publication Publication Date Title
CN107369153B (zh) 磁共振脂肪抑制图像获取方法、存储介质及扫描***
US9081073B2 (en) System for suppression of artifacts in MR imaging
CN108209918B (zh) 磁共振成像方法及磁共振***
JP2011517983A (ja) 患者の安全性、及び走査性能の改良のためのリアルタイム性の局所及び大局のsar推定
KR20130018584A (ko) 자기 공명 노출 생성 방법
JP6356809B2 (ja) 水/脂肪分離を有するゼロエコー時間mrイメージング
JP2008514259A (ja) 磁気共鳴デバイス及び方法
US20160223634A1 (en) Magnetization transfer contrast technique for chemical exchange saturation transfer (cest) mri by localized steam and method of oeration thereof
US9846213B2 (en) Optimization of the noise development of a 3D gradient echo sequence in a magnetic resonance system
KR20130033473A (ko) 자기 공명 시스템을 제어하기 위한 방법
US7358731B2 (en) Determination of spatial sensitivity profiles of RF coils in magnetic resonance imaging
US20150309139A1 (en) Method and magnetic resonance apparatus for slice-selective magnetic resonance imaging
US10310043B2 (en) Method and apparatus for the correction of magnetic resonance scan data
US10909729B2 (en) Image processing apparatus
JP2013121437A (ja) 磁気共鳴イメージング装置及びT1ρイメージング法
US10393843B2 (en) Method and apparatus for accelerated acquisition of magnetic resonance data
US10254367B2 (en) Magnetic resonance imaging method and apparatus with motion-corrected model-based acceleration of parameter mapping
US10145926B2 (en) Phase-cycled steady-state free precession imaging
US20170231523A1 (en) Method and magnetic resonance apparatus scar quantification in the myocardium
EP1461632B1 (en) Parallel mr imaging method
US11547317B2 (en) TI scout for inversion recovery sequence
US9506998B2 (en) Method and magnetic resonance apparatus to generate raw data sets from double echo data acquisitions
US11698432B2 (en) Magnetic resonance imaging system, and main magnetic field correction method therefor and storage medium
US11650276B2 (en) Method for acquiring magnetic resonance (MR) data
JP7089964B2 (ja) 磁気共鳴イメージング装置

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
CB02 Change of applicant information
CB02 Change of applicant information

Address after: 201807 Shanghai City, north of the city of Jiading District Road No. 2258

Applicant after: Shanghai Lianying Medical Technology Co., Ltd

Address before: 201807 Shanghai City, north of the city of Jiading District Road No. 2258

Applicant before: SHANGHAI UNITED IMAGING HEALTHCARE Co.,Ltd.

GR01 Patent grant
GR01 Patent grant