CN111273207B - 磁共振成像方法、装置、电子设备、存储介质 - Google Patents

磁共振成像方法、装置、电子设备、存储介质 Download PDF

Info

Publication number
CN111273207B
CN111273207B CN202010075000.3A CN202010075000A CN111273207B CN 111273207 B CN111273207 B CN 111273207B CN 202010075000 A CN202010075000 A CN 202010075000A CN 111273207 B CN111273207 B CN 111273207B
Authority
CN
China
Prior art keywords
space
magnetic resonance
target
data
resonance imaging
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
CN202010075000.3A
Other languages
English (en)
Other versions
CN111273207A (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.)
Neusoft Medical Systems Co Ltd
Shanghai Neusoft Medical Technology Co Ltd
Original Assignee
Neusoft Medical Systems Co Ltd
Shanghai Neusoft Medical Technology 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 Neusoft Medical Systems Co Ltd, Shanghai Neusoft Medical Technology Co Ltd filed Critical Neusoft Medical Systems Co Ltd
Priority to CN202010075000.3A priority Critical patent/CN111273207B/zh
Publication of CN111273207A publication Critical patent/CN111273207A/zh
Application granted granted Critical
Publication of CN111273207B publication Critical patent/CN111273207B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • 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/4818—MR characterised by data acquisition along a specific k-space trajectory or by the temporal order of k-space coverage, e.g. centric or segmented coverage of k-space
    • 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/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
    • 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/5608—Data processing and visualization specially adapted for MR, e.g. for feature analysis and pattern recognition on the basis of measured MR data, segmentation of measured MR data, edge contour detection on the basis of measured MR data, for enhancing measured MR data in terms of signal-to-noise ratio by means of noise filtering or apodization, for enhancing measured MR data in terms of resolution by means for deblurring, windowing, zero filling, or generation of gray-scaled images, colour-coded images or images displaying vectors instead of pixels

Landscapes

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

Abstract

本发明公开了磁共振成像方法、装置、电子设备、存储介质。方法包括：在对目标对象进行至少两个扫描周期的扫描过程中，基于径向轨迹采集针对所述目标对象的磁共振数据，其中，每个扫描周期包含多个时间段；针对每个扫描周期，将相同时间段的磁共振数据映射于一个K空间内，在每个K空间中形成至少一条径向线；从填补K空间中获取第一预设区域内的采样数据点，并填补于所述目标K空间的数据缺失区域；基于经过数据填补的目标K空间重建针对所述目标对象的磁共振图像。从而提高了图像重建的精度。

Description

