CN114076912A - 抑制静态组织的方法、磁共振成像方法及系统 - Google Patents

抑制静态组织的方法、磁共振成像方法及系统 Download PDF

Info

Publication number
CN114076912A
CN114076912A CN202010838070.XA CN202010838070A CN114076912A CN 114076912 A CN114076912 A CN 114076912A CN 202010838070 A CN202010838070 A CN 202010838070A CN 114076912 A CN114076912 A CN 114076912A
Authority
CN
China
Prior art keywords
tissue
magnetic resonance
static
interest
region
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
CN202010838070.XA
Other languages
English (en)
Other versions
CN114076912B (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 CN202010838070.XA priority Critical patent/CN114076912B/zh
Priority to CN202311457609.7A priority patent/CN117518053A/zh
Publication of CN114076912A publication Critical patent/CN114076912A/zh
Application granted granted Critical
Publication of CN114076912B publication Critical patent/CN114076912B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/05—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves 
    • A61B5/055—Detecting, 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
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/0033—Features or image-related aspects of imaging apparatus classified in A61B5/00, e.g. for MRI, optical tomography or impedance tomography apparatus; arrangements of imaging apparatus in a room
    • A61B5/004—Features or image-related aspects of imaging apparatus classified in A61B5/00, e.g. for MRI, optical tomography or impedance tomography apparatus; arrangements of imaging apparatus in a room adapted for image acquisition of a particular organ or body part
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/0033—Features or image-related aspects of imaging apparatus classified in A61B5/00, e.g. for MRI, optical tomography or impedance tomography apparatus; arrangements of imaging apparatus in a room
    • A61B5/004—Features or image-related aspects of imaging apparatus classified in A61B5/00, e.g. for MRI, optical tomography or impedance tomography apparatus; arrangements of imaging apparatus in a room adapted for image acquisition of a particular organ or body part
    • A61B5/0042—Features or image-related aspects of imaging apparatus classified in A61B5/00, e.g. for MRI, optical tomography or impedance tomography apparatus; arrangements of imaging apparatus in a room adapted for image acquisition of a particular organ or body part for the brain
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/70—Means for positioning the patient in relation to the detecting, measuring or recording means
    • A61B5/704—Tables
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/74—Details of notification to user or communication with user or patient ; user input means
    • A61B5/7475—User input or interface means, e.g. keyboard, pointing device, joystick
    • A61B5/748—Selection of a region of interest, e.g. using a graphics tablet
    • 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/483—NMR imaging systems with selection of signals or spectra from particular regions of the volume, e.g. in vivo spectroscopy
    • G01R33/4838—NMR imaging systems with selection of signals or spectra from particular regions of the volume, e.g. in vivo spectroscopy using spatially selective suppression or saturation of MR signals
    • 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/5607—Image 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
    • 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/563—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution of moving material, e.g. flow contrast angiography
    • G01R33/56308—Characterization of motion or flow; Dynamic imaging

Landscapes

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

Abstract

本申请涉及一种抑制静态组织的方法、磁共振成像方法及系统。上述磁共振成像方法，包括将检测对象置于静态磁场中。确定所述检测对象的感兴趣区域。所述感兴趣区域包括多个片层，且每个片层中包括静态组织和流动组织。向所述检测对象发射射频脉冲，以同时激发所述感兴趣区域中静态组织和流动组织的核自旋。施加目标梯度场，进行空间编码，获取所述感兴趣区域的磁共振信号，其中，所述感兴趣区域的磁共振信号中所述静态组织对应的信号被部分抑制或者完全抑制。重建磁共振信号，以获取所述感兴趣区域的磁共振图像。上述方法中，通过修改目标梯度场，抑制成像区域的静态组织，实现对流动物体进行成像，而不需要通过额外的预处理脉冲抑制背景组织。

Description

抑制静态组织的方法、磁共振成像方法及系统
技术领域
本申请涉及成像领域，特别是涉及一种抑制静态组织的方法、磁共振成像方法及系统。
背景技术
磁共振成像技术(MRI)在临床诊断和科学研究中得到了越来越广泛的应用，具有安全，多对比度，以及对软组织有很好的分辨能力等优点。在非增强血管成像中，时间飞越法(Time Of Flight，TOF)，基于GRE序列，利用流入增强效应进行血管成像，最后通过最大强度投影(Maximum Intensity Projection，MIP)进行血管显示，这里背景信号会给MIP处理带来干扰，尤其是高亮的脂肪信号。
传统TOF序列时序中，选择方向和读出方向的梯度零阶矩带来的相位为0，一阶矩带来的相位为0，从而进行血管成像。传统TOF序列一般通过设置回波时间(Time Echo，TE)为反相位来抑制脂肪，这就限制了使用的回波时间(Time Echo，TE)，从而带来了长TE导致的流动信号损失。
发明内容
基于此，针对传统的脂肪抑制技术会导致流动信号损失问题，本申请提供一种抑制静态组织的方法、磁共振成像方法及系统。
一种磁共振成像方法，包括：
将检测对象置于静态磁场中；
确定所述检测对象的感兴趣区域，所述感兴趣区域包括多个片层，且每个片层中包括静态组织和流动组织；
向所述检测对象发射射频脉冲，以同时激发所述感兴趣区域中静态组织和流动组织的核自旋；
施加目标梯度场，进行空间编码，获取所述感兴趣区域的磁共振信号，其中，所述感兴趣区域的磁共振信号中所述静态组织对应的信号被部分抑制或者完全抑制；
重建磁共振信号，以获取所述感兴趣区域的磁共振图像。
在其中一个实施例中，所述目标梯度场包括选层方向的梯度和读出方向的梯度，所述目标梯度场通过如下步骤确定：
获取所述流动组织相对所述静态组织的相对速度和所述流动组织的初始位置；
根据所述相对速度和所述初始位置，设定选层方向的梯度和读出方向的梯度。
在其中一个实施例中，所述获取所述流动组织相对所述静态组织的相对速度和所述流动组织的初始位置的步骤包括：
施加流速编码梯度，以获取所述初始位置和在选层方向的相对速度以及在读出方向的相对速度。
其中一个实施例中，所述感兴趣区域为血管、脑室或者脊髓中的一种。
一种磁共振成像方法，包括：
将检测对象置于静态磁场中；
确定所述检测对象的感兴趣区域，所述感兴趣区域包括多个片层，且每个片层中包括静态组织和流动组织；
向所述检测对象施加成像序列，获取第一组磁共振信号，所述成像序列包括射频脉冲和目标梯度场，所述射频脉冲用于同时激发所述感兴趣区域中静态组织和流动组织的核自旋，所述目标梯度场用于对所述核自旋进行空间编码，以获取所述第一组磁共振信号，所述第一组磁共振信号中所述静态组织对应的信号被部分抑制或者完全抑制；
根据所述第一组磁共振信号，获取所述感兴趣区域的磁共振图像。
在其中一个实施例中，所述目标梯度场包括选层方向的梯度和读出方向的梯度，所述目标梯度场通过如下步骤确定：
获取所述流动组织相对所述静态组织的相对速度和所述流动组织的初始位置；
根据所述相对速度和所述初始位置，设定选层方向的梯度和读出方向的梯度。
在其中一个实施例中，所述第一组磁共振信号沿所述流动组织的流动方向获取。
一种抑制静态组织的方法，包括：
获取流动组织的初始位置、在选层方向所述流动组织相对静态组织的相对速度以及在读出方向所述流动组织相对静态组织的相对速度；
根据所述初始位置、所述在选层方向所述流动组织相对静态组织的相对速度以及所述在读出方向所述流动组织相对静态组织的相对速度，确定选层方向的梯度和读出方向的梯度，以确定施加于检测对象的目标梯度场；所述目标梯度场用于部分抑制或者完全抑制所述静态组织对应的磁共振信号。
在其中一个实施例中，所述获取流动组织的初始位置、在选层方向所述流动组织相对静态组织的相对速度以及在读出方向所述流动组织相对静态组织的相对速度的步骤包括：
发送位置采集指令，所述位置采集指令用于指示梯度线圈向所述检测对象施加流速编码梯度，以获取反馈信号；
根据所述反馈信号，获取所述初始位置、所述在选层方向所述流动组织相对静态组织的相对速度以及所述在读出方向所述流动组织相对静态组织的相对速度。
一种磁共振成像系统，包括：
扫描床，用于至少将检测对象的感兴趣区域置于扫描腔中，所述感兴趣区域包括多个片层，且每个片层中包括静态组织和流动组织；
射频线圈，用于向所述检测对象发射射频脉冲，以同时激发所述感兴趣区域中静态组织和流动组织的核自旋；
梯度线圈，用于施加目标梯度场，进行空间编码，获取所述感兴趣区域的磁共振信号，其中，所述感兴趣区域的磁共振信号中所述静态组织对应的信号被部分抑制或者完全抑制；
存储器，用于存储计算机程序；
处理器，所述处理器执行所述计算机程序时实现上述实施例中任一项所述的抑制静态组织的方法的步骤，或者上述实施例中任一项所述的磁共振成像方法的步骤。
上述磁共振成像方法，包括将检测对象置于静态磁场中。确定所述检测对象的感兴趣区域。所述感兴趣区域包括多个片层，且每个片层中包括静态组织和流动组织。向所述检测对象发射射频脉冲，以同时激发所述感兴趣区域中静态组织和流动组织的核自旋。施加目标梯度场，进行空间编码，获取所述感兴趣区域的磁共振信号，其中，所述感兴趣区域的磁共振信号中所述静态组织对应的信号被部分抑制或者完全抑制。重建磁共振信号，以获取所述感兴趣区域的磁共振图像。上述方法中，通过修改目标梯度场，抑制成像区域的静态组织，实现对流动物体进行成像，而不需要通过额外的预处理脉冲抑制背景组织。
附图说明
为了更清楚地说明本申请实施例或传统技术中的技术方案，下面将对实施例或传统技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本申请的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1为本申请一个实施例提供的磁共振成像方法的流程图；
图2为本申请一个实施例提供的2D时序图；
图3为本申请一个实施例提供的3D时序图；
图4为本申请一个实施例提供的成像结果对比图；
图5为本申请另一个实施例提供的抑制静态组织的方法的流程图。
具体实施方式
为使本申请的上述目的、特征和优点能够更加明显易懂，下面结合附图对本申请的具体实施方式做详细的说明。在下面的描述中阐述了很多具体细节以便于充分理解本申请。但是本申请能够以很多不同于在此描述的其它方式来实施，本领域技术人员可以在不违背本申请内涵的情况下做类似改进，因此本申请不受下面公开的具体实施的限制。
可以理解，本申请所使用的术语“第一”、“第二”等可在本文中用于描述各种元件，但这些元件不受这些术语限制。这些术语仅用于将第一元件与另一个元件区分。举例来说，在不脱离本申请的范围的情况下，可以将第一获取模块称为第二获取模块，且类似地，可将第二获取模块称为第一获取模块。第一获取模块和第二获取模块两者都是获取模块，但其不是同一个获取模块。
需要说明的是，当元件被称为“设置于”另一个元件，它可以直接在另一个元件上或者也可以存在居中的元件。当一个元件被认为是“连接”另一个元件，它可以是直接连接到另一个元件或者可能同时存在居中元件。
除非另有定义，本文所使用的所有的技术和科学术语与属于本申请的技术领域的技术人员通常理解的含义相同。本文中在本申请的说明书中所使用的术语只是为了描述具体的实施例的目的，不是旨在于限制本申请。本文所使用的术语“及/或”包括一个或多个相关的所列项目的任意的和所有的组合。
请参见图1，本申请提供一种磁共振成像方法。所述磁共振成像方法包括步骤S10-S50。
S10，将检测对象置于静态磁场中。
步骤S10中，将包含待进行磁共振成像的对象记为检测对象。所述检测对象可以是健康试者，也可以是患者，或者还可以是动物。扫描设备中的主磁体可以产生施加到检测对象的静态磁场。所述静态磁场也可以成为主磁场。所述主磁体还可以控制所述静态磁场的均匀性。
S20，确定所述检测对象的感兴趣区域。所述感兴趣区域包括多个片层，且每个片层中包括静态组织和流动组织。
步骤S20中，所述检测对象的感兴趣区域可以为是任何部位或组织，如心脏、血管或其他存在搏动区域的器官或者组织。可以通过计算机设定所述检测对象的感兴趣区域。每一个感兴趣区域可以为一个三维分块/三维体积。所述三维分块包括多个二维片层。所述静态组织可以为脂肪、肌肉、脑白质、脑灰质等任何静止的组织。所述流动组织可以为血液、脑脊液。在其中一个实施例中，所述感兴趣区域为血管、脑室或者脊髓中的一种。
S30，向所述检测对象发射射频脉冲，以同时激发所述感兴趣区域中静态组织和流动组织的核自旋。
步骤S30中，所述射频脉冲可以与所述静态磁场垂直。波形发生器可以产生所述射频脉冲序列。所述射频脉冲序列可以被射频功率放大器放大，由射频电子器件处理，并且被施加到射频发射线圈以产生响应于强大电流的第三磁场，强大电流由射频电子器件基于经放大的所述射频脉冲产生。
S40，施加目标梯度场，进行空间编码，获取所述感兴趣区域的磁共振信号，其中，所述感兴趣区域的磁共振信号中所述静态组织对应的信号被部分抑制或者完全抑制。
步骤S40中，在磁共振成像中，当所述检测对象俯卧或仰卧于扫描床上时，可以利用前后方向(即y方向)的梯度场对磁共振信号进行相位编码，利用左右方向(即x方向)的梯度场进行选层(或层选)编码，利用上下方向(即z方向)的梯度场进行频率编码/频率读出编码。
在三维磁共振成像中，首先需要选层，在层面内进行频率编码和相位编码，将磁共振信号分配到不同的像素位置上，从而形成磁共振图像。例如可以是利用左右方向(相对人体的左右方向)的梯度场进行层面选择编码，相位编码方向与层面选择编码方向垂直，即相位编码方向为前后方向/沿人体的正面指向侧面。当然，层面选择编码方向还为前后方向，与之对应地，相位编码方向则为左右方向。优选的，所述根据不同方向的梯度场分别对所述磁共振信号进行相位编码和频率编码，获取所述磁共振信号对应的编码数据的方法为，磁共振信号利用相对于肢体前后方向的梯度场对所述磁共振信号进行相位编码，利用沿肢体内血流的方向进行频率读出编码，获取所述磁共振信号对应的编码数据。
可以通过修改梯度场的序列实现所述感兴趣区域的磁共振信号中所述静态组织对应的信号被部分抑制或者完全抑制。梯度场的序列施加能够使得静态组织的磁共振信号被散相，流动组织的磁共振信号被聚相。在其中一个实施例中，通过设置序列中的选层和读出方向的梯度零阶矩带来的相位为π的奇数倍，序列中的选层和读出方向的梯度一阶矩带来的相位为0，实现对静态组织的抑制，对流动组织进行成像。该方法不需要额外的准备模块就可以完成静止组织的抑制。可选的，在TOF序列中，对2D序列设置选层和读出方向的梯度零阶矩带来的相位为π的奇数倍，一阶矩带来的相位为0；对3D序列设置读出方向的梯度零阶矩带来的相位为π的奇数倍，一阶矩带来的相位为0，从而修改后的2D和3D时序图如图2和图3所示。
S50，重建磁共振信号，以获取所述感兴趣区域的磁共振图像。
步骤S50中，将磁共振信号对应的编码数据填充入K空间中。对K空间进行重建处理，获取检测对象的感兴趣区域的磁共振图像。所述磁共振信号为梯度回波信号。
上述磁共振成像方法，包括将检测对象置于静态磁场中。确定所述检测对象的感兴趣区域。所述感兴趣区域包括多个片层，且每个片层中包括静态组织和流动组织。向所述检测对象发射射频脉冲，以同时激发所述感兴趣区域中静态组织和流动组织的核自旋。施加目标梯度场，以获取所述感兴趣区域的磁共振信号，其中，所述感兴趣区域的磁共振信号中所述静态组织对应的信号被部分抑制或者完全抑制。重建磁共振信号，以获取所述感兴趣区域的磁共振图像。上述方法中，通过修改目标梯度场，抑制成像区域的静态组织，实现对流动物体进行成像，而不需要通过额外的预处理脉冲抑制背景组织。
在其中一个实施例中，所述目标梯度场包括选层方向的梯度和读出方向的梯度，确定所述目标梯度场的方法为，获取所述流动组织相对所述静态组织的相对速度和所述流动组织的初始位置。根据所述相对速度和所述初始位置，设定选层方向的梯度和读出方向的梯度。
其中，相位与流速关系如下公式所示：
Figure BDA0002640419190000091
Figure BDA0002640419190000092
Figure BDA0002640419190000093
Figure BDA0002640419190000094
Figure BDA0002640419190000095
其中，公式1表示感兴趣区域的位置随时间变化的关系；公式2表示感兴趣区域的磁共振信号相位的表达式；公式3表示经变换后感兴趣区域的磁共振信号相位的表达式，公式中γ表示磁旋比，G表示梯度，x表示位置，υ表示运动速度，φ表示相位，mn表示n阶矩，t表示时间，μ表示梯度对应的时间轴方向的变量，且0≤μ≤t。
由上式可知，选层和读出方向的梯度零阶矩带来的相位为π的奇数倍，一阶矩带来的相位为0，即，
Figure BDA0002640419190000096
γm1(t)vx0＝0。所述流动组织相对所述静态组织的相对速度和所述流动组织的初始位置确定时，即可得到选层方向的梯度波形和读出方向的梯度波形。
在其中一个实施例中，所述获取所述流动组织相对所述静态组织的相对速度和所述流动组织的初始位置的方法可以为，施加流速编码梯度，以获取所述初始位置和在选层方向的相对速度以及在读出方向的相对速度。具体的，可以通过处理器发送位置采集指令。梯度线圈接收所述位置采集指令后，向所述检测对象施加流速编码梯度。接收线圈采集所述检测对象的反馈信号，并发送给处理器。通过所述处理器分析该反馈信号后即可获取所述流动组织相对所述静态组织的相对速度和所述流动组织的初始位置。
进而，根据所述初始位置和在选层方向的相对速度，可以得到所述选层方向的梯度需要满足：
Figure BDA0002640419190000101
γm1(t)vx0＝0 (公式5)
其中，γ为磁旋比，m0(t)为选层方向的梯度的零阶矩，x0为所述流动组织的初始位置，P为正整数；m1(t)为选层方向的梯度的一阶矩，vx0为所述流动组织在选层方向的相对速度。
根据所述初始位置和在读出方向的相对速度，可以得到所述读出方向的梯度需要满足：
Figure BDA0002640419190000102
γM1(t)vy0＝0 (公式6)
其中，γ为磁旋比，M0(t)为读出方向的梯度的零阶矩，x0为所述流动组织的初始位置，P为正整数；M1(t)为读出方向的梯度的一阶矩，vy0为所述流动组织在读出方向的相对速度。
本实施例中，根据上述方法，设置序列中的选层和读出方向的梯度零阶矩带来的相位为π的奇数倍，一阶矩带来的相位为0。请参见图4，为结果对比图，根据图4可知，本申请实现了静态组织、背景信号得到更好的抑制，而对流动组织(图中的血管区域)的成像没有影响。
请参见图5，本申请一个实施例中提供一种抑制静态组织的方法。所述抑制静态组织的方法包括：
S60，获取流动组织的初始位置、在选层方向所述流动组织相对静态组织的相对速度以及在读出方向所述流动组织相对静态组织的相对速度。
步骤S60中，相位与流速关系如下公式所示：
Figure BDA0002640419190000111
Figure BDA0002640419190000112
Figure BDA0002640419190000113
Figure BDA0002640419190000114
Figure BDA0002640419190000115
其中，相位与流速关系如下公式所示，其中，下述公式中γ表示磁旋比，G表示梯度，x表示位置，υ表示运动速度，φ表示相位，mn表示n阶矩。
由上式可知，选层和读出方向的梯度零阶矩带来的相位为π的奇数倍，一阶矩带来的相位为0，即，
Figure BDA0002640419190000116
γm1(t)vx0＝0。所述流动组织相对所述静态组织的相对速度和所述流动组织的初始位置确定时，即可得到选层方向的梯度波形和读出方向的梯度波形。
在其中一个实施例中，所述获取流动组织的初始位置、在选层方向所述流动组织相对静态组织的相对速度以及在读出方向所述流动组织相对静态组织的相对速度的方法可以为，发送位置采集指令。所述位置采集指令用于指示梯度线圈向所述检测对象施加流速编码梯度，以获取反馈信号。根据所述反馈信号，获取所述初始位置、所述在选层方向所述流动组织相对静态组织的相对速度以及所述在读出方向所述流动组织相对静态组织的相对速度。具体的，可以通过处理器发送位置采集指令。梯度线圈接收所述位置采集指令后，向所述检测对象施加流速编码梯度。接收线圈采集所述检测对象的反馈信号，并发送给处理器。通过所述处理器分析该反馈信号后即可获取所述流动组织相对所述静态组织的相对速度和所述流动组织的初始位置。
S70，根据所述初始位置、所述在选层方向所述流动组织相对静态组织的相对速度以及所述在读出方向所述流动组织相对静态组织的相对速度，确定选层方向的梯度和读出方向的梯度，以确定施加于检测对象的目标梯度场；所述目标梯度场用于部分抑制或者完全抑制所述静态组织对应的磁共振信号。
步骤S70中，根据所述初始位置和在选层方向的相对速度，可以得到所述选层方向的梯度需要满足：
Figure BDA0002640419190000121
γm1(t)vx0＝0 (公式11)
其中，γ为磁旋比，m0(t)为选层方向的梯度的零阶矩，x0为所述流动组织的初始位置，P为正整数；m1(t)为选层方向的梯度的一阶矩，vx0为所述流动组织在选层方向的相对速度。
根据所述初始位置和在读出方向的相对速度，可以得到所述读出方向的梯度需要满足：
γM0(t)x0＝(2p-1)π；
γM1(t)vy0＝0 (公式12)
其中，γ为磁旋比，M0(t)为读出方向的梯度的零阶矩，x0为所述流动组织的初始位置，P为正整数；M1(t)为读出方向的梯度的一阶矩，vy0为所述流动组织在读出方向的相对速度。
本实施例中，根据上述方法，设置序列中的选层和读出方向的梯度零阶矩带来的相位为π的奇数倍，一阶矩带来的相位为0。请参见图4，图中分别显示了使用现有血管成像序列采集的头部图像和采用本申请的成像方法采集得到的头部图像，图中的高亮区域为血液、流动组织对应的区域，灰度值较小的区域为静态组织、人体内的空腔等区域。根据图4可知，本申请实现了静态组织得到更好的抑制，而对流动组织的成像没有影响。
本申请一个实施例中提供一种磁共振成像系统。所述磁共振成像系统包括扫描床、射频线圈、梯度线圈、一个或多个处理器以及存储器。
所述扫描床用于至少将测对象的感兴趣区域置于扫描腔中，所述感兴趣区域包括多个片层，且每个片层中包括静态组织和流动组织。所述射频线圈用于向所述检测对象发射射频脉冲，以同时激发所述感兴趣区域中静态组织和流动组织的核自旋。所述梯度线圈用于施加目标梯度场，以获取所述感兴趣区域的磁共振信号，其中，所述感兴趣区域的磁共振信号中所述静态组织对应的信号被部分抑制或者完全抑制。所述存储器用于存储一个或多个程序，所述处理器执行程序时实现任一项所述的抑制静态组织的方法的步骤。
可以理解的是，将有待进行磁共振成像的对象记为检测对象。所述检测对象可以是健康试者，也可以是患者，或者还可以是动物。扫描设备中的主磁体可以产生施加到检测对象的静态磁场。所述静态磁场也可以成为主磁场。所述主磁体还可以控制所述静态磁场的均匀性。
所述检测对象的感兴趣区域可以为是任何部位或组织，如心脏、血管或其他存在搏动区域的器官或者组织。每一个感兴趣区域可以为一个三维分块。所述三维分块包括多个二维片层。所述静态组织可以为脂肪。所述流动组织可以为血液。在其中一个实施例中，所述感兴趣区域为血管、脑室或者脊髓中的一种。
所述射频脉冲可以与所述静态磁场垂直。波形发生器可以产生所述射频脉冲序列。所述射频脉冲序列可以被射频功率放大器放大，由射频电子器件处理，并且被施加到射频发射线圈以产生响应于强大电流的第三磁场，强大电流由射频电子器件基于经放大的所述射频脉冲产生。
在磁共振成像中，当所述检测对象俯卧或仰卧于扫描床上时，可以利用前后方向(即y方向)的梯度场对磁共振信号进行相位编码(phase encoding，PE)，利用左右方向(即x方向)的梯度场进行选层(或层选)编码(s1ice phase encoding，SPE)，利用上下方向(即z方向)的梯度场进行频率编码/频率读出编码(readout encoding，RE)。
在三维磁共振成像中，首先需要选层，在层面内进行频率编码和相位编码，将磁共振信号分配到不同的像素位置上，从而形成磁共振图像。例如可以是利用左右方向(相对人体的左右方向)的梯度场进行层面选择编码，相位编码方向与层面选择编码方向垂直，即相位编码方向为前后方向/沿人体的正面指向侧面。当然，层面选择编码方向还为前后方向，与之对应地，相位编码方向则为左右方向。优选的，所述根据不同方向的梯度场分别对所述磁共振信号进行相位编码和频率编码，获取所述磁共振信号对应的编码数据的方法为，磁共振信号利用相对于肢体前后方向的梯度场对所述磁共振信号进行相位编码，利用沿肢体内血流的方向进行频率读出编码，获取所述磁共振信号对应的编码数据。
可以通过修改梯度场的序列实现所述感兴趣区域的磁共振信号中所述静态组织对应的信号被部分抑制或者完全抑制。在其中一个实施例中，通过设置序列中的选层和读出方向的梯度零阶矩带来的相位为π的奇数倍，序列中的选层和读出方向的梯度一阶矩带来的相位为0，实现对静态组织的抑制，对流动组织进行成像。该方法不需要额外的准备模块就可以完成静止组织的抑制。可选的，在TOF序列中，对2D序列设置选层和读出方向的梯度零阶矩带来的相位为π的奇数倍，一阶矩带来的相位为0；对3D序列设置读出方向的梯度零阶矩带来的相位为π的奇数倍，一阶矩带来的相位为0。
如图2为本申请一实施例的TOF 2D序列示意图，其中，RF表示射频线圈发射的射频脉冲；Gss表示选层方向的梯度场；Gpe表示相位编码方向的梯度场；Gro表示读出方向的梯度场。在此实施例中，TOF 2D序列具体为在RF脉冲210施加的同时在Gss方向施加层面选择梯度220，且在220之后施加沿Gss方向的失相梯度230；随后分别沿着Gpe和Gro方向先后施加失相梯度240、270，并且在失相梯度270施加的同时沿Gpe方向施加相位编码梯度250，在失相梯度270之后紧随施加频率编码梯度280以采集得到梯度回波信号290。当然为了减小梯度回波信号对下一射频脉冲激发信号的影响，梯度回波信号290采集后还可施加沿Gpe方向的失相梯度260。在此实施例中，将选层和读出方向的梯度零阶矩带来的相位为π的奇数倍，将选层和读出方向的梯度一阶矩带来的相位为0，能够使得感兴趣区域的磁共振信号中所述静态组织对应的信号被部分抑制或者完全抑制。本实施例中，沿Gpe方向施加的失相梯度240、260内部所列的相同箭头表示两者属于同类型的梯度，相位编码梯度250的箭头与失相梯度240、260相反，表示与失相梯度240、260属于不同的类型。
如图3为本申请一实施例的TOF 3D序列示意图。在此实施例中，TOF 3D序列具体为在RF脉冲310施加的同时在Gss方向施加层面选择梯度320，且在320之后先后施加沿Gss方向的编码梯度330、失相梯度340；随后分别沿着Gpe和Gro方向先后施加失相梯度350、370，且在350之后紧随施加沿Gpe方向施加编码梯度360，在370之后紧随施加频率编码梯度380以采集得到梯度回波信号390。当然为了减小梯度回波信号对下一射频脉冲激发信号的影响，梯度回波信号390采集后还可分别施加沿Gpe方向、Gss方向施加的失相梯度341、351。在此实施例中，设置沿Gro方向的梯度零阶矩带来的相位为PI的奇数倍，设置沿Gro方向的梯度一阶矩带来的相位为0，能够使得感兴趣区域的磁共振信号中所述静态组织对应的信号被部分抑制或者完全抑制。在此实施例中使用三维TOF成像能够获得空间分辨力刚搞的图像，而且由于体素较小，流动失相位相对较轻，受湍流的影响相对较小。
在一个实施例中，检测对象选择血管，流动组织为血流，该磁共振成像方法包括：首先，采用TOF 3D序列对检测对象的3D容积执行顺血流采集得到第一组磁共振信号；然后采用TOF 3D序列对检测对象的3D容积执行逆血流采集得到第二组磁共振信号；重建第一组磁共振信号和第二组磁共振信号获得检测对象得磁共振图像。在此实施例中，顺血流采集具体为先采集血流近端层面磁共振信号，然后向血流远侧逐层采集，且3D容积的分层方向与血流方向一致。逆血流采集具体为先采集血流远端层面磁共振信号，然后向血流近侧逐层采集，且3D容积的分层方向与血流方向相反。本实施例中，通过两次血流信号得采集，可以有效减少血流饱和效应。
将磁共振信号对应的编码数据填充入K空间中。对K空间进行重建处理，获取3D容积的磁共振图像。所述磁共振信号可以是梯度回波信号。
在一个实施例中，检测对象选择血管，流动组织为血流，磁共振成像方法选择如图3所示的TOF 3D成像序列对检测对象执行扫描，如前述内容，成像序列包括射频脉冲和目标梯度场，射频脉冲用于同时激发3D容积中静态组织和流动组织的核自旋，目标梯度场用于对核自旋进行空间编码，以获取磁共振信号，第一组磁共振信号中所述静态组织对应的信号被部分抑制或者完全抑制。进一步的，TOF 3D成像序列通过调节沿Gss方向的编码梯度330、沿Gpe方向的编码梯度360的施加时序，使得在磁共振信号的采集过程中沿Gss方向以连续的方式采集，在层面内Gpe方向以隔行扫描的方式采集。本实施例中通过隔行采集方式有利于减少血流饱和效应，使得整个3D容积的血流信号强度均一化，有利于减小血管内信号强度的波动，并能够显示慢血流和小血管的显示。此外，本实施例的技术方案还改变了对血流速度和方向的敏感性，改善了对血管狭窄和异常血管的显示率。
存储器作为一种计算机可读存储介质，可用于存储软件程序、计算机可执行程序以及模块，如本申请实施例中的磁共振成像方法对应的程序指令/模块。处理器通过运行存储在存储器中的软件程序、指令以及模块，从而执行设备的各种功能应用以及数据处理，即实现上述磁共振成像方法。
存储器可主要包括存储程序区和存储数据区，其中，存储程序区可存储操作系统、至少一个功能所需的应用程序。存储数据区可存储根据终端的使用所创建的数据等。此外，存储器可以包括高速随机存取存储器，还可以包括非易失性存储器，例如至少一个磁盘存储器件、闪存器件、或其他非易失性固态存储器件。在一些实例中，存储器可进一步包括相对于处理器远程设置的存储器，这些远程存储器可以通过网络连接至设备。上述网络的实例包括但不限于互联网、企业内部网、局域网、移动通信网及其组合。
上述磁共振成像系统，包括将检测对象置于静态磁场中。确定所述检测对象的感兴趣区域。所述感兴趣区域包括多个片层，且每个片层中包括静态组织和流动组织。向所述检测对象发射射频脉冲，以同时激发所述感兴趣区域中静态组织和流动组织的核自旋。施加目标梯度场，以获取所述感兴趣区域的磁共振信号，其中，所述感兴趣区域的磁共振信号中所述静态组织对应的信号被部分抑制或者完全抑制。重建磁共振信号，以获取所述感兴趣区域的磁共振图像。上述方法中，通过修改目标梯度场，抑制成像区域的静态组织，实现对流动物体进行成像，而不需要通过额外的预处理脉冲抑制背景组织。
以上所述实施例的各技术特征可以进行任意的组合，为使描述简洁，未对上述实施例中的各个技术特征所有可能的组合都进行描述，然而，只要这些技术特征的组合不存在矛盾，都应当认为是本说明书记载的范围。
以上所述实施例仅表达了本申请的几种实施方式，其描述较为具体和详细，但并不能因此而理解为对申请专利范围的限制。应当指出的是，对于本领域的普通技术人员来说，在不脱离本申请构思的前提下，还可以做出若干变形和改进，这些都属于本申请的保护范围。因此，本申请专利的保护范围应以所附权利要求为准。

Claims (10)

1.一种磁共振成像方法，其特征在于，包括：
将检测对象置于静态磁场中；
确定所述检测对象的感兴趣区域，所述感兴趣区域包括多个片层，且每个片层中包括静态组织和流动组织；
向所述检测对象发射射频脉冲，以同时激发所述感兴趣区域中静态组织和流动组织的核自旋；
施加目标梯度场，进行空间编码，获取所述感兴趣区域的磁共振信号，其中，所述感兴趣区域的磁共振信号中所述静态组织对应的信号被部分抑制或者完全抑制；
重建磁共振信号，以获取所述感兴趣区域的磁共振图像。
2.根据权利要求1所述的磁共振成像方法，其特征在于，所述目标梯度场包括选层方向的梯度和读出方向的梯度，所述目标梯度场通过如下步骤确定：
获取所述流动组织相对所述静态组织的相对速度和所述流动组织的初始位置；
根据所述相对速度和所述初始位置，设定选层方向的梯度和读出方向的梯度。
3.根据权利要求2所述的磁共振成像方法，其特征在于，所述获取所述流动组织相对所述静态组织的相对速度和所述流动组织的初始位置的步骤包括：
施加流速编码梯度，以获取所述初始位置和在选层方向的相对速度以及在读出方向的相对速度。
4.根据权利要求1所述的磁共振成像方法，其特征在于，所述感兴趣区域为血管、脑室或者脊髓中的一种。
5.一种磁共振成像方法，其特征在于，包括：
将检测对象置于静态磁场中；
确定所述检测对象的感兴趣区域，所述感兴趣区域包括多个片层，且每个片层中包括静态组织和流动组织；
向所述检测对象施加成像序列，获取第一组磁共振信号，所述成像序列包括射频脉冲和目标梯度场，所述射频脉冲用于同时激发所述感兴趣区域中静态组织和流动组织的核自旋，所述目标梯度场用于对所述核自旋进行空间编码，以获取所述第一组磁共振信号，所述第一组磁共振信号中所述静态组织对应的信号被部分抑制或者完全抑制；
根据所述第一组磁共振信号，获取所述感兴趣区域的磁共振图像。
6.根据权利要求5所述的磁共振成像方法，其特征在于，所述目标梯度场包括选层方向的梯度和读出方向的梯度，所述目标梯度场通过如下步骤确定：
获取所述流动组织相对所述静态组织的相对速度和所述流动组织的初始位置；
根据所述相对速度和所述初始位置，设定选层方向的梯度和读出方向的梯度。
7.根据权利要求5所述的磁共振成像方法，其特征在于，所述第一组磁共振信号沿所述流动组织的流动方向获取。
8.一种抑制静态组织的方法，其特征在于，包括：
获取流动组织的初始位置、在选层方向所述流动组织相对静态组织的相对速度以及在读出方向所述流动组织相对静态组织的相对速度；
根据所述初始位置、所述在选层方向所述流动组织相对静态组织的相对速度以及所述在读出方向所述流动组织相对静态组织的相对速度，确定选层方向的梯度和读出方向的梯度，以确定施加于检测对象的目标梯度场；所述目标梯度场用于部分抑制或者完全抑制所述静态组织对应的磁共振信号。
9.根据权利要求8所述的抑制静态组织的方法，其特征在于，所述获取流动组织的初始位置、在选层方向所述流动组织相对静态组织的相对速度以及在读出方向所述流动组织相对静态组织的相对速度的步骤包括：
发送位置采集指令，所述位置采集指令用于指示梯度线圈向所述检测对象施加流速编码梯度，以获取反馈信号；
根据所述反馈信号，获取所述初始位置、所述在选层方向所述流动组织相对静态组织的相对速度以及所述在读出方向所述流动组织相对静态组织的相对速度。
10.一种磁共振成像系统，其特征在于，包括：
扫描床，用于至少将检测对象的感兴趣区域置于扫描腔中，所述感兴趣区域包括多个片层，且每个片层中包括静态组织和流动组织；
射频线圈，用于向所述检测对象发射射频脉冲，以同时激发所述感兴趣区域中静态组织和流动组织的核自旋；
梯度线圈，用于施加目标梯度场，进行空间编码，获取所述感兴趣区域的磁共振信号，其中，所述感兴趣区域的磁共振信号中所述静态组织对应的信号被部分抑制或者完全抑制；
存储器，用于存储计算机程序；
处理器，所述处理器执行所述计算机程序时实现权利要求1至7任一项所述的磁共振成像方法的步骤，或者利要求8或9所述的抑制静态组织的方法的步骤。
CN202010838070.XA 2020-08-19 2020-08-19 抑制静态组织的方法、磁共振成像方法及系统 Active CN114076912B (zh)

Priority Applications (2)

Application Number Priority Date Filing Date Title
CN202010838070.XA CN114076912B (zh) 2020-08-19 2020-08-19 抑制静态组织的方法、磁共振成像方法及系统
CN202311457609.7A CN117518053A (zh) 2020-08-19 2020-08-19 磁共振成像方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010838070.XA CN114076912B (zh) 2020-08-19 2020-08-19 抑制静态组织的方法、磁共振成像方法及系统

Related Child Applications (1)

Application Number Title Priority Date Filing Date
CN202311457609.7A Division CN117518053A (zh) 2020-08-19 2020-08-19 磁共振成像方法

Publications (2)

Publication Number Publication Date
CN114076912A true CN114076912A (zh) 2022-02-22
CN114076912B CN114076912B (zh) 2023-09-22

Family

ID=80281639

Family Applications (2)

Application Number Title Priority Date Filing Date
CN202311457609.7A Pending CN117518053A (zh) 2020-08-19 2020-08-19 磁共振成像方法
CN202010838070.XA Active CN114076912B (zh) 2020-08-19 2020-08-19 抑制静态组织的方法、磁共振成像方法及系统

Family Applications Before (1)

Application Number Title Priority Date Filing Date
CN202311457609.7A Pending CN117518053A (zh) 2020-08-19 2020-08-19 磁共振成像方法

Country Status (1)

Country Link
CN (2) CN117518053A (zh)

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5386190A (en) * 1991-07-31 1995-01-31 Takeuchi; Hiroyuki Selective excitation method of nuclear spins in magnetic resonance imaging and MRI apparatus using the same
CN101023866A (zh) * 2005-12-22 2007-08-29 Ge医疗系统环球技术有限公司 磁共振成像设备
CN101129262A (zh) * 2006-08-21 2008-02-27 美国西门子医疗解决公司 在磁共振成像中的对比度增强
CN101259018A (zh) * 2006-12-07 2008-09-10 Ge医疗系统环球技术有限公司 磁共振成像设备和磁共振成像方法
EP2148212A2 (en) * 2008-07-24 2010-01-27 Kabushiki Kaisha Toshiba Magnetic resonance imaging apparatus for contrast enhancement of images

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5386190A (en) * 1991-07-31 1995-01-31 Takeuchi; Hiroyuki Selective excitation method of nuclear spins in magnetic resonance imaging and MRI apparatus using the same
CN101023866A (zh) * 2005-12-22 2007-08-29 Ge医疗系统环球技术有限公司 磁共振成像设备
CN101129262A (zh) * 2006-08-21 2008-02-27 美国西门子医疗解决公司 在磁共振成像中的对比度增强
CN101259018A (zh) * 2006-12-07 2008-09-10 Ge医疗系统环球技术有限公司 磁共振成像设备和磁共振成像方法
EP2148212A2 (en) * 2008-07-24 2010-01-27 Kabushiki Kaisha Toshiba Magnetic resonance imaging apparatus for contrast enhancement of images

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
胡丽丽: "抑制MRI运动伪影的双脉冲预饱和技术的实现" *

Also Published As

Publication number Publication date
CN114076912B (zh) 2023-09-22
CN117518053A (zh) 2024-02-06

Similar Documents

Publication Publication Date Title
JP6084573B2 (ja) マルチポイントディクソン技術を用いるｍｒ撮像
CN107072592B (zh) 磁共振成像装置以及定量性磁化率匹配方法
US5771893A (en) Method and apparatus for nuclear magnetic resonance imaging of physiological function information
US9746539B2 (en) MR imaging with suppresion of flow artifacts
JP5127841B2 (ja) 磁気共鳴イメージング装置及び磁化率強調画像撮影方法
JP5599893B2 (ja) ナビゲータを使用するｍｒイメージング
KR101629165B1 (ko) 자기공명영상장치 및 그 제어방법
US10203387B2 (en) MR imaging with enhanced susceptibility contrast
RU2739479C2 (ru) Магнитно-резонансная томография с разделением воды и жира по методу диксона
KR20130018584A (ko) 자기 공명 노출 생성 방법
WO2015190508A1 (ja) 磁気共鳴イメージング装置及び水脂肪分離画像作成方法
JP3556176B2 (ja) Ｍｒ投影画像の自動修正
US10132902B2 (en) Intrinsic navigation from velocity-encoding gradients in phase-contrast MRI
US11033199B2 (en) Echo-planar imaging magnetic resonance elastography pulse sequence
JP2017153965A (ja) 磁気共鳴イメージング装置及び磁気共鳴イメージング方法
US8143891B2 (en) System for image acquisition with fast magnetic resonance gradient echo sequences
CN109716155A (zh) 具有迪克逊型水/脂肪分离的mr成像
JP6202761B2 (ja) 磁気共鳴イメージング装置及びその処理方法
US10444316B2 (en) Reduction of eddy currents during flow encoded magnetic resonance imaging
KR101629162B1 (ko) 자기공명영상장치 및 그 제어방법
CN114076912B (zh) 抑制静态组织的方法、磁共振成像方法及系统
CN113466768B (zh) 磁共振成像方法及磁共振成像系统
US10782375B2 (en) Multi-contrast images from a magnetic resonance imaging scan
EP3859366A1 (en) Mr imaging using dixon-type water/fat separation with suppression of flow-induced leakage and/or swapping artifacts

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