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

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

Info

Publication number
CN111292388A
CN111292388A CN202010075253.0A CN202010075253A CN111292388A CN 111292388 A CN111292388 A CN 111292388A CN 202010075253 A CN202010075253 A CN 202010075253A CN 111292388 A CN111292388 A CN 111292388A
Authority
CN
China
Prior art keywords
space
data points
magnetic resonance
spaces
target
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.)
Pending
Application number
CN202010075253.0A
Other languages
English (en)
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
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 filed Critical Neusoft Medical Systems Co Ltd
Priority to CN202010075253.0A priority Critical patent/CN111292388A/zh
Publication of CN111292388A publication Critical patent/CN111292388A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T11/00—2D [Two Dimensional] image generation
    • G06T11/003—Reconstruction from projections, e.g. tomography
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2210/00—Indexing scheme for image generation or computer graphics
    • G06T2210/41—Medical

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Magnetic Resonance Imaging Apparatus (AREA)

Abstract

本发明公开了磁共振成像方法、装置、电子设备、存储介质，方法包括：在一个扫描周期内，基于径向轨迹采集针对目标对象的磁共振数据，扫描周期包含多个时间段；将相同时间段的磁共振数据映射于相同的K空间内，在每个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空间的时间段相邻的1个或2个其他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个或2个其他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空间进行加权处理的结果示意图；
图3是本发明一示例性实施例示出的一种磁共振成像方法的流程图；
图4a是对图2a中的采样数据点进行插值处理后的结果示意图；
图4b是对图2b中的采样数据点进行插值处理后的结果示意图；
图4c是对图2c中的采样数据点进行插值处理后的结果示意图；
图4d是对图2d中的采样数据点进行插值处理后的结果示意图；
图5是本发明一示例性实施例示出的使用图2b～2d的K空间的部分采样数据点对图4a的K空间进行数据填补的结果示意图；
图6是本发明一示例性实施例示出的使用图4b的K空间的采样数据点和插值数据点对图4a的K空间进行数据填补的结果示意图。
图7是本发明一示例性实施例示出的另一种磁共振成像方法的流程图；
图8是本发明一示例性实施例示出的一种磁共振成像装置的模块示意图；
图9是本发明一示例实施例示出的一种电子设备的结构示意图。
具体实施方式
这里将详细地对示例性实施例进行说明，其示例表示在附图中。下面的描述涉及附图时，除非另有表示，不同附图中的相同数字表示相同或相似的要素。以下示例性实施例中所描述的实施方式并不代表与本发明相一致的所有实施方式。相反，它们仅是与如所附权利要求书中所详述的、本发明的一些方面相一致的装置和方法的例子。
在本发明使用的术语是仅仅出于描述特定实施例的目的，而非旨在限制本发明。在本发明和所附权利要求书中所使用的单数形式的“一种”、“所述”和“该”也旨在包括多数形式，除非上下文清楚地表示其他含义。还应当理解，本文中使用的术语“和/或”是指并包含一个或多个相关联的列出项目的任何或所有可能组合。
应当理解，尽管在本发明可能采用术语第一、第二、第三等来描述各种信息，但这些信息不应限于这些术语。这些术语仅用来将同一类型的信息彼此区分开。例如，在不脱离本发明范围的情况下，第一信息也可以被称为第二信息，类似地，第二信息也可以被称为第一信息。取决于语境，如在此所使用的词语“如果”可以被解释成为“在……时”或“当……时”或“响应于确定”。
图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空间中的空白区域(没有径向线的区域)表示未被采样的区域，也即数据缺失区域。
在图像重建时，一个K空间中的采样数据较少，为了重建出分辨率高的图像，需要对K空间中的数据缺失区域(空白区域)进行数据填补。相关技术中，将一个扫描周期内获得的所有K空间进行加权处理，得到加权后的K空间，进而进行图像重建。以对图2a～图2d示出的K空间进行加权处理为例，得到的结果参见图2e所示。该方式的缺点是，K空间的中心区域A的信息是不同时间段采集的数据的加权结果，时间间隔越大，由于被检测对象的运动(肌肉收缩、呼吸运动等)带来的信息差异也越大，致使K空间的高频信息的加权结果会存在较大误差，影响最终的图像重建精度。
在K空间中，每条径向线含有低频到高频信息。越靠近K空间的中心区域的数据采样点具有越低的空间频率值；越靠近K空间的边缘区域的数据采样点具有越高的空间频率值。一个扫描周期内，所获得的每个K空间可利用前后时间段的时间分辨率信息，实现在不牺牲空间分辨率的同时，更精确地对K空间的数据缺失区域进行数据填补。基于此，本发明提供了一种磁共振成像方法，以提高对基于径向轨迹采集磁共振数据进行图像重建的精确度。
以下结合图1所示的MRI系统以及图2a～图2d，以一个扫描层面为例，对本发明的磁共振成像方法实施例进行详细描述。
图3是本发明一示例性实施例示出的一种磁共振成像方法的流程图，方法包括以下步骤：
步骤301、在一个扫描周期内，基于径向轨迹采集针对目标对象的磁共振数据，并将磁共振数据映射于多个K空间。
若将扫描周期划分为n个时间段，将该扫描周期内对一个扫描层面进行扫描得到的扫描数据，按照相同时间段进行K空间填补，即可得到n个K空间。以一个扫描周期包括4个时间段为例，即可得到4个K空间，具体可参见图2a～图2d所示。需要说明的是，每个K空间中形成的径向线数量根据实际所需要的空间分辨率而定，不限于2条，可以是1条，3条等，本实施例中为形象表示只画出两条径向线。
步骤302、基于并行成像算法对目标K空间中的采样数据点进行插值处理，得到至少一条由多个插值数据点形成的平行于径向线的填补线。
其中，目标K空间为步骤301获得的K空间中的一个。
步骤302中，对目标K空间中的采样数据点进行插值处理的目的是为了填补由于沿时间维度的欠采样导致的一些缺失数据。参见图4a～4d，图中分别示出了对图2a～2d的K空间中的采样数据点进行插值处理的结果，图中的实线表示径向线(实际采样数据)，虚线表示由插值数据点形成的填补线(插值结果)。与采样数据点类似的，插值数据点由kx和ky表示。
经过插值处理，将每条径向线扩展为由m条相互平行的K空间线(包括1条径向线和m-1条填补线)组成的K空间区段。经过插值处理，目标K空间中对应的高低频信息增多了。
其中，填补线的条数可根据实际需求自行确定，可以理解的，填补线的条数越多，图像重建的计算量越大，且离径向线的距离越远，与径向线的相关性越小，因此没有必要插值出很多的填补线。根据试验可知，在径向线的两侧插值出共4条填补线就能达到比较好的建像效果。
径向线与填补线之间的间距(或填补线与填补线之间的间距)与建像视野FOV有关，间距Δk与建像视野FOV的关系表示如下：
Figure BDA0002378342650000091
步骤303、从其他K空间获取预设区域内的采样数据点，并填补于经过插值处理的目标K空间的数据缺失区域。
其中，其他K空间为步骤301获得的K空间中除目标K空间之外的K空间。
K空间边缘区域的数据含有更多的不会随时间变化而变化的信息，因此可以获取其他K空间的边缘区域的全部或部分数据对目标K空间的空白区域进行数据填补。预设区域也即K空间的边缘区域，预设区域中的采样数据点的频率较低，针对不同的K空间，预设区域可以设置成相同，也可以设置成不同。
参见图5，图5示出了使用图2b～2d的K空间(其他K空间)的部分采样数据点对图4a的K空间(经过插值处理后的目标K空间)进行填补的结果。从图中可以看出，中心区域B使用的是对目标K空间的采样数据插值得到的数据，保留了目标K空间的时间变化信息，目标K空间的其他数据缺失区域采用其他K空间边缘区域的全部或部分数据进行填补，从而实现了在不牺牲空间分辨率和时间分辨率的同时，更精确地对K空间的数据缺失区域进行数据填补。
可以理解的，采样时间段越接近，数据相关性越好，因此可以选择与目标K空间的时间段相邻的1个或2个其他K空间对目标K空间进行数据填补。假设图2a～2d分别为根据相邻时间段采集的磁共振数据获得的K空间，若选择图2b作为目标图像，则可从图2a和图2c中获取采样数据点填补于图2b中。
步骤304、基于经过数据缺失点填补的目标K空间重建针对目标对象的磁共振图像。
具体的，对步骤303获得的K空间进行傅里叶变换，即可重建出针对目标对象的磁共振图像。
本实施例中，基于并行成像算法对目标K空间中的采样数据点进行插值处理后，目标K空间所对应的高低频信息增多了，又借助相同扫描周期的其他K空间的边缘区域数据对目标K空间进行缺失数据填补，该过程使得更多的K空间高低频真实信息被保留，在不牺牲空间分辨率和时间分辨率的同时，更精确地对K空间的数据缺失区域进行数据填补，进而重建出的图像精确度大大提高。
在另一个实施例中，步骤303之前，可以基于并行成像算法对其他K空间的采样数据点进行插值处理，具体实现过程与对目标K空间的实现过程类似，此处不再赘述。步骤303中，则从其他K空间获取预设区域内的采样数据点和插值数据点，并填补于经过插值处理的目标K空间中的数据缺失区域。参见图6，图6示出了使用图4b的K空间(其他K空间)的采样数据点和插值数据点对图4a的K空间(目标K空间)进行数据填补的结果。
其他K空间经过数据插值处理后，对应的高低频信息有所增多。从而可仅使用与目标K空间相邻的时间段的其他K空间对目标K空间进行数据填补，即可使目标K空间具有足够多的数据，重建出分辨率高的图像。
图7是本发明一示例性实施例示出的一种磁共振成像方法的流程图，方法包括以下步骤：
步骤701、在一个扫描周期内，基于径向轨迹采集针对目标对象的磁共振数据，并将磁共振数据映射于多个K空间。
步骤702、基于并行成像算法对每个K空间的采样数据点进行插值处理，得到由多个插值数据点形成的至少一条平行于径向线的填补线。
其中，步骤701和步骤702的具体实现方式与步骤301和步骤302类似此处不再赘述。
步骤703、从其他K空间获取预设区域内的采样数据点和插值数据点，并填补于经过插值处理的目标K空间的数据缺失区域。
其中，目标K空间为步骤702中经过插值处理的多个K空间中的一个，其他K空间为除目标K空间之外的其他K空间。
在另一个实施例中，也可以只选择与目标K空间对应的时间段相邻的其他K空间中获取采样数据点和插值数据点。
步骤704、基于经过数据填补的目标K空间重建针对所述目标对象的磁共振图像。
从而，从步骤702获得的多个K空间中依次选择一个作为目标K空间,即可重建出扫描周期内获得的所有K空间的磁共振图像。
与前述磁共振成像方法实施例相对应，本发明还提供了磁共振成像装置的实施例。
图8是本发明一示例性实施例示出的一种磁共振成像装置的模块示意图，装置包括：采集模块81、映射模块82、插值模块83、填补模块84和重建模块85。
采集模块81用于在一个扫描周期内，基于径向轨迹采集针对目标对象的磁共振数据，所述扫描周期包含多个时间段；
映射模块82用于将相同时间段的磁共振数据映射于相同的K空间内，在每个K空间中形成至少一条径向线，所述径向线包含多个表征频率的采样数据点；
插值模块83用于基于并行成像算法对目标K空间中的采样数据点进行插值处理，得到多个插值数据点，所述多个插值数据点形成多条平行于所述径向线的填补线，其中，所述目标K空间为多个K空间中的一个；
填补模块84用于从其他K空间获取预设区域内的采样数据点，并填补于经过插值处理的目标K空间的数据缺失区域，其中，所述其他K空间为所述多个K空间中除所述目标K空间之外的K空间；
重建模块85用于基于经过数据填补的目标K空间重建针对所述目标对象的磁共振图像。
可选地，所述插值模块还用于基于并行成像算法对所述其他K空间中的采样数据点进行插值处理，得到多个插值数据点；
所述填补模块具体用于从其他K空间获取预设区域内的采样数据点和插值数据点，并填补于经过插值处理的目标K空间中的数据缺失区域。
可选地，所述填补模块具体用于从与所述目标K空间的时间段相邻的1个或2个其他K空间中获取所述采样数据点。
可选地，所述预设区域为所述K空间的边缘区域。
本发明实施例还提供一种磁共振成像装置，所述磁共振成像装置包括：
采集模块，用于在一个扫描周期内，基于径向轨迹采集针对目标对象的磁共振数据，所述扫描周期包含多个时间段；
映射模块，用于将相同时间段的磁共振数据映射于相同的K空间内，在每个K空间中形成至少一条径向线，所述径向线包含多个表征频率的采样数据点；
插值模块，用于针对每个K空间，基于并行成像算法对每个K空间中的采样数据点进行插值处理，得到多个插值数据点，所述多个插值数据点形成多条平行于所述径向线的填补线；
填补模块，用于从其他K空间获取预设区域内的采样数据点和/或插值数据点，并填补于经过插值处理的目标K空间的数据缺失区域，其中，所述目标K空间为经过插值处理的多个K空间中的一个，所述其他K空间为经过插值处理的多个K空间中除所述目标K空间之外的K空间；
重建模块，用于基于经过数据填补的目标K空间重建针对所述目标对象的磁共振图像。
图9是本发明一示例实施例示出的一种电子设备的结构示意图，示出了适于用来实现本发明实施方式的示例性电子设备90的框图。图9显示的电子设备90仅仅是一个示例，不应对本发明实施例的功能和使用范围带来任何限制。
如图9所示，电子设备90可以以通用计算设备的形式表现，例如其可以为服务器设备。电子设备90的组件可以包括但不限于：上述至少一个处理器91、上述至少一个存储器92、连接不同系统组件(包括存储器92和处理器91)的总线93。
总线93包括数据总线、地址总线和控制总线。
存储器92可以包括易失性存储器，例如随机存取存储器(RAM)921和/或高速缓存存储器922，还可以进一步包括只读存储器(ROM)923。
存储器92还可以包括具有一组(至少一个)程序模块924的程序工具925(或实用工具)，这样的程序模块924包括但不限于：操作系统、一个或者多个应用程序、其它程序模块以及程序数据，这些示例中的每一个或某种组合中可能包括网络环境的实现。
处理器91通过运行存储在存储器92中的计算机程序，从而执行各种功能应用以及数据处理，例如上述任一实施例所提供的方法。
电子设备90也可以与一个或多个外部设备94(例如键盘、指向设备等)通信。这种通信可以通过输入/输出(I/O)接口95进行。并且，模型生成的电子设备90还可以通过网络适配器96与一个或者多个网络(例如局域网(LAN)，广域网(WAN)和/或公共网络，例如因特网)通信。如图所示，网络适配器96通过总线93与模型生成的电子设备90的其它模块通信。应当明白，尽管图中未示出，可以结合模型生成的电子设备90使用其它硬件和/或软件模块，包括但不限于：微代码、设备驱动器、冗余处理器、外部磁盘驱动阵列、RAID(磁盘阵列)系统、磁带驱动器以及数据备份存储系统等。
应当注意，尽管在上文详细描述中提及了电子设备的若干单元/模块或子单元/模块，但是这种划分仅仅是示例性的并非强制性的。实际上，根据本发明的实施方式，上文描述的两个或更多单元/模块的特征和功能可以在一个单元/模块中具体化。反之，上文描述的一个单元/模块的特征和功能可以进一步划分为由多个单元/模块来具体化。
以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本发明保护的范围之内。

