CN114076911B - 磁共振成像的方法及系统 - Google Patents

磁共振成像的方法及系统 Download PDF

Info

Publication number
CN114076911B
CN114076911B CN202110912132.1A CN202110912132A CN114076911B CN 114076911 B CN114076911 B CN 114076911B CN 202110912132 A CN202110912132 A CN 202110912132A CN 114076911 B CN114076911 B CN 114076911B
Authority
CN
China
Prior art keywords
value
acquisition
magnetic resonance
pulse sequence
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
CN202110912132.1A
Other languages
English (en)
Other versions
CN114076911A (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
Publication of CN114076911A publication Critical patent/CN114076911A/zh
Application granted granted Critical
Publication of CN114076911B publication Critical patent/CN114076911B/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/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
    • G01R33/56—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
    • G01R33/5602—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution by filtering or weighting based on different relaxation times within the sample, e.g. T1 weighting using an inversion pulse
    • 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/24—Arrangements or instruments for measuring magnetic variables involving magnetic resonance for measuring direction or magnitude of magnetic fields or magnetic flux
    • G01R33/246—Spatial mapping of the RF magnetic field B1
    • 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/448—Relaxometry, i.e. quantification of relaxation times or spin density
    • 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/483—NMR imaging systems with selection of signals or spectra from particular regions of the volume, e.g. in vivo spectroscopy
    • G01R33/4831—NMR imaging systems with selection of signals or spectra from particular regions of the volume, e.g. in vivo spectroscopy using B1 gradients, e.g. rotating frame techniques, use of surface 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/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/565—Correction of image distortions, e.g. due to magnetic field inhomogeneities
    • G01R33/5659—Correction of image distortions, e.g. due to magnetic field inhomogeneities caused by a distortion of the RF magnetic field, e.g. spatial inhomogeneities of the RF magnetic field

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Condensed Matter Physics & Semiconductors (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Engineering & Computer Science (AREA)
  • Health & Medical Sciences (AREA)
  • General Health & Medical Sciences (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Radiology & Medical Imaging (AREA)
  • Signal Processing (AREA)
  • Artificial Intelligence (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Optics & Photonics (AREA)
  • Magnetic Resonance Imaging Apparatus (AREA)

Abstract

本申请提供了一种磁共振成像的方法及系统。该方法及系统可以获取与对象相关的第一组磁共振成像数据，第一组磁共振成像数据是第一采集中对象达到第一T1加权水平时由磁共振扫描仪获取的，获取与对象相关的第二组磁共振成像数据，第二组磁共振成像数据是第二采集中当对象达到第二T1加权水平时由磁共振扫描仪获取的，第二T1加权水平和第一T1加权水平不同。该方法及系统还可以基于所述第一组磁共振成像数据和所述第二组磁共振成像数据，确定参考系数的目标值，参考系数与第一采集中的第一B1不均匀性和第二采集中的第二B1不均匀性相关联。

Description

磁共振成像的方法及系统
优先权声明
本申请要求于2020年08月11日递交的美国申请US16/989,905的优先权，其内容以引用的方式被包含于此。
技术领域
本申请一般涉及磁共振成像(MRI)，更具体地，涉及MRI中定量测量的方法和系统。
背景技术
磁共振成像系统利用强磁场和射频(RF)技术广泛应用于医学诊断和/或治疗。例如，通过MRI系统，可以对对象进行各种测量，例如纵向弛豫时间(T1)测量、横向弛豫时间(T2)测量、质子密度(PD)测量，以提供疾病诊断和/或治疗的基础。
发明内容
根据本申请的一个方面，提供了一种MRI系统。该系统可以包括包含一组指令的至少一个存储设备，以及配置成与至少一个存储设备通信的至少一个处理器。当执行指令时，所述至少一个处理器可被配置为指示系统执行以下操作。所述至少一个处理器可被配置为获取与对象相关的第一组MRI数据，该第一组MRI数据是第一采集中所述对象达到第一T1加权水平时由磁共振(MR)扫描仪获取的，根据一个或多个第一值实施所述第一采集，所述一个或多个第一值中的每个第一值对应于一个或多个脉冲序列参数中的一个脉冲序列参数；获取与所述对象相关的第二组MRI数据，该第二组MRI数据是第二采集中当所述对象达到第二T1加权水平时由所述MR扫描仪获取的，该第二T1加权水平和所述第一T1加权水平不同，根据一个或多个第二值实施第二采集，所述一个或多个第二值中的每个第二值对应于所述一个或多个脉冲序列参数中的一个脉冲序列参数，所述一个或多个脉冲序列参数中的至少一个脉冲序列参数的第一值和第二值不同；以及基于所述第一组MRI数据和所述第二组MRI数据，确定参考系数的目标值，该参考系数与所述第一采集中的第一B1不均匀性和所述第二采集中的第二B1不均匀性相关联。
在一些实施例中，所述对象可以是受试体的物理点，并且所述参考系数可以测量在所述第一采集中所述物理点处的第一信号强度与在所述第二采集中所述物理点处的第二信号强度之间的差别。
在一些实施例中，所述参考系数可以为所述第一信号强度与所述第二信号强度之比或所述第二信号强度与所述第一信号强度之比。
在一些实施例中，所述第一组MRI数据和所述第二组MRI数据可以是在所述受试体的MRI扫描中获取的。并且所述基于所述第一组MRI数据和所述第二组MRI数据，确定参考系数的目标值，该参考系数与所述第一采集中的第一B1不均匀性和所述第二采集中的第二B1不均匀性相关联，可以包括：基于所述第一组MRI数据重建所述受试体的第一图像；基于所述第二组MRI数据重建所述受试体的第二图像；以及基于所述第一图像和所述第二图像确定所述参考系数的目标值。
在一些实施例中，基于所述第一图像和所述第二图像确定所述参考系数的目标值可以包括：在所述第一图像中识别与所述物理点对应的第一元素；在所述第二图像中识别与所述物理点对应的第二元素；以及基于所述第一元素的值和所述第二元素的值来确定所述参考系数的目标值。
在一些实施例中，所述操作还可以包括至少部分地基于所述参考系数的目标值来确定所述对象的T1值。
在一些实施例中，所述至少部分地基于所述参考系数的目标值来确定所述对象的T1值可以包括：基于所述一个或多个脉冲序列参数的所述一个或多个第一值和所述一个或多个第二值确定所述参考系数的参考值；基于所述参考值和所述参考系数的目标值，确定关于T1测量的校正系数；以及至少部分地基于所述校正系数来确定所述对象的T1值。
在一些实施例中，所述一个或多个脉冲序列参数的所述一个或多个第一值包括第一翻转角，以及所述一个或多个脉冲序列参数的所述一个或多个第二值包括与所述第一翻转角不同的第二翻转角。
在一些实施例中，所述至少部分地基于所述校正系数来确定所述对象的T1值可以包括基于所述第一翻转角、所述第二翻转角和所述参考系数的参考值来确定所述对象的初始T1值，以及通过根据所述校正系数校正所述初始T1值来确定所述对象的T1值。
在一些实施例中，所述至少部分地基于校正系数来确定所述对象的T1值可以包括：通过所述校正系数校正所述第一翻转角来确定校正后的第一翻转角；以及基于所述校正后的第一翻转角、所述第二翻转角和所述参考系数的目标值来确定所述对象的T1值。
在一些实施例中，所述至少部分地基于所述参考系数的目标值来确定所述对象的T1值可以是根据解析解或驱动平衡单脉冲观测(DESPOT)算法来进行的。
根据本申请的一个方面，提供了一种MRI方法。该方法可包括获取与对象相关的第一组MRI数据，该第一组MRI数据是第一采集中所述对象达到第一T1加权水平时由MR扫描仪获取的，根据一个或多个第一值实施所述第一采集，所述一个或多个第一值中的每个第一值对应于一个或多个脉冲序列参数中的一个脉冲序列参数；获取与所述对象相关的第二组MRI数据，该第二组MRI数据是第二采集中当所述对象达到第二T1加权水平时由所述MR扫描仪获取的，该第二T1加权水平和所述第一T1加权水平不同，根据一个或多个第二值实施第二采集，所述一个或多个第二值中的每个第二值对应于所述一个或多个脉冲序列参数中的一个脉冲序列参数，所述一个或多个脉冲序列参数中的至少一个脉冲序列参数的第一值和第二值不同；以及基于所述第一组MRI数据和所述第二组MRI数据，确定参考系数的目标值，该参考系数与所述第一采集中的第一B1不均匀性和所述第二采集中的第二B1不均匀性相关联。
在一些实施例中，所述对象可以是受试体的物理点，并且所述参考系数可以测量在所述第一采集中所述物理点处的第一信号强度与在所述第二采集中所述物理点处的第二信号强度之间的差别。
在一些实施例中，所述参考系数可以为所述第一信号强度与所述第二信号强度之比或所述第二信号强度与所述第一信号强度之比。
在一些实施例中，所述操作还可以包括至少部分地基于所述参考系数的目标值来确定所述对象的T1值。
在一些实施例中，所述至少部分地基于所述参考系数的目标值来确定所述对象的T1值可以包括基于所述一个或多个脉冲序列参数的所述一个或多个第一值和所述一个或多个第二值确定所述参考系数的参考值；基于所述参考值和所述参考系数的目标值，确定关于T1测量的校正系数；以及至少部分地基于所述校正系数来确定所述对象的T1值。
在一些实施例中，所述一个或多个脉冲序列参数的所述一个或多个第一值可以包括第一翻转角，以及所述一个或多个脉冲序列参数的所述一个或多个第二值可以包括与所述第一翻转角不同的第二翻转角。
在一些实施例中，至少部分基于校正系数来确定对象的T1值可以包括基于所述第一翻转角、所述第二翻转角和所述参考系数的参考值来确定所述对象的初始T1值，以及通过根据所述校正系数校正所述初始T1值来确定所述对象的T1值。
在一些实施例中，所述至少部分地基于所述校正系数来确定所述对象的T1值可以包括通过所述校正系数校正所述第一翻转角来确定校正后的第一翻转角，以及基于所述校正后的第一翻转角、所述第二翻转角和所述参考系数的目标值来确定所述对象的T1值。
根据本申请的另一方面，提供了一种非暂时性计算机可读介质。非暂时性计算机可读介质可包括用于MRI的一组指令，其中当由至少一个处理器执行时，该组指令可指示至少一个处理器实现方法。该方法可包括获取与对象相关的第一组MRI数据，该第一组MRI数据是第一采集中所述对象达到第一T1加权水平时由MR扫描仪获取的，根据一个或多个第一值实施所述第一采集，所述一个或多个第一值中的每个第一值对应于一个或多个脉冲序列参数中的一个脉冲序列参数；获取与所述对象相关的第二组MRI数据，该第二组MRI数据是第二采集中当所述对象达到第二T1加权水平时由所述MR扫描仪获取的，该第二T1加权水平和所述第一T1加权水平不同，根据一个或多个第二值实施第二采集，所述一个或多个第二值中的每个第二值对应于所述一个或多个脉冲序列参数中的一个脉冲序列参数，所述一个或多个脉冲序列参数中的至少一个脉冲序列参数的第一值和第二值不同；以及基于所述第一组MRI数据和所述第二组MRI数据，确定参考系数的目标值，该参考系数与所述第一采集中的第一B1不均匀性和所述第二采集中的第二B1不均匀性相关联。
本申请的一部分附加特性可以在以下描述中进行说明。通过对以下描述和相应附图的研究或者对实施例的生产或操作的了解，本申请的一部分附加特性对于本领域技术人员是明显的。本申请的特征可以通过实践或使用下文讨论的具体实施例中阐述的方法、手段和组合的各个方面来实现和达到。
附图说明
本申请将通过示例性实施例进行进一步描述。这些示例性实施例将通过附图进行详细描述。这些实施例是非限制性的示例性实施例，在这些实施例中，各图中相同的编号表示相似的结构，其中：
图1是根据本申请的一些实施例所示的示例性MRI系统的示意图；
图2是根据本申请的一些实施例所示的示例性MR扫描仪的框图；
图3是根据本申请的一些实施例所示的计算设备的示例性硬件和/或软件组件的示意图；
图4是根据本申请的一些实施例所示的移动设备的示例性硬件和/或软件组件的示意图；
图5是根据本申请的一些实施例所示的示例性处理设备的框图；
图6是根据本申请的一些实施例所示的用于确定对象的T1值的示例性过程的流程图；
图7是根据本申请的一些实施例所示的用于至少部分地基于参考系数的目标值来确定物理点的T1值的示例性过程的流程图；
图8是根据本申请的一些实施例所示的模拟实验的结果的示意图；
图9A到9C是根据本申请一些实施例所示的模拟实验的结果的示意图；以及
图10A到10D是根据本申请一些实施例所示的在水模上执行的MRI扫描的结果的示意图。
具体实施方式
在下面的详细描述中，通过示例阐述了许多具体细节，为了更清楚地说明本申请实施例的技术方案。然而，对于本领域技术人员显而易见的是可以在没有这些细节的情况下实施本申请。在其它情况下，为了避免不必要地使本申请的各方面变得晦涩难懂，已经在较高的层次上描述了众所周知的方法、过程、系统、组件和/或电路。对于本领域的普通技术人员来讲，对所披露的实施例做出的各种改变是显而易见的，并且在不偏离本申请的精神和范围的情况下，本申请中所定义的普遍原则可以适用于其他实施例和应用场景。因此，本申请不限于所示的实施例，而是被赋予与权利要求一致的最宽范围。
本申请中所使用的术语仅出于描述特定示例实施例的目的，而非限制性的。如本申请使用的单数形式“一”、“一个”及“该”同样可以包括复数形式，除非上下文明确提示例外情形。还应当理解，在本申请说明书中使用的术语“包括”、“包含”仅提示存在所述特征、整体、步骤、操作、组件和/或部件，但并不排除存在或添加一个或多个其它特征、整数、步骤、操作、组件、部件和/或其组合的情况。
应当理解，这里使用的术语“系统”、“设备”、“单位”、“模块”和/或“块”是以升序区分不同级别的不同组件、元素、零件、部分或组装的一种方法。然而，如果可以达到相同的目的，这些术语也可以被其他表达替换。
通常，这里使用的词语“模块”、“单元”或“块”是指包括在硬件或固件中的逻辑或者是软件指令的集合。本文描述的模块，单元或块可以被实现为软件和/或硬件，并且可以存储在任何类型的非暂时性计算机可读介质或其他存储设备中。在一些实施例中，可以编译软件模块/单元/块并将其链接到可执行程序中。应当理解，软件模块可以从其他模块/单元/块或从其自身调用，和/或可以响应于检测到的事件或中断而被调用。配置用于在计算设备上执行的软件模块/单元/块(例如，如图3所示的处理器310)可以设置在计算机可读介质上，例如光盘、数字视频光盘、闪存驱动器、磁盘或任何其他有形介质、或作为数字下载(最初可以以压缩或可安装的格式存储，在执行之前需要安装、解压缩或解密)。这里的软件代码可以被部分的或全部的储存在执行操作的计算设备的存储设备中，并应用在计算设备的操作之中。软件指令可以嵌入固件中，例如EPROM。还应当理解，硬件模块/单元/块可以包括在连接的逻辑组件中，例如门和触发器，和/或可以包括在可编程单元中，例如可编程门阵列或处理器。这里描述的模块/单元/块或计算设备功能可以实现为软件模块/单元/块，但是可以用硬件或固件表示。通常，这里描述的模块/单元/块指的是逻辑模块/单元/块，其可以与其他模块/单元/块组合或者分成子模块/子单元/子块，不管其物理组织或存储器件如何。该描述可适用于系统、引擎或其一部分。
可以理解的是，除非上下文另有明确说明，当单元、引擎、模块或块被称为在另一单元、引擎、模块或块“上”、“连接至”或“耦合至”另一单元、引擎、模块或块时，其可以直接在其它单元、引擎、模块或块上，与之连接或耦合，或与之通信，或者可能存在中间单元、引擎、模块或块。在本申请中，术语“和/或”可包括任何一个或多个相关所列条目的任意或所有组合。本申请中的术语“像素”和“体素”可交换地用于指代图像中的元素。本申请中的术语“图像”用于指各种形式的图像，包括二维图像、三维图像、四维图像等。
根据以下对附图的描述，本申请的这些和其它的特征、特点以及相关结构元件的功能和操作方法，以及部件组合和制造经济性，可以变得更加显而易见，这些附图都构成本申请说明书的一部分。然而，应当理解的是，附图仅仅是为了说明和描述的目的，并不旨在限制本申请的范围。应当理解的是，附图并不是按比例绘制的。
本申请中使用的流程图说明了根据本申请的一些实施例的系统实现的操作。应明确理解，流程图的操作可以不按顺序执行。相反，这些操作可以以相反的顺序实现，或者同时实现。此外，可以将一个或多个其他操作添加到流程图中。可以从流程图中删除一个或多个操作。
此外，虽然本申请中公开的系统和方法主要是关于在MRI系统中对对象(例如，患者、患者的物理点)进行测量来描述的。应该理解，这只是为了说明的目的。本申请的系统和方法可以应用于任何其他种类的医学成像系统。在一些实施例中，成像系统可以包括单模态成像系统和/或多模态成像系统。单模态成像系统可以包括例如MRI系统。多模态成像系统可以包括例如磁共振成像-计算机断层扫描(MRI-CT)系统、正电子发射断层扫描-磁共振成像(PET-MRI)系统、单光子发射计算机断层扫描-磁共振成像(SPECT-MRI)系统，数字减影血管造影-磁共振成像(DSA-MRI)系统等。
本申请的一个方面涉及用于MRI中定量测量的系统和方法。为了说明的目的，参考MRI中对对象的T1测量的实现来描述本申请。对象可以是生物的，也可以是非生物的。例如，对象可以包括患者、人造受试体等。又例如，对象可以包括患者的特定部分、器官、组织和/或物理点。例如，对象可包括头部、大脑、颈部、身体、肩部、手臂、胸部、心脏、胃、血管、软组织、膝盖、脚等，或其组合。在一些实施例中，对象可以是受试体(例如，患者)的物理点。如本文所使用的，受试体的物理点是指具有物理特性(例如重量、惯性、动量等)的部分，例如受试体的点或区域。在一些实施例中，受试体的物理点可对应于受试体的MRI图像的像素或体素，该MRI图像是通过对受试体的MRI扫描而生成的。物理点是指图像中对应的物体空间位置上的每个点(相应部位/位置)。应当注意，这并非旨在限制，而且可以应用所述系统和方法来对受试体进行其他定量测量(例如，T2测量)。
传统地，为了确定受试体的物理点的T1值，可以通过MR扫描仪将磁化恢复脉冲序列(例如，反转恢复脉冲序列、饱和恢复脉冲序列)施加在受试体上以从受试体获取回波信号，并且物理点的T1值可以通过将回波信号的信号强度与T1弛豫效应相关来确定。T1测量的另一方法可通过使用MR扫描仪对受试体实施具有不同脉冲序列参数值(例如，不同翻转角(FAs)或不同重复时间)的多个采集来执行。可以在多个获取中在物理点达到不同的T1加权水平时，获取多组MRI数据。因为物理点的T1值可能影响对应于不同脉冲序列参数值的一组MRI数据，所以可以基于一组MRI数据来确定对象的T1值。然而，采集的实现可能受到各种因素的影响，例如MR扫描仪的硬件限制、B0场(即，主磁场)的波长、B1场(即，射频场)的波长、受试体的介电特性和导电特性等。
由于这些因素，在采集中，施加于物理点的实际B1场可能偏离采集的预设B1场。另外，实际B1场可以具有不均匀的空间分布(即，实际B1场可以在受试体的不同物理点处具有不同的值)，这被称为B1不均匀性。B1不均匀性可能导致在物理点处的实际FA和采集的预设FA之间的差异，这可能影响在采集过程中采集的一组MRI数据和基于该组MRI数据进行的T1测量的精度。T1测量可能需要考虑采集中B1的不均匀性。例如，传统的T1测量方法可以在每次采集期间确定物理点处的实际B1场(或实际FA或实际FA与预设FA的比率)，并且基于采集中的实际B1场确定物理点的调整后的T1值。通常，传统的T1测量方法假设B1不均匀性对于不同的采集是一致的，例如，通过使用相同的校正系数来确定不同采集中的实际FAs。然而，不同的采集可能具有不同的B1不均匀性，并且在每个采集中确定的实际B1场可能仅仅是估计值。根据该假设确定的T1值可能具有低信噪比(SNR)。期望提供通过考虑到不同采集中不同的B1不均匀性来进行T1测量的系统和方法，从而提高T1测量的准确度。
根据本申请的一些实施例，当对象达到第一T1加权水平时，系统和方法可以获取与对象相关的第一组MRI数据，该第一组MRI数据是第一采集中对象达到第一T1加权水平时由MR扫描仪获取的，可以根据一个或多个第一值实施所述第一采集，该一个或多个第一值中的每个第一值对应于一个或多个脉冲序列参数中的一个脉冲序列参数。该系统和方法可以获取与对象相关的第二组MRI数据，该第二组MRI数据是第二采集中当对象达到第二T1加权水平时由MR扫描仪获取的，第二T1加权水平和第一T1加权水平不同，根据一个或多个第二值实施第二采集，该一个或多个第二值中的每个第二值对应于一个或多个脉冲序列参数中的一个脉冲序列参数，一个或多个脉冲序列参数中的至少一个脉冲序列参数的第一值和第二值不同。该系统和方法可基于第一组MRI数据和第二组MRI数据进一步确定参考系数的目标值，该参考系数与第一采集中的第一B1不均匀性和第二采集中的第二B1不均匀性相关联。
在一些实施例中，该系统和方法还可以基于参考系数的目标值来确定定量参数，例如对象的T1值。参考系数的目标值可用于消除或减少第一和第二B1不均匀性对确定对象的T1值的影响。与传统的T1测量方法相比，该系统和方法可以同时考虑B1的第一和第二不均匀性，提高了T1测量的精度和可靠性。
图1是根据本申请的一些实施例所示的示例性MRI(磁共振成像)系统100的示意图。如图1所示，MRI系统100可以包括MR扫描仪110、处理设备120、存储设备130、一个或多个终端140和网络150。在一些实施例中，扫描仪110、处理设备120、存储设备130和终端140可经由无线连接、有线连接或其组合彼此连接和/或通信。MRI系统100中的组件之间的连接可以是可变的。例如，MR扫描仪110可以通过网络150连接到处理设备120。又例如，MR扫描仪110可以直接连接到处理设备120。
MR扫描仪(磁共振扫描仪)110可被配置为扫描对象(或对象的一部分)以获取图像数据，例如与对象相关联的回波信号(或MRI数据)。例如，MR扫描仪110可以通过在对象上施加MR脉冲序列来检测多组MRI数据。在一些实施例中，MR扫描仪110可包括例如磁性体、梯度线圈、RF线圈等，如图2所述。在一些实施例中，根据磁性体的类型，MR扫描仪110可为永磁MR扫描仪、超导电磁铁MR扫描仪或电阻电磁铁MR扫描仪等。在一些实施例中，根据磁场强度，MR扫描仪110可以是高场MR扫描仪、中场MR扫描仪和低场MR扫描仪等。
为便于说明，图1中提供了包括X轴、Y轴和Z轴的坐标系。图1中所示的X轴和Z轴可以是水平的，Y轴可以是垂直的。如图所示，从面向MR扫描仪110正面的方向看，沿着X轴的正X方向可以是从MR扫描仪110的右侧到左侧的方向；沿图1所示的Y轴的正Y方向可以是从MR扫描仪110的下部到上部；图1所示的沿Z轴的正Z方向可以是指对象移出MR扫描仪110的扫描通道(或称为孔)的方向。可以在本申请的其他地方找到MR扫描仪110的更多描述。例如，参见图2及其描述。
处理设备120可以处理从MR扫描仪110、存储设备130和/或终端140获取的数据和/或信息。例如，处理设备120可以通过处理由MR扫描仪110收集的与对象有关的图像数据(例如，MRI数据)，来确定参考系数的目标值，其中，参考系数可与多个采集中的B1不均匀性相关联。又例如，处理设备120可以至少部分地基于参考系数的目标值对对象进行T1测量。在一些实施例中，处理设备120可以是单个服务器或服务器组。服务器组可以是集中式的，也可以是分布式的。在一些实施例中，处理设备120可以是本地的或远程的。例如，处理设备120可以经由网络150访问来自MR扫描仪110、存储设备130和/或终端140的信息和/或数据。又例如，处理设备120可以直接连接到MR扫描仪110、终端140和/或存储设备130，以访问信息和/或数据。在一些实施例中，处理设备120可以在云平台上实现。例如，云平台可以包括私有云、公共云、混合云、社区云、分布式云、云间云、多云等，或者其组合。在一些实施例中，处理设备120可以由具有如图3所述的一个或多个组件的计算设备300来实现。
存储设备130可以存储数据、指令和/或任何其他信息。在一些实施例中，存储设备130可以存储从MR扫描仪110、处理设备120和/或终端140获取的数据。在一些实施例中，存储设备130可以存储处理设备120可以执行的，或处理设备120可以用来执行本申请中描述的示例性方法的数据和/或指令。在一些实施例中，存储设备130可以包括大容量存储设备、可移动存储设备、易失性读写存储器、只读存储器(ROM)等或其组合。示例性大容量存储设备可以包括磁盘、光盘、固态驱动器等。示例性可移动存储设备可以包括闪存驱动器、软盘、光盘、存储卡、压缩磁盘、磁带等。示例性易失性读写存储器可以包括随机存取存储器(RAM)。示例性RAM可包括动态RAM(DRAM)、双数据速率同步动态RAM(DDR-SDRAM)、静态RAM(SRAM)、晶闸管RAM(T-RAM)、零电容RAM(Z-RAM)等。示例性ROM可包括掩模ROM(MROM)、可编程ROM(PROM)、可擦除可编程ROM(EPROM)、电可擦除可编程ROM(EEPROM)、光盘ROM(CD-ROM)、数字多功能磁盘ROM等。在一些实施例中，存储设备130可以在本申请其他地方描述的云平台上实现。
在一些实施例中，存储设备130可以连接到网络150以与MRI系统100中的一个或多个其他组件(例如，MR扫描仪110、处理设备120和/或终端140)通信。MRI系统100的一个或多个组件可以经由网络150访问存储设备130中存储的数据或指令。在一些实施例中，存储设备130可以是处理设备120或终端140的一部分。
终端140可被配置为实现用户与MRI系统100之间的用户交互。例如，终端140可从用户接收指令以使MR扫描仪110扫描对象。又例如，终端140可以从处理设备120接收处理结果(例如，与对象有关的定量参数的值或定量地图)，并将处理结果显示给用户。在一些实施例中，终端140可以连接到MR扫描仪110、处理设备120和/或存储设备130和/或与之通信。在一些实施例中，终端140可以包括移动设备140-1、平板计算机140-2、膝上型计算机140-3等或其组合。例如，移动设备140-1可以包括移动电话、个人数字助理(PDA)、游戏设备、导航设备、销售点(POS)设备、膝上型计算机、平板计算机、台式机等或其组合。在一些实施例中，终端140可包括输入设备、输出设备等。输入设备可包括可通过键盘输入的字母数字键和其他键、触摸屏(例如，具有触觉或触觉反馈)、语音输入、眼球跟踪输入、大脑监控系统输入、或任何其他类似的输入机制。通过输入设备接收到的输入信息可以经由例如总线发送到处理设备120以进行进一步处理。其他类型的输入设备可包括光标控制设备，例如鼠标、轨迹球或光标方向键等。输出设备可包括显示器、扬声器、打印机等，或其组合。在一些实施例中，终端140可以是处理设备120或MR扫描仪110的一部分。
网络150可以包括可促进MRI系统100的信息和/或数据交换的任何合适的网络。在一些实施例中，MRI系统100的一个或多个组件(例如，MR扫描仪110、处理设备120、存储设备130、终端140等)可经由网络150与MRI系统100的一个或多个其它组件传输信息和/或数据。例如，处理设备120可经由网络150从MR扫描仪110获取图像数据(例如，回波信号)。又例如，处理设备120可经由网络150从终端140获取用户指令。网络150可包括公用网络(例如，因特网)、专用网络(例如，局域网(LAN)、广域网(WAN))等)、有线网络(例如，以太网)、无线网络(例如，802.11网络、，Wi-Fi网络等)、蜂窝网络(例如，长期演进(LTE)网络)、帧中继网络、虚拟专用网络(“VPN”)、卫星网络、电话网络、路由器、集线器、交换机、服务器计算机等，或其组合。例如，网络150可以包括电缆网络、有线网络、光纤网络、电信网络、内部网、无线局域网(WLAN)、城域网(MAN)、公共电话交换网(PSTN)、蓝牙TM网络，紫蜂TM网络、近场通信(NFC)网络等，或其组合。在一些实施例中，网络150可以包括一个或多个网络接入点。例如，网络150可以包括有线和/或无线网络接入点，例如基站和/或因特网交换点，通过这些接入点，MRI系统100的一个或多个组件可以连接到网络150以交换数据和/或信息。
该描述旨在说明，而不是限制本申请的范围。许多替代、修改和变化对于本领域技术人员来说是显而易见的。本文描述的示例性实施例的特征、结构、方法和特性可以以各种方式组合以获取附加和/或替代示例性实施例。在一些实施例中，MRI系统100可以包括一个或多个附加组件和/或可以省略上面描述的一个或多个组件。另附加地或替代地，MRI系统100的两个或多个组件可以集成到单个组件中。例如，处理设备120可以集成到MR扫描仪110中。又例如，MRI系统100的组件可以被能够实现该组件的功能的另一组件替换。在一些实施例中，存储设备130可以是包括云计算平台的数据存储器，例如公共云、私有云、社区和混合云等。然而，这些变化和修改不脱离本申请的范围。
图2是根据本申请的一些实施例所示的示例性MR扫描仪110的框图。如图2所示，MR扫描仪110可以包括磁体220、梯度线圈230、RF线圈240和脉冲序列模块250。
磁体220可在扫描对象210的至少一部分期间产生静态磁场。磁体220可为各种类型，包括例如永磁体、超导电磁铁、电阻电磁铁等。
梯度线圈230可以在X方向、Y方向和/或Z方向上向主磁场提供磁场梯度。如本文所使用的，X方向、Y方向和Z方向可以表示坐标系(例如，与图1中描述的相同或类似的坐标系)中的X轴、Y轴和Z轴。例如，Z轴可以沿着磁体220的轴，X轴和Z轴可以形成水平面，并且X轴和Y轴可以形成垂直面。在一些实施例中，梯度线圈230可包括用于在X方向上向主磁场提供磁场梯度的X方向线圈、用于在Y方向上向主磁场提供磁场梯度的Y方向线圈，和/或用于向Z方向上的主磁场提供磁场梯度Z方向线圈。在一些实施例中，X方向线圈、Y方向线圈和/或Z方向线圈可以具有各种形状或配置。例如，Z方向线圈可以基于圆形(Maxwell)线圈来设计。又例如，X方向线圈和Y方向线圈可以基于鞍形(Golay)线圈配置来设计。
RF线圈240可以向被检查对象210发射RF脉冲信号和/或接收来自被检查对象210的回波信号。在一些实施例中，RF线圈240可以包括发射线圈和接收线圈。发射线圈可发射信号(例如，RF脉冲)，该信号可激发对象210中的原子核以提供共振。接收线圈可以接收从对象210发射的回波信号。在一些实施例中，RF发射线圈和RF接收线圈可以集成到一个相同的线圈中。在一些实施例中，RF线圈240可以是各种类型的，例如，包括正交检测(QD)正交线圈、相控阵线圈、特定元件频谱线圈等。在一些实施例中，RF线圈240可以是包括多个线圈单元的相控阵线圈，每个线圈单元可以独立地检测回波信号。
在一些实施例中，RF线圈240可用于检测由MR脉冲序列产生的信号。MR脉冲序列可以是各种类型，例如自旋回波(SE)脉冲序列、GRE脉冲序列、反转恢复(IR)脉冲序列、多回波MR脉冲序列、T1ρ准备脉冲序列、T2准备脉冲序列、扩散加权成像(DWI)脉冲序列、饱和恢复(SR)脉冲序列、稳态脉冲序列等。
在一些实施例中，MR脉冲序列可以由一个或多个脉冲序列参数定义，该参数包括例如MR脉冲序列的类型、施加MR脉冲序列的时间、MR脉冲序列的持续时间、MR脉冲序列中激励脉冲的翻转角，MR脉冲序列中RF脉冲的计数(或数量)、TR、重复计数、反转时间(TI)、MR脉冲序列中采集的计数(或数量)、b值、T1ρ-准备持续时间、T2准备持续时间、回波序列长度、回波间隔、速度编码(VENC)值等。如本文所用，激励脉冲的FA可指激励脉冲相对于主磁场的净磁化矢量的旋转。TR可指MR脉冲序列中两个重复且连续的RF脉冲之间的时间间隔(例如，SE脉冲序列中两个连续的激发RF脉冲之间的时间间隔，IR脉冲序列中两个连续的180°反转脉冲之间的时间间隔)。重复计数可指MR脉冲序列中的重复计数(或次数)。TI可指在IR脉冲序列中180°反转脉冲和随后的90°激发脉冲之间的时间跨度。b值可指反映DWI脉冲序列中扩散敏化梯度的强度和时间的因子。T1ρ准备持续时间是指T1ρ准备脉冲序列中自旋锁脉冲的持续时间。T2准备持续时间可指T2准备脉冲序列中T2准备脉冲的持续时间。
在一些实施例中，RF线圈240可以检测(或接收)与MR脉冲序列激发的一个或多个回波相对应的一个或多个回波信号。在一些实施例中，回波信号(或回波)可由一个或多个参数来定义，例如，回波信号类型(自旋回波、快速自旋回波(FSE)、快速恢复FSE、单次激发FSE、梯度回波、具有稳态进动的快速成像)、回波时间(TE)、回波信号强度、检测回波信号的线圈单元(例如，表示为检测回波信号的标识(ID)或线圈单元的序列号)、在其中检测回波信号的重复(例如，表示为重复序列号)、在其中检测回波信号的采集(例如，表示为采集序列号)等。TE可指应用激励射频脉冲和由激励射频脉冲激励的回波峰值之间的时间。
脉冲序列模块250可配置为在扫描对象210之前和/或期间定义与MR扫描仪110相关的参数和布置。在一些实施例中，与MR扫描仪110相关的参数可包括与MR脉冲序列相关的由MR扫描仪110应用的一个或多个参数(例如，MR脉冲序列的类型、TR、重复计数、TI等)、与梯度线圈230产生的梯度场或射频场(例如RF中心频率、翻转角)相关的一个或多个参数、与本申请其他地方所述的RF线圈240检测到的回波信号(例如TE、自旋回波类型)有关的一个或多个参数等，或其任何组合。在一些实施例中，与MR扫描仪110相关的参数可以包括一个或多个其他成像参数，例如RF通道的计数(或数量)、图像对比度和/或比率、切片厚度、成像类型(例如，T1加权成像、T2加权成像、质子密度加权成像等)、图像的视野(FOV)、MR扫描仪110的偏心频移等，或其组合。
在一些实施例中，脉冲序列模块250可以连接到处理设备120和/或与处理设备120通信。例如，在MRI扫描处理之前，处理设备120可以根据临床需求或扫描协议来设计和/或确定与MR扫描仪110相关的参数和布置的至少一部分，并将其发送到脉冲序列模块250，MR扫描仪110可以基于由脉冲序列模块250定义的参数和布置来扫描对象210。例如，MR扫描仪110可以应用具有与由脉冲序列模块250定义的MR脉冲序列相关的特定参数的MR脉冲序列，并且RF线圈240可以根据与由脉冲序列模块250定义的回波信号相关的特定参数来接收回波信号。
关于以上所提供关于MR扫描仪110的描述旨在是说明性的，而不是限制本申请的范围。许多替代、修改和变化对于本领域技术人员来说是显而易见的。本文描述的示例性实施例的特征、结构、方法和特性可以以各种方式组合以获取附加和/或替代示例性实施例。例如，可以将脉冲序列模块250集成到处理设备120中。然而，这些变化和修改不脱离本申请的范围。
图3是根据本申请的一些实施例所示的计算设备300的示例性硬件和/或软件组件的示意图。在一些实施例中，MRI系统100的一个或多个组件可以实现在计算设备300的一个或多个组件上。仅作为示例，处理设备120和/或终端140可以分别实现在计算设备300的一个或多个组件上。
如图3所示，计算设备300可以包括处理器310、存储器320、输入/输出(I/O)330和通信端口340。处理器310可以根据本文描述的技术执行计算机指令(例如，程序代码)和执行处理设备120的功能。计算机指令可以包括例如例程、程序、对象、组件、数据结构、过程、模块和功能，其可以执行本申请所述的特定功能。例如，处理器310可以处理从MR扫描仪110、存储设备130、终端140和/或MRI系统100的任何其他组件获取的对象的图像数据。在一些实施例中，处理器310可以基于MR扫描仪110获取的与对象相关的MRI数据来确定对象的T1值。
在一些实施例中，处理器310可以包括一个或多个硬件处理器，例如微控制器、微处理器、精简指令集计算机(RISC)、专用集成电路(ASICs)、专用指令集处理器(ASIP)、中央处理单元(CPU)、图形处理单元(GPU)，物理处理单元(PPU)、微控制器单元、数字信号处理器(DSP)、现场可编程门阵列(FPGA)、高级RISC机器(ARM)、可编程逻辑器件(PLD)、能够执行一个或多个功能的任何电路或处理器等，或其组合。
仅用于说明，在计算设备300中仅描述了一个处理器。然而，应当注意，本申请中的计算设备300还可以包括多个处理器。因此，如本申请所述的由一个处理器执行的操作和/或方法步骤也可以由多个处理器联合或单独执行。例如，如果在本申请中，计算设备300的处理器同时执行步骤A和步骤B，则应当理解，步骤A和步骤B也可以由计算设备300中的两个或多个不同的处理器共同或单独地执行(例如，第一处理器执行步骤A，第二处理器执行步骤B，或者第一和第二处理器联合执行步骤A和B)。
存储器320可以存储从MR扫描仪110、存储设备130、终端140和/或MRI系统100的任何其他组件获取的数据/信息。在一些实施例中，存储器320可以包括大容量存储设备、可移动存储设备、易失性读写存储器、只读存储器(ROM)，或类似的，或其组合。例如，大容量存储设备可以包括磁盘、光盘、固态驱动器等。可移动存储设备可以包括闪存驱动器、软盘、光盘、存储卡、压缩盘、磁带等。易失性读写存储器可以包括随机存取存储器(RAM)。RAM可以包括动态RAM(DRAM)、双数据速率同步动态RAM(DDR-SDRAM)、静态RAM(SRAM)、晶闸管RAM(T-RAM)和零电容RAM(Z-RAM)等。ROM可以包括掩模ROM(MROM)、可编程ROM(PROM)、可擦除可编程ROM(EPROM)、电可擦除可编程ROM(EEPROM)、光盘ROM(CD-ROM)和数字多功能磁盘ROM等。在一些实施例中，存储器320可以存储一个或多个程序和/或指令，以执行本申请中描述的示例性方法。例如，存储器320可以存储用于处理设备120的用于确定对象的T1值的程序。
输入/输出330可以输入和/或输出信号、数据、信息等。在一些实施例中，输入/输出330可以允许与计算设备300(例如，处理设备120)的用户交互。在一些实施例中，输入/输出330可以包括输入设备和输出设备。输入设备的示例可包括键盘、鼠标、触摸屏、麦克风等，或其组合。输出设备的示例可以包括显示设备、扬声器、打印机、投影仪等或其组合。显示装置的实例可包括液晶显示器(LCD)、基于发光二极管(LED)的显示器、平板显示器、曲面屏幕、电视装置、阴极射线管(CRT)、触摸屏等，或其组合。
通信端口340可以连接到网络(例如，网络150)以促进数据通信。通信端口340可以在计算设备300(例如，处理设备120)和MRI系统100的一个或多个组件(例如，MR扫描仪110、存储设备130和/或终端140)之间建立连接。该连接可以是有线连接、无线连接、能够实现数据传输和/或接收的任何其他通信连接、和/或这些连接的组合。有线连接可以包括例如电缆、光缆、电话线等或其组合。无线连接可以包括例如蓝牙TM连接、Wi-FiTM链接、WiMaxTM链路、WLAN链路、紫蜂链路、移动网络链路(例如3G、4G、5G等)等，或其组合。在一些实施例中，通信端口340可以是和/或包括标准化的通信端口，例如RS232、RS485等。在一些实施例中，通信端口340可以是专门设计的通信端口。例如，可以根据数字成像和医学通信(DICOM)协议来设计通信端口340。
图4是根据本申请一些实施例所示的可以实现的移动设备400的示例性硬件和/或软件组件的示意图。在一些实施例中，MRI系统100的一个或多个组件可以在移动设备400的一个或多个组件上实现。仅作为示例，终端140可以在移动设备400的一个或多个组件上实现。
如图4所示，移动设备400可以包括通信平台410、显示器420、图形处理单元(GPU)430、中央处理单元(CPU)440、输入/输出450、内存460和存储器490。在一些实施例中，任何其他合适的组件，包括但不限于系统总线或控制器(未示出)，也可以包括在移动设备400中。在一些实施例中，操作系统470(例如iOSTM、安卓TM、Windows PhoneTM)和一个或多个应用程序480可以从存储器490加载到内存460中，以便由CPU 440执行。应用程序480可以包括浏览器或任何其他合适的移动应用，用于接收和呈现与MRI系统100有关的信息。用户与信息流的交互可以经由输入/输出450实现，并且经由网络150提供给处理设备120和/或MRI系统100的其他组件。
为了实现本申请中描述的各种模块、单元及其功能，可以将计算机硬件平台用作本文所描述的一个或多个元件的硬件平台。具有用户接口元件的计算机可用于实现个人计算机(PC)或任何其它类型的工作站或终端设备。如果程控得当，计算机也可以充当服务器。
图5是根据本申请的一些实施例所示的示例性处理设备120的框图。如图5所示，处理设备120可以包括获取模块502、参考系数确定模块504和T1测量模块506。
获取模块502可以被配置为获取与MRI系统100相关的信息，当获取模块502可以获取与物理点(本申请披露的对象的示例性实施例)相关的第一组MRI数据，以及获取与物理点相关的第二组MRI数据，第一组MRI数据是第一采集中对象达到第一T1加权水平时由MR扫描仪获取的，第二组MRI数据是第二采集中当对象达到第二T1加权水平时由MR扫描仪获取的。可通过对患者施加由一个或多个脉冲序列参数(例如，TR、FA)定义的MR脉冲序列来实施第一采集和第二采集。在一些实施例中，可以根据一个或多个第一值实施第一采集，每个第一值对应于一个或多个脉冲序列参数中的一个脉冲序列参数。可以根据一个或多个第二值实施第二采集，每个第二值对应于一个或多个脉冲序列参数中的一个脉冲序列参数。一个或多个脉冲序列参数中至少一个脉冲序列的第一值和第二值可以不同。关于获取第一组和第二组MRI数据的更多描述可以在本申请的其他地方找到。例如，参见步骤601和602及其相关描述。
参考系数确定模块504可以被配置为基于第一组MRI数据和第二组MRI数据确定参考系数的目标值，该参考系数与第一采集中的第一B1不均匀性和第二采集中的第二B1不均匀性相关联。B1采集过程中可能出现不均匀性，并导致采集过程中物理点处的实际FA和预设FA之间存在差异。由于第一B1不均匀性，在第一采集期间物理点处的第一实际FA可能不同于第一采集的第一FA。由于第二B1不均匀性，在第二采集期间物理点处的第二实际FA可能不同于第二采集的第二FA。本申请可以利用与第一B1不均匀性和第二B1不均匀性均相关联的参考系数来进行T1测量，从而改进T1测量的精度。关于参考系数的目标值的确定的更多描述可以在本申请的其他地方找到。例如，参见步骤603及其相关描述。
T1测量模块506可被配置为至少部分地基于参考系数的目标值来确定物理点的T1值。在一些实施例中，T1测量模块506可以至少部分地基于参考系数的目标值来确定关于T1测量的校正系数，并且至少部分基于校正系数来确定物理点的T1值。关于物理点的T1值的确定的更多描述可以在本申请的其他地方找到。例如，参见步骤604及其相关描述。
在一些实施例中，这些模块可以是处理设备120的全部或部分的硬件电路。这些模块还可以实现为由处理设备120读取和执行的应用程序或指令集。此外，这些模块可以是硬件电路和应用程序/指令的任意组合。例如，当处理设备120正在执行应用/指令集时，这些模块可以是处理设备120的一部分。
应当注意，关于处理设备120的上述描述是为了说明的目的而提供的，并不旨在限制本申请的范围。对于本领域的普通技术人员，可以根据本申请的教导进行多种变化和修改。然而，这些变化和修改并不脱离本申请的范围。在一些实施例中，处理设备120还可以包括一个或多个附加模块，例如存储模块。另外或可替换地，可以省略上述模块中的一个或多个。此外，上述任何模块可以在两个或多个独立单元中实现。
图6是根据本申请的一些实施例所示的用于确定对象的T1值的示例性过程的流程图。在一些实施例中，过程600的一个或多个步骤可以在图1所示的MRI系统100中实现。例如，过程600可以以指令的形式存储在MRI系统100的存储设备(例如，存储设备130、存储器320和/或存储器490)中，并由处理设备120(例如，如图3所示的计算设备300的处理器310、如图4所示的移动设备400的CPU 440、如图5所示的一个或多个模块)调用和/或执行。
如本文所使用的，对象的T1值是指纵向磁化(例如，沿着主磁场)再生的时间常数。出于说明目的，以下描述将患者的物理点作为示例性对象。应当注意，对象可以是另一受试体的物理点，例如患者的一部分、动物、人造受试体(例如，模型)等。此外，应当理解，下文中物理点的T1值的确定仅作为示例提供，并且过程600可用于确定对象的一个或多个其他定量参数(例如，T2值、T2*值、R2值)的值。
如本申请其他部分所述，为了在物理点上进行T1测量，当物理点达到不同的T1加权水平时，可以应用具有不同脉冲序列参数值(例如，不同的翻转角或不同的重复时间)的多次采集来获取多组MRI数据。由于在每次采集中发生的B1不均匀性，在在采集期间，物理点处的实际FA和预设FA之间可能存在差异，这可能影响在采集中采集的一组MRI数据和基于该组MRI数据进行的T1测量的精度。此外，不同的采集可能具有不同的B1不均匀性。因此，可能需要通过考虑不同采集中的不同B1不均匀性来确定物理点的T1值，从而提高所确定的T1值的精度。
在601中，处理设备120(例如，获取模块502)可以获取与物理点相关的第一组MRI数据，该第一组MRI数据是第一采集中物理点达到第一T1加权水平时由MR扫描仪获取的。
在602中，处理设备120(例如，获取模块502)可以获取与物理点相关的第二组MRI数据，该第二组MRI数据时第二采集中物理点达到第二T1加权水平时由MR扫描仪获取的。
在采集(例如，第一或第二采集)中获取的一组MRI数据可以包括MR扫描仪通过扫描患者获取的数据(例如，一组回波信号)。例如，第一组MRI数据可以包括在第一采集中获取的第一组回波信号，第二组MRI数据可以包括在第二采集中获取的第二组回波信号。在一些实施例中，MR扫描仪可包括多个线圈单元，每个线圈单元可被配置为在患者扫描期间独立地检测回波信号。第一组和第二组回波信号可以包括分别在第一采集和第二采集期间由一个或多个线圈单元检测到的回波信号。
可通过在患者上施加MR脉冲序来列对患者进行第一采集和第二采集。第一采集和第二采集可以在患者扫描期间以任何顺序执行。MR脉冲序列可以由一个或多个脉冲序列参数定义，例如，MR脉冲序列的类型、采集的计数(或数目)、采集的TR、MR脉冲序列中激励脉冲的FA、反转时间等，或其任何组合。用于实现第一和第二采集的示例性MR脉冲序列可以包括可变FA(VFA)脉冲序列(例如，在不同采集中具有不同FAs的GRE脉冲序列)、可变TR(VTR)脉冲序列(例如，在不同采集中具有不同TR的GRE脉冲序列)等。
在一些实施例中，可以根据一个或多个第一值实施第一采集，每个第一值对应于一个或多个脉冲序列参数中的一个脉冲序列参数。可以根据一个或多个第二值实施第二采集，每个第二值对应于一个或多个脉冲序列参数中的一个脉冲序列参数。一个或多个脉冲序列参数中的至少一个脉冲序列参数的第一值和第二值可以不同。在一些实施例中，在采集中获取的一组MRI数据可能受到与采集相对应的脉冲序列参数(例如，TR、FA)的值、物理点的T1值和/或物理点的一个或多个其他定量参数的影响。在第一和第二采集中使用至少一个脉冲序列参数的不同值可以在第一和第二采集中实现不同的T1加权水平。采集中的T1加权水平可以反映T1值相对于在采集中采集的一组MRI数据的影响。例如，采集中的T1加权水平可以反映T1值相对于采集中获取的MRI信号的强度的影响。
出于说明的目的，以下描述将GRE脉冲序列作为示例性的脉冲序列，但并不意在限制。在由GRE脉冲序列实现的采集中，在采集期间的物理点处的信号强度可根据如下等式(1)(或称为Block等式)来表示：
Figure BDA0003204049210000181
其中S指采集期间物理点处的信号强度，α是指采集中激励脉冲的预设FA(翻转角)。E1可以等于e-TR/T1，其中TR指采集中两个重复且连续的RF脉冲之间的时间间隔，T1是指物理点的T1值。在一些实施例中，由于TR远小于T1值，因此E1可近似为
Figure BDA0003204049210000182
根据等式(1)，物理点处的信号强度可与采集中GRE脉冲序列的FA和TR以及物理点的T1值相关联。通过在第一采集和第二采集中使用不同的FAs和/或TRs，第一采集中物理点的T1加权水平(即，第一T1加权水平)可能不同于第二采集中物理点的T1加权水平(即，第二T1加权水平)。仅作为示例，可以根据第一FA(表示为α1)实施第一采集，可以根据不同于第一FA(α1)的第二FA(表示为α2)实施第二采集，并且除FA之外的脉冲序列参数(例如TR)的第一值和第二值相同。关于获取第一组和第二组MRI数据的更多描述可以在本申请的其他地方找到。参见示例1和示例2。
在603中，处理设备120(例如，参考系数确定模块504)可以基于第一组MRI数据和第二组MRI数据确定参考系数的目标值，该参考系数与第一采集中的第一B1不均匀性和第二采集中的第二B1不均匀性相关联。
如本申请的其他地方所述，B1不均匀性可在采集期间发生，并导致在采集期间的物理点处的实际FA和预设FA之间的差异。例如，当采集的预设FA相对较小(例如，小于阈值角)时，采集期间物理点处的实际FA可近似等于γB1t，其中t指在采集期间施加激励脉冲的时间段，B1是采集过程中的实际B1场，γ是旋磁比。由于第一B1不均匀性，在第一采集期间的物理点处的第一实际FA(表示为α′1)可偏离第一采集的第一FAα1。例如，α′1可等于r1*α1，其中r1可反映第一B1不均匀性的影响因子。类似地，由于第二B1不均匀性，在第二采集期间的物理点处的第二实际FA(表示为α′2)可偏离第二采集的第二FAα2。例如，α′ 2可等于r2*α2，其中r2可反映第二B1不均匀性的影响因子。
传统的T1测量方法可以假设第一B1不均匀性与第二B1不均匀性一致，并确定第一实际FAα′1(或r1)和/或第二个实际FAα′2(或r2)。然而，由于第一B1不均匀性可能不同于第二B1不均匀性(即，r1可能不同于r2)，因此传统的T1测量方法可能具有降低的精度。本申请可以利用与第一B1不均匀性和第二B1不均匀性两者相关联的参考系数来进行T1测量，从而改进T1测量的精度。
在一些实施例中，参考系数可以测量(描述或者代表)第一采集期间物理点处的第一信号强度与第二采集期间物理点处的第二信号强度之间的差别(difference)。仅作为示例，参考系数可以是第二信号强度与第一信号强度的比率，其可以根据如下等式(2)来表示：
Figure BDA0003204049210000191
其中S(α′1)指第一采集期间物理点处的第一信号强度，S(α′2)指第二采集期间物理点处的第二信号强度，R′指参考系数的目标值。
如等式(2)所示，参考系数可以与α′1(例如r1*α1)和α′2(例如r2*α2)均相关联。换句话说，参考系数可与第一采集中的第一B1不均匀性(例如，由r1表示)和第二采集中的第二B1不均匀性(例如，由r2表示)两者相关联。应当理解，提供等式(2)仅仅是为了说明的目的，参考系数可以用任何其他形式表示。例如，可以通过从第二信号强度减去第一信号强度来确定参考系数。又例如，参考系数可以是第一信号强度与第二信号强度的比率。
参考系数的目标值是指基于第一和第二组MRI数据确定的参考系数的实际值。例如，可以在患者的MRI扫描中获取第一组MRI数据和第二组MRI数据。处理设备120可以使用诸如傅立叶变换(FFT)算法的MR图像重建算法基于第一组MRI数据重建患者的第一图像。处理设备120还可以使用图像重建算法基于第二组MRI数据重建患者的第二图像。处理设备120还可以基于第一图像和第二图像确定参考系数的目标值。
例如，与第一图像中的物理点对应的像素值或体素值可以表示在第一采集期间物理点处的第一信号强度，并且与第二图像中的物理点对应的像素值或体素值可以表示在第二采集期间物理点处的第二信号强度。处理设备120可以通过例如将第二图像除以第一图像来生成第一图像和第二图像之间的差值图像。处理设备120可以从差值图像识别与物理点相对应的元素(例如，像素或体素)，并将该元素的值指定为参考系数的目标值。又例如，处理设备120可以在第一图像中识别与物理点相对应的第一元素(例如，像素或体素)，并且在第二图像中识别与物理点相对应的第二元素(例如，像素或体素)。处理设备120可以基于第一元素的值和第二元素的值来确定参考系数的目标值。仅作为示例，可以将第二元素的值与第一元素的值的比率确定为参考系数的目标值。
在604中，处理设备120(例如，T1测量模块506)可以至少部分地基于参考系数的目标值来确定物理点的T1值。
在一些实施例中，处理设备120可以至少部分地基于参考系数的目标值来确定关于T1测量的校正系数，并且至少部分地基于校正系数来确定物理点的T1值。物理点的T1值的确定可涉及基于校正系数而进行的B1场校正，使得所确定的T1值可被视为物理点的校正T1值。关于物理点的T1值的确定的更多描述可以在本申请的其他地方找到。例如，参见图7及其相关描述。在一些实施例中，可以省略步骤604。在一些实施例中，参考系数的目标值可用于对物理点进行其它定量测量。
应当注意，关于过程600的上述描述仅仅是为了说明的目的而提供的，并不意在限制本申请的范围。对于本领域的普通技术人员，可以根据本申请的教导进行多种变化和修改。然而，这些变化和修改并不脱离本申请的范围。在一些实施例中，过程600可以用一个或多个未描述的附加步骤和/或不用上面讨论的一个或多个步骤来完成。例如，过程600可以包括将物理点的T1值发送到终端设备(例如，医生的终端140)以进行显示的附加操作。
在一些实施例中，可以对患者进行三次或多次采集。处理设备120可以获取与物理点相关的两组以上的MRI数据。至少两组MRI数据中的每组MRI数据可以在其中一次采集中获取，该采集中物理点达到区别性的T1加权水平。例如，上述采集可以应用不同的FAs。在一些实施例中，上述三次或多次采集可以形成至少两对采集。例如，具有最小FA的采集可以与其余采集的每次采集配对。又例如，上述三次或多次采集可以随机配对。对于每一次采集，处理设备120可以基于与该对采集相对应的两组MRI数据来确定物理点的初始T1值。例如，对于每一对采集，该对采集的一次采集可被视为第一采集，另一次采集可被视为第二采集，并且处理设备120可对该对采集执行过程600以确定初始T1值。处理设备120可以通过组合与该对采集对相对应的初始T1值来进一步确定物理点的T1值。例如，初始T1值的平均值可以被确定为物理点的T1值。
在一些实施例中，患者可以包括多个物理点。对于患者的每个物理点，处理设备120可以通过执行过程600来确定物理点的T1值。可选地，可以为医疗诊断生成T1图，该T1图包括每个物理点的T1值。
图7是根据本申请的一些实施例所示的用于至少部分地基于参考系数的目标值来确定物理点的T1值的示例性过程的流程图。在一些实施例中，过程700的一个或多个操作可以在图1所示的MRI系统100中实现。例如，过程700可以以指令的形式存储在MRI系统100的存储设备(例如，存储设备130、存储器320和/或存储器490)中，并由处理设备120(例如，如图3所示的计算设备300的处理器310、如图4所示的移动设备400的CPU 440、如图5所示的一个或多个模块)调用和/或执行。在一些实施例中，可以执行过程700的一个或多个操作以实现如图6所描述的操作604的至少一部分。
如图6所述，可以对患者执行第一采集和第二采集，以获取患者的第一组MRI数据和第二组MRI数据。第一采集和第二采集可对应于至少一个脉冲序列参数(例如FA)的不同值。基于第一组MRI数据和第二组MRI数据，可以确定参考系数的目标值，该参考系数与第一采集中的第一B1不均匀性和第二采集中的第二B1不均匀性相关联。出于说明目的，假设通过应用GRE脉冲序列来获取第一组和第二组MRI数据。GRE脉冲序列可以包括根据第一翻转角α1执行的第一采集以及根据第二翻转角α2执行的第二采集。第一和第二采集可以具有相同的TR。在下文中，将第二采集期间物理点处的第二信号强度与第一采集期间物理点处的第一强度的比率作为示例性参考系数。应当注意，提供以下描述仅仅是为了说明的目的，并不打算限制本申请。参考系数可以是任何其他合适的参数。可通过应用任何其它MR脉冲序列来获取第一组和第二组MRI数据。
在701中，处理设备120(例如，T1测量模块506)可以基于一个或多个脉冲序列参数的一个或多个第一值和一个或多个第二值来确定参考系数的参考值。
参考系数的参考值是指当不考虑第一B1不均匀性和第二B1不均匀性时参考系数的估计值(或理想值)。例如，如果不考虑第一B1不均匀性和第二B1不均匀性，则第一采集中的第一实际FAα1 ′可以等于第一FAα1，第二采集中的第二实际FAα′2可以等于第二FAα2。根据α1、α2和其它脉冲序列参数的值，可将第二采集中物理点处的第二估计(或理想)信号强度与第一采集中物理点处的第一估计(或理想)信号强度的比率确定为参考系数的参考值。仅作为示例，参考系数的参考值可以表示为
Figure BDA0003204049210000221
其中S(α1)可以表示当α′1等于α1时在第一采集中的物理点处的第一估计(或理想)信号强度，并且S(α2)可表示当α′2等于α2时在第二采集中的物理点处的第二估计(或理想)信号强度。
在一些实施例中，处理设备120可以确定物理点的参考T1值，并且基于参考T1值和如图6所述的等式(1)来确定参考系数的参考值。例如，物理点的参考T1值可根据物理点所在的器官或组织、物理点的材料类型(例如，水、脂肪、肌肉)等或其任何组合来确定。仅作为示例，如果物理点位于患者的大脑中，则参考T1值可以是人脑的平均T1值。可通过将参考T1值、α1和TR输入到等式(1)来确定S(α1)。可通过将参考T1值、α2和TR输入到等式(1)中来确定S(α2)。然后，处理设备120可以通过将S(α2)除以S(α1)来确定参考系数(即R)的参考值。
在702中，处理设备120(例如，T1测量模块506)可以基于参考值和参考系数的目标值确定关于T1测量的校正系数。
在一些实施例中，校正系数可以测量参考系数的目标值(即R’)和参考系数的参考值(即R)之间的相似度。例如，校正系数可以是R与R’的比率或R’与R的比率。与参考系数类似，校正系数(例如，R/R’或R’/R)也可以同时考虑第一B1不均匀性和第二B1不均匀性。
根据如结合图6所描述的等式(1)，在采集期间在物理点处的信号强度可以受到物理点处的实际FA(其反映了在采集期间B1的不均匀性)和物理点的T1值的影响。在这种情况下，校正系数可以与第一B1不均匀性、第二B1不均匀性以及物理点的T1值相关联。根据如图8所述的实验结果，与第一和第二B1不均匀性相比，T1值对校正系数的影响相对较小。当T1值超过阈值时，T1值对校正系数的影响可能非常小，并且可以忽略T1值的影响。换句话说，校正系数可以主要由第一B1不均匀性和第二B1不均匀性确定。因此，可以使用校正系数来消除或减小T1测量中的第一B1不均匀性和第二B1不均匀性的影响。
在703中，处理设备120(例如，T1测量模块506)可以至少部分地基于修正系数来确定物理点的T1值(表示为T″1)。
在一些实施例中，处理设备120可以基于第一翻转角α1(或用第一翻转角α1获取的第一组MRI数据)、第二翻转角α2(或用第二翻转角α2获取的第二组MRI数据)，以及参考系数的参考值R来确定物理点的初始T1值(表示为T′1)。例如，物理点的初始T1值可以通过将R、TR、α1和α2输入到等式(3)中来确定，如下所示：
Figure BDA0003204049210000222
处理设备120可以通过根据校正系数校正初始T1值来进一步确定物理点的T1值。例如，物理点的T1值可以等于T′1*r，其中r可以表示校正系数。
在一些实施例中，处理设备120可以根据校正系数校正第一翻转角α1和第二翻转角α2中的一个。例如，如果参考系数等于
Figure BDA0003204049210000231
则可以校正第一翻转角α1来确定校正后的第一翻转角α″1。校正后的第一翻转角α″1可以等于α1*r。然后，处理设备120可以基于校正后的第一翻转角α″1、第二翻转角α2和参考系数的目标值R'来确定物理点的T1值。例如，物理点的T1值(T″1)可以通过将α″1、α2、和R’输入到等式(3)中来确定，即，根据如下等式(4)：
Figure BDA0003204049210000232
如果参考系数等于
Figure BDA0003204049210000233
则可以校正第二翻转角α2来确定校正后的第二翻转角α″2。校正后的第二个翻转角α″2可以等于α2*r。然后，处理设备120可以基于第一翻转角α1、校正后的第二翻转角α″2和参考系数的目标值R’来确定物理点的T1值。例如，物理点的T1值(T″1)可以通过将α1、α″2、和R’输入到等式(3)中来确定，即，根据等式(5)如下：
Figure BDA0003204049210000234
根据本申请的一些实施例，可以确定参考系数的目标值和参考值。参考系数可与第一采集中的第一B1不均匀性和第二采集中的第二B1不均匀性均相关联。基于目标值和参考系数的参考值，可以确定校正系数，而非第一B1不均匀性或第二B1不均匀性的。校正系数可用于消除或减少在确定物理点的T1值时第一和第二B1不均匀性的影响，从而提高所确定的T1值的精度。
应当注意，关于过程700的上述描述仅仅是为了说明的目的而提供的，并不意在限制本申请的范围。对于本领域的普通技术人员，可以根据本申请的教导进行多种变化和修改。然而，这些变化和修改并不脱离本申请的范围。在一些实施例中，过程700可以用一个或多个未描述的附加操作和/或不用上面讨论的一个或多个操作来完成。此外，可以根据实际需要修改上述等式。例如，可以省略等式的一个或多个系数，和/或等式可以包括一个或多个附加系数。又例如，根据如上所述的等式(1)到(3)实现的T1测量的解析解可以被用于T1测量的一个或多个其他算法(例如，驱动平衡单脉冲观测(DESPOT)算法)代替。又例如，如果第一采集和第二采集是根据不同于GRE脉冲序列的脉冲序列进行的，则采集过程中物理点的信号强度可以用不同于等式(1)的另一个等式来表示。
示例
提供这些示例是为了说明的目的，而不是为了限制本申请的范围。
例1
进行模拟实验，以评估图6和7所述的T1测量系统和方法的准确性。在该模拟实验中，以患者大脑的一个物理点作为示例性对象。物理点的参考T1值(表示为T1ref)的范围设定在300到3500毫秒(ms)之间，这是典型3TMRI扫描中脑组织的典型范围。假设通过应用GRE脉冲序列在患者的大脑上进行MRI扫描，所述GRE脉冲序列包括具有第一FA的第一采集和具有第二FA的第二采集。第一FA和第二FA分别设置为等于4°和32°。第一采集和第二采集的TRs均等于10ms。
对于物理点的特定参考T1值，物理点的理想T1值(表示为T1)、初步未校正的T1值(表示为t′1)、以及校正后的T1值(表示为T″1)根据以下描述确定。
参考系数的参考值R基于TR(即10ms)、第一FA(即4°)、第二FA(即32°)和特定参考值来确定。具体而言，参考值R等于S(α2)到S(α1)，其中S(α1)表示在第一采集期间物理点处的第一估计信号强度，并且S(α2)表示在第二采集期间物理点处的第二估计信号强度。根据操作701所描述的方法来确定参考值R。然后，如703所述通过将参考值R、第一FA(即4°)和第二FA(即32°)来输入到等式3中来确定与特定参考T1值相对应的理想T1值。理想T1值的确定是在没有考虑第一采集中的第一B1不均匀性和第二采集中的第二个B1不均匀性的情况下进行的。
根据蒙特卡洛模拟算法进行多次模拟扫描。在每次模拟扫描中，添加随机噪声，使关于第一B1不均匀性的第一随机影响因子乘以第一实际FA，以生成第一采集的第一实际FA，以及使第二B1不均匀性的第二随机影响因子乘以第二FA，以产生第二采集的第二实际FA。第一随机影响因子和第二随机影响因子都在0.8到1.2之间。对于每次模拟扫描，参考系数的目标值R’被确定为
Figure BDA0003204049210000241
其中S(α′1)表示在第一采集期间物理点处的实际第一信号强度，并且S(α′2)表示在第二采集期间物理点处的实际第二信号强度。S(α′1)以及S(α′2)根据模拟扫描获取的模拟MRI数据确定。然后，如703所述，通过将R′、第一FA(即，4°)和第二FA(即，32°)输入到等式(3)中来确定与特定参考T1值相对应的物理点的初始未校正T1值。初始未校正T1值的确定考虑了第一和第二B1不均匀性(即，通过使用参考系数的目标值)，但是在没有B1场校正的情况下执行。
此外，对于每次模拟扫描，R与R’的比率被确定为校正系数(表示为r)。通过将校正系数乘以第二FA来确定校正后的第二FA。如703所述的，通过将R’、第一FA(即，4°)和校正后的第二FA输入等式(3)中来确定与特定参考T1值相对应的物理点的校正后的T1值。校正后的T1值的确定考虑了第一和第二B1的不均匀性(即利用目标值(参考系数)，并涉及B1场校正。
图8是示出的基于模拟实验的结果确定的物理点对应的校正系数和参考T1值之间的关系示意图。在图8中，点的水平坐标表示特定参考T1值，点的垂直坐标表示基于与特定参考T1值对应的多次模拟扫描而确定的校正系数的平均值，并且通过该点的虚线表示校正系数的方差值。如图8所示，与不同的参考T1值相对应的校正系数大致相同，这表明物理点的T1值对物理点的校正系数几乎没有影响。
图9A是示出的基于模拟实验的结果确定的理想T1值和参考T1值之间的关系的示意图。在图9A中，点的水平坐标表示特定参考T1值，点的垂直坐标表示基于与特定参考T1值对应的多次模拟扫描确定的理想T1值与特定参考T1值的第一比率的平均值，并且通过该点的虚线表示第一比率的方差值。如图9A所示，当参考T1值小(例如，小于200ms)时，理想T1值大于参考T1值，这表明误差是由在T1测量中进行的数学运算(例如，由等式(1)和(3))引起的。当参考T1值大于200ms时，理想T1值接近参考T1值，数学运算引起的误差可以忽略不计。
图9B是示出的基于模拟实验的结果确定的初步未校正T1值和参考T1值之间的关系的示意图。在图9B中，点的水平坐标表示特定参考T1值，点的垂直坐标表示基于与特定参考T1值对应的多次模拟扫描确定的初始未校正T1值与参考T1值的第二比率的平均值，并且通过该点的虚线表示第二比率的方差值。如图9B所示，初始未校正T1值与参考T1值有较大的偏差，特别是当参考T1值小于200ms或大于2000ms时。
图9C是示出的基于模拟实验的结果确定的校正T1值和参考T1值之间的关系的示意图。在图9C中，点的水平坐标表示特定参考T1值，点的垂直坐标表示校正后的T1值与基于与特定参考T1值对应的多次模拟扫描确定的特定参考T1值的第三比率的平均值，并且通过该点的虚线表示第三比率的方差值。如图9C所示，当参考T1值大于200ms时，校正后的T1值基本上等于参考T1值并且具有相对较小的方差。
根据图9A、9B和9C，本文公开的系统和方法可以提高T1测量的精度。
例2
根据3D-GRE脉冲序列对ISMRM水模进行MRI扫描。在ISMRM水模的扫描过程中，第一采集的第一FA为4°，第二采集的第二FA为16°。ISMRM水模的每个物理点都有一个已知的参考T1值。对于ISMRM水模的每个物理点，根据第一FA和第二FA确定初始T1值。另外，对于ISMRM水模的每个物理点，通过执行本文公开的T1测量方法(例如，图6和7)来确定物理点的校正系数的值和校正后的T1值。
图10A示出了显示ISMRM水模的每个物理点的初始T1值的初始T1图。图10B示出了ISMRM水模的每个物理点的校正系数的校正系数图。图10C示出了显示ISMRM水模的每个物理点的校正T1值的校正后的T1图。与初始T1图相比，由本文公开的系统和方法生成的校正后的T1图具有更高的质量(例如，更均匀)。
图10D是示出的基于ISMRM水体模的MRI扫描生成的T1测量结果的图。横轴表示参考T1值，纵轴表示估计T1。空心点的垂直坐标表示对应于特定参考T1值的初始T1值。实心点的垂直坐标表示对应于特定参考T1值的校正后的T1值。虚线中各点的纵坐标等于该点的横坐标，即虚线的斜率等于1，如图10D所示，对于一个特定的参考T1值，相应的校正后的T1值比相应的初始T1值更接近于特定的参考T1值。这表明本申请的系统和方法可确定具有提高的精度和可靠性的物理点的T1值。
对于本领域的普通技术人员，可以根据本申请的教导进行多种变化和修改。然而，这些变化和修改并不脱离本申请的范围。
上文已对基本概念做了描述，显然，对于阅读此申请后的本领域的普通技术人员来说，上述发明披露仅作为示例，并不构成对本申请的限制。虽然此处并未明确说明，但本领域的普通技术人员可能会对本申请进行各种修改、改进和修正。该类修改、改进和修正在本申请中被建议，所以该类修改、改进、修正仍属于本申请示范实施例的精神和范围。
上文已对基本概念做了描述，显然，对于阅读此申请后的本领域的普通技术人员来说，上述发明披露仅作为示例，并不构成对本申请的限制。虽然此处并未明确说明，但本领域的普通技术人员可能会对本申请进行各种修改、改进和修正。这些改变，改进和修改旨在通过本申请建议，并且在本申请的示例性实施例的精神和范围内。
同时，本申请使用了特定词语来描述本申请的实施例。例如“一个实施例”、“一实施例”、和/或“一些实施例”意指与本申请至少一个实施例相关的某一特征、结构或特性。因此，应当强调并注意的是，本说明书中在不同位置两次或多次提及的“一实施例”或“一个实施例”或“一替代性实施例”并不一定是指同一实施例。此外，本申请的一个或多个实施例中的某些特征、结构或特点可以进行适当的组合。
此外，本领域的普通技术人员可以理解，本申请的各方面可以通过若干具有可专利性的种类或情况进行说明和描述，包括任何新的和有用的过程、机器、产品或物质的组合，或对其任何新的和有用的改进。相应地，本申请的各个方面可以完全由硬件执行、可以完全由软件(包括韧体、常驻软件、微代码等)执行、也可以由硬件和软件组合执行。以上硬件或软件均可被称为“单元”、“模块”或“系统”。此外，本申请的各方面可以采取体现在一个或多个计算机可读介质中的计算机程序产品的形式，其中计算机可读程序代码包含在其中。
计算机可读信号介质可能包含一个内含有计算机程序代码的传播数据信号，例如在基带上或作为载波的一部分。此类传播信号可以有多种形式，包括电磁形式、光形式等或任何合适的组合。计算机可读信号介质可以是除计算机可读存储介质之外的任何计算机可读介质，该介质可以通过连接至一个指令执行系统、装置或设备以实现通信、传播或传输供使用的程序。位于计算机可读信号介质上的程序代码可以通过任何合适的介质进行传播，包括无线电、电缆、光纤电缆、RF等，或任何上述介质的组合。
本申请各部分操作所需的计算机程序代码可以用任意一种或以上程序设计语言编写，包括面向对象程序设计语言如Java、Scala、Smalltalk、Eiffel、JADE、Emerald、C++、C#、VB.NET、Python等，常规程序化程序设计语言如C程序设计语言、Visual Basic、Fortran2103、Perl、COBOL 2102、PHP、ABAP，动态程序设计语言如Python、Ruby，和Groovy，或其他程序设计语言等。该程序代码可以完全在用户计算机上运行、或作为独立的软件包在用户计算机上运行、或部分在用户计算机上运行部分在远程计算机运行、或完全在远程计算机或服务器上运行。在后种情况下，远程计算机可以通过任何网络形式与用户计算机连接，比如局域网(LAN)或广域网(WAN)，或连接至外部计算机(例如通过因特网)，或在云计算环境中，或作为服务使用如软件即服务(SaaS)。
此外，除非权利要求中明确说明，本申请所述处理元素和序列的顺序、数字字母的使用、或其他名称的使用，并非用于限定本申请流程和方法的顺序。尽管上述披露中通过各种示例讨论了一些目前认为有用的发明实施例，但应当理解的是，该类细节仅起到说明的目的，附加的权利要求并不仅限于披露的实施例，相反，权利要求旨在覆盖所有符合本申请实施例实质和范围的修正和等价组合。例如，尽管上述各种组件的实现可以体现在硬件设备中，但是它也可以实现为纯软件解决方案，例如，在现有服务器或移动设备上的安装。
同理，应当注意的是，为了简化本申请披露的表述，从而帮助对一个或多个发明实施例的理解，前文对本申请的实施例的描述中，有时会将多种特征归并至一个实施例、附图或对其的描述中。然而，本申请的该方法不应被解释为反映所申明的客体需要比每个权利要求中明确记载的更多特征的意图。相反，发明的主体应具备比上述单一实施例更少的特征。
一些实施例中使用了描述成分、属性数量的数字，应当理解的是，此类用于实施例描述的数字，在一些示例中使用了修饰词“大约”、“近似”或“大体上”来修饰。例如，除非另有说明，否则“大约”、“近似”或“基本上”可以指示其所描述的值的±1％、±5％、±10％或±20％的变化。相应地，在一些实施例中，说明书和权利要求中使用的数值参数均为近似值，该近似值根据个别实施例所需特点可以发生改变。在一些实施例中，数值参数应考虑规定的有效数位并采用一般位数保留的方法。尽管本申请一些实施例中用于确认其范围广度的数值域和参数为近似值，在具体实施例中，此类数值的设定在可行范围内尽可能精确。在一些实施例中，分类中使用的分类条件是为了说明的目的而提供的，并且根据不同的情况进行修改。例如，“概率值大于阈值”的分类条件可以进一步包括或排除“概率值等于阈值”的条件。

Claims (9)

1.一种磁共振成像的方法，包括：
获取与对象相关的第一组磁共振成像数据，所述第一组磁共振成像数据是第一采集中所述对象达到第一T1加权水平时由磁共振扫描仪获取的，根据一个或多个第一值实施所述第一采集，所述一个或多个第一值中的每个第一值对应于一个或多个脉冲序列参数中的一个脉冲序列参数；
获取与所述对象相关的第二组磁共振成像数据，所述第二组磁共振成像数据是第二采集中当所述对象达到第二T1加权水平时由所述磁共振扫描仪获取的，所述第二T1加权水平和所述第一T1加权水平不同，根据一个或多个第二值实施第二采集，所述一个或多个第二值中的每个第二值对应于所述一个或多个脉冲序列参数中的一个脉冲序列参数，所述一个或多个脉冲序列参数中的至少一个脉冲序列参数的第一值和第二值不同；其中，所述一个或多个脉冲序列参数的所述一个或多个第一值包括第一翻转角，所述一个或多个脉冲序列参数的所述一个或多个第二值包括与所述第一翻转角不同的第二翻转角；以及
基于所述第一组磁共振成像数据和所述第二组磁共振成像数据，确定参考系数的目标值，该参考系数与所述第一采集中的第一B1不均匀性和所述第二采集中的第二B1不均匀性相关联。
2.根据权利要求1所述的方法，其特征在于，所述对象是受试体的物理点，并且所述参考系数测量在所述第一采集中所述物理点处的第一信号强度与在所述第二采集中所述物理点处的第二信号强度之间的差别。
3.根据权利要求2所述的方法，其特征在于，所述参考系数为所述第一信号强度与所述第二信号强度之比或所述第二信号强度与所述第一信号强度之比。
4.根据权利要求1所述的方法，其特征在于，还包括：
基于所述参考系数的目标值来确定所述对象的T1值。
5.根据权利要求4所述的方法，其特征在于，所述基于所述参考系数的目标值来确定所述对象的T1值包括：
基于所述一个或多个脉冲序列参数的所述一个或多个第一值和所述一个或多个第二值确定所述参考系数的参考值；
基于所述参考值和所述参考系数的目标值，确定关于T1测量的校正系数；以及
基于所述校正系数来确定所述对象的T1值。
6.根据权利要求5所述的方法，其特征在于，基于校正系数来确定对象的T1值包括：
基于所述第一翻转角、所述第二翻转角和所述参考系数的参考值来确定所述对象的初始T1值；以及
通过根据所述校正系数校正所述初始T1值来确定所述对象的T1值。
7.根据权利要求5所述的方法，其特征在于，所述基于所述校正系数来确定所述对象的T1值包括：
通过所述校正系数校正所述第一翻转角来确定校正后的第一翻转角；以及
基于所述校正后的第一翻转角、所述第二翻转角和所述参考系数的目标值来确定所述对象的T1值。
8.一种非暂时性计算机可读存储介质，包括至少一组指令，当由至少一个处理器执行时，所述至少一组指令指示所述至少一个处理器执行如权利要求1-7任一所述的磁共振成像方法。
9.一种磁共振成像的系统，包括：
至少一个存储设备，包括一组指令；和
至少一个处理器，被配置为与所述至少一个存储设备通信，其中，当执行所述一组指令时，所述至少一个处理器被配置为指示所述系统执行包括如下内容的操作：
获取与对象相关的第一组磁共振成像数据，所述第一组磁共振成像数据是第一采集中所述对象达到第一T1加权水平时由磁共振扫描仪获取的，根据一个或多个第一值实施所述第一采集，所述一个或多个第一值中的每个第一值对应于一个或多个脉冲序列参数中的一个脉冲序列参数；
获取与所述对象相关的第二组磁共振成像数据，所述第二组磁共振成像数据是第二采集中当所述对象达到第二T1加权水平时由所述磁共振扫描仪获取的，所述第二T1加权水平和所述第一T1加权水平不同，根据一个或多个第二值实施第二采集，所述一个或多个第二值中的每个第二值对应于所述一个或多个脉冲序列参数中的一个脉冲序列参数，所述一个或多个脉冲序列参数中的至少一个脉冲序列参数的第一值和第二值不同；其中，所述一个或多个脉冲序列参数的所述一个或多个第一值包括第一翻转角，所述一个或多个脉冲序列参数的所述一个或多个第二值包括与所述第一翻转角不同的第二翻转角；以及
基于所述第一组磁共振成像数据和所述第二组磁共振成像数据，确定参考系数的目标值，该参考系数与所述第一采集中的第一B1不均匀性和所述第二采集中的第二B1不均匀性相关联。
CN202110912132.1A 2020-08-11 2021-08-10 磁共振成像的方法及系统 Active CN114076911B (zh)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
US16/989,905 2020-08-11
US16/989,905 US11287504B2 (en) 2020-08-11 2020-08-11 Systems and methods for magnetic resonance imaging

Publications (2)

Publication Number Publication Date
CN114076911A CN114076911A (zh) 2022-02-22
CN114076911B true CN114076911B (zh) 2023-07-14

Family

ID=73456566

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110912132.1A Active CN114076911B (zh) 2020-08-11 2021-08-10 磁共振成像的方法及系统

Country Status (2)

Country Link
US (1) US11287504B2 (zh)
CN (1) CN114076911B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US11841413B2 (en) * 2019-12-30 2023-12-12 Case Western Reserve University System and method for B1-selective spatial encoding using magnetic resonance
WO2024020858A1 (zh) * 2022-07-27 2024-02-01 维沃移动通信有限公司 表面构建方法、装置、电子设备及介质

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101919696A (zh) * 2009-06-15 2010-12-22 通用电气公司 用于磁共振rf场测量的系统、方法和设备
CN105308469A (zh) * 2013-06-06 2016-02-03 皇家飞利浦有限公司 具有B0失真校正的并行MR成像以及使用规则化的SENSE重建的多点Dixon水-脂肪分离
GB201713444D0 (en) * 2017-08-22 2017-10-04 Siemens Healthcare Ltd A method of performing magnetic resonance imaging and a magnetic resonance apparatus
CN109477877A (zh) * 2016-07-14 2019-03-15 上海联影医疗科技有限公司 磁共振成像系统和方法
CN110346741A (zh) * 2018-04-06 2019-10-18 凯斯西储大学 使用多个脉冲序列类型的磁共振指纹的系统和方法
CN111505553A (zh) * 2019-07-12 2020-08-07 上海联影医疗科技有限公司 磁共振成像系统和方法

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7906964B2 (en) * 2007-01-29 2011-03-15 New York University Method and system for determining acquisition parameters associated with magnetic resonance imaging for a particular measurement time
JP6615594B2 (ja) * 2015-11-27 2019-12-04 キヤノンメディカルシステムズ株式会社 画像処理方法、画像処理装置及び磁気共鳴イメージング装置
US10261152B2 (en) * 2017-03-22 2019-04-16 Wisconsin Alumni Research Foundation System and method for confounder-corrected T1 measures using MRI
CN107369153B (zh) * 2017-07-18 2020-11-27 上海联影医疗科技股份有限公司 磁共振脂肪抑制图像获取方法、存储介质及扫描系统
US10641853B2 (en) * 2017-09-08 2020-05-05 The Regents Of The University Of California Systems and methods for ultrashort echo time actual flip angle imaging and variable repetition time magnetic resonance imaging

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101919696A (zh) * 2009-06-15 2010-12-22 通用电气公司 用于磁共振rf场测量的系统、方法和设备
CN105308469A (zh) * 2013-06-06 2016-02-03 皇家飞利浦有限公司 具有B0失真校正的并行MR成像以及使用规则化的SENSE重建的多点Dixon水-脂肪分离
CN109477877A (zh) * 2016-07-14 2019-03-15 上海联影医疗科技有限公司 磁共振成像系统和方法
GB201713444D0 (en) * 2017-08-22 2017-10-04 Siemens Healthcare Ltd A method of performing magnetic resonance imaging and a magnetic resonance apparatus
CN110346741A (zh) * 2018-04-06 2019-10-18 凯斯西储大学 使用多个脉冲序列类型的磁共振指纹的系统和方法
CN111505553A (zh) * 2019-07-12 2020-08-07 上海联影医疗科技有限公司 磁共振成像系统和方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
吴宗山，彭传勇，叶永泉，陈艳玲，张树恒，徐启兰，李运运，项杨，陈琪.基于多维度集成算法的T2*定量成像对颅内出血的诊断.放射学实践.2019,第34卷(第10期),第1071-1076页. *

Also Published As

Publication number Publication date
US20200371184A1 (en) 2020-11-26
US11287504B2 (en) 2022-03-29
CN114076911A (zh) 2022-02-22

Similar Documents

Publication Publication Date Title
CN111505553B (zh) 磁共振成像系统和方法
CN107126211B (zh) 磁共振成像方法和系统
CN114076911B (zh) 磁共振成像的方法及系统
CN111012349B (zh) 磁共振成像中信号表征确定的系统与方法
US20240264257A1 (en) Systems and methods for magnetic resonance imaging
CN112763957B (zh) 磁共振成像系统和方法
CN113126013B (zh) 图像处理系统和方法
WO2019041607A1 (en) SYSTEM AND METHOD FOR DETERMINING ROCK ANGLES IN MAGNETIC RESONANCE IMAGING
US11592506B2 (en) System and method for magnetic resonance imaging
US11940516B2 (en) Systems and methods for magnetic resonance imaging
US20220409084A1 (en) Systems and methods for magnetic resonance imaging
US11069063B2 (en) Systems and methods for noise analysis
US11980457B2 (en) Systems and methods for simultaneous multi-slice multitasking imaging
US11353527B2 (en) Systems and methods for waveform determination in magnetic resonance imaging
EP3451009B1 (en) System and method for determining flip angles of refocusing radio frequency pulses in magnetic resonance imaging
CN113050010B (zh) 噪音分析的系统、方法
US20240053419A1 (en) Systems and methods for magnetic resonance imaging
US11982728B2 (en) Systems and methods for magnetic resonance imaging
US20230019828A1 (en) Systems and methods for magnetic resonance imaging
US11946993B2 (en) Imaging method and system
CN117518051A (zh) 一种磁共振成像方法和系统
CN116205999A (zh) 磁共振t1映射系统和方法

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