JP2003135418A - ０次位相検出方法およびｍｒｉ装置 - Google Patents
０次位相検出方法およびｍｒｉ装置Info
- Publication number
- JP2003135418A JP2003135418A JP2001328447A JP2001328447A JP2003135418A JP 2003135418 A JP2003135418 A JP 2003135418A JP 2001328447 A JP2001328447 A JP 2001328447A JP 2001328447 A JP2001328447 A JP 2001328447A JP 2003135418 A JP2003135418 A JP 2003135418A
- Authority
- JP
- Japan
- Prior art keywords
- phase detection
- zsum
- pulse sequence
- phase
- order phase
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
- G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
- G01R33/48—NMR imaging systems
- G01R33/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
- G01R33/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
-
- 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/565—Correction of image distortions, e.g. due to magnetic field inhomogeneities
- G01R33/56554—Correction of image distortions, e.g. due to magnetic field inhomogeneities caused by acquiring plural, differently encoded echo signals after one RF excitation, e.g. correction for readout gradients of alternating polarity in EPI
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
(57)【要約】
【課題】ＭＲ信号の０次位相を正しく表す０次位相を求
める。 【解決手段】ＭＲ信号をフーリエ変換して得た全サンプ
リング点の複素ベクトルの合成ベクトルの位相を０次位
相とする。 【効果】正しい０次位相および０次位相差が得られるた
め、イメージング用データの補正やそれ以外の種々の処
理に利用されたときに正しい結果が得られる。
める。 【解決手段】ＭＲ信号をフーリエ変換して得た全サンプ
リング点の複素ベクトルの合成ベクトルの位相を０次位
相とする。 【効果】正しい０次位相および０次位相差が得られるた
め、イメージング用データの補正やそれ以外の種々の処
理に利用されたときに正しい結果が得られる。
Description
【０００１】
【発明の属する技術分野】本発明は、０次位相検出方法
およびＭＲＩ（Magnetic Resonance Imaging）装置に関
し、さらに詳しくは、ＭＲ信号の０次位相を正しく検出
することが出来る０次位相検出方法およびＭＲＩ装置に
関する。
およびＭＲＩ（Magnetic Resonance Imaging）装置に関
し、さらに詳しくは、ＭＲ信号の０次位相を正しく検出
することが出来る０次位相検出方法およびＭＲＩ装置に
関する。
【０００２】
【従来の技術】図５に、マルチショット拡散強調ＥＰＩ
（Echo Plannar Imaging）法によるイメージング用パル
スシーケンスの基本例を示す。このイメージング用パル
スシーケンスでは、励起パルスRF90とスライス勾配SG90
とを印加する。次にＭＰＧ（Motion Probing Gradien
t）パルスMPGを印加する。次に反転ＲＦパルスRF180と
スライス勾配SG180を印加する。次にＭＰＧパルスMPGを
印加する。次に位相エンコード勾配pdnを印加する。次
に交互に正負に反転するデータ収集用リード勾配ｒ１，
…，ｒｍを連続的に印加し且つ反転時に位相エンコード
勾配ｐ２，…，ｐＭを印加し、第１エコーｅ１から第Ｍ
エコーｅＭが順に集束するのとタイミングを合せてサン
プリングし、各エコーｅ１，…，ｅＭに対応したイメー
ジング用データＦ(n，1)，…，Ｆ(n，M)をそれぞれ収集
する。これを位相エンコード勾配pdnの大きさを変えな
がらｎ＝１，…，Ｎについて繰り返して、ｋ空間を埋め
るイメージング用データＦ(1，1)〜Ｆ(N，M)を収集す
る。これをＮショット・Ｍエコーという。また、ショッ
トに対して実行時間順に付けた番号ｎをショット番号と
いう。また、あるショットのエコー列のエコーに対して
集束時間順に付けた番号をエコー番号という。
（Echo Plannar Imaging）法によるイメージング用パル
スシーケンスの基本例を示す。このイメージング用パル
スシーケンスでは、励起パルスRF90とスライス勾配SG90
とを印加する。次にＭＰＧ（Motion Probing Gradien
t）パルスMPGを印加する。次に反転ＲＦパルスRF180と
スライス勾配SG180を印加する。次にＭＰＧパルスMPGを
印加する。次に位相エンコード勾配pdnを印加する。次
に交互に正負に反転するデータ収集用リード勾配ｒ１，
…，ｒｍを連続的に印加し且つ反転時に位相エンコード
勾配ｐ２，…，ｐＭを印加し、第１エコーｅ１から第Ｍ
エコーｅＭが順に集束するのとタイミングを合せてサン
プリングし、各エコーｅ１，…，ｅＭに対応したイメー
ジング用データＦ(n，1)，…，Ｆ(n，M)をそれぞれ収集
する。これを位相エンコード勾配pdnの大きさを変えな
がらｎ＝１，…，Ｎについて繰り返して、ｋ空間を埋め
るイメージング用データＦ(1，1)〜Ｆ(N，M)を収集す
る。これをＮショット・Ｍエコーという。また、ショッ
トに対して実行時間順に付けた番号ｎをショット番号と
いう。また、あるショットのエコー列のエコーに対して
集束時間順に付けた番号をエコー番号という。
【０００３】図６は、ｋ空間ＫＳにおけるイメージング
用データＦ(1,1)〜Ｆ(N,M)の収集軌跡（trajectory）を
示す模式図である。但し、Ｎ＝４，Ｍ＝４としている。
ｋ空間ＫＳを位相エンコード軸方向に第１行から第Ｎ×
Ｍ行（図６では第１６行）まで分割したとき、第ｎショ
ットの第ｍエコーで第（ｎ＋（ｍ−１）Ｎ）行のイメー
ジング用データＦ(n,m)を収集するように、位相エンコ
ードpdn，ｐ2，…，ｐＭを印加している。
用データＦ(1,1)〜Ｆ(N,M)の収集軌跡（trajectory）を
示す模式図である。但し、Ｎ＝４，Ｍ＝４としている。
ｋ空間ＫＳを位相エンコード軸方向に第１行から第Ｎ×
Ｍ行（図６では第１６行）まで分割したとき、第ｎショ
ットの第ｍエコーで第（ｎ＋（ｍ−１）Ｎ）行のイメー
ジング用データＦ(n,m)を収集するように、位相エンコ
ードpdn，ｐ2，…，ｐＭを印加している。
【０００４】さて、ＥＰＩ法は位相エラーに敏感である
ため、ＭＲ信号の位相を検出し、位相補正する必要があ
る。図７は、位相検出用パルスシーケンスの一例であ
る。この位相検出用パルスシーケンスは、図５のイメー
ジング用パルスシーケンスから位相エンコード勾配を省
略したパルスシーケンスである。順に集束する第１位相
検出用エコーＥ１から第Ｍ位相検出用エコーＥＭを基に
して第１位相検出用データＤ_1から第Ｍ位相検出用デー
タＤ_Mをそれぞれ収集する。
ため、ＭＲ信号の位相を検出し、位相補正する必要があ
る。図７は、位相検出用パルスシーケンスの一例であ
る。この位相検出用パルスシーケンスは、図５のイメー
ジング用パルスシーケンスから位相エンコード勾配を省
略したパルスシーケンスである。順に集束する第１位相
検出用エコーＥ１から第Ｍ位相検出用エコーＥＭを基に
して第１位相検出用データＤ_1から第Ｍ位相検出用デー
タＤ_Mをそれぞれ収集する。
【０００５】次に、第ｍ位相検出用データＤ_mをフーリ
エ変換し、複素ベクトルＺ(n)を得る。但し、１≦ｍ≦
Ｍ、１≦ｎ≦Ｎである。次に、第ｍ位相検出用データＤ
_mの１次位相φ１_mを次式により算出する。 φ１_m＝arg｛n=1ΣN-1（Ｚ(n+1)／Ｚ(n)）｝ …(1) ここで、arg｛｝は、複素数の偏角を表す関数である。
そして、次式により１次位相の補正を行う。 Ｚcor1(n)＝Ｚ(n)・exp｛−ｉ・φ１_m・（ｎ−１）｝ …(2) ここで、exp｛｝は、指数関数である。
エ変換し、複素ベクトルＺ(n)を得る。但し、１≦ｍ≦
Ｍ、１≦ｎ≦Ｎである。次に、第ｍ位相検出用データＤ
_mの１次位相φ１_mを次式により算出する。 φ１_m＝arg｛n=1ΣN-1（Ｚ(n+1)／Ｚ(n)）｝ …(1) ここで、arg｛｝は、複素数の偏角を表す関数である。
そして、次式により１次位相の補正を行う。 Ｚcor1(n)＝Ｚ(n)・exp｛−ｉ・φ１_m・（ｎ−１）｝ …(2) ここで、exp｛｝は、指数関数である。
【０００６】次に、第ｍ位相検出用データＤ_mの０次位
相φ０_mを次式により算出する。 φ０_m＝arg｛n=1ΣNＺcor1(n)｝ …(3)
相φ０_mを次式により算出する。 φ０_m＝arg｛n=1ΣNＺcor1(n)｝ …(3)
【０００７】上記１次位相φ１_mおよび０次位相φ０_m
は、イメージング用データの補正やそれ以外の種々の処
理に利用される。
は、イメージング用データの補正やそれ以外の種々の処
理に利用される。
【０００８】また、第ｍ位相検出用データＤ_mの０次位
相φ０_mと第（ｍ＋１）位相検出用データＤ_m+1の０次
位相φ０_m+1の０次位相差Δφ０_mが、次式により求め
られる。 Δφ０_m＝φ０_m+1−φ０_m …(5) 上記０次位相差Δφ０_mも、種々の処理に利用される。
相φ０_mと第（ｍ＋１）位相検出用データＤ_m+1の０次
位相φ０_m+1の０次位相差Δφ０_mが、次式により求め
られる。 Δφ０_m＝φ０_m+1−φ０_m …(5) 上記０次位相差Δφ０_mも、種々の処理に利用される。
【０００９】図８は、Ｚ(n)，Ｚcor1(n)，Ｚcor(n)の位
相を示す概念図である。１次位相の補正および０次位相
の補正を行うことにより、第１位相検出用データＤ_1〜
第Ｍ位相検出用データＤ_Mに対する位相エラーの影響を
なくすことが出来る。
相を示す概念図である。１次位相の補正および０次位相
の補正を行うことにより、第１位相検出用データＤ_1〜
第Ｍ位相検出用データＤ_Mに対する位相エラーの影響を
なくすことが出来る。
【００１０】
【発明が解決しようとする課題】上記(3)式により算出
した０次位相φ０_mは、上記(1)式により算出した１次
位相φ１_mと共に用いて、１次位相の補正と０次位相の
補正とを続けて行う場合には、問題がなかった。しか
し、上記(3)式により算出した０次位相φ０_mを単独で
補正等の処理に利用する場合には、正しい結果が得られ
ない問題点があった。この理由は、上記(3)式により算
出した０次位相φ０_mは、第１サンプリング点の複素ベ
クトルＺ(1)の位相そのものであり、第１サンプリング
点の複素ベクトルＺ(1)〜第Ｎサンプリング点の複素ベ
クトルＺ(N)の０次位相を正しく表していないからであ
る。
した０次位相φ０_mは、上記(1)式により算出した１次
位相φ１_mと共に用いて、１次位相の補正と０次位相の
補正とを続けて行う場合には、問題がなかった。しか
し、上記(3)式により算出した０次位相φ０_mを単独で
補正等の処理に利用する場合には、正しい結果が得られ
ない問題点があった。この理由は、上記(3)式により算
出した０次位相φ０_mは、第１サンプリング点の複素ベ
クトルＺ(1)の位相そのものであり、第１サンプリング
点の複素ベクトルＺ(1)〜第Ｎサンプリング点の複素ベ
クトルＺ(N)の０次位相を正しく表していないからであ
る。
【００１１】また、上記(5)式により算出した０次位相
差Δφ０_mは、ＥＰＩ法やＧＲＡＳＥ（GRadient And S
pin Echo）法において、第ｍ位相検出用データＤ_mに対
応するリード勾配の極性と第（ｍ＋１）位相検出用デー
タＤ_m+1に対応するリード勾配の極性とが反転している
場合、図９に示すようになり、正しい０次位相差を表さ
ない問題点があった。なお、図９で「正極性のリード勾
配」とは、図７の奇数のリード勾配ｒ１，ｒ３，…を指
す。また、図９で「負極性のリード勾配」とは、図７の
偶数のリード勾配ｒ２，ｒ４，…を指す。
差Δφ０_mは、ＥＰＩ法やＧＲＡＳＥ（GRadient And S
pin Echo）法において、第ｍ位相検出用データＤ_mに対
応するリード勾配の極性と第（ｍ＋１）位相検出用デー
タＤ_m+1に対応するリード勾配の極性とが反転している
場合、図９に示すようになり、正しい０次位相差を表さ
ない問題点があった。なお、図９で「正極性のリード勾
配」とは、図７の奇数のリード勾配ｒ１，ｒ３，…を指
す。また、図９で「負極性のリード勾配」とは、図７の
偶数のリード勾配ｒ２，ｒ４，…を指す。
【００１２】そこで、本発明の第１の目的は、第１サン
プリング点の複素ベクトルＺ(1)〜第Ｎサンプリング点
の複素ベクトルＺ(N)の０次位相を正しく表す０次位相
を求めることが出来る０次位相検出方法およびＭＲＩ装
置を提供することにある。また、本発明の第２の目的
は、連続する位相検出用エコーに対するリード勾配の極
性が反転する場合でも、正しい０次位相差を求めること
が出来る０次位相検出方法およびＭＲＩ装置を提供する
ことにある。
プリング点の複素ベクトルＺ(1)〜第Ｎサンプリング点
の複素ベクトルＺ(N)の０次位相を正しく表す０次位相
を求めることが出来る０次位相検出方法およびＭＲＩ装
置を提供することにある。また、本発明の第２の目的
は、連続する位相検出用エコーに対するリード勾配の極
性が反転する場合でも、正しい０次位相差を求めること
が出来る０次位相検出方法およびＭＲＩ装置を提供する
ことにある。
【００１３】
【課題を解決するための手段】第１の観点では、本発明
は、イメージング用パルスシーケンスから位相エンコー
ド勾配を省略したパルスシーケンスである位相検出用パ
ルスシーケンスにより集束させた位相検出用エコーから
収集した位相検出用データをフーリエ変換して得られる
第ｎサンプリング点の複素ベクトルＺ(n)を、 Ｚ(n)＝ｘ(n)＋ｉ・ｙ(n) …(6) で表すとき、 Ｚsum＝n=1ΣN｛ｘ(n)｝＋ｉ・n=1ΣN｛ｙ(n)｝ …(7) により、合成ベクトルＺsumを求め、その合成ベクトル
Ｚsumを用いて、 φ０＝arg｛Ｚsum｝ …(8) により、０次の位相φ０を求めることを特徴とする０次
位相検出方法を提供する。上記第１の観点による０次位
相検出方法では、第１サンプリング点の複素ベクトルＺ
(1)〜第Ｎサンプリング点の複素ベクトルＺ(N)の合成ベ
クトルＺsumの位相を０次位相とするから、従来技術の
ように第１サンプリング点の複素ベクトルＺ(1)の位相
を０次位相とする場合に較べて、第１サンプリング点の
複素ベクトルＺ(1)〜第Ｎサンプリング点の複素ベクト
ルＺ(N)の０次位相を正しく表す０次位相を求めること
が出来る。
は、イメージング用パルスシーケンスから位相エンコー
ド勾配を省略したパルスシーケンスである位相検出用パ
ルスシーケンスにより集束させた位相検出用エコーから
収集した位相検出用データをフーリエ変換して得られる
第ｎサンプリング点の複素ベクトルＺ(n)を、 Ｚ(n)＝ｘ(n)＋ｉ・ｙ(n) …(6) で表すとき、 Ｚsum＝n=1ΣN｛ｘ(n)｝＋ｉ・n=1ΣN｛ｙ(n)｝ …(7) により、合成ベクトルＺsumを求め、その合成ベクトル
Ｚsumを用いて、 φ０＝arg｛Ｚsum｝ …(8) により、０次の位相φ０を求めることを特徴とする０次
位相検出方法を提供する。上記第１の観点による０次位
相検出方法では、第１サンプリング点の複素ベクトルＺ
(1)〜第Ｎサンプリング点の複素ベクトルＺ(N)の合成ベ
クトルＺsumの位相を０次位相とするから、従来技術の
ように第１サンプリング点の複素ベクトルＺ(1)の位相
を０次位相とする場合に較べて、第１サンプリング点の
複素ベクトルＺ(1)〜第Ｎサンプリング点の複素ベクト
ルＺ(N)の０次位相を正しく表す０次位相を求めること
が出来る。
【００１４】第２の観点では、本発明は、上記構成の０
次位相検出方法において、第１の時に集束した位相検出
用エコーを基に得た合成ベクトルをＺsum_1とし、第２
の時に集束した位相検出用エコーを基に得た合成ベクト
ルをＺsum_2とするとき、 Δφ０＝arg｛Ｚsum_1｝−arg｛Ｚsum_2｝ …(9) により、０次の位相差Δφ０を求めることを特徴とする
０次位相検出方法を提供する。上記第２の観点による０
次位相検出方法では、第１の位相検出用エコーにかかる
複素ベクトルの合成ベクトルＺsumの位相arg｛Ｚsum_
1｝と第２の位相検出用エコーにかかる複素ベクトルの
合成ベクトルＺsumの位相arg｛Ｚsum_2｝の差を０次位
相差Δφ０とするから、第１の位相検出用エコーに対す
るリード勾配の極性と第２の位相検出用エコーに対する
リード勾配の極性とが反転する場合でも、正しい０次位
相差を求めることが出来る。
次位相検出方法において、第１の時に集束した位相検出
用エコーを基に得た合成ベクトルをＺsum_1とし、第２
の時に集束した位相検出用エコーを基に得た合成ベクト
ルをＺsum_2とするとき、 Δφ０＝arg｛Ｚsum_1｝−arg｛Ｚsum_2｝ …(9) により、０次の位相差Δφ０を求めることを特徴とする
０次位相検出方法を提供する。上記第２の観点による０
次位相検出方法では、第１の位相検出用エコーにかかる
複素ベクトルの合成ベクトルＺsumの位相arg｛Ｚsum_
1｝と第２の位相検出用エコーにかかる複素ベクトルの
合成ベクトルＺsumの位相arg｛Ｚsum_2｝の差を０次位
相差Δφ０とするから、第１の位相検出用エコーに対す
るリード勾配の極性と第２の位相検出用エコーに対する
リード勾配の極性とが反転する場合でも、正しい０次位
相差を求めることが出来る。
【００１５】第３の観点では、本発明は、上記構成の０
次位相検出方法において、第１の時に集束した位相検出
用エコーを基に得た合成ベクトルをＺsum_1とし、第２
の時に集束した位相検出用エコーを基に得た合成ベクト
ルをＺsum_2とするとき、 Δφ０＝arg｛Ｚsum_1／Ｚsum_2｝ …(10) により、０次の位相差Δφ０を求めることを特徴とする
０次位相検出方法を提供する。上記第３の観点による０
次位相検出方法は、上記第２の観点による０次位相検出
方法と等価であり、正しい０次位相差を求めることが出
来る。
次位相検出方法において、第１の時に集束した位相検出
用エコーを基に得た合成ベクトルをＺsum_1とし、第２
の時に集束した位相検出用エコーを基に得た合成ベクト
ルをＺsum_2とするとき、 Δφ０＝arg｛Ｚsum_1／Ｚsum_2｝ …(10) により、０次の位相差Δφ０を求めることを特徴とする
０次位相検出方法を提供する。上記第３の観点による０
次位相検出方法は、上記第２の観点による０次位相検出
方法と等価であり、正しい０次位相差を求めることが出
来る。
【００１６】第４の観点では、本発明は、上記構成の０
次位相検出方法において、前記イメージング用パルスシ
ーケンスが、リード勾配の極性を反転させてエコーを集
束させるパルスシーケンスであることを特徴とする０次
位相検出方法を提供する。上記第４の観点による０次位
相検出方法では、リード勾配の極性を反転させてエコー
を集束させるパルスシーケンスを対象とするが、このよ
うなパルスシーケンスは位相エラーに敏感であるため、
０次位相や０次位相差を正しく求めることは特に有用で
ある。
次位相検出方法において、前記イメージング用パルスシ
ーケンスが、リード勾配の極性を反転させてエコーを集
束させるパルスシーケンスであることを特徴とする０次
位相検出方法を提供する。上記第４の観点による０次位
相検出方法では、リード勾配の極性を反転させてエコー
を集束させるパルスシーケンスを対象とするが、このよ
うなパルスシーケンスは位相エラーに敏感であるため、
０次位相や０次位相差を正しく求めることは特に有用で
ある。
【００１７】第５の観点では、本発明は、上記構成の０
次位相検出方法において、前記イメージング用パルスシ
ーケンスが、ＥＰＩ法またはＧＲＡＳＥ法のパルスシー
ケンスであることを特徴とする０次位相検出方法を提供
する。上記第５の観点による０次位相検出方法では、Ｅ
ＰＩ法またはＧＲＡＳＥ法のパルスシーケンスを対象と
するが、このようなパルスシーケンスは位相エラーに敏
感であるため、０次位相や０次位相差を正しく求めるこ
とは特に有用である。
次位相検出方法において、前記イメージング用パルスシ
ーケンスが、ＥＰＩ法またはＧＲＡＳＥ法のパルスシー
ケンスであることを特徴とする０次位相検出方法を提供
する。上記第５の観点による０次位相検出方法では、Ｅ
ＰＩ法またはＧＲＡＳＥ法のパルスシーケンスを対象と
するが、このようなパルスシーケンスは位相エラーに敏
感であるため、０次位相や０次位相差を正しく求めるこ
とは特に有用である。
【００１８】第６の観点では、本発明は、上記構成の０
次位相検出方法において、前記第１の時の位相検出用エ
コーと前記第２の時の位相検出用エコーとが連続したエ
コーであることを特徴とする０次位相検出方法を提供す
る。上記第６の観点による０次位相検出方法では、連続
するエコーを対象とするためリード勾配の極性が反転し
ているが、この場合でも０次位相や０次位相差を正しく
求めることが出来る。
次位相検出方法において、前記第１の時の位相検出用エ
コーと前記第２の時の位相検出用エコーとが連続したエ
コーであることを特徴とする０次位相検出方法を提供す
る。上記第６の観点による０次位相検出方法では、連続
するエコーを対象とするためリード勾配の極性が反転し
ているが、この場合でも０次位相や０次位相差を正しく
求めることが出来る。
【００１９】第７の観点では、本発明は、上記構成の０
次位相検出方法において、イメージング用パルスシーケ
ンスによってイメージング用データを収集するスキャン
とは別のレファレンススキャンで位相検出用パルスシー
ケンスによって位相検出用データを収集することを特徴
とする０次位相検出方法を提供する。上記第７の観点に
よる０次位相検出方法では、イメージング用パルスシー
ケンスとは別に、レファレンススキャンにより位相検出
用データを収集するから、データ収集の時間的制約が少
なくなる。
次位相検出方法において、イメージング用パルスシーケ
ンスによってイメージング用データを収集するスキャン
とは別のレファレンススキャンで位相検出用パルスシー
ケンスによって位相検出用データを収集することを特徴
とする０次位相検出方法を提供する。上記第７の観点に
よる０次位相検出方法では、イメージング用パルスシー
ケンスとは別に、レファレンススキャンにより位相検出
用データを収集するから、データ収集の時間的制約が少
なくなる。
【００２０】第８の観点では、本発明は、上記構成の０
次位相検出方法において、イメージング用パルスシーケ
ンスの前に位相検出用パルスシーケンスを実行すること
を特徴とする０次位相検出方法を提供する。上記第８の
観点による０次位相検出方法では、イメージング用パル
スシーケンスの前に位相検出用パルスシーケンスを実行
するから、イメージング用パルスシーケンスでデータ収
集を行う際に、位相検出用パルスシーケンスで求めた０
次位相や０次位相差を利用可能となる。
次位相検出方法において、イメージング用パルスシーケ
ンスの前に位相検出用パルスシーケンスを実行すること
を特徴とする０次位相検出方法を提供する。上記第８の
観点による０次位相検出方法では、イメージング用パル
スシーケンスの前に位相検出用パルスシーケンスを実行
するから、イメージング用パルスシーケンスでデータ収
集を行う際に、位相検出用パルスシーケンスで求めた０
次位相や０次位相差を利用可能となる。
【００２１】第９の観点では、本発明は、ＲＦパルス送
信手段と、勾配パルス印加手段と、ＭＲ信号受信手段
と、前記各手段を制御してイメージング用パルスシーケ
ンスから位相エンコード勾配を省略したパルスシーケン
スである位相検出用パルスシーケンスにより位相検出用
エコーから位相検出用データを収集する位相検出用デー
タ収集手段と、前記位相検出用データをフーリエ変換し
て複素ベクトルを得るフーリエ変換手段と、第ｎサンプ
リング点の前記複素ベクトルＺ(n)を、 Ｚ(n)＝ｘ(n)＋ｉ・ｙ(n) …(6) で表すとき、 Ｚsum＝n=1ΣN｛ｘ(n)｝＋ｉ・n=1ΣN｛ｙ(n)｝ …(7) により合成ベクトルＺsumを求め、その合成ベクトルＺs
umを用いて、 φ０＝arg｛Ｚsum｝ …(8) により０次の位相φ０を求める０次位相算出手段とを具
備したことを特徴とするＭＲＩ装置を提供する。上記第
９の観点によるＭＲＩ装置では、上記第１の観点による
０次位相検出方法を好適に実施できる。
信手段と、勾配パルス印加手段と、ＭＲ信号受信手段
と、前記各手段を制御してイメージング用パルスシーケ
ンスから位相エンコード勾配を省略したパルスシーケン
スである位相検出用パルスシーケンスにより位相検出用
エコーから位相検出用データを収集する位相検出用デー
タ収集手段と、前記位相検出用データをフーリエ変換し
て複素ベクトルを得るフーリエ変換手段と、第ｎサンプ
リング点の前記複素ベクトルＺ(n)を、 Ｚ(n)＝ｘ(n)＋ｉ・ｙ(n) …(6) で表すとき、 Ｚsum＝n=1ΣN｛ｘ(n)｝＋ｉ・n=1ΣN｛ｙ(n)｝ …(7) により合成ベクトルＺsumを求め、その合成ベクトルＺs
umを用いて、 φ０＝arg｛Ｚsum｝ …(8) により０次の位相φ０を求める０次位相算出手段とを具
備したことを特徴とするＭＲＩ装置を提供する。上記第
９の観点によるＭＲＩ装置では、上記第１の観点による
０次位相検出方法を好適に実施できる。
【００２２】第１０の観点では、本発明は、上記構成の
ＭＲＩ装置において、第１の時に集束した位相検出用エ
コーを基に得た合成ベクトルをＺsum_1とし、第２の時
に集束した位相検出用エコーを基に得た合成ベクトルを
Ｚsum_2とするとき、 Δφ０＝arg｛Ｚsum_1｝−arg｛Ｚsum_2｝ …(9) により０次の位相差Δφ０を求める０次位相差算出手段
を具備したことを特徴とするＭＲＩ装置を提供する。上
記第１０の観点によるＭＲＩ装置では、上記第２の観点
による０次位相検出方法を好適に実施できる。
ＭＲＩ装置において、第１の時に集束した位相検出用エ
コーを基に得た合成ベクトルをＺsum_1とし、第２の時
に集束した位相検出用エコーを基に得た合成ベクトルを
Ｚsum_2とするとき、 Δφ０＝arg｛Ｚsum_1｝−arg｛Ｚsum_2｝ …(9) により０次の位相差Δφ０を求める０次位相差算出手段
を具備したことを特徴とするＭＲＩ装置を提供する。上
記第１０の観点によるＭＲＩ装置では、上記第２の観点
による０次位相検出方法を好適に実施できる。
【００２３】第１１の観点では、本発明は、上記構成の
ＭＲＩ装置において、第１の時に集束した位相検出用エ
コーを基に得た合成ベクトルをＺsum_1とし、第２の時
に集束した位相検出用エコーを基に得た合成ベクトルを
Ｚsum_2とするとき、 Δφ０＝arg｛Ｚsum_1／Ｚsum_2｝ …(10) により０次の位相差Δφ０を求める０次位相差算出手段
を具備したことを特徴とするＭＲＩ装置を提供する。上
記第１１の観点によるＭＲＩ装置では、上記第３の観点
による０次位相検出方法を好適に実施できる。
ＭＲＩ装置において、第１の時に集束した位相検出用エ
コーを基に得た合成ベクトルをＺsum_1とし、第２の時
に集束した位相検出用エコーを基に得た合成ベクトルを
Ｚsum_2とするとき、 Δφ０＝arg｛Ｚsum_1／Ｚsum_2｝ …(10) により０次の位相差Δφ０を求める０次位相差算出手段
を具備したことを特徴とするＭＲＩ装置を提供する。上
記第１１の観点によるＭＲＩ装置では、上記第３の観点
による０次位相検出方法を好適に実施できる。
【００２４】第１２の観点では、本発明は、上記構成の
ＭＲＩ装置において、前記イメージング用パルスシーケ
ンスが、リード勾配の極性を反転させてエコーを集束さ
せるパルスシーケンスであることを特徴とするＭＲＩ装
置を提供する。上記第１２の観点によるＭＲＩ装置で
は、上記第４の観点による０次位相検出方法を好適に実
施できる。
ＭＲＩ装置において、前記イメージング用パルスシーケ
ンスが、リード勾配の極性を反転させてエコーを集束さ
せるパルスシーケンスであることを特徴とするＭＲＩ装
置を提供する。上記第１２の観点によるＭＲＩ装置で
は、上記第４の観点による０次位相検出方法を好適に実
施できる。
【００２５】第１３の観点では、本発明は、上記構成の
ＭＲＩ装置において、前記イメージング用パルスシーケ
ンスが、ＥＰＩ法またはＧＲＡＳＥ法のパルスシーケン
スであることを特徴とするＭＲＩ装置を提供する。上記
第１３の観点によるＭＲＩ装置では、上記第５の観点に
よる０次位相検出方法を好適に実施できる。
ＭＲＩ装置において、前記イメージング用パルスシーケ
ンスが、ＥＰＩ法またはＧＲＡＳＥ法のパルスシーケン
スであることを特徴とするＭＲＩ装置を提供する。上記
第１３の観点によるＭＲＩ装置では、上記第５の観点に
よる０次位相検出方法を好適に実施できる。
【００２６】第１４の観点では、本発明は、上記構成の
ＭＲＩ装置において、前記第１の時の位相検出用エコー
と前記第２の時の位相検出用エコーとが連続したエコー
であることを特徴とするＭＲＩ装置を提供する。上記第
１４の観点によるＭＲＩ装置では、上記第６の観点によ
る０次位相検出方法を好適に実施できる。
ＭＲＩ装置において、前記第１の時の位相検出用エコー
と前記第２の時の位相検出用エコーとが連続したエコー
であることを特徴とするＭＲＩ装置を提供する。上記第
１４の観点によるＭＲＩ装置では、上記第６の観点によ
る０次位相検出方法を好適に実施できる。
【００２７】第１５の観点では、本発明は、上記構成の
ＭＲＩ装置において、イメージング用パルスシーケンス
によってイメージング用データを収集するスキャンとは
別のレファレンススキャンで位相検出用パルスシーケン
スによって位相検出用データを収集することを特徴とす
るＭＲＩ装置を提供する。上記第１５の観点によるＭＲ
Ｉ装置では、上記第７の観点による０次位相検出方法を
好適に実施できる。
ＭＲＩ装置において、イメージング用パルスシーケンス
によってイメージング用データを収集するスキャンとは
別のレファレンススキャンで位相検出用パルスシーケン
スによって位相検出用データを収集することを特徴とす
るＭＲＩ装置を提供する。上記第１５の観点によるＭＲ
Ｉ装置では、上記第７の観点による０次位相検出方法を
好適に実施できる。
【００２８】第１６の観点では、本発明は、上記構成の
ＭＲＩ装置において、イメージング用パルスシーケンス
の前に位相検出用パルスシーケンスを実行することを特
徴とするＭＲＩ装置を提供する。上記第１６の観点によ
るＭＲＩ装置では、上記第８の観点による０次位相検出
方法を好適に実施できる。
ＭＲＩ装置において、イメージング用パルスシーケンス
の前に位相検出用パルスシーケンスを実行することを特
徴とするＭＲＩ装置を提供する。上記第１６の観点によ
るＭＲＩ装置では、上記第８の観点による０次位相検出
方法を好適に実施できる。
【００２９】
【発明の実施の形態】以下、図に示す実施形態により本
発明をさらに詳しく説明する。なお、これにより本発明
が限定されるものではない。
発明をさらに詳しく説明する。なお、これにより本発明
が限定されるものではない。
【００３０】図１は、本発明の一実施形態にかかるＭＲ
Ｉ装置のブロック図である。このＭＲＩ装置１００にお
いて、マグネットアセンブリ１は、内部に被検体を挿入
するための空間部分（孔）を有し、この空間部分を取り
まくようにして、被検体に一定の主磁場を印加する永久
磁石１ｐと、勾配磁場を発生するための勾配磁場コイル
１ｇ（勾配磁場コイルはｘ軸，ｙ軸，ｚ軸の各コイルを
備えており、これらの組み合わせによりスライス軸，リ
ード軸，位相エンコード軸が決まる）と、被検体内の原
子核のスピンを励起するためのＲＦパルスを送信する送
信コイル１ｔと、被検体からのＭＲ信号を受信する受信
コイル１ｒ等が配置されている。勾配磁場コイル１ｇ，
送信コイル１ｔおよび受信コイル１ｒは、それぞれ勾配
磁場駆動回路３，ＲＦ電力増幅器４および前置増幅器５
に接続されている。なお、永久磁石１ｐの代わりに、超
電導磁石を用いてもよい。
Ｉ装置のブロック図である。このＭＲＩ装置１００にお
いて、マグネットアセンブリ１は、内部に被検体を挿入
するための空間部分（孔）を有し、この空間部分を取り
まくようにして、被検体に一定の主磁場を印加する永久
磁石１ｐと、勾配磁場を発生するための勾配磁場コイル
１ｇ（勾配磁場コイルはｘ軸，ｙ軸，ｚ軸の各コイルを
備えており、これらの組み合わせによりスライス軸，リ
ード軸，位相エンコード軸が決まる）と、被検体内の原
子核のスピンを励起するためのＲＦパルスを送信する送
信コイル１ｔと、被検体からのＭＲ信号を受信する受信
コイル１ｒ等が配置されている。勾配磁場コイル１ｇ，
送信コイル１ｔおよび受信コイル１ｒは、それぞれ勾配
磁場駆動回路３，ＲＦ電力増幅器４および前置増幅器５
に接続されている。なお、永久磁石１ｐの代わりに、超
電導磁石を用いてもよい。
【００３１】計算機７は、パルスシーケンスを作成し、
シーケンス記憶回路８に渡す。シーケンス記憶回路８
は、パルスシーケンスを記憶し、そのパルスシーケンス
に基づいて勾配磁場駆動回路３を操作し、マグネットア
センブリ１の勾配磁場コイル１ｇから勾配磁場を発生さ
せると共に、ゲート変調回路９を操作し、ＲＦ発振回路
１０の搬送波出力信号を所定タイミング・所定包絡線形
状のパルス状信号に変調し、それをＲＦパルスとしてＲ
Ｆ電力増幅器４に加え、ＲＦ電力増幅器４でパワー増幅
した後、前記マグネットアセンブリ１の送信コイル１ｔ
に印加する。前置増幅器５は、マグネットアセンブリ１
の受信コイル１ｒで受信したＭＲ信号を増幅し、位相検
波器１２に入力する。位相検波器１２は、ＲＦ発振回路
１０の搬送波出力信号を参照信号とし、ＭＲ信号を位相
検波して、ＡＤ変換器１１に与える。ＡＤ変換器１１
は、アナログ信号のＭＲ信号をディジタル信号のデータ
に変換し、計算機７に入力する。計算機７は、ＡＤ変換
器１１からデータを読み込み、位相検出処理，位相補正
処理，画像再構成処理等を実行し、イメージを作成す
る。このイメージは、表示装置６にて表示される。ま
た、計算機７は、操作卓１３から入力された情報を受け
取るなどの全体的な制御を受け持つ。
シーケンス記憶回路８に渡す。シーケンス記憶回路８
は、パルスシーケンスを記憶し、そのパルスシーケンス
に基づいて勾配磁場駆動回路３を操作し、マグネットア
センブリ１の勾配磁場コイル１ｇから勾配磁場を発生さ
せると共に、ゲート変調回路９を操作し、ＲＦ発振回路
１０の搬送波出力信号を所定タイミング・所定包絡線形
状のパルス状信号に変調し、それをＲＦパルスとしてＲ
Ｆ電力増幅器４に加え、ＲＦ電力増幅器４でパワー増幅
した後、前記マグネットアセンブリ１の送信コイル１ｔ
に印加する。前置増幅器５は、マグネットアセンブリ１
の受信コイル１ｒで受信したＭＲ信号を増幅し、位相検
波器１２に入力する。位相検波器１２は、ＲＦ発振回路
１０の搬送波出力信号を参照信号とし、ＭＲ信号を位相
検波して、ＡＤ変換器１１に与える。ＡＤ変換器１１
は、アナログ信号のＭＲ信号をディジタル信号のデータ
に変換し、計算機７に入力する。計算機７は、ＡＤ変換
器１１からデータを読み込み、位相検出処理，位相補正
処理，画像再構成処理等を実行し、イメージを作成す
る。このイメージは、表示装置６にて表示される。ま
た、計算機７は、操作卓１３から入力された情報を受け
取るなどの全体的な制御を受け持つ。
【００３２】図２は、本発明の一実施形態に係る位相検
出用パルスシーケンスのパルスシーケンス図である。こ
の位相検出用パルスシーケンスは、図５のマルチショッ
ト拡散強調ＥＰＩ法によるイメージング用パルスシーケ
ンスから位相エンコード勾配を省略したパルスシーケン
スである。すなわち、励起パルスRF90とスライス勾配SG
90とを印加する。次にＭＰＧパルスMPGを印加する。次
に反転ＲＦパルスRF180とスライス勾配SG180を印加す
る。次にＭＰＧパルスMPGを印加する。次に交互に正負
に反転するデータ収集用リード勾配ｒ１，…，ｒＭを連
続的に印加し且つ位相エンコード勾配は印加せず、順に
集束する第１位相検出用エコーＥ１から第Ｍ位相検出用
エコーＥＭを基にして第１位相検出用データＤ_1から第
Ｍ位相検出用データＤ_Mをそれぞれ収集する。
出用パルスシーケンスのパルスシーケンス図である。こ
の位相検出用パルスシーケンスは、図５のマルチショッ
ト拡散強調ＥＰＩ法によるイメージング用パルスシーケ
ンスから位相エンコード勾配を省略したパルスシーケン
スである。すなわち、励起パルスRF90とスライス勾配SG
90とを印加する。次にＭＰＧパルスMPGを印加する。次
に反転ＲＦパルスRF180とスライス勾配SG180を印加す
る。次にＭＰＧパルスMPGを印加する。次に交互に正負
に反転するデータ収集用リード勾配ｒ１，…，ｒＭを連
続的に印加し且つ位相エンコード勾配は印加せず、順に
集束する第１位相検出用エコーＥ１から第Ｍ位相検出用
エコーＥＭを基にして第１位相検出用データＤ_1から第
Ｍ位相検出用データＤ_Mをそれぞれ収集する。
【００３３】この位相検出用パルスシーケンスは、レフ
ァレンススキャンで実行される。このレファレンススキ
ャンの後、イメージング用パルスシーケンスによってイ
メージング用データを収集するスキャンを行う。
ァレンススキャンで実行される。このレファレンススキ
ャンの後、イメージング用パルスシーケンスによってイ
メージング用データを収集するスキャンを行う。
【００３４】図３は、本発明の一実施形態に係る０次位
相検出処理を示すフロー図である。ステップＳ１では、
エコー番号カウンタｍ＝１に初期化する。
相検出処理を示すフロー図である。ステップＳ１では、
エコー番号カウンタｍ＝１に初期化する。
【００３５】ステップＳ２では、第ｍ位相検出用データ
Ｄ_mに対してリード軸方向に１次元フーリエ変換を施
し、複素ベクトルＺ(n)_mを得る。ｎは、サンプリング
点番号で、１≦ｎ≦Ｎである。
Ｄ_mに対してリード軸方向に１次元フーリエ変換を施
し、複素ベクトルＺ(n)_mを得る。ｎは、サンプリング
点番号で、１≦ｎ≦Ｎである。
【００３６】ステップＳ３では、第ｎサンプリング点の
複素ベクトルＺ(n)_mを、 Ｚ(n)_m＝ｘ(n)_m＋ｉ・ｙ(n)_m …(6') で表すとき、 Ｚsum_m＝n=1ΣN｛ｘ(n)_m｝＋ｉ・n=1ΣN｛ｙ(n)_m｝ …(7') により、合成ベクトルＺsum_mを求める。ステップＳ４
では、合成ベクトルＺsum_mを用いて、 φ０_m＝arg｛Ｚsum_m｝ …(8') により、０次の位相φ０_mを求める。
複素ベクトルＺ(n)_mを、 Ｚ(n)_m＝ｘ(n)_m＋ｉ・ｙ(n)_m …(6') で表すとき、 Ｚsum_m＝n=1ΣN｛ｘ(n)_m｝＋ｉ・n=1ΣN｛ｙ(n)_m｝ …(7') により、合成ベクトルＺsum_mを求める。ステップＳ４
では、合成ベクトルＺsum_mを用いて、 φ０_m＝arg｛Ｚsum_m｝ …(8') により、０次の位相φ０_mを求める。
【００３７】ステップＳ５では、エコー番号カウンタｍ
を“１”だけインクリメントする。ステップＳ６では、
ｍ＝２なら前記ステップＳ２に戻り、ｍ≧３ならステッ
プＳ７へ進む。ステップＳ７では、 Δφ０_m-2＝arg｛Ｚsum_m-2｝−arg｛Ｚsum_m-1｝ …(9') により、０次の位相差Δφ０_m-2を求める。
を“１”だけインクリメントする。ステップＳ６では、
ｍ＝２なら前記ステップＳ２に戻り、ｍ≧３ならステッ
プＳ７へ進む。ステップＳ７では、 Δφ０_m-2＝arg｛Ｚsum_m-2｝−arg｛Ｚsum_m-1｝ …(9') により、０次の位相差Δφ０_m-2を求める。
【００３８】ステップＳ８では、ｍ＝３〜Ｍなら前記ス
テップＳ２に戻り、ｍ＝Ｍ＋１なら処理を終了する。
テップＳ２に戻り、ｍ＝Ｍ＋１なら処理を終了する。
【００３９】なお、前記ステップＳ７で、
Δφ０_m-2＝arg｛Ｚsum_m-1／Ｚsum_m-2｝ …(10')
により、０次の位相差Δφ０_m-2を求めてもよい。
【００４０】図４に示すように、得られた０次位相φ０
は、全サンプリング点の複素ベクトルの０次位相を表し
ている。従って、連続する位相検出用エコーに対するリ
ード勾配の極性が反転する場合でも、正しい０次位相差
Δφ０が得られる。よって、イメージング用データの補
正やそれ以外の種々の処理に利用されたときに正しい結
果が得られる。
は、全サンプリング点の複素ベクトルの０次位相を表し
ている。従って、連続する位相検出用エコーに対するリ
ード勾配の極性が反転する場合でも、正しい０次位相差
Δφ０が得られる。よって、イメージング用データの補
正やそれ以外の種々の処理に利用されたときに正しい結
果が得られる。
【００４１】
【発明の効果】本発明の０次位相検出方法およびＭＲＩ
装置によれば、ＭＲ信号をフーリエ変換して得た全サン
プリング点の複素ベクトルの合成ベクトルの位相を０次
位相とするから、従来技術のように第１サンプリング点
の複素ベクトルの位相を０次位相とする場合に較べて、
ＭＲ信号の０次位相を正しく表す０次位相を求めること
が出来る。
装置によれば、ＭＲ信号をフーリエ変換して得た全サン
プリング点の複素ベクトルの合成ベクトルの位相を０次
位相とするから、従来技術のように第１サンプリング点
の複素ベクトルの位相を０次位相とする場合に較べて、
ＭＲ信号の０次位相を正しく表す０次位相を求めること
が出来る。
【図１】本発明の一実施形態に係るＭＲＩ装置のブロッ
ク図である。
ク図である。
【図２】位相検出用パルスシーケンスのパルスシーケン
ス図である。
ス図である。
【図３】本発明の一実施形態に係る０次位相検出処理を
示すフロー図である。
示すフロー図である。
【図４】本発明に係る０次位相および０次位相差を示す
概念図である。
概念図である。
【図５】マルチショット拡散強調ＥＰＩ法によるイメー
ジング用パルスシーケンスのパルスシーケンス図であ
る。
ジング用パルスシーケンスのパルスシーケンス図であ
る。
【図６】イメージング用データの収集軌跡を示す概念図
である。
である。
【図７】位相検出用パルスシーケンスのパルスシーケン
ス図である。
ス図である。
【図８】従来の１次位相補正と０次位相補正の説明図で
ある。
ある。
【図９】従来の０次位相および０次位相差を示す概念図
である。
である。
１００ ＭＲＩ装置
１ マグネットアセンブリ
３ 勾配磁場駆動回路
７ 計算機
８ シーケンス記憶回路
フロントページの続き
(72)発明者 池崎 吉和
東京都日野市旭ケ丘４丁目７番地の127
ジーイー横河メディカルシステム株式会社
内
Ｆターム(参考） 4C096 AA17 AD06 AD07 AD13 BA41
BA42 BB32 DB10 DB20
Claims (16)
- 【請求項１】 イメージング用パルスシーケンスから位
相エンコード勾配を省略したパルスシーケンスである位
相検出用パルスシーケンスにより集束させた位相検出用
エコーから収集した位相検出用データをフーリエ変換し
て得られる第ｎサンプリング点の複素ベクトルＺ(n)
を、 Ｚ(n)＝ｘ(n)＋ｉ・ｙ(n) で表すとき、 Ｚsum＝n=1ΣN｛ｘ(n)｝＋ｉ・n=1ΣN｛ｙ(n)｝ により、合成ベクトルＺsumを求め、その合成ベクトル
Ｚsumを用いて、 φ０＝arg｛Ｚsum｝ により、０次の位相φ０を求めることを特徴とする０次
位相検出方法。 - 【請求項２】 請求項１に記載の０次位相検出方法にお
いて、第１の時に集束した位相検出用エコーを基に得た
合成ベクトルをＺsum_1とし、第２の時に集束した位相
検出用エコーを基に得た合成ベクトルをＺsum_2とする
とき、 Δφ０＝arg｛Ｚsum_1｝−arg｛Ｚsum_2｝ により、０次の位相差Δφ０を求めることを特徴とする
０次位相検出方法。 - 【請求項３】 請求項１に記載の０次位相検出方法にお
いて、第１の時に集束した位相検出用エコーを基に得た
合成ベクトルをＺsum_1とし、第２の時に集束した位相
検出用エコーを基に得た合成ベクトルをＺsum_2とする
とき、 Δφ０＝arg｛Ｚsum_1／Ｚsum_2｝ により、０次の位相差Δφ０を求めることを特徴とする
０次位相検出方法。 - 【請求項４】 請求項１から請求項３のいずれかに記載
の０次位相検出方法において、前記イメージング用パル
スシーケンスが、リード勾配の極性を反転させてエコー
を集束させるパルスシーケンスであることを特徴とする
０次位相検出方法。 - 【請求項５】 請求項４に記載の０次位相検出方法にお
いて、前記イメージング用パルスシーケンスが、ＥＰＩ
法またはＧＲＡＳＥ法のパルスシーケンスであることを
特徴とする０次位相検出方法。 - 【請求項６】 請求項４または請求項５に記載の０次位
相検出方法において、前記第１の時の位相検出用エコー
と前記第２の時の位相検出用エコーとが連続したエコー
であることを特徴とする０次位相検出方法。 - 【請求項７】 請求項１から請求項６のいずれかに記載
の０次位相検出方法において、イメージング用パルスシ
ーケンスによってイメージング用データを収集するスキ
ャンとは別のレファレンススキャンで位相検出用パルス
シーケンスによって位相検出用データを収集することを
特徴とする０次位相検出方法。 - 【請求項８】 請求項７に記載の０次位相検出方法にお
いて、イメージング用パルスシーケンスの前に位相検出
用パルスシーケンスを実行することを特徴とする０次位
相検出方法。 - 【請求項９】 ＲＦパルス送信手段と、勾配パルス印加
手段と、ＭＲ信号受信手段と、前記各手段を制御してイ
メージング用パルスシーケンスから位相エンコード勾配
を省略したパルスシーケンスである位相検出用パルスシ
ーケンスにより位相検出用エコーから位相検出用データ
を収集する位相検出用データ収集手段と、前記位相検出
用データをフーリエ変換して複素ベクトルを得るフーリ
エ変換手段と、第ｎサンプリング点の前記複素ベクトル
Ｚ(n)を、 Ｚ(n)＝ｘ(n)＋ｉ・ｙ(n) で表すとき、 Ｚsum＝n=1ΣN｛ｘ(n)｝＋ｉ・n=1ΣN｛ｙ(n)｝ により合成ベクトルＺsumを求め、その合成ベクトルＺs
umを用いて、 φ０＝arg｛Ｚsum｝ により０次の位相φ０を求める０次位相算出手段とを具
備したことを特徴とするＭＲＩ装置。 - 【請求項１０】 請求項９に記載のＭＲＩ装置におい
て、第１の時に集束した位相検出用エコーを基に得た合
成ベクトルをＺsum_1とし、第２の時に集束した位相検
出用エコーを基に得た合成ベクトルをＺsum_2とすると
き、 Δφ０＝arg｛Ｚsum_1｝−arg｛Ｚsum_2｝ により０次の位相差Δφ０を求める０次位相差算出手段
を具備したことを特徴とするＭＲＩ装置。 - 【請求項１１】 請求項９に記載のＭＲＩ装置におい
て、第１の時に集束した位相検出用エコーを基に得た合
成ベクトルをＺsum_1とし、第２の時に集束した位相検
出用エコーを基に得た合成ベクトルをＺsum_2とすると
き、 Δφ０＝arg｛Ｚsum_1／Ｚsum_2｝ により０次の位相差Δφ０を求める０次位相差算出手段
を具備したことを特徴とするＭＲＩ装置。 - 【請求項１２】 請求項９から請求項１１のいずれかに
記載のＭＲＩ装置において、前記イメージング用パルス
シーケンスが、リード勾配の極性を反転させてエコーを
集束させるパルスシーケンスであることを特徴とするＭ
ＲＩ装置。 - 【請求項１３】 請求項１２に記載のＭＲＩ装置におい
て、前記イメージング用パルスシーケンスが、ＥＰＩ法
またはＧＲＡＳＥ法のパルスシーケンスであることを特
徴とするＭＲＩ装置。 - 【請求項１４】 請求項１２または請求項１３に記載の
ＭＲＩ装置において、前記第１の時の位相検出用エコー
と前記第２の時の位相検出用エコーとが連続したエコー
であることを特徴とするＭＲＩ装置。 - 【請求項１５】 請求項９から請求項１４のいずれかに
記載のＭＲＩ装置において、イメージング用パルスシー
ケンスによってイメージング用データを収集するスキャ
ンとは別のレファレンススキャンで位相検出用パルスシ
ーケンスによって位相検出用データを収集することを特
徴とするＭＲＩ装置。 - 【請求項１６】 請求項１５に記載のＭＲＩ装置におい
て、イメージング用パルスシーケンスの前に位相検出用
パルスシーケンスを実行することを特徴とするＭＲＩ装
置。
Priority Applications (4)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2001328447A JP3720752B2 (ja) | 2001-10-26 | 2001-10-26 | ０次位相検出方法およびｍｒｉ装置 |
| US10/279,803 US20030083571A1 (en) | 2001-10-26 | 2002-10-24 | 0-Order phase detecting method and MRI system |
| CNB021470553A CN1231180C (zh) | 2001-10-26 | 2002-10-25 | 0阶相位检测方法和mri*** |
| KR1020020065514A KR100677017B1 (ko) | 2001-10-26 | 2002-10-25 | 0차 위상 검출 방법 및 ｍｒｉ 시스템 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2001328447A JP3720752B2 (ja) | 2001-10-26 | 2001-10-26 | ０次位相検出方法およびｍｒｉ装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2003135418A true JP2003135418A (ja) | 2003-05-13 |
| JP3720752B2 JP3720752B2 (ja) | 2005-11-30 |
Family
ID=19144505
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2001328447A Expired - Fee Related JP3720752B2 (ja) | 2001-10-26 | 2001-10-26 | ０次位相検出方法およびｍｒｉ装置 |
Country Status (4)
| Country | Link |
|---|---|
| US (1) | US20030083571A1 (ja) |
| JP (1) | JP3720752B2 (ja) |
| KR (1) | KR100677017B1 (ja) |
| CN (1) | CN1231180C (ja) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111239657A (zh) * | 2020-01-20 | 2020-06-05 | 上海东软医疗科技有限公司 | 谱图的相位校正方法、装置及设备 |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2006075380A (ja) * | 2004-09-10 | 2006-03-23 | Ge Medical Systems Global Technology Co Llc | Ｍｒ画像生成方法およびｍｒｉ装置 |
| JP4192139B2 (ja) * | 2004-10-27 | 2008-12-03 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | Ｓｓｆｐパルスシーケンスにおけるフェーズサイクリング法および磁気共鳴撮影装置 |
| US20120274322A1 (en) * | 2011-04-27 | 2012-11-01 | Sangwoo Lee | Magnetic resonance imaging apparatus |
Family Cites Families (14)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| NL9002842A (nl) * | 1990-12-21 | 1992-07-16 | Philips Nv | Magnetische resonantie werkwijze en inrichting ter reductie van beeldfouten in een magnetisch resonantie beeld. |
| JPH05245123A (ja) * | 1992-03-03 | 1993-09-24 | Toshiba Corp | 磁気共鳴イメージング装置 |
| DE4415393B4 (de) * | 1993-06-01 | 2010-10-14 | Siemens Ag | Verfahren zur Erzeugung von Bildern in einem Kernspintomographiegerät mit einer Spinecho-Pulssequenz |
| JP3505294B2 (ja) * | 1995-03-28 | 2004-03-08 | ジーイー横河メディカルシステム株式会社 | Ｍｒｉ装置 |
| JP3384944B2 (ja) * | 1996-07-11 | 2003-03-10 | ジーイー横河メディカルシステム株式会社 | Ｍｒｉ装置 |
| US5825185A (en) * | 1996-11-27 | 1998-10-20 | Picker International, Inc. | Method for magnetic resonance spin echo scan calibration and reconstruction |
| US5923168A (en) * | 1997-06-17 | 1999-07-13 | General Electric Company | Correction of artifacts caused by Maxwell terms in slice offset echo planar imaging |
| US6275038B1 (en) * | 1999-03-10 | 2001-08-14 | Paul R. Harvey | Real time magnetic field mapping using MRI |
| US6285187B1 (en) * | 1999-04-28 | 2001-09-04 | General Electric Company | Method and apparatus for reducing artifacts in echo planar imaging |
| US6259250B1 (en) * | 1999-04-28 | 2001-07-10 | General Electric Company | Method and apparatus for reducing artifacts in echo planar imaging |
| DE10015068C2 (de) * | 2000-03-25 | 2002-06-27 | Bruker Medical Gmbh | Verfahren zur Erzeugung von Bildern der magnetischen Resonanz |
| US6700374B1 (en) * | 2000-03-29 | 2004-03-02 | Koninklijke Philips Electronics, N.V. | EPI calibration method to minimize ghosting in reconstructed images |
| US6320380B1 (en) * | 2000-10-03 | 2001-11-20 | Marconi Medical Systems, Inc. | MRI method and apparatus for increasing the efficiency of echo lanar imaging and other late echo techniques |
| US6933720B2 (en) * | 2001-12-11 | 2005-08-23 | Toshiba America Mri, Inc. | Sequence preconditioning for ultra-fast magnetic resonance imaging |
-
2001
- 2001-10-26 JP JP2001328447A patent/JP3720752B2/ja not_active Expired - Fee Related
-
2002
- 2002-10-24 US US10/279,803 patent/US20030083571A1/en not_active Abandoned
- 2002-10-25 CN CNB021470553A patent/CN1231180C/zh not_active Expired - Fee Related
- 2002-10-25 KR KR1020020065514A patent/KR100677017B1/ko not_active IP Right Cessation
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111239657A (zh) * | 2020-01-20 | 2020-06-05 | 上海东软医疗科技有限公司 | 谱图的相位校正方法、装置及设备 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN1413555A (zh) | 2003-04-30 |
| KR20030035963A (ko) | 2003-05-09 |
| KR100677017B1 (ko) | 2007-01-31 |
| JP3720752B2 (ja) | 2005-11-30 |
| US20030083571A1 (en) | 2003-05-01 |
| CN1231180C (zh) | 2005-12-14 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP3878429B2 (ja) | Ｍｒｉ装置 | |
| EP1635184B1 (en) | MR image production method and MRI apparatus | |
| US6842001B2 (en) | MRI systems with parallel receivers for phase correction | |
| US6512372B1 (en) | MRI apparatus which has reduced ghost artifacts due to Maxwell term phase error | |
| US20050168221A1 (en) | MR imaging method and MRI system | |
| US5051699A (en) | Magnetic resonance imaging system | |
| US6285187B1 (en) | Method and apparatus for reducing artifacts in echo planar imaging | |
| JP3117670B2 (ja) | マルチスライスｍｒイメージング方法およびｍｒｉ装置 | |
| US20040064035A1 (en) | Magnetic resonance imaging apparatus | |
| JP2003210431A (ja) | Ｆｉｄ信号の位相ずれ測定方法、ｓｅ・ｓｔｅ信号の位相ずれ測定方法、ｍｒイメージング方法およびｍｒｉ装置 | |
| US7405566B2 (en) | RF pulse applying method and MRI apparatus | |
| JP2755125B2 (ja) | Ｍｒイメージング装置 | |
| JP2003010148A (ja) | Ｍｒイメージング方法およびｍｒｉ装置 | |
| US6528996B1 (en) | Diffusion-weighted imaging method and apparatus for fast pulse sequence with MPG pulses | |
| US6940280B2 (en) | Systems and methods for manipulating changes in phase encodes | |
| JP2007090001A (ja) | Ｍｒスキャン方法およびｍｒｉ装置 | |
| JP2009160051A (ja) | 磁気共鳴イメージング装置 | |
| JP3720752B2 (ja) | ０次位相検出方法およびｍｒｉ装置 | |
| US6777935B2 (en) | MR method for generating navigator pulses | |
| JPH0810239A (ja) | 拡散強調イメージング方法およびｍｒｉ装置 | |
| JP3576641B2 (ja) | Ｍｒｉ装置 | |
| JP3447099B2 (ja) | Ｍｒｉ装置 | |
| EP4043902A1 (en) | Dixon-type water/fat separation mr imaging | |
| JP3647444B2 (ja) | Ｍｒｉ装置 | |
| JPH04309330A (ja) | 核磁気共鳴撮像装置 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A625 | Written request for application examination (by other person) |
Free format text: JAPANESE INTERMEDIATE CODE: A625 Effective date: 20031211 |
|
| A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20050414 |
|
| TRDD | Decision of grant or rejection written | ||
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20050823 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20050908 |
|
| R150 | Certificate of patent or registration of utility model |
Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20080916 Year of fee payment: 3 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090916 Year of fee payment: 4 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100916 Year of fee payment: 5 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100916 Year of fee payment: 5 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100916 Year of fee payment: 5 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110916 Year of fee payment: 6 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110916 Year of fee payment: 6 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120916 Year of fee payment: 7 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120916 Year of fee payment: 7 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120916 Year of fee payment: 7 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120916 Year of fee payment: 7 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130916 Year of fee payment: 8 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| LAPS | Cancellation because of no payment of annual fees |