[go: up one dir, main page]
More Web Proxy on the site http://driver.im/

CN1196930C - 相位误差测量方法、校正方法和装置、磁共振成像装置 - Google Patents

相位误差测量方法、校正方法和装置、磁共振成像装置 Download PDF

Info

Publication number
CN1196930C
CN1196930C CNB011357940A CN01135794A CN1196930C CN 1196930 C CN1196930 C CN 1196930C CN B011357940 A CNB011357940 A CN B011357940A CN 01135794 A CN01135794 A CN 01135794A CN 1196930 C CN1196930 C CN 1196930C
Authority
CN
China
Prior art keywords
phase
phase error
excitations
gradient magnetic
period
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
CNB011357940A
Other languages
English (en)
Other versions
CN1350176A (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 CN1350176A publication Critical patent/CN1350176A/zh
Application granted granted Critical
Publication of CN1196930C publication Critical patent/CN1196930C/zh
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Images

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/05Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves 
    • A61B5/055Detecting, 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
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00Arrangements or instruments for measuring magnetic variables
    • G01R33/20Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48NMR imaging systems
    • G01R33/54Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
    • G01R33/56Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
    • G01R33/565Correction of image distortions, e.g. due to magnetic field inhomogeneities
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00Arrangements or instruments for measuring magnetic variables
    • G01R33/20Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48NMR imaging systems
    • G01R33/54Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
    • G01R33/56Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
    • G01R33/565Correction of image distortions, e.g. due to magnetic field inhomogeneities
    • G01R33/56518Correction of image distortions, e.g. due to magnetic field inhomogeneities due to eddy currents, e.g. caused by switching of the gradient magnetic field
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00Arrangements or instruments for measuring magnetic variables
    • G01R33/20Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48NMR imaging systems
    • G01R33/54Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
    • G01R33/56Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
    • G01R33/565Correction of image distortions, e.g. due to magnetic field inhomogeneities
    • G01R33/56563Correction of image distortions, e.g. due to magnetic field inhomogeneities caused by a distortion of the main magnetic field B0, e.g. temporal variation of the magnitude or spatial inhomogeneity of B0

Landscapes

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

Abstract

为了精确地测量和校正相位轴方向上自旋的相位误差,在从第一个180°激励到第二个180°激励之间的时段期间在相位轴方向上施加一个具有零的积分值的梯度磁场,以读出第一自旋回波SE1;在从第二个180°激励到第三个180°激励之间的时段期间在相位轴方向上施加一个具有零的积分值的梯度磁场,以读出第二自旋回波信号;根据这些自旋回波确定由于残磁的效应而产生的相位误差。

Description

相位误差测量方法、校正方法和装置、磁共振成像装置
发明领域
本发明涉及一种相位误差测量方法和装置、相位误差校正方法和装置、记录媒体和磁共振成像装置,更具体来说,涉及一种用于测量因残磁(residua1 magnetization)等等的效应而引起的自旋中的相位误差的方法和装置、一种用于校正这种误差的方法和装置、一种记录有用于操作计算机来执行一个相位误差测量功能的程序的记录媒体,一种记录有用于操作计算机来执行一个相位误差校正功能的程序的记录媒体,一种包含相位误差测量装置的磁共振成像装置,和一种包含相位误差校正装置的磁共振成像装置。
背景技术
在磁共振成像(MRI)装置中,将要成像的物体装入磁系统的内部空间中,即其中生成静态磁场的空间中;施加梯度磁场和高频磁场,以从该物体内的自旋中生成磁共振信号;然后根据所接收信号重构一个层析X射线摄影图像。
梯度磁场被施加在三个互相垂直的轴向上。这三个互相垂直的轴是限幅(slice)轴、相位轴和频率轴。限幅轴方向上的梯度磁场的作用是选择性地由一个RF(无线频率)激励信号激励限幅轴上的一个期望的限幅,有时被称作限幅梯度(slice gradient)。相位轴方向上的梯度磁场的用作自旋的相位编码,有时被称作相位编码梯度(phaseencoding gradient)。频率轴方向上的梯度磁场的用作磁共振信号的读出,有时被称作读出梯度。磁共振信号是以回波信号的形式被读出的。
一种磁共振成像方法是快速自旋回波(FSE)技术。这个技术要激励90°的自旋,接着多次地重复180°激励的自旋的倒置,以获得每个90°激励的多个视象(views)的自旋回波。
图1中表示了FSE技术的一个脉冲序列。在图1中,(1)是90°和180°激励的序列,(2)、(3)和(4)分别是限幅梯度序列Gs,相位编码梯度Gp和读出梯度Gr;(5)是自旋回波SE序列。这些序列沿一个时间轴t发生。
如图所示,在施加一个限幅梯度Gs1的同时,实施一个90°激励。然后,在过了时段U1后,在施加一个限幅梯度Gs2的同时,实施第一个180°激励。然后,在过了时段U2后,在施加一个限幅梯度Gs3的同时,实施第二个180°激励。自此之后,以类似的方式,在每个时段U2,分别在施加限幅梯度Gs4、Gs5、…的同时,实施第三个、第四个、…180°激励。
在从90°激励到第一个180°激励的时段期间,施加一个读出梯度Gr1,以进行自旋的相位分散(phase dispersion,即dephasing,下文中有时称“移相”)。然后,在第一个180°激励到第二个180°激励的时段期间,施加一个读出梯度Gr2,以进行自旋的相位聚焦(phasefocusing,即rephasing,下文中有时称“再定相”),生成第一个自旋回波SE1。生成了该自旋回波的读出梯度Gr2在其后半部分分散自旋的相位(dephases the spins)。
在读出梯度Gr2的施加之前,施加一个相位编码梯度Gp1,以进行相位编码,在读出梯度Gr2的施加被完成后,在相反方向上施加一个相位编码梯度Gp1’,以消除该相位编码。
自此之后,以类似的方式在180°激励之间的每个时段期间施加读出梯度Gr3、Gr4、…,以生成各自的自旋回波SE2、SE3、…。此外,通过相位编码梯度Gp2、Gp3、…获得相位编码。使相位编码每一次都有差异。
自旋回波是一个在回波中心有最大幅值的RF信号。第一个自旋回波SE1的最大幅值或峰值在从90°激励起的一个时段TE(回波时间)后发生。第二个自旋回波SE2的峰值在从第一个自旋回波SE1的峰值起的时段TE后发生。自此之后,自旋回波SE3、SE4、…的峰值类似地按时段TE的间隔发生。峰值的生成有时被称作自旋回波的聚焦(图像形成)。
在通过永久磁体实现静态磁场生成的磁系统中,由于例如永久磁体的极片被梯度磁场的磁化,会发生残磁。由于频率轴方向上的残磁影响在90°激励与180°激励之间的时段的自旋的移相,自旋回波聚焦的定时或峰值生成的定时(timing),受到误差。类似的现象是由涡流产生的。
例如,如果第一个自旋回波SE1在比合适的TE短的时段TE’内聚焦,则由于聚焦中的定时误差,第二个自旋回波SE2在比合适的TE长的时段TE”内聚焦,第三个自旋回波SE3在比合适的TE短的时段TE’内聚焦。
由于这种聚焦误差是重构的图像中噪声因素生成的根源,人们曾试图通过调整读出梯度来消除频率轴方向上残磁、涡流等的效应,以校正聚焦误差。
当在相位轴方向上也存在残磁时,由于该残磁而发生聚焦误差。然而该聚焦误差却不能通过调整读出梯度的方法而校正,因为梯度的轴是不同的。此外,相位轴方向上残磁对自旋的相位的效应,本来就不可能准确地测量。
发明内容
因此,本发明的一个目的是提供一种用于精确地测量相位轴方向上的自旋中的相位误差的方法和装置、用于校正这种相位误差的方法和装置、一种记录有用于操作计算机来执行一个相位误差测量功能的程序的记录媒体,一种记录有用于操作计算机来执行一个相位误差校正功能的程序的记录媒体,一个包含相位误差测量装置的磁共振成像装置,和一个包含相位误差校正装置的磁共振成像装置。
(1)按照用于解决前述问题的一个方面,本发明是一种相位误差测量方法,特征在于包含:对物体自旋实施一个90°激励;在从所述90°激励起的第一时段过后实施第一个180°激励;在从所述第一个180°激励起的第二时段过后实施第二个180°激励;在从所述第二个180°激励起的第二时段过后实施第三个180°激励;在从所述第一个180°激励到所述第二个180°激励的时段期间在一个相位轴方向上施加一个零的时间积分值的梯度磁场,以读出一个第一自旋回波信号;在从所述第二个180°激励到所述第三个180°激励的时段期间在该相位轴方向上施加一个零的时间积分值的梯度磁场,以读出一个第二自旋回波信号;和根据所述第一和第二自旋回波信号确定在从所述90°激励到所述第一个180°激励的时段期间在该相位轴方向上自旋的相位误差。
(2)按照用于解决前述问题的另一个方面,本发明是一种相位误差测量装置,特征在于包含:自旋激励设备,用于对物体自旋实施一个90°激励,在从所述90°激励起的第一时段过后实施第一个180°激励,在从所述第一个180°激励起的第二时段过后实施第二个180°激励,以及在从所述第二个180°激励起的第二时段过后实施第三个180°激励;回波读出设备,用于在从所述第一个180°激励到所述第二个180°激励的时段期间在一个相位轴方向上施加一个零的时间积分值的梯度磁场,以读出一个第一自旋回波信号,以及在从所述第二个180°激励到所述第三个180°激励的时段期间在该相位轴方向上施加一个零的时间积分值的梯度磁场,以读出一个第二自旋回波信号;和相位误差计算设备,用于根据所述第一和第二自旋回波信号确定在从所述90°激励到所述第一个180°激励的时段期间在该相位轴方向上自旋的相位误差。
(3)按照用于解决前述问题的另一个方面,本发明是一种记录媒体,特征在于以计算机可读的方式记录一个程序,用于使计算机执行下述功能:对物体自旋实施一个90°激励;在从所述90°激励起的第一时段过后实施第一个180°激励;在从所述第一个180°激励起的第二时段过后实施第二个180°激励;在从所述第二个180°激励起的第二时段过后实施第三个180°激励;在从所述第一个180°激励到所述第二个180°激励的时段期间在一个相位轴方向上施加一个零的时间积分值的梯度磁场,以读出一个第一自旋回波信号;在从所述第二个180°激励到所述第三个180°激励的时段期间在该相位轴方向上施加一个零的时间积分值的梯度磁场,以读出一个第二自旋回波信号;和根据所述第一和第二自旋回波信号确定在从所述90°激励到所述第一个180°激励的时段期间在该相位轴方向上自旋的相位误差。
(4)按照用于解决前述问题的另一个方面,本发明是一种磁共振成像装置,用于用一个静态磁场、梯度磁场和一个高频磁场获得来自一个物体的磁共振信号,并根据所获得的所述磁共振信号构造一个图像,特征在于包含:自旋激励设备,用于对物体自旋实施一个90°激励,在从所述90°激励起的第一时段过后实施第一个180°激励,在从所述第一个180°激励起的第二时段过后实施第二个180°激励,以及在从所述第二个180°激励起的第二时段过后实施第三个180°激励;回波读出设备,用于在从所述第一个180°激励到所述第二个180°激励的时段期间在一个相位轴方向上施加一个零的时间积分值的梯度磁场,以读出一个第一自旋回波信号,以及在从所述第二个180°激励到所述第三个180°激励的时段期间在该相位轴方向上施加一个零的时间积分值的梯度磁场,以读出一个第二自旋回波信号;和相位误差计算设备,用于相据所述第一和第二自旋回波信号确定在从所述90°激励到所述第一个180°激励的时段期间在该相位轴方向上自旋的相位误差。
按照如就(1)-(4)所述的各方面的本发明,由于在90°与180°激励之间的时段期间在相位轴方向上没有施加梯度磁场,并且在各180°激励之间的时段期间在相位轴方向上施加各自具有零的时间积分值的梯度磁场,以读出各自的自旋回波信号,所以能获得两个自旋回波信号,这两个信号受到在90°与180°激励之间的时段期间在相位轴方向上的残磁的互相相反的影响。因此,通过计算就能从这些自旋回波信号确定自旋的一个相位误差。
在如就(1)-(4)所述的各方面的本发明中,通过在每个180°激励之前和之后在频率轴方向上施加分相器(crusher)梯度磁场,就能消除由FID信号产生的干扰,并且也能防止受激回波的生成,以提高测量的精确性。
(5)按照用于解决前述问题的另一个方面,本发明是一种相位误差测量方法,特征在于包含:对物体自旋实施一个90°激励;在从所述90°激励起的第一时段过后实施第一个180°激励;在从所述第一个180°激励起的第二时段过后实施第二个180°激励;在从所述第二个180°激励起的第二时段过后实施第三个180°激励;在从所述第一个180°激励到所述第二个180°激励的时段期间在一个相位轴方向上施加一个零的时间积分值的梯度磁场,以读出一个第一自旋回波信号;在从所述第二个180°激励到所述第三个180°激励的时段期间在该相位轴方向上施加一个零的时间积分值的梯度磁场,以读出一个第二自旋回波信号;根据所述第一和第二自旋回波信号确定在从所述90°激励到所述第一个180°激励的时段期间在该相位轴方向上自旋的相位误差;试探在所述梯度磁场的极性颠倒的情况下重复两次所述的步骤;和确定前者所确定的相位误差与后者所确定的相位误差的平均值。
(6)按照用于解决前述问题的另一个方面,本发明是一种相位误差测量装置,特征在于包含:自旋激励设备,用于对物体自旋实施一个90°激励,在从所述90°激励起的第一时段过后实施第一个180°激励,在从所述第一个180°激励起的第二时段过后实施第二个180°激励,以及在从所述第二个180°激励起的第二时段过后实施第三个180°激励;回波读出设备,用于在从所述第一个180°激励到所述第二个180°激励的时段期间在一个相位轴方向上施加一个零的时间积分值的梯度磁场,以读出一个第一自旋回波信号,以及在从所述第二个180°激励到所述第三个180°激励的时段期间在该相位轴方向上施加一个零的时间积分值的梯度磁场,以读出一个第二自旋回波信号;相位误差计算设备,用于根据所述第一和第二自旋回波信号确定在从所述90°激励到所述第一个180°激励的时段期间在该相位轴方向上自旋的相位误差;控制设备,用于试探在所述梯度磁场的极性颠倒的情况下重复两次由所述自旋激励设备实施的自旋激励、由所述回波读出设备进行的回波读出,和由所述相位误差计算设备进行的相位误差计算;和平均计算设备,用于确定前者所确定的相位误差与后者所确定的相位误差的平均值。
(7)按照用于解决前述问题的另一个方面,本发明是一种记录媒体,特征在于以计算机可读的方式记录一个程序,用于使计算机执行下述功能:对物体自旋实施一个90°激励;在从所述90°激励起的第一时段过后实施第一个180°激励;在从所述第一个180°激励起的第二时段过后实施第二个180°激励;在从所述第二个180°激励起的第二时段过后实施第三个180°激励;在从所述第一个180°激励到所述第二个180°激励的时段期间在一个相位轴方向上施加一个零的时间积分值的梯度磁场,以读出一个第一自旋回波信号;在从所述第二个180°激励到所述第三个180°激励的时段期间在该相位轴方向上施加一个零的时间积分值的梯度磁场,以读出一个第二自旋回波信号;根据所述第一和第二自旋回波信号确定在从所述90°激励到所述第一个180°激励的时段期间在该相位轴方向上自旋的相位误差;试探在所述梯度磁场的极性颠倒的情况下重复两次所述各步骤;和确定前者所确定的相位误差与后者所确定的相位误差的平均值。
(8)按照用于解决前述问题的另一个方面,本发明是一种磁共振成像装置,用于用一个静态磁场、梯度磁场和一个高频磁场获得来自一个物体的磁共振信号,并根据所获得的所述磁共振信号构造一个图像,特征在于包含:自旋激励设备,用于对物体自旋实施一个90°激励,在从所述90°激励起的第一时段过后实施第一个180°激励,在从所述第一个180°激励起的第二时段过后实施第二个180°激励,以及在从所述第二个180°激励起的第二时段过后实施第三个180°激励;回波读出设备,用于在从所述第一个180°激励到所述第二个180°激励的时段期间在一个相位轴方向上施加一个零的时间积分值的梯度磁场,以读出一个第一自旋回波信号,以及在从所述第二个180°激励到所述第三个180°激励的时段期间在该相位轴方向上施加一个零的时间积分值的梯度磁场,以读出一个第二自旋回波信号;相位误差计算设备,用于根据所述第一和第二自旋回波信号确定在从所述90°激励到所述第一个180°激励的时段期间在该相位轴方向上自旋的相位误差;控制设备,用于试探在所述梯度磁场的极性颠倒的情况下重复两次由所述自旋激励设备实施的自旋激励、由所述回波读出设备进行的回波读出,和由所述相位误差计算设备进行的相位误差计算;和平均计算设备,用于确定前者所确定的相位误差与后者所确定的相位误差的平均值。
按照如就(5)-(8)所述的各方面的本发明,由于在90°与180°激励之间的时段期间在相位轴方向上没有施加梯度磁场,并且在各180°激励之间的时段期间在相位轴方向上施加各自具有零的时间积分值的梯度磁场,以读出各自的自旋回波信号,所以能获得两个自旋回波信号,这两个信号受到在90°与180°激励之间的时段期间在相位轴方向上的残磁的互相相反的影响。因此,通过计算就能从这些自旋回波信号确定自旋的一个相位误差。
通过试探在所述梯度磁场的极性颠倒的情况下重复两次这样的相位误差测量,然后求所得到的测量值的平均值,就能消除在各180°激励之间的时段期间施加的梯度磁场对所测量值的影响,获得高精确度的测量值。
在如就(5)-(8)所述的各方面的本发明中,通过在每个180°激励之前和之后在频率轴方向上施加分相器梯度磁场,就能消除由FID信号产生的干扰,并且也能防止受激回波的生成,以提高测量的精确性。
在这种情况下,通过在两次试探中施加极性相反的分相器梯度磁场,就能消除作为分相器施加的梯度磁场对测量值的影响,由此便于进行更高精确度的测量。
(9)按照用于解决前述问题的另一个方面,本发明是一种相位误差校正方法,特征在于包含:对物体自旋实施一个90°激励;在从所述90°激励起的第一时段过后实施第一个180°激励;在从所述第一个180°激励起的第二时段过后实施第二个180°激励;在从所述第二个180°激励起的第二时段过后实施第三个180°激励;在从所述第一个180°激励到所述第二个180°激励的段期间在一个相位轴方向上施加一个零的时间积分值的梯度磁场,以读出一个第一自旋回波信号;在从所述第二个180°激励到所述第三个180°激励的时段期间在该相位轴方向上施加一个零的时间积分值的梯度磁场,以读出一个第二自旋回波信号;根据所述第一和第二自旋回波信号确定在从所述90°激励到所述第一个180°激励的时段期间在该相位轴方向上自旋的相位误差;和在从90°激励到第一个180°激励的时段期间在相位轴方向上施加一个使所述相位误差为零的梯度磁场。
(10)按照用于解决前述问题的另一个方面,本发明是一种相位误差校正装置,特征在于包含:自旋激励设备,用于对物体自旋实施一个90°激励,在从所述90°激励起的第一时段过后实施第一个180°激励,在从所述第一个180°激励起的第二时段过后实施第二个180°激励,以及在从所述第二个180°激励起的第二时段过后实施第三个180°激励;回波读出设备,用于在从所述第一个180°激励到所述第二个180°激励的时段期间在一个相位轴方向上施加一个零的时间积分值的梯度磁场,以读出一个第一自旋回波信号,以及在从所述第二个180°激励到所述第三个180°激励的时段期间在该相位轴方向上施加一个零的时间积分值的梯度磁场,以读出一个第二自旋回波信号;相位误差计算设备,用于根据所述第一和第二自旋回波信号确定在从所述90°激励到所述第一个180°激励的时段期间在该相位轴方向上自旋的相位误差;和校正梯度磁场施加设备,用于在从90°激励到第一个180°激励的时段期间在相位轴方向上施加一个使所述相位误差为零的梯度磁场。
(11)按照用于解决前述问题的另一个方面,本发明是一种记录媒体,特征在于以计算机可读的方式记录一个程序,用于使计算机执行下述功能:对物体自旋实施一个90°激励;在从所述90°激励起的第一时段过后实施第一个180°激励;在从所述第一个180°激励起的第二时段过后实施第二个180°激励;在从所述第二个180°激励起的第二时段过后实施第三个180°激励;在从所述第一个180°激励到所述第二个180°激励的时段期间在一个相位轴方向上施加一个零的时间积分值的梯度磁场,以读出一个第一自旋回波信号;在从所述第二个180°激励到所述第三个180°激励的时段期间在该相位轴方向上施加一个零的时间积分值的梯度磁场,以读出一个第二自旋回波信号;根据所述第一和第二自旋回波信号确定在从所述90°激励到所述第一个180°激励的时段期间在该相位轴方向上自旋的相位误差;和在从一个90°激励到一个180°激励的时段期间施加一个使该相位轴方向上的所述相位误差为零的梯度磁场。
(12)按照用于解决前述问题的另一个方面,本发明是一种磁共振成像装置,用于用一个静态磁场、梯度磁场和一个高频磁场获得来自一个物体的磁共振信号,并根据所获得的所述磁共振信号构造一个图像,特征在于包含:自旋激励设备,用于对物体自旋实施一个90°激励,在从所述90°激励起的第一时段过后实施第一个180°激励,在从所述第一个180°激励起的第二时段过后实施第二个180°激励,以及在从所述第二个180°激励起的第二时段过后实施第三个180°激励;回波读出设备,用于在从所述第一个180°激励到所述第二个180°激励的时段期间在一个相位轴方向上施加一个零的时间积分值的梯度磁场,以读出一个第一自旋回波信号,以及在从所述第二个180°激励到所述第三个180°激励的时段期间在该相位轴方向上施加一个零的时间积分值的梯度磁场,以读出一个第二自旋回波信号;相位误差计算设备,用于根据所述第一和第二自旋回波信号确定在从所述90°激励到所述第一个180°激励的时段期间在该相位轴方向上自旋的相位误差;和校正梯度磁场施加设备,用于在从90°激励到第一个180°激励的时段期间在相位轴方向上施加一个使所述相位误差为零的梯度磁场。
按照如就(9)-(12)所述的各方面的发明,由于与(1)-(4)的各方面的发明类似地测量一个相位误差,并且在从90°激励到第一个180°激励的时段期间在相位轴方向上施加一个使该相位误差为零的梯度磁场,就能校正由于相位轴方向上残留梯度磁场而产生的相位误差。
(13)按照用于解决前述问题的另一个方面,本发明是一种相位误差校正方法,特征在于包含:对物体自旋实施一个90°激励;在从所述90°激励起的第一时段过后实施第一个180°激励;在从所述第一个180°激励起的第二时段过后实施第二个180°激励;在从所述第二个180°激励起的第二时段过后实施第三个180°激励;在从所述第一个180°激励到所述第二个180°激励的时段期间在一个相位轴方向上施加一个零的时间积分值的梯度磁场,以读出一个第一自旋回波信号;在从所述第二个180°激励到所述第三个180°激励的时段期间在该相位轴方向上施加一个零的时间积分值的梯度磁场,以读出一个第二自旋回波信号;根据所述第一和第二自旋回波信号确定在从所述90°激励到所述第一个180°激励的时段期间在该相位轴方向上自旋的相位误差;试探在所述梯度磁场的极性颠倒的情况下重复两次所述的步骤;确定前者所确定的相位误差与后者所确定的相位误差的平均值;和在从90°激励到第一个180°激励的时段期间在相位轴方向上施加一个使所述相位误差的平均为零的梯度磁场。
(14)按照用于解决前述问题的另一个方面,本发明是一种相位误差测量装置,特征在于包含:自旋激励设备,用于对物体自旋实施一个90°激励,在从所述90°激励起的第一时段过后实施第一个180°激励,在从所述第一个180°激励起的第二时段过后实施第二个180°激励,以及在从所述第二个180°激励起的第二时段过后实施第三个180°激励;回波读出设备,用于在从所述第一个180°激励到所述第二个180°激励的时段期间在一个相位轴方向上施加一个零的时间积分值的梯度磁场,以读出一个第一自旋回波信号,以及在从所述第二个180°激励到所述第三个180°激励的时段期间在该相位轴方向上施加一个零的时间积分值的梯度磁场,以读出一个第二自旋回波信号;相位误差计算设备,用于根据所述第一和第二自旋回波信号确定在从所述90°激励到所述第一个180°激励的时段期间在该相位轴方向上自旋的相位误差;控制设备,用于试探在所述梯度磁场的极性颠倒的情况下重复两次由所述自旋激励设备实施的自旋激励、由所述回波读出设备进行的回波读出,和由所述相位误差计算设备进行的相位误差计算;平均计算设备,用于确定前者所确定的相位误差与后者所确定的相位误差的平均值;和校正梯度磁场施加设备,用于在从90°激励到第一个180°激励的时段期间在相位轴方向上施加一个使所述相位误差的平均为零的梯度磁场。
(15)按照用于解决前述问题的另一个方面,本发明是一种记录媒体,特征在于以计算机可读的方式记录一个程序,用于使计算机执行下述功能:对物体自旋实施一个90°激励;在从所述90°激励起的第一时段过后实施第一个180°激励;在从所述第一个180°激励起的第二时段过后实施第二个180°激励;在从所述第二个180°激励起的第二时段过后实施第三个180°激励;在从所述第一个180°激励到所述第二个180°激励的时段期间在一个相位轴方向上施加一个零的时间积分值的梯度磁场,以读出一个第一自旋回波信号;在从所述第二个180°激励到所述第三个180°激励的时段期间在该相位轴方向上施加一个零的时间积分值的梯度磁场,以读出一个第二自旋回波信号;根据所述第一和第二自旋回波信号确定在从所述90°激励到所述第一个180°激励的时段期间在该相位轴方向上自旋的相位误差;试探在所述梯度磁场的极性颠倒的情况下重复两次所述各步骤;确定前者所确定的相位误差与后者所确定的相位误差的平均值;和在从90°激励到第一个180°激励的时段期间在相位轴方向上施加一个使所述相位误差的平均为零的梯度磁场。
(16)按照用于解决前述问题的另一个方面,本发明是一种磁共振成像装置,用于用一个静态磁场、梯度磁场和一个高频磁场获得来自一个物体的磁共振信号,并根据所获得的所述磁共振信号构造一个图像,特征在于包含:自旋激励设备,用于对物体自旋实施一个90°激励,在从所述90°激励起的第一时段过后实施第一个180°激励,在从所述第一个180°激励起的第二时段过后实施第二个180°激励,以及在从所述第二个180°激励起的第二时段过后实施第三个180°激励;回波读出设备,用于在从所述第一个180°激励到所述第二个180°激励的时段期间在一个相位轴方向上施加一个零的时间积分值的梯度磁场,以读出一个第一自旋回波信号,以及在从所述第二个180°激励到所述第三个180°激励的时段期间在该相位轴方向上施加一个零的时间积分值的梯度磁场,以读出一个第二自旋回波信号;相位误差计算设备,用于根据所述第一和第二自旋回波信号确定在从所述90°激励到所述第一个180°激励的时段期间在该相位轴方向上自旋的相位误差;控制设备,用于试探在所述梯度磁场的极性颠倒的情况下重复两次由所述自旋激励设备实施的自旋激励、由所述回波读出设备进行的回波读出,和由所述相位误差计算设备进行的相位误差计算;平均计算设备,用于确定前者所确定的相位误差与后者所确定的相位误差的平均值;和校正梯度磁场施加设备,用于在从90°激励到第一个180°激励的时段期间在相位轴方向上施加一个使所述相位误差的平均为零的梯度磁场。
按照如就(13)一(16)所述的各方面的发明,由于与(5)-(8)的各方面的发明类似地测量一个相位误差,并且在从90°激励到第一个180°激励的时段期间在相位轴方向上施加一个使该相位误差为零的梯度磁场,就能更精确地校正由于相位轴方向上残留梯度磁场而产生的相位误差。
因此,本发明能提供一种用于测量在相位轴方向上的自旋中的相位误差的方法和装置、一种用于校正这种误差的方法和装置、一种记录有用于操作计算机来执行一个相位误差测量功能的程序的记录媒体,一种记录有用于操作计算机来执行一个相位误差校正功能的程序的记录媒体,一个包含相位误差测量装置的磁共振成像装置,和一个包含相位误差校正装置的磁共振成像装置。
本发明的其它目的和优点将在以下对如附图所示的本发明最佳实施例的说明中显而易见。
附图说明
图1是表示一种FSE技术的一个脉冲序列的图;
图2是按照本发明的一个实施例的装置的框图;
图3是表示用于相位误差测量的示例性脉冲序列的图;
图4A-4C是相位误差计算的概念图;
图5是表示用于相位误差测量的示例性脉冲序列的图;
图6A-6C是相位误差计算的概念图;
图7是表示用于相位误差测量的示例性脉冲序列的图;
图8是表示用于相位误差测量的示例性脉冲序列的图;
图9是表示用于磁共振成像的示例性脉冲序列的图。
其中附图中的符号说明如下:
100磁系统
102主磁场磁部分
106梯度线圈部分
108RF线圈部分
300物体
500托板
130梯度驱动部分
140驱动部分
150数据收集部分
160控制部分
170数据处理部分
180显示部分
190操作部分
(图3、5、7、8、9和1)(右端)时间
优选实施例说明
现在将结合附图详细地说明本发明的几个实施例。图2表示一种磁共振成像装置的框图,该装置是本发明的一个实施例。该装置的配置代表按照本发明的装置的一个实施例。该装置的操作代表按照本发明的方法的一个实施例。
如图2中所示,本装置有磁系统100。磁系统100有一个主磁场磁体部分102、一个梯度线圈(gradient coil)部分106和一个RF线圈部分108。主磁场磁体部分102和线圈部分各包含一对跨一个空间互相面对的构件。这些部分有一个大致像盘一样的形状,并被安置得具有共同的中心轴。物体300被放在托板500上,由一个图中没有表示的载物架装置运进运出磁系统100的内部空间(腔室)。
主磁场磁体部分102在磁系统100的内部空间中生成一个静态磁场。该静态磁场的方向大致上与物体300的体轴的方向垂直。就是说,生成一个“纵向”磁场。主磁场磁体部分102例如是用永久磁体构造的。
梯度线圈部分106生成三个梯度磁场,用于把各自的梯度给予三个互相垂直的轴(即限幅轴、相位轴和频率轴)的方向上的静态磁场强度。
当把静态磁场空间中互相垂直的坐标轴表示为x、y和z时,其中任何一个轴都可以是限幅轴。在这种情况中,其余两个轴中的一个是相位轴,另一个是频率轴。此外,在保持限幅、相位和频率轴的相互垂直性的同时可以赋予它们相对于x、y和z轴的任意倾角。
在限幅轴方向上的梯度磁场被称作限幅梯度磁场。在相位轴方向上的梯度磁场被称作相位编码梯度磁场。在频率轴方向上的梯度磁场被称作读出梯度磁场。为了能生成这种梯度磁场,梯度线圈部分106有三个梯度线圈,图中没有表示出它们。自此以下,梯度磁场有时被简称为梯度(gradient)。
RF线圈部分108发送一个激励信号,用于激励静态磁场空间中物体300内的自旋。该RF线圈部分108也接收由被激励的自旋生成的磁共振信号。
RF线圈部分108有发送和接收线圈,图中未显示。发送和接收线圈可以用同一个线圈或者单独的专用线圈。
梯度线圈部分106与梯度驱动部分130连接。梯度驱动部分130向梯度线圈部分106提供驱动信号,以生成梯度磁场。梯度驱动部分130有三个驱动电路,图中没有显示它们,它们对应于梯度线圈部分106中的三个梯度线圈。
RF线圈部分108与RF驱动部分140连接。RF驱动部分140向RF线圈部分108提供驱动信号,以发送RF信号,由此激励物体300内的自旋。
RF线圈部分108连接到数据收集部分150。数据收集部分150收集由RF线圈部分108通过采样而接收的接收信号,并且把这些信号收集为数字数据。
梯度驱动部分130、RF驱动部分140和数据收集部分150连接到控制部分160。控制部分160控制梯度驱动部分130、RF驱动部分140和数据收集部分150完成成像。
控制部分160是用例如计算机构造的。控制部分160有一个存储器,它在图中没有表示。存储器存储用于该控制部分160的程序和几种数据。控制部分160的功能是通过计算机执行在存储器中存储的程序而实现的。
数据收集部分150的数据的输出,连接到数据处理部分170。数据收集部分150所收集的数据,被输入到数据处理部分170。数据处理部分170是用例如计算机构造的。数据处理部分170有一个存储器,它在图中没有表示。存储器存储用于该数据处理部分170的程序和几种数据。
数据处理部分170与控制部分160连接。数据处理部分在控制部分160之上,控制控制部分160。本装置的功能是通过数据处理部分170执行在存储器中存储的程序而实现的。
数据处理部分170把由数据收集部分所收集的数据存储到存储器中。在存储器中形成一个数据空间。该数据空间构成一个二维的傅立叶(Fourier)空间。自此以下有时把傅立叶空间称作k空间。数据处理部分170对k空间中的数据进行一个二维傅立叶逆变换,以重构物体300的图像。
数据处理部分170与显示部分180和操作部分190连接。显示部分180包含一个图形显示等等。操作部分190包含一个配置有指示设备的键盘等等。
显示部分180显示从数据处理部分170输出的被重构的图像和几种信息。操作部分190被用户操作,该部分190向数据处理部分170输入几个命令、信息等等。用户通过显示部分180和操作部分190互动地操作本装置。
采用本装置的这种配置,执行对在相位轴方向上由于残磁而产生的相位误差的测量和校正。一个由磁系统100、梯度驱动部分130、RF驱动部分140、数据收集部分150、控制部分160和数据处理部分170组成的部分,是本发明的相位误差测量装置的一个实施例。该部分也是本发明的相位误差校正装置的一个实施例。该装置的配置代表按照本发明的装置的一个实施例。该装置的操作代表按照本发明的方法的一个实施例。
图3中表示一个由本装置实施的一个用于相位误差测量的示例性脉冲序列。在图3中,(1)是一个RF激励序列;(2)是一个限幅梯度序列Gs;(3)是相位轴方向上的一个梯度序列Gp;(4)是频率轴方向上的一个梯度序列Gr;和(5)是一个自旋回波序列SE。这些序列沿一个时间轴t发生。
(1)和(2)处所示的用于选择性的激励的序列与图1中所示的相同。一个由磁系统100和在实施这些序列中涉及的梯度驱动部分130和RF驱动部分140组成的部分,是本发明的自旋激励装置的一个实施例。
如(3)处所示,梯度Gp11、Gp12和Gp13的极性变换两次,它们在第一个与第二个180°激励之间的时段期间被施加在相位轴方向上;梯度Gp21、Gp22和Gp23的极性变换两次,它们在第二个与第三个180°激励之间的时段期间被施加。在90°与第一个180°激励之间的时段期间没有梯度被施加。此外,在本例中,在频率轴方向上没有梯度被施加,如(4)所示的那样。
在第一个与第二个180°激励之间的时段期间,首先施加负极性的梯度Gp11;然后,极性转换,施加正极性的梯度Gp12;然后,极性再次转换,施加负极性的梯度Gp13。
同样,在第二个与第三个180°激励之间的时段期间,首先施加负极性的梯度Gp21;然后,极性转换,施加正极性的梯度Gp22;然后,极性再次转换,施加负极性的梯度Gp23。
通过这种梯度,在相位轴方向上交替地实现了自旋的移相和再定相。具体来说,由梯度Gp11实现移相,由梯度Gp12实现再定相,由梯度Gp13实现移相。类似地,由梯度Gp21、Gp22和Gp23分别实现移相、再定相和移相。
梯度Gp11、Gp12和Gp13各自被设计得在整个梯度上有一个零的时间积分值。具体来说,梯度Gp12的前一半的时间积分值等于梯度Gp11的时间积分值,后一半的时间积分值等于梯度Gp13的时间积分值。对于梯度Gp21、Gp22和Gp23来说也是一样。自此以下,时间积分值将被简称为积分值。
梯度Gp11、Gp12和Gp13生成一个自旋回波SE1,如(5)处所示。自旋回波SE1在由该梯度产生的再定相与前面的移相平衡(equilibrate)时的点聚焦,在这个时刻具有一个最大幅值即峰值。
尽管自旋回波SE1的峰值应当在梯度Gp12的积分值等于梯度Gp11的积分值时发生,即应当在从90°激励起的时段TE过后发生,当在相位轴方向上存在残磁Gp0时,由于在90°与180°激励之间由梯度Gp0产生的移相,在自旋回波SE1的峰值出现的时刻发生一个误差Δt1。生成涡流的时,经历类似的效应。下面将对残磁作出说明。
当由残磁Gp0产生的移相的方向与梯度Gp11产生的相同时,误差Δt1是正的,方向相反时,误差Δt1是负的。尽管以下说明是对正误差的情况作说明的,如果误差是负的,以下说明同样成立,只是要把正的和负的互换。
梯度Gp21、GP22和Gp23生成一个自旋回波SE2。自旋回波SE2在由梯度Gp22产生的再定相与前面的移相平衡时的点聚焦。就是说,它在自旋回波SE1的峰值生成之后,在梯度Gp21和Gp22的积分值等于梯度Gp12和Gp13的积分值时聚焦。所以,自旋回波SE2有峰值的时刻,不同于其正确的聚焦时刻,误差是-Δt2。
有关自旋回波SE1和SE2的数据,由数据收集部分150收集。一个由磁系统100和在自旋回波SE1和SE2的生成中涉及的梯度驱动部分130和数据收集部分150组成的部分,是本发明的回波读出装置的一个实施例。
误差+Δt1与-Δt2之间的差,代表自旋的一个相位误差,这是由于在90°与180°激励之间的时段期间的残磁Gp0而产生的。该相位误差与该时段期间的残磁的积分值成正比。
由于在这个脉冲序列中从90°激励到180°激励时在相位轴方向上没有残磁Gp0以外的梯度,这个+Δt1与-Δt2之间的差精确地代表了由于残磁Gp0而产生的相位误差。
有关自旋回波SE1和SE2的数据,被输入到数据处理部分170。数据处理部分170用该输入数据进行计算,以确定该相位误差的一个测量值。
图4中表示了一个由数据处理部分170进行相位误差计算的概念图。自旋回波SE1的数据有一个与时间轴t上的原点偏离+Δt1的峰值位置,如图4中示意图(1)的左示意图所示的那样。
当把这种回波数据进行傅立叶变换时,就得到一个复数数据序列。该复数数据序列有一个对称于原点的线性相位特征曲线(profile),如右示意图所示的那样。应当注意的是,假设在Δt=0处的相位特征曲线中,所有相位是零。在以下说明作同样的假设。
相位特征曲线的斜率代表自旋回波SE1的峰值位置偏移量+Δt1。斜率的大小与峰值位置的偏移量+Δt1成正比。斜率的符号是正的,与峰值位置偏移量的符号一致。
自旋回波SE2的数据有一个与时间轴t上的原点偏离-Δt2的峰值位置,如图4中示意图(2)的左示意图所示的那样。当把这种回波数据进行傅立叶变换时,就得到一个复数数据序列。相位特征曲线变成一个对称于原点的线性特征曲线,如右示意图所示的那样。
该相位特征曲线的斜率与峰值位置的偏移量Δt2成正比。斜率的符号是负的,与峰值位置偏移量的符号一致。
通过确定这些相位特征曲线之间的差,就获得一个如示意图(3)所示的相位差特征曲线。该相位差特征曲线代表这两个自旋回波的峰值位置偏移量的差,就是说,代表由于在90°与180°激励之间的时段期间的残磁Gp0而产生的自旋的一个相位误差。该相位差特征曲线由其斜率的大小来代表相位误差的大小,由其斜率的符号来代表相位误差的符号。
如上所述的相位差特征曲线的计算是在数据处理部分170执行的。数据处理部分170是本发明的相位误差计算装置的一个实施例。
如果通过施加梯度Gp11、Gp12、Gp13、Gp21、Gp22和Gp23而生成伴生残磁,则相位误差的测量值含有一个由于该伴生残磁的效应而产生的误差。
作为针对这种情况的措施,采用如图5中示意性表示的序列来收集自旋回波SE1’和SE2’,其中,使梯度Gp11、Gp12、Gp13、Gp21、GP22和Gp23的极性全部翻转。
图5中的梯度Gp11’、Gp12’、Gp13’、Gp21’、Gp22’和Gp23’的极性与图3中的梯度Gp11、Gp12、Gp13、Gp21、Gp22和Gp23的极性相反。除极性外,图3和图5的梯度之间的条件都是相同的。
通过根据这种自旋回波SE1’和SE2’的聚焦时间的误差+Δt1’与-Δt2’之间的差确定一个相位误差,就能获得一个测量值,其中,由梯度Gp11、Gp12、Gp13、Gp21、Gp22和Gp23产生的伴生残磁的效应是相反的。
图6表示一个根据自旋回波SE1’和SE2’的相位误差计算的概念图。自旋回波SE1’的数据有一个与时间轴t上的原点偏离-Δt1’的峰值位置,如图6中(1)的左示意图所示的那样。
通过对该回波数据进行傅立叶变换而得到的一个复数数据序列的相位特征曲线,如右示意图所示的那样。该相位特征曲线的斜率与峰值位置的偏移量Δt1’成正比。斜率的符号是负的,对应于峰值位置偏移的方向。
自旋回波SE2’的数据有一个与时间轴t上的原点偏离+Δt2’的峰值位置,如图6中示意图(2)的左示意图所示的那样。通过对该回波数据进行傅立叶变换而得到的相位特征曲线,如右示意图所示的那样。该相位特征曲线的斜率与峰值位置的偏移量Δt2’成正比。斜率的符号是正的,对应于峰值位置偏移的方向。
通过确定这些相位特征曲线之间的差,就获得一个如示意图(3)所示的相位差特征曲线。应当注意的是,该相位差特征曲线是通过(2)-(1)的减操作获得的,因为梯度Gp11’、Gp12’、Gp13’、Gp21’、GP22’和Gp23’的极性已经被翻转(inverted)。
该相位差特征曲线也代表由于在90°与180°激励之间的时段期间的残磁Gp0而产生的自旋的一个相位误差。然而,与梯度Gp11、Gp12、Gp13、Gp21、GP22和Gp23伴生的残磁的效应是相反的。
图3中所示的脉冲序列的实施和根据两个自旋回波SE1和SE2的相位误差计算,以及图5中所示的脉冲序列的实施和根据两个自旋回波SE1’和SE2’的相位误差计算,是在控制部分160的控制下顺序地进行的。控制部分160是本发明的控制装置的一个实施例。
通过计算图4中(3)的相位差特征曲线和图6中(3)的相位差特征曲线的平均,就能消除伴生残磁引起的误差,由此获得一个具有良好的精确性的相位差特征曲线或相位误差测量值。
这两个相位差特征曲线的平均,是由数据处理部分170计算的。数据处理部分170是本发明的求平均值装置的一个实施例。
伴随第一个180°激励而生成的一个FID(自由感应衰减)信号,作为噪声被混杂在自旋回波SE1中,由此不利地影响相位误差测量的精确性。
相应地,为了去除基于FID的噪声,在第一个180°激励之前和之后在频率方向上分别施加梯度Gr11和Gr12,在第二个180°激励之前和之后分别施加梯度Gr21和Gr22,如图7中的(4)示例性地表示的那样。
在第一个180°激励之后施加的梯度Gr12是一个所谓的分相器,它分散被激励自旋的相位,使FID信号的信号强度为零,避免噪声影响到自旋回波SE1。
在第一个180°激励之前施加的梯度Gr11,是一个用于预先分散自旋的相位的梯度,以便防止自旋回波SE1受该分相器的影响。被梯度Gr11分散的相位被在该180°激励之后施加的梯度Grl2再定相(rephased),返回到合适的状态。
在第二个180°激励之后施加的梯度Gr22也是一个分相器,它分散被激励自旋的相位,使FID信号的信号强度为零,避免噪声影响到自旋回波SE2。
在第二个180°激励之前施加的梯度Gr21,是一个用于预先分散自旋的相位的梯度,以便防止自旋回波SE2受该分相器的影响。被梯度Gr21分散的相位被在该180°激励之后施加的梯度Gr22进行再定相,返回到合适的状态。
在第二个180°激励之后施加的梯度Gr22也是受激回波(stimulated echo)的一个分相器。生成该受激回波的自旋,从第一个180°激励到第二个180°激励以纵向磁化的形式存在,它们在这个时段期间不受梯度磁场的影响。
该纵向磁化被第二个180°激励转变成横向磁化,变成受激回波生成的源。通过由分相器Gr22移相,该回波的聚焦被禁止,由此避免受激回波生成。这样,就能防止受激回波混入自旋回波SE2中,干扰相位误差测量。
一个由磁系统和在梯度Gr11、Gr12、Gr21、Gr22的施加中涉及的梯度驱动部分130组成的部分,是本发明的分相器梯度磁场施加装置的一个实施例。
如果通过施加梯度Gr11、Gr12、Gr21和Gr22而生成伴生残磁,则相位误差的测量值含有一个由于该伴生残磁的效应而产生的误差。
作为针对这种情况的措施,采用如图8中示意性表示的序列来收集自旋回波SE1’和SE2’,其中,使梯度Gp11、Gp12、Gp13、Gp21、GP22和Gp23的极性以及梯度Gr11、Gr12、Gr21和Gr22的极性,全部是翻转的。
图8中的梯度Gp11’、Gp12’、Gp13’、Gp21’、Gp22’和Gp23’以及梯度Gr11’、Gr12’、Gr21’和Gr22’的极性,与图7中的梯度Gp11、Gp12、Gp13、Gp21、Gp22和Gp23以及梯度Gr11、Gr12、Gr21和Gr22的极性相反。除极性外,图7和图8的梯度之间的条件都是相同的。
通过根据这种自旋回波SE1’和SE2’的聚焦时间的误差+Δt1’与-Δt2’之间的差确定一个相位误差,就能获得一个测量值,其中,由梯度Gp11、Gp12、Gp13、Gp21、Gp22和Gp23以及梯度Gr11、Gr12、Gr21和Gr22产生的伴生残磁的效应是相反的。
这样,通过对从图7中的SE1和SE2获得的相位误差或相位差特征曲线和从图8中的SE1’和SE2’获得的相位误差或相位差特征曲线求平均值,就能获得不受由梯度Gp11、Gp12、Gp13、Gp21、Gp22和Gp23以及梯度Gr11、Gr12、Gr21和Gr22产生的伴生残磁的影响的相位误差测量值。
根据如此获得的相位误差,进行相位误差的校正。由于相位误差是因为在90°激励和180°激励之间的时段期间的残磁Gp0的移相而生成的,可以在该时段内反向地施加一个具有与残磁Gp0的积分值相等的积分值的梯度。
由于一个梯度的积分值与移相的量之间的关系理论上是确定的,通过计算,根据相位误差测量值就能获得用于相位误差校正的梯度的积分值。
在通过FSE技术进行磁共振成像时,由一个脉冲序列进行扫描,该脉冲序列在相位轴方向上附加一个具有这样一个积分值的梯度。图8中表示一例这样的脉冲序列。
如图9中(3)所示,在90°和180°激励之间的时段期间在相位轴方向上施加一个梯度Gp0’。梯度Gp0’的积分值等于残磁在该时段中的积分值,并且具有相反的极性。这使该时段中的移相为零,由此消除相位误差。一个由磁系统100和在梯度Gp0’的施加中涉及的梯度驱动部分130组成的部分,是本发明的校正梯度磁场施加装置的一个
实施例。
除了梯度Gp0’外,这个脉冲序列与图1中所示的普通FSE技术的脉冲序列是相同的。通过这个脉冲序列,收集到一个90°激励的例如16个视象的自旋回波数据。通过重复收集例如32次,获得512个视面的自旋回波数据。通过对自旋回波数据进行二维傅立叶变换,重构一个图像。
由于不存在相位误差,所有自旋回波的峰值都发生在适当的时刻。因此,就能获得一个高质量的重构图像,它不含有由于相位误差的噪声因素。
相位误差的测量和校正,最好每次在用FSE技术进行成像之前进行。这使扫描总是能适合最近的状态-即使残磁的状态被前面的扫描改变,由此便于总是能获得高质量的重构图像。
应当注意的是,以上所述的相位误差测量和校正,不但能与采用FSE技术的成像术结合,也能与由FSE技术衍生的采用脉冲序列的成像术结合,这类衍生的技术诸如是IR-FSE(反相恢复FSE)和快速FLAIR(快速流体细化的反相恢复)技术。
将一个用于操作计算机执行如上所述的相位误差测量功能的程序,以计算机可读的方式记录在一个记录媒体上。该记录媒体可采用例如磁记录媒体、光记录媒体和磁-光记录媒体中的任何一种,或者任何其它合适的记录媒体。该记录媒体可以是半导体存储媒体。在本说明书中存储媒体与记录媒体是同义的。
此外,将一个用于操作计算机执行如上所述的相位误差校正功能的程序,以计算机可读的方式记录在一个记录媒体上。该记录媒体可采用例如磁记录媒体、光记录媒体和磁-光记录媒体中的任何一种,或者任何其它合适的记录媒体。该记录媒体可以是半导体存储媒体。在本说明书中存储媒体与记录媒体是同义的。
在不偏离本发明的精神和范围的情况下,可以设置本发明的许多不同的实施例。应当明白的是,本发明除了由后附的权利要求书界定外,并不仅限于本说明书中所描述的特定实施例。

Claims (30)

1.一种相位误差测量方法,包含下列步骤:
对物体自旋实施一个90°激励;
在从所述90°激励起的第一时段过后实施第一个180°激励;
在从所述第一个180°激励起的第二时段过后实施第二个180°激励;
在从所述第二个180°激励起的第二时段过后实施第三个180°激励;
在从所述第一个180°激励到所述第二个180°激励的时段期间在一个相位轴方向上施加一个具有零的时间积分值的第一梯度磁场,该第一梯度磁场包括第一极性的第一脉冲、紧随其后的与所述第一极性相对的第二极性的第二脉冲、紧随其后的所述第一极性的第三脉冲,以产生一个第一自旋回波信号,具有由从所述第一自旋回波的波峰到所述第一脉冲和所述第二脉冲的时间积分相等时的时间点之间的距离决定的第一相位误差;
在从所述第二个180°激励到所述第三个180°激励的时段期间在该相位轴方向上施加一个具有零的时间积分值的第二梯度磁场,该第二梯度磁场包括第一极性的第一脉冲、紧随其后的与所述极性相对的第二极性的第二脉冲、紧随其后的所述极性的第三脉冲,以产生一个第二自旋回波信号,具有由从所述第二自旋回波的波峰到所述第一脉冲和所述第二脉冲的时间积分相等时的时间点之间的距离决定的第二相位误差;和
根据利用所述第一和第二自旋回波信号的所述第一相位误差和所述第二相位误差的计算,确定在从所述90°激励到所述第一个180°激励的时段期间在该相位轴方向上由于残磁或涡流或其结合造成的自旋的相位误差。
2.权利要求1的相位误差测量方法,特征在于进一步包含下列步骤:
在每个所述180°激励之前和之后在一个频率轴方向上施加分相器梯度磁场。
3.权利要求1的相位误差测量方法,特征在于进一步包含下列步骤:
试探在所述梯度磁场的极性颠倒的情况下重复两次所述的步骤;和
确定前者所确定的相位误差与后者所确定的相位误差的平均值。
4.权利要求3的相位误差测量方法,特征在于进一步包含下列步骤:
在每个所述180°激励之前和之后在一个频率轴方向上施加分相器梯度磁场。
5.权利要求4的相位误差测量方法,特征在于进一步包含下列步骤:
在所述的其极性颠倒的两次试探中施加所述分相器梯度磁场。
6.权利要求3的相位误差测量方法,特征在于进一步包含下列步骤:
在从一个90°激励到一个180°激励的时段期间在相位轴方向上施加一个使所述相位误差的平均为零的梯度磁场。
7.权利要求6的相位误差测量方法,特征在于进一步包含下列步骤:
在每个所述180°激励之前和之后在一个频率轴方向上施加分相器梯度磁场。
8.权利要求7的相位误差测量方法,特征在于进一步包含下列步骤:
在所述的其极性颠倒的两次试探中施加所述分相器梯度磁场。
9.权利要求1的相位误差测量方法,特征在于进一步包含下列步骤:
在从90°激励到第一个180°激励的时段期间在相位轴方向上施加一个使所述相位误差为零的梯度磁场。
10.权利要求9的相位误差测量方法,特征在于进一步包含下列步骤:
在每个所述180°激励之前和之后在一个频率轴方向上施加分相器梯度磁场。
11.一种相位误差测量装置,包含:
第一装置,用于对物体自旋实施一个90°激励,在从所述90°激励起的第一时段过后实施第一个180°激励,在从所述第一个180°激励起的第二时段过后实施第二个180°激励,以及在从所述第二个180°激励起的第二时段过后实施第三个180°激励;
第二装置,用于在从所述第一个180°激励到所述第二个180°激励的时段期间在一个相位轴方向上施加一个具有零的时间积分值的第一梯度磁场,该第一梯度磁场包括第一极性的第一脉冲、紧随其后的与所述第一极性相对的第二极性的第二脉冲、紧随其后的所述第一极性的第三脉冲,以产生一个第一自旋回波信号,具有由从所述第一自旋回波信号的波峰到所述第一脉冲和所述第二脉冲的时间积分相等时的时间点之间的距离决定的第一相位误差,以及在从所述第二个180°激励到所述第三个180°激励的时段期间在该相位轴方向上施加一个具有零的时间积分值的第二梯度磁场,该第二梯度磁场包括第一极性的第一脉冲、紧随其后的与所述第一极性相对的第二极性的第二脉冲、紧随其后的所述第一极性的第三脉冲,以产生一个第二自旋回波信号,具有由从所述第二自旋回波的波峰到所述第一脉冲和所述第二脉冲的时间积分相等时的时间点之间的距离决定的第二相位误差;和
第三装置,用于根据利用所述第一和第二自旋回波信号的所述第一相位误差和所述第二相位误差的计算,确定在从所述90°激励到所述第一个180°激励的时段期间在该相位轴方向上由于残磁或涡流或其结合造成的自旋的相位误差。
12.权利要求11的相位误差测量装置,特征在于进一步包含:
第四装置,用于在每个所述180°激励之前和之后在一个频率轴方向上施加分相器梯度磁场。
13.权利要求11的相位误差测量装置,特征在于进一步包含:
第五装置,用于试探在所述梯度磁场的极性颠倒的情况下重复两次由所述第一装置实施的自旋激励、由所述第二装置进行的自旋回波读出、和由所述第三装置进行的相位误差计算;和
第六装置,用于确定前者所确定的相位误差与后者所确定的相位误差的平均值。
14.权利要求13的相位误差测量装置,特征在于进一步包含:
第四装置,用于在每个所述180°激励之前和之后在一个频率轴方向上施加分相器梯度磁场。
15.权利要求14的相位误差测量装置,其中,所述第五装置使所述第四装置在所述的其极性颠倒的两次试探中施加所述分相器梯度磁场。
16.权利要求13的相位误差测量装置,特征在于进一步包含:
第七装置,用于在从90°激励到180°激励的时段期间在相位轴方向上施加一个使所述相位误差的平均为零的梯度磁场。
17.权利要求16的相位误差测量装置,特征在于进一步包含:
第四装置,用于在每个所述180°激励之前和之后在一个频率轴方向上施加分相器梯度磁场。
18.权利要求17的相位误差测量装置,其中,所述第四装置在所述的其极性颠倒的两次试探中施加所述分相器梯度磁场。
19.权利要求11的相位误差测量装置,特征在于进一步包含:
第六装置,用于在从90°激励到180°激励的时段期间在相位轴方向上施加一个使所述相位误差为零的梯度磁场。
20.权利要求19的相位误差测量装置,特征在于进一步包含:第四装置,用于在每个所述180°激励之前和之后在一个频率轴方向上施加分相器梯度磁场。
21.一种磁共振成像装置,用于用一个静态磁场、梯度磁场和一个高频磁场获得来自一个物体的磁共振信号,并根据所获得的所述磁共振信号构造一个图像,该装置包含:
第一装置,用于对物体自旋实施一个90°激励,在从所述90°激励起的第一时段过后实施第一个180°激励,在从所述第一个180°激励起的第二时段过后实施第二个180°激励,以及在从所述第二个180°激励起的第二时段过后实施第三个180°激励;
第二装置,用于在从所述第一个180°激励到所述第二个180°激励的时段期间在一个相位轴方向上施加一个具有零的时间积分值的第一梯度磁场,该第一梯度磁场包括第一极性的第一脉冲、紧随其后的与所述第一极性相对的第二极性的第二脉冲、紧随其后的所述第一极性的第三脉冲,以产生一个第一自旋回波信号,具有由从所述第一自旋回波信号的波峰到所述第一脉冲和所述第二脉冲的时间积分相等时的时间点之间的距离决定的第一相位误差,以及在从所述第二个180°激励到所述第三个180°激励的时段期间在该相位轴方向上施加一个具有零的时间积分值的第二梯度磁场,该第二梯度磁场包括第一极性的第一脉冲、紧随其后的与所述第一极性相对的第二极性的第二脉冲、紧随其后的所述第一极性的第三脉冲,以产生一个第二自旋回波信号,具有由从所述第二自旋回波的波峰到所述第一脉冲和所述第二脉冲的时间积分相等时的时间点之间的距离决定的第二相位误差;和
第三装置,用于根据利用所述第一和第二自旋回波信号的所述第一相位误差和所述第二相位误差的计算,确定在从所述90°激励到所述第一个180°激励的时段期间在该相位轴方向上由于残磁或涡流或其结合造成的自旋的相位误差。
22.权利要求21的磁共振成像装置,特征在于进一步包含:
第四装置,用于在每个所述180°激励之前和之后在一个频率轴方向上施加分相器梯度磁场。
23.权利要求21的磁共振成像装置,特征在于进一步包含:
第五装置,用于试探在所述梯度磁场的极性颠倒的情况下重复两次由所述第一装置实施的自旋激励、由所述第二装置进行的自旋回波读出、和由所述第三进行的相位误差计算;和
第六装置,用于确定前者所确定的相位误差与后者所确定的相位误差的平均值。
24.权利要求23的磁共振成像装置,特征在于进一步包含:
第四装置,用于在每个所述180°激励之前和之后在一个频率轴方向上施加分相器梯度磁场。
25.权利要求24的磁共振成像装置,特征在于,所述第五装置使所述第四装置在所述的其极性颠倒的两次试探中施加所述分相器梯度磁场。
26.权利要求23的磁共振成像装置,特征在于进一步包含:
第六装置,用于在从90°激励到180°激励的时段期间在相位轴方向上施加一个使所述相位误差的平均为零的梯度磁场。
27.权利要求26的磁共振成像装置,特征在于进一步包含:
第四装置,用于在每个所述180°激励之前和之后在一个频率轴方向上施加分相器梯度磁场。
28.权利要求27的磁共振成像装置,特征在于,所述第五装置使所述第四装置在所述的其极性颠倒的两次试探中施加所述分相器梯度磁场。
29.权利要求21的磁共振成像装置,特征在于进一步包含:
第六装置,用于在从90°激励到180°激励的时段期间在相位轴方向上施加一个使所述相位误差为零的梯度磁场
30.权利要求29的磁共振成像装置,特征在于进一步包含:
第四装置,用于在每个所述180°激励之前和之后在一个频率轴方向上施加分相器梯度磁场。
CNB011357940A 2000-10-20 2001-10-19 相位误差测量方法、校正方法和装置、磁共振成像装置 Expired - Fee Related CN1196930C (zh)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
JP2000320827A JP3875479B2 (ja) 2000-10-20 2000-10-20 磁気共鳴撮影装置
JP320827/2000 2000-10-20
JP320827/00 2000-10-20

Publications (2)

Publication Number Publication Date
CN1350176A CN1350176A (zh) 2002-05-22
CN1196930C true CN1196930C (zh) 2005-04-13

Family

ID=18799048

Family Applications (1)

Application Number Title Priority Date Filing Date
CNB011357940A Expired - Fee Related CN1196930C (zh) 2000-10-20 2001-10-19 相位误差测量方法、校正方法和装置、磁共振成像装置

Country Status (5)

Country Link
US (1) US6664787B2 (zh)
EP (1) EP1199578A3 (zh)
JP (1) JP3875479B2 (zh)
KR (1) KR100458774B1 (zh)
CN (1) CN1196930C (zh)

Families Citing this family (24)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP3453089B2 (ja) * 1999-08-20 2003-10-06 ジーイー横河メディカルシステム株式会社 Mri装置
GB2374672A (en) * 2001-04-20 2002-10-23 Marconi Medical Systems Uk Ltd Magnetic Resonance Imaging
JP2005512648A (ja) * 2001-12-19 2005-05-12 コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ 周期的に振幅が変調されたクラッシャーグラジエントを有する拡散重み付けされた複数スピンエコー(rare)シーケンス
CN100375904C (zh) * 2002-12-24 2008-03-19 深圳安科高技术股份有限公司 一种射频脉冲直角相位设置方法及电路
JP2005152114A (ja) * 2003-11-21 2005-06-16 Ge Medical Systems Global Technology Co Llc Mri方法およびmri装置
JP3968352B2 (ja) 2004-02-03 2007-08-29 ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー Mri装置
JP3968353B2 (ja) * 2004-02-18 2007-08-29 ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー Mri装置
JP4619674B2 (ja) * 2004-03-24 2011-01-26 株式会社東芝 磁気共鳴イメージング装置
JP4347788B2 (ja) * 2004-12-01 2009-10-21 ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー Mri装置
US7187171B2 (en) * 2005-02-10 2007-03-06 Kaburhiki Kaisha Toshiba Magnetic resonance imaging apparatus and its control method
DE102005015069B4 (de) * 2005-04-01 2008-03-27 Siemens Ag Verfahren zur Vermeidung linearer Phasenfehler in der Magnetresonanz-Spektroskopie
WO2008005409A2 (en) * 2006-06-29 2008-01-10 Vassol Inc. Automatic segmentation of stationary tissue in pcmr imaging
US8085041B2 (en) * 2008-04-10 2011-12-27 General Electric Company Three-point method and system for fast and robust field mapping for EPI geometric distortion correction
DE102012205587B4 (de) 2012-04-04 2013-12-24 Siemens Aktiengesellschaft Schichtspezifische Phasenkorrektur bei Schicht-Multiplexing
KR101474757B1 (ko) 2013-07-08 2014-12-19 삼성전자주식회사 자장 측정 방법 및 장치
EP3060116B1 (en) * 2013-10-21 2022-06-15 Koninklijke Philips N.V. Mr imaging with temperature mapping
CN104458785B (zh) * 2014-12-12 2016-09-07 中国科学院武汉物理与数学研究所 一种核磁共振波谱谱峰对齐及谱峰提取方法
US10126402B2 (en) * 2015-06-02 2018-11-13 Case Western Reserve University Magnetic resonance imaging (MRI) with artifact-free T2 mapping
CN106814339B (zh) * 2015-11-27 2019-11-12 西门子(深圳)磁共振有限公司 相位信息获取方法、系统及磁共振成像系统
CN107024670B (zh) * 2016-01-29 2020-08-04 上海联影医疗科技有限公司 磁共振系统的校正方法及装置
CN106353706B (zh) * 2016-10-08 2019-07-23 上海东软医疗科技有限公司 磁共振成像序列相位校正方法及磁共振设备
CN110988765B (zh) * 2019-12-25 2022-03-04 东软医疗系统股份有限公司 一种磁共振相位校正方法和装置
DE102020212281B4 (de) * 2020-09-29 2022-05-12 Siemens Healthcare Gmbh Verfahren zur zeitsparenden Erzeugung einer B0-Karte basierend auf einer Doppelecho-Sequenz mit stimulierten Echos und Magnetresonanzvorrichtung
CN115685032B (zh) * 2022-11-02 2023-06-23 佛山瑞加图医疗科技有限公司 一种幅条k空间下的快速自旋回波的校正方法及系统

Family Cites Families (26)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4714884A (en) * 1986-06-13 1987-12-22 General Electric Company Method of eliminating effects of spurious NMR signals caused by imperfect 180 degree RF pulses
US4978918A (en) * 1988-04-24 1990-12-18 Mitsubishi Denki Kabushiki Kaisha Magnetic resonance imaging method
DE4027252A1 (de) * 1990-08-29 1992-03-05 Philips Patentverwaltung Verfahren fuer die zweidimensionale kernresonanzspektroskopie
US5270654A (en) * 1991-07-05 1993-12-14 Feinberg David A Ultra-fast multi-section MRI using gradient and spin echo (grase) imaging
US5493224A (en) * 1992-03-03 1996-02-20 Hitachi, Ltd. Ultra high-speed magnetic resonance imaging method and apparatus
EP0577188B1 (en) * 1992-06-29 1999-09-01 Koninklijke Philips Electronics N.V. Method and apparatus for magnetic resonance imaging
US5281916A (en) * 1992-07-29 1994-01-25 General Electric Company NMR angiography using fast spin echo pulse sequences
JP3246020B2 (ja) * 1992-12-31 2002-01-15 株式会社島津製作所 Mrイメージング装置
JPH06217959A (ja) * 1993-01-26 1994-08-09 Toshiba Corp 磁気共鳴イメージングにおける位相エンコード量調整法
EP0633480A1 (en) * 1993-07-09 1995-01-11 Koninklijke Philips Electronics N.V. Artefact suppression in GRASE MR imaging
US5378985A (en) * 1993-07-15 1995-01-03 General Electric Company Fast spin echo prescan for MRI system
WO1995005610A1 (en) * 1993-08-13 1995-02-23 The Government Of The United States Of America As Represented By The Secretary, Department Of Health And Human Services Method for magnetic resonance spectroscopic imaging with multiple spin-echoes
JP3688773B2 (ja) * 1995-10-31 2005-08-31 株式会社東芝 Mri装置
JP3573570B2 (ja) * 1996-04-26 2004-10-06 株式会社日立メディコ 磁気共鳴イメージング装置
JP3384944B2 (ja) * 1996-07-11 2003-03-10 ジーイー横河メディカルシステム株式会社 Mri装置
US6078176A (en) * 1996-11-08 2000-06-20 General Electric Company Fast spin echo pulse sequence for diffusion weighted imaging
US6008647A (en) * 1997-02-11 1999-12-28 General Electric Company Method for reducing Maxwell term artifacts in fast spin echo MR images
US6011392A (en) * 1997-04-10 2000-01-04 General Electric Company Method for reducing Maxwell term artifacts in fast spin echo MR images
US5899858A (en) * 1997-04-10 1999-05-04 Mayo Foundation For Medical Education And Research MR imaging with enhanced sensitivity of specific spin motion
DE19826864A1 (de) * 1998-06-17 1999-12-23 Philips Patentverwaltung MR-Verfahren
US6369568B1 (en) * 1999-06-03 2002-04-09 Ge Medical Systems Global Technology Company, Llc Fast spin echo phase correction for MRI system
JP3365983B2 (ja) * 1999-09-28 2003-01-14 ジーイー横河メディカルシステム株式会社 Mri装置
KR100331449B1 (ko) * 1999-11-12 2002-04-09 윤종용 패스트 스핀 에코 영상법에서 위상 부호화 경사 자계펄스의 위상 오류 보정방법
US6392411B1 (en) * 2000-01-03 2002-05-21 Ge Yokogawa Medical Systems, Limited MR imaging method, phase shift measuring method and MR imaging system
JP3440049B2 (ja) * 2000-02-25 2003-08-25 ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー Mri装置
US6541971B1 (en) * 2001-06-28 2003-04-01 Koninklijke Philips Electronics, N.V. Multi-dimensional spatial NMR excitation

Also Published As

Publication number Publication date
US6664787B2 (en) 2003-12-16
KR20020031061A (ko) 2002-04-26
JP3875479B2 (ja) 2007-01-31
EP1199578A3 (en) 2004-06-02
CN1350176A (zh) 2002-05-22
EP1199578A2 (en) 2002-04-24
KR100458774B1 (ko) 2004-12-03
JP2002136498A (ja) 2002-05-14
US20020047708A1 (en) 2002-04-25

Similar Documents

Publication Publication Date Title
CN1196930C (zh) 相位误差测量方法、校正方法和装置、磁共振成像装置
CN1299121C (zh) 磁共振成像设备
CN1231177C (zh) 磁共振成像方法及实施该方法的装置
CN100339047C (zh) 磁共振成像设备和磁共振成像方法
JP3348572B2 (ja) Mrイメージング装置
CN1265762C (zh) 线圈灵敏图生成方法和平行成像方法以及磁共振成像装置
CN1248648C (zh) 磁共振成像设备
CN1683939A (zh) 在磁共振断层造影中融合加权成像的运动校正的多点方法
JP6285669B2 (ja) 磁気共鳴イメージング装置
CN105473067A (zh) 用于降低mr成像中声音噪声水平的系统和方法
JP3369688B2 (ja) 核磁気共鳴を用いた検査装置
CN105676155A (zh) 磁共振成像系统的磁场不均匀性值获取方法和失真校正方法
CN1196940C (zh) 磁共振信号获取方法及装置、记录媒体和磁共振成像装置
CN1237487C (zh) 图像处理方法和装置、成像装置
CN1920592A (zh) 相位分布及相位校正方法和装置以及磁共振成像方法和装置
US6528996B1 (en) Diffusion-weighted imaging method and apparatus for fast pulse sequence with MPG pulses
CN1258001A (zh) Mri系统中剩余磁化强度的修正
CN1192743C (zh) Mr成像方法,相位误差测量方法和mri系统
JP5971683B2 (ja) 磁気共鳴イメージング装置
CN1244821C (zh) 多重回波成像方法
CN1262238C (zh) 磁共振信号采集的方法和装置及磁共振成像装置
JP2009160051A (ja) 磁気共鳴イメージング装置
CN1346622A (zh) 磁谐振信号获取方法和装置、记录媒体和磁谐振图象装置
CN113406546B (zh) 一种抗涡流效应的跨膜水交换磁共振成像序列设计方法
CN109752681A (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
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: 20050413

Termination date: 20141019

EXPY Termination of patent right or utility model