CN1247151C - 磁共振成像系统 - Google Patents

磁共振成像系统 Download PDF

Info

Publication number
CN1247151C
CN1247151C CNB021272611A CN02127261A CN1247151C CN 1247151 C CN1247151 C CN 1247151C CN B021272611 A CNB021272611 A CN B021272611A CN 02127261 A CN02127261 A CN 02127261A CN 1247151 C CN1247151 C CN 1247151C
Authority
CN
China
Prior art keywords
data set
data
root
mean
operator
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
CNB021272611A
Other languages
English (en)
Other versions
CN1395904A (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.)
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 CN1395904A publication Critical patent/CN1395904A/zh
Application granted granted Critical
Publication of CN1247151C publication Critical patent/CN1247151C/zh
Anticipated expiration legal-status Critical
Expired - Fee Related 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/561—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution by reduction of the scanning time, i.e. fast acquiring systems, e.g. using echo-planar pulse sequences
    • G01R33/5613—Generating steady state signals, e.g. low flip angle sequences [FLASH]
    • 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

Landscapes

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

Abstract

本发明涉及一种磁共振(MR)成像法和一种磁共振成像(MRI)系统。本发明的目的是产生无伪象带的图像。在从组成K空间的视图中获取从数据fv(0)到fv(N-1)范围中的数据fv(k)前，扫描设备通过顺序改变相位编码的相位对稳定状态自由进动(SSFP)中的数据重复采集N次(这里N表示2的幂)。此时，射频(RF)脉冲的相位根据一个表达式360°□v□k/N而变化。如果操作者指定傅里叶变换(FT)成像，数据处理设备就对数据fv(k)相应于用RF脉冲表示的每个相位进行傅里叶变换，因此产生数据Fv(n)。如果操作者不指定傅里叶变换成像，数据fv(k)实际上就被认为是Fv(n)。此后，为了产生数据Av，就对数据Fv(n)执行操作者选择的至少加权加法、最大强度投影(MIP)处理和均方根转换中的任一过程。图像根据数据Av进行重建。

Description

