CN110623666B - 磁共振指纹方法 - Google Patents

磁共振指纹方法 Download PDF

Info

Publication number
CN110623666B
CN110623666B CN201910540740.7A CN201910540740A CN110623666B CN 110623666 B CN110623666 B CN 110623666B CN 201910540740 A CN201910540740 A CN 201910540740A CN 110623666 B CN110623666 B CN 110623666B
Authority
CN
China
Prior art keywords
parameter
determined
comparison
value
signal
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
CN201910540740.7A
Other languages
English (en)
Other versions
CN110623666A (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.)
Siemens Healthineers AG
Original Assignee
Siemens Healthineers AG
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 Siemens Healthineers AG filed Critical Siemens Healthineers AG
Publication of CN110623666A publication Critical patent/CN110623666A/zh
Application granted granted Critical
Publication of CN110623666B publication Critical patent/CN110623666B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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/4828—Resolving the MR signals of different chemical species, e.g. water-fat imaging
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00—Arrangements or instruments for measuring magnetic variables
    • G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48—NMR imaging systems
    • G01R33/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
    • G01R33/56—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
    • G01R33/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
    • 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
    • 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/28—Details of apparatus provided for in groups G01R33/44 - G01R33/64
    • G01R33/38—Systems for generation, homogenisation or stabilisation of the main or gradient magnetic field
    • G01R33/385—Systems for generation, homogenisation or stabilisation of the main or gradient magnetic field using gradient magnetic field coils
    • 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/443—Assessment of an electric or a magnetic field, e.g. spatial mapping, determination of a B0 drift or dosimetry
    • 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/50—NMR imaging systems based on the determination of relaxation times, e.g. T1 measurement by IR sequences; T2 measurement by multiple-echo 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
    • 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

Landscapes

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

Abstract

一种磁共振指纹(MRF)方法，包括：建立的体素时间序列的信号特性与第一比较信号特性(字典)的第一信号比较以及基于第一比较信号特性和第一信号比较中确定的值生成另外(合成的)比较信号特性。使用所生成的另外比较信号特性来执行通过其可以确定至少第一另外参数和第二另外参数的值的另一信号比较。通过生成另外比较信号特性，可以确定在第一信号比较中不一定已经确定的至少一个另外参数的值。在此期间，所需的处理能力保持较低。通过在另一信号比较的过程中确定第一另外参数和第二另外参数，该方法使得能够整体地确定另外参数的值，因此可以避免由其他另外参数引起的效应造成的另外参数的值的伪造。

Description