磁共振成像方法、装置、电子设备、存储介质
技术领域
本发明涉及医学成像技术领域，特别涉及磁共振成像方法、装置、电子设备、存储介质。
背景技术
磁共振成像(Magnetic Resonance Imaging，MRI)是现代医学影像中主要的成像方式之一，其基本工作原理是利用磁共振现象，采用激发被检测对象中的氢质子，利用梯度场进行位置编码，随后使用接收线圈接收带位置信息的信号，最终通过傅里叶变换重建出图像信息。
为了重建出高分辨的图像，对于径向轨迹采样得到的K空间数据，目前一般采用KWIC(k-space weighted image contrast，K空间图像加权)重建算法进行图像重建，KWIC重建算法将来源于不同采样时间段的多个K空间进行加权处理。随着采样时间间隔越大，因被检测对象的运动带来的高频信息差异也越大，单纯地将多个K空间进行加权处理，进而进行图像重建，得到的图像精度较低。
发明内容
本发明提供磁共振成像方法、装置、电子设备、存储介质，以提高图像重建的精度。
具体地，本发明是通过如下技术方案实现的：
第一方面，提供一种磁共振成像方法，所述磁共振成像方法包括：
在对目标对象进行至少两个扫描周期的扫描过程中，基于径向轨迹采集针对所述目标对象的磁共振数据，其中，每个扫描周期包含多个时间段，在所述扫描过程中，对于不同扫描周期中对应的时间段，所述目标对象处于相同的呼吸状态；
针对每个扫描周期，将相同时间段的磁共振数据映射于一个K空间内，在每个K空间中形成至少一条径向线，所述径向线包含多个表征频率的采样数据点；
从填补K空间中获取第一预设区域内的采样数据点，并填补于所述目标K空间的数据缺失区域，其中，所述目标K空间为多个K空间中的一个，所述填补K空间为不同扫描周期中与所述目标K空间的时间段对应的时间段的K空间；
基于经过数据填补的目标K空间重建针对所述目标对象的磁共振图像。
可选地，所述磁共振成像方法还包括：
从其他K空间中获取第二预设区域内的采样数据点，并填补于所述目标K空间的数据缺失区域，其中，所述其他K空间为所述目标K空间的扫描周期中除所述目标K空间之外的K空间。
可选地，所述磁共振成像方法还包括：
从相邻K空间中获取第二预设区域内的采样数据点，并填补于所述目标K空间的数据缺失区域，其中，所述相邻K空间为所述目标K空间的扫描周期中与所述目标K空间的时间段相邻的K空间。
可选地，所述K空间包含中心区域、环绕所述中心区域的中间区域以及环绕所述中间区域的边缘区域；
所述第一预设区域包含所述边缘区域；
所述第二区域包含所述中间区域和所述边缘区域。
可选地，采集针对所述目标对象的磁共振数据，包括：
基于呼吸门控技术采集所述磁共振数据。
第二方面，提供一种磁共振成像装置，所述磁共振成像装置包括：
采集模块，用于在对目标对象进行至少两个扫描周期的扫描过程中，基于径向轨迹采集针对所述目标对象的磁共振数据，其中，每个扫描周期包含多个时间段，在所述扫描过程中，对于不同扫描周期中对应的时间段，所述目标对象处于相同的呼吸状态；
映射模块，用于针对每个扫描周期，将相同时间段的磁共振数据映射于一个K空间内，在每个K空间中形成至少一条径向线，所述径向线包含多个表征频率的采样数据点；
填补模块，用于从填补K空间中获取第一预设区域内的采样数据点，并填补于所述目标K空间的数据缺失区域，其中，所述目标K空间为多个K空间中的一个，所述填补K空间为不同扫描周期中与所述目标K空间的时间段对应的时间段的K空间；
重建模块，用于基于经过数据填补的目标K空间重建针对所述目标对象的磁共振图像。
可选地，所述填补模块还用于：
从其他K空间中获取第二预设区域内的采样数据点，并填补于所述目标K空间的数据缺失区域，其中，所述其他K空间为所述目标K空间的扫描周期中除所述目标K空间之外的K空间。
可选地，所述填补模块还用于：
从相邻K空间中获取第二预设区域内的采样数据点，并填补于所述目标K空间的数据缺失区域，其中，所述相邻K空间为所述目标K空间的扫描周期中与所述目标K空间的时间段相邻的K空间。
可选地，所述K空间包含中心区域、环绕所述中心区域的中间区域以及环绕所述中间区域的边缘区域；
所述第一预设区域包含所述边缘区域；
所述第二区域包含所述中间区域和所述边缘区域。
可选地，所述采集模块具体用于：
基于呼吸门控技术采集所述磁共振数据。
第三方面，提供一种电子设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现上述任一项所述的磁共振成像方法。
第四方面，提供一种计算机可读存储介质，其上存储有计算机程序，所述计算机程序被处理器执行时实现上述任一项所述的磁共振成像方法的步骤。
本发明的实施例提供的技术方案可以包括以下有益效果：
本发明实施例中，充分地利用了采集的磁共振数据在多个时间维度的高低频信息的相似性，对K空间进行数据填补，可以在维持图像时间分辨率和空间分辨率的前提下，提高现有重建算法的对比度和保真度，进而有效地提高了重建图像的精确度，并具有较强的鲁棒性。
应当理解的是，以上的一般描述和后文的细节描述仅是示例性和解释性的，并不能限制本发明。
附图说明
此处的附图被并入说明书中并构成本说明书的一部分，示出了符合本发明的实施例，并与说明书一起用于解释本发明的原理。
图1是本发明一示例性实施例示出的磁共振成像系统的结构示意图；
图2a是本发明一示例性实施例示出的将基于径向轨迹采集的磁共振数据映射于K空间的第一结果示意图；
图2b是本发明一示例性实施例示出的将基于径向轨迹采集的磁共振数据映射于K空间的第二结果示意图；
图2c是本发明一示例性实施例示出的将基于径向轨迹采集的磁共振数据映射于K空间的第三结果示意图；
图2d是本发明一示例性实施例示出的将基于径向轨迹采集的磁共振数据映射于K空间的第四结果示意图；
图2e是现有技术中对图2a～图2d示出的K空间进行加权处理的结果示意图；
图3a是本发明一示例性实施例示出的一种磁共振成像方法的流程图；
图3b是本发明一示例性实施例示出的表征采用呼吸门控技术获取的3个扫描周期的磁共振数据的图像示意图；
图3c是使用图3a示出的磁共振成像方法对一个K空间进行数据填补的结果示意图；
图4a是本发明一示例性实施例示出的另一种磁共振成像方法的流程图；
图4b是使用图4a示出的磁共振成像方法对一个K空间进行数据填补的结果示意图；
图5a是本发明一示例性实施例示出的另一种磁共振成像方法的流程图；
图5b是使用图5a示出的磁共振成像方法对一个K空间进行数据填补的结果示意图；
图6是本发明一示例性实施例示出的一种磁共振成像装置的模块示意图；
图7是本发明一示例实施例示出的一种电子设备的结构示意图。
具体实施方式
这里将详细地对示例性实施例进行说明，其示例表示在附图中。下面的描述涉及附图时，除非另有表示，不同附图中的相同数字表示相同或相似的要素。以下示例性实施例中所描述的实施方式并不代表与本发明相一致的所有实施方式。相反，它们仅是与如所附权利要求书中所详述的、本发明的一些方面相一致的装置和方法的例子。
在本发明使用的术语是仅仅出于描述特定实施例的目的，而非旨在限制本发明。在本发明和所附权利要求书中所使用的单数形式的“一种”、“所述”和“该”也旨在包括多数形式，除非上下文清楚地表示其他含义。还应当理解，本文中使用的术语“和/或”是指并包含一个或多个相关联的列出项目的任何或所有可能组合。
应当理解，尽管在本发明可能采用术语第一、第二、第三等来描述各种信息，但这些信息不应限于这些术语。这些术语仅用来将同一类型的信息彼此区分开。例如，在不脱离本发明范围的情况下，第一信息也可以被称为第二信息，类似地，第二信息也可以被称为第一信息。取决于语境，如在此所使用的词语“如果”可以被解释成为“在……时”或“当……时”或“响应于确定”。
图1是本发明一示例性实施例示出的磁共振成像(MRI)系统的结构示意图，如图1所示，MRI系包括MRI设备110、显示操作设备120、数据处理器130。MRI设备110包括磁体111、梯度线圈112、射频线圈113。
磁体111产生用于将被检测对象(可以是人或动物，或者人或动物的部位)中的原子核的磁偶极矩的方向调整为恒定方向的静磁场。
梯度线圈112包括用于沿着彼此相交成直角的X轴方向、Y轴方向和Z轴方向产生磁场梯度的X线圈、Y线圈和Z线圈。梯度线圈112可通过根据被检测对象的区域而不同地诱导共振频率来提供被检测对象的空间定位信息。具体的，用一个方向的梯度磁场作为层面选择梯度，确定扫描层面，然后用另外两个方向的梯度磁场来确定层面内的坐标位置，通过三个梯度磁场的不同组合，可以实现任意层面断层成像。
射频线圈113可向被检测对象发送RF(射频)信号，并采集从检测对象发射的MR(磁共振)数据。具体的，射频线圈113产生具有与原子核的类型相应的RF信号，并且将RF信号施加到被检测对象，使得检测对象的原子核从低能态跃迁到高能态。当由射频线圈113产生的RF信号消失时，原子核从高能态跃迁到低能态，由此发射具有拉莫尔频率的电磁波(MR信号)。换句话讲，当停止向检测对象施加RF信号时，检测对象的原子核能级从高能级改变到低能级，因此发射具有拉莫尔频率的电磁波，此时射频线圈113可采集MR数据。
数据处理器130可对射频线圈113采集的MR数据进行处理，以重建成图像。显示操作设备120可显示数据处理器130重建的图像。
对于MR数据的采集，传统采用笛卡尔轨迹，它的磁共振成像方法虽然简单，但是逐行采样对运动(肌肉收缩、呼吸运动等)造成的伪影很敏感。为了解决伪影问题，基于径向轨迹的磁共振成像技术发展越来越迅速，其能够减少伪影数量和数据采集时间，在动态加速成像、运动伪影消除、超短回波时间成像等方面发挥着重要的优势。
磁共振成像过程中，需要将采集的磁共振数据映射到K空间，也即将原始采集的MR数据(时间域数据)映射到频率域。图像重建时，对K空间进行傅里叶变换，即可得到重建图像。
图2a～图2d示出了一个扫描周期内获得的4个K空间，每个K空间中均有2条投影角度交错的径向线。若一个扫描周期包含4个时间段，在该扫描周期内基于径向轨迹采集MR数据，将同一时间段采集的MR数据进行编码后填补于对应的K空间，即可得到4个K空间。针对一个K空间进行数据填补时，按照数据采样轨迹将对应的MR数据填补于该K空间，即可得到至少一条径向线。每条径向线由多个采样数据点组成，每个采样数据点表征与MR数据对应的空间频率。每个采样数据点由kx和ky表示，kx和ky分别代表两个相互垂直的空间频率。其中，K空间中的空白区域(没有径向线的区域)表示未被采样的区域，也即数据缺失区域。
参见图2a，为了便于下文描述，将K空间划分为中心区域(圈a围成的区域)、环绕中心区域的中间区域(圈b围成的区域)以及环绕中间区域的边缘区域(圈b的外围区域)。其中，K空间的划分方式可根据实际需求进行划分。
在K空间中，每条径向线含有低频到高频信息。越靠近K空间的中心区域的数据采样点具有越低的空间频率值；越靠近K空间的边缘区域的数据采样点具有越高的空间频率值。高频信息决定了重建图像的空间分辨率和细节信息，例如，被检测对象的位移、形态变化；低频信息决定了重建图像的对比度信息。
在图像重建时，由于K空间中的采样数据较少，重建得到的图像往往空间分辨率较低，为了重建出分辨率高的图像，需要对K空间中的数据缺失区域(空白区域)进行数据填补。相关技术中，将一个扫描周期内获得的所有K空间进行加权处理，得到加权后的K空间，进而进行图像重建。以对图2a～图2d示出的K空间进行加权处理为例，得到的结果参见图2e所示，不同时间段的K空间中心区域的信息是采用各自对应时间段采集的信息，而K空间外部区域的信息则是将同一周期所有时间段的信息进行填充。该方式的缺点是，同一周期中时间间隔越大，由于被检测对象的运动(肌肉收缩、呼吸运动等)带来的信息差异也越大，致使不同时间段K空间的高频信息的填充结果会存在较大误差，影响最终的图像重建精度。
基于此，本发明提供了一种磁共振成像方法，以提高对基于径向轨迹采集磁共振数据进行图像重建的精确度。以下结合图1所示的MRI系统以及图2a～图2d，对本发明的磁共振成像方法实施例进行详细描述。
图3a是本发明一示例性实施例示出的一种磁共振成像方法的流程图，方法包括以下步骤：
步骤301、在对目标对象进行至少两个扫描周期的扫描过程中，基于径向轨迹采集针对目标对象的磁共振数据。
其中，每个扫描周期包含多个时间段。在扫描过程中，对于不同扫描周期中相对应的时间段，目标对象(被检测对象)处于相同的呼吸状态，对应的形态结构相同。
步骤301中，基于径向轨迹采集磁共振数据的过程中，可以但不限于利用呼吸门控技术，以确保不同扫描周期中对应的时间段采集的磁共振数据，是在目标对象处于相同的呼吸状态下采集的。
参见图3b，为采用呼吸门控技术获取的3个扫描周期的磁共振数据，图中每个6边形表征一个时间段采集的磁共振数据，填充图案相同的六边形表征目标对象处于相同的呼吸状态(形态结构相同)下采集的磁共振数据。
步骤302、针对每个扫描周期，将相同时间段的磁共振数据映射于一个K空间内。
假设针对目标对象进行了m个扫描周期的扫描，每个扫描周期划分为n个时间段，将针对目标对象的磁共振数据填充于K空间，即可得到m*n个K空间。以图3b示出的磁共振数据为例，一个扫描周期包括4个时间段，进行了3个扫描周期的扫描，即可得到3*4个，共12个欠采样的K空间。
需要说明的是，每个K空间中形成的径向线数量根据实际所需要的空间分辨率而定，不限于2条，可以是1条，3条等，本实施例中为形象表示只画出两条径向线。
步骤303、从填补K空间中获取第一预设区域内的采样数据点，并填补于目标K空间的数据缺失区域。
其中，目标K空间为多个K空间中的一个，填补K空间为不同扫描周期中与目标K空间的时间段对应的时间段的K空间。以图3b示出的磁共振数据为例，若目标K空间为第1个扫描周期中第1时间段的磁共振数据映射得到的，填补K空间即为第2个扫描周期和第3个扫描周期中第1时间段的磁共振数据映射得到的K空间；若目标K空间为第2个扫描周期中第2时间段的磁共振数据映射得到的，填补K空间即为第1个扫描周期和第3个扫描周期中第2时间段的磁共振数据映射得到的K空间；以此类推。
本实施例中，第一预设区域包含K空间的边缘区域。根据不同扫描周期对应的时间段所采集的磁共振数据得到的K空间具有类似的高频信息，反映在K空间的边缘区域。由此步骤303中，可获取目标K空间的相邻周期同一时间段K空间的边缘区域的数据采样点，填补于目标K空间的欠采边缘区域。参见图3b，图中填充图案相同的六边形表征的磁共振数据得到的3个K空间，具有相似的高频信息。针对不同的K空间，预设区域可以设置成相同，也可以设置成不同。
参见图3c，图中示出了使用第1个扫描周期的第1时间段的K空间(填补K空间)和第3个扫描周期的第1时间段的K空间(填补K空间)对第2个扫描周期的第1时间段的K空间(目标K空间)进行数据填补的结果。从图中可以看出，中心区域和中间区域采用的是目标K空间的采样数据，保留了目标K空间的对比度信息，目标K空间中的边缘区域的数据缺失部分采用填补K空间的边缘区域的全部或部分数据进行填补，从而实现了在不牺牲空间分辨率和时间分辨率的同时，更精确地对K空间的数据缺失区域进行数据填补。
步骤304、基于经过数据填补的目标K空间重建针对目标对象的磁共振图像。
具体的，对步骤304获得的K空间进行傅里叶变换，即可重建出针对目标对象的磁共振图像。
从而，将步骤302获得的每个K空间依次选取为目标K空间，并进行步骤303和步骤304的数据填补及图像重建，既能得到多幅磁共振图像，实现动态增强磁共振成像，且每磁共振图具有高时间分辨率和空间分辨率，去除了因呼吸运动出现的伪影。
本实施例中，借助相邻扫描周期中与目标K空间的呼吸状态相同的时间段的K空间的边缘区域数据对目标K空间进行缺失数据填补，该过程使得目标K空间所对应的高频信息增多了，更多的K空间高频真实信息被保留，在不牺牲空间分辨率和时间分辨率的同时，充分利用了k空间的多维时间信息，更精确地对K空间的数据缺失区域进行数据填补，进而重建出的图像精确度大大提高。
图4a是本发明一示例性实施例示出的另一种磁共振成像方法的流程图，本实施例与图3a示出的磁共振成像方法基本相同，不同之处在于，本实施例中还从与目标K空间的扫描周期相同的其他K空间获取采样数据点，以进一步增加目标K空间的信息。参见图4a，方法包括以下步骤：
步骤401、在对目标对象进行至少两个扫描周期的扫描过程中，基于径向轨迹采集针对目标对象的磁共振数据。
步骤402、针对每个扫描周期，将各自时间段的磁共振数据映射于对应的K空间内。
步骤403、从填补K空间中获取第一预设区域内的采样数据点，并填补于目标K空间的数据缺失区域。
步骤404、从其他K空间中获取第二预设区域内的采样数据点，并填补于目标K空间的数据缺失区域。
其中，其他K空间为目标K空间的扫描周期中除目标K空间之外的K空间。以图3b示出的磁共振数据为例，若目标K空间为第1个扫描周期中第1时间段的磁共振数据映射得到的，其他K空间即为第1个扫描周期中第2时间段、第3时间段和第4时间段的磁共振数据映射得到的K空间；若目标K空间为第2个扫描周期中第2时间段的磁共振数据映射得到的，其他K空间即为第2个扫描周期中第1时间段、第3时间段和第4时间段的磁共振数据映射得到的K空间；以此类推。
本实施例中，第二预设区域包含K空间的中间区域和边缘区域。步骤404中，也即获取同一周期时间段K空间的中间区域和边缘区域的数据采样点并填补于目标K空间的数据缺失区域。同一周期不同时间段采集的数据具有相似的对比度信息，同一周期相邻的时间段采集的数据具有相似的高频信息，因此，可以从同一周期时间段K空间的中间区域和边缘区域中获取数据采样点并对目标K空间进行数据填补。针对不同的K空间，预设区域可以设置成相同，也可以设置成不同。
参见图4b，图4b为在图3c的基础上，使用第2个扫描周期的第2时间段～第4时间段的K空间(其他K空间)对目标K空间进一步进行填补的结果。从图中可以看出，中心区域采用的是目标K空间的采样数据，保留了目标K空间的对比度信息；目标K空间中的中间区域的数据缺失部分，采用同一周期其他时间段的K空间的中间区域的全部或部分数据进行填补；目标K空间中的边缘区域的数据缺失部分，采用同一周期其他时间段以及相邻周期同一时间段的K空间的边缘区域的全部或部分数据以及填补K空间的边缘区域的全部或部分数据进行填补。从而，实现了对图3c示出的K空间的数据缺失区域进行了进一步的数据填补。
步骤405、基于经过数据填补的目标K空间重建针对目标对象的磁共振图像。
本实施例的步骤401～步骤403、步骤405与步骤301～步骤304的具体实现方式类似，此处不再赘述。
需要说明的是，步骤404不限于在步骤403之后执行，步骤404可以在步骤403之前执行，步骤404也可以与步骤403同步执行。
在另一个实施例中，若目标K空间的中间区域的采样数据点足够多，能够重建出质量较好的图像，步骤404中，也可从其他K空间中仅获取第一预设区域(边缘区域)内的采样数据点，并填补于目标K空间的数据缺失区域。
图5a是本发明一示例性实施例示出的另一种磁共振成像方法的流程图，本实施例与图4a示出的磁共振成像方法基本相同，不同之处在于，本实施例中从目标K空间的扫描周期中与目标K空间的时间段相邻的K空间获取采样数据点。参见图5a，方法包括以下步骤：
步骤501、在对目标对象进行至少两个扫描周期的扫描过程中，基于径向轨迹采集针对目标对象的磁共振数据。
步骤502、针对每个扫描周期，将相同时间段的磁共振数据映射于一个K空间内。
步骤503、从填补K空间中获取第一预设区域内的采样数据点，并填补于目标K空间的数据缺失区域。
步骤504、从相邻K空间中获取第二预设区域内的采样数据点，并填补于目标K空间的数据缺失区域。
其中，相邻K空间为目标K空间的扫描周期中与目标K空间的时间段相邻的K空间。以图3b示出的磁共振数据为例，若目标K空间为第1个扫描周期中第1时间段的磁共振数据映射得到的，相邻K空间即为第1个扫描周期中第2时间段的磁共振数据映射得到的K空间；若目标K空间为第2个扫描周期中第2时间段的磁共振数据映射得到的，相邻K空间即为第2个扫描周期中第1时间段和第3时间段的磁共振数据映射得到的K空间；以此类推。
可以理解的，相同扫描周期中，与目标K空间的时间段越相近的时间段对应的K空间的对比度信息越相似，也即采样时间段越接近，数据相关性越好，因此可以选择与目标K空间的时间段相邻的1个或2个相邻K空间对目标K空间进行数据填补。本实施例中，在目标K空间的数据足够用于图像重建的情况下，仅从目标K空间的相邻K空间中获取采样数据点，以避免K空间边缘区域中较大的时间窗口会导致高空间频率数据的平均。
参见图5b，图5b为在图3c的基础上，使用第2个扫描周期的第2时间段的K空间(相邻K空间)对目标K空间进一步进行填补的结果。从图中可以看出，中心区域采用的是目标K空间的采样数据，保留了目标K空间的对比度信息；目标K空间中的中间区域的数据缺失部分，采用同一周期相邻K空间的中间区域的全部或部分数据进行填补；目标K空间中的边缘区域的数据缺失部分，采用相邻周期同一时间段的K空间的边缘区域的全部或部分数据以及填补K空间的边缘区域的全部或部分数据进行填补。
步骤505、基于经过数据填补的目标K空间重建针对目标对象的磁共振图像。
本实施例的步骤501～步骤503、步骤505与步骤401～步骤403、步骤405的具体实现方式类似，此处不再赘述。
本实施例中，不仅利用同一扫描周期相邻时间段的K空间，还利用了相邻扫描周期同一呼吸状态的K空间对目标K空间进行数据填补，从而一方面增加了目标K空间的高低频信息，另一方面减小了K空间的边缘区域的时间窗口，减少了高频信息重建的不准确性，提高了高频信息重建的准确性。
与前述磁共振成像方法实施例相对应，本发明还提供了磁共振成像装置的实施例。
图6是本发明一示例性实施例示出的一种磁共振成像装置的模块示意图，该磁共振成像装置包括：采集模块61、映射模块62、填补模块63和重建模块64。
采集模块61用于在对目标对象进行至少两个扫描周期的扫描过程中，基于径向轨迹采集针对所述目标对象的磁共振数据，其中，每个扫描周期包含多个时间段，在所述扫描过程中，对于不同扫描周期中对应的时间段，所述目标对象处于相同的呼吸状态；
映射模块62用于针对每个扫描周期，将相同时间段的磁共振数据映射于一个K空间内，在每个K空间中形成至少一条径向线，所述径向线包含多个表征频率的采样数据点；
填补模块63用于从填补K空间中获取第一预设区域内的采样数据点，并填补于所述目标K空间的数据缺失区域，其中，所述目标K空间为多个K空间中的一个，所述填补K空间为不同扫描周期中与所述目标K空间的时间段对应的时间段的K空间；
重建模块64用于基于经过数据填补的目标K空间重建针对所述目标对象的磁共振图像。
可选地，所述填补模块还用于：
从其他K空间中获取第二预设区域内的采样数据点，并填补于所述目标K空间的数据缺失区域，其中，所述其他K空间为所述目标K空间的扫描周期中除所述目标K空间之外的K空间。
可选地，所述填补模块还用于：
从相邻K空间中获取第二预设区域内的采样数据点，并填补于所述目标K空间的数据缺失区域，其中，所述相邻K空间为所述目标K空间的扫描周期中与所述目标K空间的时间段相邻的K空间。
可选地，所述K空间包含中心区域、环绕所述中心区域的中间区域以及环绕所述中间区域的边缘区域；
所述第一预设区域包含所述边缘区域；
所述第二区域包含所述中间区域和所述边缘区域。
可选地，所述采集模块具体用于：
基于呼吸门控技术采集所述磁共振数据。
对于装置实施例而言，由于其基本对应于方法实施例，所以相关之处参见方法实施例的部分说明即可。以上所描述的装置实施例仅仅是示意性的，其中所述作为分离部件说明的单元可以是或者也可以不是物理上分开的，作为单元显示的部件可以是或者也可以不是物理单元，即可以位于一个地方，或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部模块来实现本发明方案的目的。本领域普通技术人员在不付出创造性劳动的情况下，即可以理解并实施。
图7是本发明一示例实施例示出的一种电子设备的结构示意图，示出了适于用来实现本发明实施方式的示例性电子设备70的框图。图7显示的电子设备70仅仅是一个示例，不应对本发明实施例的功能和使用范围带来任何限制。
如图7所示，电子设备70可以以通用计算设备的形式表现，例如其可以为服务器设备。电子设备70的组件可以包括但不限于：上述至少一个处理器71、上述至少一个存储器72、连接不同系统组件(包括存储器72和处理器71)的总线73。
总线73包括数据总线、地址总线和控制总线。
存储器72可以包括易失性存储器，例如随机存取存储器(RAM)721和/或高速缓存存储器722，还可以进一步包括只读存储器(ROM)723。
存储器72还可以包括具有一组(至少一个)程序模块724的程序工具725(或实用工具)，这样的程序模块724包括但不限于：操作系统、一个或者多个应用程序、其它程序模块以及程序数据，这些示例中的每一个或某种组合中可能包括网络环境的实现。
处理器71通过运行存储在存储器72中的计算机程序，从而执行各种功能应用以及数据处理，例如上述任一实施例所提供的方法。
电子设备70也可以与一个或多个外部设备74(例如键盘、指向设备等)通信。这种通信可以通过输入/输出(I/O)接口75进行。并且，模型生成的电子设备70还可以通过网络适配器76与一个或者多个网络(例如局域网(LAN)，广域网(WAN)和/或公共网络，例如因特网)通信。如图所示，网络适配器76通过总线73与模型生成的电子设备70的其它模块通信。应当明白，尽管图中未示出，可以结合模型生成的电子设备70使用其它硬件和/或软件模块，包括但不限于：微代码、设备驱动器、冗余处理器、外部磁盘驱动阵列、RAID(磁盘阵列)系统、磁带驱动器以及数据备份存储系统等。
应当注意，尽管在上文详细描述中提及了电子设备的若干单元/模块或子单元/模块，但是这种划分仅仅是示例性的并非强制性的。实际上，根据本发明的实施方式，上文描述的两个或更多单元/模块的特征和功能可以在一个单元/模块中具体化。反之，上文描述的一个单元/模块的特征和功能可以进一步划分为由多个单元/模块来具体化。
本发明实施例还提供一种计算机可读存储介质，其上存储有计算机程序，其特征在于，所述计算机程序被处理器执行时实现上述任一项实施例提供的方法的步骤。
以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本发明保护的范围之内。