磁共振成像系统
本发明的背景技术
本发明涉及一种磁共振(MR)成像法和一种磁共振成像(MRI)系统。更具体地说，本发明涉及一种能够重建良好质量图像的MR成像法和MRI系统。
专利No.2,898,329已揭示了一种根据下面步骤的MR成像法：
(1)在从组成K空间的视图v中获得数据fv(0)前通过顺序改变相位编码的相位对稳定状态自由进动(SSFP)中的数据重复采集；
(2)通过顺序改变相位编码的相位和交替地对RF脉冲移相180度SSFP中的数据重复采集，从而从组成K空间的视图v中获得数据fv(1)；
(3)对fv(0)和fv(1)进行相加和相减，以得到如下表达的数据Av：
Av＝0.5*Fv(0)+0.5*Fv(1)或
Av＝0.5*Fv(0)-0.5*Fv(1)；和
(4)根据所产生的数据Av重建图像。
根据专利No.2,898,329揭示的MR成像法，在某些情况下能够产生良好质量的图像。但是，在其它情况下能仅产生质量差的图像(例如，由不均匀磁场引起的具有伪象带的图像)。
本发明的概述
因此，本发明的目的是一种在任何常用MR成像法仅提供质量差的图像的情况下提供一种能够重建良好质量图像的MR成像法和MRI系统。
根据本发明的第一种方案，提供一种具有下面描述步骤的MR成像法。也就是说，在第一步骤(1)中，在从组成K空间的视图v中获得从数据fv(0)到数据fv(N-1)范围中的数据fv(k)前通过顺序改变相位编码的相位，对稳定状态自由进动(SSFP)中的数据重复采集N次(这里N等于2的幂)。此时，射频(RF)脉冲的相位根据一个表达式360°vk/N而变化。在第二步骤(2)中，如果操作者指定傅里叶变换(FT)成像，数据fv(k)相应于RF脉冲表示的相位进行相位编码，然后再进行傅里叶变换。这样就会产生数据Fv(n)。相反，如果操作者不指定傅里叶变换成像，数据fv(k)实际上就被认为是Fv(n)。在第三步骤(3)中，为了产生数据Av，就对数据Fv(n)执行操作者选择的至少加权加法、最大强度投影(MIP)处理和均方根转换中的任一过程。
根据本发明第一个方案中提供的MR成像法，操作者能够指定是否对数据fv(k)相应于用RF脉冲表示的每个相位进行傅里叶变换(FT)。傅里叶变换能够根据下面所述的原理指定数据的自由感应衰减(FID)分量和自旋回波分量或其激励回波分量中任一个应该占主导地位。
例如，N＝4时，如果k＝0，对所有数据采集的次数而言，RF脉冲的相位就设定为0。数据fv(0)的FID分量的极性与正Y方向(Y轴的正方向)一致，而数据的自旋回波分量或其激励回波分量与负Y方向(Y轴的负方向)一致。如果k＝1，RF脉冲的相位就顺序地设定为0，π/2，π，3π/2等。数据fv(1)的FID分量的极性与正Y方向一致，而数据的自旋回波分量或其激励回波分量与负X方向(X轴的负方向)一致。如果k＝2，RF脉冲的相位就交替地设定为0和π。数据fv(2)的FID分量的极性与正Y方向一致，而数据的自旋回波分量或其激励回波分量也与正Y方向一致。如果k＝3，RF脉冲的相位就顺序地设定为0，3π/2，π/2，π等。数据fv(3)的FID分量的极性与正Y方向一致，而数据的自旋回波分量或其激励回波分量与负X方向(X轴的负方向)一致。
由于数据Fv(0)＝fv(0)+fv(1)+fv(2)+fv(3)，FID分量被完整保留，因为由于上述极性而使自旋回波或激励回波分量被消除。实际上，情形与理想不一致。然而，在合成数据Fv(0)中，其FID分量是占主导地位。而且，由于数据Fv(1)＝fv(0)-j.fv(1)-fv(2)+j.fv(3)，自旋回波或激励回波分量被完整保留，因为FID分量由于上述极性被消除。因此，数据Fv(1)的自旋回波或激励回波分量占主导地位。通常，如果在数据Fv(n)中n假定为奇数值，FID分量就占主导地位。如果n假定为偶数值，自旋回波或激励回波分量就占主导地位。因此，能够指定FID分量和自旋回波或激励回波分量中任一个占主导地位。
此外，根据本发明第一个方案中提供的MR成像法，操作者能够从至少加权加法、MIP和均方根转换中选择任一个对数据Fv(n)执行处理。如果对数据执行加权加法处理，就可以指定FID分量和自旋回波或激励回波中任一个占主导地位。如果执行MIP处理，就能提高信噪比。而且，如果执行均方根转换，就能提高信噪比。
根据本发明第一个方案中提供的MR成像法，至少可从四种处理中选择执行处理。这样，能够在任何常用MR成像法仅能产生质量差的图像的情况下产生质量好的图像。
本发明者所作的研究表明：根据专利No.2,898,329揭示的MR成像法，如果对表示具有伪象带的可视图像的数据fv(0)和fv(1)执行均方根转换，就可产生无伪象带的图像。而且，如果增加重复次数N(例如，为8或更大)，并选择傅里叶变换和均方根转换，从而可在许多情况下产生质量好的图像。
根据本发明的第二个方案，提供一种包括下列步骤的MR成像法。也就是说，在第一步骤(1)中，在从组成K空间的视图v中获得从数据fv(0)到数据fv(N-1)范围中的数据fv(k)前通过顺序改变相位编码的相位，对SSFP中的数据重复采集N次(这里N等于2的幂)。此时，射频(RF)脉冲的相位根据一个表达式360°□v□k/N而变化。在第二步骤(2)中，为了产生数据Fv(n)，对数据fv(k)相应于用RF脉冲表示的每个相位进行傅里叶变换。在第三步骤(3)中，为了产生数据Av，就对数据Fv(n)执行操作者选择的至少加权加法、MIP处理和均方根转换中的任一过程。在第四步骤(4)中。根据所产生的数据Av重建图像。
根据本发明第二个方案中提供的MR成像法，对数据fv(k)相应于用RF脉冲表示的每个相位进行傅里叶变换。傅里叶变换能够指定数据的FID分量和自旋回波分量或其激励回波分量中任一个应该占主导地位。此外，操作者能够从至少加权加法、MIP和均方根转换中选择任一个对数据Fv(n)执行处理。如果执行加权加法处理，就指定FID分量和自旋回波或激励回波中任一个占主导地位。如果执行MIP处理，就能提高信噪比。如果执行均方根转换，就能提高信噪比。如上所述，可选择至少两种处理中的一个。这样，能够在任何常用MR成像法仅能产生质量差的图像的情况下产生质量好的图像。
本发明者所作的研究表明：增加重复次数N(例如，为8或更大)时，如果执行傅里叶变换和均方根转换，就可在许多情况下产生质量好的图像。
根据本发明的第三个方案，提供一种主要由下列步骤组成的MR成像法。也就是说，在第一步骤(1)中，在从组成K空间的视图v中获得从数据fv(0)到数据fv(N-1)范围中的数据fv(k)前通过顺序改变相位编码的相位，对SSFP中的数据重复采集N次(这里N等于2的幂)。此时，射频(RF)脉冲的相位根据一个表达式360°□v□k/N而变化。在第二步骤(2)中，数据fv(k)实际上可被认为是数据Fv(n)。在第三步骤(3)中，为了产生数据Av，就对数据Fv(n)执行操作者选择的至少加权加法、MIP处理和均方根转换中的任一过程。在第四步骤(4)中。根据所产生的数据Av重建图像。
根据本发明第三个方案中提供的MR成像法，操作者能够从至少加权加法、MIP和均方根转换中选择任一个对数据Fv(n)执行处理。如果执行加权加法处理，就指定FID分量和自旋回波或激励回波中任一个占主导地位。如果执行MIP处理，就能提高信噪比。如果执行均方根转换，就能提高信噪比。因此，操作者能够选择至少两种处理中的一个。这样，能够在任何常用MR成像法仅能产生质量差的图像的情况下产生质量好的图像。
本发明者所作的研究表明：根据专利No.2,898,329揭示的MR成像法，如果对显示示具有伪象带的可视图像的数据fv(0)和fv(1)执行均方根转换，在某些情况下，就可产生无伪象带的图像。
根据本发明的第四种方案，在任一前述的MR成像法基础上提供一种MR成像法。特别是，一种用于获取SSFT中数据的脉冲序列能够启动FID信号和回波的同时采集。而且，每个磁场梯度表示随时间变化的强度，在一个重复时间TR过程中检测的其强度值被积分为0。
众所周知，SSFT中的数据采集采用不同的脉冲序列。例如，大家知道用于稳定状态采集(FIESTA)快速成像所使用的一种脉冲序列和用于TrueSSFT的-种脉冲序列。
根据本发明第四种方案提供的MR成像法，可采用FIESTA所使用的脉冲序列。
根据本发明的第五种方案，在任一前述的MR成像法基础上提供一种MR成像法。这里，从相应于用RF脉冲表示的每个相位执行傅里叶变换中产生的数据Fv(n)表达如下：
Fv ( n ) = Σ K = 0 N - 1 fv ( k ) . exp { - j . n . 2 π . k / M }
根据本发明第五种方案提供的MR成像法，就能对数据fv(k)相应于用RF脉冲表示每个相位执行傅里叶变换。
根据本发明的第六种方案，在任一前述的MR成像法基础上提供一种MR成像法。这里，从N＝2时执行加权加法中产生的数据Av表达如下：
Av＝0.5*Fv(0)+0.5*Fv(1)
根据本发明第六种方案提供的MR成像法，可使数据fv(k)的FID分量占主导地位。
根据本发明的第七种方案，在任一前述的MR成像法基础上提供一种MR成像法。这里，从N＝2时执行加权加法中产生的数据Av表达如下：
Av＝0.5*Fv(0)-0.5*Fv(1)
根据本发明第七种方案提供的MR成像法，可使数据fv(k)的自旋回波或激励回波分量占主导地位。
根据本发明的第八种方案，在任一前述的MR成像法基础上提供一种MR成像法。这里，从解函数最大值{}以提供一最大值的MIP处理中产生的数据Av表达如下：
Av＝max{Fv(0)，etc，Fv(N-1)}
根据本发明第八种方案提供的MR成像法，可采用所有N个信号中最大值的一个信号。因此，在许多情况下能够提高信噪比。
根据本发明的第九种方案，在任一前述的MR成像法基础上提供一种MR成像法。这里，从均方根转换中产生的数据Av表达如下：
Av＝{(Fv(0)2+…+Fv(N-1)2)/N}
根据本发明第九种方案提供的MR成像法，N数据项都可采用，将不会被消除。因此，在许多情况下能够提高信噪比。
根据本发明的第十种方案，提供一种主要由发射机线圈，梯度线圈单元，接收机线圈，扫描设备和数据处理设备组成的MRI系统。发射机线圈发射射频(RF)脉冲。梯度线圈单元应用磁场梯度。接收机线圈接收NMR信号。扫描设备驱动发射机线圈，梯度线圈单元和接收机线圈以获取数据。数据处理设备对获取的数据执行算术操作以产生图像。这里，在从组成K空间的视图中获取从数据fv(0)到fv(N-1)范围中的数据fv(k)前，扫描设备通过顺序改变相位编码的相位对稳定状态自由进动(SSFP)中的数据重复采集N次(这里N表示2的幂)。此时，射频(RF)脉冲的相位根据一个表达式360°□v□k/N而变化。如果操作者指定傅里叶变换(FT)成像，数据处理设备就对数据fv(k)相应于用RF脉冲表示的每个相位进行傅里叶变换，因此产生数据Fv(n)。如果操作者不指定傅里叶变换成像，数据fv(k)实际上就被认为是Fv(n)。此后，为了产生数据Av，就对数据Fv(n)执行操作者选择的至少加权加法、最大强度投影(MIP)处理和均方根转换中的任一过程。因此，根据所产生的数据Av重建图像。
在根据本发明第十种方案提供的MRI系统中，能够较佳地实现本发明第一种方案中提供的MR成像法。
根据本发明的第十一种方案，提供一种主要由发射机线圈，梯度线圈单元，接收机线圈，扫描设备和数据处理设备组成的MRI系统。发射机线圈发射射频(RF)脉冲。梯度线圈单元应用磁场梯度。接收机线圈接收NMR信号。这里，在从组成K空间的视图中获取从数据fv(0)到fv(N-1)范围中的数据fv(k)前，扫描设备通过顺序改变相位编码的相位对SSFP中的数据重复采集N次(这里N表示2的幂)。此时，RF脉冲的相位根据一个表达式360°□v口k/N而变化。数据处理设备对数据fv(k)相应于用RF脉冲表示的每个相位进行傅里叶变换以产生数据Fv(n)。此后，为了产生数据Av，就对数据Fv(n)执行操作者选择的至少加权加法、最大强度投影(MIP)处理和均方根转换中的任一过程。之后，根据所产生的数据Av重建图像。
在根据本发明第十一种方案提供的MRI系统中，能够较佳地实现本发明第二种方案中提供的MR成像法。
根据本发明的第十二种方案，提供一种主要由发射机线圈，梯度线圈单元，接收机线圈，扫描设备和数据处理设备组成的MRI系统。发射机线圈发射射频(RF)脉冲。梯度线圈单元应用磁场梯度。接收机线圈接收NMR信号。这里，在从组成K空间的视图中获取从数据fv(0)到fv(N-1)范围中的数据fv(k)前，扫描设备通过顺序改变相位编码的相位对SSFP中的数据重复采集N次(这里N表示2的幂)。此时，RF脉冲的相位根据一个表达式360°□v□k/N而变化。数据处理设备认为数据fv(k)实际上就是数据Fv(n)。此后，为了产生数据Av，对数据Fv(n)执行操作者选择的至少加权加法、最大强度投影(MIP)处理和均方根转换中的任一过程。因此，根据所产生的数据Av重建图像。
在根据本发明第十二种方案提供的MRI系统中，能够较佳地实现本发明第三种方案中提供的MR成像法。
根据本发明的第十三种方案，在任一前述的MRI系统基础上提供一种MRI系统。一种用于获取SSFP中数据的脉冲序列能够启动FID信号和回波的同时采集。而且，每个磁场梯度表示随时间变化的强度，在一个重复时间TR过程中检测的其强度值被积分为0。
在根据本发明第十三种方案提供的MRI系统中，能够较佳地实现本发明第四种方案中提供的MR成像法。
根据本发明的第十四种方案，在任一前述的MRI系统基础上提供一种MRI系统。这里，从相应于用RF脉冲表示的每个相位执行傅里叶变换中产生的数据Fv(n)表达如下：
Fv ( n ) = Σ K = 0 N - 1 fv ( k ) . exp { - j . n . 2 π . k / M }
在根据本发明第十四种方案提供的MRI系统中，能够较佳地实现本发明第五种方案中提供的MR成像法。
根据本发明的第十五种方案，在任一前述的MRI系统基础上提供一种MRI系统。这里，从N＝2时执行加权加法中产生的数据Av表达如下：
Av＝0.5*Fv(0)+0.5*Fv(1)
在根据本发明第十五种方案提供的MRI系统中，能够较佳地实现本发明第六种方案中提供的MR成像法。
根据本发明的第十六种方案，在任一前述的MRI系统基础上提供一种MRI系统。这里，从N＝2时执行加权加法中产生的数据Av表达如下：
Av＝0.5*Fv(0)-0.5*Fv(1)
在根据本发明第十六种方案提供的MRI系统中，能够较佳地实现本发明第七种方案中提供的MR成像法。
根据本发明的第十七种方案，在任一前述的MRI系统基础上提供一种MRI系统。这里，从解决函数最大值{}以获得一最大值的MIP处理中产生的数据Av表达如下：
Av＝max{Fv(0)，.....，Fv(N-1)}
在根据本发明第十七种方案提供的MRI系统中，能够较佳地实现本发明第八种方案中提供的MR成像法。
根据本发明的第十八种方案，在任一前述的MRI系统基础上提供一种MRI系统。这里，从均方根转换中产生的数据Av表达如下：
Av＝{(Fv(0)2+…+Fv(N-1)2)/N}
在根据本发明第十八种方案提供的MRI系统中，能够较佳地实现本发明第九种方案中提供的MR成像法。
因此，在本发明中实现的MR成像法和MRI系统能够产生无伪象带的质量好的图像。
本发明的其它目的和优点从下面如附图所示本发明最佳实施例的描述中将会变得更清楚。
附图的简述
图1是根据本发明一个实施例的MRI系统的方框图。
图2是在本发明实施例中实现MR成像的流程图。
图3是FIESTA使用的脉冲序列的实例视图。
图4是表示K空间和视图的解释图。
图5表示数据采集次序实例和RF脉冲相位的实例的视图。
图6表示数据采集次序另一个实例和RF脉冲相位的另一个实例的视图。
本发明的详细描述
下面将进一步描述本发明有关附图中所示的一个实施例。
图1是根据本发明一个实施例的MRI系统的方框图。
在MRI系统100中，一磁组件1具有-物体插入其中的钻孔。一静磁场线圈1p，一梯度线圈单元1g，一发射机线圈1t，和一接收机线圈1r排列成似乎能包围钻孔。静磁场线圈1p将一定强度的静磁场作用于一物体。梯度线圈单元1g在一空间的X，Y和Z轴方向(定义为限幅(slice)选择，读出和相位编码的方向)产生磁场梯度。发射机线圈1t应用射频(RF)脉冲以激励物体的原子核自旋。接收机线圈1r检测物体感应的NMR信号。静磁场线圈1p，梯度线圈单元1g，发射机线圈1t，和接收机线圈1r分别与静磁场电源2，梯度驱动电路3，RF功率放大器4和前端放大器5相连接。
顺便提及，一永磁铁可代替静磁场线圈1p。
根据一存储在时序存储器电路6中的脉冲序列，相应于计算机7发出的指令可致动梯度驱动电路3。这就使包含在磁组件1中的梯度线圈单元1g产生磁场梯度。同时，一门控调制电路8将由RF振荡电路9提供的载波输出信号调制为表示设定定时和用一设定包线检测的脉冲信号。该合成信号作为RF脉冲传输给RF功率放大器4。该功率放大器4对RF脉冲的功率进行放大。然后生成的脉冲作用于包含在磁组件1中的发射机线圈1t上。之后，将被成像的所需限幅(slice)就能被选择和激励。
前端放大器5放大物体感应和经过包含在磁组件1中的接收机线圈1r检测的NMR信号，并将该生成信号传输给相位检测器10。相位检测器10使用RF振荡电路9的载波输出信号作为参考信号来检测前端放大器5输送的NMR信号的相位，然后将该信号传输给A/D转换器11。A/D转换器11将相位经过检测的模拟信号转换为数字信号，并将该数字信号传输给计算机7。
计算机7接收于操作者控制台12输入的信息，并负责控制整个系统。而且，计算机7从A/D转换器中读取数字数据，对该数据执行算术操作，从而产生图像。
之后，图像显示在显示设备13上或信息呈现在显示设备13上。
图2是描述用MRI系统100实现MR成像的流程图。
在步骤S1中，重复计数器k初始化为0。
在步骤S2中，在从组成K空间的视图v中获得数据fv(k)前通过顺序改变相位编码的相位，稳定状态自由进动(SSFP)中的数据被重复采集。此时，射频(RF)脉冲的相位根据一个表达式360°口v□k/N而变化。
在步骤S3和S4中，步骤S2的处理从计数器k表示为1的时刻到计数器表示为N-1的时刻重复N-1次。
为了简短起见，重复的次数N应该为8。此时，RF脉冲的相位根据表达式v.k.π/4而变化。
图3表示在SSFP中数据进行采集使用的其中一个脉冲序列的实例，即，使用稳定状态采集(FIESTA)快速成像所使用的脉冲序列。
根据FIESTA脉冲序列，RF脉冲在每隔一比将被测量物体显示的横弛豫时间短的重复时间TR内被重复使用。为了获取数据fv(k)，应检测SSFP中产生的自由感应衰减(FID)信号和回波(自旋回波或激励回波)。将被采用的磁场梯度是这样的：每个磁场梯度显示出随时间变化的强度，其在一个重复时间TR过程中检测的强度值被积分为0。而且，相位编码的磁场强度顺序相应于每个视图V改变强度。RF脉冲的相位相应于每个视图v和重复计数器k的值被顺序改变。
图4是表示一k空间和一k空间被扫描以获取数据所沿的轨迹的总体图。
K空间是一用读出方向和相位编码方向定义的二维空间。
这里，假定视图V#0至#31沿相位编码磁场梯度方向并置。
假定k空间被分为由视图#0至#7组成的第一磁场，由视图#8至#15组成的第二磁场，由#视图#16至#23组成的第三磁场和由视图#24至#31组成的第四磁场。
图5表示k＝0的数据采集次序的一个实例和RF脉冲相位的一个实例。
首先，数据没有被获取，但是为了获得SSFP重复应用FIESTA脉冲序列。这种应用应该称之为无效果应用。
在SSEP中，为了从组成第一磁场的视图#0至#7中顺序获取数据f0(0)至数据f7(0)就应用FIESTA脉冲。此时，对所有应用RF脉冲的相位设定为0。
此后，SSFP通过无效果应用获得。然后从组成第二磁场的视图#8至#15中顺序获取数据f8(0)至数据f15(0)。此时，对所有应用RF脉冲的相位设定为0。
此后，SSFP通过无效果应用获得。然后从组成第三磁场的视图#16至#23中顺序获取数据f16(0)至数据f23(0)。此时，对所有应用RF脉冲的相位设定为0。
此后，SSFP通过无效果应用获得。然后从组成第四磁场的视图#24至#31中顺序获取数据f24(0)至数据f31(0)。此时，对所有应用RF脉冲的相位设定为0。
图6表示数据采集次序k＝1的一个实例和RF脉冲相位的一个实例的视图。
首先，SSFP通过无效果应用获得。然后从组成第一磁场的视图#0至#7中顺序获取数据f0(1)至数据f7(1)。此时，RF脉冲的相位依次变化为0，π/4，π/2，3π/4，π，5π/4，3π/2，和7π/4。
此后，SSFP通过无效果应用获得。然后从组成第二磁场的视图#8至#15中顺序获取数据f8(1)至数据f15(1)。此时，RF脉冲的相位依次变化为0，π/4，π/2，3π/4，π，5π/4，3π/2，和7π/4。
此后，SSFP通过无效果应用获得。然后从组成第三磁场的视图#16至#23中顺序获取数据f16(1)至数据f23(1)。此时，RF脉冲的相位依次变化为0，π/4，π/2，3π/4，π，5π/4，3π/2，和7π/4。
此后，SSFP通过无效果应用获得。然后从组成第四磁场的视图#24至#31中顺序获取数据f24(1)至数据f31(1)。此时，RF脉冲的相位依次变化为0，π/4，π/2，3π/4，π，5π/4，3π/2，和7π/4。
回来参考图2，在步骤S5中，如果操作者指定傅里叶变换成像，控制就转到步骤S6。否则，控制就转到步骤S7。
在步骤S6中，为了产生数据Fv(n)，对数据fv(k)相应于用RF脉冲表示每个相位进行傅里叶变换。数据Fv(n)表达如下：
Fv ( n ) = Σ k = 0 7 fv ( k ) . exp { - j . n . π . k / 4 }
然后控制转到步骤S8。
在步骤S7中，数据fv(k)实际上被认为是数据Fv(n)。然后控制就转到步骤S8。
在步骤S8中，如果操作者指定均方根转换成像，控制就转到步骤S9。否则，控制就转到步骤S10。
在步骤S9中，为了产生如下表达的数据Av，对数据Fv(n)进行均方根转换：
Av = { ( Fv ( 0 ) 2 + . . . . . + Fv ( 7 ) 2 / 8 )
然后控制就转到步骤S11。
在步骤S10中，为了产生数据Av，对数据Fv(n)执行操作者选择的加权加法或最大强度投影(MIP)处理。加权加法产生的数据Av表达如下：
Av＝C0.Fv(0)+C1.Fv(1)+...+C7.Fv(7)
这里C0，C1，....C7表示权重。
顺便提及，N＝2时，C0＝C1＝0.5或C0＝0.5和C1＝-0.5。
从MIP处理中产生的数据Av表达如下：
Av＝max{Fv(0)，...，Fv(7)}
这里max{}表示提供最大值的函数。
然后控制就转到步骤S11。
图像根据数据Av进行重建。然后终止处理。
根据前述MRI系统100，操作者能够指定是否应该对数据fv(k)相应于用RF脉冲表示的每个相位进行傅里叶变换。而且，操作者能够从至少加权加法、MIP和均方根转换中选择任一个对数据Fv(n)执行处理。这样，能够在任何常用MR成像法仅能产生质量差的图像的情况下产生质量好的图像。
在前述实施例中，k空间被划分为四个场。相反，k空间可不必划分或可划分为五个或更多场。
本发明可由许多广泛不同的实施例构成而不脱离本发明的精神和变化范围。应该明白除了后面权利要求限定的外，本发明并不局限于本说明书中描述的这些
具体实施例。

Claims (13)

1.一种磁共振成像系统，包括：
用于发射射频脉冲的发射机线圈；
用于施加磁场梯度的梯度线圈单元；
用于接收核磁共振信号的接收机线圈；
一个扫描设备，用于：
通过驱动上述发射机线圈、梯度线圈单元和接收机线圈来获取数据，
通过顺序改变相位编码的相位，在稳定状态自由进动中重复信息采集N次，直到从视图v获取在从数据fv(0)到fv(N-1)范围中的第一数据组fv(k)，这里N表示2的幂，且在上述信息采集的同时，根据一个表达式360°vk/N来变化射频脉冲的相位；和一个数据处理设备，用于：
对所述信息执行算术操作；
如果操作者指定傅里叶变换成像，上述数据处理设备就对上述第一数据组fv(k)相应于用射频脉冲表示的每个相位进行傅里叶变换以产生第二数据组Fv(n)；如果操作者不指定傅里叶变换成像，则使用所述第一数据组fv(k)作为第三数据组Fv1(n)；从操作者接收一个选择：是否对所述第二数据组和第三数据组之一执行均方根转换；
如果操作者选择了所述均方根转换，则对所述第二数据组和第三数据组之一执行均方根转换，以产生一个第四数据组；和
如果操作者没有选择所述均方根转换，则对所述第二数据组和第三数据组之一执行加权加法和最大强度投影处理之一，以产生一个第五数据组；和
从上述第四数据组和第五数据组之一重建一个图像。
2.一种磁共振成像系统，包括：
用于发射射频脉冲的发射机线圈；
用于施加磁场梯度的梯度线圈单元；
用于接收核磁共振信号的接收机线圈；
一个扫描设备，用于：
通过驱动上述发射机线圈、梯度线圈单元和接收机线圈来获取数据，通过顺序改变相位编码的相位在稳定状态自由进动中重复信息采集N次，直到从视图v获取在从数据fv(0)到fv(N-1)范围中的第一数据组fv(k)，这里N表示2的幂，在所述信息采集的同时，根据一个表达式360°vk/N变化射频脉冲的相位；和一个数据处理设备，用于：
对所述信息执行算术操作；
通过对第一数据组fv(k)相应于用射频脉冲表示的每个相位进行傅里叶变换以产生第二数据组Fv(n)，从操作者接收一个是否对所述第二数据组执行均方根转换的选择；
如果操作者选择了均方根转换，则对所述第二数据组执行均方根转换以产生一个第三数据组；和
如果操作者未选择均方根转换，则对所述第二数据组执行加权加法和最大强度投影之一以产生一个第四数据组；和
根据所述第三数据组和第四数据组之一重建一个图像。
3.一种磁共振成像系统，包括：
用于发射射频脉冲的发射机线圈；
用于施加磁场梯度的梯度线圈单元；
用于接收核磁共振信号的接收机线圈；
一个扫描设备，用于：
通过驱动上述发射机线圈、梯度线圈单元和接收机线圈来获取数据；
通过顺序改变相位编码的相位，在稳定状态自由进动中重复信息采集N次，直到在从视图v中获取在数据fv(0)到fv(N-1)范围中的第一数据组fv(k)，这里N表示2的幂，在所述信息采集的同时，根据一个表达式360°vk/N变化射频脉冲的相位；
一个数据处理设备，用于：
对所述信息执行算术操作；
如果操作者没有指定傅里叶变换成象，则将所述第一数据组fv(k)作为第二数据组数据Fv(n)；
从操作者接收一个是否对所述第二数据组执行均方根转换的选择；
如果操作者选择了均方根转换，则对所述第二数据组执行均方根转换以产生一个第三数据组；和
如果操作者未选择均方根转换，则对所述第二数据组执行加权加法和最大强度投影之一以产生一个第四数据组；和
根据所述第三数据组和第四数据组之一重建一个图像。
4.根据权利要求1至3任一项所述的磁共振成像系统，其特征在于：用于在稳定状态自由进动中采集数据的一个脉冲序列用于启动一个自由感应衰减信号和一个回波的同时采集，且所采用的磁场梯度是：每个磁场梯度显示出随时间变化的强度，在一个重复时间的过程中检测的强度值被积分为0。
5.根据权利要求1或2所述的磁共振成像系统，其特征在于：从相应于用射频脉冲表示的每个相位所执行的傅里叶变换中产生的第二数据组Fv(n)表达如下：
Fv ( n ) = Σ k = 0 N - 1 fv ( k ) . exp { - j . n . 2 π . k / N } .
6.根据权利要求1所述的磁共振成像系统，其特征在于：从N＝2时执行的加权加法中产生的第五数据组Av表达如下：
Av＝0.5×Fv(0)+0.5×Fv(1)。
7.根据权利要求1所述的磁共振成像系统，其特征在于：从N＝2时执行的加权加法中产生的第五数据组Av表达如下：
Av＝0.5×Fv(0)-0.5×Fv(1)。
8.根据权利要求1所述的磁共振成像系统，其特征在于：从解函数最大值以提供一最大值的最大强度投影处理中产生的第五数据组Av表达如下：
Av＝max{Fv(0)，.....，Fv(N-1)}。
9.根据权利要求1所述的磁共振成像系统，其特征在于：从均方根转换中产生的第五数据组Av表达如下：
Av＝{(Fv(0)2+…+Fv(N-1)2)/N}。
10.根据权利要求2或3所述的磁共振成像系统，其特征在于：从N＝2时执行的加权加法中产生的第四数据组Av表达如下：
Av＝0.5×Fv(0)+0.5×Fv(1)。
11.根据权利要求2或3所述的磁共振成像系统，其特征在于：从N＝2时执行的加权加法中产生的第四数据组Av表达如下：
Av＝0.5×Fv(0)-0.5×Fv(1)。
12.根据权利要求2或3所述的磁共振成像系统，其特征在于：从解函数最大值以提供一最大值的最大强度投影处理中产生的第四数据组Av表达如下：
Av＝max{Fv(0)，.....，Fv(N-1)}。
13.根据权利要求2或3所述的磁共振成像系统，其特征在于：从均方根转换中产生的第三数据组Av表达如下：
Av＝{(Fv(0)2+…+Fv(N-1)2)/N}。
CNB021272611A 2001-06-21 2002-06-21 磁共振成像系统 Expired - Fee Related CN1247151C (zh)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
JP187956/2001 2001-06-21
JP2001187956A JP3891799B2 (ja) 2001-06-21 2001-06-21 Ｍｒｉ装置
JP187956/01 2001-06-21

Publications (2)

Publication Number Publication Date
CN1395904A CN1395904A (zh) 2003-02-12
CN1247151C true CN1247151C (zh) 2006-03-29

Family

ID=19027165

Family Applications (1)

Application Number Title Priority Date Filing Date
CNB021272611A Expired - Fee Related CN1247151C (zh) 2001-06-21 2002-06-21 磁共振成像系统

Country Status (4)

Country Link
US (1) US6819104B2 (zh)
JP (1) JP3891799B2 (zh)
KR (1) KR100459099B1 (zh)
CN (1) CN1247151C (zh)

Families Citing this family (21)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7558612B2 (en) * 2002-04-16 2009-07-07 The Board Of Trustees Of The Leland Stanford Junior University Motion compensated spiral FISP MRI
US7009395B2 (en) * 2002-05-15 2006-03-07 General Electric Company Method for MRI with steady-state free precession pulse sequences
DE10353342B4 (de) * 2003-11-14 2008-07-17 Siemens Ag Verbesserte MRT-Bildgebung auf Basis konventioneller PPA-Rekonstruktionsverfahren
US7587231B2 (en) * 2004-01-09 2009-09-08 Toshiba America Mri, Inc. Water fat separated magnetic resonance imaging method and system using steady-state free-precession
JP3968352B2 (ja) * 2004-02-03 2007-08-29 ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー Ｍｒｉ装置
JP4192139B2 (ja) * 2004-10-27 2008-12-03 ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー Ｓｓｆｐパルスシーケンスにおけるフェーズサイクリング法および磁気共鳴撮影装置
US7132828B2 (en) * 2005-03-28 2006-11-07 The Board Of Trustees Of The Leland Stanford Junior University Artifact reduction in SSFP MRI using super field view reconstruction
JP4327171B2 (ja) * 2006-04-06 2009-09-09 ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー 最大値投影方法および装置
JP4262737B2 (ja) 2006-10-06 2009-05-13 ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー 磁気共鳴イメージング装置、スキャン装置、プログラムおよび記憶媒体
CN101174331B (zh) * 2006-11-01 2011-07-27 深圳市蓝韵实业有限公司 一种医学影像最大密度投影生成方法
EP1959397B1 (en) * 2007-02-19 2019-08-07 Wisconsin Alumni Research Foundation Iterative HYPR medical image reconstruction
CN101281241B (zh) * 2007-04-06 2012-06-20 Ge医疗系统环球技术有限公司 Mri设备及其控制方法
JP4585539B2 (ja) * 2007-04-12 2010-11-24 ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー Ｍｒｉ装置
JP5068606B2 (ja) * 2007-08-29 2012-11-07 ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー 磁気共鳴イメージング装置，プログラム
JP5361311B2 (ja) * 2007-09-26 2013-12-04 シーメンス アクチエンゲゼルシヤフト 高周波送信装置、高周波送信アンテナ装置、磁気共鳴システムおよび磁気共鳴画像を生成する方法
WO2014153775A1 (zh) * 2013-03-29 2014-10-02 深圳先进技术研究院 一种快速磁共振成像方法和系统
CN103278786B (zh) * 2013-03-29 2015-08-19 深圳先进技术研究院 一种快速磁共振成像方法和系统
KR101650891B1 (ko) * 2015-07-09 2016-08-24 성균관대학교산학협력단 자기 공명 전기 임피던스 영상기술을 이용한 자기 공명 영상시스템 및 전도율 분포 영상 생성방법
CN107192970A (zh) * 2017-05-05 2017-09-22 中国科学院电工研究所 一种磁共振成像系统的阵列梯度线圈驱动装置
CN107861078A (zh) * 2017-12-15 2018-03-30 冯友珍 一种mri影像诊断系统
CN108152769B (zh) * 2017-12-22 2019-08-16 中国科学院武汉物理与数学研究所 一种气体波谱的角度和弛豫时间常数t1同时测量方法

Family Cites Families (31)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4639671A (en) * 1983-05-31 1987-01-27 General Electric Company Simultaneous NMR imaging system
US4766382A (en) * 1984-05-17 1988-08-23 Jeol Ltd. Two-dimensional nuclear magnetic resonance spectrometry
US4752734A (en) * 1985-08-06 1988-06-21 The General Hospital Corporation Flow imaging by means of nuclear magnetic resonance
US4728933A (en) * 1986-03-24 1988-03-01 U.S. Holding Company, Inc. Miniature dual cavity ringer
US4706024A (en) * 1986-05-05 1987-11-10 General Electric Company Suppression of non-coupled spins in NMR imaging and spectroscopy
US4800889A (en) * 1987-04-06 1989-01-31 General Electric Company Rapid-scan NMR angiography
US4912413A (en) * 1987-08-14 1990-03-27 Picker International, Inc. Conjugate symmetry magnetic resonance imaging
US4851779A (en) * 1987-08-14 1989-07-25 Picker International, Inc. Three dimensional volume imaging with conjugate symmetrization
US4780675A (en) * 1987-08-14 1988-10-25 Picker International, Inc. Conjugate symmetry magnetic resonance imaging
US4849697A (en) * 1988-06-27 1989-07-18 General Electric Company Three-dimensional magnetic resonance flow-contrast angiography with suppression of stationary material
US4918386A (en) * 1988-12-23 1990-04-17 General Electric Company Method for simultaneously obtaining three-dimensional NMR angiograms and stationary tissue NMR images
DE4004185C2 (de) 1989-02-24 1997-08-07 Siemens Ag Verfahren zur Gewinnung von flußkompensierten, T¶2¶- gewichteten Bildern mittels der kernmagnetischen Resonanz
JPH03268742A (ja) * 1990-03-19 1991-11-29 Hitachi Ltd イメージング装置
JPH05184555A (ja) * 1991-03-20 1993-07-27 Hitachi Ltd 磁気共鳴イメージング方法及び装置
US5347216A (en) * 1992-06-23 1994-09-13 General Electric Company Fast NMR image acquisition with spectrally selective inversion pulse
US5256967A (en) * 1992-10-01 1993-10-26 General Electric Company Fast NMR image acquisition with spectrally selective inversion pulses
US5729140A (en) * 1996-10-04 1998-03-17 Mayo Foundation For Medical Education And Research Suppression of artifacts in NMR images by correlation of two NEX acquistions
US6577128B1 (en) * 1998-10-15 2003-06-10 Btg International Limited NQR method and apparatus for testing a sample by applying multiple excitation blocks with different delay times
US6392408B1 (en) * 1998-05-06 2002-05-21 Quamtum Magnetics, Inc. Method and system for cancellation of extraneous signals in nuclear quadrupole resonance spectroscopy
US6043656A (en) * 1998-11-23 2000-03-28 General Electric Company Method for compensating an MRI system for residual magnetization
US6289233B1 (en) * 1998-11-25 2001-09-11 General Electric Company High speed tracking of interventional devices using an MRI system
DE19901763B4 (de) * 1999-01-18 2005-12-01 Siemens Ag Impulssequenz für ein Kernspintomographiegerät
US6275038B1 (en) * 1999-03-10 2001-08-14 Paul R. Harvey Real time magnetic field mapping using MRI
US6307368B1 (en) * 1999-05-14 2001-10-23 Board Of Trustees Of The Leland Stanford Junior University Linear combination steady-state free precession MRI
US6235264B1 (en) * 1999-12-01 2001-05-22 General Electric Company Medical imaging method for characterizing tumor angiogenesis using polymeric contrast agents
US6292684B1 (en) * 1999-12-14 2001-09-18 Ge Medical Systems Global Technology Company, Llc Respiratory displacement and velocity measurement using navigator MRI echo signals
KR100375922B1 (ko) * 2000-01-25 2003-03-15 주식회사 메디슨 Ｍｒａ의 영상구현시 발생하는 아티펙트제거방법
US6493569B2 (en) * 2000-12-30 2002-12-10 Ge Medical Systems Global Technology Company, Llc Method and apparatus using post contrast-enhanced steady-state free precession in MR imaging
JP3814157B2 (ja) * 2001-04-17 2006-08-23 ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー Ｍｒｉ装置
US6552542B1 (en) * 2001-09-28 2003-04-22 The Board Of Trustees Of The Leland Stanford Junior University Oscillating dual-equilibrium steady state angiography using magnetic resonance imaging
JP3796446B2 (ja) * 2002-01-23 2006-07-12 ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー Ｍｒｉ装置

Also Published As

Publication number Publication date
JP3891799B2 (ja) 2007-03-14
CN1395904A (zh) 2003-02-12
US20020196019A1 (en) 2002-12-26
KR20020097017A (ko) 2002-12-31
US6819104B2 (en) 2004-11-16
KR100459099B1 (ko) 2004-12-03
JP2003010148A (ja) 2003-01-14

Similar Documents

Publication Publication Date Title
CN1247151C (zh) 磁共振成像系统
CN1085342C (zh) 磁共振成像方法和装置
CN1251644C (zh) 多片层磁共振成象方法和装置
CN1745705A (zh) 磁共振影像生成方法和磁共振成像装置
CN102621510B (zh) 用于抑制mr成像中伪影的系统
CN1130570C (zh) 磁共振成象方法和装置
CN1851494A (zh) 改进发送侧加速的体积选择磁共振成像的方法和设备
CN101051074A (zh) 磁共振成像装置和磁共振成像方法
CN1265762C (zh) 线圈灵敏图生成方法和平行成像方法以及磁共振成像装置
CN1943510A (zh) 磁共振成像设备
CN1289920A (zh) 核磁共振成像装置
CN1765320A (zh) Ssfp脉冲序列中的相位循环方法和磁共振成像设备
CN1683939A (zh) 在磁共振断层造影中融合加权成像的运动校正的多点方法
CN1217622C (zh) 减少快速自旋回波磁共振图象中麦克斯韦项假象的方法
CN1378816A (zh) 相位校正方法和磁共振成像系统
CN101034151A (zh) 基于梯度回波序列进行磁共振成像的方法和设备
CN1356091A (zh) 图像处理方法与装置,记录介质及成像装置
CN1989900A (zh) Rf脉冲施加方法和mri装置
CN1769921A (zh) 在磁共振断层造影中避免皱褶伪影的优化方法
CN1237487C (zh) 图像处理方法和装置、成像装置
CN1441258A (zh) Mri设备和mra成像方法
CN1939213A (zh) 磁共振扫描方法和磁共振成像设备
CN1192743C (zh) Mr成像方法,相位误差测量方法和mri系统
JP3668076B2 (ja) Ｍｒｉ装置
CN1262238C (zh) 磁共振信号采集的方法和装置及磁共振成像装置

Legal Events

Date Code Title Description
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
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

Granted publication date: 20060329

Termination date: 20140621

EXPY Termination of patent right or utility model