CN109975343B - 磁共振成像方法、***及存储介质 - Google Patents

磁共振成像方法、***及存储介质 Download PDF

Info

Publication number
CN109975343B
CN109975343B CN201910232085.9A CN201910232085A CN109975343B CN 109975343 B CN109975343 B CN 109975343B CN 201910232085 A CN201910232085 A CN 201910232085A CN 109975343 B CN109975343 B CN 109975343B
Authority
CN
China
Prior art keywords
space data
data
magnetic resonance
fov
generate
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
CN201910232085.9A
Other languages
English (en)
Other versions
CN109975343A (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 CN201910232085.9A priority Critical patent/CN109975343B/zh
Publication of CN109975343A publication Critical patent/CN109975343A/zh
Priority to US16/813,881 priority patent/US11592506B2/en
Application granted granted Critical
Publication of CN109975343B publication Critical patent/CN109975343B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N24/00—Investigating or analyzing materials by the use of nuclear magnetic resonance, electron paramagnetic resonance or other spin effects
    • G01N24/08—Investigating or analyzing materials by the use of nuclear magnetic resonance, electron paramagnetic resonance or other spin effects by using nuclear magnetic resonance
    • 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
    • G01R33/4824—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 using a non-Cartesian trajectory
    • 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/561—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution by reduction of the scanning time, i.e. fast acquiring systems, e.g. using echo-planar pulse sequences
    • 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/561—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution by reduction of the scanning time, i.e. fast acquiring systems, e.g. using echo-planar pulse sequences
    • G01R33/5619—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution by reduction of the scanning time, i.e. fast acquiring systems, e.g. using echo-planar pulse sequences by temporal sharing of data, e.g. keyhole, block regional interpolation scheme for k-Space [BRISK]
    • 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/58—Calibration of imaging systems, e.g. using test probes, Phantoms; Calibration objects or fiducial markers such as active or passive RF coils surrounding an MR active material
    • G01R33/583—Calibration of signal excitation or detection systems, e.g. for optimal RF excitation power or frequency
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00—Image analysis
    • G06T7/0002—Inspection of images, e.g. flaw detection
    • G06T7/0012—Biomedical image inspection
    • 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/4806—Functional imaging of brain activation
    • 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
    • G01R33/4824—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 using a non-Cartesian trajectory
    • G01R33/4826—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 using a non-Cartesian trajectory in three dimensions
    • 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/5601—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution involving use of a contrast agent for contrast manipulation, e.g. a paramagnetic, super-paramagnetic, ferromagnetic or hyperpolarised contrast agent
    • 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
    • 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/561—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution by reduction of the scanning time, i.e. fast acquiring systems, e.g. using echo-planar pulse sequences
    • G01R33/5611—Parallel magnetic resonance imaging, e.g. sensitivity encoding [SENSE], simultaneous acquisition of spatial harmonics [SMASH], unaliasing by Fourier encoding of the overlaps using the temporal dimension [UNFOLD], k-t-broad-use linear acquisition speed-up technique [k-t-BLAST], k-t-SENSE
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/10—Image acquisition modality
    • G06T2207/10072—Tomographic images
    • G06T2207/10088—Magnetic resonance imaging [MRI]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/30—Subject of image; Context of image processing
    • G06T2207/30004—Biomedical image processing
    • 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/12—Bounding box
    • 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
    • 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/52—Parallel processing

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Engineering & Computer Science (AREA)
  • Condensed Matter Physics & Semiconductors (AREA)
  • General Health & Medical Sciences (AREA)
  • Health & Medical Sciences (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Radiology & Medical Imaging (AREA)
  • Signal Processing (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Quality & Reliability (AREA)
  • Medical Informatics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Pathology (AREA)
  • Immunology (AREA)
  • Biochemistry (AREA)
  • Analytical Chemistry (AREA)
  • Chemical & Material Sciences (AREA)
  • Magnetic Resonance Imaging Apparatus (AREA)

Abstract

本发明实施例公开了一种磁共振成像方法、***及存储介质，该方法包括：获取受检者的第一K空间数据，第一K空间数据基于非笛卡尔采样方式获取；对所述第一K空间数据进行重栅格化以生成第二K空间数据；对所述第二K空间数据进行校正处理以生成第三K空间数据，以使所述第三K空间数据对应的FOV小于所述第二K空间数据对应的FOV；基于并行成像方法和/或压缩感知方法确定所述第三K空间数据对应的磁共振图像。解决了现有技术的对非笛卡尔坐标系下的磁共振数据进行图像重建的速度较慢的技术问题，以达到实时成像的技术效果。

Description

磁共振成像方法、***及存储介质
技术领域
本发明实施例涉及磁共振成像领域，尤其涉及一种磁共振成像方法、***及存储介质。
背景技术
磁共振成像(MRI)可以在无损伤、无电离辐射的情况下得到样品/组织内部的高对比度清晰图像，在各个领域，尤其是临床医学诊断中得到了广泛应用。
受传统奈奎斯特采样定理的限制，基于笛卡尔采样获得磁共振数据的时间较长，成像速度较慢，无法满足临床高端实时成像的要求。相反地，非笛卡尔采样具有很多优点，比如成像速度快、对运动和流动不敏感等。然而，由于非笛卡尔采样的磁共振数据没有落在笛卡尔网格点上，因此不能利用传统的快速傅立叶变换对其进行图像重建，而必须使用其他的图像重建算法，但现有技术所使用的图像重建算法普遍较为复杂，执行效率低，进而导致图像重建速度较慢。
因此，有必要提供一种磁共振成像方法，以提高对非笛卡尔采样所获得的磁共振数据的图像重建速度，从而实现实时成像。
发明内容
本发明实施例提供了一种磁共振成像方法、***及存储介质，解决了现有技术对非笛卡尔坐标系下的磁共振数据进行图像重建的速度较慢的技术问题。
第一方面，本发明实施例提供了一种磁共振成像方法，包括：
采用获取受检者的第一K空间数据，所述第一K空间数据基于非笛卡尔采样方式获取；
对所述第一K空间数据进行重栅格化以生成第二K空间数据；
对所述第二K空间数据进行校正处理以生成第三K空间数据，以使所述第三K空间数据对应的FOV小于所述第二K空间数据对应的FOV；
基于并行成像方法和/或压缩感知方法确定所述第三K空间数据对应的磁共振图像。
第二方面，本发明实施例还提供了一种磁共振成像***，包括：
MR扫描器，用于通过获取受检者的第一K空间数据，所述第一K空间数据基于非笛卡尔采样方式获取；
处理器，用于获取所述第一K空间数据，以及对所述第一K空间数据进行重栅格化以生成第二K空间数据；以及对所述第二K空间数据进行校正处理以生成第三K空间数据，以使所述第三K空间数据对应的FOV小于所述第二K空间数据对应的FOV；以及基于并行成像方法和/或压缩感知方法确定所述第三K空间数据对应的磁共振图像。
第三方面，本发明实施例还提供了一种包含计算机可执行指令的存储介质，所述计算机可执行指令在由计算机处理器执行时用于执行如第一方面所述的磁共振成像方法。
本发明实施例提供的磁共振成像方法的技术方案，包括采用获取受检者的第一K空间数据，第一K空间数据基于非笛卡尔采样方式获取；对第一K空间数据进行重栅格化以生成第二K空间数据，从而使第二K空间数据的各个数据位于笛卡尔坐标的坐标点上；对第二K空间数据进行校正处理以生成第三K空间数据，以使第三K空间数据对应的FOV小于第二K空间数据对应的FOV；基于并行成像方法和/或压缩感知方法确定第三K空间数据对应的磁共振图像。本发明实施例使非笛卡尔采样结合并行成像方法和/或压缩感知方法进行图像重建成为了可能，从而可以显著提高图像重建的速度，满足临床实时成像的需求。
附图说明
为了更清楚地说明本发明实施例中的技术方案，下面将对实施例描述中所需要使用的附图做一简单地介绍，显而易见地，下面描述中的附图是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1是本发明实施例一提供的磁共振成像方法的流程图；
图2A是本发明实施例一提供的螺旋扫描轨迹的示意图；
图2B是本发明实施例一提供的径向扫描轨迹的示意图；
图3A是本发明实施例一提供的第二K空间数据对应的FOV的示意图；
图3B是本发明实施例一提供的第三K空间数据对应的FOV的示意图；
图4是本发明实施例二提供的磁共振成像方法的流程图；
图5是本发明实施例二提供的权重参数确定方法的流程图；
图6是本发明实施例三提供的磁共振成像方法的流程图；
图7是本发明实施例四提供的磁共振成像装置的结构框图；
图8是本发明实施例五提供的磁共振成像***的结构框图；
图9是本发明实施例五提供的又一磁共振成像***的结构框图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚，以下将参照本发明实施例中的附图，通过实施方式清楚、完整地描述本发明的技术方案，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
实施例一
图1是本发明实施例一提供的磁共振成像方法的流程图。本实施例的技术方案适用于对基于非笛卡尔采样方式获取的磁共振数据进行处理，从而使处理后的磁共振数据能够通过并行成像和/或压缩感知方法进行图像重建的情况。该方法可以由本发明实施例提供的磁共振成像装置来执行，该装置可以采用软件和/或硬件的方式实现，并配置在处理器中应用。该方法具体包括如下步骤：
S101、获取受检者的第一K空间数据，第一K空间数据基于非笛卡尔采样方式获取。
磁共振数据的采样方式通常包括笛卡尔采样方式和非笛卡尔采样方式，前者采集的磁共振数据均匀地分布在笛卡尔坐标的整数格点上，但采集过程中每次相位编码的耗时都比较长，因此随着相位编码步数的增加会显著增加扫描时间；后者的编码梯度不再是线性的，而是随着采样时间连续变化，每条采样轨迹也不再是笛卡尔坐标系下横平竖直的直线，因此能够有效地缩短采样时间，而且该采样方式还具有固有的变密度采样特性(K空间不同区域采样点的分布密度不同)，以及可有效地填充K空间等优点。
常用的非笛卡尔采样方式有螺旋扫描轨迹(参见图2A)和径向扫描轨迹(参见图2B)。图2B示出的径向扫描轨迹在kx-ky平面之内采用径向采集(沿极坐标的径向方向，每条线都通过原点，但是角度不同)，在kz方向依然采用笛卡尔坐标采集，同时还结合了等分圆周采集技术和采用黄金分割角(Golden Angle)采集技术。需要说明的是，以上仅是对非笛卡尔采样方式的示例性说明，本实施例对非笛卡尔采样方式的具体实现方式不予限定。
对于径向扫描轨迹，由于K空间中心在扫描过程中被重复多次采样，因此其对运动伪影不敏感，并且在K空间严重欠采样的情况下仍能获得较好的重建效果，因此可以用于对成像速度要求较高的场合，比如功能磁共振成像，螺旋采样可以有效利用梯度***硬件，常用于实时快速成像。
S102、对第一K空间数据进行重栅格化以生成第二K空间数据。
由于第一K空间数据的数据点不在笛卡尔坐标点(即：整数格点)上，直接对其进行图像重建的运算量非常大，图像重建速度特别慢，实践中难以直接使用。而要提高图像重建速度以满足临床实时成像的要求，通常需要先将其转换到整数格点上。为此本实施例先对第一K空间数据进行重栅格化以生成第二K空间数据，从而使第二K空间数据的各个数据点位于笛卡尔坐标点上。
S103、对第二K空间数据进行校正处理以生成第三K空间数据，以使第三K空间数据对应的FOV小于第二K空间数据对应的FOV。
由于重栅格化后的第一K空间数据，即第二K空间数据的数据量通常比较大，直接对其进行图像处理的运算量仍然比较高，因此本实施例对第二K空间数据进行校正处理以生成第三K空间数据，以使第三K空间数据对应的FOV小于第二K空间数据对应的FOV。
可以理解的是，原始磁共振数据对应的原始MRI图像(原始FOV)通常包括位于中心的有效图像部和位于周围的无效图像部，无效图像部理论上是没有信号的(通常对应磁共振图像外周的黑色区域)，实际上可能有一些没有价值的信号。考虑到重栅格化会使MRI图像的FOV进一步增大，但是位于中间的有效图像部不会发生变化，只是会增大无效图像部。因此本实施例的校正处理优选为图像域的校正处理，即对FOV进行处理，使处理后的FOV面积小于第二K空间数据对应的FOV面积，因此第三K空间数据的数据信息量小于现有技术的图像重建对象(第二K空间数据)的数据信息量。
可选地，校正处理为去掉第二K空间数据所对应FOV(参见图3A)的部分无效图像部以生成矩形FOV(参见图3B)，并将该矩形FOV对应的K空间数据作为第三K空间数据。
由于正方形也是矩形，因此本实施例在对第二K空间数据对应的正方形FOV进行处理时，可以将该正方形FOV(参见图3A)变成小一点的正方形FOV，也可以将其裁剪成长方形FOV(参见图3B)，无论裁剪后的FOV的形状如何，该FOV均包括全部的有效图像部。需要说明的是，通过裁剪改变FOV的形状只是其中的一种方法，本实施例不限于此。优选地，该矩形FOV的面积小于第二K空间数据对应的FOV面积，那么其对应的第三K空间数据的数据信息量小于第二K空间数据的数据信息量，即实现了通过裁剪FOV的无效图像部减少了最终参与图像重建的数据信息量。
S104、基于并行成像方法和/或压缩感知方法确定第三K空间数据对应的磁共振图像。
得到第三K空间数据后，采用并行成像方法和/或压缩感知方法对第三K空间数据进行图像重建以生成磁共振图像。
本发明实施例提供的磁共振成像方法的技术方案，包括采用获取受检者的第一K空间数据，第一K空间数据基于非笛卡尔采样方式获取；对第一K空间数据进行重栅格化以生成第二K空间数据，从而使第二K空间数据的各个数据位于笛卡尔坐标的坐标点上；对第二K空间数据进行校正处理以生成第三K空间数据，以使第三K空间数据对应的FOV小于第二K空间数据对应的FOV；基于并行成像方法和/或压缩感知方法确定第三K空间数据对应的磁共振图像。K空间数据量的减少使采用并行成像方法和/或压缩感知方法进行图像重建成为了可能，从而可以显著提高图像重建的速度，满足临床实时成像的需求。
实施例二
图4是本发明实施例二提供的磁共振成像方法的流程图。本发明实施例在上述实施例的基础上，增加了确定第三K空间数据的权重参数的步骤。
相应地，本实施例的方法包括：
S201、获取受检者的第一K空间数据，第一K空间数据基于非笛卡尔采样方式获取。
S202、对第一K空间数据进行重栅格化以生成第二K空间数据。
S203、对第二K空间数据进行校正处理以生成第三K空间数据，以使第三K空间数据对应的FOV小于第二K空间数据对应的FOV。
S204、确定与第三K空间数据对应的权重参数。
为了提高基于第三K空间数据重建的磁共振图像的准确性，本实施例在使用第三K空间数据进行图像重建之前，还通过权重参数对其进行校正。如图5所示，权重参数的确定方法可选为：
S2041、将第一K空间数据中的每个数据重置为预设数值以得到初始权重参数。
其中，预设数值通常根据实际需求进行设置，本实施例优选将其设置为1。
S2042、对初始权重参数进行重栅格化以生成中间权重参数。
为了使权重参数与第三K空间数据相对应，本实施例对初始权重参数也进行重栅格化以生成中间权重参数。
S2043、对中间权重参数对应的FOV进行处理以生成矩形FOV，并将矩形FOV对应的各个权重参数的模值作为权重参数。
在得到中间权重参数之后，对该中间数据对应的FOV进行处理以生成矩形FOV，且该矩形FOV与第三K空间数据对应的矩形FOV一致，然后将该矩形FOV对应的各个权重参数的模值作为权重参数，由于权重参数对应的FOV与第三K空间数据对应的FOV相同，因此权重参数与第三K空间数据一一对应。
S205、基于并行成像方法和/或压缩感知方法确定第三K空间数据和权重参数对应的磁共振图像。
相较于现有技术，本发明实施例在前述实施例的基础上，通过改变MRI图像的FOV减少了参与图像重建的数据量，使基于非笛卡尔采样方式采集的磁共振数据可以通过并行成像和/或压缩感知方法进行图像重建，通过权重参数来校正参与图像重建的第三K空间数据，可以提高重建出的磁共振图像的图像质量。
实施例三
图6是本发明实施例三提供的磁共振成像方法的流程图。本发明实施例在上述任意实施例的基础上，增加了密度校正函数的确定步骤，包括：
S301、获取受检者的第一K空间数据，第一K空间数据基于非笛卡尔采样方式获取。
S302、对第一K空间数据进行重栅格化以生成第二K空间数据。
S303、对第二K空间数据进行校正处理以生成第三K空间数据，以使第三K空间数据对应的FOV小于第二K空间数据对应的FOV。
S304、确定第三K空间数据的权重参数。
S305、确定第三K空间数据的密度校正函数。
为了提高磁共振图像的质量，通常需要通过密度校正函数对磁共振数据进行密度校正，基于此，本实施例基于现有技术的密度校正函数方法确定第一K空间数据的初始密度校正函数。该初始密度校正函数与第一K空间数据一一对应，而第三K空间数据与第一K空间数据并不对应，为了得到第三K空间的密度校正函数，本实施例以预设插值方法对初始密度校正函数进行插值，以得到与第三K空间数据对应的密度校正函数。其中，预设插值方法采用现有技术的插值方法即可，本实施例在此不对其进行限定。
另外，第一K空间数据通常包括中心满采区和外周的非满采区，本实施例也可以将第一K空间数据的中心满采区除以权重参数以更新第一K空间数据。基于更新后的第一K空间数据得到的第三K空间数据的中心满采区的权重参数和密度校正函数均为1。
需要说明的是，本实施例不对S304和S305的执行顺序进行限定，即二者的执行顺序可以调换。
S306、基于并行成像方法和/或压缩感知方法确定第三K空间数据、权重参数和密度校正函数对应的磁共振图像。
第三K空间数据、权重参数和密度校正函数得到之后，即可基于并行成像方法和/或压缩感知方法确定第三K空间数据、权重参数和密度校正函数对应的磁共振图像，比如：根据第三K空间数据、权重参数和密度校正函数确定线圈敏感度数据；然后根据第三K空间数据、权重参数、密度校正函数和线圈敏感度数据确定磁共振图像。
其中，根据第三K空间数据、权重参数、密度校正函数和线圈敏感度数据确定磁共振图像，可通过求解以下目标函数实现：
Figure BDA0002007037520000101
其中，w1为权重参数；w2为密度校正函数；F为傅里叶变换，C为线圈敏感度数据，d为重建后的图像或者多幅动态图像；Sc为第三K空间数据；T为稀疏变换。
在前述实施例的基础上引入密度校正函数，可以显著提高基于并行成像和/或压缩感知方法确定的第三K空间数据对应的磁共振图像的图像质量。
实施例四
图7是本发明实施例四提供的磁共振成像装置的结构框图。该装置用于执行上述任意实施例所提供的磁共振成像方法，该装置可选为软件或硬件实现。
该装置包括：
获取模块1201，用于获取受检者的第一K空间数据，第一K空间数据基于非笛卡尔采样方式获取；
第二K空间数据生成模块1202，用于对第一K空间数据进行重栅格化以生成第二K空间数据；
第三K空间数据生成模块1203，用于对第二K空间数据进行校正处理以生成第三K空间数据，以使第三K空间数据对应的FOV小于第一K空间数据对应的FOV；
图像重建模块1204，用于基于并行成像方法和/或压缩感知方法确定第三K空间数据对应的磁共振图像。
本发明实施例提供的磁共振成像装置的技术方案，通过获取模块获取受检者的第一K空间数据，第一K空间数据基于非笛卡尔采样方式获取；通过第二K空间数据生成模块对第一K空间数据进行重栅格化以生成第二K空间数据，以将磁共振数据转换至笛卡尔坐标的整数格点上；通过第三K空间数据生成模块对第二K空间数据进行校正处理以生成第三K空间数据，以使第三K空间数据对应的FOV小于第二K空间数据对应的FOV；通过图像重建模块基于并行成像方法和/或压缩感知方法确定第三K空间数据对应的磁共振图像。K空间数据量的减少使采用并行成像方法和/或压缩感知方法进行图像重建成为了可能，从而可以显著提高图像重建的速度，以满足临床实时成像的需求。
本发明实施例所提供的磁共振成像装置可执行本发明任意实施例所提供的磁共振成像方法，具备执行方法相应的功能模块和有益效果。
实施例五
图8是本发明实施例五提供的磁共振成像***的结构框图。该磁共振成像***100包括MR扫描器110和处理器120，MR扫描器110用于获取受检者的第一K空间数据，第一K空间数据基于非笛卡尔采样方式获取；处理器120用于获取第一K空间数据，以及对第一K空间数据进行重栅格化以生成第二K空间数据；以及对第二K空间数据进行校正处理以生成第三K空间数据，以使第三K空间数据对应的FOV小于第二K空间数据对应的FOV；以及基于并行成像方法和/或压缩感知方法确定第三K空间数据对应的磁共振图像。
其中，对第二K空间数据进行数据校正以生成第三K空间数据，包括：对第二K空间数据对应的FOV进行处理以生成矩形FOV，并将该矩形FOV对应的K空间数据作为第三K空间数据。
相对于现有技术来说，本实施例处理器可以对非笛卡尔坐标系下的K空间数据进行重栅格化和校正处理，使经过重栅格化和校正处理后的磁共振数据即第三K空间数据，能够通过并行成像方法和/或压缩感知方法进行图像重建。从而可以显著提高图像重建的速度，以满足临床实时成像的需求。
需要说明的是，本发明实施例的处理器被配置为实现前述实施例所述的磁共振成像方法，具体的磁共振成像方法可参见前述实施例，本实施例在此不予赘述。
如图9所示，该磁共振成像***100还包括控制器130和输出装置140，其中，处理器130可同时监测或控制MR扫描器110、处理器120和输出装置140。控制器130可以包括中央处理器(Central Processing Unit，CPU)、专门应用集成电路(Application-SpecificIntegrated Circuit，ASIC)、专用指令处理器(Application Specific Instruction SetProcessor，ASIP)、图形处理单元(Graphics Processing Unit，GPU)、物理处理器(PhysicsProcessing Unit，PPU)、数字信号处理器(Digital Processing Processor，DSP)、现场可编程逻辑门阵列(Field-Programmable Gate Array，FPGA)、ARM处理器等中的一种或几种的组合。
输出装置140，比如显示器，可显示感兴趣区域的磁共振图像。进一步地，输出装置140还可显示受检者的身高、体重、年龄、成像部位、以及MR扫描器110的工作状态等。输出装置140的类型可以是阴极射线管(CRT)输出装置、液晶输出装置(LCD)、有机发光输出装置(OLED)、等离子输出装置等中的一种或几种的组合。
磁共振成像***100可连接一个局域网(Local Area Network，LAN)、广域网(WideArea Network，WAN)、公用网络、私人网络、专有网络、公共交换电话网(Public SwitchedTelephone Network，PSTN)、互联网、无线网络、虚拟网络、或者上述网络的任何组合。
MR扫描器110包括MR信号获取模块、MR控制模块和MR数据存储模块。其中，MR信号获取模块包括磁体单元和射频单元。磁体单元主要包括产生B0主磁场的主磁体和产生梯度的梯度组件。磁体单元中包含的主磁体可以是永磁体或超导磁体，梯度组件主要包含梯度电流放大器(AMP)、梯度线圈，梯度组件还可包含三个独立通道Gx、Gy、Gz，每个梯度放大器激发梯度线圈组中对应的一个梯度线圈，产生用于生成相应空间编码信号的梯度场，以对磁共振信号进行空间定位。射频单元主要包括射频发射线圈和射频接收线圈，射频发射线圈用于向受检者或人体发射射频脉冲信号，射频接收线圈用于接收从人体采集的磁共振信号，且根据功能的不同，组成射频单元的射频线圈可分为体线圈和局部线圈。在一个实施例中，体线圈或局部线圈的种类可以是鸟笼形线圈、螺线管形线圈、马鞍形线圈、亥姆霍兹线圈、阵列线圈、回路线圈等。在一个具体实施例中，局部线圈设置为阵列线圈，且该阵列线圈可设置为4通道模式、8通道模式或16通道模式。磁体单元和射频单元可组成开放性低场磁共振装置或者封闭型超导磁共振装置。
MR控制模块可监测包含磁体单元和射频单元的MR信号获取模块、MR数据处理模块。具体地，MR控制模块可接收MR信号获取模块发送的信息或者脉冲参数；此外，MR控制模块还可控制MR数据处理模块的处理过程。在一个实施例中，MR控制模块还连接有包含脉冲序列发生器、梯度波形发生器、发射机和接收机等，在接受用户从控制台发出的指令后，控制磁场模块执行相应扫描序列。
示例性地，本发明MR扫描器110产生MR数据的具体过程包括：主磁体产生B0主磁场，受检者体内的原子核在主磁场作用下产生进动频率，该进动频率与主磁场强度呈正比；MR控制模块存储和发送需要执行的扫描序列(scan sequence)的指令，脉冲序列发生器根据扫描序列指令对梯度波形发生器和发射机进行控制，梯度波形发生器输出具有预定时序和波形的梯度脉冲信号，该信号经过Gx、Gy和Gz梯度电流放大器，再通过梯度组件中的三个独立通道Gx、Gy、Gz，每个梯度放大器激发梯度线圈组中对应的一个梯度线圈，产生用于生成相应空间编码信号的梯度场，以对磁共振信号进行空间定位；脉冲序列发生器还执行扫描序列，输出包括射频发射的射频脉冲的计时、强度、形状等数据以及射频接收的计时和数据采集窗口的长度到发射机，同时发射机将相应射频脉冲发送至射频单元中的体发射线圈产生B1场，在B1场作用下病人体内被激发的原子核发出的信号被射频单元中的接收线圈感知到，然后通过发送/接收开关传输到MR数据处理模块，经过放大、解调、过滤、AD转换等数字化处理，然后传输到MR数据存储模块。当MR数据存储模块获取一组原始的k-空间数据后，扫描结束。原始的k-空间数据被重新整理成与每个将被重建的图像对应的单独的k-空间数据组，每个k-空间数据组被输入到阵列处理器，进行图像重建后结合磁共振信号，形成一组图像数据。
实施例六
本发明实施例六还提供一种包含计算机可执行指令的存储介质，所述计算机可执行指令在由计算机处理器执行时用于执行磁共振成像方法，该方法包括：
采用获取受检者的第一K空间数据，所述第一K空间数据基于非笛卡尔采样方式获取；
对所述第一K空间数据进行重栅格化以生成第二K空间数据；
对所述第二K空间数据进行校正处理以生成第三K空间数据，以使第三K空间数据对应的FOV小于第二K空间数据对应的FOV；
基于并行成像方法和/或压缩感知方法确定所述第三K空间数据对应的磁共振图像。
当然，本发明实施例所提供的一种包含计算机可执行指令的存储介质，其计算机可执行指令不限于如上所述的方法操作，还可以执行本发明任意实施例所提供的磁共振成像方法中的相关操作。
通过以上关于实施方式的描述，所属领域的技术人员可以清楚地了解到，本发明可借助软件及必需的通用硬件来实现，当然也可以通过硬件实现，但很多情况下前者是更佳的实施方式。基于这样的理解，本发明的技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来，该计算机软件产品可以存储在计算机可读存储介质中，如计算机的软盘、只读存储器(Read-Only Memory，简称ROM)、随机存取存储器(RandomAccess Memory，简称RAM)、闪存(FLASH)、硬盘或光盘等，包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等)执行本发明各个实施例所述的磁共振成像方法。
值得注意的是，上述磁共振成像装置的实施例中，所包括的各个单元和模块只是按照功能逻辑进行划分的，但并不局限于上述的划分，只要能够实现相应的功能即可；另外，各功能单元的具体名称也只是为了便于相互区分，并不用于限制本发明的保护范围。
注意，上述仅为本发明的较佳实施例及所运用技术原理。本领域技术人员会理解，本发明不限于这里所述的特定实施例，对本领域技术人员来说能够进行各种明显的变化、重新调整和替代而不会脱离本发明的保护范围。因此，虽然通过以上实施例对本发明进行了较为详细的说明，但是本发明不仅仅限于以上实施例，在不脱离本发明构思的情况下，还可以包括更多其他等效实施例，而本发明的范围由所附的权利要求范围决定。

Claims (7)

1.一种磁共振成像方法，其特征在于，包括：
获取受检者的第一K空间数据，所述第一K空间数据基于非笛卡尔采样方式获取；
对所述第一K空间数据进行重栅格化以生成第二K空间数据；
对所述第二K空间数据进行校正处理以生成第三K空间数据，以使所述第三K空间数据对应的FOV小于所述第二K空间数据对应的FOV；
确定所述第三K空间数据的权重参数；
在得到所述第三K空间数据之后，还包括：
确定第一K空间数据的初始密度校正函数；
基于预设插值法，根据第一K空间数据对应的初始密度校正函数和所述第三K空间数据的各个坐标点确定所述第三K空间数据对应的密度校正函数；
基于并行成像方法和/或压缩感知方法，确定所述第三K空间数据、权重参数和密度校正函数对应的磁共振图像，包括：
根据所述第三K空间数据、权重参数和密度校正函数确定线圈敏感度数据；
根据所述第三K空间数据、权重参数、密度校正函数和所述线圈敏感度数据确定磁共振图像；
所述根据第三K空间数据、权重参数、密度校正函数和线圈敏感度数据确定磁共振图像，包括：
根据以下目标函数确定磁共振图像：
Figure FDA0003632624260000011
其中，w1为权重参数；w2为密度校正函数；F为傅里叶变换，C为线圈敏感度数据，d为重建后的图像或者多幅动态图像；Sc为第三K空间数据；T为稀疏变换。
2.根据权利要求1所述的方法，其特征在于，所述对所述第二K空间数据进行校正处理以生成第三K空间数据，包括：
对所述第二K空间数据对应的FOV进行处理以生成矩形FOV，并将所述矩形FOV对应的K空间数据作为第三K空间数据。
3.根据权利要求1所述的方法，其特征在于，所述确定第三K空间数据对应的权重参数，包括：
将第一K空间数据中的每个数据重置为预设数值以得到初始权重参数；
对所述初始权重参数进行重栅格化以生成中间权重参数；
对所述中间权重参数对应的FOV进行处理以生成矩形FOV，并将所述矩形FOV对应的各个权重参数的模值作为权重参数。
4.根据权利要求1所述的方法，其特征在于，所述第一K空间数据包括中心满采区和外周非满采区，相应的，在所述确定所述第三K空间数据的权重参数之后，还包括：
将所述第一K空间数据的中心满采区的数据除以权重参数，以更新所述第一K空间数据。
5.一种磁共振成像***，其特征在于，包括：
MR扫描器，用于获取受检者的第一K空间数据，所述第一K空间数据基于非笛卡尔采样方式获取；
处理器，用于获取所述第一K空间数据，以及对所述第一K空间数据进行重栅格化以生成第二K空间数据；以及对所述第二K空间数据进行校正处理以生成第三K空间数据，以使所述第三K空间数据对应的FOV小于所述第二K空间数据对应的FOV；
所述处理器，还用于确定所述第三K空间数据的权重参数；
所述处理器，还用于确定第一K空间数据的初始密度校正函数；
基于预设插值法，根据第一K空间数据对应的初始密度校正函数和所述第三K空间数据的各个坐标点确定所述第三K空间数据对应的密度校正函数；
所述处理器，还用于根据所述第三K空间数据、权重参数和密度校正函数确定线圈敏感度数据；
根据所述第三K空间数据、权重参数、密度校正函数和所述线圈敏感度数据确定磁共振图像；
所述处理器，还用于根据以下目标函数确定磁共振图像：
Figure FDA0003632624260000031
其中，w1为权重参数；w2为密度校正函数；F为傅里叶变换，C为线圈敏感度数据，d为重建后的图像或者多幅动态图像；Sc为第三K空间数据；T为稀疏变换。
6.根据权利要求5所述的***，其特征在于，所述对所述第二K空间数据进行校正处理以生成第三K空间数据，包括：
对所述第二K空间数据对应的FOV进行处理以生成矩形FOV，并将所述矩形FOV对应的K空间数据作为第三K空间数据。
7.一种包含计算机可执行指令的存储介质，其特征在于，所述计算机可执行指令在由计算机处理器执行时用于执行如权利要求1-4中任一所述的磁共振成像方法。
CN201910232085.9A 2019-03-26 2019-03-26 磁共振成像方法、***及存储介质 Active CN109975343B (zh)

Priority Applications (2)

Application Number Priority Date Filing Date Title
CN201910232085.9A CN109975343B (zh) 2019-03-26 2019-03-26 磁共振成像方法、***及存储介质
US16/813,881 US11592506B2 (en) 2019-03-26 2020-03-10 System and method for magnetic resonance imaging

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910232085.9A CN109975343B (zh) 2019-03-26 2019-03-26 磁共振成像方法、***及存储介质

Publications (2)

Publication Number Publication Date
CN109975343A CN109975343A (zh) 2019-07-05
CN109975343B true CN109975343B (zh) 2022-06-28

Family

ID=67080649

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910232085.9A Active CN109975343B (zh) 2019-03-26 2019-03-26 磁共振成像方法、***及存储介质

Country Status (2)

Country Link
US (1) US11592506B2 (zh)
CN (1) CN109975343B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US11348291B2 (en) 2019-11-29 2022-05-31 Shanghai United Imaging Intelligence Co., Ltd. System and method for reconstructing magnetic resonance images
CN111257809B (zh) * 2020-01-23 2022-05-17 上海东软医疗科技有限公司 磁共振成像方法、装置、存储介质及医疗设备

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104965184A (zh) * 2015-05-25 2015-10-07 厦门大学 基于时空编码单扫描磁共振成像的螺旋采样及重建方法
WO2017007279A1 (ko) * 2015-07-09 2017-01-12 성균관대학교산학협력단 자기 공명 영상 장치 및 방법

Family Cites Families (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8040328B2 (en) * 2000-10-11 2011-10-18 Peter Smith Books, papers, and downloaded information to facilitate human interaction with computers
DE102004024459B4 (de) * 2004-05-14 2007-09-20 Siemens Ag Verfahren zur Magnetresonanz-Bildgebung mittels einer partiellen parallelen Akquisitionstechnik bei nicht-kartesischer Belegung des k-Raums
DE102005046732B4 (de) * 2005-04-18 2010-04-15 Siemens Ag Verbessertes Rekonstruktionsverfahren bei der Propellerbildgebung in der Magnetresonanztomographie
EP2521987B1 (en) * 2010-01-04 2024-05-08 Mayo Foundation For Medical Education And Research System and method for combined time-resolved magnetic resonance angiography and perfusion imaging
US20130289912A1 (en) * 2012-03-30 2013-10-31 Siemens Aktiengesellschaft Eigen-vector approach for coil sensitivity maps estimation
WO2013159044A1 (en) 2012-04-19 2013-10-24 New York University System, method and computer-accessible medium for highly-accelerated dynamic magnetic resonance imaging using golden-angle radial samplng and compressed sensing
WO2015019970A1 (ja) * 2013-08-08 2015-02-12 株式会社東芝 画像処理装置及び磁気共鳴イメージング装置
KR101857795B1 (ko) * 2013-12-18 2018-05-14 삼성전자주식회사 자기 공명 영상 장치 및 그 동작방법
US9983283B2 (en) * 2015-03-16 2018-05-29 Toshiba Medical Systems Corporation Accelerated MRI using radial strips and undersampling of k-space
JP6811010B2 (ja) * 2015-10-26 2021-01-13 ミヨシ油脂株式会社 樹脂に添加されるフィラー、樹脂材料とその弾性率を調整する方法および製造方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104965184A (zh) * 2015-05-25 2015-10-07 厦门大学 基于时空编码单扫描磁共振成像的螺旋采样及重建方法
WO2017007279A1 (ko) * 2015-07-09 2017-01-12 성균관대학교산학협력단 자기 공명 영상 장치 및 방법

Also Published As

Publication number Publication date
US20200309883A1 (en) 2020-10-01
US11592506B2 (en) 2023-02-28
CN109975343A (zh) 2019-07-05

Similar Documents

Publication Publication Date Title
US10670678B2 (en) MR imaging using stack-of stars acquisition
JP5048052B2 (ja) 磁気共鳴イメージング装置および磁気共鳴イメージング方法
JP6513398B2 (ja) 事前情報に制約される正則化を使用するｍｒ画像再構成
US6448771B1 (en) Magnetic resonance method for forming a fast dynamic image
CN109188326B (zh) 磁共振成像方法和磁共振***
WO2020019412A1 (zh) 磁共振成像方法、装置、设备及存储介质
US7372269B2 (en) Magnetic resonance imaging method and apparatus
US9535146B2 (en) Magnetic resonance imaging apparatus and magnetic resonance imaging method
CN109907759B (zh) 磁共振成像方法和***
CN107209238B (zh) 具有对边带伪迹的抑制的并行多切片mr成像
WO2012077543A1 (ja) 磁気共鳴イメージング装置及びコントラスト強調画像取得方法
CN110664378B (zh) 磁共振成像方法、装置、***及存储介质
US10317487B2 (en) Method and magnetic resonance apparatus for acquisition of MR data of a two-dimensional volume segment
JP2005144173A (ja) Ｍｒｔ画像形成方法
CN109975343B (zh) 磁共振成像方法、***及存储介质
JP2017529960A (ja) アーチファクト抑制を有するプロペラｍｒイメージング
US9274197B2 (en) Magnetic resonance imaging data sampling methods and systems
CN112272535A (zh) 磁共振弥散成像的方法和***
US20160124065A1 (en) Method and apparatus for correction of magnetic resonance image recordings with the use of a converted field map
JP2000287949A (ja) 磁気共鳴イメージング装置
JPH11262479A (ja) 連続ｍｒｉ画像再構成方法及び装置
US10254367B2 (en) Magnetic resonance imaging method and apparatus with motion-corrected model-based acceleration of parameter mapping
US11269037B2 (en) MR imaging using motion-dependent radial or spiral k-space sampling
CN114325524B (zh) 磁共振图像重建方法、装置、***及存储介质
CN108577841B (zh) 一种propeller技术中抑制非刚性运动的权重计算方法

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