(57)【要約】【課題】空間分解能を低下させないで画像化に必要な磁気共鳴情報を短時間で取得可能な２次元磁気共鳴イメージング方法を提供する。【解決手段】画像を取得しようとする測定断面上に並ぶ複数の異なるラインについて、測定断面と非平行な平面を選択励起することによって各ラインを励起し、励起された各ラインの発する信号を観測することによって、励起面の重複を回避しながら複数の異なるラインの１次元磁気共鳴情報を取得する。そして、取得した複数のラインの１次元磁気共鳴情報から、前記測定断面の画像を形成する。【選択図】        図６−３
2006026311.tif 000002
【特許請求の範囲】
【請求項１】
測定対象物の断面の画像を取得するための２次元磁気共鳴イメージング方法であって、
画像を取得しようとする測定断面上に並ぶ複数の異なるラインについて、１次元測定法を用いて、各ラインの１次元磁気共鳴情報を取得する１次元磁気共鳴情報取得工程と、
前記複数のラインの１次元磁気共鳴情報から、前記測定断面の画像を形成する画像形成工程と、を含み、
前記１次元磁気共鳴情報取得工程において、
前記測定断面と非平行で、かつ、１次元磁気共鳴情報を取得しようとするライン上で前記測定断面と交わる平面を選択励起することによって、各ラインを励起し、励起された各ラインの発する信号を観測することによって、
励起面の重複を回避しながら複数の異なるラインの１次元磁気共鳴情報を測定することを特徴とする、２次元磁気共鳴イメージング方法。
【請求項２】
前記１次元磁気共鳴情報取得工程において、
前記測定断面と非平行で、かつ、１次元磁気共鳴情報を取得しようとするライン上で前記測定断面と交わる第１の平面を選択励起するとともに、
前記測定断面および前記第１の平面と非平行で、かつ、前記１次元磁気共鳴情報を取得しようとするライン上で前記測定断面および前記第１の平面と交わる第２の平面を選択励起し、
前記第１の平面と第２の平面との交線部分の磁気共鳴情報を、当該ラインの１次元磁気共鳴情報として取得することを特徴とする、請求項１に記載の２次元磁気共鳴イメージング方法。
【請求項３】
前記第１の平面を、９０°パルスおよび前記第１の平面に直交する勾配磁場によって励起するとともに、
前記第２の平面を、１８０°パルスおよび前記第２の平面に直交する勾配磁場によって励起することを特徴とする、請求項２に記載の２次元磁気共鳴イメージング方法。
【請求項４】
前記１次元磁気共鳴情報は、プロトン核の磁気共鳴情報であることを特徴とする、請求項１〜３のいずれか一項に記載の２次元磁気共鳴イメージング方法。
【請求項５】
前記１次元磁気共鳴情報取得工程において、マルチプレーン法により、前記各ラインの１次元磁気共鳴情報を取得することを特徴とする、請求項１〜４のいずれか一項に記載の２次元磁気共鳴イメージング方法。
【請求項６】
請求項１〜５のいずれか一項に記載の２次元磁気共鳴イメージング方法を用いて、前記測定対象物について複数の断面の画像を取得し、
取得した複数の断面の画像から、前記測定対象物の３次元磁気共鳴情報を形成することを特徴とする、３次元磁気共鳴イメージング方法。
【請求項７】
請求項１〜６のいずれか一項に記載された方法をコンピュータで実行させることを特徴とする、磁気共鳴イメージングプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、空間分解能を低下させることなく画像化に必要な磁気共鳴情報を短時間で取得することができる２次元磁気共鳴イメージング方法、および、これを用いた３次元磁気共鳴イメージング方法、磁気共鳴イメージングプログラムに関する。
【背景技術】
【０００２】
これまで主に医療分野で発展してきたＭＲＩ（Magnetic Resonance Imaging）は、生体に対して非侵襲的な測定方法であることを特徴とする。ＭＲＩは、多くの場合、生体内の水（プロトン）の二次元あるいは三次元的分布を画像化する技術である。そのため、もっぱら画像診断に利用されている。
【０００３】
生きたままの植物について、内部の水分布やその水移動を探るのに決定的な手法を欠いていた植物学の分野にも、この手法は応用が可能である。この際、空間分解能が高く、良好な画質の像を得るには、描画の対象となる範囲を表現するために必要な画素数を増す必要がある。ところが、画素数が増すにつれ、各画素に応じて膨大な測定時間を要し、このことが医療現場や、研究などの技術に応用する際の制限となることが多い。
【０００４】
従来の画像化法では２次元ＳＥ法（Spin Echo method）が一般的に用いられる。２次元ＳＥ法は画像化法の基本とも言うべき手法である。図２０に２次元ＳＥ法のパルスシーケンスを示す。２次元ＳＥ法は、まず９０°パルスを用いて磁化をＸＹ平面に倒す。このとき同時にＺ方向に勾配磁場をかけ、Ｚ方向のスライス部位を選択的に励起する。その後、ＸＹ平面に展開した磁化にＹ軸から勾配磁場をかける（位相エンコード）。このとき同時にＸ軸方向にリフォーカス勾配磁場をかけ、正常なエコー信号を生じさせる。その後、１８０°パルスにより展開した磁化を再結合させ、同時にＸ軸方向に位置情報を得る勾配磁場をかける。得られた情報をフーリエ変換することで２次元の画像を得ることができる。図２１に示すように、２次元ＳＥ法では、位置情報（画像）はＸとＹの２つの勾配磁場を用いた面として捉えられる。
【０００５】
２次元ＳＥ法において、測定時間（ＡＴ：acquisition time）は、一般的にＡＴ＝ＴＲ＊Ｎp＊Ａv で与えられる（ＴＲ：繰り返し時間（repetition time）、Ｎp：位相エンコードステップ数、Ａv：積算回数）。したがって、この３つのパラメータのうちのどれかを少なくすれば、撮像時間の短縮が可能である。しかし、繰り返し時間ＴＲを極端に短縮すると、測定対象によっては十分な信号を得づらくなる。また、位相エンコードステップ数Ｎpは画像分解能を決定する重要なパラメータであるから、Ｎpの短縮は結果として画像分解能の低下を意味する。積算回数Ａvは、Ｓ／Ｎ比を向上させるために同じ信号を複数回取得する場合の回数である。同じ信号を２回取れば単純に測定時間は２倍となる（通常積算回数は４〜１６回）。測定時間を短縮するため積算回数を少なくすることは一般によく取られる方法である。しかし、部位による磁化率の違いが顕著に見られる生体試料では、Ｓ／Ｎ比の低下が著しく、積算回数を少なくするとノイズの影響を強く受けてしまう。
【０００６】
画像化にかかる時間を短縮する方法はいくつか考案され、すでに実用化されているものもある（非特許文献１〜２１）。しかし、いずれも画像自体にかなりの妥協（劣化）を強いるものである。特に試料自体が小さいものの測定を試みた場合、細かい組織の解析をするうえで空間分解能の優劣が大きな問題となる。このことが短時間で起こる生体内の水動態の変化を調査する際の大きな障害となっていた。
【０００７】
【非特許文献１】Purcell EM, Torrey HC, Pound RV. Phys Rev 69:37, 1946
【非特許文献２】Bloch F, Hansen WW, Packard ME. Phys Rev 69:127, 1946
【非特許文献３】Hahn EL. Phys Rev 80:580, 1950
【非特許文献４】Lauterbur PC. Nature 242:190, 1973
【非特許文献５】Damadian R. US patent 3789832, 1972
【非特許文献６】Kumar A, Welti D, Ernst RR. J Magn Reson 18:69, 1975
【非特許文献７】Manfield P, Pykett LL. J Magn Reson 29:355, 1978
【非特許文献８】Hounsfield GN Brit J Radiol 46:1016, 1973
【非特許文献９】R. Damadian: Science, 171, 1151, 1971.
【非特許文献１０】W.A. EdelsteJn et al. : Phys. Mad. Biol.. 25. 751. 1980.
【非特許文献１１】AoN. Garroway, et al.: JoPhys. C, 7, L458, 1974.
【非特許文献１２】P. Mansfield, et al.. Phys. E, 9, 271, 1976.
【非特許文献１３】P. Mansfield, et al.: Brit. J. Radiol.,51, 921, 1978.
【非特許文献１４】P. Mansfield: J. Phys. C, 10, L55, 1977.
【非特許文献１５】W.S. Hinshaw: Phys. Letters, 48A, 87, 1974.
【非特許文献１６】W.S. Hinshaw: J. Appl. Phys., 47, 3709, 1976.
【非特許文献１７】W.S. Hinshaw, et al. : Nature, 270, 722, 1977.
【非特許文献１８】W.S. Hinshaw,et al.: Brit. jr. Radiol.,51, 273, 1978.
【非特許文献１９】W.S. Hinshaw,et al.: Brit. J. Radiol., 52, 36, 1979.
【非特許文献２０】P.Mansfield, P.R.Harvey.: Magnetic Resonance in Medicine. 29(6). 1993. 746-758.
【非特許文献２１】石田信昭・狩野広美・小川秀次郎  分析化学、Vol41、 561-566、 1992
【発明の開示】
【発明が解決しようとする課題】
【０００８】
本発明は、上記に鑑みてなされたものであって、空間分解能を低下させないで画像化に必要な磁気共鳴情報を短時間で取得することが可能な２次元磁気共鳴イメージング方法、３次元磁気共鳴イメージング方法、および、磁気共鳴イメージングプログラムを提供することを目的とする。
【課題を解決するための手段】
【０００９】
本発明者らは、上記課題を解決するため、鋭意研究の結果、従来から画像化前のキャリブレーション目的で用いられている１次元測定法に着目し、１次元測定法において励起面を調整して１次元磁気共鳴情報を取得することにより、画像を取得しようとする測定断面上に並ぶ複数のラインの１次元磁気共鳴情報を短時間で測定することが可能となることを見出し、本発明に想到した。
【００１０】
すなわち、本発明は下記の通りである。
〔１〕  測定対象物の断面の画像を取得するための２次元磁気共鳴イメージング方法であって、
画像を取得しようとする測定断面上に並ぶ複数の異なるラインについて、１次元測定法を用いて、各ラインの１次元磁気共鳴情報を取得する１次元磁気共鳴情報取得工程と、
前記複数のラインの１次元磁気共鳴情報から、前記測定断面の画像を形成する画像形成工程と、を含み、
前記１次元磁気共鳴情報取得工程において、
前記測定断面と非平行で、かつ、１次元磁気共鳴情報を取得しようとするライン上で前記測定断面と交わる平面を選択励起することによって、各ラインを励起し、励起された各ラインの発する信号を観測することによって、
励起面の重複を回避しながら複数の異なるラインの１次元磁気共鳴情報を測定することを特徴とする、２次元磁気共鳴イメージング方法。
【００１１】
〔２〕  前記１次元磁気共鳴情報取得工程において、
前記測定断面と非平行で、かつ、１次元磁気共鳴情報を取得しようとするライン上で前記測定断面と交わる第１の平面を選択励起するとともに、
前記測定断面および前記第１の平面と非平行で、かつ、前記１次元磁気共鳴情報を取得しようとするライン上で前記測定断面および前記第１の平面と交わる第２の平面を選択励起し、
前記第１の平面と第２の平面との交線部分の磁気共鳴情報を、当該ラインの１次元磁気共鳴情報として取得することを特徴とする、〔１〕に記載の２次元磁気共鳴イメージング方法。
【００１２】
〔３〕  前記第１の平面を、９０°パルスおよび前記第１の平面に直交する勾配磁場によって励起するととともに、
前記第２の平面を、１８０°パルスおよび前記第２の平面に直交する勾配磁場によって励起することを特徴とする、〔２〕に記載の２次元磁気共鳴イメージング方法。
【００１３】
〔４〕  前記１次元磁気共鳴情報は、プロトン核の磁気共鳴情報であることを特徴とする、〔１〕〜〔３〕のいずれか一項に記載の２次元磁気共鳴イメージング方法。
【００１４】
〔５〕  前記１次元磁気共鳴情報取得工程において、マルチプレーン法により、前記各ラインの１次元磁気共鳴情報を取得することを特徴とする、〔１〕〜〔４〕のいずれか一項に記載の２次元磁気共鳴イメージング方法。
【００１５】
〔６〕  〔１〕〜〔５〕のいずれか一項に記載の２次元磁気共鳴イメージング方法を用いて、前記測定対象物について複数の断面の画像を取得し、
取得した複数の断面の画像から、前記測定対象物の３次元磁気共鳴情報を形成することを特徴とする、３次元磁気共鳴イメージング方法。
【００１６】
〔７〕  〔１〕〜〔６〕のいずれか一項に記載された方法をコンピュータで実行させることを特徴とする、磁気共鳴イメージングプログラム。
【発明の効果】
【００１７】
本発明によれば、励起面の重複を回避しながら複数の異なるラインの１次元磁気共鳴情報を測定するため、Ｔ１緩和時間を待たずに、すぐに次のラインを励起することが可能となる。このため、空間分解能を低下させないで、撮像時間を大幅に短縮することが可能となる。
【００１８】
また、得られた複数の１次元磁気共鳴情報をそのまま利用して画像を形成することができるため、作像処理時間を短縮できる。
【発明を実施するための最良の形態】
【００１９】
以下に、本発明の実施形態について図面に基づいて詳細に説明する。まず、本発明の概要について説明し、その後、具体的な実施の形態を示しながら、本発明の構成および処理等について詳細に説明する。
【００２０】
（本発明の概要）
図１は、本発明の２次元磁気共鳴イメージング方法の概要を説明するためのフローチャートである。本発明の２次元磁気共鳴イメージング方法は、測定対象物の断面の画像を取得するための方法であって、本発明者らが開発した新規な１次元測定法を用いて、測定断面上に並ぶ複数の異なるライン上の１次元磁気共鳴情報を取得した後(ステップＳ１０)、取得した１次元磁気共鳴情報を用いて測定断面の画像を形成するものである(ステップＳ２０)。
【００２１】
本発明において、２次元磁気共鳴イメージングの対象となる測定対象物は、核磁気共鳴現象を起こさせることができる原子核を有する物体である。通常は、人、動物、植物等が測定対象物となる。
【００２２】
一般的に用いられている２次元ＳＥ法等の２次元測定法が、位置情報を面で捉えるのに対し、本発明で用いる１次元測定法は、位置情報を線（ライン）として捉える測定方法である。１次元測定法では、位相エンコードを用いないので、２次元測定法と比較して非常に短い測定時間で１次元磁気共鳴情報を得ることができる。例えば、観測核が1Ｈの場合、通常、数秒から数十秒で測定が可能である。
【００２３】
以下、代表的なスピン系列であるスピンエコー法を例にして、１次元測定法について説明する。図２は、１次元スピンエコー法のパルスシーケンスの一例である。スピンエコー法の場合、最初の９０°選択励起パルスとＺ軸勾配磁場の組合せで、Ｚ＝０の平面（第１の平面）が選択的に励起され、次の１８０°選択励起パルスとＹ軸勾配磁場でＹ＝０の平面（第２の平面）が選択的に励起される。従って、最終的にスピンエコーを結んで観測されるのは、Ｚ＝０（第１の平面）とＹ＝０（第２の平面）の交線部分となる。ここでスピンエコー信号観測時にＸ軸勾配磁場をかければ、Ｚ＝０とＹ＝０の交線部分のＸ方向の１次元磁気共鳴情報が得られる。選択励起される平面（以下、励起面と言う場合がある）は、実際には若干の厚みを有しており、この厚みは印加する勾配磁場の強さによって決まり、強い勾配磁場を用いるほど励起面を薄くできる。１次元スピンエコー法では、スピンエコー信号を発する領域は、図３に示すような、第１の平面と第２の平面とが重なりあった棒状の領域となる。空間分解能を向上させるためには、強い勾配磁場を用いることによって、この棒状の励起領域の断面を絞り込めばよい。この棒状領域の断面の大きさは、測定条件によって異なるが、通常、縦・横ともに１０μｍ〜５００μｍ程度である。
【００２４】
本発明者らが開発した１次元測定法の特徴は、測定断面上に並ぶ各ラインの１次元磁気共鳴情報を測定する際、測定断面と非平行で、かつ、１次元磁気共鳴情報を取得しようとするライン上で前記測定断面と交わる平面を選択励起することによって各ラインを励起することによって、励起面の重複を回避しながら複数の異なるラインの１次元磁気共鳴情報を測定する点にある。例えば、上述のスピンエコー法は、２つの平面（Ｚ＝０、Ｙ＝０）を選択励起する測定方法であるが、本発明者らが開発した１次元測定法をスピンエコー法に適用する場合は、測定断面と非平行な２つの異なる平面を選択励起する。
【００２５】
本発明者らが開発した１次元測定法と従来技術の違いについて、スピンエコー法によってＺ＝０平面の断面のイメージングを行う場合を例にして、説明する。１次元測定法によって１次元磁気共鳴情報を取得するラインは、図４に示すように、測定断面（Ｚ＝０）上に並ぶ複数のラインＸ1〜Ｘnである。図４では、ラインＸ1〜Ｘnは、測定断面上において、一定の間隔で平行に並んでいる。
【００２６】
従来の１次元測定法では、まずラインＸ1について１次元磁気共鳴情報を取得するため、図５−１に示すように、測定断面と同一面である第１の平面Ａ1（Ｚ＝０）、および、測定断面と直交する第２の平面Ｂ1（Ｙ＝ｎ1：ｎ1は実数）を選択励起し、第１の平面Ａ1と第２の平面Ｂ1の交線部分であるラインＸ1の発するスピンエコー信号を観測する。そして、次のラインＸ2についても同様に、図５−２に示すように、測定断面と同一面である第１の平面Ａ2（Ｚ＝０）、および、測定断面と直交する第２の平面Ｂ2（Ｙ＝ｎ2：ｎ2は実数）を選択励起して、第１の平面Ａ2と第２の平面Ｂ2の交線部分であるラインＸ2の発するスピンエコー信号を観測する。
【００２７】
しかし、ここで注意すべきことは、一度励起した面を再び励起するには、Ｔ１緩和時間を待つ必要があることである。従来の１次元測定法では、各ラインの１次元磁気共鳴情報を測定する度に、第１の平面であるＺ＝０平面を繰り返し励起していた。このため、次のラインを測定する前に、前回のラインの測定の際に励起されたＺ＝０平面上の原子核を基底の安定な状態に戻さなくてはならないため、毎回Ｔ１緩和時間待つことになり、このことが測定時間の長時間化の原因となっていた。
【００２８】
これに対し、本発明者らが開発した新規な１次元測定法によれば、測定断面と非平行な平面を選択励起することにより、励起面の重複を回避しながら複数の異なるラインの１次元磁気共鳴情報を測定することが可能である。このため、次のラインを測定する前にＴ１時間待つ必要がなくなり、測定時間を大幅に短縮することができる。
【００２９】
図６−１〜図６−３を参照しながら、本発明者らが開発した新規な１次元測定法について説明する。まず、ラインＸ1について１次元磁気共鳴情報を取得するため、図６−１に示すように、測定断面と非平行で、かつ、ラインＸ1上で測定断面と交わる第１の平面Ａ1、および、測定断面および第１の平面Ａ1と非平行で、かつ、ラインＸ1で測定断面および第１の平面Ａ1と交わる第２の平面Ｂ1を選択励起する。そして、第１の平面Ａ1と第２の平面Ｂ1の交線部分であるラインＸ1が発するスピンエコー信号を観測する。ここで、第１の平面Ａ1は、９０°パルスを照射する際、第１の平面Ａ1に対して垂直方向の勾配磁場を印加することによって選択励起することができる。また、第２の平面Ｂ1は、１８０°パルスを照射する際、第２の平面Ｂ1に対して垂直方向の勾配磁場を印加することによって選択励起することができる。
【００３０】
次いで、ラインＸ2についても１次元磁気共鳴情報を取得する。すなわち、図６−２に示すように、測定断面と非平行で、かつ、ラインＸ2上で測定断面と交わる第１の平面Ａ2、および、測定断面および第１の平面Ａ2と非平行で、かつ、ラインＸ2で測定断面および第１の平面Ａ2と交わる第２の平面Ｂ2を選択励起する。そして、第１の平面Ａ2と第２の平面Ｂ2の交線部分であるラインＸ2が発するスピンエコー信号を観測する。励起平面を、第１の平面Ａ1から第１の平面Ａ2に平行移動し、かつ、第２の平面Ｂ1からＢ2に平行移動するには、ＲＦ周波数を変えるだけでよい。すなわち、９０°パルスの周波数を、第１の平面Ａ2に対応する共鳴周波数に変えるとともに、１８０°パルスの周波数を、第２の平面Ｂ2に対応する共鳴周波数に変えることで、励起平面を平行移動することができる。
【００３１】
以降、同様にして、ＲＦ周波数を変えながら、図６−３に示すように、ラインＸnについても、測定断面と非平行で、かつ、ラインＸn上で測定断面と交わる第１の平面Ａnおよび第２の平面Ｂnを選択励起し、各ラインの発するスピンエコー信号を観測する。
【００３２】
このように、測定断面と非平行な平面を選択励起することによって、各ラインの測定時に生じる励起面の重複を回避することが可能となる。本発明者らが開発した１次元測定法によれば、異なるラインの測定において、励起面（第１の平面Ａと第２の平面Ｂ）が、まったく重ならないので、Ｔ１緩和時間を待たずに、すぐに次のラインを励起することが可能となり、測定時間を短縮することができる。
【００３３】
かくして得られた各ラインＸ1〜Ｘnのスピンエコー信号をフーリエ変換することによって、各ラインＸ1〜Ｘnの１次元磁気共鳴情報を得る。そして、各ラインＸ1〜Ｘnの１次元磁気共鳴情報を、測定断面上における各ラインの位置と対応付けて並べることにより、測定断面の画像を形成することができる。
【００３４】
本発明者らが開発した１次元測定法は、選択励起パルスを用いるパルスシーケンスには全て適用可能である。すなわち、線形傾斜磁場をかけた状態で周波数にある特定の幅を持たせたＲＦパルスを印加するパルスシーケンスであれば、特に限定なく適用することができる。したがって、上述のスピンエコー法の他、例えば、グラジェントエコー法、スティミュレーテッドエコー法等にも適用できる。グラジェントエコー法およびスティミュレーテッドエコー法は、スピンエコー法と異なり、９０°選択励起パルスのみのシークエンスであるが、９０°選択励起パルスで、測定断面と非平行な平面を選択励起し、励起された平面と測定断面の交線部分の発する信号を観測することにより、励起面の重複を回避しながら複数の異なるラインの１次元磁気共鳴情報を測定することができる。
【００３５】
また、本発明者らが開発した１次元測定法と、多断層法（マルチプレーン法）と組み合わせることによって、多数の測定ラインをほぼ同時に測定することも可能になり、測定時間をより高速化できる。図７にマルチプレーン法のパルスシーケンスの概念図を示す。マルチプレーン法によれば、１つ目のラインの励起が終了し、磁化の回復を待つＴＲの間に次のラインを励起することによって、ほぼ同時に多数のラインの１次元磁気共鳴情報を取得することができる。マルチプレーン法は、例えばＸＹ平面の２次元イメージをＺ軸方向に数枚から数十枚をほぼ同時に計測する際などに用いられていた手法であり、従来の１次元測定法では、前述のように励起面の重複があったため、マルチプレーン法を適用できなかった。しかし、本発明者らが開発した１次元測定法によれば、励起面の重複を回避することができるため、マルチプレーン法との併用が可能となる。
【００３６】
また、本発明の２次元磁気共鳴イメージング方法を用いて、測定対象物について複数の断面の画像を取得し、取得した複数の断面の画像を組み合わせて、測定対象物の３次元磁気共鳴イメージングを行うこともできる。
【００３７】
以上説明したように、本発明者らが開発した１次元測定法を用いることにより、撮像時間を大幅に短縮することができる。このため、従来は事実上測定困難であった核種を観測核とした２次元イメージングが可能となる。従来の方法では測定時間が長かったため、観測核は測定時間の比較的短いプロトン核が主流であった。しかし、本発明では高速測定が可能であるため、観測核は、核磁気共鳴現象を起こす原子核であれば特に限定されず、プロトン核の他、例えば、11Ｂ、13Ｃ、14Ｎ、15Ｎ、17Ｏ、19Ｆ、31Ｐ、33Ｓ、129Ｘｅ等を観測核とした測定が可能となる。このため、測定対象となる生体内物質の幅が広がり、各種疾患に対する診断・治療に役立つ。また、時間分解能を改善することができるため、例えば、生体内での代謝動態をリアルタイムでイメージングすることも可能となる。
【００３８】
また、本発明者らが開発した１次元測定法を用いることにより、画像解析に要する処理時間を短縮することが可能となる。従来の２次元測定法では、得られた信号からＺ＝０断面の画像を形成するに際し、観測した信号をＸ軸およびＹ軸ともにフーリエ変換する必要があったが、本発明によれば、各ラインについて観測した信号をＸ軸のみにフーリエ変換して得たデータを、そのまま画像解析に利用することができる。すなわち、本発明の方法は、２次元測定法と比較して、フーリエ変換処理が一軸少なくて済むので、作像処理に要する処理時間を短縮できる。
【００３９】
（実施の形態１）
実施の形態１は、スピンエコー法によって測定断面Ｚ＝０上の各ラインの１次元磁気共鳴情報を取得し、測定断面Ｚ＝０の２次元磁気共鳴イメージングを行う実施例である。
【００４０】
実施の形態１では、測定断面Ｚ＝０上において一定のピッチでＸ軸方向に平行に並ぶ各ラインｙ＝ａn（ａnは実数）について、１ラインづつ１次元磁気共鳴情報を測定していく。図８は、ラインｙ＝ａnを測定する場合に励起される平面の模式図である。実施の形態１は、測定断面Ｚ＝０に対して４５°傾いた平面Ｚ＋Ｙ＝ａn、および、９０°傾いた平面Ｙ＝ａnを励起することによって、ラインｙ＝ａnの１次元磁気共鳴情報を取得することを特徴とする。ここで、Ｚ＋Ｙ＝ａn平面は、９０°パルスとＺＹ勾配磁場の組合せによって励起することができ、Ｙ＝ａn平面は、１８０°パルスとＹ軸勾配磁場の組合せによって励起することができる。
【００４１】
１次元磁気共鳴情報は、公知のＭＲＩ装置を用いて取得することができる。１次元磁気共鳴測定に用いるＭＲＩ装置（図示せず）は、制御装置３０によって制御されている。図９は、ＭＲＩ装置を制御するための制御装置３０の構成を示す図である。制御装置３０は、ＣＰＵを含む制御部３０１と、メモリ３０２と、表示部３０３と、入力部３０４と、フーリエ変換部３０５と、作像部３０６と、ＣＤ−ＲＯＭドライブユニット３０７と、ディスクユニット３０８と、を備え、これらは、それぞれシステムバスＡを介して接続されている。
【００４２】
制御部３０１は、ＭＲＩ装置の動作を制御して、本発明者らが開発した１次元測定法に基づくパルスシーケンスを実行させる。メモリ３０２は、ＲＡＭ，ＲＯＭ等のメモリを含み、制御部３０１が実行すべきプログラム、および、測定の過程で得られた必要な情報を記憶する。表示部３０３は、ＣＲＴやＬＣＤ（液晶表示パネル）等で構成され、ユーザに対して、パルスシーケンスの設定画面、測定断面の画像等、各種画面を表示する。入力部３０４は、キーボード、マウス等で構成され、ユーザが、パルスシーケンスの測定条件等の各種情報の入力を行うために使用する。フーリエ変換部３０５は、ＭＲＩ装置で観測した各ラインの発する信号をフーリエ変換して各ラインの１次元磁気共鳴情報に変換する。作像部３０６は、各ラインの１次元磁気共鳴情報から測定断面の画像を形成する。
【００４３】
また、図示のＣＤ−ＲＯＭ４０には、本実施の形態の磁気共鳴イメージングプログラムを記述したプログラムが格納されている。まず、ＣＤ−ＲＯＭドライブユニット３０５にセットされたＣＤ−ＲＯＭ４０からプログラムがディスクユニット３０６にインストールされる。そして、ＭＲＩ制御装置を立ち上げるときにディスクユニット１０６から読み出されたプログラムが、メモリ３０２に格納される。この状態で、制御部３０１（ＣＰＵ）は、メモリ３０２に格納されたプログラムにしたがって、測定対象物に印加する静磁場および勾配磁場の制御、および、照射する電磁波の制御を行い、ＭＲＩ装置に本発明者らが開発した１次元測定法に基づくパルスシーケンスを実行させる。
【００４４】
なお、本実施の形態においては、ＣＤ−ＲＯＭ４０にて本発明の磁気共鳴イメージングプログラムを提供しているが、このプログラムの記録媒体は、これに限定されることなく、システムを構成するコンピュータに応じて、たとえば、フレキシブルディスク等の磁気ディスク，光磁気ディスク，磁気テープ等の他の記録媒体を用いることも可能である。
【００４５】
次に本実施の形態の処理について説明する。まず、撮影前に、ユーザは、入力部３０４から、各種測定条件を設定する。実施の形態１では、観測核はプロトン核であり、測定条件は下記の通りに設定した。なお、下記に示す測定条件の設定値は、あくまで一例であり、測定目的に応じて、設定値を適宜調整することができることは言うまでもない。
【００４６】
観測核の９０°パルス幅：48 [us]
パルスのアッテネータ値
周波数の変調を行うdomain：x_domain
周波数の変調を行うoffset：x_offset
選択励起パルスのパルス幅：1000[us]
選択励起パルスのアッテネータ値：23.1[dB]
１８０°パルスの定義：17.1[dB]
１８０°パルスの移動：可能に設定
ＲＦパルスの変調1ppmあたりの移動距離：0.25[mm]
選択励起パルスの関数：ガウス関数
スライスの高さ方向を指定・移動可能：5[ppm]
Ｘ軸磁場勾配のリフォーカス値（grad_1）：29.5[%]
Ｚ軸磁場勾配（選択励起）のリフォーカス値（grad_2）：11[%], 0[%]→60[%] : 0.1[%];
Ｚ軸磁場勾配のリフォーカス値（grad_3）：-6[%]
Ｙ軸磁場勾配（選択励起）のリフォーカス値（grad_4）：-3.5[%]
指定時間（grad_recover時間）：0.3[ms]
Ｔ２時間およびその制限値：0.2[ms], 0.1[ms] →300[s] : 0.1[ms]
測定開始までの待ち時間：1.0[s]
緩和時間およびその制限値：1.0[s], 0.1[s] →300[s] : 0.1[s]
【００４７】
測定条件を定義した後、ｙ＝ａ1、ｙ＝ａ2、ｙ＝ａ3、・・・ｙ＝ａnの順に、各ラインの１次元磁気共鳴情報を取得していく。図１０は、各ラインの１次元磁気共鳴情報測定のパルスシーケンスを実行するためのフローチャートである。まず、測定開始まで所定の待ち時間を待つ(ステップＳ１１１)。
【００４８】
次いで、９０°パルスとＺＹ勾配磁場でＺ＋Ｙ＝ａn平面を励起する。まず、Ｚ軸勾配磁場をgrad2で印加し、100μsec待った後(ステップＳ１２１)、Ｙ軸勾配磁場をgrad2で印加し、grad_recover待ち(ステップＳ１２２)、選択励起パルス（９０°パルス）を印加する(ステップＳ１２３)。ここで、選択励起パルスのＲＦ周波数は、Ｚ＋Ｙ＝ａn平面を励起するために所定の共鳴周波数に設定されている。ついで、Ｚ軸勾配磁場を切り、100μsec待った後(ステップＳ１２４)、Ｙ軸勾配磁場を切り、Ｔ２時間待つ(ステップＳ１２５)。
【００４９】
Ｚ＋Ｙ＝ａn平面を励起した後、正常なスピンエコー信号を得るため、リフォーカス勾配磁場をかける。Ｘ軸勾配磁場をgrad1で、Ｙ軸勾配磁場をgrad4で、Ｚ軸勾配磁場をgrad3で印加し、信号取り込み時間の半分の時間待つ(ステップＳ１３１)。次いで、Ｘ，Ｙ，Ｚ軸勾配磁場を切り、grad_recover時間待つ(ステップＳ１３２)。
【００５０】
リフォーカス勾配磁場をかけた後、１８０°パルスとＹ軸勾配磁場でＹ＝ａn平面を励起する。まず、Ｙ軸勾配磁場をgrad2で印加して、grad_recover時間待ち(ステップＳ１４１)、選択励起パルス（９０°パルス）を印加する(ステップＳ１４２)。ここで照射する選択励起パルスのＲＦ周波数は、Ｙ＝ａn平面を励起するために所定の共鳴周波数に設定されている。ついで、Ｙ軸勾配磁場を切り、Ｔ２時間待つ(ステップＳ１４３)。そして、grad_recover時間待つ(ステップＳ１４４)。
【００５１】
Ｚ＋Ｙ＝ａn平面およびＹ＝ａn平面を励起した後、２つの平面の交線部分となるｙ＝ａ1が発するスピンエコー信号を観測する。まず、データを取り込むためＸ軸勾配磁場をかけ、40μsec待つ(ステップＳ１５１)。そして、ｙ＝ａnの発するスピンエコー信号を取り込み、Ｘ軸勾配磁場をきる(ステップＳ１５２)。
【００５２】
このように、測定断面Ｚ＝０上の各ラインのスピンエコー信号を、ｙ＝ａ1、ｙ＝ａ2、ｙ＝ａ3、・・・ｙ＝ａnの順に取得していく。測定ラインを１ピッチずらして、次のラインを測定するには、選択励起のＲＦ周波数を変化させればよい。各ラインのピッチ（スライス厚）は、適宜設定することができるが、本実施の形態では、０．５ｍｍに設定されている。
【００５３】
フーリエ変換部３０５は、観測した各ラインのスピンエコー信号をフーリエ変換し、各ラインの１次元磁気共鳴情報として、メモリ３０２に保存する。
【００５４】
作像部３０６は、メモリ３０２から各ラインの１次元磁気共鳴情報を読み出し、各ラインの位置と各ラインの１次元磁気共鳴情報とを対応付けて並べることにより、測定断面（Ｚ＝０）の画像を形成する。得られた画像は、表示部３０３に表示され、ユーザが確認することができる。
【００５５】
以上説明したように、実施の形態１は、測定断面に対して、９０°パルスを４５°傾け、１８０°パルスを９０°傾けることにより、励起面の重複を回避しながら複数の異なるラインの１次元磁気共鳴情報測定を可能としている。
【００５６】
（実施の形態２）
実施の形態２は、実施の形態１の応用例であり、９０°パルスおよび１８０°パルスをそれぞれ４５°傾けて、Ｚ＝０の断面画像を取得する実施例である。なお、実施の形態２において、パルスシーケンス以外の構成および処理については、実施の形態１の２次元磁気共鳴イメージング方法と同様であるので、説明を省略する。
【００５７】
図１１は、実施の形態２で励起される平面の模式図である。実施の形態２では、測定断面Ｚ＝０に対して４５°傾いた平面Ｚ＋Ｙ＝ａn、および、４５°傾いた平面Ｚ−Ｙ＝ａnを励起することによって、測定断面Ｚ＝０上の各ラインｙ＝ａnの１次元磁気共鳴情報を取得することを特徴とする。ここで、Ｚ＋Ｙ＝ａn平面は、９０°パルスとＺＹ勾配磁場の組合せによって励起することができ、Ｙ＝ａn平面は、１８０°パルスとＺＹ勾配磁場の組合せによって励起することができる。
【００５８】
図１２は、実施の形態２において１次元磁気共鳴情報の測定に用いるパルスシーケンスである。また、図１３は、実施の形態２のパルスシーケンスを実行するためのフローチャートである。実施の形態１と同様にして、まず、測定開始まで所定の待ち時間を待ち、(ステップＳ１１１)、９０°パルスでＺ＋Ｙ＝ａn平面を励起し(ステップＳ１２１〜１２５)、リフォーカス勾配磁場をかける(ステップＳ１３１〜１３２)。
【００５９】
実施の形態２では、リフォーカス勾配磁場をかけた後、以下のように、１８０°パルスでＺ−Ｙ＝ａn平面を励起する。まず、Ｙ軸勾配磁場をgrad2で印加し、100μsec待ち(ステップＳ２０１)、Ｚ軸勾配磁場を−grad2で印加し、grad_recover時間待ち(ステップＳ２０２)、選択励起パルス（１８０°パルス）を印加する(ステップＳ２０３)。次いで、Ｙ軸勾配磁場を切り、100μsec待った後(ステップＳ２０４)、Ｚ軸勾配磁場を切って、Ｔ２時間待つ(ステップＳ２０５)。そして、grad_recover時間待つ(ステップＳ２０６)。
【００６０】
Ｚ＋Ｙ＝ａn平面およびＺ−Ｙ＝ａn平面を励起した後、実施の形態１と同様にして、２つの平面の交線部分となるｙ＝ａ1が発するスピンエコー信号を観測する(ステップＳ１５１、１５２)。
【００６１】
このようにして、測定断面Ｚ＝０上の各ラインの１次元磁気共鳴情報を、ｙ＝ａ1、ｙ＝ａ2、ｙ＝ａ3、・・・ｙ＝ａnの順に取得していくことができる。
【００６２】
以上説明したように、実施の形態２は、測定断面に対して、９０°パルスおよび１８０°パルスをそれぞれ４５°傾けることにより、励起面の重複を回避しながら複数の異なるラインの１次元磁気共鳴情報測定を可能としている。
【００６３】
（実施の形態３）
実施の形態３は、９０°パルスおよび１８０°パルスをそれぞれ４５°傾けて、Ｙ＝０の断面の画像を取得する実施例である。なお、実施の形態３において、測定断面およびパルスシーケンス以外の構成および処理については、実施の形態１の２次元磁気共鳴イメージング方法と同様であるので、説明を省略する。
【００６４】
図１４は、実施の形態３で励起される平面の模式図である。実施の形態３では、測定断面Ｙ＝０上において一定のピッチでＸ軸方向に平行に並ぶ各ラインｚ＝ａn（ａnは実数）について、測定断面Ｙ＝０に対してそれぞれ４５°傾いた平面Ｚ＋Ｙ＝ａnおよび平面Ｚ−Ｙ＝ａnを励起することによって、１次元磁気共鳴情報を取得することを特徴とする。ここで、Ｚ＋Ｙ＝ａn平面は、９０°パルスとＺＹ勾配磁場の組合せによって励起することができ、Ｚ−Ｙ＝ａn平面は、１８０°パルスとＺＹ勾配磁場の組合せによって励起することができる。
【００６５】
図１５は、実施の形態３において１次元磁気共鳴情報の測定に用いるパルスシーケンスである。また、図１６は、実施の形態３のパルスシーケンスを実行するためのフローチャートである。
【００６６】
実施の形態１と同様にして、まず、測定開始まで所定の待ち時間を待ち、(ステップＳ１１１)、９０°パルスでＺ＋Ｙ＝ａn平面を励起し(ステップＳ１２１〜１２５)、リフォーカス勾配磁場をかける(ステップＳ１３１〜１３２)。
【００６７】
実施の形態３では、リフォーカス勾配磁場をかけた後、以下のように、１８０°パルスでＺ−Ｙ＝ａn平面を励起する。まず、Ｙ軸勾配磁場を-grad2で印加し、100μsec待ち(ステップＳ３０１)、Ｚ軸勾配磁場をgrad2で印加し、grad_recover時間待ち(ステップＳ３０２)、選択励起パルス（１８０°パルス）を印加する(ステップＳ３０３)。次いで、Ｙ軸勾配磁場を切り、100μsec待った後(ステップＳ３０４)、Ｚ軸勾配磁場を切って、Ｔ２時間待つ(ステップＳ３０５)。そして、grad_recover時間待つ(ステップＳ３０６)。
【００６８】
Ｚ＋Ｙ＝ａn平面およびＺ−Ｙ＝ａn平面を励起した後、実施の形態１と同様にして、２つの平面の交線部分となるｚ＝ａnが発するスピンエコー信号を観測する(ステップＳ１５１、１５２)。
【００６９】
このようにして、測定断面Ｙ＝０上の各ラインの１次元磁気共鳴情報を、ｚ＝ａ1、ｚ＝ａ2、ｚ＝ａ3、・・・ｚ＝ａnの順に取得していくことができる。
【００７０】
以上説明したように、実施の形態３は、測定断面Ｙ＝０の場合についても、励起面の重複を回避しながら複数の異なるラインの１次元磁気共鳴情報測定を可能としている。
【００７１】
（実施の形態４）
実施の形態４は、実施の形態２の応用例であり、マルチプレーン法により、測定断面Ｚ＝０上の各ラインｙ＝ａn（ａnは実数）の１次元磁気共鳴情報をほぼ同時に測定する実施例である。実施の形態４において、マルチプレーン法を用いる以外の構成および処理については、実施の形態２の２次元磁気共鳴イメージング方法と同様であるので、説明を省略する。
【００７２】
図１７は、実施の形態４のパルスシーケンスを実行するためのフローチャートである。実施の形態４では、マルチプレーン法を適用することにより、測定断面Ｚ＝０上の各ラインｙ＝ａ1、ｙ＝ａ2、ｙ＝ａ3、・・・ｙ＝ａnの１次元磁気共鳴情報を、連続的に測定することを特徴とする。
【００７３】
まず、測定開始まで所定の待ち時間を待つ(ステップＳ１１１)。
【００７４】
実施の形態４では、周波数を変調しながら選択励起パルスを印加して、Ｚ＋Ｙ＝ａn平面を連続的に励起する。まず、Ｚ軸勾配磁場をgrad2で印加し、100μsec待った後(ステップＳ４０１)、Ｙ軸勾配磁場をgrad2で印加し、grad_recover待ち(ステップＳ４０２)、周波数を変調しながら選択励起パルス（９０°パルス）を印加する(ステップＳ４０３)。周波数を変調しながら選択励起パルスを印加することにより、Ｚ＋Ｙ＝ａ1、Ｚ＋Ｙ＝ａ2、・・・、Ｚ＋Ｙ＝ａn平面をほぼ同時に励起することができる。ここでは、例えば、周波数を−20ppmから1ppmステップで変調するよう設定する。ついで、Ｚ軸勾配磁場を切り、100μsec待った後(ステップＳ４０４)、Ｙ軸勾配磁場を切り、Ｔ２時間待つ(ステップＳ４０５)。
【００７５】
Ｚ＋Ｙ＝ａn平面を連続的に励起した後、実施の形態２と同様に、リフォーカス勾配磁場をかけ(ステップＳ１３１〜１３２)、１８０°パルスでＺ−Ｙ＝ａn平面を励起し(ステップＳ２０１〜２０６)、Ｚ＋Ｙ＝ａn平面およびＺ−Ｙ＝ａn平面の交線部分となるｙ＝ａnが発するスピンエコー信号を観測する(ステップＳ１５１、１５２)。
【００７６】
このように、マルチプレーン法を用いたパルスシークエンスにより、測定断面Ｚ＝０上の各ラインｙ＝ａ1、ｙ＝ａ2、ｙ＝ａ3、・・・ｙ＝ａnの１次元磁気共鳴情報を、一気に取得することができる。
【００７７】
以上説明したように、実施の形態４は、マルチプレーン法を適用することにより、測定断面上の複数の異なるラインの１次元磁気共鳴情報測定をほぼ同時に測定することができるため、測定時間をより短縮することができる。
【００７８】
（実施の形態５）
実施の形態５は、実施の形態３の応用例であり、マルチプレーン法により、測定断面Ｙ＝０上の各ラインｚ＝ａn（ａnは実数）の１次元磁気共鳴情報をほぼ同時に測定する実施例である。実施の形態５において、マルチプレーン法を用いる以外の構成および処理については、実施の形態３の２次元磁気共鳴イメージング方法と同様であるので、説明を省略する。
【００７９】
図１８は、実施の形態５のパルスシーケンスを実行するためのフローチャートである。実施の形態５では、マルチプレーン法を適用することにより、測定断面Ｙ＝０上の各ラインｚ＝ａ1、ｚ＝ａ2、ｚ＝ａ3、・・・ｚ＝ａnの１次元磁気共鳴情報を、連続的に測定することを特徴とする。
【００８０】
まず、測定開始まで所定の待ち時間を待つ(ステップＳ１１１)。
【００８１】
そして、実施の形態４と同様に、周波数を変調しながら選択励起パルスを印加して、Ｚ＋Ｙ＝ａn平面を連続的に励起する(ステップＳ４０１〜４０５)。
【００８２】
Ｚ＋Ｙ＝ａn平面を連続的に励起した後、実施の形態３と同様に、リフォーカス勾配磁場をかけ(ステップＳ１３１〜１３２)、１８０°パルスでＺ−Ｙ＝ａn平面を励起し(ステップＳ３０１〜３０６)、Ｚ＋Ｙ＝ａn平面およびＺ−Ｙ＝ａn平面の交線部分となるｙ＝ａnが発するスピンエコー信号を観測する(ステップＳ１５１、１５２)。
【００８３】
このように、マルチプレーン法を用いたパルスシークエンスにより、測定断面Ｙ＝０上の各ラインｚ＝ａ1、ｚ＝ａ2、ｚ＝ａ3、・・・ｚ＝ａnの１次元磁気共鳴情報を、一気に取得することができる。
【００８４】
以上説明したように、実施の形態５は、測定断面Ｙ＝０の場合についても、マルチプレーン法を適用することにより、測定断面Ｙ＝０の場合についても、測定時間の短縮を可能としている。
【００８５】
以上、本発明の２次元磁気共鳴イメージング方法の実施の形態の一例を示したが、測定断面およびパルスシークエンスは、これらの実施の形態に限定されるものではない。
【００８６】
すなわち、測定断面は、Ｚ＝０、Ｙ＝０に限定されず、平面状である限り、あらゆる断面を測定することができ、例えば、Ｚ＝１やＸ＋Ｙ＋Ｚ＝３平面などについても、本発明による２次元イメージングが可能である。
【００８７】
また、ＲＦ周波数と勾配磁場の組合せにより、あらゆる平面を選択励起することが可能であるため、パルスシークエンスについて、様々なバリエーションが考えられる。例えば、９０°パルスにＺＹ勾配磁場をかける時、Ｙ軸方向の勾配磁場強度をＺ軸方向の勾配磁場強度の半分にすれば、Ｚ＋(Ｙ／２)＝ａn平面を選択励起できる。このとき、Ｙ軸方向の勾配磁場強度は半分になり、Ｙ軸方向への選択幅（励起面のＹ軸方向の厚み）は２倍になる。１８０°パルスについても同様の操作が可能である。
【実施例】
【００８８】
以下、実施例に基づき、本発明について、さらに詳細に説明する。なお、本発明は下記実施例に限定されるものではない。
【００８９】
＜実施例１＞
実施の形態２で説明したパルスシーケンスを用いて、１次元測定法により１次元磁気共鳴情報を取得し、２次元磁気共鳴イメージングを行った。測定装置、測定試料、測定条件は下記の通りである。実施例１において測定に要した時間は２２秒であった。
【００９０】
(a)測定装置日本電子株式会社製  JNM-ECA500WB
イメージング用プローブ  DOTY製  20ｍｍ用プローブ
(b)測定試料Phantom CuSO4 水溶液
(c)測定条件
FOV=20*20
測定ポイント数：256
スライス厚：0.5ｍｍ
RG （レシーバゲイン）＝36
TE（エコータイム）=0.2msec
TR（繰り返し時間）=5sec
90°パルス＝49μsec
積算回数：2回
【００９１】
＜実施例２＞
実施の形態４で説明したパルスシーケンスを用いて、マルチプレーン法を適用した１次元測定法により１次元磁気共鳴情報を取得し、２次元磁気共鳴イメージングを行った。測定条件は下記の通りである。測定装置、測定試料は、実施例１と同じものを使用した。実施例２において測定に要した時間は２０秒であった。測定結果を図１９−１に示す。図１９−１の上図は、得られた断面画像であり、下図は、前記断面画像におけるＡ−Ａ´ライン上のシグナル強度である。
【００９２】
FOV=20*20
測定ポイント数：256
スライス厚：0.5ｍｍ
RG（レシーバゲイン）＝36
TE（エコータイム）=0.2msec
TR（繰り返し時間）=約5sec
90°パルス＝49μsec 
積算回数：2回
【００９３】
＜比較例１＞
２次元ＳＥ法により、２次元磁気共鳴イメージングを行った。測定条件は下記の通りである。測定装置、測定試料は、実施例１と同じものを使用した。実施例２において測定に要した時間は４２分４５秒であった。測定結果を図１９−２に示す。図１９−２の上図は、得られた断面画像であり、下図は、前記断面画像から、画像解析ソフト（Scion Image By Scion Corporation URL: http://www.scioncorp.com/index.htm）を用いて解析したＢ−Ｂ´ライン上のシグナル強度である。
【００９４】
FOV=20*20
測定ポイント数：256*256（2次元）
スライス厚：0.5ｍｍ
RG（レシーバゲイン）＝36
TE（エコータイム）=0.2msec
TR（繰り返し時間）=5sec
90°パルス＝49μsec
積算回数：2回
【００９５】
実施例１および２により、２次元ＳＥ法を用いた比較例１と比較して、１００分の１以下の測定時間で撮像することができた。また、実施例１および２において、１次元測定法を用いて得られたシグナルは、比較例１の２次元画像から抽出されたシグナルとまったく遜色のない結果となった。図１９−１と図１９−２を比較すると、１次元測定法と２次元ＳＥ法では、シグナル強度の差がわずかに見られる。これは、比較例１は、Ｘ軸およびＹ軸ともにフーリエ変換を行った後に解析されるのに対し、実施例２ではＸ軸のみにフーリエ変換をして得られるデータを即シグナル解析に利用できることに起因すると考えられる。すなわち、比較例１では、フーリエ変換を１軸多く行っているためにデータに差が出ていると考えられる。したがって、比較例１のように、２次元ＳＥ法によって画像を得てからのシグナル強度を解析するよりも、本発明のように、１次元測定法によって直接シグナル強度を解析する方が、正確なデータを取得できると考えられる。
【産業上の利用可能性】
【００９６】
以上のように、本発明にかかる核磁気イメージング方法は、空間分解能を低下させないで画像化に必要な磁気共鳴情報を短時間で取得することができるため、医療現場や研究等、測定スピードが要求される分野で幅広く利用することができる。
【図面の簡単な説明】
【００９７】
【図１】本発明の２次元磁気共鳴イメージング方法の概要を説明するためのフローチャートである。
【図２】１次元スピンエコー法のパルスシーケンスである。
【図３】１次元スピンエコー法による励起面および信号観測領域を示す模式図である。
【図４】１次元磁気共鳴情報を取得するラインの一例を示す模式図である。
【図５−１】従来の１次元測定法によるラインＸ1の励起時の状態を示す模式図である。
【図５−２】従来の１次元測定法によるラインＸ2の励起時の状態を示す模式図である。
【図６−１】本発明の１次元測定法によるラインＸ1の励起時の状態を示す模式図である。
【図６−２】本発明の１次元測定法によるラインＸ2の励起時の状態を示す模式図である。
【図６−３】本発明の１次元測定法によるラインＸnの励起時の状態を示す模式図である。
【図７】マルチプレーン法のパルスシーケンスの概念図である。
【図８】実施の形態１で励起される平面の模式図である。
【図９】ＭＲＩ装置を制御する制御装置３０の構成図である。
【図１０】実施の形態１のパルスシーケンスを実行するためのフローチャートである。
【図１１】実施の形態２で励起される平面の模式図である。
【図１２】実施の形態２において１次元磁気共鳴情報の測定に用いるパルスシーケンスである。
【図１３】実施の形態２のパルスシーケンスを実行するためのフローチャートである。
【図１４】実施の形態３で励起される平面の模式図である。
【図１５】実施の形態３において１次元磁気共鳴情報の測定に用いるパルスシーケンスである。
【図１６】実施の形態３のパルスシーケンスを実行するためのフローチャートである。
【図１７】実施の形態４のパルスシーケンスを実行するためのフローチャートである。
【図１８】実施の形態５のパルスシーケンスを実行するためのフローチャートである。
【図１９−１】実施例２によって得られた断面画像および当該断面画像におけるＡ−Ａ´ライン上のシグナル強度を示す図表である。
【図１９−２】比較例１によって得られた断面画像および当該断面画像におけるＢ−Ｂ´ライン上のシグナル強度を示す図表である。
【図２０】２次元ＳＥ法のパルスシーケンスである。
【図２１】２次元ＳＥ法の概念図である。
【図１】
2006026311.tif 000003
【図２】
2006026311.tif 000004
【図３】
2006026311.tif 000005
【図４】
2006026311.tif 000006
【図５−１】
2006026311.tif 000007
【図５−２】
2006026311.tif 000008
【図６−１】
2006026311.tif 000009
【図６−２】
2006026311.tif 000010
【図６−３】
2006026311.tif 000011
【図７】
2006026311.tif 000012
【図８】
2006026311.tif 000013
【図９】
2006026311.tif 000014
【図１０】
2006026311.tif 000015
【図１１】
2006026311.tif 000016
【図１２】
2006026311.tif 000017
【図１３】
2006026311.tif 000018
【図１４】
2006026311.tif 000019
【図１５】
2006026311.tif 000020
【図１６】
2006026311.tif 000021
【図１７】
2006026311.tif 000022
【図１８】
2006026311.tif 000023
【図１９−１】
2006026311.tif 000024
【図１９−２】
2006026311.tif 000025
【図２０】
2006026311.tif 000026
【図２１】
2006026311.tif 000027