CN110215210B - 磁共振成像方法、系统、计算机设备和可读存储介质 - Google Patents

磁共振成像方法、系统、计算机设备和可读存储介质 Download PDF

Info

Publication number
CN110215210B
CN110215210B CN201910452080.7A CN201910452080A CN110215210B CN 110215210 B CN110215210 B CN 110215210B CN 201910452080 A CN201910452080 A CN 201910452080A CN 110215210 B CN110215210 B CN 110215210B
Authority
CN
China
Prior art keywords
radio frequency
pulse
magnetic resonance
bandwidth
magnetic field
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
CN201910452080.7A
Other languages
English (en)
Other versions
CN110215210A (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 CN201910452080.7A priority Critical patent/CN110215210B/zh
Publication of CN110215210A publication Critical patent/CN110215210A/zh
Application granted granted Critical
Publication of CN110215210B publication Critical patent/CN110215210B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • 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/7203Signal processing specially adapted for physiological signals or for diagnostic purposes for noise prevention, reduction or removal
    • A61B5/7217Signal processing specially adapted for physiological signals or for diagnostic purposes for noise prevention, reduction or removal of noise originating from a therapeutic or surgical apparatus, e.g. from a pacemaker
    • 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/7235Details of waveform analysis
    • A61B5/7253Details of waveform analysis characterised by using transforms
    • A61B5/7257Details of waveform analysis characterised by using transforms using Fourier transforms
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02ATECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/30Assessment of water resources

Landscapes

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

Abstract

本发明涉及一种磁共振成像方法、系统、计算机设备和可读存储介质,该方法包括:控制磁体产生主磁场;检测对象位于主磁场形成的视场区域内,检测对象携带金属植入物,且金属植入物使主磁场非均匀;控制射频发射线圈发射射频脉冲;射频脉冲的整体带宽覆盖范围是由金属植入物引起的场不均匀性程度确定的;采用预设的读出带宽控制射频接收线圈接收回波信号,将回波信号对应的回波数据保存在K空间中,对K空间中的回波数据进行傅里叶变换,得到磁共振图像。该方法能够使得发射的射频脉冲覆盖不同的进动频率范围,最终使得成像区域的质子都能够被有效激发,并能够更加准确的反映扫描体的有效信息,减少因植入物带来的场不均匀性造成的图像伪影。

Description

磁共振成像方法、系统、计算机设备和可读存储介质
技术领域
本发明涉及磁共振成像领域,特别是涉及一种磁共振成像方法、系统、计算机设备和可读存储介质。
背景技术
磁共振成像技术中,氢原子核在磁体提供的主磁场中由杂乱无序的热平衡状态转为部分顺主磁场方向、部分逆主磁场方向,二者之差形成净磁化矢量,射频发射线圈将氢原子核由主磁场方向翻转到横向平面,并绕主磁场进动,射频接收线圈感应出电流信号,经信号处理和图像重建单元得到被成像的组织的图像,而为了能够有效定位来自不同空间位置的信号,会施加额外的选层梯度、相位编码梯度和读出梯度,建立起接收信号频率与空间位置之间的线性对应关系。但当主磁场不均匀时,该对应关系不再成立,会造成选层不准、图像混叠错位等问题。
鉴于此,有必要提出一种改进型的能够对磁场均匀性不敏感的磁共振成像方法。
发明内容
基于此,有必要针对当主磁场不均匀时,接收信号频率与空间位置之间的线性对应关系不再成立,造成选层不准、图像混叠错位等问题,提供一种磁共振成像方法、系统、计算机设备和可读存储介质。
第一方面,本发明实施例提供一种磁共振成像方法,所述方法包括:
控制磁体产生主磁场;检测对象位于所述主磁场形成的视场区域内,所述检测对象携带金属植入物,且所述金属植入物使所述主磁场非均匀;
控制射频发射线圈发射射频脉冲;所述射频脉冲的整体带宽覆盖范围是由所述金属植入物引起的场不均匀性程度确定的;
采用预设的读出带宽控制射频接收线圈接收回波信号,将所述回波信号对应的回波数据保存在K空间中,对所述K空间中的回波数据进行傅里叶变换,得到磁共振图像。
在其中一个实施例中,所述射频脉冲包括一个中心频率,相应的,所述回波信号对应一个数据集。
在其中一个实施例中,所述射频脉冲包括至少两个中心频率,相应的,所述回波信号对应至少两个数据集,每个中心频率对应一个数据集。
在其中一个实施例中,所述射频脉冲包括激发射频脉冲和回聚射频脉冲,所述激发射频脉冲包括多个中心频率,所述回聚射频脉冲包括一个中心频率,且所述激发射频脉冲的各中心频率对应不同的频带,各频带整体范围覆盖所述回聚射频脉冲所对应的带宽范围;
或者,所述回聚射频脉冲包括多个中心频率,所述激发射频脉冲包括一个中心频率,且所述回聚射频脉冲的各中心频率对应不同的频带,各频带整体范围能够覆盖所述激发射频脉冲所对应的带宽范围。
在其中一个实施例中,所述至少两个数据集分别填充至至少两个K空间,对所述K空间回波数据进行傅里叶变换,得到磁共振图像,包括:
分别对所述各K空间进行傅里叶变换,得到至少两个三维子图像;
对所述各三维子图像进行复合,获取所述磁共振图像。
在其中一个实施例中,所述方法还包括:
控制所述射频发射线圈多次发射射频脉冲,并采用所述预设的读出带宽控制所述射频接收线圈多次接收回波信号,得到多个数据集。
在其中一个实施例中,所述预设的读出带宽大于预设的带宽阈值。
第二方面,本发明实施例提供一种磁共振成像系统,所述系统包括:磁体、射频发射线圈、射频接收线圈和控制系统;
所述磁体,用于产生主磁场,并形成容纳检测对象的孔腔,且所述检测对象携带金属植入物;
所述射频发射线圈,用于发射射频脉冲;
所述射频接收线圈,用于接收磁共振信号;
所述控制系统,用于执行上述实施例所述的磁共振成像方法。
第三方面,本发明实施例提供一种计算机设备,包括存储器和处理器,所述存储器存储有计算机程序,所述处理器执行所述计算机程序时实现以下步骤:
控制磁体产生主磁场;检测对象位于所述主磁场形成的视场区域内,所述检测对象携带金属植入物,且所述金属植入物使所述主磁场非均匀;
控制射频发射线圈发射射频脉冲;所述射频脉冲的整体带宽覆盖范围是由所述金属植入物引起的场不均匀性程度确定的;
采用预设的读出带宽控制射频接收线圈接收回波信号,将所述回波信号对应的回波数据保存在K空间中,对所述K空间中的回波数据进行傅里叶变换,得到磁共振图像。
第四方面,本发明实施例提供一种计算机可读存储介质,其上存储有计算机程序,所述计算机程序被处理器执行时实现以下步骤:
控制磁体产生主磁场;检测对象位于所述主磁场形成的视场区域内,所述检测对象携带金属植入物,且所述金属植入物使所述主磁场非均匀;
控制射频发射线圈发射射频脉冲;所述射频脉冲的整体带宽覆盖范围是由所述金属植入物引起的场不均匀性程度确定的;
采用预设的读出带宽控制射频接收线圈接收回波信号,将所述回波信号对应的回波数据保存在K空间中,对所述K空间中的回波数据进行傅里叶变换,得到磁共振图像。
上述实施例提供的磁共振成像方法、系统、计算机设备和可读存储介质中,射频发射线圈发射的射频脉冲的整体带宽覆盖范围是由检测对象携带的金属植入物引起的场不均匀性程度确定的,当主磁场不均匀时,会使得质子的进动频率偏置分布,由于发射的射频脉冲的整体带宽范围足够大,能够使得发射的射频脉冲覆盖不同的进动频率范围,最终使得成像区域的质子都能够被有效激发,并采用预设的读出带宽接收相应的回波信号,并保存在对应的K空间中,进而对k空间中的回波数据进行傅里叶变换得到的磁共振图像,能够更加准确的反映扫描体的有效信息,减少因植入物带来的场不均匀性造成的图像伪影。
附图说明
图1为一个实施例提供的磁共振成像方法的应用环境示意图;
图2为一个实施例提供的磁共振成像方法的流程示意图;
图3为一个实施例提供的发射的射频脉冲的示意图;
图4为另一个实施例提供的发射的射频脉冲的示意图;
图5为另一个实施例提供的磁共振成像方法的流程示意图;
图6为一个实施例提供的磁共振成像系统结构示意图;
图7为一个实施例提供的计算机设备的内部结构示意图。
具体实施方式
为了使本申请的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本申请进行进一步详细说明。应当理解,此处描述的具体实施例仅仅用以解释本申请,并不用于限定本申请。
本申请实施例提供的磁共振成像方法,可以适用于如图1所示的应用环境中。人体中氢原子核自旋,可等效为一个小磁针,在主磁体提供的强磁场中,氢原子核由杂乱无序的热平衡状态转为部分顺、部分逆主磁场方向,二者的差形成净磁化矢量,射频发射线圈将氢原子核由主磁场方向翻转到横向平面,并绕主磁场进动,梯度线圈发射梯度脉冲定位来自不同位置的信号,射频接收线圈接收磁共振信号,计算机设备根据接收到的磁共振信号得到被成像的组织的图像。
下面以具体的实施例对本发明的技术方案以及本发明的技术方案如何解决上述技术问题进行详细说明。下面这几个具体的实施例可以相互结合,对于相同或相似的概念或过程可能在某些实施例中不再赘述。
图2为一个实施例提供的磁共振成像方法的流程示意图。本实施例涉及的是磁共振系统发射射频脉冲,且发射的射频脉冲的整体带宽覆盖范围是由检测对象携带的金属植入物引起的场不均匀性程度确定的,并采用预设的读出带宽接收回波信号,将回波信号对应的回波数据保存在K空间中,对K空间中的回波数据进行傅里叶变换,得到磁共振图像的具体实现过程。如图2所示,该方法可以包括:
S201,控制磁体产生主磁场;检测对象位于主磁场形成的视场区域内,检测对象携带金属植入物,且金属植入物使主磁场非均匀。
具体的,磁体可环绕形成具有容纳空间的孔腔,检测对象可被已送至孔腔内,计算机设备控制磁体产生主磁场,检测对象位于主磁场形成的视场区域内,且检测对象携带有金属植入物,该金属植入物使主磁场非均匀。其中,主磁场的场强越高,检测对象携带的金属植入物会使得主磁场的非均匀性越严重。
S202,控制射频发射线圈发射射频脉冲;射频脉冲的整体带宽覆盖范围是由金属植入物引起的场不均匀性程度确定的。
具体的,计算机设备控制射频发射线圈发射射频脉冲,其中,发射的射频脉冲的整体带宽覆盖范围是由上述金属植入物引起的场不均匀性程度确定的。可选的,发射的射频脉冲可以包括激发射频脉冲和回聚射频脉冲,发射的激发射频脉冲和回聚射频脉冲中至少一种脉冲的带宽大于预设的第一带宽阈值。可选的,预设的第一带宽阈值可以为3KHz。需要说明的是,本实施例发射的射频激发脉冲的扫描序列为快速自旋回波序列(Fast SpinEcho,FSE),FSE序列是一种常用的磁共振扫描序列,FSE序列在激发脉冲作用的T时间后,会在与激发脉冲相差90度相位的x-y平面施加回聚脉冲,在回聚脉冲的180度分量的作用下,所有的质子会产生相位翻转,在x-y平面内翻转180度,而其进动频率和方向不发生改变,在经过相同的时间T后,各质子自旋将会再次位于同相位。可选的,射频发射线圈可以为容积线圈,其可可包括鸟笼线圈、横电磁线圈、表面线圈、鞍形线圈等,也可以为局部线圈,其可包括鸟笼线圈、螺线管线圈、鞍形线圈、柔性线圈等。
S203,采用预设的读出带宽控制射频接收线圈接收回波信号,将回波信号对应的回波数据保存在K空间中,对K空间中的回波数据进行傅里叶变换,得到磁共振图像。
具体的,计算机设备采用预设的读出带宽控制射频接收线圈接收回波信号,将回波信号对应的回波数据保存在K空间中,对K空间中的回波数据进行傅里叶变换,得到磁共振图像。可选地,在接收回波信号之前可控制梯度线圈产生空间梯度场,该梯度场可包括选层梯度、相位编码梯度和读出梯度,通过该梯度场可建立起接收信号频率与空间位置之间的线性对应关系。在一些实施例中,所捕获到的回波信号可根据采样算法来被采样以生成要被填充到k空间中的一组离散数据。可选的,可以有各种用于以所捕获到的信号来填充k空间的采样算法,包括笛卡尔采样(逐行)、径向采样、螺旋采样、回旋形采样等。其中,预设的读出带宽根据扫描序列的模数转换(Analog to Digital Converter,ADC)采集窗口确定,回波信号为受检者对射频脉冲响应后生成的携带受检者相关信息的信号。可选的,采用的预设的读出带宽的带宽可以为大于500HZ的带宽。需要说明的是,K空间是一般空间在傅里叶转换下的对偶空间,是以频率空间为坐标系承载磁共振图像原始回波数据的空间,计算机设备可以对保存在K空间中的回波数据进行傅里叶变换,得到时间域函数,解析出各质子的空间定位,从而得到磁共振图像。可以理解的是,磁共振成像中采集的数据是频率域函数,而为了空间定位,需要用到几个方向的梯度对空间进行编码,因此需要对保存在K空间中的数据进行傅里叶变化,得到时间域函数,并解析出每个质子的空间定位,得到磁共振图像。可选地,磁体、射频发射线圈、梯度线圈和射频接收线圈都可用以控制系统控制。仅作为示例,梯度线圈可包括第一方向上的第一磁场梯度、第二方向上的第二磁场梯度,以及第三方向上的第三磁场梯度。在一些实施例中,第一、第二和第三方向可分别沿X轴、Y轴和Z轴。在一些实施例中,沿X轴、Y轴和/或Z轴的磁场梯度可对应于k空间中的不同编码/读出方向(例如,kx轴的方向、ky轴的方向、kz轴的方向、或任何其他方向),控制系统可控制X方向、Y方向和Z方向上的磁场梯度,从而进行选层、位置定位等。在一些实施例中,控制系统可接收由例如用户提供的命令,并根据收到的命令来调节磁体和/或射频发射线圈、梯度线圈以拍摄感兴趣对象的图像。仅作为示例,该命令可涉及磁场梯度(例如,散相梯度)的极性、波形、强度和/或定时。在一些实施例中,射频发射线圈发射射频脉冲的时序、ADC采集窗的时序以及梯度线圈产生梯度脉冲的时序可由控制系统发出的脉冲序列确定。脉冲序列可以是自旋回波(SE)序列、快速自旋回波(FSE)序列、超短回波时间(UTE)序列等。
在本实施例中,射频发射线圈发射的射频脉冲的整体带宽覆盖范围是由检测对象携带的金属植入物引起的场不均匀性程度确定的,当主磁场不均匀时,会使得质子的进动频率偏置分布,由于发射的射频脉冲的整体带宽范围足够大,能够使得发射的射频脉冲覆盖不同的进动频率范围,最终使得成像区域的质子都能够被有效激发,并采用预设的读出带宽接收相应的回波信号,并保存在对应的K空间中,进而对k空间中的回波数据进行傅里叶变换得到的磁共振图像,能够更加准确的反映扫描体的有效信息,减少因植入物带来的场不均匀性造成的图像伪影。
在上述发射射频脉冲的场景中,发射的射频脉冲的带宽与角度具有一定的相关关系。在上述实施例的基础上,作为一种可选的实施方式,射频脉冲的带宽越高,射频脉冲的角度越小。
具体的,上述发射的射频脉冲中,射频脉冲的带宽越高,该射频脉冲的角度越小。可选的,发射的射频脉冲可以包括激发射频脉冲和回聚射频脉冲。例如,若激发射频脉冲的带宽越高,则激发射频脉冲的角度越小;若回聚射频脉冲的带宽越高,则回聚射频脉冲的角度越小。
在本实施例中,发射的射频脉冲的带宽越高,射频脉冲的角度越小,而带宽越高的射频脉冲能够覆盖质子不同的进动频率范围,从而使得采用预设的读出带宽接收的回波信号为扫描体不同区域的完整的信号,使得保存在K空间中的回波数据为不同区域的完整的回波信号对应的回波数据。
在上述发射射频脉冲的场景中,在上述实施例的基础上,作为一种可选的实施方式,射频脉冲包括一个中心频率,相应的,在该中心频率下采集得到的回波信号对应一个数据集。具体的,在上述发射的射频脉冲中,射频脉冲包括一个中心频率,该中心频率下激发得到的回波数据对应一个数据集。可选的,发射的射频脉冲可以包括均为高带宽、小角度的激发射频脉冲和回聚射频脉冲,即激发射频脉冲的角度和回聚射频脉冲的角度均小于预设的角度阈值,激发射频脉冲的带宽和回聚射频脉冲的带宽均大于预设的第一带宽阈值。可选的,预设的角度阈值可以为60度,预设的第一带宽阈值可以为3kHZ。
在本实施例中,发射的射频脉冲包括一个中心频率,发射的射频脉冲可以包括均为小角度、高带宽的激发射频脉冲和回聚射频脉,能够使得发射的射频脉冲覆盖不同的进动频率范围,从而使得采用预设的读出带宽接收的回波信号为扫描体不同区域的完整的信号,使得保存在K空间中的回波数据为不同区域的完整的回波信号对应的回波数据。
图3为一个实施例提供的发射的射频脉冲的示意图。图4为一个实施例提供的发射的射频脉冲的示意图。在上述发射射频脉冲的场景中,发射的射频脉冲可以包括多个低带宽的激发射频脉冲和一个小角度、高带宽的回聚射频脉冲或者,可以包括一个小角度、高带宽的激发射频脉冲和多个低带宽的回聚射频脉冲。在上述实施例的基础上,作为一种可选的实施方式,射频脉冲包括至少两个中心频率,相应的,该两个中心频率下激发核自旋产生的回波信号对应至少两个数据集,每个中心频率对应一个数据集。
具体的,在上述发射的射频脉冲中,射频脉冲包括至少两个中心频率,相对应的,检测对象被激发产生的回波信号对应至少两个数据集,且每个中心频率对应一个数据集。可选的,如图3所示,可以是激发射频脉冲包括多个中心频率,回聚射频脉冲包括一个中心频率,也就是发射的射频脉冲可以是多个低带宽的激发射频脉冲和一个小角度、高带宽的回聚射频脉冲,即多个激发射频脉冲的带宽均小于预设的第一带宽阈值,回聚射频脉冲的角度小于预设的角度阈值,回聚射频脉冲的带宽大于预设的第一带宽阈值,且激发射频脉冲的各中心频率对应不同的频带,各频带整体范围覆盖回聚射频脉冲所对应的带宽范围;如图4所示,也可以是回聚射频脉冲包括多个中心频率,激发射频脉冲包括一个中心频率,也就是发射的射频脉冲可以是一个小角度、高带宽的激发射频脉冲和多个低带宽的回聚射频脉冲,即激发射频脉冲的角度小于预设的角度阈值,激发射频脉冲的带宽大于预设的第一带宽阈值,多个回聚射频脉冲的带宽均小于预设的第一带宽阈值,且回聚射频脉冲的各中心频率对应不同的频带,各频带整体范围能够覆盖激发射频脉冲所对应的带宽范围。可选的,预设的角度阈值可以为60度,预设的第一带宽阈值可以为3kHZ。可选的,在本实施例中,计算机设备控制射频发射线圈多次发射射频脉冲,并控制射频接收线圈多次采用预设的读出带宽接收回波信号,得到多个数据集。可选的,预设的读出带宽大于预设的带宽阈值。可选的,预设的带宽阈值可以为500Hz。
在本实施例中,发射的射频脉冲包括至少两个中心频率,相对应的回波信号对应至少两个数据集,每个中心频率下对应产生一个数据集,能够使得回聚射频脉冲的整体带宽范围覆盖激发脉冲的带宽范围,保证发射的射频脉冲覆盖不同的进动频率范围,从而使得采用预设的读出带宽接收的回波信号为扫描体不同区域的完整的信号,使得保存在K空间中的回波数据为不同区域的完整的回波信号对应的回波数据。
在上述实施例的基础上,作为一种可选的实施方式,上述S203中的将回波信号对应的回波数据保存在K空间中之前,上述方法还包括:对回波数据进行磁化处理和/或压脂处理,得到处理后的回波数据;磁化处理用于对回波数据进行T2加权;压脂处理用于消除回波数据中的脂肪数据。
具体的,计算机设备将回波信号对应的回波数据保存在K空间中之前,计算机设备会对回波数据进行磁化处理和/或压脂处理,得到处理后的回波数据,其中,磁化处理用于对回波数据进行T2加权,压脂处理用于消除回波数据中的脂肪数据,对回波数据进行处理后,计算机设备将处理后的回波数据保存在K空间中。
在本实施例中,计算机设备将回波信号对应的回波数据保存在K空间中之前,对回波数据进行的磁化处理和/或压脂处理能够消除回波数据中的干扰数据,使得经过处理的回波数据更加地准确,进而使得保存在K空间中的回波数据更加的准确。
图5为另一个实施例提供的磁共振成像方法的流程示意图。在上述射频脉冲包括至少两个中心频率,回波信号对应至少两个数据集的场景中。如图5所示,将至少两个数据集分别填充至至少两个K空间,对K空间回波数据进行傅里叶变换,得到磁共振图像,包括:
S501,分别对至少两个K空间进行傅里叶变换,得到至少两个三维子图像。
具体的,计算机设备将至少两个数据集分别填充至至少两个K空间后,计算机设备分别对各K空间进行傅里叶变换,得到时间域函数,并解析出每个质子的空间定位,得到至少两个三维子图像。
S502,对各三维子图像进行复合,获取磁共振图像。
具体的,计算机设备分别对各K空间进行傅里叶变换,得到至少两个三维子图像后,对各三维子图像进行复合,得到磁共振图像。可选的,对各三维子图像进行复合可以是进行图像拼接处理。
在本实施例中,计算机设备将至少两个数据集分别填充至至少两个K空间后,分别对各K空间进行傅里叶变换,得到至少两个三维子图像,由于得到的数据集中包含的数据较多,进而提高了得到的三维子图像的准确度,而磁共振图像是对各三维子进行复合处理得到的,进而提高了得到的磁共振图像的准确度。
应该理解的是,虽然图2和图5的流程图中的各个步骤按照箭头的指示依次显示,但是这些步骤并不是必然按照箭头指示的顺序依次执行。除非本文中有明确的说明,这些步骤的执行并没有严格的顺序限制,这些步骤可以以其它的顺序执行。而且,图2和图5中的至少一部分步骤可以包括多个子步骤或者多个阶段,这些子步骤或者阶段并不必然是在同一时刻执行完成,而是可以在不同的时刻执行,这些子步骤或者阶段的执行顺序也不必然是依次进行,而是可以与其它步骤或者其它步骤的子步骤或者阶段的至少一部分轮流或者交替地执行。
图6为一个实施例提供的磁共振成像系统结构示意图。如图6所示,该系统可以包括:磁体、射频发射线圈、射频接收线圈和控制系统;磁体,用于产生主磁场,并形成容纳检测对象的孔腔,且所述检测对象携带金属植入物;射频发射线圈,用于发射射频脉冲;射频接收线圈,用于接收磁共振信号;控制系统,用于执行上述实施例所述的磁共振成像方法,其实现原理和技术效果类似,在此不再赘述。
本发明实施例提供的磁共振成像方法,可以适用于如图7所示的计算机设备。该计算机设备包括通过系统总线连接的处理器、存储器,该存储器中存储有计算机程序,处理器执行该计算机程序时可以执行下述方法实施例的步骤。可选的,该计算机设备还可以包括网络接口、显示屏和输入装置。其中,该计算机设备的处理器用于提供计算和控制能力。该计算机设备的存储器包括非易失性存储介质、内存储器,该非易失性存储介质存储有操作系统和计算机程序。该内存储器为非易失性存储介质中的操作系统和计算机程序的运行提供环境。该计算机设备的网络接口用于与外部的终端通过网络连接通信。可选的,该计算机设备可以是个人计算机,个人数字助理,平板电脑、手机等等,还可以是云端或者远程服务器,本发明实施例对计算机设备的具体形式并不做限定。
在一个实施例中,提供了一种计算机设备,包括存储器和处理器,存储器中存储有计算机程序,该处理器执行计算机程序时实现以下步骤:
控制磁体产生主磁场;检测对象位于主磁场形成的视场区域内,检测对象携带金属植入物,且金属植入物使主磁场非均匀;
控制射频发射线圈发射射频脉冲;射频脉冲的整体带宽覆盖范围是由金属植入物引起的场不均匀性程度确定的;
采用预设的读出带宽控制射频接收线圈接收回波信号,将回波信号对应的回波数据保存在K空间中,对K空间中的回波数据进行傅里叶变换,得到磁共振图像。
上述实施例提供的计算机设备,其实现原理和技术效果与上述方法实施例类似,在此不再赘述。
在一个实施例中,提供了一种计算机可读存储介质,其上存储有计算机程序,计算机程序被处理器执行时实现以下步骤:
控制磁体产生主磁场;检测对象位于主磁场形成的视场区域内,检测对象携带金属植入物,且金属植入物使主磁场非均匀;
控制射频发射线圈发射射频脉冲;射频脉冲的整体带宽覆盖范围是由金属植入物引起的场不均匀性程度确定的;
采用预设的读出带宽控制射频接收线圈接收回波信号,将回波信号对应的回波数据保存在K空间中,对K空间中的回波数据进行傅里叶变换,得到磁共振图像。
上述实施例提供的计算机可读存储介质,其实现原理和技术效果与上述方法实施例类似,在此不再赘述。
本领域普通技术人员可以理解实现上述实施例方法中的全部或部分流程,是可以通过计算机程序来指令相关的硬件来完成,所述的计算机程序可存储于一非易失性计算机可读取存储介质中,该计算机程序在执行时,可包括如上述各方法的实施例的流程。其中,本申请所提供的各实施例中所使用的对存储器、存储、数据库或其它介质的任何引用,均可包括非易失性和/或易失性存储器。非易失性存储器可包括只读存储器(ROM)、可编程ROM(PROM)、电可编程ROM(EPROM)、电可擦除可编程ROM(EEPROM)或闪存。易失性存储器可包括随机存取存储器(RAM)或者外部高速缓冲存储器。作为说明而非局限,RAM以多种形式可得,诸如静态RAM(SRAM)、动态RAM(DRAM)、同步DRAM(SDRAM)、双数据率SDRAM(DDRSDRAM)、增强型SDRAM(ESDRAM)、同步链路(Synchlink)DRAM(SLDRAM)、存储器总线(Rambus)直接RAM(RDRAM)、直接存储器总线动态RAM(DRDRAM)、以及存储器总线动态RAM(RDRAM)等。
以上所述实施例的各技术特征可以进行任意的组合,为使描述简洁,未对上述实施例中的各个技术特征所有可能的组合都进行描述,然而,只要这些技术特征的组合不存在矛盾,都应当认为是本说明书记载的范围。
以上所述实施例仅表达了本发明的几种实施方式,其描述较为具体和详细,但并不能因此而理解为对发明专利范围的限制。应当指出的是,对于本领域的普通技术人员来说,在不脱离本发明构思的前提下,还可以做出若干变形和改进,这些都属于本发明的保护范围。因此,本发明专利的保护范围应以所附权利要求为准。

Claims (8)

1.一种磁共振成像方法,其特征在于,所述方法包括:
控制磁体产生主磁场;检测对象位于所述主磁场形成的视场区域内,所述检测对象携带金属植入物,且所述金属植入物使所述主磁场非均匀;
控制射频发射线圈发射射频脉冲;所述射频脉冲的整体带宽覆盖范围是由所述金属植入物引起的场不均匀性程度确定的;所述射频脉冲包括激发射频脉冲和回聚射频脉冲,所述激发射频脉冲包括多个中心频率,所述回聚射频脉冲包括一个中心频率,且所述激发射频脉冲的各中心频率对应不同的频带,各频带整体范围覆盖所述回聚射频脉冲所对应的带宽范围;或者,所述回聚射频脉冲包括多个中心频率,所述激发射频脉冲包括一个中心频率,且所述回聚射频脉冲的各中心频率对应不同的频带,各频带整体范围能够覆盖所述激发射频脉冲所对应的带宽范围;
采用预设的读出带宽控制射频接收线圈接收回波信号,将所述回波信号对应的回波数据保存在K空间中,对所述K空间中的回波数据进行傅里叶变换,得到磁共振图像。
2.根据权利要求1所述的方法,其特征在于,所述射频脉冲包括至少两个中心频率,相应的,所述回波信号对应至少两个数据集,每个中心频率对应一个数据集。
3.根据权利要求2所述的方法,其特征在于,所述至少两个数据集分别填充至至少两个K空间,对所述K空间回波数据进行傅里叶变换,得到磁共振图像,包括:
分别对所述各K空间进行傅里叶变换,得到至少两个三维子图像;
对所述各三维子图像进行复合,获取所述磁共振图像。
4.根据权利要求2所述的方法,其特征在于,所述方法还包括:
控制所述射频发射线圈多次发射射频脉冲,并采用所述预设的读出带宽控制所述射频接收线圈多次接收回波信号,得到多个数据集。
5.根据权利要求4所述的方法,其特征在于,所述预设的读出带宽大于预设的带宽阈值。
6.一种磁共振成像系统,其特征在于,所述系统包括:磁体、射频发射线圈、射频接收线圈和控制系统;
所述磁体,用于产生主磁场,并形成容纳检测对象的孔腔,且所述检测对象携带金属植入物;
所述射频发射线圈,用于发射射频脉冲;
所述射频接收线圈,用于接收磁共振信号;
所述控制系统,用于执行如权利要求1-5任一项所述的磁共振成像方法。
7.一种计算机设备,包括存储器和处理器,所述存储器存储有计算机程序,其特征在于,所述处理器执行所述计算机程序时实现权利要求1-5中任一项所述方法的步骤。
8.一种计算机可读存储介质,其上存储有计算机程序,其特征在于,所述计算机程序被处理器执行时实现权利要求1-5中任一项所述方法的步骤。
CN201910452080.7A 2019-05-28 2019-05-28 磁共振成像方法、系统、计算机设备和可读存储介质 Active CN110215210B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910452080.7A CN110215210B (zh) 2019-05-28 2019-05-28 磁共振成像方法、系统、计算机设备和可读存储介质

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910452080.7A CN110215210B (zh) 2019-05-28 2019-05-28 磁共振成像方法、系统、计算机设备和可读存储介质

Publications (2)

Publication Number Publication Date
CN110215210A CN110215210A (zh) 2019-09-10
CN110215210B true CN110215210B (zh) 2023-08-22

Family

ID=67818175

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910452080.7A Active CN110215210B (zh) 2019-05-28 2019-05-28 磁共振成像方法、系统、计算机设备和可读存储介质

Country Status (1)

Country Link
CN (1) CN110215210B (zh)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111257810B (zh) * 2020-01-19 2021-11-23 上海联影医疗科技股份有限公司 磁共振信号采集系统、方法
CN112806982A (zh) * 2020-12-30 2021-05-18 深圳市联影高端医疗装备创新研究院 磁共振成像方法、装置、设备、存储介质及系统
CN112986399B (zh) * 2021-03-15 2022-06-28 南昌航空大学 一种电磁超声sh导波换能器及在线检测系统、方法
CN114236439A (zh) * 2021-11-09 2022-03-25 深圳市联影高端医疗装备创新研究院 线圈定位方法、装置、磁共振设备及存储介质

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN85103309A (zh) * 1985-05-14 1986-11-12 菲利浦光灯制造公司 确定一物体中的一个区域中的核磁共振(nmr)分布的快速方法与设备
JP2006000265A (ja) * 2004-06-16 2006-01-05 Hitachi Medical Corp 磁気共鳴イメージング装置
CN105531597A (zh) * 2013-09-10 2016-04-27 皇家飞利浦有限公司 耐金属的mr成像
CN105662415A (zh) * 2016-03-03 2016-06-15 哈尔滨医科大学 一种多核磁共振成像系统
CN107850649A (zh) * 2015-06-26 2018-03-27 皇家飞利浦有限公司 用于检测金属植入物并且为高效mri工作流选择磁共振脉冲序列的方法和检测单元

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7821264B2 (en) * 2008-01-18 2010-10-26 General Electric Company System and method for MR imaging in inhomogeneous magnetic fields
US8274286B2 (en) * 2008-01-18 2012-09-25 General Electric Company System and method for multi-spectral MR imaging near metal
US9594134B2 (en) * 2014-01-22 2017-03-14 Wisconsin Alumni Research Foundation System and method for fully phase-encoded magnetic resonance imaging using multiband radio frequency excitation
US10725133B2 (en) * 2016-04-13 2020-07-28 The Board Of Trustees Of The Leland Stanford Junior University Field-mapping and artifact correction in multispectral imaging
US10613171B2 (en) * 2017-03-14 2020-04-07 Siemens Healthcare Gmbh Multi-banded RF-pulse enhanced magnetization imaging

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN85103309A (zh) * 1985-05-14 1986-11-12 菲利浦光灯制造公司 确定一物体中的一个区域中的核磁共振(nmr)分布的快速方法与设备
JP2006000265A (ja) * 2004-06-16 2006-01-05 Hitachi Medical Corp 磁気共鳴イメージング装置
CN105531597A (zh) * 2013-09-10 2016-04-27 皇家飞利浦有限公司 耐金属的mr成像
CN107850649A (zh) * 2015-06-26 2018-03-27 皇家飞利浦有限公司 用于检测金属植入物并且为高效mri工作流选择磁共振脉冲序列的方法和检测单元
CN105662415A (zh) * 2016-03-03 2016-06-15 哈尔滨医科大学 一种多核磁共振成像系统

Also Published As

Publication number Publication date
CN110215210A (zh) 2019-09-10

Similar Documents

Publication Publication Date Title
CN110215210B (zh) 磁共振成像方法、系统、计算机设备和可读存储介质
KR101663367B1 (ko) 물의 t1 시간 및 지방의 t1 시간을 결정하는 자기 공명 시스템 및 방법
US10274569B2 (en) Magnetic resonance imaging system with ghost artifact reduction and method operation thereof
US20090273345A1 (en) Method for determining the spatial distribution of magnetic resonance signals with use of local spatially encoding magnetic fields
US9594141B2 (en) System for accelerated magnetic resonance imaging using parallel coils
US10317487B2 (en) Method and magnetic resonance apparatus for acquisition of MR data of a two-dimensional volume segment
US10359487B2 (en) Zero echo time MR imaging
US20150268322A1 (en) Method and magnetic resonance apparatus for the acquisition of mr data of a slice within a subject
US20050083054A1 (en) Method and magnetic resonance imaging apparatus for compensating contrast inhomogeneities in magnetic resonance images
US8310234B2 (en) System for fat suppression in MR image acquisition
US8456163B2 (en) Variable image resolution MR imaging system
US10794980B2 (en) Cross-term spatiotemporal encoding for magnetic resonance imaging
US20060244450A1 (en) Method for shimming a main magnetic field in magnetic resonance
US10502804B2 (en) Magnetic resonance imaging apparatus and method of controlling the same
US20170261583A1 (en) Method and apparatus for sectional optimization of radial mr pulse sequences
US10151811B2 (en) Method for producing an MR image and correspondingly designed magnetic resonance scanner
US10371777B2 (en) MRI propeller with motion correction, water-fat separation and estimation of magnetic field inhomogeneity information
US20160054411A1 (en) Method and apparatus for acquiring magnetic resonance data
CN112826491B (zh) 磁共振成像方法及系统
US11169236B2 (en) Phase correction for echo-planar imaging
US10145926B2 (en) Phase-cycled steady-state free precession imaging
US20160124065A1 (en) Method and apparatus for correction of magnetic resonance image recordings with the use of a converted field map
CN110488210B (zh) 磁共振系统的扫描方法、计算机设备和存储介质
US10578697B2 (en) Method and magnetic resonance apparatus using a spin-echo sequence for the spatially selective acquisition of magnetic resonance data
US20160349343A1 (en) Method and apparatus for the correction of magnetic resonance scan data

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

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.

CB02 Change of applicant information
GR01 Patent grant
GR01 Patent grant