Claims (12)

1.一种磁共振成像方法，其特征在于，所述磁共振成像方法包括：
在对目标对象进行至少两个扫描周期的扫描过程中，基于径向轨迹采集针对所述目标对象的磁共振数据，其中，每个扫描周期包含多个时间段，对于不同扫描周期中对应的时间段，所述目标对象在所述扫描过程中处于相同的呼吸状态；
针对每个扫描周期，将相同时间段的磁共振数据映射于一个K空间内，在每个K空间中形成至少一条径向线，所述径向线包含多个表征频率的采样数据点；
从填补K空间中获取第一预设区域内的采样数据点，并填补于目标K空间的数据缺失区域，其中，所述目标K空间为多个K空间中的一个，所述填补K空间为不同扫描周期中与所述目标K空间的时间段对应的时间段的K空间；
基于经过数据填补的目标K空间重建针对所述目标对象的磁共振图像。
2.如权利要求1所述的磁共振成像方法，其特征在于，所述磁共振成像方法还包括：
从其他K空间中获取第二预设区域内的采样数据点，并填补于所述目标K空间的数据缺失区域，其中，所述其他K空间为所述目标K空间的扫描周期中除所述目标K空间之外的K空间。
3.如权利要求1所述的磁共振成像方法，其特征在于，所述磁共振成像方法还包括：
从相邻K空间中获取第二预设区域内的采样数据点，并填补于所述目标K空间的数据缺失区域，其中，所述相邻K空间为所述目标K空间的扫描周期中与所述目标K空间的时间段相邻的K空间。
4.如权利要求2或3所述的磁共振成像方法，其特征在于，所述K空间包含中心区域、环绕所述中心区域的中间区域以及环绕所述中间区域的边缘区域；
所述第一预设区域包含所述边缘区域；
所述第二预设区域包含所述中间区域和所述边缘区域。
5.如权利要求1所述的磁共振成像方法，其特征在于，采集针对所述目标对象的磁共振数据，包括：
基于呼吸门控技术采集所述磁共振数据。
6.一种磁共振成像装置，其特征在于，所述磁共振成像装置包括：
采集模块，用于在对目标对象进行至少两个扫描周期的扫描过程中，基于径向轨迹采集针对所述目标对象的磁共振数据，其中，每个扫描周期包含多个时间段，在所述扫描过程中，对于不同扫描周期中对应的时间段，所述目标对象处于相同的呼吸状态；
映射模块，用于针对每个扫描周期，将相同时间段的磁共振数据映射于一个K空间内，在每个K空间中形成至少一条径向线，所述径向线包含多个表征频率的采样数据点；
填补模块，用于从填补K空间中获取第一预设区域内的采样数据点，并填补于目标K空间的数据缺失区域，其中，所述目标K空间为多个K空间中的一个，所述填补K空间为不同扫描周期中与所述目标K空间的时间段对应的时间段的K空间；
重建模块，用于基于经过数据填补的目标K空间重建针对所述目标对象的磁共振图像。
7.如权利要求6所述的磁共振成像装置，其特征在于，所述填补模块还用于：
从其他K空间中获取第二预设区域内的采样数据点，并填补于所述目标K空间的数据缺失区域，其中，所述其他K空间为所述目标K空间的扫描周期中除所述目标K空间之外的K空间。
8.如权利要求6所述的磁共振成像装置，其特征在于，所述填补模块还用于：
从相邻K空间中获取第二预设区域内的采样数据点，并填补于所述目标K空间的数据缺失区域，其中，所述相邻K空间为所述目标K空间的扫描周期中与所述目标K空间的时间段相邻的K空间。
9.如权利要求7或8所述的磁共振成像装置，其特征在于，所述K空间包含中心区域、环绕所述中心区域的中间区域以及环绕所述中间区域的边缘区域；
所述第一预设区域包含所述边缘区域；
所述第二预设区域包含所述中间区域和所述边缘区域。
10.如权利要求6所述的磁共振成像装置，其特征在于，所述采集模块具体用于：
基于呼吸门控技术采集所述磁共振数据。
11.一种电子设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，所述处理器执行所述计算机程序时实现权利要求1至5中任一项所述的磁共振成像方法。
12.一种计算机可读存储介质，其上存储有计算机程序，其特征在于，所述计算机程序被处理器执行时实现权利要求1至5中任一项所述的磁共振成像方法的步骤。
CN202010075000.3A 2020-01-22 2020-01-22 磁共振成像方法、装置、电子设备、存储介质 Active CN111273207B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010075000.3A CN111273207B (zh) 2020-01-22 2020-01-22 磁共振成像方法、装置、电子设备、存储介质

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010075000.3A CN111273207B (zh) 2020-01-22 2020-01-22 磁共振成像方法、装置、电子设备、存储介质

