CN110226098A - 动态对比增强的磁共振成像中的图像质量控制 - Google Patents
动态对比增强的磁共振成像中的图像质量控制 Download PDFInfo
- Publication number
- CN110226098A CN110226098A CN201780084811.6A CN201780084811A CN110226098A CN 110226098 A CN110226098 A CN 110226098A CN 201780084811 A CN201780084811 A CN 201780084811A CN 110226098 A CN110226098 A CN 110226098A
- Authority
- CN
- China
- Prior art keywords
- magnetic resonance
- resonance imaging
- calibration
- imaging system
- estimation
- 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.)
- Granted
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/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
- G01R33/543—Control of the operation of the MR system, e.g. setting of acquisition parameters prior to or during MR data acquisition, dynamic shimming, use of one or more scout images for scan plane prescription
-
- 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/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/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
Landscapes
- Physics & Mathematics (AREA)
- High Energy & Nuclear Physics (AREA)
- Condensed Matter Physics & Semiconductors (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Signal Processing (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Radiology & Medical Imaging (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
本发明提供了一种磁共振成像系统(100),包括用于存储机器可执行指令(140)和脉冲序列命令(142)的存储器(134)。脉冲序列命令被配置用于根据DCE磁共振成像协议控制磁共振成像系统。所述磁共振成像系统还包括用户界面(200)和用于控制所述控制磁共振成像系统的处理器(130)。机器可执行指令的执行使处理器:使用脉冲序列命令控制(500)磁共振成像系统,以针对变化的翻转角两次或更多次采集校准磁共振数据(144);将校准磁共振数据的每个采集重建(502)成校准图像(146)以创建一组可变翻转角图像(148);使用所述一组可变翻转角图像计算(504)T1映射(150);至少部分地使用T1映射计算(506)预定磁共振成像造影剂的造影剂校准(152);使用校准准确度模型计算(508)估计的校准误差(154),其描述造影剂校准和/或T1映射中的估计的误差,其中,校准准确度模型被配置用于使用所述一组可变翻转角图像计算估计的校准误差;并且如果估计的校准误差在预定校准误差范围之外,则在用户界面上显示(510)校准警告消息(202)。
Description
技术领域
本发明涉及磁共振成像,具体地涉及动态对比增强的磁共振成像。
背景技术
大的静态磁场由磁共振成像(MRI)扫描器使用以对准原子的核自旋,作为用于产生患者的身体内的图像的流程的部分。该大的静态磁场被称为B0场或主磁场。对象的各种数量或性质可以使用MRI来空间地测量。
一种MRI成像技术是T1加权成像,在不同翻转角度下采集的一系列T1加权图像可以用于计算T1弛豫时间的空间映射。存在可注射的造影剂,其可以改变对象内的局部T1弛豫时间。在注射造影剂之前,可以执行对象内的T1弛豫时间的映射。该初始T1弛豫时间图可以用于校准已知造影剂对对象的影响。然后可以以不同的时间间隔在注射造影剂之后采集一系列T1加权图像。该一系列T1加权图像和根据初始T1弛豫时间图创建的校准可以用于完成对对象内的造影剂的定量测量结果。然后可以使用示踪动力学建模依据生理组织特性来解释这些定量测量结果。一组技术是动态对比度增强(DCE)磁共振成像。期刊文章Sourbron,S.P.和Buckley,D.L.(2013),“Classic models for dynamic contrast-enhanced MRI”(NMR Biomed.,26:1004-1027,doi:10.1002/nbm.2940)评论了各种DCE磁共振成像技术和数据分析方法。
美国专利申请US2014/0010432公开了一种方法,其包括获得患者的图像数据。图像数据对应于来自成像采集的采集数据,所述成像采集来自患者的检查计划中的一组规划的图像采集。方法还包括基于成像实践指南利用处理器分析图像数据并产生指示分析的电子格式化数据。处理器基于指示分析的数据来生成指示对检查计划的改变的推荐的信号。
发明内容
本发明在独立权利要求中提供磁共振成像系统、方法和计算机程序产品。在从属权利要求中给出了实施例。
DCE MRI的困难在于在执行利用造影剂的成像之前的患者特异性协议,其中,选择各种扫描参数并首先执行校准。如果在校准或扫描参数的选择中出现误差,则操作者可能不会意识到误差,直到整个DCE MRI协议已经完成。这能够要求对象在很晚之后或甚至在不同的一天返回以进行重复检查。
实施例可以提供通过估计T1映射和/或造影剂校准中的校准误差来确保DCE图像的质量的手段。可以以变化的翻转角采集一系列T1加权图像,以形成一组可变翻转角图像。该组可变翻转角图像用于构建T1映射。比较所述一组可变翻转角图像的个体图像可以例如用于实现T1图中的噪声的估计。T1图中的该噪声估计可以用作一个准则或“校准误差”,其可以用于警告用户采取校正动作。可以将误差的传播应用于造影剂的校准方程,以计算造影剂校准其自身的“校准误差”的另一变型。这可以具有能够警告操作者应该重复校准或者应该改变特定扫描参数的优点。
在其他实施例中,可以添加额外的改善。例如,可以使用用于DCE MRI协议的扫描参数和校准误差的知识来计算预测的DCE图像质量。例如,扫描参数和校准误差可以与用于基于梯度回波的磁共振成像的稳态信号方程一起使用,以对预测的DCE图像质量进行前向建模。这样做的一个优点是在已经注射造影剂之前估计该DCE图像质量。如果DCE图像质量不足,仍可采取校正动作。
在其他实施例中,可以添加甚至更多改善。另一个困难是DCE MRI研究的结果能够取决于所选择的扫描参数。这通常很难正确进行,因为当一个扫描参数被改变时,其能够影响其他值。可以将当前扫描参数的值与中央数据库或甚至患者特异性数据库进行比较,以确保DCE图像是一致的。如果值在预定操作范围之外,则可以警告操作者和/或为操作者提供校正以实现扫描参数。
作为具体范例,扫描参数的两个范例是视场(FOV)和脉冲重复时间(TR)。操作者选择不正确扫描参数的常见情况会是对于需要增加的视场(FOV)或感兴趣区域的大患者的情况。操作者可以增加FOV,其可以引导用户增加采集矩阵以便保持空间分辨率。然而,这导致增加的TR和更差的时间分辨率,从而潜在地违反规定的范围。如果用户选择比通常更多的成像切片以覆盖更大的解剖区域-这也导致更差的时间分辨率,则能够出现类似的问题。
在一个方面中,本发明提供磁共振成像系统。磁共振成像系统包括用于存储机器可执行指令和脉冲序列命令的存储器。本文使用的脉冲序列命令包括命令或数据,所述命令或数据可以被转换成用于使磁共振成像系统采集某类磁共振数据的命令。脉冲序列命令实际上可以包含用于多于一种类型的磁共振成像数据采集的指令,以满足用于特定成像协议的数据。
脉冲序列命令被配置用于根据DCE或动态对比度增强磁共振成像协议来控制磁共振成像系统。在DCE磁共振成像中,校准磁共振成像系统对造影剂的响应,使得可以测量诸如血管外细胞外空间或其他组织参数(诸如转移常数)的数量的绝对量度。这些可以例如有用于检测肿瘤组织的代谢。
磁共振成像系统还包括用户界面。磁共振成像系统还包括用于控制磁共振成像系统的处理器。
机器可执行指令的执行使处理器使用脉冲序列命令来控制磁共振成像系统,以针对变化的翻转角两次或更多次采集校准磁共振数据。术语“校准磁共振数据”是指示用于完成T1映射的特定磁共振数据的标签,所述T1映射可用于定量校准造影剂。机器可执行指令的执行还使处理器将校准磁共振数据的每个采集重建为校准图像,以创建一组可变翻转角图像。对于变化的翻转角度,多次采集校准磁共振数据。然后将该校准磁共振数据重建为一组磁共振图像,其在本文中称为一组可变翻转角图像。
机器可执行指令的执行还使处理器使用所述一组可变翻转角图像计算T1映射。机器可执行指令的执行还使处理器至少部分地使用T1映射来计算针对预定磁共振成像造影剂的造影剂校准。T1映射可以用于计算造影剂校准,因为所述一组可变翻转角图像提供针对各种翻转角的基线测量结果。机器可执行指令的执行还使处理器使用校准准确度模型计算估计的校准误差,所述估计的校准误差描述造影剂校准和/或T1映射中的估计的误差。
校准准确度模型被配置用于使用一组可变翻转角图像来计算估计的校准误差。校准准确度模型是使用可变翻转角图像中的个体图像来估计造影剂校准或T1映射中的误差的模型。例如,可以将各种图像彼此进行比较,或者可以将其他统计量度应用于各种图像以估计测量结果中的信噪比。这可以用于提供造影剂和/或T1映射中的信噪比或其他误差的准确估计。
如果估计的校准误差在预定校准误差范围之外,则机器可执行指令的执行还使处理器在用户界面上显示校准警告消息。这可能是有益的,因为DCE磁共振图像的采集能够涉及大量时间,其中,首先校准磁共振数据被采集并用于进行校准。在执行校准之后,例如可以利用对象中的预定磁共振成像造影剂来执行成像。在利用主题中的造影剂采集图像上浪费时间之前估计校准中的误差能够是有益的。
在另一实施例中,机器可执行指令的执行还使处理器接收一组扫描参数,以用于在脉冲序列命令的执行期间控制磁共振成像。例如,扫描参数可以指示对象的检查的特定细节。这例如可以指示从其采集数据的感兴趣的位置或区域,或者可以存在在采集某类磁共振数据之前由用户设置的其他参数。机器可执行指令的执行还使处理器使用所述一组扫描参数、估计的校准误差并且使用用于基于梯度回波的磁共振成像的稳态信号方程来计算预测的DCE图像质量。如果预测的DCE图像质量在预定图像质量范围之外,则机器可执行指令的执行还使处理器在用户界面上显示采集警告消息。
在该实施例中,其能够是有益的,因为可以甚至在采集DCE磁共振图像之前预测预测的DCE图像质量。这是通过使用特定于特定检查的所述一组扫描参数和估计的校准误差和通过使用稳态信号方程的建模来完成的。这可以有用于预测DC图像质量,例如预测的信噪比和/或根据造影剂测量结果导出的一个或多个特定造影剂测量结果的预测准确度。这可以提供在DC图像中的噪声水平或图像质量为差的情况下防止重新检查的有效手段。
在另一个实施例中,机器可执行指令的执行还使处理器利用脉冲序列命令控制磁共振成像系统,以采集DCE磁共振数据的序列。机器可执行指令的执行还使处理器使用DCE磁共振数据和造影剂校准将DCE磁共振数据的序列重建为一系列DCE磁共振图像。
该实施例能够是有益的,因为预测的DCE图像质量和估计的校准误差的组合检查可以帮助确保一系列DCE磁共振图像具有更高的图像质量和/或在造影剂测量结果或根据造影剂测量结果导出的数量中具有更高的预测准确度。
在另一实施例中,机器可执行指令的执行还使处理器将扫描参数存储在存储器中的日志文件中。
机器可执行指令的执行还使处理器重复读取日志文件以获得一组当前扫描参数。机器可执行指令的执行还使处理器将当前扫描参数与一组预定扫描参数进行比较,以确定当前扫描参数是否在预定操作范围之外。如果当前扫描参数在预定操作范围之外,则机器可执行指令的执行还使处理器在用户界面上显示操作参数警告。
该实施例能够是有益的,因为其可以提供检查扫描参数是否合适或者在用于采集DCE磁共振图像或校准磁共振数据的预定止指南内的自动手段。通常,磁共振成像参数是相互关联的。对于操作磁共振成像系统的技术人员而言选择适当的操作参数可能是困难的或是不可能的。通过从日志文件中读取扫描参数,机器可执行指令可以有效地监测扫描参数并自动将其与预定扫描参数进行比较。
在另一实施例中,预定扫描参数由以下中的任一项确定:公布的标准、医师推荐或命令、患者简档、患者数据库及其组合。该实施例能够是有益的,因为其可以提供可以比较扫描参数的多个位置。
例如,这可以有益于使用特定患者简档或患者数据库提供质量控制,以使DCE磁共振图像与先前或后续检查相比能够更好。
在另一实施例中,由处理器作为背景过程执行日志文件的重复读取和当前扫描参数的比较。这可以具有以下优点:其自动且透明地起作用并且不需要由磁共振成像系统的操作者监测。
在另一个实施例中,机器可执行指令的执行还使处理器响应于以下中的任一项而在用户界面上显示一组建议的扫描参数改变:显示操作者参数警告、采集警告消息、校准警告消息及其组合。该实施例能够是有益的,因为其可以帮助磁共振成像系统的操作者更好地顺应DCE扫描的标准操作条件,或者确保当前采集的DCE图像和导出的数据与先前采集的DCE磁共振图像和/或导出的数据相比可以更好。
在另一实施例中,机器可执行指令的执行还使处理器至少部分地使用以下中的任一项来确定所述一组建议的扫描参数改变:查找表、先前扫描数据库、扫描参数之间的相互关系的模型及其组合。该实施例能够是有益的,因为其全部可以提供有效地向磁共振成像系统的操作者提供建议的参数的手段。
在另一实施例中,校准警告包括以下中的任一项:估计的校准误差在预定校准范围之外的警告、建议的校正动作及其组合。该实施例能够是有益的,因为其可以为磁共振成像系统的操作者提供一种在流程已经进展很远之前校正磁共振DCE图像的采集的手段。
在另一实施例中,校准准确度模型被配置用于使用以下中的任一项来估计估计的校准误差:直方图基础噪声估计模型、基于小波变换的噪声估计模型、或使用由所述一组可变翻转角图像的不同组合的相减构建的数据集来估计噪声的噪声估计模型。在这些范例的每个中,T1映射中的信噪比和/或造影剂校准中的估计的误差是通过使用所述一组可变翻转角产生和估计信噪比来进行的。然后使用图像中的噪声的该估计来计算估计的校准误差。
在另一实施例中,机器可执行指令的执行还使处理器使用脉冲序列命令来控制磁共振成像系统以采集B1映射磁共振数据。例如,可以使用制作B1磁共振图的任何常用方法。B1磁共振图的知识可以有益于确定更准确的T1映射或造影剂校准。机器可执行指令的执行还使处理器使用B1映射磁共振数据计算B1图。还至少部分地使用B1图来计算造影剂校准。
在另一方面中,本发明提供了一种计算机程序产品,包括用于由控制磁共振成像系统的处理器执行的机器可执行指令。磁共振成像系统包括用户界面。机器可执行指令的执行使处理器使用脉冲序列命令来控制磁共振成像系统,以针对变化的翻转角两次或更多次采集校准磁共振数据。脉冲序列命令被配置用于根据DCE磁共振成像协议控制磁共振成像系统。机器可执行指令的执行还使处理器将校准磁共振数据的每个采集重建为校准图像,以计算一组可变翻转角图像。
所述一组可变翻转角图像由个体校准图像组成。机器可执行指令的执行还使处理器使用所述一组可变翻转角图像计算T1映射。机器可执行指令的执行还使处理器至少部分地使用T1映射来计算预定磁共振成像造影剂的造影剂校准。机器可执行指令的执行还使处理器使用校准准确度模型计算估计的校准误差,所述估计的校准误差描述造影剂校准和/或T1映射中的估计的误差。校准准确度模型被配置用于使用所述一组可变翻转角图像计算估计的校准误差。如果估计的校准误差在预定校准误差范围之外,则机器可执行指令的执行还使处理器在用户界面上显示校准警告消息。
之前已经讨论了该系统的优点。
在另一方面中,本发明提供了一种操作磁共振成像系统的方法。所述磁共振成像系统包括用户界面。所述方法包括将对象放置在磁共振成像系统中。所述方法还包括使用脉冲序列命令控制磁共振成像系统,以针对变化的翻转角两次或更多次采集校准磁共振数据。脉冲序列命令被配置用于根据DCE磁共振成像协议控制磁共振成像系统。所述方法还包括将校准磁共振数据的每个采集重建为校准图像以创建一组可变翻转角图像。所述方法还包括使用所述一组可变翻转角图像计算T1映射。
所述方法还包括至少部分地使用T1映射计算针对预定磁共振成像造影剂的造影剂校准。所述方法还包括使用校准准确度模型计算估计的校准误差,所述估计的校准误差描述造影剂和/或T1映射中的估计的误差。校准准确度模型被配置用于使用所述一组可变翻转角图像计算估计的校准误差。所述方法还包括如果估计的校准误差在预定校准误差范围之外,则在用户界面上显示校准警告消息。之前已经讨论过该系统的优点。
在另一实施例中,所述方法还包括接收一组扫描参数,以用于在执行脉冲序列命令期间控制磁共振成像。所述方法还包括使用所述一组扫描参数、估计的校准误差和使用用于创建基于回波的磁共振成像的稳态信号方程来计算预定DCE图像质量。所述方法还包括:如果预测的DCE图像质量在预定图像质量范围之外,则在用户界面上显示采集警告消息。之前已经讨论了这些额外的步骤的优点。
在范例中,所述方法还包括向对象注射预定磁共振成像造影剂。所述方法还包括利用脉冲序列命令控制磁共振成像系统,以根据DCE磁共振成像协议采集DCE磁共振数据的序列。例如,可以以特定时间间隔采集DCE磁共振数据的图像。所述方法还包括使用DCE磁共振数据和造影剂校准将DCE磁共振数据的序列重建为一系列DCE磁共振图像。可以存在另外的后处理步骤,其中,一系列DCE磁共振图像用于计算导出的数量,例如血管外细胞外空间或组织参数,例如转移常数。
应当理解,可以组合本发明的前述实施例中的一个或多个,只要组合的实施例不是相互排斥的。
如本领域的技术人员将认识到的,本发明的各个方面可以实现为装置、方法或计算机程序产品。相应地,本发明的各个方面可以采取完全硬件实施例、完全软件实施例(包括固件、常驻软件、微代码等)或组合软件和硬件方面的实施例(在本文中总体上全部可以被称为“电路”、“模块”或“系统”)的形式。此外,本发明的各个方面可以采取实现在一个或多个计算机可读介质中的计算机程序产品的形式,所述一个或多个计算机可读介质具有实现在其上的计算机可执行代码。
可以利用一个或多个计算机可读介质的任何组合。所述计算机可读介质可以是计算机可读信号介质或计算机可读存储介质。本文使用的“计算机可读存储介质”涵盖任何可以存储可由计算设备的处理器执行的指令的有形存储介质。可以将计算机可读存储介质称为计算机可读非暂态存储介质。也可以将计算机可读存储介质称为有形计算机可读介质。在一些实施例中,计算机可读存储介质还可以能够存储可以由计算设备的处理器访问的数据。计算机可读存储介质的范例包括但不限于:软盘、磁硬盘驱动器、固态硬盘、闪速存储器、USB拇指驱动器、随机存取存储器(RAM)、只读存储器(ROM)、光盘、磁光盘以及处理器的寄存器文件。光盘的范例包括压缩盘(CD)和数字通用盘(DVD),例如,CD-ROM、CD-RW、CD-R、DVD-ROM、DVD-RW或DVD-R盘。术语计算机可读存储介质还指能够经由网络或通信链路由计算机设备访问的各种类型的记录介质。例如,可以在调制调解器、因特网或局域网上检索数据。可以使用任何适当介质发送实现在计算机可读介质上的计算机可执行代码,所述任何适当介质包括但不限于无线的、有线的、光纤线缆的、RF等或者前面的任何合适的组合。
计算机可读信号介质可以包括具有实现在其中的计算机可执行代码的传播的数据信号,例如,在基带中或作为载波的部分。这样的传播的信号可以采取任何各种形式,包括但不限于电磁的、光学的或它们的任何合适的组合。计算机可读信号介质可以是这样的任何计算机可读介质:不是计算机可读存储介质,并且能够传达、传播或传输由指令执行系统、装置或设备使用或与指令执行系统、装置或设备结合使用的程序。
“计算机存储器”或“存储器”是计算机可读存储介质的范例。计算机存储器是可由处理器直接访问的任何存储器。“计算机存储设备”或“存储设备”是计算机可读存储介质的另外的范例。计算机存储设备可以是任何易失性或非易失性计算机可读存储介质。
如本文中所使用的“处理器”涵盖能够执行程序或机器可执行指令或计算机可执行代码的电子部件。对包括“处理器”的计算设备的引用应当被解读为能够包含多于一个的处理器或处理核。所述处理器可以例如是多核处理器。处理器也可以指在单个计算机系统之内的或分布在多个计算机系统之间的处理器的集合。术语计算设备也应当被解读为能够指每个包括一个或多个处理器的计算设备的集合或网络。计算机可执行代码可以由可以在相同的计算设备之内或甚至可以分布在多个计算设备之间的多个处理器来执行。
计算机可执行代码可以包括令处理器执行本发明的方面的机器可执行指令或程序。用于执行针对本发明的方面的操作的计算机可执行代码可以以一个或多个编程语言的任何组合来编写并且被编译为机器可执行指令,所述一个或多个编程语言包括诸如Java、Smalltalk、C++等的面向对象的编程语言以及诸如“C”编程语言或相似编程语言的常规过程性编程语言。在一些实例中,所述计算机可执行代码可以采取高级语言的形式或者采取预编译的形式并且结合在工作时生成机器可执行指令的解读器一起被使用。
所述计算机可执行代码可以完全在用户的计算机上、部分在用户的计算机上(作为独立的软件包)、部分在用户的计算机上并且部分在远程计算机上、或完全在远程计算机或服务器上执行。在后一种情形下,所述远程计算机可以通过包括局域网(LAN)或广域网(WAN)的任何类型的网络连接到用户的计算机,或者可以(例如,通过使用因特网服务提供商的因特网)对外部计算机进行连接。
参考根据本发明的实施例的方法、装置(系统)和计算机程序产品的流程图、图示和/或方框图来描述本发明的方面。应理解,当可应用时,能够通过采取计算机可执行代码的形式的计算机程序指令来实施流程图、图示和/或方框图的方框的每个方框或部分。还应理解,当互不排斥时,可以组合不同流程图、图示和/或方框图中的方框的组合。这些计算机程序指令可以被提供给通用计算机、专用计算机或产生机器的其他可编程数据处理装置的处理器,使得经由计算机或其他可编程数据处理装置的处理器执行的指令创建用于实施在流程图和/或一个或多个方框图框中指定的功能/动作的单元。
这些计算机程序指令还可以存储在计算机可读介质中,所述计算机可读介质能够指引计算机、其他可编程数据处理装置或其他设备以特定的方式来工作,使得在计算机可读介质中存储的指令产生包括实施在流程图和/或一个或多个方框图框中指定的功能/动作的指令的制品。
所述计算机程序指令还可以加载到计算机、其他可编程数据处理装置或其他设备上,以令在计算机、其他可编程装置或其他设备上执行一系列操作步骤,从而产生计算机实施的过程,使得在计算机或其他可编程装置上执行的指令提供用于在流程图和/或一个或多个方框图框中指定的功能/动作的过程。
如本文所使用的“用户接口”是允许用户或操作人员与计算机或计算机系统交互的接口。“用户接口”还可以被称为“人机接口设备”。用户接口可以向操作人员提供信息或数据和/或从操作人员接收信息或数据。用户接口可以使得来自操作人员的输入能够被计算机接收并且可以从计算机向用户提供输出。换言之,所述用户接口可以允许操作人员控制或操控计算机,并且所述接口可以允许计算机指示操作人员的控制或操控的效果。显示器或图形用户接口上的数据或信息的显示是向操作人员提供信息的范例。通过键盘、鼠标、跟踪球、触摸板、指点杆、图形输入板、操纵杆、游戏手柄、网络摄像头、耳机、踏板、有线手套、遥控器和加速度计对数据的接收全都是实现对来自操作人员的信息或数据的接收的用户接口部件的范例。
如本文所使用的“硬件接口”涵盖使得计算机系统的处理器能够与外部计算设备和/或装置交互和/或控制外部计算设备和/或装置的接口。硬件接口可以允许处理器将控制信号或指令发送到外部计算设备和/或装置。硬件接口也可以使得处理器能够与外部计算设备和/或装置交换数据。硬件接口的范例包括但不限于:通用串行总线、IEEE 1394端口、并行端口、IEEE 1284端口、串行端口、RS-232端口、IEEE-488端口、蓝牙连接、无线局域网连接、TCP/IP连接、以太网连接、控制电压接口、MIDI接口、模拟输入接口以及数字输入接口。
如本文所使用的“显示器”或“显示设备”涵盖适于显示图像或数据的输出设备或用户接口。显示器可以输出视觉、音频和/或触觉数据。显示器的范例包括但不限于:计算机监测器、电视屏幕、触摸屏、触觉电子显示器、盲文屏幕、阴极射线管(CRT)、存储管、双稳态显示器、电子纸、矢量显示器、平板显示器、真空荧光显示器(VF)、发光二极管(LED)显示器、电致发光显示器(ELD)、等离子体显示板(PDP)、液晶显示器(LCD)、有机发光二极管显示器(OLED)、投影仪和头戴式显示器。
磁共振(MR)数据在本文中被定义为在磁共振成像扫描期间使用磁共振装置的天线所记录的由原子自旋发射的射频信号的测量结果。磁共振数据是医学成像数据的范例。磁共振(MR)图像在本文中被定义为磁共振成像数据内包含的解剖数据的经重建的二维或三维可视化。
附图说明
在下文中将仅通过范例并且参考附图描述本发明的优选实施例,在附图中:
图1图示了磁共振成像系统的范例;
图2图示了用户接口的范例;
图3图示了用户接口的另外的范例;
图4图示了用户接口的另外的范例;
图5示出了图示使用图1的磁共振成像系统的方法的流程图;
图6示出了图示使用图1的磁共振成像系统的另外的方法的流程图;
图7示出了图示使用图1的磁共振成像系统的另外的方法的流程图;并且
图8包含针对其功能的至少部分来表示磁共振成像系统的图。
附图标记列表
100 磁共振成像系统
104 磁体
106 磁体的膛
108 成像区
109 感兴趣区域
110 磁场梯度线圈
112 磁场梯度线圈电源
114 射频线圈
116 收发器
118 对象
120 对象支撑件
126 计算机系统
128 硬件接口
130 处理器
132 显示器
134 计算机存储器
140 机器可执行指令
142 脉冲序列命令
144 校准磁共振数据
146 校准图像
148 一组可变翻转角图像
150 T1映射
152 造影剂校准
154 估计的校准误差
156 一组扫描参数
158 预测的DCE图像质量
160 日志文件
162 一组预定扫描参数
200 图形用户界面
202 校准警告消息
300 采集警告消息
400 操作参数警告消息
402 建议的校正动作
800 质量控制和工作流程系统
802 B1图
804 注射
806 DCE磁共振图像
808 协议和质量@t1
810 准确度模型和比较器
812 当前扫描参数
具体实施方式
在这些附图中相似编号的元件或为等价元件或执行相同的功能。如果功能等价,则先前已经论述的元件将不必要在后面的附图中论述。
图1示出了具有磁体104的磁共振成像系统100的范例。磁体104是具有穿过其的膛106的超导圆柱型磁体。使用不同类型的磁体也是可能的;例如也能够使用分裂式圆柱形磁体和所谓的开放式磁体两者。分裂式圆柱形磁体类似于标准圆柱形磁体,除了低温恒温器已经被分裂成两部分以允许进入磁体的等平面,这样的磁体可以例如与带电粒子束治疗结合使用。开放式磁体具有两个磁体部分,一个在另一个上方,其间具有足够大的空间以接收对象:两个部分区的布置类似于亥姆霍兹线圈的布置。开放式磁体是受欢迎的,因为对象较少受限制。在圆柱形磁体的低温恒温器内部存在超导线圈的集合。在圆柱形磁体104的膛106内,存在成像区108,其中,磁场足够强且均匀以执行磁共振成像。感兴趣区域109被示出在成像区108内。对象118被示出为由对象支撑体120支撑,使得对象118的至少部分在成像区108和感兴趣区域109内。
在磁体的膛106内,还存在一组磁场梯度线圈110,所述一组磁场梯度线圈用于采集磁共振数据以在空间上编码磁体104的成像区108内的磁自旋。磁场梯度线圈110连接到磁场梯度线圈电源112。磁场梯度线圈110旨在是代表性的。通常,磁场梯度线圈110包含三个独立的线圈组,以用于在三个正交的空间方向上进行空间编码。磁场梯度电源向磁场梯度线圈供应电流。供应给磁场梯度线圈110的电流被控制为时间的函数并且可以是斜变或脉冲的。
与成像区108相邻的是射频线圈114,以用于操纵成像区108内的磁自旋的取向并且用于从也在成像区108内的自旋接收无线电发射。射频天线可以包含多个线圈元件。射频天线也可以被称为通道或天线。射频线圈114连接到射频收发器116。射频线圈114和射频收发器116可以由独立的发射线圈和接收线圈以及独立的发射器和接收器代替。应理解,射频线圈114和射频收发器116是代表性的。射频线圈114也旨在表示专用发射天线和专用接收天线。类似地,收发器116也可以表示独立的发射器和接收器。射频线圈114也可以具有多个接收/发射元件,并且射频收发器116可以具有多个接收/发射通道。例如,如果诸如SENSE的并行成像技术被执行,则射频能够114将具有多个线圈元件。
收发器116和梯度控制器112被示出为连接到计算机系统126的硬件接口128。计算机系统还包括与硬件系统128、存储器134和显示器132通信的处理器130。存储器134可以是处理器130可访问的存储器的任何组合。这可以包括诸如主存储器、高速缓存存储器以及非易失性存储器(诸如闪存RAM、硬盘驱动器或其他存储设备)的事物。在一些范例中,存储器130可以被认为是非瞬态计算机可读介质。
计算机存储器被示出为包含使得处理器130能够控制磁共振成像系统100的操作和功能的一组机器可执行指令140。计算机存储器134还被示出为包含脉冲序列命令142,脉冲序列命令142使得处理器130能够控制磁共振成像系统100根据DCE磁共振成像协议来采集磁共振数据。计算机存储器134还被示出为包含通过利用脉冲序列命令142的部分控制磁共振成像系统所采集的校准磁共振数据144。校准磁共振数据144被采集多次并且针对各种翻转角被采集。
计算机存储器134被示出为包含一幅校准图像146,校准图像146是根据针对特定翻转角的校准磁共振数据144的一次采集重建的。计算机存储器134还被示出为包含一组可变翻转角图像148,其由针对不同翻转角采集的一组校准图像146组成。计算机存储器134还被示为包含根据所述一组可变翻转角图像计算的T1映射150。计算机存储器134还被示出为包含使用T1映射150的预定磁共振成像造影剂的造影剂校准152。
计算机存储器134还被示出为包含估计的校准误差154,估计的校准误差154通过将所述一组可变翻转角图像148中的个体图像彼此比较来计算。估计的校准误差可以例如是T1映射或造影剂校准中的信噪比的估计,或者可以是例如从其导出的数量。计算机存储器134还被示出为包含一组扫描参数156,其可以用于使用如由脉冲序列命令142定义的DCE磁共振成像协议来定制磁共振数据的采集。
计算机存储器134还被示为包含使用所述一组扫描参数156、估计的校准误差154以及例如可以是用于磁共振成像的稳态信号方程的模型计算的预测的DCE图像质量158,所述模型用于对基于梯度回波的磁共振信号进行建模。计算机存储器134还被示为包含日志文件160,日志文件160包含磁共振成像系统的当前操作状态。例如,当特定扫描参数156正用于采集磁共振数据或规划的磁共振数据时,其被存储在日志文件160中。
计算机存储器还被示出为包含一组预定扫描参数162,可以将存储在日志文件160中的所述一组扫描参数156与所述一组预定扫描参数进行比较。
图2示出了具有图形用户界面200的显示器132的范例。当估计的校准误差154高于预定量时,图形用户界面200可以显示校准警告消息202。
图3示出了如图2所示的显示器132和图形用户界面200。在这种情况下,图形用户界面200响应于预测的DCE图像质量158在预定图像质量范围之外而显示采集警告消息300。
图4示出了图形用户界面200中的显示的另一视图。在该范例中,图形用户界面200显示操作参数警告消息400,操作参数警告消息400响应于日志文件160中的当前扫描参数在预定操作范围之外而显示。图形用户界面200还可以显示一组建议的校正动作402或一组建议的扫描参数,以校正引起操作参数警告消息400的误差。
图5示出了图示操作图1的磁共振成像系统100的方法的流程图。首先,在步骤500中,使用脉冲序列命令142来控制磁共振成像系统100,以采集校准磁共振数据144。针对变化的翻转角度多次采集校准磁共振数据144。接下来在步骤502中,校准磁共振数据144的每次采集被重建为校准图像146。这被完成以创建一组可变翻转角图像148。接下来在步骤504中,使用所述一组可变翻转角图像148来计算T1映射150。
然后在步骤506中,至少部分地使用T1映射150计算针对预定磁共振成像造影剂的造影剂校准152。接下来在步骤508中,使用校准准确度模型计算估计的校准误差154,估计的校准误差154描述造影剂校准152和/或T1映射150中的估计的误差。校准准确度模型被配置用于使用所述一组可变翻转角图像148来计算估计的校准误差。最后,在步骤510中,如果估计的校准误差在预定校准误差范围之外,则在用户界面200上显示校准警告消息202。
图6示出了图示操作图1的磁共振成像系统100的另一方法的流程图。图6的方法类似于图5,具有若干额外的步骤。图6中的方法从步骤600开始。在步骤600中,接收一组扫描参数156,以用于控制在执行脉冲序列命令142期间的磁共振成像。然后方法行进到步骤500,如图5所图示的。在执行步骤510之后,执行步骤602。在步骤602中,使用所述一组扫描参数156、估计的校准误差154和用于对稳态信号方程建模以创建基于回波的磁共振成像的模型来计算预测的DCE图像质量158。最后在步骤604中,如果预测的DCE图像质量158在预定图像质量范围之外,则在用户界面200上显示采集警告消息300。
图7示出了图示操作图1的磁共振成像系统100的另一方法的流程图。图7中所图示的方法类似于图6中所图示的方法。图7中的方法开始于步骤600,再次如图6所描述的。方法然后行进到步骤700。在步骤700中,扫描参数被存储在计算机存储器134中的日志文件中。接下来在步骤702中,在操作和采集磁共振成像数据期间重复读取日志文件以获得一组当前扫描参数。步骤702可以重复多次。在执行步骤702之后,执行步骤704,并且将如从日志文件读取的当前扫描参数与所述一组预定扫描参数162进行比较,以确定当前扫描参数是否在预定操作范围之外。接下来在步骤706中,如果当前扫描参数在预定操作范围之外,则在用户界面200上显示操作者参数警告400。方法然后行进到步骤500,如图6所图示的。然后执行步骤500、502、504、506、508、510、602和604,如图6所图示的。
范例可以提供在标准化的定量DCE-MRI测量的定义、采集和后处理期间用于实时主动质量控制和工作流程支持的软件或机器可执行指令。其可以包括三个质量控制水平中的一个或多个:1)其确保顺应预定义的标准化指南以及后续检查的可比性。2)其便于图像质量的实时自动分析和预测。3)其使用正向建模预测定量后处理的准确度。流程期间或测量或预测的质量中的显着偏差及时报告给操作者,包括校正动作的建议。这些特征便于DCE-MRI数据的标准化、鲁棒和可重复的采集,以及具有受控定量准确度的导出参数图的随后的产生。
定量磁共振成像(MRI)被广泛认为是朝向个性化精准医学的范式转变,因为其使用绝对组织性质而不是定性对比度加权图像的视觉解释来促进病理学的评价。一项重要的技术是动态对比增强的MRI(DCE-MRI),其可用于量化组织参数,诸如转移常数Ktrans或血管外细胞外空间的分数体积,ve。这些参数已经被示出为有用,作为大量的研究中的肿瘤代谢的生物标记。定量DCE-MRI的主要挑战是由于在数据采集和后处理两者中标准化和质量控制的缺乏的结果的高度变异性。因此,不同研究的结果难以比较。因此,定量DCE-MRI目前主要在研究医院中执行,并且在标准临床设置中的采用仍然是有限的。
为了克服这些挑战,如定量成像生物标记联盟(QIBA)这样的组织寻求通过提供以不同流程的详细推荐为特征的概况来减少定量DCE-MRI的变异性。
即使利用标准化的协议推荐或指南,获得定量参数图的整个流程仍然很复杂,并且在采集和后处理期间包括多个关键且容易出错的步骤。因此,DCE-MRI在临床研究或临床常规中的应用面临着若干挑战:
大多数标准化指南包括MR采集参数的允许范围的详细列表。在临床常规中,确保根据指南执行所有测量是具有挑战性的。例如,特定患者解剖结构能够需要调节MR序列参数。在这种情况下,用户面对调节参数而不违反指南的困难任务。类似地,在后续检查的情况下,理想情况下,所有采集参数针对所有测量应该相同,以确保可比性。然而,先前检查中这些参数的精确值通常不能直接对用户可用。
在临床设置中,测量结果的准确度和精度能够受到各种事件的影响,例如硬件缺陷或线圈错位。在这些情况下,即使根据标准化指南执行采集,也能够损害所采集的图像的质量。这样的情况目前需要对图像的仔细视觉检查。
在测量的时间处,用户不知道特定参数设置或图像质量对组织参数量化的影响。由于次优参数设置的问题通常仅将在后处理期间展现,即在测量之后很久展现。这通常导致患者重新排程以获得足够的定量结果。
因此,范例可以实现主动质量控制和工作流程支持系统,其可以包括以下三个水平中的一个或多个:
在第一水平中,连续监测MR测量协议的状态。检测到根据标准化指南与预定义的协议的任何显着偏差,并向用户发布警告。在对治疗控制的后续检查的情况下,如果检测到与先前检查的设置的偏差,则通过发布警告来确保结果的可比性。
在第二水平中,通过实时分析来自于已经完成的扫描的图像来确保足够的图像质量。在低图像质量的情况下,向用户发布警告,包括修改采集的建议。由于定量DCE-MRI检查包括2次或更多次扫描,因此在开始后续扫描之前,即在修改仍然可能时,可以有利地使用在第一扫描中获得的信息来查找问题。具体地,可以在造影剂注射之前完成必要的修改,之后扫描操作是时间关键的并且不能重复。
在第三水平中,使用正向建模在检查期间连续预测定量后处理的预期准确度(DCE图像质量)。在低预期准确度的情况下,显示警告以及协议修改的建议。
实时主动质量控制基于两个软件(SW)部件,并且在下面进一步描述。
操作MR系统(“MR控制台”)和/或执行后处理步骤(“MR工作站”)的SW的修改,其自动地将关于当前、规划和完成的扫描的所有所需信息供应给质量控制SW。
质量和工作流程控制SW(在MR控制台或MR工作站上运行或在单独的监测系统,例如平板计算机上运行),其实现主动质量控制,包括对采集的数据的实时评价和对用户的反馈。
作为第一步骤,MR控制台和/或MR工作站上的额外的SW将关于当前测量协议的信息(例如MR序列参数、扫描状态、处理参数、应用的数据模型和患者信息)供应给质量和工作流程控制SW。该信息可以例如包括在MR控制台/工作站上创建的标准状态记录机构中,并且对应的日志文件可以由质量和工作流程控制SW传递和评价。此外,从完成的扫描获得的重建图像数据被传递到质量控制SW(例如,通过临时磁盘存储),以允许诸如信噪比(SNR)的图像质量度量的实时分析。
最后,质量和工作流程控制SW可以使用本地患者特异性数据库来存储关于测量协议的信息,以允许比较相同患者的后续检查的条件。
在第二步骤中,质量和工作流程控制SW被实施为旨在与MR采集或图像后处理并行运行的额外的程序。这可以在MR控制台上或在MR工作站上实施,或者在单独的监测系统(例如平板计算机)上实施。图8中示出了预期功能的可能实施方式的示意性概述。
图8从功能的角度表示磁共振成像系统100。框100表示由磁共振成像系统100采集的多个测量结果。框800表示质量控制和工作流程系统800,其可以例如由机器可执行指令实施。框126表示由计算机系统126执行的工作流程。质量控制工作流程系统800的输出可以例如是多个警告和报告,其可以等价于校准警告消息202、采集警告消息300和操作者参数警告消息400。
磁共振成像系统可以采集各种数据。这可以包括B1图802,构建成T1映射的可变翻转角图像150还可以表示造影剂804到患者或对象118中的注射和一系列DCE磁共振图像806的采集。然而,应注意,在执行任何注射之前执行各种警告消息202、300、400的计算和显示。质量控制和工作流程支持系统800可以具有多个部件。其可以具有患者特异性数据库162、标准和指南以及协调,并且协议162可以全部组合成一组预定扫描参数。
例如,可以将实际协议或当前扫描参数812输入到准确度模型810中,并且使用各种协议和图像质量808的比较器810具有输入。准确度模型810和比较器可以对应于校准准确度模型和用于创建基于回波的磁共振成像的稳态信号方程,以及在当前扫描参数与预定扫描参数之间进行比较的代码。在一些范例中,质量控制和工作流程支持系统作为背景过程运行,使得操作者不知道其存在,除非需要发布警告或报告202、300、404。
在典型的定量DCE-MRI检查中,示意性地描绘在图8的顶部:其包括B1映射扫描、稍后用于T1映射(造影剂浓度的校准所需要的)的一系列可变翻转角(VFA)扫描,以及用于在采集动态系列(DCE扫描)期间的造影剂注射。
由SW使用的实际协议使用预定义的外部标准、推荐或指南或通过研究特异性协议来初始化,所述研究特异性协议定义质量目标,例如允许的序列参数范围和期望的图像质量度量,例如,SNR、T1/B1映射中的容忍误差的边界或基于模型的定量后处理中的预期变化的边界。
SW从MR控制台或工作站接收连续数据流(例如,通过从日志文件中连续提取数据)以构建当前测量协议及其质量的内部实时表示。该信息被供应给不同的准确度模型,其结果与SW的实际概况进行比较。在准则顺应的情况下,准确度模型简单地产生当前一组序列参数。备选地,如果检测到与先前检查的设置的偏差,则通过发布警告来确保结果的可比性。可以采用不同的用户配置文件来显示警告,这取决于其严重性和测量目的。
用于定义合适和及时警告的中心指南是将用户交互减少到最小,以免分散标准化工作流程的注意力。重要的是,在扫描和后处理条件仍在预期范围内时,或者当前没有可能的校正动作时,质量和工作流程支持SW应保持不可见(无用户交互)。没有实时缓解的质量问题应在随后的性能分析中报告(例如,用于医院管理)。应在背景中执行用户动作的序列、图像质量和后处理分析结果以及实际发布的警告的大量记录,以便允许针对不同利益相关者(例如技术专家、放射科医师、医院经理)的利益在多个水平上进行性能分析。
作为实时图像质量分析的范例,一旦完成扫描的执行,MR控制台上的额外的SW就自动将重建图像供应给主动质量控制SW。在VFA扫描(其通常涉及信号平均(信号平均数(NSA)≥2))的情况下,在平均之前重建个体图像并将其发送到质量控制SW。然后,准确度模型实现这些个体图像的相减以获得可以用于确定平均图像中的噪声水平的数据集。备选地,可以应用基于自动直方图或基于小波变换的用于噪声估计的方法。同时,使用天然T1图以及所检查的感兴趣器官的典型T1值的查找表两者来确定针对感兴趣组织的近似信号水平。将这两个结果组合产生感兴趣器官的SNR。然后使用该VFA-SNR以及选定的采集参数(NSA、翻转角、重复时间等)和用于基于梯度回波的序列的稳态信号方程来预测DCE扫描的预期SNR。将这两个计算的SNR值与预定义的概况进行比较,以检测异常低的图像质量。除了该基于SNR的评估之外,可以实现基于与整体伪影水平相关的度量的其他准确度模型以评价图像质量。
最后,为了确保后处理的质量控制,准确度模型预测根据各种图像导出的定量参数的准确度和精度。使用蒙特卡罗方法确定T1映射的准确度和精度,所述蒙特卡罗方法采用不同VFA扫描的先前确定的SNR。闭式表达式可以用于浓度测量准确度和精度的时间有效计算。这样,可以在不可重复的造影剂注射之前校正不充足的序列参数设置或图像伪影。
在由用户对协议的每次修改之后,SW重建其测量协议的内部表示并再次执行所有准确度模型。如果检测到偏差(例如,如果参数值在其允许范围之外,或者检测到的图像SNR低于期望),则向用户显示视觉警告,并向操作者示出针对校正的建议。
尽管已经在附图和前面的描述中详细说明和描述了本发明,但这样的说明和描述被认为是说明性或示范性的而非限制性的;本发明不限于公开的实施例。
本领域技术人员通过研究附图、说明书和权利要求书,在实践要求保护的本发明时能够理解和实现所公开实施例的其他变型。在权利要求书中,“包括”一词不排除其他元件或步骤,词语“一”或“一个”不排除多个。单个处理器或其他单元可以满足权利要求中记载的若干项目的功能。尽管在互不相同的从属权利要求中记载了特定元件,但是这并不指示不能有利地使用这些元件的组合。计算机程序可以存储和/或分布在适当的介质上,例如与其他硬件一起提供或作为其他硬件的一部分提供的光学存储介质或固态介质,但计算机程序也可以以其他形式分布,例如经由因特网或其他有线或无线的远程通信系统分布。权利要求书中的任何附图标记都不得被解释为对范围的限制。
Claims (13)
1.一种磁共振成像系统(100),其中,所述磁共振成像系统包括:
-存储器(134),其用于存储机器可执行指令(140)和脉冲序列命令(142),其中,所述脉冲序列命令被配置用于根据DCE磁共振成像协议来控制所述磁共振成像系统,
-用户界面(200);
-处理器(130),其用于控制所述磁共振成像系统,其中,所述机器可执行指令的执行使所述处理器:
-使用所述脉冲序列命令来控制(500)所述磁共振成像系统,以针对变化的翻转角两次或更多次采集校准磁共振数据(144);
-将所述校准磁共振数据的每个采集重建(502)为校准图像(146),以创建一组可变翻转角图像(148);
-使用所述一组可变翻转角图像来计算(504)T1映射(150);
-至少部分地使用所述T1映射来计算(506)针对预定磁共振成像造影剂的造影剂校准(152);
-使用校准准确度模型来计算(508)估计的校准误差(154),所述估计的校准误差描述所述造影剂校准和/或所述T1映射中的估计的误差,其中,所述校准准确度模型被配置用于使用所述一组可变翻转角图像来计算所述估计的校准误差;并且
-在所述估计的校准误差在预定校准误差范围之外的情况下在所述用户界面上显示(510)校准警告消息(202),并且其中,所述机器可执行指令的执行还使所述处理器:
-接收(600)一组扫描参数(156),所述一组扫描参数用于在所述脉冲序列命令的执行期间控制磁共振成像;
-使用所述一组扫描参数、所述估计的校准误差和用于基于梯度回波的磁共振成像的稳态信号方程来计算(602)预测的DCE图像质量(158);
-在所述预测的DCE图像质量在预定图像质量范围之外的情况下在所述用户界面上显示(604)采集警告消息(300)。
2.如权利要求1所述的磁共振成像系统,其中,所述机器可执行指令的执行还使所述处理器:
-利用所述脉冲序列命令来控制所述磁共振成像系统,以采集DCE磁共振数据的序列;并且
-使用所述DCE磁共振数据和所述造影剂校准将所述DCE磁共振数据的序列重建为一系列DCE磁共振图像(806)。
3.如权利要求1或2所述的磁共振成像系统,其中,所述机器可执行指令的执行还使所述处理器:
-将所述扫描参数存储(700)在所述存储器中的日志文件(160)中;
-重复(702)读取所述日志文件以获得一组当前扫描参数;
-将所述当前扫描参数与一组预定扫描参数(162)进行比较(704),以确定所述当前扫描参数是否在预定操作范围之外;并且
-在所述当前扫描参数在所述预定操作范围之外的情况下显示(706)所述用户界面的操作参数警告(400)。
4.如权利要求3所述的磁共振成像系统,其中,所述预定扫描参数是根据以下中的任一项确定的:公布的标准、医师推荐或命令、患者简档、患者数据库以及其组合。
5.如权利要求3或4所述的磁共振成像系统,其中,对所述日志文件的所述重复读取和所述当前扫描参数的所述比较是作为背景过程由所述处理器执行的。
6.如权利要求3、4或5所述的磁共振成像系统,其中,所述机器可执行指令的执行还使所述处理器响应于以下中的任一项而在所述用户界面上显示一组建议的扫描参数改变:显示所述操作者参数警告、所述采集警告消息、所述校准警告消息以及其组合。
7.如权利要求6所述的磁共振成像系统,其中,所述机器可执行指令的所述执行还使所述处理器至少部分地使用以下中的任一项来确定所述一组建议的扫描参数改变:查找表、先前扫描数据库、扫描参数之间的相互关系的模型以及其组合。
8.如前述权利要求中的任一项所述的磁共振成像系统,其中,所述校准警告包括以下中的任一项:所述估计的校准误差在所述预定校准误差范围之外的警告、建议的校正动作以及其组合。
9.如前述权利要求中的任一项所述的磁共振成像系统,其中,所述校准准确度模型被配置用于使用以下中的一项来计算所述估计的校准误差:基于直方图的噪声估计模型、基于小波变换的噪声估计模型、或者通过使用根据所述一组可变翻转角图像的不同组合的相减构建的数据集来估计噪声的噪声估计模型。
10.如前述权利要求中的任一项所述的磁共振成像系统,其中,所述机器可执行指令的执行还使所述处理器:
-使用所述脉冲序列命令来控制所述磁共振成像系统,以采集B1映射磁共振数据;并且
-使用所述B1映射磁共振数据来计算B1图,其中,还至少部分地使用所述B1图来计算所述造影剂校准。
11.一种包括用于由控制磁共振成像系统(100)的处理器执行的机器可执行指令(140)的计算机程序产品,其中,所述磁共振成像系统包括用户界面(200),其中,所述机器可执行指令的执行使所述处理器:
-使用脉冲序列命令(142)来控制(500)所述磁共振成像系统,以针对变化的翻转角两次或更多次采集校准磁共振数据(144),其中,所述脉冲序列命令被配置用于根据DCE磁共振成像协议来控制所述磁共振成像系统;
-将所述校准磁共振数据的每个采集重建(502)为校准图像(146),以创建一组可变翻转角图像(148);
-使用所述一组可变翻转角图像来计算(504)T1映射(150);
-至少部分地使用所述T1映射来计算(506)针对预定磁共振成像造影剂的造影剂校准(152);
-使用校准准确度模型来计算(508)估计的校准误差(154),所述估计的校准误差描述所述造影剂校准和/或所述T1映射中的估计的误差,其中,所述校准准确度模型被配置用于使用所述一组可变翻转角图像来计算所述估计的校准误差;并且
-在所述估计的校准误差在预定校准误差范围之外的情况下在所述用户界面上显示(510)校准警告消息(202),并且
其中,所述机器可执行指令的执行还使所述处理器:
-接收(600)一组扫描参数(156),所述一组扫描参数用于在所述脉冲序列命令的执行期间控制磁共振成像;
-使用所述一组扫描参数、所述估计的校准误差和用于基于梯度回波的磁共振成像的稳态信号方程来计算(602)预测的DCE图像质量(158);
-在所述预测的DCE图像质量在预定图像质量范围之外的情况下在所述用户界面上显示(604)采集警告消息(300)。
12.一种操作磁共振成像系统(100)的方法,其中,所述磁共振成像系统包括用户界面,其中,所述方法包括:
-使用脉冲序列命令(142)来控制(500)所述磁共振成像系统,以针对变化的翻转角两次或更多次采集校准磁共振数据,其中,所述脉冲序列命令被配置用于根据DCE磁共振成像协议来控制所述磁共振成像系统;
-将所述校准磁共振数据的每个采集重建(502)为校准图像(146),以创建一组可变翻转角图像(148);
-使用所述一组可变翻转角图像来计算(504)T1映射(150);
-至少部分地使用所述T1映射来计算(506)针对预定磁共振成像造影剂的造影剂校准(152);
-使用校准准确度模型来计算(508)估计的校准误差(154),所述估计的校准误差描述所述造影剂校准和/或所述T1映射中的估计的误差,其中,所述校准准确度模型被配置用于使用所述一组可变翻转角图像来计算所述估计的校准误差;并且
-在所述估计的校准误差在预定校准误差范围之外的情况下在用户界面(200)上显示(510)校准警告消息(202),并且
其中,所述方法还包括:
-接收(600)一组扫描参数(156),所述一组扫描参数用于在所述脉冲序列命令的执行期间控制磁共振成像;
-使用所述一组扫描参数、所述估计的校准误差和用于基于梯度回波的磁共振成像的稳态信号方程来计算(602)预测的DCE图像质量(158);
-在所述预测的DCE图像质量在预定图像质量范围之外的情况下在所述用户界面上显示(604)采集警告消息(300)。
13.如权利要求12所述的方法,其中,所述方法还包括:
-向对象注射所述预定磁共振成像造影剂;
-利用所述脉冲序列命令来控制所述磁共振成像系统,以采集DCE磁共振数据的序列;并且
-使用所述DCE磁共振数据和所述造影剂校准将所述DCE磁共振数据的序列重建为一系列DCE磁共振图像。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
EP16200952.6 | 2016-11-28 | ||
EP16200952 | 2016-11-28 | ||
PCT/EP2017/080135 WO2018096006A1 (en) | 2016-11-28 | 2017-11-23 | Image quality control in dynamic contrast enhanced magnetic resonance imaging |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110226098A true CN110226098A (zh) | 2019-09-10 |
CN110226098B CN110226098B (zh) | 2021-11-23 |
Family
ID=57406174
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201780084811.6A Active CN110226098B (zh) | 2016-11-28 | 2017-11-23 | 动态对比增强的磁共振成像中的图像质量控制 |
Country Status (6)
Country | Link |
---|---|
US (1) | US10935617B2 (zh) |
EP (1) | EP3545327A1 (zh) |
JP (1) | JP7199352B2 (zh) |
CN (1) | CN110226098B (zh) |
RU (1) | RU2019119770A (zh) |
WO (1) | WO2018096006A1 (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111257811A (zh) * | 2020-04-29 | 2020-06-09 | 山东奥新医疗科技有限公司 | 磁共振成像方法、装置及系统 |
CN115965621A (zh) * | 2023-02-15 | 2023-04-14 | 中国医学科学院阜外医院 | 一种基于磁共振成像的主要心脏不良事件的预测装置 |
WO2023134116A1 (zh) * | 2022-01-13 | 2023-07-20 | 浙江大学 | 测量水分子跨细胞膜流出速率的方法及应用、胶质瘤磁共振影像学标志物的测量方法及系统 |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10775464B2 (en) * | 2017-01-31 | 2020-09-15 | Regents Of The University Of Minnesota | System and method for dynamic, cardiac phase-resolved quantitative longitudinal relaxation parameter mapping |
EP3584596A1 (en) * | 2018-06-19 | 2019-12-25 | Koninklijke Philips N.V. | Tracking of image quality in magnetic resonance imaging |
JP7301724B2 (ja) * | 2019-11-29 | 2023-07-03 | キヤノンメディカルシステムズ株式会社 | 医用画像処理装置、方法及びプログラム |
US11982728B2 (en) | 2022-02-16 | 2024-05-14 | Shanghai United Imaging Healthcare Co., Ltd. | Systems and methods for magnetic resonance imaging |
Citations (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2005342511A (ja) * | 2004-06-01 | 2005-12-15 | General Electric Co <Ge> | 自動プロトコル選択のための方法及び装置 |
CN101912262A (zh) * | 2010-07-22 | 2010-12-15 | 中国科学院深圳先进技术研究院 | 磁共振成像定量参数计算方法及系统 |
US20110026799A1 (en) * | 2008-03-27 | 2011-02-03 | Koninklijke Philips Electronics N.V. | Flip angle imaging with improved b1 mapping for multi-rf transmit systems |
CN102077108A (zh) * | 2008-04-28 | 2011-05-25 | 康奈尔大学 | 分子mri中的磁敏度精确量化 |
US20120128132A1 (en) * | 2009-05-14 | 2012-05-24 | University Health Network | Phantom for contrast imaging calibration |
US20120277833A1 (en) * | 2011-04-29 | 2012-11-01 | Medtronic, Inc. | Electrical stimulation therapy based on head position |
CN103282788A (zh) * | 2010-10-13 | 2013-09-04 | 皇家飞利浦电子股份有限公司 | 用于t1校准的多隔室mri体模 |
CN103400376A (zh) * | 2013-07-19 | 2013-11-20 | 南方医科大学 | 一种乳腺动态增强磁共振图像序列的配准方法 |
CN103424726A (zh) * | 2012-05-21 | 2013-12-04 | 西门子公司 | 连续校正磁共振测量序列的相位误差的方法和磁共振设备 |
CN104155623A (zh) * | 2013-04-17 | 2014-11-19 | 西门子公司 | 用于自动确定组织种类的磁场反转时间的方法和系统 |
CN104345289A (zh) * | 2013-08-08 | 2015-02-11 | 西门子公司 | 用于确定水的t1 时间和脂肪的t1 时间的方法和磁共振设备 |
CN104519794A (zh) * | 2012-08-06 | 2015-04-15 | 皇家飞利浦有限公司 | 基于动态对比增强成像的渗透率度量 |
US20150160320A1 (en) * | 2013-12-11 | 2015-06-11 | Siemens Aktiengesellschaft | Accuracy and Off-Resonance Performance in Quantitative Magnetization Prepared Magnetic Resonance Imaging |
US20150309141A1 (en) * | 2014-04-25 | 2015-10-29 | New York University | System, method, and computer-accessible medium for determining at least one characteristic of at least one tissue or at least one mri measurement condition of the at least one tissue using active contrast encoding magnetic resonance imaging procedure(s) |
US20160291107A1 (en) * | 2015-04-02 | 2016-10-06 | The General Hospital Corporation | 3D Balanced EPI Magnetic Resonance Fingerprinting |
Family Cites Families (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5600574A (en) * | 1994-05-13 | 1997-02-04 | Minnesota Mining And Manufacturing Company | Automated image quality control |
US20060025667A1 (en) | 2004-07-29 | 2006-02-02 | Edward Ashton | Method for tumor perfusion assessment in clinical trials using dynamic contrast enhanced MRI |
DE102005028873B4 (de) | 2005-06-22 | 2014-07-24 | Siemens Aktiengesellschaft | Verfahren zur Durchführung einer Untersuchung und Diagnosevorrichtung hierfür |
US7528365B2 (en) * | 2006-02-07 | 2009-05-05 | Applera Corporation | Chemical noise reduction for mass spectrometry |
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 |
US8705819B2 (en) | 2007-06-21 | 2014-04-22 | Koninklijke Philips N.V. | Adjusting acquisition protocols for dynamic medical imaging using dynamic models |
DE102007036564B4 (de) * | 2007-08-03 | 2009-07-09 | Siemens Ag | Verfahren zur gewebespezifischen Optimierung von Bildgebungsparametern |
DE102009038139A1 (de) | 2009-08-19 | 2011-02-24 | Siemens Aktiengesellschaft | Verfahren zur Steuerung einer Magnetresonanzanlage |
US8731736B2 (en) * | 2011-02-22 | 2014-05-20 | Honda Motor Co., Ltd. | System and method for reducing driving skill atrophy |
US9715575B2 (en) | 2011-03-29 | 2017-07-25 | Koninklijke Philips N.V. | Image acquisition and/or image related parameter recommender |
US9013182B2 (en) | 2011-12-16 | 2015-04-21 | Rajiv Gandhi Cancer Institute & Research Centre | Method for computing pharmacokinetic parameters in MRI |
JP5932519B2 (ja) | 2012-06-27 | 2016-06-08 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | スキャン条件調整装置および撮影装置並びにプログラム |
US8995738B2 (en) * | 2013-02-18 | 2015-03-31 | Wisconsin Alumni Research Foundation | System and method for magnetic resonance imaging parametric mapping using confidence maps |
US9339239B2 (en) * | 2013-09-10 | 2016-05-17 | Ohio State Innovation Foundation | Methods and devices for optimization of magnetic resonance imaging protocols |
DE102013220288B4 (de) | 2013-10-08 | 2016-09-01 | Siemens Healthcare Gmbh | Optimieren von Steuerbefehlen zum Aufnehmen von Magnetresonanz-Bilddaten |
EP3125765A4 (en) * | 2014-03-04 | 2018-04-18 | The Trustees of Columbia University in the City of New York | Regularization of images |
-
2017
- 2017-11-23 JP JP2019528079A patent/JP7199352B2/ja active Active
- 2017-11-23 US US16/463,922 patent/US10935617B2/en active Active
- 2017-11-23 EP EP17804867.4A patent/EP3545327A1/en active Pending
- 2017-11-23 CN CN201780084811.6A patent/CN110226098B/zh active Active
- 2017-11-23 WO PCT/EP2017/080135 patent/WO2018096006A1/en active Application Filing
- 2017-11-23 RU RU2019119770A patent/RU2019119770A/ru not_active Application Discontinuation
Patent Citations (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2005342511A (ja) * | 2004-06-01 | 2005-12-15 | General Electric Co <Ge> | 自動プロトコル選択のための方法及び装置 |
US20110026799A1 (en) * | 2008-03-27 | 2011-02-03 | Koninklijke Philips Electronics N.V. | Flip angle imaging with improved b1 mapping for multi-rf transmit systems |
CN102077108A (zh) * | 2008-04-28 | 2011-05-25 | 康奈尔大学 | 分子mri中的磁敏度精确量化 |
US20120128132A1 (en) * | 2009-05-14 | 2012-05-24 | University Health Network | Phantom for contrast imaging calibration |
CN101912262A (zh) * | 2010-07-22 | 2010-12-15 | 中国科学院深圳先进技术研究院 | 磁共振成像定量参数计算方法及系统 |
CN103282788A (zh) * | 2010-10-13 | 2013-09-04 | 皇家飞利浦电子股份有限公司 | 用于t1校准的多隔室mri体模 |
US20120277833A1 (en) * | 2011-04-29 | 2012-11-01 | Medtronic, Inc. | Electrical stimulation therapy based on head position |
CN103424726A (zh) * | 2012-05-21 | 2013-12-04 | 西门子公司 | 连续校正磁共振测量序列的相位误差的方法和磁共振设备 |
CN104519794A (zh) * | 2012-08-06 | 2015-04-15 | 皇家飞利浦有限公司 | 基于动态对比增强成像的渗透率度量 |
CN104155623A (zh) * | 2013-04-17 | 2014-11-19 | 西门子公司 | 用于自动确定组织种类的磁场反转时间的方法和系统 |
CN103400376A (zh) * | 2013-07-19 | 2013-11-20 | 南方医科大学 | 一种乳腺动态增强磁共振图像序列的配准方法 |
CN104345289A (zh) * | 2013-08-08 | 2015-02-11 | 西门子公司 | 用于确定水的t1 时间和脂肪的t1 时间的方法和磁共振设备 |
US20150160320A1 (en) * | 2013-12-11 | 2015-06-11 | Siemens Aktiengesellschaft | Accuracy and Off-Resonance Performance in Quantitative Magnetization Prepared Magnetic Resonance Imaging |
US20150309141A1 (en) * | 2014-04-25 | 2015-10-29 | New York University | System, method, and computer-accessible medium for determining at least one characteristic of at least one tissue or at least one mri measurement condition of the at least one tissue using active contrast encoding magnetic resonance imaging procedure(s) |
US20160291107A1 (en) * | 2015-04-02 | 2016-10-06 | The General Hospital Corporation | 3D Balanced EPI Magnetic Resonance Fingerprinting |
Non-Patent Citations (1)
Title |
---|
PING WANG: "Effects of flip angle uncertainty and noise on the accuracy of DCE-MRI metrics:compasion between standard concentration-based and signal difference methods", 《MAGNETIC RESONANCE IMAGING》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111257811A (zh) * | 2020-04-29 | 2020-06-09 | 山东奥新医疗科技有限公司 | 磁共振成像方法、装置及系统 |
WO2023134116A1 (zh) * | 2022-01-13 | 2023-07-20 | 浙江大学 | 测量水分子跨细胞膜流出速率的方法及应用、胶质瘤磁共振影像学标志物的测量方法及系统 |
US11988734B2 (en) | 2022-01-13 | 2024-05-21 | Zhejiang University | Method and application for measuring the intracellular water transmembrane efflux rate, and measurement method and system for magnetic resonance imaging biomarker of glioma |
CN115965621A (zh) * | 2023-02-15 | 2023-04-14 | 中国医学科学院阜外医院 | 一种基于磁共振成像的主要心脏不良事件的预测装置 |
Also Published As
Publication number | Publication date |
---|---|
JP2020503093A (ja) | 2020-01-30 |
RU2019119770A (ru) | 2020-12-28 |
US10935617B2 (en) | 2021-03-02 |
EP3545327A1 (en) | 2019-10-02 |
US20190285711A1 (en) | 2019-09-19 |
WO2018096006A1 (en) | 2018-05-31 |
CN110226098B (zh) | 2021-11-23 |
JP7199352B2 (ja) | 2023-01-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110226098A (zh) | 动态对比增强的磁共振成像中的图像质量控制 | |
US11249160B2 (en) | Image segmentation using reference gray scale values | |
EP3270178B1 (en) | A system and method for determining optimal operating parameters for medical imaging | |
CN110476075A (zh) | 针对解剖区域的磁共振指纹识别词典的选择 | |
JP2014064889A (ja) | 医療映像装置及びその制御方法 | |
WO2015140277A1 (en) | Control of magnetic resonance imaging acquisition using modeling | |
US20200041594A1 (en) | Sub voxel resolution magnetic resonance fingerprinting imaging | |
CN112292606A (zh) | 在磁共振成像中跟踪图像质量 | |
JP2021526403A (ja) | Cest磁気共鳴イメージングにおける動き検出 | |
JP6636676B1 (ja) | 強度補正された磁気共鳴画像 | |
JP7420834B2 (ja) | 運動アーチファクト検出の方法 | |
US20220301687A1 (en) | Uncertainty maps for deep learning eletrical properties tomography | |
EP4306983A1 (en) | Making anatomical measurements using magnetic resonance imaging | |
CN111033290B (zh) | 傅里叶空间中的磁共振指纹识别 | |
US20240156346A1 (en) | Staged reconstruction of planning images for cardiac magnetic resonance imaging | |
EP3936881A1 (en) | Magnetic resonance fingerprinting quality assurance | |
EP4220212A1 (en) | Simulated computer tomographic imaging | |
JP2016209336A (ja) | 磁気共鳴イメージング装置 |
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 |