CN102138792A - 磁共振成像设备 - Google Patents

磁共振成像设备 Download PDF

Info

Publication number
CN102138792A
CN102138792A CN2011100373484A CN201110037348A CN102138792A CN 102138792 A CN102138792 A CN 102138792A CN 2011100373484 A CN2011100373484 A CN 2011100373484A CN 201110037348 A CN201110037348 A CN 201110037348A CN 102138792 A CN102138792 A CN 102138792A
Authority
CN
China
Prior art keywords
blood vessel
imaging apparatus
magnetic resonance
data
vessel
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.)
Pending
Application number
CN2011100373484A
Other languages
English (en)
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.)
GE Medical Systems Global Technology Co LLC
Original Assignee
GE Medical Systems Global Technology Co LLC
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 GE Medical Systems Global Technology Co LLC filed Critical GE Medical Systems Global Technology Co LLC
Publication of CN102138792A publication Critical patent/CN102138792A/zh
Pending legal-status Critical Current

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/5635—Angiography, e.g. contrast-enhanced angiography [CE-MRA] or time-of-flight angiography [TOF-MRA]
    • 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/56308—Characterization of motion or flow; Dynamic imaging
    • G01R33/56316—Characterization of motion or flow; Dynamic imaging involving phase contrast techniques
    • 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/56509—Correction of image distortions, e.g. due to magnetic field inhomogeneities due to motion, displacement or flow, e.g. gradient moment nulling

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)
  • Radiology & Medical Imaging (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Biophysics (AREA)
  • Molecular Biology (AREA)
  • Veterinary Medicine (AREA)
  • Condensed Matter Physics & Semiconductors (AREA)
  • Signal Processing (AREA)
  • Pathology (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Medical Informatics (AREA)
  • General Physics & Mathematics (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • Public Health (AREA)
  • Vascular Medicine (AREA)
  • Hematology (AREA)
  • Cardiology (AREA)
  • Physiology (AREA)
  • Magnetic Resonance Imaging Apparatus (AREA)

Abstract

磁共振成像设备(1)执行用于生成每个自旋的与该自旋的流速对应的相移的脉冲序列、由此从受检者(13)获取磁共振信号，并且根据所述磁共振信号的每个来确定所述受检者(13)的各血管的位置。磁共振成像设备(1)包括血管位置指定装置(102)，用于根据所述磁共振信号的信号强度随时间的变化以及根据所述自旋的流速随时间的变化来指定各血管的位置。

Description

磁共振成像设备
技术领域
本发明涉及用于确定受检者的各血管的位置的磁共振成像设备及其程序。
背景技术
在执行取决于血流速度的脉冲序列时，可预先测量血流速度。作为测量血流速度的一种方法，已知的是一种执行扫描的方法，用于测量血流速度，从而使操作人员从每个所获取磁共振图像中找出血管以及围绕这个血管以作为ROI(感兴趣区域)。
[专利文献1]日本未审查专利公开No.2005-305151
发明内容
技术问题
但是出现以下问题：因为当例如血管较小时，设置感兴趣区域ROI的工作变得很麻烦。因此希望解决此问题。
问题的解决方案
本发明的一个方面是一种磁共振成像设备，它执行用于生成每个自旋的与该自旋的流速对应的相移的脉冲序列，由此从受检者获取磁共振信号，并且根据磁共振信号的每个来确定受检者的各血管的位置，包括：血管位置指定装置，用于根据磁共振信号的信号强度随时间的变化以及自旋的流速随时间的变化来指定各血管的位置。
本发明的另一个方面是磁共振成像设备的一种程序，它执行用于生成每个自旋的与该自旋的流速对应的相移的脉冲序列，由此从受检者获取磁共振信号，并且根据磁共振信号的每个来确定受检者的各血管的位置，其中，提供程序以执行血管位置指定过程，用于根据磁共振信号的信号强度随时间的变化以及自旋的流速随时间的变化来指定各血管的位置。
发明的有利效果
本发明能够易于根据各磁共振信号的信号强度随时间的变化以及流速随时间的变化来判定各血管的位置。
通过以下对附图所示的本发明的优选实施例进行的描述，本发明的其它目的和优点将会非常明显。
附图说明
图1是示出根据本发明的第一实施例的磁共振成像设备1的简图。
图2是示出MRI设备1的处理流程的简图。
图3是示出受检者13的层片SL的位置和通过相衬度(phasecontrast)方法所得到的电影图像的简图。
图4是确定最大值c_max(x，y)时使用的说明图。
图5是示出表示等式(2)或等式(3)对于层片SL的平面的各位置(x，y)是否成立的二值图像的一个示例的简图。
图6是示意示出所提取血管区域的简图。
图7是用于描述用于确定数据c(x，y，t)的时间方向的相关性关于彼此相邻的像素是否较高的方法的一个示例的简图。
图8是示出根据第三实施例的流程的简图。
图9是用于说明根据第三实施例的流程的简图。
图10是示出根据第四实施例的处理流程的简图。
图11是用于说明根据第四实施例的处理流程的简图。
图12是用于说明根据第四实施例的另一个处理流程的简图。
具体实施方式
虽然下面将说明用于执行本发明的模式，但是用于执行本发明的模式并不局限于以下模式或实施例。
(1)第一实施例
图1是示出根据本发明的第一实施例的磁共振成像设备1的简图。
磁共振成像(MRI(磁共振成像))设备1具有磁场发生器2、台架3、托架4、接收线圈5等等。
磁场发生器2具有保持受检者13的膛21、超导线圈22、梯度线圈23和发射线圈24。超导线圈22施加静态磁场B0，并且梯度线圈23在频率编码方向、相位编码方向和层片选择方向施加梯度磁场。发射线圈24发射RF脉冲。顺便提到，虽然当前实施例中使用超导线圈22，但是永磁体可用于代替超导线圈22。
托架4配置成使得能够从台架3移动到膛21。受检者13由托架4传送到膛21。
接收线圈5附连到受检者13的各腿部13a。接收线圈5接收从受检者13所生成的每个磁共振信号。
MRI设备1还具有定序器6、发送器7、梯度磁场电源8、接收器9、中央处理器10、输入装置11和显示装置12。
在中央处理器10的控制下，定序器6向发送器7传送与脉冲序列的各RF脉冲有关的信息(中心频率、带宽等)，并且向梯度磁场电源8发送与梯度磁场有关的信息(梯度磁场的强度等)。
发送器7根据从定序器6所传送的信息来输出用于驱动发射线圈24的驱动信号。
梯度磁场电源8根据从定序器6所发送的信息来输出用于驱动梯度线圈23的驱动信号。
接收器9对接收线圈5所接收的每个磁共振信号进行例如数字转换等信号处理，并且将其输出给中央处理器10。
中央处理器10控制MRI设备1的相应部分的操作，以便实现MRI设备1的各种操作，例如向定序器6和显示装置12传送必要信息以及根据从接收器9所接收的各信号来重构图像。中央处理器10例如由计算机来配置。中央处理器10具有图像生成装置101和血管位置指定装置102。图像生成装置101生成图像CIk(其中k＝1至m)，稍后进行描述。血管位置指定装置102根据图像CIk(其中k＝1至m)来指定血管位置。顺便提到，中央处理器10通过执行预定程序来用作图像生成装置101和血管位置指定装置102。
输入装置11响应操作人员14的操纵而向中央处理器10输入各种指令。显示装置12在其上显示各种信息。
按以上所述来配置MRI设备1。接下来将说明MRI设备1的处理流程。
图2是示出MRI设备1的处理流程的简图。将说明处理流程，同时在描述图2时根据需要参照图3至图6。顺便提到，将对于确定受检者13的各腿部13a中的每个血管的位置的示例来进行以下描述。但是，本发明可适用于确定受检者13的任意部位或区域中的血管、如受检者13的腹部区域中的血管的位置的情况。
在步骤S1，操作人员14将层片SL设置到受检者13的腿部13a(参见图3(a))。顺便提到，虽然图3(a)中仅设置一片层片SL，但是可设置多片层片。
在已经设置层片SL之后，执行使用相衬度方法的脉冲序列，以便从层片SL获取磁共振信号，并且由此生成电影图像，其中的每个取决于磁共振信号的强度和各自旋的流速。在相衬度方法中，自旋的相移的幅值能够按照自旋的流速来改变。相应地，与各自旋的流速有关的信息能够通过经由相衬度方法获取的磁共振信号来获得。在第一实施例中，在以脉冲序列改变梯度磁场的极性的同时执行成像两次，由此获取复合数据f1和f2。图像生成装置101(参见图1)根据这些复合数据f1和f2(参见图3(b))来生成电影图像，其中的每个取决于各磁共振信号的强度和各自旋的流速。
图3(b)是示意示出取决于每个磁共振信号的强度和各自旋的流速的各电影图像的简图。
例如，图像CIk(其中k＝1至m)能够通过求复合数据f1和f2的微分来获得。图像CIk(其中k＝1至m)中的各像素的位置和时间以(x，y，t)来表达。各像素所表示的数据以c(x，y，t)来表达。在当前实施例中，数据c(x，y，t)由下式(1)定义：
c(x，y，t)
＝a(x，y，t)*sin(π*v(x，y，t)/VENC/2)        …(1)
其中，a(x，y，t)：在各像素的位置和时间(x，y，t)的信号强度，
v(x，y，5)：在各像素的位置和时间(x，y，t)的自旋的流速，以及
VENC：速度编码的梯度量。
关于在x＝xi和y＝yj处的像素的数据c(xi，yj，t1)、c(xi，yj，t2)、c(xi，yj，tm)在图3(b)中作为代表示出。在已经生成图像CI1至CIm之后，操作人员14进行步骤S2。
在步骤S2，血管位置指定装置102(参见图1)首先计算数据c(x，y，t)的绝对值|c(x，y，t)|，并且确定层片SL的平面的各位置(x，y)的绝对值|c(x，y，t)|在时间方向的最大值c_max(x，y)(参见图4)。
图4是确定最大值c_max(x，y)时使用的说明图。
图4(a)是示出图像CI1至CIm的简图。
当例如确定层片SL的平面的位置(xi，yj)处的最大值c_max(x，y)时，血管位置指定装置102使用位置(xi，yj)处的数据c(xi，yj，t)(参见图4(b))。
图4(b)是示出数据序列Cij的简图，其中位置(xi，yj)处的数据c(xi，yj，t)按照时间序列来排列。
血管位置指定装置102确定关于数据序列Cij的各数据c(xi，yj，t)的绝对值|c(xi，yj，t)|，并且确定各绝对值|c(xi，yj，t)|在时间方向的最大值c_max(xi，yj)。图4(b)中，在时间tα的数据c(xi，yj，tα)的绝对值|c(xi，yj，tα)|取时间方向的最大值。相应地，最大值c_max(xi，yj)按照下式(2)来表达：
c_max(xi，yj)＝|c(xi，yj，tα)|    …(2)
因此，位置(xi，yj)处的绝对值|c(xi，yj，t)|在时间方向的最大值c_max(xi，yj)能够根据等式(2)来计算。
以上描述中已经示出用于计算层片SL的平面的位置(xi，yj)处的最大值c_max(xi，yj)的过程。但是，能够通过相同过程来确定层片SL的平面的任何其它位置(x，y)处的最大值c_max(x，y)。例如，层片SL的平面的位置(xp，yq)处的最大值c_max(xp，yq)(参见图4(a))能够从其中位置(xp，yq)处的数据c(xp，yq，t)按时间序列排列的数据序列Cpq(参见图4(c))来计算。血管位置指定装置102确定关于数据序列Cpq的数据c(xp，yq，t)的绝对值|c(xp，yq，t)|，并且确定绝对值|c(xp，yq，t)|在时间方向的最大值c_max(xp，yq)。在图4(c)的数据序列Cpq中，在时间tβ的数据c(xp，yq，tβ)的绝对值|(xp，yq，tβ)|取时间方向的最大值。相应地，最大值c_max(xp，yq)按照下式(3)来表达：
c_max(xp，yq)＝|c(xp，yq，tβ)|    …(3)
因此，位置(xp，yq)处的绝对值|c(xp，yq，t)|在时间方向的最大值c_max(xp，yq)能够根据等式(3)来计算。
在按照上述过程对层片SL的平面的各位置(x，y)确定了数据c(x，y，t)的绝对值|c(x，y，t)|在时间方向的最大值c_max(x，y)之后，操作人员14进行步骤S3。
在步骤S3，血管位置指定装置102确定在步骤S2所确定的最大值c_max(x，y)是否小于阈值c_limit。一般来说，最大值c_max(x，y)趋向于在产生于血管的每个磁共振信号的情况下变大，而最大值c_max(x，y)趋向于在产生于固定组织的每个磁共振信号以及来自受检者13体外的信号的情况下变小。因此，能够确定当下式(4)成立时，它表示产生于固定组织或噪声的磁共振信号。另一方面，能够确定当下式(5)成立时，存在它表示产生于血管的磁共振信号的高可能性。
c_max(x，y)＜c_imit     …(4)
c_max(x，y)≥c_limit    …(5)
顺便提到，c_limit能够通过迭代计算等等来优化。
例如，在层片SL的平面的位置(xi，yj)，最大值c_max(xi，yj)大于阈值c_limit，如图4(b)所示。因此，由于等式(5)在层片SL的平面的位置(xi，yj)成立，所以认为血管的可能性较高。
另一方面，在图4(c)所示的层片SL的平面的位置(xp，yq)处最大值c_max(xp，yq)小于阈值c_limit。因此，由于等式(4)在层片SL的平面的位置(xp，yq)成立，所以认为固定组织或噪声的可能性较高(即，血管的可能性较低)。
类似地，确定等式(4)或(5)关于层片SL的平面的任何其它位置(x，y)是否成立(参见图5)。
图5是示出表示等式(4)或(5)对于层片SL的平面的各位置(x，y)是否成立的二值图像的一个示例的简图。
图5(a)是示出层片SL的平面的简图，以及图5(b)示出表示等式(4)或(5)在位于图5(a)所示的层片SL的部分区域R之内的各位置处是否成立的二值图像。
图5(b)中，所示加对角阴影的各像素指示在层片SL的区域R中等式(4)成立的位置(即，固定组织的可能性或者受检者体外的可能性较高的位置)。例如，由于在位置(xp，yq)处根据数据序列Cpq所计算的最大值c_max(xp，yq)在等式(4)中成立，所以像素P(xp，yq)指示固定组织或者受检者体外的可能性较高的位置。
另一方面，每个开放像素(open pixel)指示在层片SL的区域R中等式(5)成立的位置(x，y)(即血管的可能性较高的位置)。例如，由于在位置(xi，yj)处根据数据Cij所计算的最大值c_max(xi，yj)满足等式(5)，所以像素P(xi，yj)指示血管的可能性较高的位置。
因此，血管的可能性较高的各像素(开放像素)可通过确定等式(4)是否成立来指定。虽然为了说明方便起见在图5的层片SL的部分区域R中示出血管的可能性较高的各像素(开放像素)，但是血管的可能性较高的像素实际上在层片SL的整个区域上来指定。在已经执行步骤S3之后，操作人员14进行步骤S4。
在步骤S4，血管位置指定装置102(参见图1)耦合血管的可能性较高的像素(图5所示的开放像素)内彼此相邻的像素，并且提取血管区域(参见图6)。
图6是示意示出所提取血管区域的简图。
将邻接像素相互耦合实现了血管区域R1和R2的提取。
顺便提到，由于像素P(xr，ys)属于图6中的开放像素，所以它对应于在步骤S3判断为血管的可能性较高的像素。但是，像素P(xr，ys)由固定组织的像素或者位于受检者13体外的像素包围(图6中所示加对角阴影的像素)。因此，甚至在判断为血管的可能性较高的像素的情况下，认为在像素由固定组织的像素或者受检者13体外的像素包围的地方，血管的可能性较低。因此，将像素判定为没有对应于血管。
该流程按照上述方式结束。
一般来说，最大值c_max(x，y)趋向于在产生于血管的每个磁共振信号的情况下变大，而最大值c_max(x，y)趋向于在位于受检者13体外的信号(噪声)的情况下变小。相应地，血管区域能够通过从取决于信号强度和流速的图像CIk计算层片SL的每一个位置(x，y)的最大值c_max(x，y)来提取。
在第一实施例中，各血管的位置根据数据c(x，y，t)的绝对值|c(xi，yj，t)|在时间方向的最大值c_max(x，y)来指定。但是，由于在没有使数据c(x，y，t)成为负值的情况下不需要确定绝对值|c(x，y，t)|，所以血管的位置可根据数据c(x，y，t)在时间方向的最大值来指定。此外，对数据c(x，y，t)进行加权，并且各血管的位置可根据加权的数据c(x，y，t)来指定。
顺便提到，数据c(x，y，t)属于取决于信号强度a(x，y，t)和流速v(x，y，t)(参见等式(1))的数据。相应地，各血管的位置可通过以下来指定：确定信号强度a(x，y，t)和流速v(x，y，t)而无需确定数据c(x，y，t)，以及分析信号强度a(x，y，t)随时间的变化以及流速v(x，y，t)随时间的变化。
(2)第二实施例
将说明第二实施例，同时参照图2所示的流程。顺便提到，由于第二实施例在步骤S1至S3方面与第一实施例相同，所以省略步骤S1至S3的描述，因此仅说明步骤S4。
在第二实施例的步骤S4，血管位置指定装置102(参见图1)确定关于在步骤S3判定为各血管的可能性较高的像素(图5所示的开放像素)中彼此相邻的像素，数据c(x，y，t)的时间方向的相关性是否较高。
图7是用于说明用于确定关于彼此相邻的像素，数据c(x，y，t)的时间方向的相关性是否较高的方法的一个示例的简图。
当关于例如邻接像素P(xi，yj)和P(xi，yj-1)来确定数据c(x，y，t)的时间方向的相关性时，可计算在像素P(xi，yj)的数据序列Cij与在像素P(xi，yj-1)的数据序列Ci，j-1之间的相关性系数COR。一般来说，数据c(x，y，t)的时间方向的相关性系数COR趋向于在各血管的像素变大。因此，当相关性系数COR较大(例如COR＞0.8)时，认为对应像素是血管像素。另一方面，当相关性系数COR较小(例如COR≤0.8)时，认为血管像素的可能性较低。因此，即使在步骤S3错误地确定非对应于血管的像素属于血管的可能性较高的像素，也能够通过在步骤S4计算相关性系数COR将其从血管像素中排除。在第二实施例中，当确定数据c(x，y，t)的时间方向的相关性较高时，邻接像素相互耦合，并且由此提取各血管区域。因此能够以高准确度来提取血管区域。
顺便提到，由于静脉的流速比动脉要慢，所以相关性系数COR的值变小，其取决于静脉的流速，而与静脉的各像素的存在无关，使得可判断不是血管像素。虽然像素P(xi，yj)已经确定为对应于其中例如血管的区域R1属于静脉区域的血管像素，可能判断像素P(xj，yj-1)不是血管像素。作为一种用于避免这类误判的方法，考虑例如其中确定在像素P(xi，yj)的数据序列Cij的平均值M1及其标准偏差σ1以及在像素P(xi，yj-1)的数据序列Ci，j-1的平均值M2及其标准偏差σ2并且由此执行F测试和T测试的情况。一般已知，当静脉像素相互比较时，数据c(x，y，t)的数据序列的平均值趋向于近似为相同值，并且数据c(x，y，t)的数据序列的平均值也趋向于近似为相同值。因此，允许F测试和T测试通过时，对应像素可判断为静脉像素。因此，当允许F测试和T测试通过时，即使从血管像素中排除像素P(xi，yj-1)，但从相关性系数COR的值进行确定的情况下，像素(P(xi，yj-1)也能够确定为对应于静脉像素，由此使得能够以较高准确度来提取血管区域。
(3)第三实施例
将参照图8所示的流程来说明第三实施例。
图8是示出根据第三实施例的流程的简图。
由于第三实施例在步骤S1和S2方面与第一实施例相同，所以将省略步骤S1和S2的描述。在步骤S2已经结束之后，操作人员14进行步骤S21。
在步骤S21，血管位置指定装置102(参见图1)消除了伪像(artifact)的可能性较高的像素。顺便提到，虽然以下描述将对像素P(xi，yj)的伪像可能性较高的一个示例来进行，但是能够通过相似方法来确定其它像素是否也是伪像。
在时间方向的数据c(xi，yj，t)之间的微分首先对在像素P(xi，yj)处的数据序列Cij(参见例如图4(b))来执行(参见图9)。
图9(a)是示意示出在像素P(xi，yj)的数据序列Cij的简图，以及图9(b)是示意示出通过在时间方向的数据序列Cij的数据c(xi，yj，t)之间微分所得到的差分数据序列Dij的简图。
差分数据序列Dij在时间tn(其中n＝1至m-1)的数据(以下称作“差分数据”)d(xi，yj，tn)在下式(6)中使用数据c(xi，yj，tn+1)和c(xi，yj，tn)来表达：
d(xi，yj，tn)
＝c(xi，yj，tn+1)-c(xi，yj，tn)    …(6)
因此，例如，在时间tk的差分数据d(xi，yj，tk)在下式(6’)中通过等式(6)来表达：
d(xi，yj，tk)
＝c(xi，yj，tk-1)-c(xi，yj，tk)    …(6’)
在已经确定差分数据d(xi，yj，tn)之后，确定差分数据d(xi，yj，tn)的绝对值|d(xi，yj，tn)|在时间方向的最大值d_max(xi，yj)。图9(b)中，在时间tα的差分数据的绝对值|d(xi，yj，tα)|取时间方向的最大值。相应地，最大值d_max(xi，yj)按照下式(7)来表达：
d_max(xi，yj)＝|d(xi，yj，tα)|
随后，将在差分数据序列Dij的最大值d_max(xi，yj)和在数据序列Cij的最大值c_max(xi，yj)相互进行比较。一般来说，数据c(x，y，t)在血流的情况下在时间方向平滑且逐渐地改变。因此，在血管像素的情况下，在差分数据序列Dij的最大值d_max(xi，yj)必然产生比在数据序列Cij的最大值c_max(xi，yj)小的值。另一方面，在例如伪像等异常信号的情况下，不必使最大值d_max(xi，yj)为小值。因此，当下式(8)成立时，能够判断为指示产生于伪像的磁共振信号。另一方面，当下式(9)成立时，能够确定产生于血管的磁共振信号的可能性较高。
d_max(xi，yj)＞const1*c_max(xi，yj)…(8)
d_max(xi，yj)≤const1*c_max(xi，yj)…(9)
顺便提到，const1具有经验值。
因此，伪像可通过确定等式(8)是否成立来排除，由此使得能够以较高准确度来提取血管。在排除伪像之后，操作人员14进行步骤S3和S4，其中提取血管区域。
顺便提到，在以上描述中，根据在差分数据序列Dij的最大值d_max(xi，yj)与在数据序列Cij的最大值c_max(xi，yj)之间的比较结果，来排除伪像。但是，确定差分数据序列Dij的标准偏差d_std(xi，yj)和数据序列Cij的标准偏差d_std(xi，yj)，并且由此可根据这些标准偏差d_std(xi，yj)和c_std(xi，yj)之间的比较结果来排除伪像。由于数据c(xi，yj，t)在血流的情况下在时间方向平滑且逐渐地改变，所以差分数据序列Dij的标准偏差d_std(xi，yj)必然产生比数据序列Cij的标准偏差c_std(xi，yj)小的值。另一方面，在例如伪像等异常信号的情况下，不必使标准偏差d_std(xi，yj)为小值。因此，当下式(10)成立时，能够判断为指示产生于伪像的磁共振信号。另一方面，当等式(11)成立时，能够判断产生于血管的磁共振信号的可能性较高。
d_std(xi，yj)＞const2*c_std(xi，yj)…(10)
d_std(xi，yj)≤const2*c_std(xi，yj)…(11)
顺便提到，const2指示经验值。
因此，伪像甚至能够通过比较标准偏差来排除，由此使得能够以较高准确度来执行血管的提取。
可考虑最大值d_max(xi，yj)与c_max(xi，yj)之间的比较结果以及标准偏差d_std(xi，yj)与c_std(xi，yj)之间的比较结果来排除伪像。
(4)第四实施例
图10是示出根据第四实施例的处理流程的简图。顺便提到，将说明图10的处理流程，同时在其描述时根据需要参照图11和图12。
在步骤S1，首先设置层片SL(参见图11(a))。在已经设置层片SL之后，执行使用相衬度方法的脉冲序列，以便从层片SL获取磁共振信号，并且由此生成各指示磁共振信号的强度的电影图像以及各取决于磁共振信号的强度和各自旋的流速的电影图像。在相衬度方法中，自旋的相移的幅值能够按照自旋的流速来改变。相应地，与各自旋的流速有关的信息能够通过经由相衬度方法获取磁共振信号来获得。在第四实施例中，在以脉冲序列改变梯度磁场的极性的同时执行成像两次，由此获取复合数据f1和f2。图像生成装置101(参见图1)根据这些复合数据f1和f2来生成电影图像，其中的每个取决于磁共振信号的强度和自旋的流速。图11(b)示出指示信号强度的强度图像AIk(其中k＝1至m)，以及图11(c)示出取决于信号强度和流速的图像CIk(其中k＝1至m)。强度图像AIk(其中k＝1至m)能够作为例如复合数据的绝对值|f1|(＝|f2|)来获得。图像CIk(其中k＝1至m)能够通过例如求复合数据f1和f2的微分来获得。
强度图像AI1至AIm的相应像素的位置和时间以(x，y，t)来表达。各像素所表示的信号强度以a(x，y，t)来表达。在x＝xi和y＝yj处的像素的信号强度a(xi，yj，t1)、a(xi，yj，t2)、c(xi，yj，tm)通常如图11(b)所示。顺便提到，由于各取决于信号强度和流速的图像CI1至CIm与第一实施例中采用的相似，所以将省略其描述。
在已经生成强度图像AI1至AIm以及各取决于信号强度和流速的图像CI1至CIm之后，操作人员14进行步骤S11。
在步骤S11，使用强度图像AI1至AIm对层片SL的平面的各位置(x，y)计算信号强度a(x，y，t)在时间方向的最大值a_max(x，y)。
图12是计算时间方向的信号强度a(x，y，t)的最大值a_max(x，y)时使用的说明图。
图12(a)是示出强度图像AI1至AIm的电影图像的简图。
当例如计算在层片SL的平面的位置(xt，yu)处的信号强度a(xt，yu，t)在时间方向的最大值a_max(xt，yu)时，在层片SL的平面的位置(xt，yu)处的信号强度a(xt，yu，t1)至a(xt，yu，tm)可取自强度图像AI1至AIm(参见图12(b))。
图12(b)示出指示信号强度a(xt，yu，t1)至a(xt，yu，tm)随时间的变化的强度数据序列Atu。确定强度数据序列Atu实现了在层片SL的平面的位置(xt，yu)处的信号强度在时间方向的最大值a_max(xt，yu)的计算。
在以上描述中已经示出了计算层片SL的平面的位置(xt，yu)处的信号强度a(xt，yu，t)在时间方向的最大值a_max(xt，yu)的过程。但是，在层片SL的平面的各其它位置(x，y)处的信号强度a(x，y，t)在时间方向的最大值也能够按照相似过程来确定。例如，在层片SL的平面的位置(xv，yw)(参见图12(a))处的信号强度a(xv，yw，t)在时间方向的最大值a_max(xy，yw)如图12(c)所示。能够从层片SL的平面的位置(xv，yw)处的信号强度a(xv，yw，t1)至a(xv，yw，tm)的强度数据序列Avw来确定最大值a_max(xv，yw)。
在上述过程中对于层片SL的平面的各位置(x，y)确定了时间方向的信号强度a(x，y，t)的最大值a_max(x，y)之后，操作人员14进行步骤S12。
在步骤S12，确定在步骤S11已经确定的时间方向的信号强度的最大值a_max(x，y)是否小于阈值a_limit。一般来说，在时间方向的信号强度的最大值a_max(x，y)趋向于在受检者13体内变大，而在时间方向的信号强度的最大值a_max(x，y)在受检者13体外变小。因此，能够确定当下式(12)成立时，能够判断噪声的可能性较高。另一方面，能够确定当下式(13)成立时，能够判断与噪声不同的信号(产生于受检者13体内的磁共振信号)的可能性较高。
a_max(x，y)＜a_limit    …(12)
a_max(x，y)≥a_limit    …(13)
顺便提到，a_limit能够通过迭代计算等等来优化。
例如，在层片SL的平面的位置(xt，yu)，信号强度在时间方向的最大值a_max(xt，yu)大于阈值a_limit，如图12(b)所示。因此，由于等式(13)在层片SL的平面的位置(xt，yu)成立，所以认为除噪声之外的信号(产生于受检者13体内的磁共振信号)的可能性较高。
另一方面，在图12(c)所示的层片SL的平面的位置(xv，yw)处，信号强度在时间方向的最大值a_max(xv，yw)小于阈值a_limit。因此，由于等式(12)在层片SL的平面的位置(xv，yu)成立，所以认为噪声(受检者体外的磁共振信号)的可能性较高。
类似地，确定等式(12)或(13)关于层片SL的平面的任何其它位置(x，y)是否成立。因此，通过确定在层片SL的整个平面上在时间方向的信号强度的最大值a_max(x，y)是否大于或等于阈值a_limit，能够有效地排除噪声。在步骤S12结束之后，操作人员14进行步骤S2。
由于步骤S2至S4与第一实施例中采用的相似，所以将省略其描述。
在第四实施例中，在步骤S12已经确定在时间方向的信号强度的最大值a_max(x，y)在层片SL的整个平面上是否大于或等于阈值a_limit。相应地，能够在步骤S4将像素相互耦合之前有效地排除噪声，由此使得能够以较高准确度来提取血管区域。
第四实施例已经说明了除了取决于信号强度和流速的图像CIk之外还生成强度图像AIk的示例。但是，除了强度图像AIk之外(或者代替强度图像AIk)还可生成指示各自旋流速的流速图像。由于静脉的流速比动脉要慢，所以能够通过生成流速图像来识别所提取血管区域是属于静脉血管还是动脉血管。
顺便提到，在第一至第四实施例的每个中使用了通过等式(1)所表达的数据c(x，y，t)。但是，由于数据取决于流速v(x，y，t)和信号强度a(x，y，t)，所以可使用与数据c(x，y，t)不同的数据。例如，可使用通过将信号强度a(x，y，t)与流速v(x，y，t)彼此相乘所得到的数据p(x，y，t)。在这种情况下，数据p(x，y，t)通过下式(14)来表达：
p(x，y，t)＝a(x，y，t)*v(x，y，t)…(14)
甚至当使用等式(14)代替等式(1)所定义的数据p(x，y，t)时，也能够指定各血管的位置能够。
可配置本发明的许多极为不同的实施例，而没有背离本发明的精神和范围。应当理解，本发明并不局限于说明书中的具体实施例，而是由所附权利要求书来定义。

Claims (10)

1.一种磁共振成像设备(1)，执行用于生成每个自旋的与该自旋的流速对应的相移的脉冲序列、由此从受检者(13)获取磁共振信号，并且根据所述磁共振信号的每个来确定所述受检者(13)的各血管的位置，包括：
血管位置指定装置(102)，用于根据所述磁共振信号的信号强度随时间的变化以及根据所述自旋的流速随时间的变化来指定各血管的位置。
2.如权利要求1所述的磁共振成像设备(1)，其中，所述血管位置指定装置(102)根据在所述受检者(13)的预定切割面的各位置处的磁共振信号的信号强度随时间的变化以及在其各位置处的自旋的流速随时间的变化来指定各血管的位置。
3.如权利要求1或2所述的磁共振成像设备(1)，其中，所述血管位置指定装置(102)根据取决于所述信号强度和所述流速的数据来指定各血管的位置。
4.如权利要求3所述的磁共振成像设备(1)，其中，对所述受检者(13)的所述预定切割面的各位置来确定指示所述数据随时间的变化的数据序列，并且各血管的位置根据所述数据序列来指定。
5.如权利要求4所述的磁共振成像设备(1)，其中，所述血管位置指定装置(102)确定所述预定切割面的各位置的所述数据序列中的所述数据的绝对值在时间方向的最大值，并且根据所述最大值来指定各血管的位置。
6.如权利要求4或5所述的磁共振成像设备(1)，其中，所述血管位置指定装置(102)根据所述预定切割面中彼此相邻的位置处的数据序列之间的相关性来指定各血管的位置。
7.如权利要求4至6中的任一项所述的磁共振成像设备(1)，其中，所述血管位置指定装置(102)确定所述预定切割面中彼此相邻的位置处的数据序列的平均值或标准偏差，并且根据所述平均值或所述标准偏差来指定所述血管的位置。
8.如权利要求4至7中的任一项所述的磁共振成像设备(1)，其中，所述血管位置指定装置(102)确定各数据序列的数据之间在时间方向的差分，并且根据所述差分来排除伪像。
9.如权利要求1至8中的任一项所述的磁共振成像设备(1)，其中，所述血管位置指定装置(102)根据所述磁共振信号的信号强度随时间的变化来排除噪声。
10.如权利要求1至9中的任一项所述的磁共振成像设备(1)，其中，所述血管位置指定装置(102)根据所述流速随时间的变化来确定所提取的血管是动脉还是静脉。
CN2011100373484A 2010-01-29 2011-01-28 磁共振成像设备 Pending CN102138792A (zh)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
JP2010019108A JP2011156078A (ja) 2010-01-29 2010-01-29 磁気共鳴イメージング装置およびプログラム
JP2010-019108 2010-01-29

Publications (1)

Publication Number Publication Date
CN102138792A true CN102138792A (zh) 2011-08-03

Family

ID=44342238

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2011100373484A Pending CN102138792A (zh) 2010-01-29 2011-01-28 磁共振成像设备

Country Status (3)

Country Link
US (1) US20110190622A1 (zh)
JP (1) JP2011156078A (zh)
CN (1) CN102138792A (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109901087A (zh) * 2017-11-10 2019-06-18 温伯格医学物理有限公司 用于提高磁共振成像中的信噪比的方法

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP6104510B2 (ja) * 2012-02-29 2017-03-29 東芝メディカルシステムズ株式会社 画像処理装置及び制御プログラム
US9014781B2 (en) * 2012-04-19 2015-04-21 General Electric Company Systems and methods for magnetic resonance angiography
CN105022719B (zh) * 2014-04-23 2019-06-28 Ge医疗系统环球技术有限公司 医学造影系统及方法

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5016637A (en) * 1988-05-27 1991-05-21 Hitachi, Ltd. MRI method
US6505064B1 (en) * 2000-08-22 2003-01-07 Koninklijke Philips Electronics, N.V. Diagnostic imaging systems and methods employing temporally resolved intensity tracing
JP2006130116A (ja) * 2004-11-08 2006-05-25 Hitachi Medical Corp 磁気共鳴イメージング装置
CN1817305A (zh) * 2005-02-08 2006-08-16 西门子公司 成像诊断系统
US20070057671A1 (en) * 2005-09-14 2007-03-15 The Government of the United States of America as represented by the Secretary of the Imaging and reconstruction of partial field of view in phase contrast MRI
CN1943515A (zh) * 2005-08-22 2007-04-11 西门子公司 在数字图像数据中识别注有造影剂的血管的方法
CN102707251A (zh) * 2011-12-12 2012-10-03 中国科学院深圳先进技术研究院 计算space序列信号的方法和系统及主动脉信号的采集方法
CN102749603A (zh) * 2011-04-21 2012-10-24 西门子公司 用于建立磁共振血管造影图像的方法

Family Cites Families (23)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4836209A (en) * 1986-08-07 1989-06-06 Stanford University NMR imaging of moving material using variable spatially selected excitation
US4718424A (en) * 1986-08-07 1988-01-12 Stanford University NMR imaging of blood flow by moment variation of magnetic gradients
US5031624A (en) * 1990-08-17 1991-07-16 Wisconsin Alumni Research Foundation Phase contrast, line-scanned method for NMR angiography
JP3940442B2 (ja) * 1995-04-28 2007-07-04 株式会社東芝 磁気共鳴イメージング装置
US5997883A (en) * 1997-07-01 1999-12-07 General Electric Company Retrospective ordering of segmented MRI cardiac data using cardiac phase
US7254437B2 (en) * 1998-04-17 2007-08-07 Kabushiki Kaisha Toshiba MR imaging providing tissue/blood contrast image
US20030053669A1 (en) * 2001-07-18 2003-03-20 Marconi Medical Systems, Inc. Magnetic resonance angiography method and apparatus
US6801800B2 (en) * 1999-11-29 2004-10-05 Kabushiki Kaisha Toshiba MR imaging using ECG-prep scan
US6718055B1 (en) * 2000-12-05 2004-04-06 Koninklijke Philips Electronics, N.V. Temporal and spatial correction for perfusion quantification system
US6611701B2 (en) * 2000-12-30 2003-08-26 Ge Medical Systems Global Technology Company, Llc Method and apparatus for fast breath-held 3D MR data acquisition using variable sampling
DE10117787A1 (de) * 2001-04-10 2002-10-17 Philips Corp Intellectual Pty MR-Verfahren zur Untersuchung eines zyklisch veränderlichen Objekts
US6968225B2 (en) * 2001-08-24 2005-11-22 General Electric Company Real-time localization, monitoring, triggering and acquisition of 3D MRI
US6650115B2 (en) * 2001-10-12 2003-11-18 The Board Of Trustees Of The Leland Stanford Junior University Variable density fourier velocity encoding MR imaging
US7024027B1 (en) * 2001-11-13 2006-04-04 Koninklijke Philips Electronics N.V. Method and apparatus for three-dimensional filtering of angiographic volume data
JP3455530B1 (ja) * 2001-12-14 2003-10-14 株式会社東芝 Ｍｒ信号受信装置及び磁気共鳴イメージング装置
JP4045143B2 (ja) * 2002-02-18 2008-02-13 テセラ・インターコネクト・マテリアルズ，インコーポレイテッド 配線膜間接続用部材の製造方法及び多層配線基板の製造方法
US6957097B2 (en) * 2002-04-17 2005-10-18 The Board Of Trustees Of The Leland Stanford Junior University Rapid measurement of time-averaged blood flow using ungated spiral phase-contrast MRI
US6922580B2 (en) * 2002-06-04 2005-07-26 Koninklijke Philips Electronics N.V. Blood flow gated MRI
US7343193B2 (en) * 2003-06-16 2008-03-11 Wisconsin Alumni Research Foundation Background suppression method for time-resolved magnetic resonance angiography
CN100479748C (zh) * 2004-03-26 2009-04-22 株式会社东芝 核磁共振成像装置
US7382132B1 (en) * 2005-04-29 2008-06-03 General Electric Company 6-channel array coil for magnetic resonance imaging
US7945308B2 (en) * 2005-12-27 2011-05-17 General Electric Company Systems, methods and apparatus for an endo-rectal receive-only probe
US20110230756A1 (en) * 2008-09-30 2011-09-22 University Of Cape Town Fluid flow assessment

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5016637A (en) * 1988-05-27 1991-05-21 Hitachi, Ltd. MRI method
US6505064B1 (en) * 2000-08-22 2003-01-07 Koninklijke Philips Electronics, N.V. Diagnostic imaging systems and methods employing temporally resolved intensity tracing
JP2006130116A (ja) * 2004-11-08 2006-05-25 Hitachi Medical Corp 磁気共鳴イメージング装置
CN1817305A (zh) * 2005-02-08 2006-08-16 西门子公司 成像诊断系统
CN1943515A (zh) * 2005-08-22 2007-04-11 西门子公司 在数字图像数据中识别注有造影剂的血管的方法
US20070057671A1 (en) * 2005-09-14 2007-03-15 The Government of the United States of America as represented by the Secretary of the Imaging and reconstruction of partial field of view in phase contrast MRI
CN102749603A (zh) * 2011-04-21 2012-10-24 西门子公司 用于建立磁共振血管造影图像的方法
CN102707251A (zh) * 2011-12-12 2012-10-03 中国科学院深圳先进技术研究院 计算space序列信号的方法和系统及主动脉信号的采集方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
ALBERT C.S.CHUNG,ET AL: "Fusing speed and phase information for vascular segmentation of phase contrast MR angiograms", 《MEDICAL IMAGE ANALYSIS》 *
ALBERT C.S.CHUNG,ET AL: "Vascular Segmentation of Phase Contrast Magnetic Resonance Angiograms Based on Statistical Mixture Modeling and Local Phase Coherence", 《IEEE TRANSCATIONS ON MEDICAL IMAGING》 *
NOAM ALPERIN，ET AL: "PUBS:Pulsatility-Based Segmentation of Lumens Conducting Non-steady Flow", 《MAGNETIC RESONANCE IN MEDICINE》 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109901087A (zh) * 2017-11-10 2019-06-18 温伯格医学物理有限公司 用于提高磁共振成像中的信噪比的方法

Also Published As

Publication number Publication date
JP2011156078A (ja) 2011-08-18
US20110190622A1 (en) 2011-08-04

Similar Documents

Publication Publication Date Title
US10168405B2 (en) Method and apparatus for quantifying properties of an object through magnetic resonance imaging (MRI)
CN105934683B (zh) 具有对主磁场变化的估计的dixon型水/脂肪分离的mri
CN111265206B (zh) 一种磁共振血管成像方法、装置及设备
JP5559848B2 (ja) 複数タイプの磁気共鳴映像を同時に生成する装置及びその方法
CN111417979A (zh) 时间分辨的血管造影图像与空间分辨血管造影图像的组合
US10175330B2 (en) Dixon MR imaging with suppression of flow artifacts
JP6417406B2 (ja) 強調磁化率コントラストによるｍｒイメージング
CN102621510A (zh) 用于抑制mr成像中伪影的系统
CN106415298A (zh) 涉及对血管造影加权图像和灌注加权图像的采集的mri
US10429472B2 (en) Magnetic resonance imaging apparatus and method for magnetic resonance imaging with copying and setting of parameter values
JP2018008114A (ja) 画像処理装置、磁気共鳴イメージング装置及び画像処理方法
CN106456046A (zh) 经改进的多时相动态对比增强磁共振成像的方法
WO2017221654A1 (ja) 磁気共鳴イメージング装置、画像処理装置、及び拡散強調画像計算方法
CN112292606A (zh) 在磁共振成像中跟踪图像质量
CN102138792A (zh) 磁共振成像设备
US8160342B2 (en) System and method for processing data signals
JP2016140683A (ja) 医用画像処理装置および磁気共鳴イメージング装置
US5531223A (en) Method and apparatus of nuclear magnetic resonance imaging with nonlinearly processed image display
JP4912808B2 (ja) 磁気共鳴イメージング装置
JP5101896B2 (ja) 画像表示装置、画像表示方法、および、磁気共鳴イメージング装置
JP2012105982A (ja) 磁気共鳴イメージング装置及び磁気共鳴イメージング方法
US20140009155A1 (en) Labeling area determining apparatus, magnetic resonance apparatus and method for determining tilt of labeling area
WO2018065618A1 (en) Combined arterial spin labeling and magnetic resonance fingerprinting
CN108896940B (zh) 用于补偿mri的梯度脉冲的系统和方法
JP2021512678A (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
C02 Deemed withdrawal of patent application after publication (patent law 2001)
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20110803