磁共振指纹方法
技术领域
本发明涉及一种用于改进对检查对象的局部参数值的确定的磁共振指纹方法。
背景技术
磁共振技术(下面的缩写MR代表磁共振)是一种已知技术，通过该技术可以生成检查对象的内部的图像。简而言之，出于该目的，检查对象位于磁共振装置中，在相对强的、静态的、均匀的主磁场(也被称为B0场)中，其中场强为0.2特斯拉至7特斯拉及以上，使得在主磁场的方向上获得其核自旋。为了触发核自旋共振，检查对象被高频激励脉冲(RF脉冲)照射，所触发的核自旋共振被测量为所谓的k空间数据，并且基于此，该MR图像被重建或者光谱数据被识别。为了对测量数据进行空间编码，将快速切换磁梯度场覆加在主磁场上，并且这些磁场铺设轨道，沿着轨道k空间中的测量数据被读出。捕获的测量数据被数字化并作为复数值存储在k空间矩阵中。可以例如通过多维傅里叶变换根据填充有值的k空间矩阵重建关联的MR图像。要发射的一系列RF脉冲、要激活的梯度以及用于该目的并以特定方式排序的读出过程被指定为序列。
已知对描述检查的检查对象中的材料的参数(例如纵向弛豫时间T1、横向弛豫时间T2和质子密度)具有不同的灵敏度水平的各种序列类型。根据使用特定序列类型捕获的测量数据重建的MR图像以根据所使用的序列类型的灵敏度加权的方式表示检查对象的图。
使用磁共振系统产生的磁共振成像可用于确定检查对象中材料的存在和/或分布。该材料可以是例如检查对象的(可能是病理性的)组织、造影剂、标记物质或代谢产物。
这里，可以以各种方式从捕获的测量数据获得关于存在的材料的信息。相对简单的信息源是例如根据测量数据重建的图像数据。然而，还存在更复杂的方法，其例如从根据连续测量的测量数据的集重建的图像数据的体素时间序列中识别关于所检查的检查对象的信息。
借助于定量MR成像技术，可以确定被测对象的绝对性质，例如在患者的情况下的组织特异性T1和T2弛豫时间。与此相反，通常在临床日常中使用的常规序列仅生成不同组织类型的相对信号强度(所谓的权重)，具有以下结果：诊断解释高度依赖于放射科医师进行的主观评估。因此，定量技术提供了客观可比较的明显优势，但由于其长测量时间，在当前的日常中很难使用。
相对较新的定量测量方法，例如磁共振指纹方法(MRF方法)，可以将所述长测量时间的缺点降低到可接受的水平。在MRF方法中，对于每个体素(voxel)或至少在图像数据中感兴趣的体素，从利用不同成像参数连续捕获的测量数据重建的图像数据的信号特性被认为是体素时间序列，其中体素时间序列的信号特性被认为是在测量期间在检查对象的体素位置中占优势的参数的“指纹”。使用模式识别方法将这些信号特性与是特定材料的特性的信号特性的先前识别的数据库(所谓的“字典”)中的信号特性进行比较，以便在绘制的检查对象中识别根据测量数据重建的图像数据中表示的材料，或者特定于组织的参数的空间分布(例如横向弛豫时间T2或纵向弛豫时间T1-所谓的T1和T2图)。在这种情况下，这种字典中的信号特性也可以通过模拟生成。
因此，该方法的原理是将测量的信号特性与多个先前已知的信号特性进行比较。这里，可以针对字典识别T1和T2弛豫时间以及其他参数的不同组合的信号特性。对于要确定的每个参数，存在所谓的字典中相应的“维度”，其中包括相应参数的不同参数值，以便使不同的比较值可用。然后，特别地通过将测量的信号特性与模拟信号特性中的全部或一些进行比较来确定图像中的体素(可以是像素或体素)的参数值，例如T1和T2。该过程被称为“匹配”。字典中最接近测量信号特性的信号特性确定相应体素的参数，例如弛豫参数T1和T2。
除了横向弛豫时间T2(也称为自旋-自旋弛豫)之外，还已知比横向弛豫时间T2更短的有效横向弛豫时间T2*。有效横向弛豫时间T2*另外考虑主磁场B0中的局部变化并因此考虑偏共振效应，这导致横向磁化的更快衰减。两个横向弛豫参数T2和T2*之间的关系可被表示为如下：
其中T2'是由主磁场B0中的局部变化引起的横向磁化的衰减，以及借此的通过相分散横向磁化的衰减。
这里在原则上，除了上面提到的被检查对象的特定于组织的参数之外，还可以识别特定于测量的参数，例如所施加的磁场的场强，或者真正被照射的高频场B1+的强度的局部分布，这是由于由MR技术捕获的信号可能取决于被研究对象中的特定于组织的参数以及描述测量期间占优势的条件的特定于测量的参数。在这种情况下使用的成像参数被选择成使得捕获的测量数据显示对要确定的期望参数的依赖性。例如，对期望参数敏感的序列类型可以用于MRF方法。由于捕获参数中的依赖性和变化以及比较信号特性中对其的考虑，可以根据已经以这种方式捕获的体素时间序列来确定期望的参数。
对于MRF方法，原则上任何回波技术可以与用于k空间扫描的任何方法(例如笛卡尔、螺旋、径向)结合使用。
考虑所使用的字典中的特定于组织的参数T1和T2并且在测量的体素时间序列中确定它们的MRF方法在例如Ma等人的文章“Magnetic Resonance Fingerprinting”，Nature，495：第187-192页(2013)中描述。在那里，基于TrueFISP(具有稳态自由进动的真实快速成像)的序列与螺旋k空间扫描结合使用。US 20160061922A1描述了另一MRF方法。
原则上，在字典中的比较信号特性中没有考虑的任何效应——例如因为它们不在根据其生成字典的信号模型中，但是对在生成体素时间序列时捕获的信号有影响——可能伪造MRF实验的结果，也就是说确定的参数值。根据生成体素时间序列时使用的序列，这种效应可能或多或少显着。
在MR实验(MR测量)中可能对体素时间序列中测量的信号具有影响的这种效应是例如由于偏共振基于因此由于体素内相分散造成的信号衰减。如上面描述的，该效应的影响通常通过另一弛豫时间——有效横向弛豫时间T2*来量化。为了识别对象内的每个体素的相分散的差异，可以使用其中信号响应取决于T2*的序列以例如生成体素时间序列。由于T2*被假定为指数衰减(例如，利用体素内的偏共振的高斯分布的另一假设)，因此可以例如使用依赖于T2*的多回波FLASH序列来扫描该信号衰减。使用以这种方式捕获的体素时间序列，以及取决于T2*的对应的比较信号特性，T2*时间的参数可以通过MRF方法中的常规的信号比较来确定。
相反，如果没有考虑参数的可能效应，例如所述相分散(也称为体素内相移)，例如在引用的Ma等人的文章中的情况，则这也可能伪造确定其他参数的值的结果。例如，在Chiu等人的文章“Effects of RF pulse profile and intra-voxel phase dispersion on MRfingerprinting with balanced SSFP readout”，Magnetic Resonance Imaging 41，第80-86页，2017中，基于引用的Ma等人的文章中的MRF实现描述了MRF结果对体素内相分散的依赖性。
为了减少这种依赖性，等人的文章“Pseudo Steady-State FreePrecession for MR Fingerprinting””，Magnetic Resonance in Medicine 77，第1151-1161页，2017中提出了一种旨在通过改变体素时间序列的生成中使用的序列来减少对体素内相移的所述灵敏度的MRF实现。
为了生成体素时间序列，其他MRF实现使用这样序列：其中尽管发生体素内相移的这种效应，但它们相对较小。示例是Jiang等人在在文章“MR Fingerprinting Using FastImaging with Steady State Precession(FISP)with Spiral Readout”，MagneticResonance in Medicine 74，第1621-1631页，2015中描述的使用FISP(具有稳态进动的快速成像)序列的MRF实现，或使用FLASH(快速小角度激发)序列的MRF实现。利用这些序列类型，在每个重复时间TR中，横向磁化被完全消灭(“破坏”)。因此，相分散的效果被限制于更快的信号衰减。以这种方式，例如通过在序列的过程中保持回波时间TE恒定，也可以保持相分散的影响恒定，因此MRF结果最多经受低水平伪造。
在Wang等人的文章“In vivo Simultaneous Measurement ofδf,T1,T2and T2*byMagnetic Resonance Fingerprinting with Quadratic RF Phase”，Proc.Intl.Soc.Mag.Reson.Med.25，第3960页，2017中描述了另一种MRF方法，其中在一种匹配类型中为每个体素分配相分散。那里，在确定T1、T2和B0(δf)的参数值之后，在另一匹配步骤中确定相分散。对于另一匹配步骤，使用不同宽度的以大约对B0确定的值为中心的洛伦兹分布对来自字典的B0维度中的信号进行加权，并求和。使用另一匹配步骤和由此获得的信号来识别“最佳拟合”宽度，并将其用作用于相分散的测量。在另一匹配步骤中，保留已确定的T1、T2和B0的参数值，其中仅确定相分散(在整个宽度上)。
发明内容
本发明的目的是使得使用MRF方法能够改进对参数值的确定。
该目的通过根据权利要求1所述的一种用于使用磁共振指纹(MRF)技术来确定检查对象的体素中的参数值的方法，根据权利要求13所述的磁共振系统，根据权利要求14所述的计算机程序以及根据权利要求15所述的电子可读数据介质来实现。
根据本发明的一种用于使用磁共振指纹(MRF)技术来确定检查对象的体素中的参数值的方法包括以下步骤：
借助于MRF成像方法建立至少一个体素时间序列(BZS)，其中，要根据所述至少一个体素时间序列来确定由相应体素表示的检查对象中的位置处的至少一个参数的相应值(P，P')，
执行所建立的体素时间序列(BZS)的相应信号特性的至少一个片段与第一比较信号特性(D)的对应片段的第一信号比较(103)，以确定要被确定的参数中的第一参数的相应的值(P)，
基于第一比较信号特性和第一信号比较中确定的值来生成另外比较信号特性(D')，
执行所建立的体素时间序列(BZS)的相应信号特性的至少一个片段与所生成的另外比较信号特性(D')的对应片段的另一信号比较，以分别(独立地)确定要被确定的参数中的至少第一另外参数和第二另外参数的值，
输出针对相应体素确定的要被确定的参数的值。
通过使用在第一信号比较中确定的参数值创造性地根据所提供的比较信号特性的集来生成另外比较信号特性，可以确定不一定在第一信号比较中已经确定的至少一个另外参数的值。因此，特别地，在第一信号比较中可能无法确定值在另一信号比较中确定的第二另外参数。这不仅增加了可确定参数的数量，而且同时使另一信号比较中的处理能力保持很小，这是由于另外比较信号特性不需要“实际”的另一维度以确定其他参数，但是，可以说，合成地将在第一信号比较中确定的参数的维度转换成另一其他参数的维度。从一类子字典生成的另外(模拟)比较信号特性，其避免了“实际”额外的参数维度(对于另一参数)并且因此与使用这样的“实际”另一参数维度相比显着地降低了所需的处理能力。
通过在另一信号比较的过程中确定第一另外参数和第二另外参数，该方法使至少这两个另外参数的值能够被整体地确定，因此可以避免根据其他参数的值引起的效应造成的一个参数的值的伪造。如果第二另外参数对第一个另外参数具有影响或依赖性，则这首先特别有用，并且视情况而定反之亦然。
在第一信号比较期间，可以确定局部主磁场的值，并且作为第一另外参数和第二另外参数，例如横向弛豫时间T2的依赖于相应的体素中占优势的相分散的参数，例如相分散本身可以被选择，因此可以避免由相分散的效应造成的所确定的横向弛豫时间T2的值的伪造。因为如此确定的体素中的相分散的值可以与体素中表示的组织的易感性相关联，所以相分散值的确定也可以用于提供诊断问题的答案。
在第一信号比较期间，还可以确定纵向弛豫时间T1的值或横向弛豫时间T2的值，并且第一另外参数可以是纵向弛豫时间T1或横向弛豫时间T2，并且第二另外参数可以是纵向弛豫与横向弛豫之间的比率。通过使用另一信号比较来确定纵向弛豫与横向弛豫的占优势的比率，可以特别地避免所谓的“部分体积”效应。
以这种方式，所描述的方法使得能够尽可能自由地使用对不同参数具有不同依赖性的不同序列类型。这使得可以使用MRF方法尽可能精确地确定尽可能多的不同参数的值，而不必接受由于因为使用序列类型对在生成体素时间序列时捕获的信号具有影响的其他参数的效应造成的某些参数值的伪造，并且不会过度增加处理能力。
这里，在另一信号比较期间可以“保持”在第一信号比较期间确定的显示对第二另外参数无依赖性的参数的值，也就是说，它们不再被确定。因为不存在所述依赖性意味着不需要预期对在第一信号比较中确定的值的影响，以这种方式，处理能力可以在经济上再次保持较低。
根据本发明的磁共振系统包括磁体单元、梯度单元、高频单元和控制装置，该控制装置旨在利用参数值确定单元来执行根据本发明的方法。
当在控制装置上执行根据本发明的方法时，根据本发明的计算机程序在控制装置上实现它。
这里，计算机程序也可以以计算机程序产品的形式提供，该计算机程序产品可以直接加载到控制装置的存储器中，该计算机程序具有程序编码装置，以在计算机程序产品在处理系统的处理器单元中被执行时执行根据本发明的方法。
根据本发明的电子可读数据介质包括存储在其上的电子可读控制信息，其中，控制信息包括至少一个根据本发明的计算机程序，并采取这样的形式：使得在数据介质用于磁共振系统的控制装置中时，所述控制信息执行根据本发明的方法。
与该方法有关的指定的优点和实施方式也类似地适用于磁共振系统、计算机程序和电子可读数据介质。
附图说明
根据下面描述的示例性实施方式和附图，本发明的其他优点和细节将变得明显。给出的示例不代表对本发明的任何限制。在附图中：
图1示出了根据本发明的方法的示意流程图，
图2示出了体素时间序列的可能的捕获方案，以及
图3示出了根据本发明的示意性示出的磁共振系统。
具体实施方式
图1是根据本发明的方法的示意流程图，该方法用于使用磁共振指纹(MRF)技术来确定检查对象的体素中的参数值。
这里，借助于MRF方法建立至少一个体素时间序列(BZS)，根据所述至少一个体素时间序列(BZS)来确定由相应体素中表示的检查对象的位置处的至少一个参数的相应值(P，P')(框101)。这里，特别地建立体素时间序列BZS，使得体素时间序列的测量信号特性显示对要被确定的参数的依赖性，也就是说要被确定的参数的值的变化将导致特性中测量的信号中的至少一些信号的变化。
建立的体素时间序列BZS的相应信号特性的至少一个片段与第一比较信号特性D的对应片段的第一信号比较被执行，以便确定要被确定的参数中的至少第一参数的相应值P(框103)。
在第一信号比较期间，例如可以确定来自包括主磁场B0、发射场B1、横向弛豫时间T2和纵向弛豫时间T1的参数组中的至少一个参数，特别是所有参数。
出于该目的，该方法可以用于例如MRF方法的背景中，用于同时确定T1、T2、B1和B0，其中参数B1+、B0、T1和T2的值在每个体素处的单个连续测量中被识别。下面将参考图2来解释用于这种MRF方法的可能捕获方案的示例。在所示的捕获方案中，其中使用TrueFISP型序列成分，这是由于这些成分的特征在于B0参数的良好编码和高信号效率。
基于第一比较信号特性D和在第一信号比较中确定的值P生成另外比较信号特性(框105)。
这里，对另外比较信号特性D'的生成可以包括：在每种情况下针对第一另外参数的值，在第一信号比较中确定的相应体素中的第一参数的值的相应数量的预定的假定可能变化(V)上，对第一信号特性进行求和。
这里，可以基于被研究对象的先验知识预先确定预定的假定可能变化。
例如，在第一信号比较中，可以确定纵向弛豫时间T1的值或横向弛豫时间T2的值，并且第一另外参数可以是纵向弛豫时间T1或横向弛豫时间T2，以及第二另外参数是纵向弛豫与横向弛豫之间的比率。这里，基于检查对象中的组织类型的先验知识和其参数的典型比率，例如纵向弛豫时间比横向弛豫时间T1/T2，组织类型可能存在的比率T1/T2可以被预先确定为假定可能的变化。因此，所生成的另外比较信号特性可以作为合成维度从而改变可能的比率T1/T2，其中T1或T2的值已经由第一信号比较确定。如果在另一信号比较中，将T1或T2的值同时确定为第一另外参数，则也可以在适当的情况下迭代地改善结果。一旦另一信号比较已经确定了“最佳匹配”并因此确定了比率T1/T2的值和T1或T2的值，则还确定了另一参数——T2或T1的值。
因此，例如在第一信号比较期间，可以确定纵向弛豫时间T1的值，并且对于第一另外参数和第二另外参数可以是横向弛豫时间T2以及纵向弛豫与横向弛豫之间的比率。
例如，在第一信号比较期间，还可以确定横向弛豫时间T2的值，并且对于第一另外参数和第二另外参数是纵向弛豫时间T1以及纵向弛豫与横向弛豫之间的比率。
此外，例如在第一信号比较期间，还可以确定横向弛豫时间T2的值，并且对于第一另外参数和第二另外参数可以是横向弛豫时间T2以及纵向弛豫与横向弛豫之间的比率。在这种情况下，纵向弛豫时间T1的值也可以从确定的T2的值和纵向弛豫与横向弛豫之间的比率中找到。
此外，例如在第一信号比较期间，这里还可以确定纵向弛豫时间T1的值，并且对于第一另外参数和第二另外参数是纵向弛豫时间T1和纵向弛豫与横向弛豫之间的比率。在这种情况下，横向弛豫时间T2的值也可以从确定的T1的值和纵向弛豫与横向弛豫之间的比率中找到。
另外或作为替选，预定的假定可能变化可以对应于统计分布，其中该数量对应于统计分布的不同的可能假定宽度的数量。在第二另外参数依赖于第一参数，其中所确定的值形成用于生成另外比较信号特性D'的基础时，这主要是有用的，因为以这种方式由所述依赖性产生的第一参数对第二另外参数的效应可以被研究。
可以例如通过对第一比较信号特性D求和来识别另外比较信号特性D'。这可以特别地在每种情况下针对第一另外参数的固定值的这些比较信号特性D中的与第一另外参数例如T2的值的固定值相对应的比较信号特性D来完成，第一参数的值中的位于由所确定的第一参数的值和统计分布的相应宽度预先确定的值范围内的那些值被求和。在这种情况下，在第一信号比较中确定的第一参数的值可以例如位于相应宽度的中心。相应的求和还可以包括根据统计分布加权。因此，对要求和的比较信号特性D进行求和，使得根据统计分布进行加权。
这里，统计分布可以特别地是高斯分布或洛伦兹分布或真的均匀分布等。可以例如通过实验或根据基于不同参数之间的关系的先验知识的评估来选择合适的分布。
在第一信号比较期间，例如可以确定局部主磁场B0的值，并且第一另外参数可以是横向弛豫时间T2，并且第二另外参数可以是依赖于体素中占优势的相分散的参数，特别是相分散本身，由主磁场的变化引起的横向弛豫时间T2'，或有效横向弛豫时间T2*，或其倒数之一。
例如，预定的假定可能变化因此可以是以大约第一信号比较中确定的值例如主磁场B0为中心的统计分布的不同宽度。通过另一信号比较确定的在体素中占优势的这种宽度的值对应于体素中的相分散(偏共振效应)。
由主磁场的局部变化引起的横向磁化T2'的衰减的局部值可以通过由另一信号比较确定的体素中占优势的这种宽度的值来描述，因为由主磁场的局部变化引起的横向磁化衰减的倒数与指定为例如半峰全宽FWHM的所述宽度成比例：
在每种情况下，对所识别的体素时间序列BZS的相应信号特性的至少一个片段与在框105中生成的另外比较信号特性D'的对应片段进行另一信号比较，以便分别确定要确定的参数中的至少第一另外参数和第二另外参数P'的值(框107)。
由于通过使用在一定程度上形成子字典地所生成的比较信号特性的另一信号比较同时确定第一另外参数和第二另外参数的值，因此第一另外参数和第二另外参数的值被整体确定，因此可以避免由(先前未知的)其他另外参数的值对一个另外参数的值的伪造。
针对相应体素确定的要确定的参数的值P和P'可以例如以参数图的形式被输出，或者例如被存储以供之后使用。
因此，所描述的方法可以用于例如最小化由例如(T2*)相分散引起的MRF参数图的伪造。因此，在生成体素时间序列时，还可以使用其中信号响应受T2*效应影响的序列类型——例如在TrueFISP序列中的情况——其生成明显依赖于阻带(或暗带)内的这样的偏共振效应的信号，这例如可能导致错误的T2值。因此，该方法可以与任何MRF方案一起使用，其中所使用的序列类型显示信号的时变T2*依赖性。
特别地，为了避免由于相分散效应的参数图的伪造，该方法可以例如在第一信号比较(框103)中确定主磁场B0的值，并且在适当的情况下还确定参数T1、T2和B1的值，其中，第一比较信号特性D各自包括值在第一信号比较中确定的参数的维度，但不包括相分散的参数的维度。由于已知相分散对根据MRF信号比较来确定参数T1、B0和B1的值的影响很小(参见引用的Chiu等人的文章)，因此之后对这些参数确定的值可以被保留。在另一信号比较(框107)中，使用专门出于该目的生成的另一比较信号特性(框105)，可以(更精确地)确定参数T2的值，并且可以将相分散(首次)确定为独立参数。生成另外比较信号特性D'(框105)，使得它们(如先前的第一比较信号特性D)在T2并考虑通过改变的第二另外参数对信号的可能改变的情况下，在改变的相分散的情况下包括第一另外参数的单独维度。出于该目的，在生成另外比较信号特性D'时，可以针对每个体素计算表示由相分散导致的信号变化的子字典，其中例如针对T2维度的每个值，对第一比较信号特性求和，在B0维度中加权，这给出了大致与由在不同频率处的自旋进动产生的比较信号特性相对应的比较信号特性。加权由例如体素内的偏共振的假定分布(例如，高斯、洛伦兹、均匀分布等)造成。这里，分布特别地以第一信号比较中确定的相应B0值的中心频率为中心。不同的分布宽度导致不同(体素内)相分散值的信号特性，这可以因此在确定参数值时被考虑在内。
这里描述的方法还可以用作MRF方法中的部分步骤，该MRF方法具有所建立的体素时间序列与比较信号特性之间的多个信号比较，用于逐步确定参数值。
图2示意性地示出了用于捕获至少一个体素时间序列的捕获方法，由于它可以用于根据本发明的方法。所示的示例示出了用于建立体素时间序列的捕获方法，其中使用来自序列类型TrueFISP(具有稳态自由进动的真实快速成像)、FISP(具有稳态自由进动的快速成像)和FLASH(快速小角度激发)的组的三种不同的序列类型。序列类型的性质特别是指各个序列类型对特定于组织和/或特定于测量的参数的变化的灵敏度。例如，FISP序列对主磁场B0的变化几乎没有灵敏度，而trueFISP序列对主磁场B0的变化更敏感。FLASH序列和FISP序列对照射的高频场B1+中的局部变化敏感。
在图2所示的示例中，时间序列中捕获的图像数据集的相应序列号在轴26上被指示，并且不同的变量在轴27上分别指示。作为第一变量，从原点处0°以°为单位至轴点28处90°的翻转角被指示。在所示的示例中，轴26从图像数据集1进行到图像数据集3000。
3000个图像数据集被分布在十二个片段29、30、31、32、33、34、35、36、37、38和39上。
在第一片段29中，通过200个图像数据集的曲线40指示在捕获期间使用的翻转角，其中，FISP序列可以用于片段29中的捕获。如关于图1所描述的，一旦施加具有特定翻转角的RF激励脉冲，则捕获了完整的图像数据集，并且然后具有下一个翻转角的下一个RF激励脉冲被施加，并捕获另一图像数据集。在片段29中，图2示出了与sin2半曲线相对应的翻转角分布。最大翻转角可以是例如24°，并且可以使用恒定相位。
仅作为示例，针对第100个图像数据集绘制线41。对应的翻转角是曲线40的最大翻转角。
在第二片段30中，在所示的示例中，使用不同的序列类型例如TrueFISP序列获取400个图像数据集。在这种情况下，使用根据曲线42和曲线43的翻转角。在曲线42的情况下，这些达到45°，而在曲线43的情况下达到72°。
仅作为示例，对于片段30，针对第四百个图像数据组再次对翻转角绘制线44。在这种情况下，翻转角为1°。
片段30中的特定特征是使用两个不同的相位周。在通过曲线42的翻转角时，使用00相位周或没有相位周，并且在通过曲线43时，使用180°的相位周。00相位周表示固定相位。
在下一个片段31中，在曲线45中，使用另一序列类型例如FLASH序列对450个图像数据集的捕获指示翻转角。这些比FISP或TrueFISP序列中更小，并达到6°。它们的分布也是sin2分布。
除了改变翻转角之外，当重复FLASH序列时，应用相位周以产生RF破坏。如上面描述的，这使相位增加了117°的倍数。
在片段29、30和31中使用的不同序列类型的序列一起形成框55。这在图2中总共使用了三次。这里，仅参考序列的类型而不是图像数据集的数量或翻转角曲线。
在片段32中，使用框55的第一序列类型即例如FISP序列再次捕获200个图像数据集。如同在片段29中，相位是恒定的，但最大翻转角是45°。使用的翻转角位于曲线46上。
这之后是片段33中使用框55的第二序列类型例如TrueFISP序列获取的200个图像数据集。这里，使用90°相位周，并且最大翻转角为50°。翻转角被指示在曲线47上。
如同在片段31中，在片段34中，使用框55的第三序列类型例如FLASH序列捕获接下来约450个图像数据集。曲线48显示具有最大值14°的sin2分布。
第三次使用框55的第一序列类型例如FISP序列时，片段35中的曲线49达到72°并且示出了高频脉冲19的翻转角。在这个过程中，相位也是恒定的。
在利用框55的第二序列类型例如TrueFISP序列获取另外200个图像数据集时，使用270°相位周。在片段36中的曲线50中绘制的翻转角达到65°。
使用框55的第三序列类型例如FLASH序列捕获片段37中的接下来约450个图像数据集。曲线51显示延伸至至多20°的翻转角曲线，再次以sin2分布。
在最后一个片段38中，存在两条曲线52和53，用于利用框55的第一序列类型例如FISP序列捕获图像数据集。这些再次表示翻转角曲线。如上在前面的片段中，例如在利用FISP序列捕获测量数据中使用恒定相位。
通过其在具有不同灵敏度的不同序列类型的片段中的呈现，所示的示例允许在匹配步骤中确定多个特定于组织和特定于测量的参数；特别地，可以确定参数T1、T2和B0以及B1+。然而，所示的示例不应被认为是限制性的。原则上，也可以仅借助于一种序列类型建立体素时间序列。然而，使用不同特征的多种不同序列类型，特别是参考对特定于组织和/或特定于测量的参数的相应灵敏度，增加了可以根据以这种方式捕获的体素时间序列确定的参数的数量，和/或以与所使用的序列类型的特征相对应的方式确定的参数的值的质量。
所示的利用一种序列类型捕获的图示数据集的数量和所示的翻转角曲线也应该仅被视为示例。
为了说明补偿过程，在检查对象中示出特定主磁场值B0的空间分布作为示例。在记录根据其已经重建图像数据的测量数据期间，应将每个体素的信号特性作为体素时间序列与比较信号特性进行比较，以至少确定主磁场B0的局部值。在检查对象中生成基本磁场B0“从上到下”的线性特性以阐明效果。
左侧所示的检查对象中的主磁场值(B0图)的空间分辨分布对应于通过体素时间序列与比较信号特性的信号比较逐个体素确定的主磁场值，其中，用于信号比较的比较信号特性是字典的比较信号特性，其比较信号特性覆盖比检查对象中出现的B0值更小的B0值范围。因此，左侧所示的B0图覆盖了例如在-40到+40的范围内的B0值，其以任何期望单位由体素时间序列下面的测量数据的记录类型和对应的比较信号特性(字典)预定。
中心所示的差分B0图对应于左侧所示的B0图的差分图像，其中B0图具有粗略空间分辨率，其由从体素时间序列的特定B0敏感部分导出的B0值组成，并且因此以无限小的步长覆盖B0值。可能使用平滑操作，通过确定已经在不同回波时间处记录的体素时间序列的部分中信号的相位差来根据体素时间序列的所述部分确定粗略分辨的B0图，其中，利用具有FISP序列的性质的序列类型已经记录了测量数据。
如果左侧所示的所述粗略空间分辨的B0图和由信号比较产生的B0图的差分B0图被四舍五入为1/TR的倍数，则这得到对应于1.TR的倍数的片段(在所示的示例中，可以识别出五个这样的片段)。如果从受到伪影的B0图(左)中减去所述四舍五入的结果，则获得补偿的、解卷积的B0图，如右边所示。右边所示的B0图是对在测量期间中人为生成的主磁场B0的线性特性的很好的表示，即使通过利用粗略B0图在较大值范围内所描述的处理，在以与左侧所示的B0图中相同的单位的从-200至+200的范围中示出的示例中。
这里，使用这样的事实：补偿数据已经可以从没有信号比较的情况下使用的体素时间序列的信号中获得，例如，通过对要被确定的参数的局部至的至少粗略分辨的确定，特别是对特定于测量的参数，例如主磁场B0，该补偿数据可用于补偿，例如，对通过信号比较确定的参数值进行解卷积。尽管由于参数值的局部分布，这样的粗略分辨的参数图的分辨率不是直接足够的，但是粗略分辨的参数图的参数值可以用作补偿数据，以便补偿通过信号比较已经确定的更高分辨的参数，例如，如上面描述的示例中的解卷积。
因此，在所创建的体素时间序列的高分辨率中，针对体素时间序列确定的参数的补偿值被确定，这可以改善先前确定的相应参数的值并因此补偿伪像，其中，与利用纯信号比较可以实现的值的范围相比，甚至可以扩大在确定补偿的参数值时可以实现的值的范围。
利用可以确定参数T1、T2、B0和B1+的MRF记录方法，用于利用局部参数的多级确定来确定参数T1、T2、B0和B1的局部参数值的根据本发明的方法的示例性可能序列可以如下出现。
在建立的体素时间序列与比较信号特性之间的第一信号比较，其中特别地将体素时间序列的整个信号特性与比较信号特性的相应整个特性进行比较之后，可以确定所有参数T1、T2、B0和B1+的局部值(图)。
图3示意性地示出了根据本发明的磁共振系统1。这包括用于生成主磁场的磁体单元3、用于生成梯度场的梯度单元5、用于发射和接收高频信号的高频单元7以及旨在执行根据本发明的方法的控制装置9。在图3中，磁共振系统1的这些部分单元仅以粗略示意方式示出。特别地，高频单元7可以包括多个子单元，例如多个线圈，例如示意性示出的线圈7.1和7.2或更多线圈，其可以旨在仅用于发射高频信号或仅用于接收高频信号，或者用于两者。
为了研究检查对象U，例如患者或实际上是体模，可以将所述对象引入到桌子L上的磁共振系统1的测量腔中。切片S代表从中捕获测量数据的检查对象的示例性目标体积。
控制装置9用于控制磁共振系统，并且可以特别地通过梯度控制器5'控制梯度单元5，并且通过高频发射/接收控制器7'控制高频单元7。在这种情况下，高频单元7可以包括多个信道，在信道上可以发射或接收信号。
高频单元7与其高频发射/接收控制器7'一起负责生成并照射(发射)高频交变场，以便操纵检查对象U的要操纵区域中的自旋(例如要测量的切片S中)。为此，高频交变场(也被指定为B1场)的中心频率被设定至最大可能程度，使得它接近于要被操纵的自旋的共振频率。中心频率与共振频率的偏差被称为偏共振。为了生成B1场，使用高频发射/接收控制器7'控制的电流被施加至高频单元7中的HF线圈。
此外，控制装置9包括参数值确定单元15，通过参数值确定单元15可以执行根据本发明的信号比较以确定参数值。控制装置9作为整体旨在执行根据本发明的方法。
包括在控制装置9中的处理器单元13旨在执行所需测量和确定所需的所有处理操作。出于该目的所需的或在此期间识别的结果和中间结果可以被存储在控制装置9的存储器单元S中。这里，所示的单元不一定被理解为物理上分离的单元，而仅仅表示细分成有意义的单元。然而，这也可以例如在更少或甚至仅在一个物理单元中实现。
控制命令可以被发送至磁共振系统以及/或者结果例如图像数据可以例如由用户通过磁共振系统1的输入/输出装置E/A显示在控制装置9上。
本文中描述的方法还可以以计算机程序产品的形式提供，该计算机程序产品包括程序并且如果其在控制装置9上执行，则在控制装置9上实现所描述的方法。同样可以提供将电子可读控制信息存储在其上的电子可读数据介质26，其中，该控制信息至少包括刚刚描述的类型的计算机程序产品，并且采用这样的形式：在数据介质26在磁共振系统1的控制装置9中使用时，它执行所描述的方法。

Claims (15)

1.一种用于使用磁共振指纹MRF技术来确定检查对象的体素中的参数值的方法，包括以下步骤：
借助于MRF捕获方法建立至少一个体素时间序列(BZS)，其中，根据所述至少一个体素时间序列由相应体素表示的所述检查对象中的位置处的至少一个参数的相应值(P，P')将被确定，
执行所建立的体素时间序列(BZS)的相应信号特性的至少一个片段与第一比较信号特性(D)的对应片段的第一信号比较(103)，以确定将被确定的参数中的至少第一参数的相应值(P)，
基于所述第一比较信号特性和所述第一信号比较中确定的值来生成另外比较信号特性(D')，
执行所建立的体素时间序列(BZS)的相应信号特性的至少一个片段与所生成的另外比较信号特性(D')的对应片段的另一信号比较，以分别确定将被确定的参数中的至少第一另外参数和第二另外参数的值，
输出针对所述相应体素确定的将被确定的参数的值。
2.根据权利要求1所述的方法，其中，建立所述体素时间序列(BZS)，使得所述体素时间序列的信号特性显示对将被确定的所述参数的依赖性。
3.根据权利要求1或2所述的方法，其中，所述第二另外参数对所述第一另外参数有影响。
4.根据权利要求1或2所述的方法，其中，所述第二另外参数对所述第一参数有依赖性，其中所述第一参数的所确定的值形成用于生成所述另外比较信号特性(D')的基础。
5.根据权利要求1或2所述的方法，其中，在所述另一信号比较期间不再确定显示不依赖于所述第二另外参数的参数的在所述第一信号比较期间确定的值。
6.根据权利要求1或2所述的方法，其中，对所述另外比较信号特性(D')的生成包括：在每种情况下针对所述第一另外参数的值，在所述第一信号比较中确定的相应体素中的第一参数的值的相应数量的预定假定可能变化(V)上，对所述第一比较信号特性求和。
7.根据权利要求6所述的方法，其中，基于关于所检查对象的先验知识来预先确定所述预定假定可能变化。
8.根据权利要求6所述的方法，其中，所述预定假定可能变化对应于统计分布，其中所述预定假定可能变化的数量对应于不同的可能假定宽度的数量。
9.根据权利要求1、2、7和8中的任一项所述的方法，其中，在所述第一信号比较期间，确定来自包括主磁场、发射场、横向弛豫时间T2和纵向弛豫时间T1的一组参数中的至少一个参数。
10.根据权利要求1、2、7和8中任一项所述的方法，其中，在所述第一信号比较期间，确定局部主磁场的值，并且所述第一另外参数是横向弛豫时间T2，并且所述第二另外参数是取决于相应体素中占优势的相分散的参数。
11.根据权利要求1、2、7和8中任一项所述的方法，其中，在所述第一信号比较期间，确定纵向弛豫时间T1的值或横向弛豫时间T2的值，并且所述第一另外参数是所述纵向弛豫时间T1或所述横向弛豫时间T2，并且所述第二另外参数是所述纵向弛豫时间T1与所述横向弛豫时间T2之间的比率。
12.根据权利要求1、2、7、8和11中任一项所述的方法，其中，值在所述另一信号比较中确定的第二另外参数不在所述第一信号比较中确定。
13.一种磁共振系统(1)，包括：磁单元(3)、梯度单元(5)、高频单元(7)以及具有高频发射/接收控制器(7')和参数值确定单元(15)的控制装置(9)，其中，所述控制装置(9)旨在在所述磁共振系统(1)上执行根据权利要求1至12之一所述的方法。
14.一种磁共振系统(1)的控制装置(9)的存储器，计算机程序能够被直接加载到所述存储器中，所述计算机程序具有编程装置，以在所述磁共振系统(1)的所述控制装置(9)中执行所述程序时执行根据权利要求1至12之一所述的方法的步骤。
15.一种电子可读数据介质，其存储有电子可读控制信息，其中，所述控制信息包括至少一个加载到根据权利要求14所述的存储器中的计算机程序，并采用这样的形式：在所述数据介质在磁共振系统(1)的控制装置(9)中使用时，所述控制信息执行根据权利要求1至11之一所述的方法。
CN201910540740.7A 2018-06-21 2019-06-21 磁共振指纹方法 Active CN110623666B (zh)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
US201862688013P 2018-06-21 2018-06-21
US62/688,013 2018-06-21

Publications (2)

Publication Number Publication Date
CN110623666A CN110623666A (zh) 2019-12-31
CN110623666B true CN110623666B (zh) 2023-07-21

Family

ID=66589395

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910540740.7A Active CN110623666B (zh) 2018-06-21 2019-06-21 磁共振指纹方法

Country Status (3)

Country Link
US (1) US11639978B2 (zh)
EP (1) EP3588121B1 (zh)
CN (1) CN110623666B (zh)

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103376432A (zh) * 2012-04-20 2013-10-30 西门子公司 快速地位置分辨地确定检查区域内磁共振弛豫参数的方法
WO2014111266A2 (en) * 2013-01-17 2014-07-24 Ruprecht-Karls-Universität Heidelberg Synthetic parameter-weigthed images based on relaxometry (spare)
DE102018209584A1 (de) * 2018-06-14 2019-12-19 Siemens Healthcare Gmbh Magnetresonanz-Fingerprinting-Verfahren

Family Cites Families (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10261154B2 (en) * 2014-04-21 2019-04-16 Case Western Reserve University Nuclear magnetic resonance (NMR) fingerprinting tissue classification and image segmentation
US9897675B2 (en) * 2014-05-28 2018-02-20 Case Western Reserve University Magnetic resonance fingerprinting (MRF) with simultaneous multivolume acquisition
DE102014217284A1 (de) 2014-08-29 2016-03-03 Siemens Aktiengesellschaft Verfahren zu einer Magnetresonanz-Untersuchung eines Untersuchungsobjekts
US10209334B2 (en) * 2015-04-08 2019-02-19 The General Hospital Corporation System and method for controlling partial volume effects in magnetic resonance fingerprinting
US10492710B2 (en) * 2016-05-06 2019-12-03 Case Western Reserve University System and method for spectroscopic magnetic resonance fingerprinting
DE102016217675A1 (de) * 2016-09-15 2018-03-15 Siemens Healthcare Gmbh Verbesserte Erzeugung von Bildpunkt-Zeit-Serien eines Untersuchungsobjektes mittels Magnetresonanztechnik
EP3382416A1 (en) * 2017-03-30 2018-10-03 Koninklijke Philips N.V. Selection of magnetic resonance fingerprinting dictionaries for anatomical regions
US10859654B2 (en) * 2018-03-13 2020-12-08 Case Western Reserve University Systems and methods for acceleration of dictionary generation and matching in perfusion analysis
CN110346741A (zh) 2018-04-06 2019-10-18 凯斯西储大学 使用多个脉冲序列类型的磁共振指纹的系统和方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103376432A (zh) * 2012-04-20 2013-10-30 西门子公司 快速地位置分辨地确定检查区域内磁共振弛豫参数的方法
WO2014111266A2 (en) * 2013-01-17 2014-07-24 Ruprecht-Karls-Universität Heidelberg Synthetic parameter-weigthed images based on relaxometry (spare)
DE102018209584A1 (de) * 2018-06-14 2019-12-19 Siemens Healthcare Gmbh Magnetresonanz-Fingerprinting-Verfahren

Also Published As

Publication number Publication date
US11639978B2 (en) 2023-05-02
US20190391220A1 (en) 2019-12-26
EP3588121B1 (en) 2022-09-14
CN110623666A (zh) 2019-12-31
EP3588121A1 (en) 2020-01-01

Similar Documents

Publication Publication Date Title
US10663545B2 (en) Method and apparatus for low-artifact magnetic resonance fingerprinting scan
JP6557710B2 (ja) 核磁気共鳴（ｎｍｒ）フィンガープリンティング
US10335049B2 (en) Magnetic resonance fingerprinting with steady state precession (MRF-FISP)
Ben‐Eliezer et al. Rapid and accurate T2 mapping from multi–spin‐echo data using Bloch‐simulation‐based reconstruction
US9766313B2 (en) MR imaging using apt contrast enhancement and sampling at multiple echo times
US10466322B2 (en) MR imaging using multi-echo k-space acquisition
EP2798364B1 (en) Mr imaging with suppression of flow artefacts
US11373392B2 (en) System and method for magnetic resonance fingerprinting using a plurality of pulse sequence types
CN107440719B (zh) 用于显示定量磁共振图像数据的方法
US10564243B2 (en) System and method for phase cycling magnetic resonance fingerprinting (phc-MRF)
US11874355B2 (en) Magnetic resonance fingerprinting method
US10921406B2 (en) Magnetic resonance fingerprinting method and apparatus
CN109791186B (zh) 对磁共振指纹期间的b0偏共振场的直接测量
RU2538421C2 (ru) Картирование градиента восприимчивости
KR20170139037A (ko) 구획 교환을 포함하는 자기 공명 핑거프린팅(mrf)을 이용한 멀티 컴포넌트 복셀 분리
US20160116560A1 (en) Mr imaging with enhanced susceptibility contrast
US20140121492A1 (en) Mri with separation of different chemical species using a spectral model
EP3413070A1 (en) Dual-echo dixon-type water/fat separation mr imaging
JP7023954B2 (ja) プロペラｍｒイメージング
US20200033432A1 (en) Method and apparatus for identifying at least one material comprised in a voxel
CN110623666B (zh) 磁共振指纹方法
US11073585B2 (en) Low-field magnetic resonance fingerprinting
US20190361086A1 (en) Method and apparatus for implementing a magnetic resonance measurement that is insensitive to off-resonance
US10725132B2 (en) Method and apparatus for mitigating the effect of magnetization transfer in model-based magnetic resonance techniques
US10317490B2 (en) Method and magnetic resonance apparatus for determining dephasing factors in MR multiecho techniques

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