(57)【特許請求の範囲】
【請求項１】
イメージング用パルスシーケンスから位相エンコード勾配を省略したパルスシーケンスである位相検出用パルスシーケンスにより集束させた位相検出用エコーから収集した位相検出用データをフーリエ変換して得られる第ｎサンプリング点の複素ベクトルＺ(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を求め、その合成ベクトルＺsumを用いて、
φ０＝arg｛Ｚsum｝
により０次の位相φ０を求める０次位相算出手段とを具備したことを特徴とするＭＲＩ装置。
【請求項１０】
請求項９に記載のＭＲＩ装置において、第１の時に集束した位相検出用エコーを基に得た合成ベクトルをＺsum_1とし、第２の時に集束した位相検出用エコーを基に得た合成ベクトルをＺsum_2とするとき、
Δφ０＝arg｛Ｚsum_1｝−arg｛Ｚsum_2｝
により０次の位相差Δφ０を求める０次位相差算出手段を具備したことを特徴とするＭＲＩ装置。
【請求項１１】
請求項９に記載のＭＲＩ装置において、第１の時に集束した位相検出用エコーを基に得た合成ベクトルをＺsum_1とし、第２の時に集束した位相検出用エコーを基に得た合成ベクトルをＺsum_2とするとき、
Δφ０＝arg｛Ｚsum_1／Ｚsum_2｝
により０次の位相差Δφ０を求める０次位相差算出手段を具備したことを特徴とするＭＲＩ装置。
【請求項１２】
請求項９から請求項１１のいずれかに記載のＭＲＩ装置において、前記イメージング用パルスシーケンスが、リード勾配の極性を反転させてエコーを集束させるパルスシーケンスであることを特徴とするＭＲＩ装置。
【請求項１３】
請求項１２に記載のＭＲＩ装置において、前記イメージング用パルスシーケンスが、ＥＰＩ法またはＧＲＡＳＥ法のパルスシーケンスであることを特徴とするＭＲＩ装置。
【請求項１４】
請求項１２または請求項１３に記載のＭＲＩ装置において、前記第１の時の位相検出用エコーと前記第２の時の位相検出用エコーとが連続したエコーであることを特徴とするＭＲＩ装置。
【請求項１５】
請求項９から請求項１４のいずれかに記載のＭＲＩ装置において、イメージング用パルスシーケンスによってイメージング用データを収集するスキャンとは別のレファレンススキャンで位相検出用パルスシーケンスによって位相検出用データを収集することを特徴とするＭＲＩ装置。
【請求項１６】
請求項１５に記載のＭＲＩ装置において、イメージング用パルスシーケンスの前に位相検出用パルスシーケンスを実行することを特徴とするＭＲＩ装置。
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】
本発明は、０次位相検出方法およびＭＲＩ（Magnetic Resonance Imaging）装置に関し、さらに詳しくは、ＭＲ信号の０次位相を正しく検出することが出来る０次位相検出方法およびＭＲＩ装置に関する。
【０００２】
【従来の技術】
図５に、マルチショット拡散強調ＥＰＩ（Echo Plannar Imaging）法によるイメージング用パルスシーケンスの基本例を示す。
このイメージング用パルスシーケンスでは、励起パルスRF90とスライス勾配SG90とを印加する。次にＭＰＧ（Motion Probing Gradient）パルスMPGを印加する。次に反転ＲＦパルスRF180とスライス勾配SG180を印加する。次にＭＰＧパルスMPGを印加する。次に位相エンコード勾配pdnを印加する。次に交互に正負に反転するデータ収集用リード勾配ｒ１，…，ｒｍを連続的に印加し且つ反転時に位相エンコード勾配ｐ２，…，ｐＭを印加し、第１エコーｅ１から第ＭエコーｅＭが順に集束するのとタイミングを合せてサンプリングし、各エコーｅ１，…，ｅＭに対応したイメージング用データＦ(n，1)，…，Ｆ(n，M)をそれぞれ収集する。これを位相エンコード勾配pdnの大きさを変えながらｎ＝１，…，Ｎについて繰り返して、ｋ空間を埋めるイメージング用データＦ(1，1)〜Ｆ(N，M)を収集する。これをＮショット・Ｍエコーという。また、ショットに対して実行時間順に付けた番号ｎをショット番号という。また、あるショットのエコー列のエコーに対して集束時間順に付けた番号をエコー番号という。
【０００３】
図６は、ｋ空間ＫＳにおけるイメージング用データＦ(1,1)〜Ｆ(N,M)の収集軌跡（trajectory）を示す模式図である。但し、Ｎ＝４，Ｍ＝４としている。
ｋ空間ＫＳを位相エンコード軸方向に第１行から第Ｎ×Ｍ行（図６では第１６行）まで分割したとき、第ｎショットの第ｍエコーで第（ｎ＋（ｍ−１）Ｎ）行のイメージング用データＦ(n,m)を収集するように、位相エンコードpdn，ｐ2，…，ｐＭを印加している。
【０００４】
さて、ＥＰＩ法は位相エラーに敏感であるため、ＭＲ信号の位相を検出し、位相補正する必要がある。
図７は、位相検出用パルスシーケンスの一例である。
この位相検出用パルスシーケンスは、図５のイメージング用パルスシーケンスから位相エンコード勾配を省略したパルスシーケンスである。
順に集束する第１位相検出用エコーＥ１から第Ｍ位相検出用エコーＥＭを基にして第１位相検出用データＤ_1から第Ｍ位相検出用データＤ_Mをそれぞれ収集する。
【０００５】
次に、第ｍ位相検出用データＤ_mをフーリエ変換し、複素ベクトルＺ(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は、イメージング用データの補正やそれ以外の種々の処理に利用される。
【０００８】
また、第ｍ位相検出用データＤ_mの０次位相φ０_mと第（ｍ＋１）位相検出用データＤ_m+1の０次位相φ０_m+1の０次位相差Δφ０_mが、次式により求められる。
Δφ０_m＝φ０_m+1−φ０_m      …(5)
上記０次位相差Δφ０_mも、種々の処理に利用される。
【０００９】
図８は、Ｚ(n)，Ｚcor1(n)，Ｚcor(n)の位相を示す概念図である。
１次位相の補正および０次位相の補正を行うことにより、第１位相検出用データＤ_1〜第Ｍ位相検出用データＤ_Mに対する位相エラーの影響をなくすことが出来る。
【００１０】
【発明が解決しようとする課題】
上記(3)式により算出した０次位相φ０_mは、上記(1)式により算出した１次位相φ１_mと共に用いて、１次位相の補正と０次位相の補正とを続けて行う場合には、問題がなかった。
しかし、上記(3)式により算出した０次位相φ０_mを単独で補正等の処理に利用する場合には、正しい結果が得られない問題点があった。この理由は、上記(3)式により算出した０次位相φ０_mは、第１サンプリング点の複素ベクトルＺ(1)の位相そのものであり、第１サンプリング点の複素ベクトルＺ(1)〜第Ｎサンプリング点の複素ベクトルＺ(N)の０次位相を正しく表していないからである。
【００１１】
また、上記(5)式により算出した０次位相差Δφ０_mは、ＥＰＩ法やＧＲＡＳＥ（GRadient And Spin Echo）法において、第ｍ位相検出用データＤ_mに対応するリード勾配の極性と第（ｍ＋１）位相検出用データＤ_m+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／Ｚsum_2｝      …(10)
により、０次の位相差Δφ０を求めることを特徴とする０次位相検出方法を提供する。
上記第３の観点による０次位相検出方法は、上記第２の観点による０次位相検出方法と等価であり、正しい０次位相差を求めることが出来る。
【００１６】
第４の観点では、本発明は、上記構成の０次位相検出方法において、前記イメージング用パルスシーケンスが、リード勾配の極性を反転させてエコーを集束させるパルスシーケンスであることを特徴とする０次位相検出方法を提供する。
上記第４の観点による０次位相検出方法では、リード勾配の極性を反転させてエコーを集束させるパルスシーケンスを対象とするが、このようなパルスシーケンスは位相エラーに敏感であるため、０次位相や０次位相差を正しく求めることは特に有用である。
【００１７】
第５の観点では、本発明は、上記構成の０次位相検出方法において、前記イメージング用パルスシーケンスが、ＥＰＩ法またはＧＲＡＳＥ法のパルスシーケンスであることを特徴とする０次位相検出方法を提供する。
上記第５の観点による０次位相検出方法では、ＥＰＩ法またはＧＲＡＳＥ法のパルスシーケンスを対象とするが、このようなパルスシーケンスは位相エラーに敏感であるため、０次位相や０次位相差を正しく求めることは特に有用である。
【００１８】
第６の観点では、本発明は、上記構成の０次位相検出方法において、前記第１の時の位相検出用エコーと前記第２の時の位相検出用エコーとが連続したエコーであることを特徴とする０次位相検出方法を提供する。
上記第６の観点による０次位相検出方法では、連続するエコーを対象とするためリード勾配の極性が反転しているが、この場合でも０次位相や０次位相差を正しく求めることが出来る。
【００１９】
第７の観点では、本発明は、上記構成の０次位相検出方法において、イメージング用パルスシーケンスによってイメージング用データを収集するスキャンとは別のレファレンススキャンで位相検出用パルスシーケンスによって位相検出用データを収集することを特徴とする０次位相検出方法を提供する。
上記第７の観点による０次位相検出方法では、イメージング用パルスシーケンスとは別に、レファレンススキャンにより位相検出用データを収集するから、データ収集の時間的制約が少なくなる。
【００２０】
第８の観点では、本発明は、上記構成の０次位相検出方法において、イメージング用パルスシーケンスの前に位相検出用パルスシーケンスを実行することを特徴とする０次位相検出方法を提供する。
上記第８の観点による０次位相検出方法では、イメージング用パルスシーケンスの前に位相検出用パルスシーケンスを実行するから、イメージング用パルスシーケンスでデータ収集を行う際に、位相検出用パルスシーケンスで求めた０次位相や０次位相差を利用可能となる。
【００２１】
第９の観点では、本発明は、ＲＦパルス送信手段と、勾配パルス印加手段と、ＭＲ信号受信手段と、前記各手段を制御してイメージング用パルスシーケンスから位相エンコード勾配を省略したパルスシーケンスである位相検出用パルスシーケンスにより位相検出用エコーから位相検出用データを収集する位相検出用データ収集手段と、前記位相検出用データをフーリエ変換して複素ベクトルを得るフーリエ変換手段と、第ｎサンプリング点の前記複素ベクトルＺ(n)を、
Ｚ(n)＝ｘ(n)＋ｉ・ｙ(n)      …(6)
で表すとき、
Ｚsum＝n=1ΣN｛ｘ(n)｝＋ｉ・n=1ΣN｛ｙ(n)｝      …(7)
により合成ベクトルＺsumを求め、その合成ベクトルＺsumを用いて、
φ０＝arg｛Ｚsum｝      …(8)
により０次の位相φ０を求める０次位相算出手段とを具備したことを特徴とするＭＲＩ装置を提供する。
上記第９の観点によるＭＲＩ装置では、上記第１の観点による０次位相検出方法を好適に実施できる。
【００２２】
第１０の観点では、本発明は、上記構成のＭＲＩ装置において、第１の時に集束した位相検出用エコーを基に得た合成ベクトルをＺsum_1とし、第２の時に集束した位相検出用エコーを基に得た合成ベクトルをＺsum_2とするとき、
Δφ０＝arg｛Ｚsum_1｝−arg｛Ｚsum_2｝      …(9)
により０次の位相差Δφ０を求める０次位相差算出手段を具備したことを特徴とするＭＲＩ装置を提供する。
上記第１０の観点によるＭＲＩ装置では、上記第２の観点による０次位相検出方法を好適に実施できる。
【００２３】
第１１の観点では、本発明は、上記構成のＭＲＩ装置において、第１の時に集束した位相検出用エコーを基に得た合成ベクトルをＺsum_1とし、第２の時に集束した位相検出用エコーを基に得た合成ベクトルをＺsum_2とするとき、
Δφ０＝arg｛Ｚsum_1／Ｚsum_2｝      …(10)
により０次の位相差Δφ０を求める０次位相差算出手段を具備したことを特徴とするＭＲＩ装置を提供する。
上記第１１の観点によるＭＲＩ装置では、上記第３の観点による０次位相検出方法を好適に実施できる。
【００２４】
第１２の観点では、本発明は、上記構成のＭＲＩ装置において、前記イメージング用パルスシーケンスが、リード勾配の極性を反転させてエコーを集束させるパルスシーケンスであることを特徴とするＭＲＩ装置を提供する。
上記第１２の観点によるＭＲＩ装置では、上記第４の観点による０次位相検出方法を好適に実施できる。
【００２５】
第１３の観点では、本発明は、上記構成のＭＲＩ装置において、前記イメージング用パルスシーケンスが、ＥＰＩ法またはＧＲＡＳＥ法のパルスシーケンスであることを特徴とするＭＲＩ装置を提供する。
上記第１３の観点によるＭＲＩ装置では、上記第５の観点による０次位相検出方法を好適に実施できる。
【００２６】
第１４の観点では、本発明は、上記構成のＭＲＩ装置において、前記第１の時の位相検出用エコーと前記第２の時の位相検出用エコーとが連続したエコーであることを特徴とするＭＲＩ装置を提供する。
上記第１４の観点によるＭＲＩ装置では、上記第６の観点による０次位相検出方法を好適に実施できる。
【００２７】
第１５の観点では、本発明は、上記構成のＭＲＩ装置において、イメージング用パルスシーケンスによってイメージング用データを収集するスキャンとは別のレファレンススキャンで位相検出用パルスシーケンスによって位相検出用データを収集することを特徴とするＭＲＩ装置を提供する。
上記第１５の観点によるＭＲＩ装置では、上記第７の観点による０次位相検出方法を好適に実施できる。
【００２８】
第１６の観点では、本発明は、上記構成のＭＲＩ装置において、イメージング用パルスシーケンスの前に位相検出用パルスシーケンスを実行することを特徴とするＭＲＩ装置を提供する。
上記第１６の観点によるＭＲＩ装置では、上記第８の観点による０次位相検出方法を好適に実施できる。
【００２９】
【発明の実施の形態】
以下、図に示す実施形態により本発明をさらに詳しく説明する。なお、これにより本発明が限定されるものではない。
【００３０】
図１は、本発明の一実施形態にかかるＭＲＩ装置のブロック図である。
このＭＲＩ装置１００において、マグネットアセンブリ１は、内部に被検体を挿入するための空間部分（孔）を有し、この空間部分を取りまくようにして、被検体に一定の主磁場を印加する永久磁石１ｐと、勾配磁場を発生するための勾配磁場コイル１ｇ（勾配磁場コイルはｘ軸，ｙ軸，ｚ軸の各コイルを備えており、これらの組み合わせによりスライス軸，リード軸，位相エンコード軸が決まる）と、被検体内の原子核のスピンを励起するためのＲＦパルスを送信する送信コイル１ｔと、被検体からのＭＲ信号を受信する受信コイル１ｒ等が配置されている。勾配磁場コイル１ｇ，送信コイル１ｔおよび受信コイル１ｒは、それぞれ勾配磁場駆動回路３，ＲＦ電力増幅器４および前置増幅器５に接続されている。なお、永久磁石１ｐの代わりに、超電導磁石を用いてもよい。
【００３１】
計算機７は、パルスシーケンスを作成し、シーケンス記憶回路８に渡す。
シーケンス記憶回路８は、パルスシーケンスを記憶し、そのパルスシーケンスに基づいて勾配磁場駆動回路３を操作し、マグネットアセンブリ１の勾配磁場コイル１ｇから勾配磁場を発生させると共に、ゲート変調回路９を操作し、ＲＦ発振回路１０の搬送波出力信号を所定タイミング・所定包絡線形状のパルス状信号に変調し、それをＲＦパルスとしてＲＦ電力増幅器４に加え、ＲＦ電力増幅器４でパワー増幅した後、前記マグネットアセンブリ１の送信コイル１ｔに印加する。
前置増幅器５は、マグネットアセンブリ１の受信コイル１ｒで受信したＭＲ信号を増幅し、位相検波器１２に入力する。位相検波器１２は、ＲＦ発振回路１０の搬送波出力信号を参照信号とし、ＭＲ信号を位相検波して、ＡＤ変換器１１に与える。ＡＤ変換器１１は、アナログ信号のＭＲ信号をディジタル信号のデータに変換し、計算機７に入力する。
計算機７は、ＡＤ変換器１１からデータを読み込み、位相検出処理，位相補正処理，画像再構成処理等を実行し、イメージを作成する。このイメージは、表示装置６にて表示される。
また、計算機７は、操作卓１３から入力された情報を受け取るなどの全体的な制御を受け持つ。
【００３２】
図２は、本発明の一実施形態に係る位相検出用パルスシーケンスのパルスシーケンス図である。
この位相検出用パルスシーケンスは、図５のマルチショット拡散強調ＥＰＩ法によるイメージング用パルスシーケンスから位相エンコード勾配を省略したパルスシーケンスである。
すなわち、励起パルスRF90とスライス勾配SG90とを印加する。次にＭＰＧパルスMPGを印加する。次に反転ＲＦパルスRF180とスライス勾配SG180を印加する。次にＭＰＧパルスMPGを印加する。次に交互に正負に反転するデータ収集用リード勾配ｒ１，…，ｒＭを連続的に印加し且つ位相エンコード勾配は印加せず、順に集束する第１位相検出用エコーＥ１から第Ｍ位相検出用エコーＥＭを基にして第１位相検出用データＤ_1から第Ｍ位相検出用データＤ_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を求める。
【００３７】
ステップＳ５では、エコー番号カウンタｍを“１”だけインクリメントする。ステップＳ６では、ｍ＝２なら前記ステップＳ２に戻り、ｍ≧３ならステップＳ７へ進む。
ステップＳ７では、
Δφ０_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を求めてもよい。
【００４０】
図４に示すように、得られた０次位相φ０は、全サンプリング点の複素ベクトルの０次位相を表している。従って、連続する位相検出用エコーに対するリード勾配の極性が反転する場合でも、正しい０次位相差Δφ０が得られる。よって、イメージング用データの補正やそれ以外の種々の処理に利用されたときに正しい結果が得られる。
【００４１】
【発明の効果】
本発明の０次位相検出方法およびＭＲＩ装置によれば、ＭＲ信号をフーリエ変換して得た全サンプリング点の複素ベクトルの合成ベクトルの位相を０次位相とするから、従来技術のように第１サンプリング点の複素ベクトルの位相を０次位相とする場合に較べて、ＭＲ信号の０次位相を正しく表す０次位相を求めることが出来る。
【図面の簡単な説明】
【図１】本発明の一実施形態に係るＭＲＩ装置のブロック図である。
【図２】位相検出用パルスシーケンスのパルスシーケンス図である。
【図３】本発明の一実施形態に係る０次位相検出処理を示すフロー図である。
【図４】本発明に係る０次位相および０次位相差を示す概念図である。
【図５】マルチショット拡散強調ＥＰＩ法によるイメージング用パルスシーケンスのパルスシーケンス図である。
【図６】イメージング用データの収集軌跡を示す概念図である。
【図７】位相検出用パルスシーケンスのパルスシーケンス図である。
【図８】従来の１次位相補正と０次位相補正の説明図である。
【図９】従来の０次位相および０次位相差を示す概念図である。
【符号の説明】
１００          ＭＲＩ装置
１              マグネットアセンブリ
３              勾配磁場駆動回路
７              計算機
８              シーケンス記憶回路
【図１】
0003720752.tif 000002
【図２】
0003720752.tif 000003
【図３】
0003720752.tif 000004
【図４】
0003720752.tif 000005
【図５】
0003720752.tif 000006
【図６】
0003720752.tif 000007
【図７】
0003720752.tif 000008
【図８】
0003720752.tif 000009
【図９】
0003720752.tif 000010