(57)【要約】【課題】曲面の断層画像の画質を向上させると共に所望の曲面をインタラクティブに設定できリアルタイムに近い時間で断層画像を容易に得る。【解決手段】磁気共鳴イメージング装置は、入力されたスライス曲線を受け付けて、被検体における柱状領域をスライス曲線に沿ってオーバーラップするように移動させるときのそれぞれの位置を算出すると共に、算出したそれぞれの位置情報およびスライス曲線の位置情報を保持する位置計算手段と、柱状領域を励起して柱状領域から発生するＭＲ信号を検出する撮像手段と、撮像手段による一連の動作を撮像手段に繰り返し行わせる制御手段と、検出された複数のＭＲ信号を用いて画像を再構成する画像再構成処理手段と、再構成された画像である曲面断層画像を表示する表示手段と、を備え、画面上でインタラクティブにスライス曲線を延長する入力を受け付けた場合、延長されたスライス曲線に応じて表示中の曲面断層画像を更新する。【選択図】図１０
2018139779.tif 000002
【特許請求の範囲】
【請求項１】
被検体の所定のスライス曲面に直交する断面画像として表示されたスカウト画像に対して重畳するように前記スライス曲面を指定するためのスライス曲線を入力する入力手段と、
入力されたスライス曲線を受け付けて、前記被検体における所定の柱状領域を、前記スライス曲線に沿ってオーバーラップするように移動させるときのそれぞれの位置を算出すると共に、算出したそれぞれの位置情報および前記スライス曲線の位置情報を保持する位置計算手段と、
前記位置計算手段で算出されたそれぞれの位置に配置される前記柱状領域を、傾斜磁場の印加および高周波磁場の照射により励起して前記柱状領域から発生する磁気共鳴信号を検出する撮像手段と、
前記柱状領域の位置を所定の順序で変化させながら前記撮像手段による一連の動作を前記撮像手段に繰り返し行わせる制御手段と、
検出された複数の前記磁気共鳴信号を用いて画像を再構成する画像再構成処理手段と、
再構成された画像である曲面断層画像を表示する表示手段と、を備え、
撮像中および画像の再構成中であっても前記スライス曲線を延長する入力を受け付け可能であり、前記スライス曲線を延長する入力を受け付けた場合、延長されたスライス曲線に応じて表示中の曲面断層画像を更新することを特徴とする磁気共鳴イメージング装置。
【請求項２】
前記柱状領域の断面の形状は菱形であり、前記菱形の一方の対角線は他方の対角線よりも短く、前記一方の対角線の向きに前記柱状領域を移動させる請求項１に記載の磁気共鳴イメージング装置。
【請求項３】
前記柱状領域を形成する前記菱形の他方の対角線の長さは、スライス厚の１．５〜２．５倍である請求項２に記載の磁気共鳴イメージング装置。
【請求項４】
前記柱状領域を形成する前記菱形の一方の対角線の長さは、画素サイズの１．５〜２．５倍である請求項２または請求項３に記載の磁気共鳴イメージング装置。
【請求項５】
前記制御手段は、前記柱状領域の位置を変化させることで前記スライス曲線に沿って形成される前記柱状領域を１以上前方または後方に飛び越えて、前記撮像手段による一連の動作を繰り返し行わせながら、前記スライス曲線で指定されるスライス曲面内のすべての前記柱状領域を励起させる請求項１から請求項４のいずれか一項に記載の磁気共鳴イメージング装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、磁気共鳴イメージング装置に関する。
【背景技術】
【０００２】
磁気共鳴イメージング装置において、一般的には、撮像して得られる画像は平面に沿った画像（以下、平面の画像）である。しかし、生体での組織の形状は、湾曲等している場合が多い。例えば、脊椎に沿ったコロナル断面（冠状断面）を得たい場合、脊髄はＳ字状に湾曲しているため、１枚の平面の画像では、脊髄に沿った画像は得られない。
【０００３】
そのため、次の２つのいずれかの手法によって脊髄の画像を読影しているのが現状である。第１の手法では、平面状のスライス断面についての複数の平面の画像を取得し、湾曲した脊髄をすべて含むような画像が得られるような撮像を行う方法である。この場合、取得した複数の平面の画像に跨って部分的に表示されている脊髄の画像で読影することになる。
第２の手法では、前記の撮像によって得られた複数の平面の画像を用いて、ＭＰＲ（multi planar reconstruction）処理を用いて、湾曲した脊髄における曲面に沿った画像（以下、曲面の画像）を計算処理により生成し、その曲面の画像で読影することとなる。
【０００４】
ここで、第１の手法の場合、湾曲した脊髄は複数の平面の画像に跨っていて、各１枚の画像には脊髄の一部だけしか描出されていない。よって、１枚の画像では脊髄の一部だけしか描出されていない画像を複数枚寄せ集めて、それらを用いて脊髄の全体像を思い浮かべ、読影しなければならない。そのため、読影が困難であるという問題があった。
【０００５】
また、第２の手法の場合、次のような一連の作業を行わなければならない。すなわち、操作者は、まず、前記の撮像によって得られた複数の平面の画像においてＭＰＲ処理に用いる画像を選択し、選択した画像を磁気共鳴イメージング装置に読み込ませ、生成しようとする曲面の位置設定を行う。その上で、この磁気共鳴イメージング装置が、操作者により設定された位置でＭＰＲ画像を計算して生成する処理を実施することで曲面の画像を得る。そのため、撮像によって複数の画像を得た後に、追加の手間と時間を要してしまう。さらに、ＭＰＲ処理で作成した画像の空間分解能は、特に撮像によって得られた画像のスライス厚に強く依存している。そのため、具体的には、スライス厚が厚いと、ＭＰＲ処理で作成した画像の空間分解能は低下し診断に耐えない画像となるという問題があり、また、スライス厚が薄いと、必要な画像の枚数が増えるため撮像時間が延長するという問題があった。
【０００６】
そこで、平面に沿って撮像するのではなく、曲面に沿って撮像を行い、画像を再構成することが考えられる。例えば、特許文献１に記載の磁気共鳴イメージング装置は、被検体の所定の柱状領域を励起し、柱状領域からの磁気共鳴信号を検出する動作を、柱状領域を移動させながら繰り返し行い、得られた複数の磁気共鳴信号を用いて画像を再構成する。この従来技術では、ディスプレイに表示されたスカウト画像であるサジタル画像上に、操作者がマウス等を操作することにより曲線を描くことで、磁気共鳴イメージング装置が、この曲線を受け付けて、曲線を含み且つスカウト画像に垂直な曲面を撮像すべき曲面として設定している。
【先行技術文献】
【特許文献】
【０００７】
【特許文献１】特開２００９−１７８２８７号公報
【発明の概要】
【発明が解決しようとする課題】
【０００８】
例えば診断のために、脊髄等の曲面に沿った撮像によって得られる断層画像について画質の向上が望まれている。具体的には、曲面に沿った撮像において、１つの柱状領域は、再構成画像の１つの画素に相当し、断面のサイズである太さが再構成画像における空間分解能に対応している。また、同じ太さであれば、曲面の方向における画像空間分解能は、柱状領域の間隔に依存する。このため、例えば柱状領域の間隔が所定値である場合、柱状領域を太くし過ぎると、画像の空間分解能が低下するので精細な画像が得られなくなる問題がある。逆に、柱状領域を細くし過ぎると、柱状領域の体積が減ることにより、柱状領域から発生する磁気共鳴信号の強度が低下するので、得られる画像のＳＮ比が低下する問題がある。
【０００９】
また、操作者が指定して撮像した断面以外の位置の断層画像を得たいと望む場合がある。このような場合、従来技術では、同じことを再度、繰り返さなければならない。すなわち、操作者は、再度、曲線を指定し、磁気共鳴イメージング装置は、再度、曲面に沿った撮像を実施しなければならない。したがって、従来は、生体の湾曲した組織の形状に沿った所望の位置における曲面の断層画像を容易に得ることはできなかった。
【００１０】
そこで、曲面の断層画像の画質を向上させると共に所望の位置における曲面の断層画像を容易に得ることのできる磁気共鳴イメージング装置を提供することを課題とする。
【課題を解決するための手段】
【００１１】
前記課題を解決するため、本発明に係る磁気共鳴イメージング装置は、被検体の所定のスライス曲面に直交する断面画像として表示されたスカウト画像に対して重畳するように前記スライス曲面を指定するためのスライス曲線を入力する入力手段と、入力されたスライス曲線を受け付けて、前記被検体における所定の柱状領域を、前記スライス曲線に沿ってオーバーラップするように移動させるときのそれぞれの位置を算出すると共に、算出したそれぞれの位置情報および前記スライス曲線の位置情報を保持する位置計算手段と、前記位置計算手段で算出されたそれぞれの位置に配置される前記柱状領域を、傾斜磁場の印加および高周波磁場の照射により励起して前記柱状領域から発生する磁気共鳴信号を検出する撮像手段と、前記柱状領域の位置を所定の順序で変化させながら前記撮像手段による一連の動作を前記撮像手段に繰り返し行わせる制御手段と、検出された複数の前記磁気共鳴信号を用いて画像を再構成する画像再構成処理手段と、再構成された画像である曲面断層画像を表示する表示手段と、を備え、撮像中および画像の再構成中であっても前記スライス曲線を延長する入力を受け付け可能であり、前記スライス曲線を延長する入力を受け付けた場合、延長されたスライス曲線に応じて表示中の曲面断層画像を更新することを特徴とする。
【００１２】
かかる構成によれば、磁気共鳴イメージング装置は、被検体における柱状領域を、スライス曲線に沿ってオーバーラップするように移動させながら柱状領域から発生する磁気共鳴信号を検出する。このように柱状領域をオーバーラップさせることにより、柱状領域の体積を大幅に減らすことなく磁気共鳴信号を検出できるため、得られる画像のＳＮ比が向上する。また、オーバーラップさせることは、柱状領域の間隔を狭く設定したことになるので、結果として画像空間分解能も向上させることができる。
また、磁気共鳴イメージング装置は、スライス曲線を延長する入力を受け付けた場合、延長されたスライス曲線に応じて表示中の曲面断層画像を更新する。これにより、操作者は、最初に入力したスライス曲線で指定されたスライス曲面に沿った曲面断層画像の表示を参照し、画面上でインタラクティブに、最初に入力したスライス曲線を延伸することができる。したがって、曲面の撮像を継続的に実施しながら追加された位置における画像の再構成処理を行うので、撮像を一旦中止して再度位置決めして撮像することなく、所望の曲面をインタラクティブに設定できリアルタイムに近い時間で断層画像が容易に得られる。
【００１３】
また、本発明に係る磁気共鳴イメージング装置は、前記柱状領域の断面の形状は菱形であり、前記菱形の一方の対角線は他方の対角線よりも短く、前記一方の対角線の向きに前記柱状領域を移動させるようにしてもよい。これにより、画像空間分解能を向上させることができる。また、この場合、前記柱状領域を形成する前記菱形の他方の対角線の長さは、スライス厚の１．５〜２．５倍であることが好ましく、前記柱状領域を形成する前記菱形の一方の対角線の長さが、画素サイズの１．５〜２．５倍であることが好ましい。
【００１４】
また、本発明に係る磁気共鳴イメージング装置は、前記制御手段が、前記柱状領域の位置を変化させることで前記スライス曲線に沿って形成される前記柱状領域を１以上前方または後方に飛び越えて、前記撮像手段による一連の動作を繰り返し行わせながら、前記スライス曲線で指定されるスライス曲面内のすべての前記柱状領域を励起させるようにしてもよい。
【００１５】
かかる構成によれば、磁気共鳴イメージング装置は、スライス曲線に沿って算出された位置に配置される各柱状領域において、隣接する柱状領域を連続的に励起するのでなく、１以上の柱状領域を飛び越えるほど離間した位置に形成される柱状領域を順番に励起する。よって、一柱状領域に着目すると、前記一柱状領域を励起した後、ある程度の時間を待ち、前記一柱状領域における磁化の緩和が進んで元の状態に戻ってから、この一柱状領域に対して隣接する他柱状領域を励起することになる。これにより、一柱状領域と、これに隣接する柱状領域との干渉を抑制し、干渉に由来する信号強度の低下が抑制された磁気共鳴信号を得ることができるので画像のＳＮ比が向上できる。
【発明の効果】
【００１６】
本発明に係る磁気共鳴イメージング装置によれば、得られる画像のＳＮ比が向上すると共に画像空間分解能を向上させることができるので、曲面の断層画像の画質を向上させることができる。
また、本発明に係る磁気共鳴イメージング装置によれば、操作者は、画面上でインタラクティブに、最初に入力したスライス曲線を延伸することができるので、撮像を一旦中止して再度位置決めして撮像することなく、所望の位置における曲面の断層画像を容易に得ることができる。
【図面の簡単な説明】
【００１７】
【図１】本発明の実施形態に係る磁気共鳴イメージング装置の構成を模式的に示すブロック図である。
【図２】本発明の実施形態に係る磁気共鳴イメージング装置のパルスシーケンスの一例を示す図である。
【図３】本発明の実施形態に係る磁気共鳴イメージング装置において励起される２平面と柱状領域とを柱状領域の軸方向から見たときの空間的な関係を模式的に示す図である。
【図４】本発明の実施形態に係る磁気共鳴イメージング装置において励起される２平面と柱状領域とを模式的に示す斜視図である。
【図５】本発明の実施形態に係る磁気共鳴イメージング装置において励起される柱状領域で形成される湾曲面を模式的に示す斜視図である。
【図６】柱状領域の断面の形状を模式的に示す図である。
【図７】柱状領域のオーバーラップを模式的に示す図である。
【図８】本発明の実施形態に係る磁気共鳴イメージング装置においてデータ収集を説明する図であって、（ａ）は柱状領域がオーバーラップする方向のＦＯＶを満たすデータ、（ｂ）はスライスの厚み方向の信号強度プロファイルを示している。
【図９】本発明の実施形態に係る磁気共鳴イメージング装置において柱状領域の励起順を模式的に示す図である。
【図１０】本発明の実施形態に係る磁気共鳴イメージング装置の処理の流れを示すフローチャートである。
【図１１】本発明の実施形態に係る磁気共鳴イメージング装置において操作者によるスライス曲線の設定手法の一例を模式的に示す図である。
【図１２】本発明の実施形態に係る磁気共鳴イメージング装置における画像再構成の手順を模式的に示す図であって、（ａ）は操作者により設定された各柱状領域、（ｂ）は各柱状領域のＭＲ信号、（ｃ）はＭＲ信号の１次元フーリエ変換後に計算された絶対値、（ｄ）は脊椎に沿った曲面のコロナル画像を示している。
【図１３】本発明の実施形態に係る磁気共鳴イメージング装置において操作者により延伸されたスライス曲線の一例を模式的に示す図である。
【図１４】本発明の実施形態に係る磁気共鳴イメージング装置において設定されたスライス曲線で指定されたスライス曲面に関する曲面断層画像をモニタ画面に再構成画像として表示した模式図であって、（ａ）はスライス曲線を延伸する前の状態、（ｂ）はスライス曲線を延伸した後の撮像途中の状態、（ｃ）はスライス曲線を延伸した後の最終状態を示している。
【発明を実施するための形態】
【００１８】
以下、図面を参照して本発明の実施形態に係る磁気共鳴イメージング装置について詳細に説明する。
図１に示すように、磁気共鳴イメージング装置１は、核磁気共鳴現象（ＮＭＲ現象）を利用して、被検体Ｐの断層画像を得るもので、静磁場発生系２と、傾斜磁場発生系３と、送信系４と、受信系５と、制御部（制御手段）６と、信号処理部７とを備えている。
【００１９】
静磁場発生系２は、図示しないテーブルに載置された被検体Ｐに静磁場を与えるものであり、静磁場発生部８を備えている。静磁場発生部８は、被検体Ｐに強く均一な静磁場を発生させるものであり、被検体Ｐの周りのある広がりをもった空間に配置された磁石からなる。この磁石は、例えば永久磁石、または常電導方式の磁石あるいは超電導方式の磁石である。ここで、静磁場発生部８が常電導方式あるいは超電導方式の磁石の場合、静磁場発生系２は、静磁場発生部８に電源を供給する静磁場電源を備える。
【００２０】
傾斜磁場発生系３は、被検体Ｐに傾斜磁場を与えるものであり、傾斜磁場コイル１０と、傾斜磁場電源１１とを備えている。傾斜磁場コイル１０は、Ｘ、Ｙ、Ｚの三軸方向（第１のスライス方向、第２のスライス方向、周波数エンコード方向）に巻かれており、それぞれのコイルは傾斜磁場電源１１に接続されている。
【００２１】
傾斜磁場電源１１は、制御部６からの命令に従って傾斜磁場コイル１０を駆動するものである。これにより、傾斜磁場コイル１０は、Ｘ、Ｙ、Ｚの三軸方向の傾斜磁場Ｇｘ、Ｇｙ、Ｇｚを被検体Ｐに印加するようになっている。傾斜磁場Ｇｘ、Ｇｙ、Ｇｚを用い、第１のスライス方向傾斜磁場、第２のスライス方向傾斜磁場、周波数エンコード方向傾斜磁場を印加する。第１のスライス方向傾斜磁場と同時に印加する高周波パルスで励起される断面と、第２のスライス方向傾斜磁場と同時に印加される高周パルスで励起される断面と、の交差する柱状領域を励起することができる。また、周波数エンコード方向の傾斜磁場（柱状領域の長手方向の傾斜磁場）により、磁気共鳴信号(エコー信号、ＮＭＲ信号、以下、単にＭＲ信号ともいう)に柱状領域の長手方向の位置情報を付与することができる。
【００２２】
送信系４は、被検体Ｐの生体組織を構成する原子の原子核に核磁気共鳴（ＮＭＲ）を起こさせる高周波パルスを照射するものであり、高周波発生器１２と、変調器１３と、高周波アンプ１４と、送信側高周波コイル１５ａとを備えている。
【００２３】
高周波発生器１２は、高周波パルスを発生する。変調器１３は、高周波発生器１２から出力された高周波パルスを制御部６の命令に従って変調する。高周波アンプ１４は、変調された高周波パルスを増幅する。送信側高周波コイル１５ａは、被検体Ｐに近接して配置されており、増幅された高周波パルス（電磁波）を被検体Ｐに照射する。
【００２４】
受信系５は、被検体Ｐの生体組織における原子核の磁気共鳴により放出されるＮＭＲ信号を検出するものであり、受信側高周波コイル１５ｂと、増幅器１６と、Ａ／Ｄ変換器１７と、検波器１８とを備えている。
【００２５】
受信側高周波コイル１５ｂは、被検体Ｐに近接して配置されており、送信側高周波コイル１５ａから照射された電磁波である高周波パルスによる被検体Ｐからの応答の電磁波であるＮＭＲ信号を検出する。増幅器１６は、検出されたＮＭＲ信号を増幅する。Ａ／Ｄ変換器１７は、増幅されたＮＭＲ信号をデジタル信号に変換する。検波器１８は、送信系４の高周波発生器１２と接続されており、デジタル変換されたＮＭＲ信号を直交位相検波し、制御部６からの命令によるタイミングでリサンプリングした二系列の収集データとする。この二系列の収集データの信号は、制御部６から信号処理部７に送られるようになっている。
【００２６】
傾斜磁場発生系３の傾斜磁場コイル１０と、送信系４の送信側高周波コイル１５ａと、受信系５の受信側高周波コイル１５ｂとは、静磁場発生系２の静磁場発生部８の発生する静磁場の磁場空間内に設置されており、本実施形態における撮像手段９を構成している。
撮像手段９は、所定の撮像条件で被検体Ｐに核磁気共鳴信号（ＮＭＲ信号）を発生させると共に当該核磁気共鳴信号を検出する。
【００２７】
制御部６は、撮像手段９を制御するものであり、ＣＰＵ（Central Processing Unit）の制御で動作する。制御部６は、高周波パルスと傾斜磁場を所定のパルスシーケンス（pulse sequence）で繰り返し印加してＮＭＲ信号を計測データとして収集する。制御部６は、被検体Ｐの断層画像のデータ収集に必要な種々の命令を送信系４、傾斜磁場発生系３、または受信系５に送り、所定のパルスシーケンスで傾斜磁場および高周波パルスを発生させるとともにＮＭＲ信号を検出し、所定のＮＭＲ信号が収集されるようにする。
【００２８】
信号処理部７は、計算機１９と、ディスプレイ（表示手段）２０と、記憶部２１と、入力デバイス（入力手段）２２とを備えている。計算機１９は、例えば、一般的なコンピュータで実現することができ、ＣＰＵと、ＲＡＭ（Random Access Memory）と、ＲＯＭ（Read Only Memory）と、ＨＤＤ（Hard Disk Drive）と、入力／出力インタフェースとを含んで構成されている。計算機１９は、受信系５で検出されたＮＭＲ信号から画像を再構成してディスプレイ２０に表示するようになっている。
【００２９】
ディスプレイ２０は、再構成された画像である断層画像を表示するものであり、例えば、液晶ディスプレイ等から構成される。再構成された画像は記憶部２１に記憶されている。
【００３０】
記憶部２１は、断層画像等の画像データや各種データを記憶するものであり、例えば磁気ディスク、光ディスク、一般的な画像メモリ等から構成される。記憶部２１は、再構成された画像と入力された撮像条件とを組み合わせて記憶する。この撮像条件は、入力デバイス２２で入力される。
【００３１】
入力デバイス２２は、信号処理部７で行う処理や撮像条件や撮像する曲面の設定など制御のために必要な情報や被検体情報を入力するものであり、例えばキーボードやマウスあるいはタッチパネルと、条件設定に必要なＧＵＩ（Graphical User Interface）を表示する表示部などを備えている。
【００３２】
入力デバイス２２を介して設定される撮像条件として、パルスシーケンスの選択や、撮像する曲面（スライス曲面）の設定等が含まれる。
スライス曲面の設定は、例えばディスプレイ２０に、被検体Ｐの所定のスライス曲面に直交する断面画像としてスカウト画像が表示された状態で行われる。ここで、スカウト画像は、短時間で撮像された分解能が低い位置決め用の画像であって、例えばトランスバース、サジタルおよびコロナルのいずれかの画像を用いる。例えば、被検体の脊髄を含む曲面の撮像を行う場合、被検体の脊髄を含むスカウト画像を事前に撮像しておく。
【００３３】
入力デバイス２２は、表示されたスカウト画像に対して重畳するようにスライス曲面を指定するための曲線（スライス曲線）を入力する。
例えばスカウト画像として表示された矢状断面の画像（サジタル画像）上に、操作者がマウス等の入力デバイス２２を操作することにより曲線を描くことで、計算機１９が、この曲線を受け付けて、曲線を含み且つコロナル断面（冠状断面）に対応した湾曲した曲面を撮像すべき曲面として設定することができる。
【００３４】
計算機１９は、画像再構成処理部（画像再構成処理手段）２３と、位置計算部（位置計算手段）２４とを備えている。
画像再構成処理部２３は、撮像手段９により検出されたＮＭＲ信号を用いて被検体の画像を作成する。再構成された画像は、その撮像条件と共に記憶部２１に格納される。画像再構成処理部２３が行う画像再構成では、フーリエ変換を行い、絶対値画像を作成する。
【００３５】
位置計算部２４は、入力デバイス２２により入力されたスライス曲線を受け付けて、被検体Ｐにおける所定の柱状領域を、スライス曲線に沿って移動させるときのそれぞれの位置を算出するものである。位置計算部２４は、算出したそれぞれの位置情報およびスライス曲線の位置情報を保持する。
本実施形態では、位置計算部２４は、柱状領域を、スライス曲線に沿ってオーバーラップするように移動させる。
以下、柱状領域の励起によるエコー信号（磁気共鳴信号）の収集の詳細について説明する。
【００３６】
［パルスシーケンス］
まず、パルスシーケンスについて図２を参照（適宜図１および図３参照）して説明する。図２において横軸は時間を示し、縦軸は、各項目の強度を示している。
制御部６の制御の下、撮像手段９により行うパルスシーケンスの概要は次の通りである。
まず、撮像手段９は、磁化を９０度励起する第１の高周波パルス２０１と同時に、第１のスライス方向傾斜磁場２０３を被検体に印加する。このとき第１のスライス方向傾斜磁場２０３と同時に第２のスライス方向傾斜磁場２０５を被検体に印加する。これにより、被検体において、第１のスライス断面２２５（図３参照）が励起される。
【００３７】
次に、撮像手段９は、周波数エンコード方向傾斜磁場２０７を被検体に印加する。周波数エンコード方向傾斜磁場２０７は、検出すべきエコー信号の最大強度が時間軸上でＡ／Ｄ変換のための所定期間における中心付近に位置するようにするために印加するものである。
【００３８】
次に、撮像手段９は、磁化を１８０度励起する第２の高周波パルス２０２と同時に、第１のスライス方向傾斜磁場２０４を被検体に印加する。このとき第１のスライス方向傾斜磁場２０４と同時に第２のスライス方向傾斜磁場２０６を被検体に印加する。これにより第２のスライス断面２２６（図３参照）が励起される。
これにより、第１のスライス断面２２５と第２のスライス断面２２６とが交差する柱状領域２２７（図３においてハッチングで示す菱形領域）だけが、磁化を９０度励起する第１の高周波パルス２０１と、磁化を１８０度励起する第２の高周波パルス２０２との２つの高周波パルスによる励起を受けることになる。
【００３９】
次に、撮像手段９は、周波数エンコード方向傾斜磁場２０８を被検体に印加して、第１のスライス断面２２５と第２のスライス断面２２６とが交差する柱状領域２２７（図３参照）からのエコー信号２０９を検出する。そして、得られたエコー信号２０９は、タイミング２１０の期間でＡ／Ｄ変換器１７によりＡ／Ｄ変換される。なお、信号処理部７は、Ａ／Ｄ変換されたデータを用いて、画像再構成処理を行い、断層画像を表示する。
【００４０】
撮像手段９は、第１の高周波パルス２０１の印加からエコー信号２０９の検出までの動作を、操作者によって撮像条件として指定されたＴＲ（繰返し時間）で、次に励起する柱状領域の位置に切り替えて繰り返し実行する。ここで、第１の高周波パルス２０１の中心からエコー信号２０９の中心までの時間はＴＥ（エコー時間）である。
【００４１】
［柱状領域の励起］
次に、柱状領域の励起について図３を参照（適宜図２参照）して説明する。図３では、柱状領域の断面がＸＹ平面に配置されているものとしており、横軸は位置座標Ｘを示し、縦軸は位置座標Ｙを示している。
【００４２】
図３において、ベクトル２２０は、第１のスライス方向傾斜磁場２０３（図２参照）の印加方向と印加強度を示すベクトルである。また、ベクトル２２１は、第１のスライス方向傾斜磁場２０３と同時に印加される第２のスライス方向傾斜磁場２０５（図２参照）の印加方向と印加強度を示すベクトルである。第１のスライス方向傾斜磁場２０３と第２のスライス方向傾斜磁場２０５とは同時に被検体に印加されるので、被検体に実質的に印加される磁場は、ベクトル２２０とベクトル２２１とが合成された合成ベクトル２２３で表される傾斜磁場となる。ここで、励起させたい所望のスライス断面に対応した高周波周波数を設定された第１の高周波パルス２０１（図２参照）を印加すると、合成ベクトル２２３に直交した第１のスライス断面２２５が励起される。ここで、第１のスライス断面２２５の厚みは例えばＤであるものとする。
【００４３】
また、第１のスライス断面２２５を励起する高周波周波数の値と、合成ベクトル２２３等との間には次の関係がある。ここでは、合成ベクトル２２３に対応する傾斜磁場の強度をＧ1とする。また、座標の原点から第１のスライス断面２２５に下した垂線を想定する。このとき原点から第１のスライス断面２２５までの距離をＲ1とする。なお、第１のスライス断面２２５には厚みＤがあるので、距離Ｒ1とは、原点から第１のスライス断面２２５の厚みの中心点までの距離を意味する。この場合、第１のスライス断面２２５を励起する高周波周波数Ｆ1は次の式（１）で表される。
Ｆ1＝γ・(Ｇ1・Ｒ1＋Ｈ0)  …  式（１）
ここで、γは水素原子核の核磁気回転比であり、具体的には４２．５ＭＨｚ/Ｔである。また、Ｈ0は、静磁場強度である。
【００４４】
前記ベクトル２２０は、図３において、第１のスライス方向傾斜磁場２０４（図２参照）の印加方向と印加強度を示すベクトルも表している。また、ベクトル２２２は、第１のスライス方向傾斜磁場２０４と同時に印加される第２のスライス方向傾斜磁場２０６（図２参照）の印加方向と印加強度を示すベクトルである。第１のスライス方向傾斜磁場２０４と第２のスライス方向傾斜磁場２０６とは同時に被検体に印加されるので、被検体に実質的に印加される磁場は、ベクトル２２０とベクトル２２２とが合成された合成ベクトル２２４で表される傾斜磁場となる。ここで、前記した所望のスライス断面に対応した高周波周波数を設定された第２の高周波パルス２０２（図２参照）を印加すると、合成ベクトル２２４に直交した第２のスライス断面２２６が励起される。ここで、第２のスライス断面２２６の厚みは例えばＤであるものとする。
【００４５】
また、第２のスライス断面２２６を励起する高周波周波数の値と、合成ベクトル２２４等との間には次の関係がある。ここでは、合成ベクトル２２４に対応する傾斜磁場の強度をＧ2とする。また、座標の原点から第２のスライス断面２２６までの距離をＲ2とする。なお、第２のスライス断面２２６には厚みＤがあるので、距離Ｒ2とは、原点から第２のスライス断面２２６の厚みの中心点までの距離を意味する。この場合、第２のスライス断面２２６を励起する高周波周波数Ｆ2は次の式（２）で表される。
Ｆ2＝γ・(Ｇ2・Ｒ2＋Ｈ0)  …  式（２）
なお、γは水素原子核の核磁気回転比であり、Ｈ0は、静磁場強度である。
【００４６】
以上より、前記したように、第１のスライス断面２２５と第２のスライス断面２２６とが交差する柱状領域２２７だけが、所望のスライス断面を特定する周波数値が設定された第１の高周波パルス２０１と第２の高周波パルス２０２との両方で励起されるため、被検体において、所望のスライス断面に配置された柱状領域２２７だけからエコー信号が得られる。
【００４７】
ここで、入力されたスライス曲線に沿って柱状領域を移動させる方法は、被検体を寝台等に載せて物理的に移動させるのではなく、被検体に印加される高周波パルスに設定される高周波周波数の値を変化させることで実行できる。例えば、柱状領域を最初の所望の位置座標（Ｘ1，Ｙ1）に配置するには、第１回目に、位置座標（Ｘ1，Ｙ1）を特定する周波数値ｆ11が設定された第１の高周波パルス２０１と周波数値ｆ12が設定された第２の高周波パルス２０２との両方で被検体を励起する。そして、柱状領域を移動先である次の所望の位置座標（Ｘ2，Ｙ2）に配置するには、撮像条件として指定されたＴＲ（繰返し時間）後の第２回目に、位置座標（Ｘ2，Ｙ2）を特定する周波数値ｆ21が設定された第１の高周波パルス２０１と周波数値ｆ22が設定された第２の高周波パルス２０２との両方で被検体を励起する。以下、同様に周波数値を変化せることで、柱状領域を所望の位置座標（Ｘn，Ｙn）に配置することができる。
【００４８】
図３では、簡便なベクトル図で説明するために柱状領域を２次元空間で説明したが、柱状領域は、図４に示すような３次元領域である。図４では、柱状領域２２７の長手方向を図面上の奥行き方向に一致させて示している。なお、柱状領域２２７の長手方向の長さをＦＯＶ（１）と表記している。
【００４９】
［スライス曲面］
次に、スライス曲面について図５を参照（適宜図２参照）して説明する。図５では、柱状領域２２７の長手方向を図面の上下方向に一致させて配置している。
図５に示す柱状領域２２７を移動させること、すなわち、入力されたスライス曲線に沿ってオーバーラップするように配置された複数の柱状領域２２７を想定することは、被検体において励起する３次元領域（柱状領域２２７）の位置を変化させることに対応する。
【００５０】
そして、入力されたスライス曲線２２８に沿うように、被検体において柱状領域２２７が励起する位置を変える動作を、高周波パルスに設定される高周波周波数の値を変化させながら、所定の回数だけ、所定の繰返し時間（ＴＲ）で繰り返すことで、湾曲したスライス曲面２２９全体のエコー信号を得ることができる。そして、それら全体のエコー信号についてフーリエ変換等を行うことで、スライス曲面２２９に沿った曲面断層画像を再構成することができる。なお、以下では、曲面断層画像のＦＯＶのうち、柱状領域２２７の長手方向の長さをＦＯＶ（１）、柱状領域２２７の連鎖方向の長さをＦＯＶ（２）とする。
【００５１】
［柱状領域の形状］
次に、柱状領域の形状について図６を参照（適宜図３及び図５参照）して説明する。図６には、柱状領域の断面形状を示している。
本実施形態では、柱状領域２２７の断面の形状は菱形としている。また、図６では、菱形の水平方向の対角線（一方の対角線）が垂直方向の対角線（他方の対角線）よりも短くなるように図示されている。柱状領域２２７を形成する菱形の垂直方向の対角線の長さは、スライス厚Ｔの１．５〜２．５倍であることが好ましく、菱形の水平方向の対角線の長さは、画素サイズの１．５〜２．５倍であることが好ましい。ここで、画素サイズは、スライス曲線に沿って配置される柱状領域２２７の個数（スライス曲線に沿ったデータ点数）Ｎと、柱状領域２２７の連鎖方向の長さをＦＯＶ（２）（図５参照）によって決定される。すなわち、画素サイズは、ＦＯＶ（２）／Ｎで表される。
【００５２】
例えば、操作者（以下、ユーザーともいう）が撮像条件として、スライス厚をＴ（ｍｍ）、柱状領域２２７の連鎖方向のＦＯＶ（２）をＦＯＶ（ｍｍ）、柱状領域２２７の連鎖方向のデータ点数をＮとした場合、図６に示すように、柱状領域２２７を形成する菱形の垂直方向の対角線の長さが２×Ｔであり、かつ、水平方向の対角線の長さが２×ＦＯＶ／Ｎあることが特に好ましい。本実施形態では、上記の長さの対角線とした菱形の柱状領域２２７で励起することとする。
【００５３】
［柱状領域のオーバーラップ］
次に柱状領域のオーバーラップについて図７および図８を参照して説明する。
本実施形態では、菱形の柱状領域２２７を、柱状領域の連鎖方向にオーバーラップするような位置で励起する。この場合、隣接する菱形の柱状領域２２７同士の距離は、柱状領域の連鎖方向における対角線の長さの１／２であるＦＯＶ／Ｎとする。これにより菱形の柱状領域２２７は図７に示すようにオーバーラップする。なお、この場合、オーバーラップした領域の長手方向の長さはスライス厚Ｔとなる。
本実施形態では、柱状領域２２７を菱形の短い対角線の向きに移動させるので、画像空間分解能を高めることができる。
【００５４】
図８（ａ）は、菱形の柱状領域２２７における柱状領域の連鎖方向の対角線が、ユーザーによって入力されたスライス曲線２２８の接線と重なるように、且つ、図７と同様にオーバーラップして移動する場合の軌跡を模式的に示す図である。図８（ａ）は、柱状領域の連鎖方向におけるＦＯＶ（その長さはＦＯＶ（２））をすべて満たすようにデータを収集した状態を示している。ここでは、柱状領域を図６および図７に示した形状および配置としていることから、柱状領域の連鎖方向のＦＯＶを満たす菱形の柱状領域２２７の個数は、ユーザーが撮像条件として設定したデータ点数と同じ値（Ｎ）となる。この場合、ユーザーがスライス曲線２２８により指定したスライス曲面２２９に沿った断層画像が得られる。
【００５５】
図８（ｂ）は、図８（ａ）のスライス曲面２２９についての信号強度プロファイルを示している。図８（ｂ）において、縦軸は図８（ａ）のスライス曲面２２９の厚み方向の位置を示している。また、横軸は、ＭＲ信号強度を示している。なお、ＭＲ信号強度の値は横軸において左の方が右よりも高くなるように表されている。
この場合、図８（ａ）では、柱状領域が図６および図７に示した形状および配置としていることから、図８（ｂ）に示すプロファイルにおいて、ＭＲ信号の強度のＦＷＨＭ（半値幅）で評価すると、ＦＷＨＭは、ユーザーが撮像条件で設定したスライス厚Ｔ（ｍｍ）となる。
【００５６】
［柱状領域の励起順］
次に、柱状領域の好ましい励起順序について説明する。菱形の柱状領域２２７の位置を移動させながら曲面を撮像する際、隣接する柱状領域を連続的に励起していくことも可能である。ただし、柱状領域の磁化の緩和が十分に進んでいないうちに隣接する柱状領域を励起すると、干渉により、信号強度が低下してしまう。そのため、隣接する柱状領域を励起するためには、時間を空ける必要がある。このような待ち時間のため、撮像時間が長くなってしまう。そこで、撮像時間を短縮できるような柱状領域の励起順序について図９を参照して説明する。
【００５７】
以下では、例えば、ユーザーが撮像条件として、柱状領域の連鎖方向のＦＯＶを満たす菱形の柱状領域２２７の個数をＮに設定したものとして説明する。また、オーバーラップしながら配置される柱状領域の連続的に図９において左から順番に、Ｐ１、Ｐ２、Ｐ３、・・・Ｐ２５６と番号付けするものとする。さらに、ここでは、一例として６本の柱状領域を１つのグループとする。具体的には、Ｐ１〜Ｐ６が第１のグループ、Ｐ７〜Ｐ１２が第２のグループに属する。
【００５８】
この場合、第１のグループでは、励起の順番を、Ｐ１、Ｐ４、Ｐ２、Ｐ５、Ｐ３、Ｐ６として、３本おきに励起することで、隣接する柱状領域を連続励起しないようにする。
また、第２のグループでは、励起の順番を、Ｐ９、Ｐ１２、Ｐ８、Ｐ１１、Ｐ７、Ｐ１０として、３本おきに励起することで、隣接する柱状領域を連続励起しないようにする。
なお、第３グループ以降も、以下同様に、すべての柱状領域の励起を行う。
【００５９】
ここで、６本の柱状領域を３本おきに励起するので、隣接する柱状領域を励起するまでには、少なくとも（６／３）×ＴＲの時間が経過した後となるので、緩和が進んだ状態で、隣接領域を励起するため、干渉によるエコー信号の強度が低下することが抑制される。
【００６０】
［磁気共鳴イメージング装置の処理の流れ］
次に、磁気共鳴イメージング装置１の処理の流れについて図１０を参照（適宜図１、図２、図９および図１１〜図１４参照）して説明する。
まず、磁気共鳴イメージング装置１は、信号処理部７によって、ディスプレイ２０の画面上にスカウト画像を表示し、ユーザーは、そのスカウト画像上に撮像したい曲面を指定するためのスライス曲線を、マウス等の入力デバイス２２を用いて描く。これにより、信号処理部７の位置計算部２４は、スカウト画像に入力されたスライス曲線を受け付ける（ステップＳ１）。
【００６１】
図１１は、ユーザーがディスプレイ２０に表示されたスカウト画像に対して、例えばマウス等の入力デバイス２２を用い、撮像を行いたいスライス曲線２２８を設定した状態を示す模式図である。ここでは、スカウト画像は、一例として、被検体Ｐである犬の胸椎３０１から腰椎３０２まで脊髄液３０３と共に写っているサジタル画像であるものとした。
【００６２】
次に、磁気共鳴イメージング装置１は、このインタラクティブな曲面撮像を開始する（ステップＳ２）。磁気共鳴イメージング装置１は、例えば図１１に示す画像上で、ユーザーがスライス曲線２２８を設定した位置に関して、図２のパルスシーケンス図に従ったシーケンスを繰り返すことで、エコー信号を収集する。具体的には、撮像手段９は、位置計算部２４で算出されたそれぞれの位置に配置される柱状領域２２７を、傾斜磁場の印加および高周波磁場の照射により励起して柱状領域２２７から発生する磁気共鳴信号を検出する。また、制御部６は、柱状領域２２７の位置を所定の順序で変化させながら撮像手段９による一連の動作を撮像手段９に繰り返し行わせる。これにより、信号処理部７は、指定された曲面からのＭＲ信号の一部（第１グループのＰ１〜Ｐ６：図９参照）を得る（ステップＳ３）。
【００６３】
そして、信号処理部７は、画像再構成処理部２３によって、その得られたＭＲ信号の一部（Ｐ１〜Ｐ６）を用いて画像再構成を行う（ステップＳ４）。
そして、信号処理部７は、一部のＭＲ信号（Ｐ１〜Ｐ６）を用いて再構成された画像である曲面断層画像をディスプレイ２０に表示する（ステップＳ５）。ここで、ディスプレイ２０には、撮像したい曲面を指定するためのスライス曲線と共にスカウト画像が表示されているが、ディスプレイ２０の同じモニタ画面上に、再構成された曲面断層画像も表示される。
【００６４】
次に、指定された曲面からのＭＲ信号の収集が終了していない場合（ステップＳ６においてＮｏ）、磁気共鳴イメージング装置１は、ステップＳ３に戻り、信号処理部７は、指定された曲面からのＭＲ信号の一部（第２グループのＰ７〜Ｐ１２）を得る（ステップＳ３）。この場合、信号処理部７は、画像再構成処理部２３によって、その得られたＭＲ信号の一部（Ｐ７〜Ｐ１２）を用いて画像再構成を行う（ステップＳ４）。
【００６５】
そして、信号処理部７は、前回までに得られた一部のＭＲ信号（Ｐ１〜Ｐ６）を用いて再構成された曲面断層画像に加えて、今回得られた一部のＭＲ信号（Ｐ７〜Ｐ１２）を用いて再構成された曲面断層画像をディスプレイ２０に表示する（ステップＳ５）。
このようにして、信号処理部７は、各グループのデータが収集できた度に、曲面断層画像を更新していき、ユーザーが、撮像とリアルタイムに近い時間で曲面断層画像を確認できるようにする。
そして、指定されたスライス曲面における柱状領域から発生するＭＲ信号の収集が終了した場合（ステップＳ６においてＹｅｓ）、信号処理部７は、徐々に更新された曲面断層画像について現在の表示状態を保つ。
【００６６】
ここで、菱形の柱状領域２２７の位置を移動させながら収集したエコー信号の画像再構成の一例について図１２を参照して説明する。例えば、ユーザーが撮像条件の１つとして、柱状領域の連鎖方向のＦＯＶを満たす菱形の柱状領域２２７の個数（Ｎ）に、Ｎ＝２５６と設定したものとする。画像再構成処理部２３は、図１２（ａ）に示す各菱形の柱状領域２２７（Ｐ１〜Ｐ２５６）からのエコー信号（図１２（ｂ）参照）をそれぞれ１次元フーリエ変換し、その結果の絶対値を得る（図１２（ｃ）参照）。なお、図１２（ｃ）において、横軸は周波数、縦軸は強度を示している。
【００６７】
図２を参照して説明したように、撮像手段９は、エコー信号２０９の収集と同時に、周波数エンコード方向傾斜磁場２０８を被検体に印加している。周波数エンコード方向とは、柱状領域２２７の長手方向に対応している。したがって、１次元フーリエ変換して絶対値を計算した結果は、まさに、菱形の柱状領域２２７の断層画像情報である。これらの、１次元フーリエ変換して絶対値を計算した結果のデータ（図１２（ｃ）参照）を、２５６データ分並べれば、図１２（ｄ）に示すように、ユーザーが設定したスライス曲線２２８によるスライス曲面２２９（図１１参照）全体を示す曲面の断層画像（胸椎３０１から腰椎３０２にかけての脊椎に沿った曲面のコロナル画像）となる。
【００６８】
図１０に戻って磁気共鳴イメージング装置１の処理の流れについて説明を続ける。
例えば、指定されたスライス曲面における柱状領域から発生するＭＲ信号の収集が終了し（ステップＳ６においてＹｅｓ）、現在の表示状態を保っている場合、ユーザーは、ディスプレイ２０の画面上に表示された曲面断層画像を参考にして、図１１のスカウト画像上のスライス曲線を、マウス等の入力デバイス２２を用いて延長させることができる。これにより、信号処理部７の位置計算部２４は、スカウト画像に入力されたスライス曲線の延伸を受け付けることができる（ステップＳ７Ａ）。なお、スライス曲線の延伸ができるタイミングはこれに限定されず、撮像中でも、再構成中でもスライス曲線の延伸は画面上で随時受け付けている。
そして、例えば、スライス曲線が延伸されず（ステップＳ７ＢにおいてＮｏ）、撮像終了の指示がなければ（ステップＳ８においてＮｏ）、磁気共鳴イメージング装置１は、待機状態となる。
【００６９】
一方、ユーザーがスライス曲線を延伸した場合（ステップＳ７ＢにおいてＹｅｓ）、磁気共鳴イメージング装置１はステップＳ３に戻る。
図１３は、ユーザーが、スカウト画像上のスライス曲線を延伸した状態を示す模式図であり、破線部分が延伸したスライス曲線３０５である。図１３は、延伸したスライス曲線３０５が追加された点が図１１と異なっている。ここでは、例えば腰椎３０２に関して、背から腹にかけての横断面（アキシャル面）を観察したくなった場合に当該位置を含むようにユーザーがスライス曲線を延伸したこととしている。
このような場合であっても、磁気共鳴イメージング装置１は、延伸したスライス曲線３０５で指定されたスライス曲面における柱状領域から発生するＭＲ信号の収集と画像再構成、画像表示を前記の手順と同様に行う。
【００７０】
ここで、スライス曲線を延伸したときに、ディスプレイ２０に表示される曲面断層画像の遷移を図１４（ａ）〜図１４（ｂ）に示す。
図１４（ａ）は、図１１に示す実線のスライス曲線２２８で指定されるスライス曲面について再構成された曲面断層画像をモニタ画面に再構成画像として表示した模式図である。
図１４（ｂ）は、図１１に示す実線のスライス曲線２２８と、図１３に示す破線のスライス曲線３０５とで指定されるスライス曲面について撮像の途中までに再構成された曲面断層画像をモニタ画面に再構成画像として表示した模式図である。ここでは、腰椎３０２の横断面（アキシャル面）の一部まで画像再構成が行われている。
図１４（ｃ）は、図１４（ｂ）と同様であるが、延伸されたスライス曲線３０５で指定されるスライス曲面における柱状領域から発生するすべてのＭＲ信号が収集されて再構成された曲面断層画像をモニタ画面に再構成画像として表示した模式図である。
このように、磁気共鳴イメージング装置１は、撮像の進捗とリアルタイムに近い時間で再構成画像が増加し、これを反映して再構成された曲面断層画像を徐々に更新しながらディスプレイ２０に表示することができる。
【００７１】
図１０に戻って説明を続ける。
次に、延伸されたスライス曲線３０５に対応したＭＲ信号の収集が終了し、スライス曲線がそれ以上延長されず（ステップＳ７ＢにおいてＮｏ）、この曲面撮像を終了する場合（ステップＳ８においてＹｅｓ）、信号処理部７は、撮像により得られた画像を記憶部２１に保存した上で、曲面撮像を終了する（ステップＳ９）。
【００７２】
以上説明したように、磁気共鳴イメージング装置１は、被検体における柱状領域を、ユーザー（操作者）が指定したスライス曲線に沿ってオーバーラップするように移動させながら柱状領域から発生するＭＲ信号（磁気共鳴信号）を検出するので、得られる画像のＳＮ比が向上すると共に画像空間分解能を向上させることができる。
【００７３】
また、磁気共鳴イメージング装置１は、生体の湾曲した組織の形状に沿って曲面の断層画像を撮像し画像再構成する際に、撮像の進捗に従いリアルタイムに近い状態で曲面断層画像を確認できるように表示することができる。また、操作者が、ディスプレイ２０に表示された曲面断層画像を参考にして、スカウト画像上でスライス曲線をインタラクティブに延伸すると、磁気共鳴イメージング装置１は、これを反映して再構成された曲面断層画像を徐々に更新しながらディスプレイ２０に表示することができる。したがって、撮像を一旦中止して再度位置決めして撮像することなく、所望の位置における曲面の断層画像を容易に得ることができる。
【００７４】
なお、本発明は、前記実施形態に限定されるものではなく、その技術的思想の範囲内で種々の改造および変更が可能であり、本発明はこれら改造および変更された発明にも及ぶことは勿論である。例えば、被検体として犬を図示したが、他の動物でも人体でもよい。
また、スライス曲面として指定する部位を脊椎としたが、これに限るものではない。
また、スカウト画像はサジタル画像であるものとしたが、これに限らず、コロナルやトランスバースの画像であってもよい。
【符号の説明】
【００７５】
１      磁気共鳴イメージング装置
２      静磁場発生系
３      傾斜磁場発生系
４      送信系
５      受信系
６      制御部（制御手段）
７      信号処理部
８      静磁場発生部
９      撮像手段
１０    傾斜磁場コイル
１５ａ  送信側高周波コイル
１５ｂ  受信側高周波コイル
１９    計算機
２０    ディスプレイ（表示手段）
２１    記憶部
２２    入力デバイス（入力手段）
２３    画像再構成処理部（画像再構成処理手段）
２４    位置計算部（位置計算手段）
【図１】
2018139779.tif 000003
【図２】
2018139779.tif 000004
【図３】
2018139779.tif 000005
【図４】
2018139779.tif 000006
【図５】
2018139779.tif 000007
【図６】
2018139779.tif 000008
【図７】
2018139779.tif 000009
【図８】
2018139779.tif 000010
【図９】
2018139779.tif 000011
【図１０】
2018139779.tif 000012
【図１１】
2018139779.tif 000013
【図１２】
2018139779.tif 000014
【図１３】
2018139779.tif 000015
【図１４】
2018139779.tif 000016