Claims (12)

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

Priority Applications (1)

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

Applications Claiming Priority (1)

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

Publications (1)

Publication Number Publication Date
CN111292388A true CN111292388A (zh) 2020-06-16

Family

ID=71022396

Family Applications (1)

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

Country Status (1)

Country Link
CN (1) CN111292388A (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112967297A (zh) * 2021-03-23 2021-06-15 中国科学院深圳先进技术研究院 磁共振氧十七代谢成像方法、装置、存储介质及终端设备

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1973211A (zh) * 2004-05-14 2007-05-30 皇家飞利浦电子股份有限公司 涉及k－空间中心过度采样的非笛卡尔轨迹的对比度预备mri
US20130119985A1 (en) * 2010-07-22 2013-05-16 Koninklijke Philips Electronics N.V. Mri method for retrospective motion correction with interleaved radial acquisition
US20150187073A1 (en) * 2014-01-02 2015-07-02 Siemens Aktiengesellschaft Method, reconstruction device, and magnetic resonance apparatus for reconstructing magnetic resonance raw data
US20160003929A1 (en) * 2014-07-03 2016-01-07 Siemens Aktiengesellschaft Magnetic resonance imaging with different frequency encoding patterns
US20160274209A1 (en) * 2015-03-16 2016-09-22 Kabushiki Kaisha Toshiba Accelerated MRI Using Radial Strips and Undersampling of K-Space

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1973211A (zh) * 2004-05-14 2007-05-30 皇家飞利浦电子股份有限公司 涉及k－空间中心过度采样的非笛卡尔轨迹的对比度预备mri
US20130119985A1 (en) * 2010-07-22 2013-05-16 Koninklijke Philips Electronics N.V. Mri method for retrospective motion correction with interleaved radial acquisition
US20150187073A1 (en) * 2014-01-02 2015-07-02 Siemens Aktiengesellschaft Method, reconstruction device, and magnetic resonance apparatus for reconstructing magnetic resonance raw data
US20160003929A1 (en) * 2014-07-03 2016-01-07 Siemens Aktiengesellschaft Magnetic resonance imaging with different frequency encoding patterns
US20160274209A1 (en) * 2015-03-16 2016-09-22 Kabushiki Kaisha Toshiba Accelerated MRI Using Radial Strips and Undersampling of K-Space

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
LI FENG ET AL: "Rapid free-breathing dynamic contrast-enhanced MRI using motion-resolved compressed sensing" *
SAMSONOV AA ET AL: "Advance in locally constrained k-space parallel MRI" *
SHEREEN M. EI-METWALLU ET AL: "Rapid cariac MRI using random radial trajectories" *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112967297A (zh) * 2021-03-23 2021-06-15 中国科学院深圳先进技术研究院 磁共振氧十七代谢成像方法、装置、存储介质及终端设备
CN112967297B (zh) * 2021-03-23 2024-05-03 中国科学院深圳先进技术研究院 磁共振氧十七代谢成像方法、装置、存储介质及终端设备

Similar Documents

Publication Publication Date Title
JP5960163B2 (ja) 校正スキャン、コイル感度マップ及びナビゲータを使用する、剛体動き補償のためのパラレルｍｒｉ方法
JP6998218B2 (ja) 動き検出を用いるｍｒ撮像
Stehning et al. Free‐breathing whole‐heart coronary MRA with 3D radial SSFP and self‐navigated image reconstruction
JP4951264B2 (ja) シェルｋ空間サンプリング軌跡を使用したアンダーサンプリング３ＤＭＲＩ
US7840045B2 (en) Method and system for parallel reconstruction in the K-space domain for application in imaging systems
US20140126796A1 (en) Mri reconstruction with motion-dependent regularization
US20140086469A1 (en) Mri reconstruction with incoherent sampling and redundant haar wavelets
US20070156045A1 (en) Diffusion tensor imaging using highly constrained image reconstruction method
US8022700B2 (en) Method and apparatus for view ordering of magnetic resonance imaging data for dynamic studies
RU2707661C2 (ru) Мр (магнитно-резонансная) томография методом propeller c подавлением артефактов
CN112370040B (zh) 磁共振成像方法、装置、存储介质及电子设备
CN114167334A (zh) 磁共振图像的重建方法、装置及电子设备
CN110542872B (zh) 磁共振成像方法及设备、相位信息的获取方法及装置
CN114913255A (zh) 用于使用磁共振图像中的高分辨率相位重建的系统和方法
CN111292388A (zh) 磁共振成像方法、装置、电子设备、存储介质
JP4343317B2 (ja) 磁気共鳴イメージング装置
CN111273207B (zh) 磁共振成像方法、装置、电子设备、存储介质
CN113466768B (zh) 磁共振成像方法及磁共振成像系统
JP7060685B2 (ja) 体動依存型ラジアル又はスパイラルｋ空間サンプリングを使用するＭＲイメージング
CN111948590B (zh) 磁共振成像方法及装置、电子设备、存储介质
Lee et al. Unsupervised motion artifact correction of turbo spin‐echo MRI using deep image prior
CN117630782A (zh) 磁共振扫描控制方法、磁共振成像方法及存储介质
CN116930839A (zh) 多对比度定量磁共振成像采样方法、系统、介质、电子设备
Nozaki et al. MRI apparatus and method with B 0 corrections using reference and object average units

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