(57)【特許請求の範囲】
【請求項１】
ＭＲ画像再生装置であって、
該装置は、
ａ）核スピンを励起パルス（ＲＦ）により励起する励起フェーズと、
ｂ）核スピンを少なくとも１つの位相符号化勾配（ＧＰ１〜ＧＰ３）により位相符号化する符号化フェーズと、
ｃ）核共鳴信号（Ｓ）を読み出す読み出しフェーズ（ＡＱ）とからなる個別測定のシーケンスを実行する手段を有しており、
ここで前記核共鳴信号は先行する位相符号化勾配（ＧＰ１〜ＧＰ３）により検出された点にｋ空間内で割り当てられたものであり、
さらに読み出しフェーズ（ＡＱ）で得られた各信号（Ｓ）がｋ空間マトリクス（ＲＤ）にプロットされ、
前記フェーズａ）からｃ）による測定シーケンスが、ｋ空間を完全に走査するため異なる位相符号化によりｎ回実行される装置において、
第１の位相符号化勾配（ＧＰ１）が、時間的に関連する１群の個別測定の間、スイッチオンのままにされ、
その強度が測定シーケンス毎に変更される
ことを特徴とするＭＲ画像再生装置。
【請求項２】
励起パルス（ＲＦ）は広帯域である
請求項１記載の装置。
【請求項３】
核共鳴信号（Ｓ）はＦＩＤ信号である
請求項１または２記載の装置。
【請求項４】
第１の位相符号化勾配（ＧＰ１）は１群の個別測定内で線形に上昇する
請求項１から３までのいずれか１項記載の装置。
【請求項５】
各群の個別測定の間、第２の位相符号化勾配（ＧＰ２）が一定の大きさでスイッチオンされ、
該第２の位相符号化勾配は第１の位相符号化勾配（ＧＰ１）に対して垂直であり、
第２の位相符号化勾配（ＧＰ２）の大きさは、個別測定の群の間で、２次元ｋ空間のサンプリングが完全に行われるように連続的に切り替えられる
請求項１から４までのいずれか１項記載の装置。
【請求項６】
各群の個別測定の間、第３の位相符号化勾配（ＧＰ３）が一定の大きさでスイッチオンされ、
該第３の位相符号化勾配は、第１および第２の位相符号化勾配（ＧＰ１，ＧＰ２）に対して垂直であり、
第１および第２の位相符号化勾配（ＧＰ１，ＧＰ２）は、個別測定の群の間で、３次元ｋ空間のサンプリングが完全に行われるように連即的に切り替えられる
請求項５記載の装置。
【請求項７】
各群の個別測定の間、ｋ空間マトリクスの完全なラインが得られる
請求項１から６までのいずれか１項記載の装置。
【請求項８】
各読み出しフェーズで各共鳴信号が、励起に対してｙ個の固定時間間隔（ｔｉ）で読み出され、
所定の時間間隔（ｔｉ）に割り当てられた核共鳴信号がそれぞれｋ空間マトリクス（ＲＤｉ）にプロットされ、
ｙ個のｋ空間マトリクスが得られる
請求項１から７までのいずれか１項記載の装置。
【請求項９】
各ｋ空間マトリクス（ＲＤｉ）に対してフーリエ変換が実行され、
変換後に得られたデータセット（ＢＤｉ）が画像データセット（ＲＤ）の作成のために平均される
請求項８記載の装置。
【請求項１０】
平均の際に、変換後に得られたデータセットの補間が、種々異なるズーム係数の補正のために実行される
請求項９記載の装置。
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】
本発明は、以下のフェーズからなる個別測定のシーケンスによりＭＲ画像再生装置であって、
ａ）核スピンを励起パルスにより励起する励起フェーズ、
ｂ）核スピンを少なくとも１つの位相符号化勾配により位相符号化する符号化フェーズ、
ｃ）核共鳴信号を読み出す読み出しフェーズ、
ここで前記核共鳴信号は先行する位相符号化勾配（ＧＰ１〜ＧＰ３）により検出された点にｋ空間内で割り当てられたものであり、
読み出しフェーズ（ＡＱ）で得られた各信号（Ｓ）がｋ空間マトリクス（ＲＤ）にプロットされ、
前記ステップａ）からｃ）による測定シーケンスが、ｋ空間を完全に走査するため異なる位相符号化によりｎ回実行される装置に関する。
【０００２】
【従来の技術】
この種の技術は例えば、ドイツ連邦共和国特許第４２１９６１０号明細書およびドイツ連邦共和国特許出願公開第４２３２７３１号から公知である。
【０００３】
基本磁界での不均一性によって、通常のＭＲ画像再生シーケンスでは公知のように画像歪みが発生するようになる。現在使用されているパルスシーケンスは一般的に、いわゆる“スピンワープ”法に基づくものである。これは例えば、米国特許第４７７７０６０２５号明細書に記載されている。ここで核共振信号は読み出す前に少なくとも１つの方向で位相符号化され、読み出しの間に別の方向で読み出し勾配により周波数符号化される。基本磁界での不均一性は、位相符号化方向では比較的クリティカルでない。なぜなら、個々の位相符号化ステップ間での信号差だけが問題となるからである。しかし読み出し勾配の方向では、基本磁界の不均一性と読み出し勾配との重畳が歪みにつながる。層読み出しによる従来の方法では、高周波励起相中に別の歪みが発生する。なぜなら、磁界不均一性が座標に依存する付加的な層読み出し勾配として作用するからである。このことは、層面の湾曲と層厚の変化に結びつく。
【０００４】
したがって典型的には、核スピントモグラフィー装置用の磁石では、測定容積内で約１ＰＰＭまでの線形偏差しか許容されない。
【０００５】
この要求を満たすには大きなコストが必要である。とりわけ磁石を、本来の測定容積と比較して比較的に長く、ないし（磁極片磁石の場合）比較的大きな面積に構成しなければならない。これにより所要の均一性が一応達成される。
【０００６】
極端にエコー時間の短い固体材料の核スピントモグラフィーに対しては、例えばＮａｕｅｒｔ  Ａ．等著、ＳＭＲＭ  Ａｂｓｔｒａｃｔｓ，  １９９３，  １２１５頁から、励起後に位相符号化だけを実行し、励起後非常に高速にＦＩＤ信号を読み出し勾配の作用なしで検出することが公知である。その際、通常のスピンワープ法とは異なり、各信号によりｋ空間で１つの点しか得られない。したがって１２８×１２８画素の画像を形成するためには、それぞれ高周波励起および位相符号化した１２８×１２８回の個別測定が必要である。前記の選択励起の問題のためこれを省略し、その代わりに３次元位相符号化を実行すれば、例えば１２８×１２８×１２８＝２’０９７．１５７回の個別測定が必要になってしまう。勾配磁界は各個別測定毎にスイッチオンし、また再びスイッチオフしなければならないから、勾配コイルに対する所要のランプ時間のため測定時間が非常に長くなり、測定の効率が極端に悪くなる。
【０００７】
【発明が解決しようとする課題】
本発明の課題は、画像再生に対する時間のかからないＭＲ画像再生装置を提供することである。
【０００８】
【課題を解決するための手段】
上記課題は本発明により、第１の位相符号化勾配が、時間的に関連する１群の個別測定の間、スイッチオンのままにされ、
その強度が測定シーケンス毎に変更されることを特徴とするＭＲ画像再生装置を構成することにより解決される。
【０００９】
本発明の有利な構成は従属請求項に記載されている。
【００１０】
【発明の実施の形態】
以下本発明を、図１から図２２に基づいて詳細に説明する。
【００１１】
従来のスピンワープシーケンスが図１から図６に示されており、これらは問題提起の説明のためにだけ用いる。図示の実施例ではまず、層読み出し勾配ＧＳの作用下で周波数選択性の高周波パルスＲＦが照射される。これにより核スピンは検査装置の層でのみ励起される。引き続き、層読み出し勾配ＧＳの正の部分パルスに惹起されたデフェージングが負の部分パルスＧＳ−により再びリフェーズされる。さらに位相符号化勾配ＧＰが照射される。引き続き、実施例ではさらに負の読み出し勾配ＧＲ−がスイッチオンされる。
【００１２】
後続の読み出しフェーズの間は、単に正の読み出し勾配ＧＲ＋のみが照射される。発生するエコー信号Ｓは軸ＡＱに矢印で示されたように、Ｍ回走査され、得られたＭ個の測定値は図７のローデータマトリクスのラインにプロットされる。
【００１３】
図示のパルスシーケンスはＮ回、異なる値の位相符号化勾配ＧＰにより繰り返される。これにより全部でＮ列の測定マトリクスが得られる。ここで通常、位相符号化勾配はパルスシーケンス毎に同じステップで、正の最高値から負の最高値へ、および反対に連続的に切り替えられる。ローデータマトリクスＲＤは測定データ空間としてみなすことができ、この空間は実施例の２次元の場合では測定データ面とみなすことができる。この測定データ空間は核スピントモグラフィでは“ｋ空間”と称される。
【００１４】
画像再生に必要な、信号寄与値の空間的起源に関する情報は位相係数で符号化される。ここでは、デカルト座標ｘ，ｙ，ｚの位置空間とｋ空間との間に数学的にフーリエ変換を介した関係が存在する。次式が当てはまる。
【００１５】
【数１】
0003624021.tif 000002
【００１６】
ここでは以下の定義が当てはまる。
【００１７】
γ          ＝磁気回転比
ＧＲ（ｔ’）  ＝読み出し勾配の（ｘ方向での）瞬時値
ＧＰ（ｔ’）  ＝位相符号化勾配の（ｙ方向での）瞬時値
ρ（ｘｙ）  ＝核スピン密度
図７に示されたρデータマトリクスＲＤでは各ラインが個々の核共鳴信号に相応する。位相符号化勾配ＧＰがステップ毎に進められる場合、ｋ空間での走査は順次連続するラインで行われる。各個別測定の開始時に、各共鳴信号Ｓのそれぞれ前に位相符号化勾配ＧＰがスイッチオンされる。この位相符号化勾配の勾配振幅は部分シーケンス毎に連続してステップ的に上昇する。例えば各核共鳴信号を１２８個の測定点により走査し、１２８個の位相符号化ステップを実行する場合、１２８列と１２８行のρデータマトリクスが得られる。すなわち、１２８×１２８の測定値がｋ空間に得られる。したがって、図１から図６のパルスシーケンスで得られたアナログ測定信号はｋ空間のラスタにデジタル化される。
【００１８】
次に、ローデータマトリクスないしｋ空間マトリクスＲＤから２次元フーリエ変換により画像マトリクスが得られる。しかし冒頭にすでに述べたように、不均一な基本磁界では２種類の歪みが発生する。
【００１９】
−励起フェーズ中に、磁界の不均一性は層面の歪みおよび層厚の変化を惹起する。
【００２０】
−読み出しフェースでは、磁界の不均一性により読み出し勾配ＧＲ歩行での歪みが惹起される。読み出し勾配ＧＲによって理想的には線形の関係が、周波数符号化方向での位置ｘと核スピンの所属する共振周波数ｆとの間に成り立たなければならない。これは図８に示されている。しかし磁界の不均一性は図９に示すようにこの関係を非線形にしてしまう。したがって、線形読み出し勾配ＧＲに、基本磁界の不均一性を表す勾配が重畳される。
【００２１】
磁界は読み出し勾配の間、一定であると仮定すれば、不均一性は位相符号化方向において歪みにはつながらない。なぜなら、ここでは実質的に順次連続する位相符号ステップ間の信号差が重要だからである。
【００２２】
いわゆるＳＰＩシーケンス（信号ポイントイメージング）では読み出し勾配が省略され、位相符号化勾配が適用される。冒頭に述べた刊行物，Ｎａｕｅｒｔｈ等著、ＳＭＲＭ  Ａｂｓｔｒａｃｔｓ  １９９３，１２１５頁から公知のこのようなシーケンスが図１０から図１４に、同じように説明のために示されている。ここでは非選択性の高周波パルスＲＦの後に核共鳴信号の位相符号化が３方向で、３つの相互に垂直の位相符号化勾配ＧＰ１，ＧＰ２，ＧＰ３により行われる。励起後の設定された時間間隔ｔ１で、発生したＦＩＤ（Ｆｒｅｅ  Ｉｎｄｕｃｔｉｏｎ  Ｄｅｃａｙ）信号が読み出される。しかし一方で上記の従来の方法では、各信号により常に完全なラインがｋ空間マトリクスで得られる。ここでは、勾配ＧＰ１からＧＰ３により定義されたｋ空間に１つのポイントが得られるだけである。したがって、１２８×１２８×１２８のｋ空間ポイントを有するデータセットを得るために相応する数の個別測定が必要である。ここでは、勾配ＧＰ１からＧＰ３の相応する順次切り替えによって所要のすべてのｋ空間ポイントが検出される。
【００２３】
所要時間を許容できる範囲に保持するためには、個別測定の繰り返し時間をできるだけ短くするようにしなければならない。しかし各個別測定の際に、３つの勾配ＧＰ１からＧＰ３をそれぞれスイッチオン、オフしなければならない。勾配コイルは相当に大きなインダクタンスを有するから、切り替え過程を任意に高速にすることはできない。それどころか、スイッチオン−オフに対してかなりのランプ時間が必要であり、このランプ時間は繰り返し時間の短縮を困難にする。そのため実際の動作に対して許容できないような全体測定時間が発生する。
【００２４】
図１５から図１８の本発明の実施例に相応して、位相符号化勾配を各個別測定毎にスイッチオン−オフするのではなく、少なくとも１群の個別測定に対してスイッチオンすれば、測定時間が格段に短縮される。
【００２５】
図１５には高周波パルスＲＦのシーケンスが示されている。高周波パルス間には、それぞれＡＱで示されたデータ捕捉フェーズがある。位相符号化勾配ＧＰ２とＧＰ３は個別測定の図示の群では図１７と図１８に相応して一定に留まり、一方、位相符号化勾配ＧＰ１は負の値から正の値に線形に変化する。このようにしてＮ個の各共鳴信号が得られる。これら信号は位相符号化勾配ＧＰ１の方向ではすべて異なって位相符号化されている。したがって典型的にはｋ空間マトリクスの完全なラインが検出される。
【００２６】
全部でＭラインを有するｋ空間マトリクスの別のラインを検出するために、図示の測定はＭ回、値の異なる位相符号化勾配ＧＰ２により繰り返される。相応して、第３の次元を検出するための全体過程はＰ回、値の異なる位相符号化勾配ＧＰ３によって繰り返される。これにより全体で、図１５から図１８によるＭ・Ｐ個の測定が３次元ｋ空間を網羅するために必要である。
【００２７】
前記のパルスシーケンスの利点は、勾配パルスに対するランプが個別測定の間なくなり、これにより個別測定を格段に高速に繰り返すことができることである。
【００２８】
高周波パルスは“ハード”なパルス、すなわち広帯域パルスとして実施され、これにより同時にスイッチオンされる勾配が不所望なことに層選択するようなことが回避される。各読み出しフェーズの間、勾配ＧＰ１がスイッチオンされるが、しかしこのことは従来のデータ捕捉の場合のようには歪みにつながらない。なぜなら、各測定の際にただ１つの測定点だけが高周波パルスＲＦに対して固定間隔で記録され、読み出し勾配に起因する周波数依存性を備えた信号全部が読み出されるわけではないからである。
【００２９】
明瞭にするために、図１５と図１８のシーケンスの一部が図１９と図２０に拡大した時間尺度で示されている。ここでは、各信号Ｓｋが所属の励起パルスＲＦｋの後に固定の時間間隔ｔ１で読み出されることがわかる。位相符号化勾配ＧＰ１の方向での位相符号化は、励起と読み出し時点との間のこの勾配の時間積分により、すなわち時間間隔ｔ１にわたる時間積分により検出される。個の時間積分は図２０にハッチングで示されている。ここでは、位相符号化勾配ＧＰ１の方向での位相符号化は個別測定毎に上昇することがわかる。
【００３０】
前記の手法によって、従来の方法程度の測定時間が得られる。例えば、６４×１２８×１２８の測定点に対して１５分の測定時間を達成することができた。知覚し得るような画像アーチファクトのない許容できる磁界不均一性は、３００ｍｍの対象伸張にわたって、励起と読み出しとの間の時間間隔が５００μｓ、最大１０ｍＴ／ｍ・ｓの勾配装置の場合、０．３Ｔの主磁界の約１％である。ここでは１０μｓの持続時間の高周波励起パルスと４゜の励起角度が使用された。各共鳴信号に対するサンプリングレートは２００ｋＨｚであった。これにより２．５３ｍｍのＶｏｘｅｌに対して線形の解像度が達成された。
【００３１】
前記の手法では核共鳴信号からデータ点のみが走査されるから、全測定時間を基準にしたＳ／Ｎ比はまったく不利である。したがって有利な実施例では各個別測定で、発生した核共鳴信号を多数回、すなわちｙ個の異なる時間間隔ｔｉで励起後に走査する。これは図２２に矢印で示されている。引き続き各サンプリング時点に対して完全に別個のｋ空間マトリクスが作成され、全体でｙ個のｋ空間マトリクスＲＤ１からＲＤＹが得られる。別の処理が図２２のブロック回路図に示されている。各ｋ空間マトリクスＲＤ１からＲＤＹは別個にＦＦＴ（高速フーリエ変換）される。これにより相応する数の画像マトリクスＢＤ１’からＢＤＹ’が得られる。これら画像マトリクスは、励起に対して異なる時間間隔で走査したので位相符号化の強度が異なっており、そのため異なるズーム係数を有する。このことは図２２に画像マトリクスＢＤ２’とＢＤＹ’内の破線によって示されている。
【００３２】
まだ複素画像マトリクスであるＢＤ１’〜ＢＤＹ’から次に絶対値形成によって実画像マトリクスＢＤ１”〜ＢＤＹ”が得られる。択一的に絶対値形成の代わりに位相補正を行うこともできる。なぜなら、このステップでは最終的に位相誤差をが除去することが重要だからである。
【００３３】
画像マトリクスＢＤ１”〜ＢＤＹ”はいぜんとして異なるズーム係数を有する。このズーム係数は次のステップで相応の回転によって補正される。画像マトリクスにより設定されたラスタ内の画素を再生するために、ここでは補間が実行される。最後に画像データがｙ個の画像データマトリクスＢＤ１”’〜ＢＤＹ”’から平均され、これにより画像データマトリクスＢＤが得られる。このようにして得られた画像は明らかに改善されたＳ／Ｎ比を有する。正確に言えば、Ｓ／Ｎ比は平均された測定点の数のルートだけ改善される。しかもこの場合でも、測定時間は１つだけのデータ点を記録するＳＰＩ法に対して延長されない。
【図面の簡単な説明】
【図１】従来のスピンワープシーケンスも問題を説明するための図である。
【図２】従来のスピンワープシーケンスも問題を説明するための図である。
【図３】従来のスピンワープシーケンスも問題を説明するための図である。
【図４】従来のスピンワープシーケンスも問題を説明するための図である。
【図５】従来のスピンワープシーケンスも問題を説明するための図である。
【図６】従来のスピンワープシーケンスも問題を説明するための図である。
【図７】同様にｋ空間マトリクスの問題を説明するための図である。
【図８】磁界不均一性の作用を説明するための線図である。
【図９】磁界不均一性の作用を説明するための線図である。
【図１０】従来のＳＰＩシーケンスの問題を説明するための線図である。
【図１１】従来のＳＰＩシーケンスの問題を説明するための線図である。
【図１２】従来のＳＰＩシーケンスの問題を説明するための線図である。
【図１３】従来のＳＰＩシーケンスの問題を説明するための線図である。
【図１４】従来のＳＰＩシーケンスの問題を説明するための線図である。
【図１５】本発明の実施例を説明するための線図である。
【図１６】本発明の実施例を説明するための線図である。
【図１７】本発明の実施例を説明するための線図である。
【図１８】本発明の実施例を説明するための線図である。
【図１９】２つの個別測定の詳細を示す線図である。
【図２０】２つの個別測定の詳細を示す線図である。
【図２１】２つの個別測定の詳細を示す線図である。
【図２２】複数の画像データセットを平均するためのブロック図である。
【符号の説明】
ｆ  共振周波数
ＡＱ  データ捕捉フェーズ
ＢＤ  画像データマトリクス
ＧＰ  位相符号化勾配
ＧＳ  層選択勾配
ＧＳ−  部分パルス
ＲＤ  ローデータマトリクス
ＲＤ１〜ＲＤＹ  ｋ空間マトリクス
ＲＦ  高周波
Ｓ  核共鳴信号
ｔ１  時間間隔
ｘ  箇所
【図１】
0003624021.tif 000003
【図２】
0003624021.tif 000004
【図３】
0003624021.tif 000005
【図４】
0003624021.tif 000006
【図５】
0003624021.tif 000007
【図６】
0003624021.tif 000008
【図７】
0003624021.tif 000009
【図８】
0003624021.tif 000010
【図９】
0003624021.tif 000011
【図１０】
0003624021.tif 000012
【図１１】
0003624021.tif 000013
【図１２】
0003624021.tif 000014
【図１３】
0003624021.tif 000015
【図１４】
0003624021.tif 000016
【図１５】
0003624021.tif 000017
【図１６】
0003624021.tif 000018
【図１７】
0003624021.tif 000019
【図１８】
0003624021.tif 000020
【図１９】
0003624021.tif 000021
【図２０】
0003624021.tif 000022
【図２１】
0003624021.tif 000023
【図２２】
0003624021.tif 000024