磁気共鳴画像法のソースを表示

<div align="right">  
<font size="+1">[http://researchmap.jp/kojifujimoto 藤本晃司]</font><br>
''京都大学医学部医学研究科脳機能総合研究センター''<br>
<font size="+1">[http://researchmap.jp/takashihanakawa 花川 隆]</font><br>
''京都大学医学部医学研究科脳統合イメージング分野''<br>
DOI：<selfdoi />　原稿受付日：2020年3月28日　原稿完成日：2020年X月XX日<br>
担当編集委員：[http://researchmap.jp/read0048432 定藤 規弘]（自然科学研究機構生理学研究所 大脳皮質機能研究系）<br>
</div>

英語名：magnetic resonance imaging　独：Magnetresonanztomographie　仏：imagerie par résonance magnétique

英略号：MRI

{{box|text=　磁気共鳴画像法(MRI)は、生体中に多数存在する水素原子核（プロトン）と、外部から与える電磁波との相互作用を利用することで、多彩な生命現象を可視化する技術であり、特にヒトを対象とする非侵襲的脳科学計測研究においては最も重要な手法の一つである。MRIによる生体内情報の可視化には、強い静磁場を形成するドーナツ型の超電導磁石、生体内のプロトンにエネルギーを送信して核磁気共鳴を生じさせ、緩和で放出される電磁波を受信するコイル、プロトンの空間分布をエンコード・デコードするための勾配磁場コイルが必要となる。電磁波で励起されたプロトンが、エネルギーを放出して定常状態に戻る際に放出する電磁波は、プロトン周囲の微小環境を反映したT<sub>1</sub>, T<sub>2</sub>, T<sub>2</sub>T<sub>2</sub><sup>*</sup>と呼ばれる時定数を持つ。RFコイルによる電磁波と勾配磁場の印加方法を巧みに操ることでこれらの時定数を強調した画像が得られる。さらに、複数の受信コイルから得られた信号を組み合わせる手法（パラレルイメージング）と併用することで、機能的磁気共鳴画像（fMRI）実験において1秒以内に全脳を撮像することも可能となっている。}}

== はじめに ==
　地球上の生命体は、[[wj:水|水]]分子、[[脂質]]や[[アミノ酸]]など、[[wj:水素|水素]][[wj:原子|原子]]を含んだ数多くの化合物から構成されている。医療や脳科学研究に広く用いられているMRIは、これら化合物中の水素原子核（[[wj:プロトン|プロトン]]：物理学による定義）が有する小さな[[wj:磁石|磁石]]としての性質（[[wj:原子核スピン|原子核スピン]]）と、これに外部から特定の[[wj:電磁波|電磁波]]を与えた際に生じる相互作用（[[wj:核磁気共鳴現象|核磁気共鳴現象]]）を利用して、脳を含む生体内の情報を非侵襲的に画像化する手法である。

　1990年代初頭の[[機能的磁気共鳴画法]]（fMRI）の原理の発見により、[[ヒト]]を対象とする非侵襲的脳科学計測研究において最も重要な手法の一つになり、その後も[[拡散強調トラクトグラフィー]]による[[白質]]線維連絡の評価法など脳科学のツールとして発展を続けている。

　アメリカ合衆国の化学者[[wj:ポール・ラウターバー|Paul Lauterbur]]とイギリスの物理学者[[wj:ピーター・マンスフィールド|Peter Mansfield]]は、MRIに関する発見により、2003年の[[wj:ノーベル生理学・医学賞|ノーベル生理学・医学賞]]を受賞している。


== MRI信号の原理：核磁気共鳴と緩和 ==
[[File:Hanakawa Fig 1.png|thumb|right|300px| '''図1. 京都大学医学部医学研究科のヒト用7-T MRI装置''']]
=== 超電導磁石による静磁場の形成 ===
　現在の医療・脳科学研究では、1.5T（テスラ）から7Tの静磁場強度を持つMRIが使用されている('''図1''')。このような磁場強度を持つMRIは、NbTi（ニオブチタン）合金による[[wj:超伝導|超伝導]]ワイヤーを何重にも巻いたドーナツ型の空芯コイルを、液体[[wj:ヘリウム|ヘリウム]]（沸点4.2K）で常時冷却可能な[[wj:魔法瓶|デュワー瓶]]（Dewar flask）中に置き、そこに電流を流すことで、電気抵抗による発熱の問題を克服した強力な[[wj:超電導磁石|超電導磁石]]を用いている。

=== 静磁場中の原子核スピンの振る舞い ===
　この超電導磁石により形成される[[wj:静磁場|静磁場]]（外部磁場またはB0とも呼ばれる）中に置かれた[[wj:原子核|原子核]]（典型的には生体内の水素原子核）は、固有の周波数（[[wj:ラーモア歳差運動|ラーモア周波数]]ω）で静磁場の方向を回転軸とする[[wj:歳差運動|歳差運動]]を行う。ラーモア周波数は静磁場の強さに比例する（<math>\omega=\gamma B_0</math>、<math>B_0</math>は静磁場の強さ、<math>\gamma</math>は磁気回転比と呼ばれる定数）。強い静磁場内では、静磁場の向きと一致（＋）した方向と逆（－）方向を向いて歳差運動を行う原子核の個数はわずかに異なることが知られている（[[wj:ゼーマン効果|ゼーマン分裂]]）。生体内の水素原子核の数が非常に多いため、強い静磁場内に置かれた生体内には＋方向を向いた巨視的磁化が形成される。
[[File:Hanakawa Fig 2.gif|thumb|right|300px|'''図2. RFパルスによる励起、および<math>T_1</math>、<math>T_2</math>緩和＞定常状態（緩和しきった状態）のスピン集団に対して、共鳴周波数でRFパルスを照射した後のスピンの緩和過程'''<br>共鳴周波数と同じ速度で回転する座標系から眺めているため、スピン集団は静止してみえる。上段の３つの図はスピン集団を各々上(xy平面)、横（xz平面およびyz平面）から眺めた図。下段左は斜め上方から眺めた図。下段右の上段には縦磁化の時間変化を、下段右の下段には横磁化の時間変化を示す。RFパルスを照射されたスピンはまず、xy平面上に倒れる（＝横磁化の出現）。xy平面上で円弧を描くように広がる（<math>T_2^*</math>緩和）その後、z軸の（＋）方向にむかって、z軸方向の磁化（＝縦磁化）が回復して（定常状態に戻って）ゆく。]]
=== 外部からの電磁波による「核磁気共鳴現象または励起現象」 ===
　例えば3T MRI装置における水素原子核のラーモア周波数は128MHzである。この周波数は[[wj:FMラジオ|FMラジオ]]が使用する周波数帯（radio frequency，RF）である。送信コイルを用いてラーモア周波数の回転磁場（RFパルスまたは<math>B_1</math>とも呼ばれる）を照射すると（通常は数ミリ秒程度のごく短時間）、水素原子核がエネルギーを吸収し、低い[[wj:エネルギー準位|エネルギー準位]]から高いエネルギー準位に遷移する（核磁気共鳴）。この際、外部から観測される磁化（巨視的磁化）は、[[wj:回転座標系|回転座標系]]において('''図2''')回転磁場および静磁場の双方に直交する方向を軸として回転する。この巨視的磁化は、静止座標系においては、静磁場と直交する平面上で、共鳴周波数で回転する磁化（横磁化）の出現および静磁場と平行な成分（縦磁化）の減少として観測される（励起）。

=== 電磁波（RF）の停止に続く緩和現象 ===
　励起された原子核は、RFパルスの照射が終わった直後から、エネルギーを放出して定常状態に戻ろうとする（緩和現象）。この緩和現象の主なものとして、縦磁化の回復（縦緩和、<math>T_1</math>緩和）と横磁化の消失（横緩和、<math>T_2</math>緩和）がある。

==== 縦緩和（''T<sub>1</sub>''緩和） ====
　励起後の巨視的磁化の縦磁化成分の回復（<math>T_1</math>緩和）は[[wj:スピン―格子緩和|スピン―格子緩和]]とも呼ばれ、原子核スピンが周囲の格子(プロトン以外の物質)にエネルギーを放出して[[wj:熱平衡状態|熱平衡状態]]に戻ることによる現象である。この際の縦磁化の時間変化は、巨視的磁化ベクトル<math>M</math>の、ある時間<math>t</math>における静磁場方向に平行な成分を<math>M_z</math>、熱平衡状態の縦磁化を<math>M_{z,eq}</math>とおくと

:::<math>
M_z (t)=M_{z,eq}-[M_{z,eq}-M_z (0)]e^{-t/T_1}
</math>

で回復する。<math>T_1</math>はこの回復の早さを決める時定数である。

==== 横緩和（''T<sub>2</sub>''緩和） ====
　励起後の巨視的磁化の横磁化成分の消失（<math>T_2</math>緩和）は[[wj:スピン―スピン緩和|スピン―スピン緩和]]とも呼ばれ、原子核スピン集団のコヒーレンス（位相同期）の消失（即ち位相分散）による。個々の原子核スピンが「感じる」静磁場はわずかながらランダムに変動しており、結果として時間の経過とともにスピン集団のコヒーレンスが失われてゆく（ある原子核スピンは位相が早くなり、別の原子核スピンは位相が遅くなる）。この際の横磁化の時間変化は、巨視的磁化ベクトル<math>M</math>の、ある時間<math>t</math>における静磁場方向に直交する成分を<math>M_{xy}(t)</math>とおくと<math>M_{xy}(t)=M_{xy}(0)e^{-t/T_2}</math>で減衰する。<math>T_2</math>はこの減衰の早さを決める時定数である。

==== ''T<sub>2</sub><sup>*</sup>''緩和 ====
　前述の''T<sub>2</sub>''緩和に加え、現実のMRI環境では、励起されたスピン集団の外部磁場は場所により少しずつ異なる（ただしここでは緩和過程中の時間変動はないものと仮定する）。結果として外部から観測される横磁化成分の減衰の早さは、静磁場の空間的分布の不均一性の強さ<math>\Delta B_0</math>および磁気回転比（gyromagnetic ratio <math>\gamma</math>）を用いて

:::<math>
\frac{1}{T_2^*}=\frac{1}{T_2}+\gamma\Delta B_0=\frac{1}{T_2}+\frac{1}{T_{B0inhomogeneity}}
</math>

とあらわされ、''T<sub>2</sub><sup>*</sup>''と呼ばれる。

== 空間情報のエンコーディングとデコーディング ==
　　MRIは、緩和過程において放出される電磁波を受信コイルにより測定し、励起された水素原子核集団の挙動の違いを画像化している。まず二次元（2D）MRIにおいては、スライス方向に直行する軸に勾配磁場（傾斜磁場とも呼ばれる）を短時間加え、その間に帯域幅の限られたRFパルスを与えることで「スライス選択励起」を行っている。三次元（3D）MRIでは厚さを持った範囲（スラブ）を励起する勾配磁場とRFパルスを与える。

　励起されたスライスやスラブ内での水素原子核集団の分布を知るために、勾配磁場を用いたエンコーディングおよび[[wj:離散フーリエ変換|離散フーリエ変換]]を用いたデコーディングが必要である。典型的には、RFパルスによる励起に加えて、撮像領域（field-of-view, FOV）の中心部からx方向あるいはy方向に線形に変化する勾配磁場を追加することで、FOV内の水素原子核の共鳴周波数が位置に依存する状態をつくりだすことが可能となる。そうすれば、FOV内の特定の部位に存在する水素原子核集団は、特定の周波数の電磁波を放出することになるため、位置情報の特定が可能となる。MRIでは勾配磁場の与え方を工夫することで、測定された電磁波から空間情報を読み取る際に離散フーリエ変換を用いる。

　MRIの測定信号（電磁波）の画像化の理解には、離散フーリエ変換（discrete Fourier transform, DFT）の原理の理解が重要であるため簡単に説明する。フランスの数学者[[wj:ジョゼフ・フーリエ|ジョゼフ・フーリエ]]は、あらゆる周期関数（や周期信号）は、三角級数の（無限の）和として表現できることを発見した。すなわち、実数xを変数とする周期2nの[[wj:周期関数|周期関数]]<math>f(x)</math>について、

:::<math>a_n=\frac{1}{\pi}\int_{-\pi}^{\pi} f(t)</math>cos<math> nt\ dt,(n=0,1,2,3,\ldots)</math>

:::<math>b_n=\frac{1}{\pi}\int_{-\pi}^{\pi} f(t)</math>sin<math> nt\ dt,(n=0,1,2,3,\ldots)</math>

と置くと、

:::<math>f(x)=\frac{a_0}{2}+\sum_{n=1}^{\infty}(a_n</math>cos<math>nx+b_n</math>sin<math>nx)</math>

と書ける。これを[[wj:複素数|複素数]]に拡張すれば、

:::<math>c_n=\frac{1}{2\pi}\int_{-\pi}^\pi f(t) e^{-\iota nt} dt,\ (n=0,\pm1,\pm2,\ldots)</math>

として

:::<math>f(x)=\sum_{n=-\infty}^{\infty}(c_n e^{-ixn})</math>

と書ける。

　f(x)がデジタル信号で離散化できる場合、上記のフーリエ級数表現を2次元、かつ複素数に拡張したものは以下のようにあらわすことができる（離散フーリエ変換, discrete Fourier transform, DFT）

:::<math>F(u,v)=\sum_{x=0}^{M-1}\sum_{y=0}^{M-1}f(x,y)e^{-\iota 2\pi(\frac {u}{M}x+\frac{v}{N} y)}</math>

　また、この逆変換にあたる逆離散フーリエ変換（inverse discrete Fourier transform, IDFT）は以下のようにかける

:::<math>f(x,y)=\frac{1}{MN}\sum_{u=0}^{M-1}\sum_{v=0}^{N-1}F(u,v)e^{-\iota 2\pi(\frac {u}{M}x+\frac{v}{N} y)}</math>


　これは画素数MxNの画像の画素値が、周波数領域(u,v)空間では波数0～M-1,0～N-1個の[[wj:三角関数|三角関数]]の係数値として表現できることを示している。
例えばu=1,v=0とした場合、

:::<math>F(1,0)=\sum_{x=0}^{M-1}\sum_{y=0}^{M-1}f(x,y)e^{-\iota 2\pi(\frac {1}{M}x)}</math>

となるが、これはMRIでは撮像領域(FOV)の右端から左端にかけて複素数の重み係数

:::<math>e^{i\theta}=</math> cos <math>\theta + i</math>⁡sin<math>\theta, \theta =[-\pi\ \pi]</math>

を画像に乗じたのちに総和をとったことに等しい。同様に、<math>u=2,\ v=0</math>とした場合には、FOVの右端から左端にかけて<math>\theta =[-2\pi\ 2\pi]</math>の重み係数を乗じたのちに総和をとったことに等しい。

　勾配磁場コイルを用いて空間内に線形の周波数変化をもたらせば、場所に応じた連続的な位相の変化としてこの重み係数を物理的につくりだすことが可能であり、結果として撮像対象にフーリエ変換を行っていることと等しい。

　画像<math>f(x,y)</math>の周波数空間での表現である<math>F(u,v)</math>はMRIにおいてはk-spaceと呼ばれる。MRIではk-spaceの信号を逆フーリエ変換することで画像を得ている。

== 主なMRI撮像法 ==
　さまざまなMRI撮像法が提案されているが主な違いはRFパルスを照射する回数、タイミングや大きさ、勾配磁場の印加法である。これら電磁波の照射の時系列制御がキーであるため撮像シークエンスとも呼ぶ。代表的なMRI撮像法を簡単に紹介する。
=== スピンエコー法 ===
Spin echo

　励起のための電磁波（90°RFパルス）を与えたのちに、再度180°RFパルス（再収束パルス、refocus pulse）を与えることで、局所の静磁場の不均一性による位相分散の影響を取り除くことが出来る。

=== グラジエントエコー法 ===
Gradient echo

　励起のための電磁波（RFパルス）を与えたのちに、再収束パルスを与えず、勾配磁場を用いて信号を取り出す手法。局所の静磁場の不均一性による位相分散の影響を取り除くことは出来ないが、再収束パルスを必要としないため高速撮像に向く。また、後述のBOLD fMRIのように、（鉄などによる）局所の静磁場の不均一性を強調したい場合にも用いられる。

=== エコープラナー法 ===
Echo planar imaging, EPI

　一度のRFパルスの後、[[グラジエントエコー法]]あるいは[[スピンエコー法]]の信号収集時間を極端に延長し、読み出し勾配磁場を急速に変動させることで連続的なグラジエントエコーを発生させ、画像化に必要なデータを全て収集してしまう方法。<math>T_2^*</math>緩和の影響が強く、かつ原理的にもっとも高速な撮像法の一つである。fMRIで利用されるblood oxygenation-level dependent （BOLD）コントラスト（後述）は<math>T_2^*</math>緩和に依存し、かつ高い時間分解能が必要とされるため、本手法が用いられる('''図3''')。

[[ファイル:Hanakawa Fig 3.png|サムネイル|右|400px| '''図3. それぞれの強調像の比較'''<br>三次元撮像（3d）T<sub>1</sub>強調画像（magnetization-prepared rapid gradient echo, MPRAGE）では脳脊髄液が黒く、白質が明るい灰色、灰白質がその中間の暗い灰色に見える。<br>T<sub>2</sub>強調画像（sampling perfection with application optimized contrasts using different flip angle evolution, SPACE）では逆に白質が最も暗く、灰白質、脳脊髄液と明るく見える。<br>機能的MRI用のmulti-band echo planar imaging（MB-EPI）ではT<sub>2</sub>に似たコントラストを示すが、副鼻腔による磁場不均一の影響を受ける<math>T_2^*</math>強調であるため前頭部の信号の一部が欠損している。]]

== 主なMRIコントラスト ==
=== T<sub>1</sub>強調像 ===
T<sub>1</sub>-weighted image, T<sub>1</sub>WI

　縦磁化が十分に回復しないうちに信号収集を行うことで、各組織における縦磁化回復の早さの違いを強調した画像が得られる。具体的にはスピンエコー法では励起から次の励起までの時間（repetition time, TR)および励起から収集までの時間(echo time, TE)を短くすることで、またグラジエントエコー法では[[wj:フリップ角|フリップ角]]（FA）を適度に大きく、TEを短くすることで[[T1強調像|T<sub>1</sub>強調像]]が得られる。

　撮像部位に流入する血管内の血液が高信号を示すことを利用して、[[造影剤]]を用いずに[[脳血管]]を可視化する手法であるTime-of-flight (TOF)法では、グラジエントエコー法によるT1強調像が用いられる。

=== T<sub>2</sub>強調像 ===
T<sub>2</sub>-weighted image, T<sub>2</sub>WI

　励起から収集までの時間（TE）を長くとれば、各組織における横磁化の減衰速度の違いを強調した画像が得られる('''図3''')。画像収集にスピンエコー法を用いた場合、得られる画像は時定数[[T2強調像|T<sub>2</sub>強調像]]となる。

=== T<sub>2</sub><sup>*</sup>強調像 ===
T<sub>2</sub>-star-weighted image, T<sub>2</sub><sup>*</sup>WI

　TEを長くとり、画像収集にグラジエントエコー法を用いた場合、得られる画像は[[T2*強調像|T<sub>2</sub><sup>*</sup>強調像]]となる('''図3''')。得られた画像に特殊な画像処理を施すことで、[[磁化率強調像]]（susceptibility-weighted image, SWI）<ref name=Haacke2004><pubmed>15334582</pubmed></ref> や定量的磁化率マップ（quantitative susceptibility map, QSM）<ref name=Li2011><pubmed>21224002</pubmed></ref> といった画像が得られる。

=== 拡散強調像 ===
Diffusion-weighted image, DWI

　撮像の際に、数msから数十ms程度のごく短時間で反転する一組の強い勾配磁場、すなわち[[wj:運動検出勾配磁場|運動検出勾配磁場]]（motion probing gradient, MPG）を追加することで水の動きを強調した画像。MPG を与えた方向に拡散する水由来の信号は低下する。[[脳梗塞]]を起こした部位では健常組織よりも水の拡散が制限されることが知られており、脳梗塞を早期に検出する手法として臨床で広く用いられている。6方向以上のMPGを用いて拡散強調像を収集し、テンソルモデルを用いて撮像単位内に存在する水の動きやすさの「方向」を推定することも可能である（[[拡散テンソルイメージング]], diffusion tensor imaging, DTI）。この推定値を用いて白質を通る神経線維の走行を推測する方法を[[拡散テンソルトラクトグラフィー]]（diffusion tensor tractography）と呼ぶ<ref name=Basser2000><pubmed>11025519</pubmed></ref><ref name=Melhem2002><pubmed>11756078</pubmed></ref> 。拡散テンソルイメージングは各ボクセル内に1種類の神経線維の存在を仮定しているため、神経線維が交叉する撮像単位での推定に限界があり、近年ではより高度なモデルを用いた解析手法も提唱されている（Diffusion spectrum imaging, DSIなど）<ref name=Wedeen2008><pubmed>18495497</pubmed></ref> 。

=== MRスペクトロスコピー ===
MR spectroscopy、MRS

　MRSは、化合物に存在するプロトンの共鳴周波数が、水の水素原子の共鳴周波数とわずかに異なることを利用して、測定対象物の化合物を推定、あるいはその濃度を計測する手法である。脳MRSにおいては、[[シナプス]]間の[[神経伝達物質]]として最もポピュラーな[[グルタミン酸]]（glutamate）をはじめ、複数の神経伝達物質が計測のターゲットとなる。

== 撮像高速化の手法 ==
=== パラレルイメージング ===
Parallel imaging

　従来、撮像対象に対して一つの受信コイルを用いた撮像が行われていたが、1990年代に感度領域が限られた複数の受信コイルを並べて配置する[[wj:アレイコイル|アレイコイル]]と呼ばれる技術が開発された<ref name=Roemer1990><pubmed>2266841</pubmed></ref> 。コイルごとに少しずつ異なる感度分布をうまく利用することで、MRIの信号取得をスキップし、撮像時間を短縮する方法が開発された。この技術は複数の受信コイルを用いて同時（並列）に得られた信号を用いていることから、パラレルイメージングと呼ばれる。コイルの空間感度分布を用いる手法をsensitivity encoding (SENSE<ref name=Pruessmann1999><pubmed>10542355</pubmed></ref> 、コイルの空間感度分布を用いず、k-spaceにおける近傍点の相対的な関係をもとに重み係数を決定し欠損データを補完する手法をGeneRalized Autocalibrating Partial Parallel Acquisition（GRAPPA）と呼ぶ<ref name=Griswold2002><pubmed>12111967</pubmed></ref> 

=== Multi-band/Simultaneous Multi-Slice撮像 ===
　2010年にDavid Feinbergらは「スライス選択励起」RFパルスを工夫することで、複数のスライスをひとつのRFパルスで同時に励起する手法を開発し<ref name=Feinberg2010><pubmed>21187930</pubmed></ref> 、multiband MRI と名付けた('''図3''')。従来はEPIを用いて全脳を撮像するのに2-3秒程度必要であったが、multiband撮像を用いると1秒以内に全脳を撮像することが可能である。

== 機能的磁気共鳴画像法 ==
　磁気共鳴機能画像法（functional MRI, fMRI）ともいい、生体の脳をMRIで数分～数十分間連続撮像する間に、脳活動（神経活動とシナプス活動の総和）に相関して変化するMRI信号変化を非侵襲的に計測する手法である。詳細については[[機能的磁気共鳴画像法]]の項を参照されたい。

==関連項目==
* [[機能的磁気共鳴画像法]]
==参考文献==
<references />