CN101329388B - 借助磁共振设备自动确定灌注的方法和装置 - Google Patents

借助磁共振设备自动确定灌注的方法和装置 Download PDF

Info

Publication number
CN101329388B
CN101329388B CN200810110254.3A CN200810110254A CN101329388B CN 101329388 B CN101329388 B CN 101329388B CN 200810110254 A CN200810110254 A CN 200810110254A CN 101329388 B CN101329388 B CN 101329388B
Authority
CN
China
Prior art keywords
data
magnetic resonance
perfusion
regressor
biosome
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.)
Expired - Fee Related
Application number
CN200810110254.3A
Other languages
English (en)
Other versions
CN101329388A (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 AG
Original Assignee
Siemens 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 AG filed Critical Siemens AG
Publication of CN101329388A publication Critical patent/CN101329388A/zh
Application granted granted Critical
Publication of CN101329388B publication Critical patent/CN101329388B/zh
Expired - Fee Related 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/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/563—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution of moving material, e.g. flow contrast angiography
    • G01R33/56341—Diffusion imaging
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/02—Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
    • A61B5/026—Measuring blood flow
    • A61B5/0263—Measuring blood flow using NMR
    • 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/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/563—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution of moving material, e.g. flow contrast angiography
    • G01R33/56366—Perfusion 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/4806—Functional imaging of brain activation

Landscapes

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

Abstract

本发明涉及借助磁共振设备(5)自动确定灌注的方法和装置。在此用灌注灵敏的成像序列从生命体的体积单元(15)获得时间段上的多个第一MR数据(21)，用控制成像序列、尤其是利用灌注不灵敏的成像序列从同一个体积单元(15)获得该时间段上的多个第二MR数据(22)。对第一MR数据和第二MR数据进行统计分析，以确定该体积单元(15)中的灌注(μ)。

Description

借助磁共振设备自动确定灌注的方法和装置
技术领域
本发明涉及利用磁共振设备自动确定灌注的方法和装置。此外，本发明还涉及装备了本发明装置的磁共振设备以及计算机程序产品，利用该计算机程序产品可以执行本发明的方法，以及涉及电子可读数据载体，在该数据载体上存储了本发明的方法。 
背景技术
MR灌注方法(借助磁共振设备确定灌注的方法)例如用于测量不同身体部位如头部的血流(脑血流(CBF))。CBF是在每分钟100g组织里流动的动脉血的体积(mL)，对人来说大脑里的CBF的典型值是60mL/(100g×min)。如果大脑的密度接近1g/mL，则人的CBF是0.6mL/(mL×min)或者0.01s-1。也就是说如果把人看作一个体积，则CBF的大小就是时间的倒数，即定义动脉血对组织体积的供应的心率常数。CBF与组织中的血量(CBV，血体积)没有因果联系与体积内的血液运动(血流速度)也没有因果联系。在借助磁共振设备确定灌注时，根据现有技术采用ASL方法(arterial spin labelling，动脉旋转标签)，其中血液中包含的水被置于一种特别的磁状态(通常是反磁化)；简而言之：对血液进行了标注，以便将流入所观察的体积单位的血液微粒与该体积单位中的其它组织区分开来。在此一方面利用灌注灵敏的图像序列产生MR图像，另一方面产生具有控制成像序列的MR控制图像。灌注信息在此仅通过图像对比度的很小的变化形成，该图像对比度存在于流入感兴趣区域的、具有特殊磁状态的标记的血液微粒和在该感兴趣区域中被采集了MR图像的组织之间。典型的，灌注信号的数量级仅为相应MR图像的总强度的很小一部分。因此对相对灌注的图像或用于计算量化灌注的图像的获取很容易受到伪影的影响。因此目前必须产生多幅MR图像来获取很多分钟内的小灌注信号，然后从中产生在这段时间上的用灌注信息标记的一系列MR图像和控制图像(没有灌注信息)。在此交替地产生具有灌注信息的图像和对应的控制图像。 
为了从这些具有灌注信息的图像和对应的控制图像中产生其中可以看见灌注的MR图像，根据现有技术分别形成具有灌注信息的图像和该图像对应的控制图像之间的差。然后将由此获得的差值关于所述图像序列取平均值。最后确定比例因子或校准因子，以获得所产生的MR图像中的相对或量化灌注信息。 
发明内容
常规计算的、显示灌注信息的MR图像的质量很差，因为根据现有技术的用于确定灌注的方法是很容易出现错误的，例如很容易出现伪影或其它干扰。因此本发明要解决的技术问题是改善显示灌注信息的MR图像的质量。 
按照本发明，该技术问题通过一种用于自动确定灌注的方法，通过一种用于自动确定灌注的磁共振设备的装置，通过一种磁共振设备，通过一种计算机程序产品和一种电子可读数据载体解决。 
在本发明的范围中，提供一种借助磁共振设备自动确定灌注的方法。在此本发明的方法连续产生多个第一MR数据，该MR数据是用灌注灵敏的成像序列从生命体的体积单元获得的。按照相似的方式利用控制成像序列、尤其是利用灌注不灵敏的成像序列从同一个体积单元获得连续的多个第二MR数据。此外，本发明的方法还包括步骤：借助统计分析对第一MR数据和第二MR数据进行分析，使得由此确定该体积单元中的灌注。 
在此灌注灵敏的成像序列应理解为：受到灌注的微粒被置于第一磁状态，从而使在借助磁共振设备采集第一MR数据时当微粒流入该体积单元时该微粒与其它微粒、例如没有受到灌注的组织的微粒不同。通过控制成像序列，使受到灌注的微粒置于第二磁状态，该第二磁状态在借助磁共振设备采集第二MR数据时应尽可能好地与第一磁状态区分开来，并优选对灌注效应不灵敏。 
与现有技术不同，在本发明的方法中没有在第一MR数据和第二MR数据之间形成差来确定灌注，而是对全部第一MR数据和第二MR数据进行统计分析。根据本发明，通过不在第一MR数据和第二MR数据之间形成差，还可以具有优点地由第一MR数据和第二MR数据计算出相应的绝对值，该绝对值可以用于确定灌注。 
在本发明的方法中，尤其是从第一MR数据和第二MR数据出发产生MR图像序列，其中具体地说由这些MR图像内的相应的图像点或体素确定一个信号分布。优选地，然后通过统计分析检查该信号分布，在此该统计分析确定系 数，然后分析这些系数以获得关于灌注的信息。 
在本发明优选的实施方式中，对体积单元交替地采集第一MR数据，接着是第二MR数据，接着又是第一MR数据，依次下去。换句话说，借助磁共振设备交替地用灌注灵敏的成像序列产生MR图像和然后用控制成像序列产生MR图像。具有灌注的体素关于由磁共振设备采集的灌注灵敏的值的信号分布具有之字形，因为例如第一MR数据中的值具有高的测量值，而第二MR数据中的值相应地具有比较低的测量值。然后对该之字形的信号分布进行统计分析，以获得相应体素的灌注信息。 
该统计分析可以根据通用线性模型或根据student t检验来进行。替换的，还可以采用其它统计方法，如交叉相关。 
在此，可以利用通用线性模型获得相对灌注信息以及量化灌注信息，而student t检验尤其是用于获得关于第一MR数据组中的像素的强度与第二MR数据组中的相应像素的强度不同的结论。此外还可以获得关于第一和第二MR数据的质量的结论，然后有利地将它们用于例如由于质量较差而不将特定的第一和/或第二MR数据用于确定灌注信息，和/或指示操作人员采集其它的第一和第二MR数据，因为还没有足够多的第一和第二MR数据达到令人满意的质量。 
特别是对于student t检验，可以计算按照本发明产生的MR图像内的对比度的统计数值，该统计数值直接与涉及MR图像中的灌注信息的对比度-噪声比相关。对该数值或对比度-噪声比的实时显示使得操作人员可以依据采集第一和第二MR数据的时间来跟踪按照本发明产生的MR图像内的灌注信息的质量。 
为了通过统计分析改善灌注信息，可以向该统计分析添加至少一个模型函数或函数模型作为至少一个回归量。在此，该至少一个回归量相当于一组回归量中的一个或多个回归量，该回归量组包括以下回归量： 
-从由磁共振设备采集的测量结果中推导出来的回归量。例如，如果可以从该测量结果中推导出磁共振设备的系统误差，则可以通过相应的模型函数来通过统计分析考虑该误差。被确定灌注的生物体的运动也可以通过该测量结果获得，只要这是刚性体的运动，即该生物体在运动时没有出现变形。 
-由不属于磁共振设备的装置采集的回归量。如果影响第一和第二MR数据的采集的干扰也被该装置采集，则可以通过相应的模型函数由统计分析来考虑这些干扰。 
-从功能变化或生物体的变化推导出的回归量。在此，功能变化在此尤其 是生物体的生理状态的变化或者生物体的功能活动。功能活动的例子是生物体部位、如生物体的手指的周期性运动。在此，可以通过灌注的变化、通过脑部血流的变化或通过BOLD(Blood Oxygenation Level Dependent，血氧水平依赖)效应中的变化来采集功能活动。 
回归量在此对应于统计分析应理解为一种解释参数，该参数具有对要解释的参数、在本发明中尤其是对灌注的解释的影响。在此有不期望的或引起干扰的回归量，如计划之外的生物体运动，但还有期望或自然的回归量，例如生物体的计划内的功能活动。虽然在确定灌注时可以由统计分析考虑不期望的以及期望的回归量，但是例如在通用线性模型中可以更好或更精确地计划期望的回归量，因为事先就知道存在相应的回归量以及该回归量何时、如何出现。 
该模型函数与用于做出关于灌注的结论的灌注模型函数一起，形成按照通用线性模型的多维统计分析的输入。 
此外，按照本发明所述回归量组可以包括以下回归量： 
-生物体的运动。在此尤其是理解为计划外的运动。 
-磁共振设备的扫描稳定性。例如，如果磁共振设备采集第一和第二MR数据时所具有的稳定性不恒定，则第一和第二MR数据的采集就会受到一定的波动。 
-生物体的呼吸。生物体的至少一部分根据呼吸运动。 
-生物体的心跳。至少生物体的心脏依据心跳进行运动，并且例如血液的流速也取决于心跳。 
本发明的方法不仅可以产生灌注信息，而且还可以获得以下结果： 
-关于通过本发明为相应的体积单元获得的结果的可靠性的说明(Angaben)。例如借助student t检验可以作出质量说明，即本发明确定的量化灌注信息有多好。 
-关于体积单元内的一个特定像素的对比度-噪声比的说明。由此可以有利地做出关于本发明方法确定的灌注信息的品质的结论。 
-关于以下伪影的说明： 
-生物体的呼吸 
-生物体的心跳 
-生物体的运动 
-生物体的功能活动 
-BOLD效应。 
在此，关于伪影的说明例如包括关于伪影的范围或大小的结论，这例如又可以推导出特定灌注信息的质量。如果例如通过统计分析获得生物体的呼吸和/或运动在采集第一和第二MR数据时过度剧烈，从而由此可能影响灌注信息的品质，则该信息对于分析灌注信息是有价值的。 
根据本发明还可以产生分开地对每个回归量显示不同回归量的强度和/或统计重要性的映像(Abbildung)。通过这些映像例如可以推导出特定MR图像内的哪些体素会受到呼吸、心跳、运动等的影响。 
借助上面描述的结果，可以对采集的第一和第二MR数据进行质量控制。为此可以依据该结果在该方法中排除或不分析已经产生的第一和/或第二MR数据。此外如果结果接近，还可以产生其它第一和/或第二MR数据。这例如可能是统计分析获得一个或多个伪影的范围或大小超过预定阈值的情况。所采集的第一和第二MR数据的质量例如通过生物体的剧烈运动而受到强烈影响，其中该运动本身通过对MR数据的分析来识别。也可能从上面描述的结果中推导出错误的积极的灌注信息或错误的负面的灌注信息，并用于质量控制。 
通过本发明的方法可以借助统计分析产生包含以下信息的MR图像： 
-所观察的体积单元的组织的血氧变化。由于BOLD效应对磁共振设备采集的结果产生影响，因此可以从第一和第二MR数据中推导出血氧变化。 
-生命体的功能活动。由于生命体的功能活动同样具有对磁共振设备获取的结果的影响，因此还可以从第一和第二MR数据中推导出功能活动的范围。 
-BOLD效应和特定功能活动之间的相关的结果。由于统计分析已经用于获得灌注信息，因此有利的是无需大的花费就可以借助统计分析来确定BOLD效应和特定功能活动之间的相关。 
根据本发明，可以始终实时地产生含有灌注信息并由第一和第二MR数据产生的MR图像，从而MR图像可以始终通过根据本发明新采集的第一和第二MR数据来更新。换句话说，显示灌注信息的第一序列MR图像尤其是在统计分析事先已经确定在该第一序列中显示的灌注信息具有相应品质时产生。然后该第一序列基于已经获得的第一和第二MR数据和新采集的其它第一和第二MR数据而始终是最新的。 
当然还可以在结束对所有第一和第二MR数据的采集之后才分析第一和第二MR数据。 
本发明方法产生的、具有灌注信息的MR图像，可以基于根据本发明所采集的全部第一和第二MR数据产生，或者还可以基于第一和第二MR数据的特定子集产生。如果MR图像由第一和第二MR数据的子集产生，则该子集尤其是不包含根据统计分析不能用于获得灌注信息的第一和第二MR数据。此外，该子集包括尽可能最新的第一和第二MR数据，从而例如从该子集中周期性地去掉最老的第一和第二MR数据。 
如果用通用线性模型进行统计分析，则可以调整该通用线性模型的系数，使得利用该系数可以对相对灌注做出说明和/或对绝对灌注做出说明。 
在本发明范围内，还提供用于自动确定灌注的磁共振设备的装置。该装置包括用于控制磁共振设备的控制单元，用于接收多个第一和第二MR数据的接收装置，该第一和第二MR数据由磁共振设备(例如本地线圈)记录，以及分析该第一和第二MR数据并从中产生MR图像的分析装置。所述装置实施为通过控制装置控制磁共振设备，使得磁共振设备记录或采集第一和第二MR数据。控制装置控制磁共振设备，使得该磁共振设备在接收第一MR数据时实施涉及体积单元的灌注灵敏的成像序列，在接收第二MR数据时实施涉及体积单元的控制成像序列，尤其是灌注不灵敏的成像序列。所述装置可以借助其分析装置对第一和第二MR数据进行统计分析，以由此确定该体积单元中的灌注。 
本发明装置的优点基本上对应于本发明方法的优点，因此不再重复。 
在本发明的范围中，还提供一种包括上述装置的磁共振设备。 
此外本发明描述一种计算机程序产品，尤其是软件，可以加载到磁共振设备的可编程控制装置的存储器中。利用程序装置和该计算机程序产品可以执行本发明方法的所有上述实施方式，如果该计算机程序产品在装置中运行的话。 
最后本发明公开了一种电子可读数据载体，例如DVD，在该DVD上存储了电子可读控制信息、尤其是软件。如果从数据载体读取该控制信息并存储到磁共振设备的控制装置中，则可以执行上述方法的本发明所有实施方式。 
利用本发明，相对于现有技术大大改善了含有灌注信息的MR图像的获取和评估。除了灌注信息之外，按照本发明，还可以确定所获得的结果的精度信息(例如量化的灌注信息)、关于对比度-噪声比的信息和关于特定伪影的范围的信息。这些信息可以用图像的形式表示。 
本发明尤其适用于借助磁共振设备确定生物体身体中液体的灌注或流动，以识别MR图像中的灌注。当然本发明不限于优选的应用领域，还可以用于获 得其它信息或结果，如组织的血氧变化、功能活动的范围或MR图像中显示的信息的品质。一般来说，本发明的方法可用于确定来自包含信息的数据组以及来自检查数据组的信息，该检查数据组主要通过不包含该信息而与该数据组区分开。 
附图说明
下面借助本发明的实施方式参照附图解释本发明。 
图1示出具有本发明装置的磁共振设备。 
图2A示出第一和第二MR数据或MR图像的时间序列， 
图2B示出第一和第二MR数据的相应通用线性模型。 
具体实施方式
图1示出磁共振设备5的实施例，利用该磁共振设备可以自动确定灌注。该磁共振设备5的核心是断层造影装置3，患者O在该断层造影装置中位于环形基本磁场磁铁(未示出)中的患者卧榻2上，该基本磁场磁铁包围测量空间4。 
患者卧榻2可在纵向、即沿着断层造影装置3的纵轴移动。在断层造影装置3中，在基本磁场磁铁内具有整体线圈(未示出)，利用该整体线圈发射高频脉冲以及接收该高频脉冲。此外，断层造影装置3具有梯度线圈(未示出)，以便能够在每个空间方向上施加磁场梯度。 
断层造影装置3由控制装置6控制，该控制装置6在此与断层造影装置3分开示出。控制装置6与终端7连接，该终端7具有显示屏8、键盘9和鼠标10。终端7尤其用作用户接口，操作人员通过该用户接口操作控制装置6并由此操作断层造影装置3。控制装置6和终端7都是磁共振设备5的部件。 
此外在图1中示出DVD 14，在该DVD上存储软件，当该软件加载到控制装置6中之后可以执行本发明的方法。 
此外，磁共振设备5可以具有所有其它常见的部件或特征，如用于连接通信网络的接口、用于连接图像信息系统的接口等等。但是为了清楚起见未在图1中示出所有这些部件。 
通过终端7，操作人员可以与控制装置6通信，从而负责执行期望的测量，其方法是例如由控制装置6控制断层造影装置3，使得通过天线发射需要的高 频脉冲序列并以合适的方式接通梯度线圈。通过控制装置6还获取来自断层造影装置3的第一MR数据21和第二MR数据22，并在分析单元13中转换为相应的图像(MR图像)，该分析单元是控制装置6的模块。这些图像然后例如在显示屏8上显示和/或存储在存储器中，或通过网络发送。 
患者卧榻2可以借助控制装置6电动地在断层造影装置3内移动。控制装置6具有控制单元11，其自动驱动患者卧榻2穿过断层造影装置3或驶过断层造影装置3内的不同位置。此外，控制装置11负责给出定义的磁场梯度，同时由整体线圈发射高频信号，该高频信号基本上对应于磁共振频率。替换的，该高频信号还可以用特别实施的局部线圈(发射/接收线圈)发射。 
然后借助控制装置6的接收信道12或测量装置，利用局部线圈1来确定和记录患者O体内相应体积单元15的第一MR数据21和第二MR数据22。然后在分析装置13中，由这些第一MR数据或MR图像21和第二MR数据或MR图像22产生其中可视灌注的MR图像。 
在图2A中示出第一MR数据21和第二MR数据22的时间序列。在此，已经借助ASL(动脉旋转标记)在拍摄第一MR图像21时标记了生命体O头部内流入被观察体素19中的动脉血，使得流入体素19中的血液与该体积单元内的组织可以区分开。相反，在拍摄第二MR图像22时不进行这样的对动脉血的标记。如图2A所示，在拍摄一张第一MR图像21之后拍摄一张第二MR图像22，而在每张第二MR图像22之后又拍摄一张第一MR图像21。换句话说在该时间段上拍摄的第一和第二MR图像21、22的序列交替地由第一MR图像21和第二MR图像22组成。 
用附图标记19标记的白色矩形对应于其中当前确定了灌注的体素。要指出图1中显示的体积单元15在患者O的腿内，而图2A中的体素19设置在患者的头部内。 
图2B中示出通用线性模型。在此，左侧向量Y的各个相应分量值(y1，y2，y3等)对应于相应多的在时间上连续的第一和第二MR图像21、22的关于体素19的各个测量值。换句话说，向量Y对应于连续的MR图像21、22的图像点或体素19在一段时间上的信号变化。 
直接在等号右侧的矩阵一方面包含灌注模型18，另一方面包含3个函数模型17。灌注模型18在此是一个具有值1、0、1、0等的向量，从而在图2B中灰色的划线对应于1，白色的划线对应于0。在此，值1意味着Y向量的相应 分量具有灌注信息，而值0则说明Y向量的相应分量不包含灌注信息。 
3个函数模型17是用于反映接收装置12的扫描特性的稳定性或不稳定性的函数模型。同样，在此可以采用用于反映患者O心跳的函数模型，和用于反映功能活动如患者O手指的周期性运动的函数模型。在此，如果存在不稳定性或心脏恰好在跳动或手指在运动，则表示相应函数模型的向量17通常作为分量值分别具有1，如果不是，则作为分量值具有0。例如如果BOLD效应在采集第一和第二MR图像21、22期间的一段时间内具有较高的值，则表示BOLD效应的函数模型的相应向量(未示出)在这段时间内对相应的第一MR图像21和相应的第二MR图像22都具有作为分量值的1。 
由灌注模型18和3个函数模型17组成的矩阵也称为设计矩阵。 
要指出，在该通用线性模型中还可以采用分量值-1来代替分量值0，从而使设计矩阵17、18的向量具有分量值-1和+1。 
与设计矩阵17、18相乘的向量包含由利用通用线性模型16的统计分析确定的系数μ、τ1、τ2、τ3或量化参数。在此，参数μ给出体素19中灌注的量化值。以同样的方式，参数τ1给出接收装置12的扫描特性的稳定性的量化值，参数τ2 给出患者O心跳的量化值，参数τ3给出患者O的关于体素19的功能活动的量化值。 
向量U(u1，u2，u3，...)包括由于噪声引起或者由没有被函数模型17反映的误差形成的剩余误差。 
在此确定灌注或参数μ以及参数τ1、τ2、τ3是借助通过与设计矩阵17、18的实施相匹配的通用线性模型，基于不同体素的相应测量值进行的，从而最后可以确定并显示关于大的体积段中的灌注的信息(以及关于扫描特性的稳定性、心跳和功能活动的信息)。可以这样调整参数μ，使得该参数代表相对或绝对的灌注值(单位：mL/(100g×min))。 
在相应选择函数模型17时，还可以利用关于相应系数τ1、τ2、τ3的通用线性模型获得代表每个体素的确定的结果(例如灌注)的可靠性或质量的数值的概率和其他统计数值。 

Claims (27)

1.一种借助磁共振设备(5)自动确定灌注的方法，其中该方法包括以下步骤：
产生在时间段上的多个第一磁共振数据(21)，该第一磁共振数据是用灌注灵敏的成像序列从体积单元(15)采集的，
产生在时间段上的多个第二磁共振数据(22)，该第二磁共振数据是利用控制成像序列从同一个体积单元(15)采集的，
其特征在于，该方法还包括步骤：对该第一磁共振数据(21)和第二磁共振数据(22)进行统计分析(16)，以确定该体积单元(15)中的灌注(μ)，
其中，向所述统计分析(16)添加作为至少一个回归量的至少一个模型函数(17)，其中该至少一个回归量包括一组回归量中的一个或多个回归量，该组回归量包括以下回归量：从由磁共振设备(5)采集的测量结果中导出的回归量，由不属于磁共振设备(5)的装置采集的回归量，以及从被执行灌注测量的生物体(O)的功能变化推导出的回归量。
2.根据权利要求1所述的方法，其特征在于，所述控制成像序列是灌注不灵敏的成像序列。
3.根据权利要求2所述的方法，其特征在于，灌注灵敏的成像序列借助灌注灵敏的磁场梯度设置产生，灌注不灵敏的成像序列借助灌注不灵敏的磁场梯度设置产生。
4.根据权利要求1所述的方法，其特征在于，从第一磁共振数据(21)和第二磁共振数据(22)中产生磁共振图像序列，其中确定各个图像点(19)的在所述时间段上的信号分布。
5.根据权利要求4所述的方法，其特征在于，通过统计分析(16)来检查该信号分布，其中确定系数(μ)，并通过分析该系数(μ)获得关于灌注的信息。
6.根据权利要求1所述的方法，其特征在于，对所述体积单元(15)交替地采集第一磁共振数据(21)的一部分和第二磁共振数据(22)的一部分。
7.根据权利要求1所述的方法，其特征在于，根据通用线性模型(16)或根据student t检验或根据交叉相关来进行所述统计分析。
8.根据权利要求1至7中任一项所述的方法，其特征在于，向所述统计分析(16)添加作为至少一个回归量的至少一个模型函数(17)，其中该至少一个回归量包括一组回归量中的一个或多个回归量，该组回归量包括以下回归量：被执行灌注测量的生物体(O)的运动，磁共振设备(5)的扫描稳定性，生物体(O)的呼吸，生物体(O)的心跳，生物体(O)的生理状态以及生物体(O)的功能活动。
9.根据权利要求1至7中任一项所述的方法，其特征在于，通过所述统计分析(16)确定来自一个组的至少一个其它结果，该组包括：
关于通过该方法为体积单元(15)确定的结果的可靠性的说明，
关于该体积单元(15)的像素(19)的对比度-噪声比的说明，
关于伪影的说明，
其中，该伪影包括被执行灌注测量的生物体(O)的呼吸，生物体(O)的心跳，生物体(O)的运动，生物体(O)的功能活动和BOLD效应。
10.根据权利要求9所述的方法，其特征在于，借助所述至少一个其它结果对产生的第一磁共振数据(21)和第二磁共振数据(22)进行质量控制，其中如果通过统计分析(16)获得的伪影超过预定阈值，则在该方法中排除已经产生的第一磁共振数据(21)和第二磁共振数据(22)和/或产生其它第一磁共振数据(21)和第二磁共振数据(22)。
11.根据权利要求1至7中任一项所述的方法，其特征在于，借助统计分析(16)产生体积单元(15)的磁共振图像，使得该磁共振图像包含一组信息的至少一个信息，该组信息包括：
该体积单元(15)的组织的血氧含量，
被执行灌注测量的生命体(O)的功能活动，
BOLD效应与该功能活动之间相关的结果。
12.根据权利要求1至7中任一项所述的方法，其特征在于，始终实时地产生体积单元(15)的磁共振图像，并且通过新的第一磁共振数据(21)和第二磁共振数据(22)来更新该磁共振图像。
13.根据权利要求1至7中任一项所述的方法，其特征在于，体积单元(15)的磁共振图像基于全部第一磁共振数据(21)和第二磁共振数据(22)产生，或者基于第一磁共振数据(21)和第二磁共振数据(22)的子集产生。
14.根据权利要求1至7中任一项所述的方法，其特征在于，用通用线性模型(16)来进行所述统计分析(16)，
调整该通用线性模型的系数(μ)，使得该系数(μ)给出体积单元(15)中的相对灌注单位或绝对灌注单位。
15.一种用于自动确定灌注的磁共振设备的装置，其中，
该装置(6)包括用于控制磁共振设备(5)的控制单元(11)，用于接收由磁共振设备(5)拍摄的多个体积单元(15)的第一磁共振数据(21)和第二磁共振数据(22)的接收装置(12)，以及用于分析该第一磁共振数据(21)和第二磁共振数据(22)并从第一磁共振数据(21)和第二磁共振数据(22)中产生该体积单元(15)的磁共振图像的分析装置(13)，在该磁共振图像中可以看见体积单元(15)中的灌注(μ)，
所述装置(6)实施为，该装置(6)通过控制单元(11)控制磁共振设备(5)，使得磁共振设备(5)在接收第一磁共振数据(21)时在体积单元(15)中实施灌注灵敏的成像序列，在接收第二磁共振数据(22)时在体积单元(15)中实施控制成像序列，
其特征在于，所述分析装置(13)实施为，该分析装置(13)对第一磁共振数据(21)和第二磁共振数据(22)进行统计分析(16)，以确定该体积单元(15)中的灌注(μ)。
16.根据权利要求15所述的装置，其特征在于，所述装置(6)实施为，该装置(6)通过控制单元(11)控制磁共振设备(5)，使得磁共振设备(5)在产生第二磁共振数据(22)时在体积单元(15)中实施灌注不灵敏的成像序列。
17.根据权利要求15所述的装置，其特征在于，所述装置(6)实施为，这样控制控制单元(11)，使得接收装置(12)交替地接收第一磁共振数据(21)的一部分和第二磁共振数据(22)的一部分。
18.根据权利要求15所述的装置，其特征在于，所述分析装置(13)实施为，该分析装置(13)根据通用线性模型(16)或根据student t检验或根据交叉相关来执行统计分析。
19.根据权利要求15所述的装置，其特征在于，所述分析装置(13)实施为，该分析装置(13)在所述统计分析中添加至少一个模型函数(17)作为至少一个回归量，
其中，该至少一个回归量包括一组回归量中的一个或多个回归量，该组回归量包括以下回归量：从由磁共振设备(5)采集的测量结果中导出的回归量，由不属于磁共振设备(5)的装置采集的回归量，以及从被执行灌注测量的生物体(O)的功能变化推导出的回归量。
20.根据权利要求15至19中任一项所述的装置，其特征在于，所述分析装置(13)实施为，该分析装置(13)在所述统计分析中添加至少一个模型函数(17)作为至少一个回归量，
其中，该至少一个回归量包括一组回归量中的一个或多个回归量，该组回归量包括以下回归量：被执行灌注测量的生物体(O)的运动，磁共振设备(5)的扫描稳定性，生物体(O)的呼吸，生物体(O)的心跳，生物体(O)的生理状态以及生物体(O)的功能活动。
21.根据权利要求15至19中任一项所述的装置，其特征在于，所述分析装置(13)实施为，该分析装置(13)通过统计分析(16)获得来自一个组的至少一个其它结果，该组包括：
关于通过该方法为体积单元(15)获得的结果的可靠性的说明，
关于体积单元(15)的像素(19)的对比度-噪声比的说明，
关于伪影的说明，
其中，该伪影包括被执行灌注测量的生物体(O)的呼吸，生物体(O)的心跳，生物体(O)的运动，生物体(O)的功能活动和BOLD效应。
22.根据权利要求21所述的装置，其特征在于，所述装置(6)实施为，该装置(6)借助所述至少一个其它结果对产生的第一磁共振数据(21)和第二磁共振数据(22)进行质量控制，其中如果分析装置(13)确定通过统计分析(16)获得的伪影超过预定阈值，则分析装置(13)在产生磁共振图像时排除已经产生的第一磁共振数据(21)和第二磁共振数据(22)和/或该装置(6)产生其它的第一磁共振数据(21)和第二磁共振数据(22)。
23.根据权利要求15至19中任一项所述的装置，其特征在于，所述分析装置(13)实施为，该分析装置(13)借助统计分析(16)产生磁共振图像，使得该磁共振图像包含一组信息的至少一个信息，该组信息包括：
该体积单元(15)的组织的血氧含量，
被执行灌注测量的生命体(O)的功能活动，
BOLD效应与该功能活动之间相关的结果。
24.根据权利要求15至19中任一项所述的装置，其特征在于，所述装置(6)实施为，所述分析装置(13)始终实时地产生磁共振图像，并通过新的第一磁共振数据(21)和第二磁共振数据(22)更新该磁共振图像。
25.根据权利要求15至19中任一项所述的装置，其特征在于，所述分析装置(13)实施为，该分析装置(13)基于全部第一磁共振数据(21)和第二磁共振数据(22)产生磁共振图像，或者基于第一磁共振数据(21)和第二磁共振数据(22)的子集产生磁共振图像。
26.根据权利要求15至19中任一项所述的装置，其特征在于，所述分析装置(13)实施为，该分析装置(13)用通用线性模型(16)执行统计分析，以及
该分析装置(13)调整该通用线性模型的系数(μ)，使得该系数(μ)给出体积单元(15)中的相对灌注单位或绝对灌注单位。
27.一种包括根据权利要求15-26中任一项所述的装置(6)的磁共振设备。
CN200810110254.3A 2007-06-22 2008-06-23 借助磁共振设备自动确定灌注的方法和装置 Expired - Fee Related CN101329388B (zh)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
DE102007028901.6 2007-06-22
DE102007028901A DE102007028901B4 (de) 2007-06-22 2007-06-22 Verfahren und Vorrichtung zur automatischen Bestimmung von Perfusion mittels einer Magnetresonanzanlage

Publications (2)

Publication Number Publication Date
CN101329388A CN101329388A (zh) 2008-12-24
CN101329388B true CN101329388B (zh) 2013-09-11

Family

ID=40075837

Family Applications (1)

Application Number Title Priority Date Filing Date
CN200810110254.3A Expired - Fee Related CN101329388B (zh) 2007-06-22 2008-06-23 借助磁共振设备自动确定灌注的方法和装置

Country Status (4)

Country Link
US (1) US8233961B2 (zh)
JP (1) JP2009000525A (zh)
CN (1) CN101329388B (zh)
DE (1) DE102007028901B4 (zh)

Families Citing this family (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE102009015516A1 (de) * 2009-04-02 2010-10-07 Mediri Gmbh Verfahren und Vorrichtung zur Ermittlung von T2 und T2* mittels Magnet- Resonanz-Tomographie
US8532353B2 (en) * 2011-11-23 2013-09-10 Vital Images, Inc. Synthetic visualization and quantification of perfusion-related tissue viability
JP5820549B2 (ja) * 2012-04-17 2015-11-24 コーニンクレッカ フィリップス エヌ ヴェＫｏｎｉｎｋｌｉｊｋｅ Ｐｈｉｌｉｐｓ Ｎ．Ｖ． 灌流イメージング
JP6645831B2 (ja) * 2012-10-31 2020-02-14 コーニンクレッカ フィリップス エヌ ヴェＫｏｎｉｎｋｌｉｊｋｅ Ｐｈｉｌｉｐｓ Ｎ．Ｖ． かん流撮像
DE102013224264B4 (de) * 2013-11-27 2018-12-27 Siemens Healthcare Gmbh Verfahren zur Verarbeitung von Magnetresonanz-Diffusionsbilddaten
DE102014205789A1 (de) * 2014-03-27 2015-10-01 Fraunhofer-Gesellschaft zur Förderung der angewandten Forschung e.V. Magnetresonanzsystem
DE102015109678A1 (de) * 2015-06-17 2016-12-22 Julius-Maximilians-Universität Würzburg Dynamische Fourier-Raum-Magnetresonanz
US10478557B2 (en) 2015-08-21 2019-11-19 Medtronic Minimed, Inc. Personalized parameter modeling methods and related devices and systems
CN109669151B (zh) 2019-02-12 2024-06-14 泰山医学院 磁共振动脉自旋标记灌注成像的质量控制体模及评价方法
DE102019214359A1 (de) * 2019-09-20 2021-03-25 Siemens Healthcare Gmbh Verfahren zu einer adaptiven Ansteuerung eines Magnetresonanzgerätes
US11295448B1 (en) * 2021-07-28 2022-04-05 Ischemaview, Inc. Concurrent display of hemodynamic parameters and damaged brain tissue

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1285044A (zh) * 1998-01-05 2001-02-21 耐克麦德英梅金公司 磁共振研究方法
CN1494871A (zh) * 2002-05-13 2004-05-12 ��ʽ���綫֥ 磁共振成像仪、流量量化仪及流量量化方法

Family Cites Families (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6980845B1 (en) * 1998-04-13 2005-12-27 The Trustees Of The University Of Pennsylvania Multi-slice cerebral blood flow imaging with continuous arterial spin labeling MRI
US6564080B1 (en) * 1999-03-31 2003-05-13 Kabushiki Kaisha Toshiba MR imaging on ASL technique
US6717405B2 (en) * 2002-04-12 2004-04-06 Beth Israel Deaconess Medical Center, Inc. Arterial spin labeling using time varying gradients
US7587233B2 (en) * 2002-05-13 2009-09-08 The Regents Of The University Of California Velocity-selective arterial spin labeling without spatial selectivity
US7243064B2 (en) * 2002-11-14 2007-07-10 Verizon Business Global Llc Signal processing of multi-channel data
JP4342809B2 (ja) * 2003-02-14 2009-10-14 株式会社東芝 Ｍｒｉ装置
US7328054B2 (en) * 2003-04-09 2008-02-05 The Mcw Research Foundation, Inc. Perfusion magnetic resonance imaging using encoded RF tagging pulses
WO2004114213A2 (en) * 2003-06-13 2004-12-29 The General Hospital Corporation Noise and cancellation for time domain reconstruction
DE102004043809B4 (de) * 2004-09-08 2008-01-24 Charité-Universitätsmedizin Berlin Verfahren zum Kalibrieren einer kontrastmittelgestützten Perfusionsbildgebung
US8131034B2 (en) * 2005-05-20 2012-03-06 Advanced Mri Technologies, Llc Efficient visualization of cortical vascular territories and functional MRI by cycled arterial spin labeling MRI
EP1978868A4 (en) * 2005-09-19 2011-05-11 Advanced Mri Technologies Llc EXTREME ACCELERATION OF PERFUSION TIME SERIES ACQUISITION BY MAGNETIC RESONANCE IMAGING OF CYCLIC ARTERIAL SPIN MARKING
US7545142B2 (en) * 2005-11-01 2009-06-09 Beth Israel Deaconess Medical Center, Inc. Arterial spin labeling with pulsed radio frequency sequences
WO2007123799A2 (en) * 2006-04-03 2007-11-01 The Trustees Of The University Of Pennsylvania Assessing subject's reactivity to psychological stress using fmri
US8195274B2 (en) * 2007-04-27 2012-06-05 The Regents Of The University Of California Mapping of vascular perfusion territories
EP2189112A1 (en) * 2008-11-24 2010-05-26 Bracco Research S.A. Real-time perfusion imaging and quantification
WO2010108161A2 (en) * 2009-03-19 2010-09-23 The Regents Of The University Of California Multi-phase pseudo-continuing arterial spin labeling

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1285044A (zh) * 1998-01-05 2001-02-21 耐克麦德英梅金公司 磁共振研究方法
CN1494871A (zh) * 2002-05-13 2004-05-12 ��ʽ���綫֥ 磁共振成像仪、流量量化仪及流量量化方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
Thomas T. Liu等."Analysis and Design of Perfusion-Based Event-Related fMRI Experiments".《NeuroImage》.2002,(第16期),269-282.

Also Published As

Publication number Publication date
CN101329388A (zh) 2008-12-24
DE102007028901B4 (de) 2010-07-22
JP2009000525A (ja) 2009-01-08
DE102007028901A1 (de) 2009-01-02
US8233961B2 (en) 2012-07-31
US20080319302A1 (en) 2008-12-25

Similar Documents

Publication Publication Date Title
CN101329388B (zh) 借助磁共振设备自动确定灌注的方法和装置
Jansz et al. Metric optimized gating for fetal cardiac MRI
Markl et al. Time‐resolved three‐dimensional phase‐contrast MRI
CN104101849B (zh) 磁共振系统控制序列的确定
US8125484B2 (en) Method, apparatus and user interface for determining an arterial input function used for calculating hemodynamic parameters
CN104095633B (zh) 用于产生图像数据的方法
CN101856228B (zh) 磁共振成像装置
CN103800006B (zh) 磁共振成像设备和获取功能图像的方法
JP6979151B2 (ja) 磁気共鳴イメージング装置及び磁気共鳴画像処理方法
US9693706B2 (en) Image processing of flowing fluids within imaged atonomical structure
CN102525466A (zh) 图像处理装置和磁共振成像装置
CN108629816A (zh) 基于深度学习进行薄层磁共振图像重建的方法
US8099149B2 (en) MRI method for quantification of cerebral perfusion
WO2011115956A1 (en) Systems and methods for detection and prediction of brain disorders based on neural network interaction
JP2007525250A (ja) 組織血流の遅延−補償型算出方法
US11647915B2 (en) Systems and methods for medical imaging incorporating prior knowledge
CN103080763A (zh) 用于给麻醉系统发送控制消息的磁共振成像系统、计算机系统以及计算机程序产品
KR20140120860A (ko) 영상 데이터를 생성하는 방법
CN103654781A (zh) 借助激励检查对象进行磁共振成像的方法和磁共振设备
CN106456046A (zh) 经改进的多时相动态对比增强磁共振成像的方法
CA2686150A1 (en) Image processing method
US11619694B2 (en) Magnetic resonance imaging apparatus and medical data processing apparatus
US20040102691A1 (en) Method and system for correcting images and measurements acquired by MRI
CN1689515A (zh) 用于对动态信号进行核自旋断层造影测量的方法
JP7337603B2 (ja) 磁気共鳴イメージング装置

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20130911

Termination date: 20200623