Publications (2)

Publication Number Publication Date
CN111273207A CN111273207A (zh) 2020-06-12
CN111273207B true CN111273207B (zh) 2022-06-10

Family

ID=70996944

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010075000.3A Active CN111273207B (zh) 2020-01-22 2020-01-22 磁共振成像方法、装置、电子设备、存储介质

Country Status (1)

Country Link
CN (1) CN111273207B (zh)

Family Cites Families (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP3699304B2 (ja) * 1999-08-13 2005-09-28 ジーイー横河メディカルシステム株式会社 磁気共鳴撮像装置
US7609058B2 (en) * 2006-11-17 2009-10-27 Siemens Medical Solutions Usa, Inc. Method and apparatus for generating a magnetic resonance data file
US8089278B1 (en) * 2008-04-11 2012-01-03 The Regents Of The University Of California Time-resolved contrast-enhanced magnetic resonance (MR) angiography
JP5591493B2 (ja) * 2008-07-17 2014-09-17 株式会社東芝 磁気共鳴イメージング装置
US20100145185A1 (en) * 2008-12-04 2010-06-10 Xiaole Hong Extending the resolution of mri data by combining subsets from plural image acquisitions
DE102009007361B4 (de) * 2009-02-04 2011-04-07 Siemens Aktiengesellschaft Aufnahme von zeitaufgelösten Magnetresonanzdaten
CN106308798B (zh) * 2013-02-04 2019-05-31 上海联影医疗科技有限公司 一种磁共振扫描短te成像方法及磁共振扫描系统
DE102018203507B4 (de) * 2018-03-08 2020-10-22 Siemens Healthcare Gmbh Verfahren zur Aufnahme und Rekonstruktion eines vierdimensionalen Dynamikbilddatensatzes mit einer Magnetresonanzeinrichtung, Magnetresonanzeinrichtung, Computerprogramm und elektronisch lesbarer Datenträger
CN110338795B (zh) * 2019-07-11 2023-09-01 中国科学院深圳先进技术研究院 一种径向黄金角磁共振心脏电影成像方法、装置和设备

Also Published As

Publication number Publication date
CN111273207A (zh) 2020-06-12

Similar Documents

Publication Publication Date Title
CN111513716B (zh) 使用扩展灵敏度模型和深度神经网络进行磁共振图像重建的方法和系统
JP6998218B2 (ja) 動き検出を用いるｍｒ撮像
JP4951264B2 (ja) シェルｋ空間サンプリング軌跡を使用したアンダーサンプリング３ＤＭＲＩ
US7358730B2 (en) Diffusion tensor imaging using highly constrained image reconstruction method
US9482732B2 (en) MRI reconstruction with motion-dependent regularization
CN103477238B (zh) 采用来自先验采集的约束的压缩感测mr图像重建
WO2020019412A1 (zh) 磁共振成像方法、装置、设备及存储介质
US20140086469A1 (en) Mri reconstruction with incoherent sampling and redundant haar wavelets
US8022700B2 (en) Method and apparatus for view ordering of magnetic resonance imaging data for dynamic studies
WO2007124444A2 (en) Method and system for parallel reconstruction in the k-space domain for application in imaging systems
JP3556176B2 (ja) Ｍｒ投影画像の自動修正
RU2707661C2 (ru) Мр (магнитно-резонансная) томография методом propeller c подавлением артефактов
WO2014154544A1 (en) Real-time motion correction for mri using fat navigators
CN112370040B (zh) 磁共振成像方法、装置、存储介质及电子设备
US20140152303A1 (en) Magnetic resonance imaging data sampling methods and systems
JP2004000593A (ja) 画像再構成するための方法及び装置
CN114913255A (zh) 用于使用磁共振图像中的高分辨率相位重建的系统和方法
CN110542872B (zh) 磁共振成像方法及设备、相位信息的获取方法及装置
CN111273207B (zh) 磁共振成像方法、装置、电子设备、存储介质
CN111292388A (zh) 磁共振成像方法、装置、电子设备、存储介质
CN113466768B (zh) 磁共振成像方法及磁共振成像系统
JP7060685B2 (ja) 体動依存型ラジアル又はスパイラルｋ空間サンプリングを使用するＭＲイメージング
CN115552272A (zh) 使用多个磁共振成像系统配置对磁共振图像的校正
CN111948590B (zh) 磁共振成像方法及装置、电子设备、存储介质
US11474183B1 (en) Motion correction systems and methods of propeller magnetic resonance images

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