(57)【要約】【課題】  スカウト画像の撮像から本撮像を行うまでに、被検体が動いたことを好適に検知することができる。【解決手段】  互いに法線の方向が異なるスカウト画像を2枚以上撮像する撮像部と、該スカウト画像を記憶する記憶部と、前記スカウト画像を表示する表示部と、該表示手段により表示されたスカウト画像上に、操作者が直線を入力することによって、本撮像を行う位置及び方向を入力可能とする操作部を備え、前記撮像制御部は、前記操作部により入力された位置及び方向で本撮像を行う磁気共鳴イメージング装置において、前記スカウト画像を撮像してから前記本撮像を撮像するまでに、前記被検体の向きが第1の角度から第2の角度へ動いたことを検出する位置ずれ検出部を備える。【選択図】          図2
2014046123.tif 000002
【特許請求の範囲】
【請求項１】
互いに法線の方向が異なる被検体のスカウト画像を２枚以上撮像する撮像制御部と、該スカウト画像を記憶する記憶部と、前記スカウト画像を表示する表示部と、該表示部により表示されたスカウト画像上に、操作者が直線を入力することによって、前記被検体の本撮像を行う位置及び方向を入力可能とする操作部を備え、
前記撮像制御部は、前記操作部により入力された位置及び方向で本撮像を行う磁気共鳴イメージング装置において、
前記スカウト画像を撮像してから前記本撮像を撮像するまでに、前記被検体の向きが第１の角度から第２の角度へ動いたことを検出する位置ずれ検出部と、該位置ずれ検出部により検出された位置ずれを知らせる報知部を備えたことを特徴とする磁気共鳴イメージング装置。
【請求項２】
前記スカウト画像の２つの法線は、前記磁気共鳴イメージング装置における直交する第１の方向と第２の方向であり、本撮像により取得する本撮像画像の法線は、前記第１の方向と第２の方向に直交する第３の方向であることを特徴とする請求項１に記載の磁気共鳴イメージング装置。
【請求項３】
前記検出部は、２枚のスカウト画像上における、前記直線上の前記被検体の領域の幅を、第１の長さ及び第２の長さとし、前記本撮像により得られる画像上での被検体領域の前記第２の方向への長さを第３の長さ、前記第１の方向への長さを第４の長さとし、前記第１の長さと前記第３の長さ、前記第２の長さと前記第４の長さを比較することにより、前記被検体が動いたことを検出することを特徴とする請求項２に記載の磁気共鳴イメージング装置。
【請求項４】
前記第１の長さ及び前記第２の長さは、横軸を前記直線上の位置、縦軸を前記直線上の画素値とするプロファイルを生成し、該プロファイルの裾野を所定の閾値に基づいて切り取り、切り取られなかった部分の前記直線に沿った長さを測ることにより、求められることを特徴とする請求項３に記載の磁気共鳴イメージング装置。
【請求項５】
前記第３の長さは、前記本撮像の際に得られた計測データより、前記第２の方向の中央１ラインをフーリエ変換してプロファイルを生成し、該プロファイルの裾野を所定の閾値に基づいて切り取り、切り取られなかった部分の前記プロファイルの横軸に沿った長さを測ることにより、求められ、
前記第４の長さは、前記本撮像の際に得られた計測データより、前記第１の方向の中央１ラインをフーリエ変換してプロファイルを生成し、該プロファイルの裾野を所定の閾値に基づいて切り取り、切り取られなかった部分の前記プロファイルの横軸に沿った長さを測ることにより、求められることを特徴とする請求項３に記載の磁気共鳴イメージング装置。
【請求項６】
前記第３の長さは、前記本撮像により得られる画像を２次元フーリエ変換して得られる計測データより、前記第２の方向の中央１ラインをフーリエ変換してプロファイルを生成し、該プロファイルの裾野を所定の閾値に基づいて切り取り、切り取られなかった部分の前記プロファイルの横軸に沿った長さを測ることにより、求められ、
前記第４の長さは、前記本撮像により得られる画像を２次元フーリエ変換して得られる計測データより、前記第１の方向の中央１ラインをフーリエ変換してプロファイルを生成し、該プロファイルの裾野を所定の閾値に基づいて切り取り、切り取られなかった部分の前記プロファイルの横軸に沿った長さを測ることにより、求められることを特徴とする請求項３に記載の磁気共鳴イメージング装置。
【請求項７】
前記第２の角度と前記第１の角度の差分を推定する差分角度推定部と、前記操作部により入力された前記被検体の所望の位置及び方向を磁気共鳴イメージング装置座標系で前記差分の角度だけ回転させて再撮像のための位置及び方向を求める計算部と、該計算部で求めた位置及び方向で再撮像する撮像制御部を備え、前記検出部が前記被検体の位置ずれを検出すると、前記差分角度推定部が前記差分を推定し、前記計算部が再撮像のための位置及び方向を計算し、前記撮像制御部が、再撮像をすることを特徴とする請求項４乃至６のいずれか一項に記載の磁気共鳴イメージング装置。
【請求項８】
前記回転は、前記被検体の前記第２の方向を中心とする回転であり、前記差分角度推定部は、前記第２の方向を法線とするスカウト画像を少しずつ回転させ、前記操作部により入力された本撮像の方向及び位置での直線上のプロファイルの縦軸が閾値以上の値を持つ横軸に沿った長さが、前記第４の長さと同じであるような角度を前記差分として求めることを特徴とする請求項７に記載の磁気共鳴イメージング装置。
【請求項９】
前記回転は、前記被検体の前記第３の方向を中心とする回転であり、前記差分角度推定部は、縦の長さが前記第２の長さ、横の長さが前記第１の長さとするマスクを少しずつ回転させ、該マスクが前記本撮像による画像と最も合致する角度を前記差分として求めることを特徴とする請求項７に記載の磁気共鳴イメージング装置。
【請求項１０】
前記回転は、前記被検体の前記第１の方向を中心とする回転であり、前記差分角度推定部は、前記第１の方向を法線とするスカウト画像を少しずつ回転させ、前記操作部により入力された本撮像の方向及び位置での直線上のプロファイルの縦軸が閾値以上の値を持つ横軸に沿った長さが、前記第３の長さと同じであるような角度を前記差分として求めることを特徴とする請求項７に記載の磁気共鳴イメージング装置。
【請求項１１】
前記スカウト画像は、法線が互いに直交する第１及び第２及び第３の３つの方向からなる第１及び第２及び第３のスカウト画像であり、本撮像による本撮像画像を取得すべき前記被検体の所望の位置及び方向は、前記被検体の斜断面であり、前記第２のスカウト画像上に、斜めの直線で前記操作部により入力され、
前記第２のスカウト画像の右側に、第３のスカウト画像が、前記第１の方向の位置が合うようにして配置され、
前記第２のスカウト画像上で、前記斜めの線が前記被検体の領域の境界と交わる点の一方を第１の点とし、前記第１のスカウト画像上の中央ラインで、前記第１の方向への直線を第１の中央線とし、該第１の中央線と前記斜めの線の交点を第２の点とし、
前記第１の点と前記第２の点の距離を幅１とし、
前記第２のスカウト画像上前記第３の方向へ第１の破線を前記第２の点から前記第３のスカウト画像上まで描き、前記第１の破線の前記第３のスカウト画像上での前記被検体領域の幅を第５の長さとし、
前記第２のスカウトを右側へ９０°回転させたものを回転した第２のスカウト画像とし、
前記第２のスカウト画像の右側に、前記第１のスカウト画像が、前記第３の方向の位置が合うように配置され、
前記回転した第２のスカウト画像上に、前記斜めの直線を描き、
該斜めの直線上に前記第１の点に相当するものを描き、
前記回転した第２のスカウト画像上の中央ラインで、前記第３の方向への直線を第２の中央線とし、前記第２の中央線と前記斜めの線の交点を第３の点とし、
前記第１の点と第３の点の距離を幅２とし、
前記回転した第２のスカウト画像上前記第１の方向へ第２の破線を前記第１のスカウト画像上まで描き、前記第２の破線の前記第１のスカウト画像上での前記被検体領域の幅を第６の長さとし、
前記本撮像により得られた画像を前記第１のスカウト画像と前記第３のスカウト画像の間に、前記第２の方向の位置が合うように配列され、
前記本撮像による画像の被検体領域上で、前記第１の点に相当する点を第４の点とし、
前記第４の点から、前記本撮像による画像上で縦方向に下側に前記幅１の距離にあり、
横方向の第２の直線を描き、その直線上の前記被検体領域の長さを第７の長さとし、
前記第４の点から、前記本撮像による画像上で縦方向に下側に前記幅２の距離にあり、
横方向の第３の直線を描き、その直線上の前記被検体領域の長さを第８の長さとし、
前記第５の長さと前記第７の長さ、前記第６の長さと前記第８の長さを比較することにより、前記被検体が第１の角度から第２の角度へ動いたことを検出することを特徴とする請求項１に記載の磁気共鳴イメージング装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、磁気共鳴イメージング(以下、MRIという。)装置に係り、特に、スカウト画像の撮像から本撮像を行うまでに、被検体が動いたことを好適に検知することが可能なMRI装置に関する。
【背景技術】
【０００２】
MRI装置は、磁場中に置かれた被検体の核磁気共鳴(以下「NMR」という)現象から得られる信号を計測し演算処理することにより、被検体中の核スピンの密度分布、緩和時間分布等を断層像として画像表示するものであり、人体を被検体として各種の診断等に使用されている。NMR現象から信号を得るためには、空間的、時間的に一様な静磁場中に被検体を置き、高周波コイルによりパルス状に電磁波を被検体に照射し、それによって発生するNMR信号を高周波コイルにより受信する。さらにNMR信号に位置情報を与えるために静磁場に傾斜磁場が重畳される。
【０００３】
MRI装置において、本撮像を行う撮像断面の位置を設定するために、スカウト画像を取得し、モニタ上に表示したスカウト画像上にROIや線分を入力することは、通常行われている(例えば、特許文献1、図8参照)。スカウト画像は、特許文献1、図8に示されているように、Axial像(横断面：人体の体軸に垂直な面)、Sagittal像(矢状断面：人体の額の面と垂直で、体軸と平行面)、Coronal像(冠状断面：人体の額と体軸の両方に平行面)の3種類があり、撮像断面の位置を設定するためには、これら3種類のスカウト画像の内少なくとも1枚以上が用いられる。特許文献1では特に、検査中に被検体が動いて、撮像目的部位が設定済みのスライス位置からずれた場合でも、予め取得された3次元画像データと、被検体が動いた後に撮像された画像の少なくとも一部のデータと、から被検体が動く前のスライス位置を求め、被検体が動く前のスライス位置を撮像する技術が開示されている。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特開2010-57532号公報
【発明の概要】
【発明が解決しようとする課題】
【０００５】
しかしながら、特許文献1では被検体が動いたことを検知するために3次元画像データを必要とするが、3次元画像データを取得するためには時間がかかる問題点があった。
【０００６】
そこで、本発明の目的は、スカウト画像の撮像から本撮像を行うまでに、被検体が動いたことを好適に検知することが可能なMRI装置を提供することにある。
【課題を解決するための手段】
【０００７】
上記の課題を解決するために、本発明の磁気共鳴イメージング装置は、スカウト画像を撮像してから本撮像を撮像するまでに、被検体の向きが第1の角度から第2の角度へ動いたことを検出する位置ずれ検出部と、該位置ずれ検出部により検出された位置ずれを知らせる報知部を備える。
【発明の効果】
【０００８】
本発明の磁気共鳴イメージング装置によれば、スカウト画像の撮像から本撮像を行うまでに、被検体が動いたことを好適に検知することが可能となる。
【図面の簡単な説明】
【０００９】
【図１】本発明に係るMRI装置の全体構成を示すブロック図
【図２】本発明の実施例1に係るSagittal方向の本撮像において、位置ずれが無い状態、および位置ずれ検出の方法を示す図
【図３】本発明の実施例1に係るSagittal方向の本撮像において、頭部が左に傾いた状態、および再撮像を示す図
【図４】本発明の実施例1に係るSagittal方向の本撮像において、顎を引いた状態、および再撮像を示す図
【図５】本発明の実施例1に係るSagittal方向の本撮像において、左を向いた状態、および再撮像を示す図
【図６】(a)本発明の実施例1のブロック図。(b)は(a)の位置ずれ検出部の構成を示すブロック図であり、(c)は(a)の再撮像位置設定部の構成を示すブロック図
【図７】本発明の実施例1に係るSagittal方向の本撮像における位置ずれ検出および再撮像について説明するフローチャート
【図８】本発明の実施例2に係る斜断面の本撮像において、位置ずれが無い状態、および位置ずれ検出の方法を示す図
【図９】本発明の実施例2に係る斜断面の本撮像において、頭部が左に傾いた状態を示す図
【図１０】本発明の実施例2に係る斜断面の本撮像において、顎を引いた状態を示す図
【図１１】本発明の実施例2に係る斜断面の本撮像において、左を向いた状態を示すブロック図
【図１２】本発明の実施例2に係る斜断面の本撮像において、再撮像を説明する図
【図１３】(a)は本発明の実施例2のブロック図を示す図。(b)は(a)の位置ずれ検出部の構成を示すブロック図であり、(c)は(a)の再撮像位置取得部の構成を示すブロック図
【図１４】本発明の実施例2に係る斜断面の本撮像における位置ずれ検出および再撮像について説明するフローチャート
【発明を実施するための形態】
【００１０】
本発明のMRI装置は、互いに法線の方向が異なるスカウト画像を2枚以上撮像する撮像制御部と、該スカウト画像を記憶する記憶部と、スカウト画像を表示する表示部と、該表示手段により表示されたスカウト画像上に、操作者が直線を入力することによって、本撮像を行う位置及び方向を入力可能とする操作部を備え、撮像制御部は、操作部により入力された位置及び方向で本撮像を行う磁気共鳴イメージング装置において、スカウト画像を撮像してから本撮像を撮像するまでに、被検体の向きが第1の角度から第2の角度へ動いたことを検出する位置ずれ検出部と、該位置ずれ検出部により検出された位置ずれを知らせる報知部を備える。
【００１１】
このように、スカウト画像を撮像してから本撮像を撮像するまでに、被検体の向きが第1の角度から第2の角度へ動いたことを検出する位置ずれ検出部を備えることにより、被検体がスカウト画像を撮像してから本撮像を撮像するまでに動いたことを検知できる。
【００１２】
また、本発明の別の態様のMRI装置では、スカウト画像の2つの法線は、磁気共鳴イメージング装置系における直交する第1の方向と第2の方向であり、本撮像により取得する本撮像画像の法線は、第1の方向と第2の方向に直交する第3の方向であることを特徴とする。
【００１３】
この第2の態様のMRI装置では、被検体の動きがより容易に検知される。
【００１４】
また、本発明の第3の態様のMRI装置では、検出部は、2枚のスカウト画像上における、直線上の被検体の領域の幅を、第1の長さ及び第2の長さとし、本撮像により得られる画像上での被検体領域の第2の方向への長さを第3の長さ、第1の方向への長さを第4の長さとし、第1の長さと第3の長さ、第2の長さと第4の長さを比較することにより、被検体が動いたことを検出することを特徴とする。
【００１５】
この第3の態様のMRI装置では、被検体の動きがより容易化される。
【００１６】
また、本発明の第4の態様のMRI装置では、第1の長さ及び第2の長さは、横軸を直線上の位置、縦軸を直線上の画素値とするプロファイルを生成し、該プロファイルの裾野を所定の閾値に基づいて切り取り、切り取られなかった部分の直線に沿った長さを測ることにより、求められることを特徴とする。
【００１７】
この第4の態様のMRI装置では、よりスカウト画像上での被検体の領域の幅がより明確に検知される。
【００１８】
また、本発明の第5の態様のMRI装置では、第3の長さは、本撮像の際に得られた計測データより、第2の方向の中央1ラインをフーリエ変換してプロファイルを生成し、該プロファイルの裾野を所定の閾値に基づいて切り取り、切り取られなかった部分のプロファイルの横軸に沿った長さを測ることにより、求められ、
第4の長さは、本撮像の際に得られた計測データより、第1の方向の中央1ラインをフーリエ変換してプロファイルを生成し、該プロファイルの裾野を所定の閾値に基づいて切り取り、切り取られなかった部分のプロファイルの横軸に沿った長さを測ることにより、求められることを特徴とする。
【００１９】
この第5の態様のMRI装置では、より本撮像により得られる画像上での被検体領域の幅がより明確に検知される。
【００２０】
また、本発明の第6の態様のMRI装置では、第3の長さは、本撮像により得られる画像を2次元フーリエ変換して得られる計測データより、第2の方向の中央1ラインをフーリエ変換してプロファイルを生成し、該プロファイルの裾野を所定の閾値に基づいて切り取り、切り取られなかった部分のプロファイルの横軸に沿った長さを測ることにより、求められ、
第4の長さは、本撮像により得られる画像を2次元フーリエ変換して得られる計測データより、第1の方向の中央1ラインをフーリエ変換してプロファイルを生成し、該プロファイルの裾野を所定の閾値に基づいて切り取り、切り取られなかった部分のプロファイルの横軸に沿った長さを測ることにより、求められる。この第6の態様のMRI装置では、より本撮像により得られる画像上での被検体領域の幅がより明確に検知される。
【００２１】
また、本発明の第7の態様のMRI装置では、第2の角度と第1の角度の差分を推定する差分角度推定部と、操作部により入力された被検体の所望の位置及び方向を磁気共鳴イメージング装置座標系で差分の角度だけ回転させて再撮像のための位置及び方向を求める計算部と、該計算部で求めた位置及び方向で再撮像する撮像制御部を備え、検出部が被検体の位置ずれを検出すると、差分角度推定部が差分を推定し、計算部が再撮像のための位置及び方向を計算し、撮像制御部が、再撮像をする。この第7の態様のMRI装置では、より再撮像のための撮像位置が明確に検知される。
【００２２】
また、本発明の第8の態様のMRI装置では、回転は、被検体の第2の方向を中心とする回転であり、差分角度推定部は、第2の方向を法線とするスカウト画像を少しずつ回転させ、操作部により入力された本撮像の方向及び位置での直線上のプロファイルの縦軸が閾値以上の値を持つ横軸に沿った長さが、第4の長さと同じであるような角度を差分として求めることを特徴とする。この第8の態様のMRI装置では、より再撮像のための撮像位置が明確に検知される。
【００２３】
また、本発明の第9の態様のMRI装置では、回転は、被検体の第3の方向を中心とする回転であり、差分角度推定部は、縦の長さが第2の長さ、横の長さが第1の長さとするマスクを少しずつ回転させ、該マスクが本撮像による画像と最も合致する角度を差分として求めることを特徴とする。この第9の態様のMRI装置によれば、より再撮像のための撮像位置が明確に検知される。
【００２４】
また、本発明の第10の態様のMRI装置では、回転は、被検体の第1の方向を中心とする回転であり、差分角度推定部は、第1の方向を法線とするスカウト画像を少しずつ回転させ、操作部により入力された本撮像の方向及び位置での直線上のプロファイルの縦軸が閾値以上の値を持つ横軸に沿った長さが、第3の長さと同じであるような角度を差分として求めることを特徴とする。この第10の態様のMRI装置によれば、より再撮像のための撮像位置が明確に検知される。
【００２５】
また、本発明の第11の態様のMRI装置では、スカウト画像は、法線が互いに直交する第1及び第2及び第3の3つの方向からなる第1及び第2及び第3のスカウト画像であり、本撮像による本撮像画像を取得すべき被検体の所望の位置及び方向は、被検体の斜断面であり、第2のスカウト画像上に、斜めの直線で操作部により入力され、第2のスカウト画像の右側に、第3のスカウト画像が、第1の方向の位置が合うようにして配置され、第2のスカウト画像上で、斜めの線が被検体の領域の境界と交わる点の一方を第1の点とし、第2のスカウト画像上の中央ラインで、第1の方向への直線を第1の中央線とし、該第1の中央線と斜めの線の交点を第2の点とし、第1の点と第2の点の距離を幅1とし、第2のスカウト画像上第3の方向へ第1の破線を第2の点から第3のスカウト画像上まで描き、第1の破線の第2のスカウト画像上での被検体領域の幅を第5の長さとし、第2のスカウト画像を右側へ90°回転させたものを回転した第2のスカウト画像とし、回転した第2のスカウト画像の右側に、第1のスカウト画像が、第3の方向の位置が合うように配置され、回転した第2のスカウト画像上に、斜めの直線を描き、該斜めの直線上に第1の点に相当するものを描き、回転した第2のスカウト画像上の中央ラインで、第3の方向への直線を第2の中央線とし、第2の中央線と斜めの線の交点を第3の点とし、第1の点と第3の点の距離を幅2とし、回転した第2のスカウト画像上第1の方向へ第2の破線を第1のスカウト画像上まで描き、第2の破線の第1のスカウト画像上での被検体領域の幅を第6の長さとし、本撮像により得られた画像を第3のスカウト画像と第1のスカウト画像の間に、第2の方向の位置が合うように配列され、本撮像による画像の被検体領域上で、第1の点に相当する点を第4の点とし、第4の点から、本撮像による画像上で縦方向に下側に幅1の距離にあり、横方向の第2の直線を描き、その直線上の被検体領域の長さを第7の長さとし、第4の点から、本撮像による画像上で縦方向に下側に幅2の距離にあり、横方向の第3の直線を描き、その直線上の被検体領域の長さを第8の長さとし、第5の長さと第7の長さ、第6の長さと第8の長さを比較することにより、被検体が第1の角度から第2の角度へ動いたことを検出する。
【００２６】
この第11の態様のMRI装置では、本撮像の方向が斜断面である場合において、被検体の動きがより容易に検知される。
【００２７】
以下、添付図面に従って本発明のMRI装置の好ましい実施形態について詳説する。なお、発明の実施形態を説明するための全図において、同一機能を有するものは同一符号を付け、その繰り返しの説明は省略する。
【００２８】
最初に、本発明に係るMRI装置の一例の全体概要を図1に基づいて説明する。図1は、本発明に係るMRI装置の一実施例の全体構成を示すブロック図である。このMRI装置は、NMR現象を利用して被検体の断層画像を得るもので、図1に示すように、MRI装置は静磁場発生系2と、傾斜磁場発生系3と、送信系5と、受信系6と、信号処理系7と、シーケンサ4と、中央処理装置(CPU)8とを備えて構成される。
【００２９】
静磁場発生系2は、垂直磁場方式であれば、被検体1の周りの空間にその体軸と直交する方向に、水平磁場方式であれば、体軸方向に均一な静磁場を発生させるもので、被検体1の周りに永久磁石方式、常電導方式あるいは超電導方式の静磁場発生源が配置されている。
【００３０】
傾斜磁場発生系3は、MRI装置の座標系(静止座標系)であるX，Y，Zの3軸方向に傾斜磁場を発生する傾斜磁場コイル9と、それぞれの傾斜磁場コイルを駆動する傾斜磁場電源10とから成り、後述のシ−ケンサ4からの命令に従ってそれぞれのコイルの傾斜磁場電源10を駆動することにより、X，Y，Zの3軸方向に傾斜磁場Gx，Gy，Gzを印加する。撮影時には、スライス面(撮影断面)に直交する方向にスライス方向傾斜磁場パルス(Gs)を印加して被検体1に対するスライス面を設定し、そのスライス面に直交して且つ互いに直交する残りの2つの方向に位相エンコード方向傾斜磁場パルス(Gp)と周波数エンコード方向傾斜磁場パルス(Gf)を印加して、エコー信号にそれぞれの方向の位置情報をエンコードする。
【００３１】
シーケンサ4は、高周波磁場パルス(以下、「RFパルス」という)と傾斜磁場パルスをある所定のパルスシーケンスで繰り返し印加する制御手段で、CPU8の制御で動作し、被検体1の断層画像のデータ収集に必要な種々の命令を送信系5、傾斜磁場発生系3、および受信系6に送る。
【００３２】
送信系5は、被検体1の生体組織を構成する原子の原子核スピンに核磁気共鳴を起こさせるために、被検体1にRFパルスを照射するもので、高周波発振器11と変調器12と高周波増幅器13と送信側の高周波コイル(送信コイル)14aとから成る。高周波発振器11から出力された高周波パルスをシーケンサ4からの指令によるタイミングで変調器12により振幅変調し、この振幅変調された高周波パルスを高周波増幅器13で増幅した後に被検体1に近接して配置された高周波コイル14aに供給することにより、RFパルスが被検体1に照射される。
【００３３】
受信系6は、被検体1の生体組織を構成する原子核スピンの核磁気共鳴により放出されるエコー信号(NMR信号)を検出するもので、受信側の高周波コイル(受信コイル)14bと信号増幅器15と直交位相検波器16と、A/D変換器17とから成る。送信側の高周波コイル14aから照射された電磁波によって誘起された被検体1の応答のNMR信号が被検体1に近接して配置された高周波コイル14bで検出され、信号増幅器15で増幅された後、シーケンサ4からの指令によるタイミングで直交位相検波器16により直交する二系統の信号に分割され、それぞれがA/D変換器17でディジタル量に変換されて、信号処理系7に送られる。
【００３４】
信号処理系7は、各種データ処理と処理結果の表示及び保存等を行うもので、光ディスク19、磁気ディスク18等の外部記憶装置と、CRT等からなるディスプレイ20とを有し、受信系6からのデータがCPU8に入力されると、CPU8が信号処理、画像再構成等の処理を実行し、その結果である被検体1の断層画像をディスプレイ20に表示すると共に、外部記憶装置の磁気ディスク18等に記録する。
【００３５】
操作部25は、MRI装置の各種制御情報や上記信号処理系7で行う処理の制御情報を入力するもので、トラックボール又はマウス23、及び、キーボード24から成る。この操作部25はディスプレイ20に近接して配置され、操作者がディスプレイ20を見ながら操作部25を通してインタラクティブにMRI装置の各種処理を制御する。
【００３６】
なお、送信側の高周波コイル14aと傾斜磁場コイル9は、被検体1が挿入される静磁場発生系2の静磁場空間内に、垂直磁場方式であれば被検体1に対向して、水平磁場方式であれば被検体1を取り囲むようにして設置されている。また、受信側の高周波コイル14bは、被検体1に対向して、或いは取り囲むように設置されている。
【００３７】
現在MRI装置の撮像対象核種は、臨床で普及しているものとしては、被検体の主たる構成物質である水素原子核(プロトン)である。プロトン密度の空間分布や、励起状態の緩和時間の空間分布に関する情報を画像化することで、人体頭部、腹部、四肢等の形態または、機能を2次元もしくは3次元的に撮像する。
【００３８】
次に、本発明の実施例の説明に入る前に、MRI装置における本撮像までの一連の流れについて説明する。MRI装置では、本撮像の前にスカウト撮像を行う。スカウト撮像とは、本撮像を行う位置を決めるために行う前計測のことであり、通常Axial(横断面)、Coronal(冠状断面)、Sagittal(矢状断面)の3断面を撮像するものである。操作者は、これら3断面のうち少なくとも1枚以上のスカウト画像を用いその画像上に、線分を入力することにより、本撮像をする位置を決める。
【実施例１】
【００３９】
次に、実施例1に係るMRI装置を図2〜図4を用い説明する。
【００４０】
ここでは、本実施例(本撮像の方向がSagittal方向)における位置ずれの態様についてその概略を説明する。ただし、本実施例における位置ずれとは、スカウト撮像から本撮像までに被検体が動くことによって、実際と異なる部位が撮像されてしまうことを言う。
【００４１】
先ず、図2は、位置決めに使用されたスカウト画像の撮像方向がAxial方向もしくはCoronal方向であり、本撮像の方向がSagittal方向である場合において、スカウト撮像から本撮像までに被検体1が動かない場合の例を示す図である。ただし、図2において、26aは、Axial方向のスカウト画像、27aは、Coronal方向のスカウト画像、28aはSagittal方向の本撮像により得られた画像を示す。なお、図2では、Axial方向のスカウト画像26aとCoronal方向のスカウト画像27aの被検体領域の横方向の位置が一致するように配置されており、Axial方向のスカウト画像26aの本撮像位置における被検体領域の縦幅とSagittal方向の本撮像画像28aの被検体領域の横幅が対比可能なように破線が描かれており、Coronal方向のスカウト画像27aの本撮像位置における被検体領域の縦幅と、Sagittal方向の本撮像画像28aの被検体領域の縦幅が比較可能なように配置されている。
【００４２】
これによれば、被検体1がスカウト撮像から本撮像までに動かない場合、Axial方向のスカウト画像26aの本撮像位置における被検体領域の縦幅とSagittal方向の本撮像画像28aの被検体領域の横幅とが一致し、Coronal方向のスカウト画像27aの本撮像位置における被検体領域の縦幅とSagittal方向の本撮像画像28aの被検体領域の縦幅とが一致していることがわかる。また、Sagittal方向の本撮像画像28aの被検体領域は、Axial方向のスカウト画像26aとCoronal方向のスカウト画像27aの本撮像位置における被検体領域の縦幅によって描かれる破線の領域の中に納まっていることがわかる。
【００４３】
次に、位置決めに使用されたスカウト画像の撮像方向がAxial方向もしくはCoronal方向であり、本撮像の方向がSagittal方向である場合において、被検体1がスカウト撮像から本撮像までに動いた場合、本撮像画像上の被検体領域がどのように変化するかを説明する。
【００４４】
先ず図3を用い、スカウト撮像から本撮像までに被検体1の頭部が左に傾いた場合(以下、Z軸を中心に回転した場合と表現する)、本撮像画像上の被検体領域がどのように変化するかを説明する。ただし、図3において26bは、Axial方向のスカウト画像、27bは、Coronal方向のスカウト画像、28bは、Sagittal方向の本撮像により得られた画像であり、それぞれの画像上に、被検体1の頭部がZ軸を中心に回転した場合、被検体1の頭部がどのように描出されるかが被検体1の頭部を表す模式図あるいは斜線で示されている。なお、図3の破線は、被検体1がスカウト撮像から本撮像までに動かない場合のSagittal方向の本撮像画像28aの被検体領域の横幅と縦幅であり、各スカウト画像内には、この場合の被検体1の頭部を表す模式図が点線で描かれている。これによれば、被検体1の頭部がZ軸を中心に回転した場合、本撮像位置によって、Sagittal方向の本撮像画像28bの被検体領域の縦幅には、長短の差ができ、被検体領域の横幅は、短くなることがわかる。
【００４５】
次に、図4を用い、スカウト撮像から本撮像までに被検体1の頭部が顎を引いた場合(以下、X軸を中心に回転した場合と表現する)、本撮像画像上の被検体領域がどのように変化するかを説明する。ただし、図4において26cは、Axial方向のスカウト画像、27cは、Coronal方向のスカウト画像、28cは、Sagittal方向の本撮像により得られた画像であり、それぞれの画像上に、被検体1の頭部がX軸を中心に回転した場合、被検体1の頭部がどのように描出されるかが被検体1の頭部を表す模式図あるいは斜線で示されている。なお、図4の破線は、被検体1がスカウト撮像から本撮像までに動かない場合のSagittal方向の本撮像画像28aの被検体領域の横幅と縦幅であり、各スカウト画像内には、この場合の被検体1の頭部を表す模式図が薄線で描かれている。これによれば、被検体1の頭部がX軸を中心に回転した場合、本撮像位置によって、Sagittal方向の本撮像画像28cの被検体領域の縦幅は、短くなり、被検体領域の横幅は、長くなることがわかる。
【００４６】
次に、図5を用い、スカウト撮像から本撮像までに被検体1の頭部が左を向いた場合(以下、Y軸を中心に回転した場合と表現する)、本撮像画像上の被検体領域がどのように変化するかを説明する。ただし、図5において26dは、Axial方向のスカウト画像、27dは、Coronal方向のスカウト画像、28dはSagittal方向の本撮像により得られた画像であり、それぞれの画像上に、被検体1の頭部がY軸を中心に回転した場合、被検体1の頭部がどのように描出されるかが被検体1の頭部を表す模式図あるいは斜線で示されている。なお、図5の破線は、被検体1がスカウト撮像から本撮像までに動かない場合のSagittal方向の本撮像画像28dの被検体領域の横幅と縦幅であり、各スカウト画像内には、この場合の被検体1の頭部を表す模式図が破線で描かれている。これによれば、被検体1の頭部がY軸を中心に回転した場合、本撮像位置によって、Sagittal方向の本撮像画像28dの被検体領域の縦幅は、短くなり、被検体領域の横幅には、長短の差ができることがわかる。
【００４７】
以上より、位置決めに使用されたスカウト画像の撮像方向がAxial方向もしくはCoronal方向であり、本撮像の方向がSagittal方向である場合において、位置ずれの特徴を説明したので、この特徴を用いた本実施例における位置ずれの検知の手法について説明する。この方法は、上述した直交座標系の3軸を中心とした回転による位置ずれを検知する方法である。ここでは、本撮像で描出される被検体領域の幅(縦方向と横方向)が、本撮像と異なる2方向から撮像されたスカウト画像における該幅と一致するかにより、位置ずれを検知する。また、2つのスカウト画像における前記幅を用いて本撮像画像上に長方形を描いた領域内に、本撮像により描出された被検体領域がどの程度含まれているかにより検知する。
【００４８】
以下、位置決めに使用されたスカウト画像の撮像方向がAxial方向もしくはCoronal方向であり、本撮像の方向がSagittal方向である場合において、位置ずれの検知方法の詳細を図2を用い説明する。
【００４９】
先ず、Axial方向のスカウト画像26aと、Coronal方向のスカウト画像27a上での本撮像位置の被検体領域の幅を各々算出する。より具体的には、まず、Axial方向のスカウト画像26a上あるいは、Coronal方向のスカウト画像27a上に、本撮像位置としてGUIを用いて直線Lが描かれると、該直線が横切る画素の画素値を縦軸とするプロファイルを生成する。続いて、そのプロファイルの裾野を所定の画素値の閾値(例えば被検体領域の画素値の1/20)に基づいて削除し、最初に被検体領域が確認できる点を点A1、点B1と定義し、この線分A1B1の長さを被検体領域幅L1aとして長さを求める。ここで、GUI上に描かれた本撮像位置を表す直線が、2つのスカウト画像のいずれか一方のみ上にしか配置されていないような場合、2つのスカウト画像を並列して配置し、一方に配置されている直線を延長させて、双方のスカウト画像上でSagittal方向の本撮像画像28aと比較する本撮像位置(直線Lの位置)の被検体領域幅を見積もるようにする。そして、得られたAxial方向のスカウト画像26a上の本撮像位置における被検体領域上の直線の長さをL1a、Coronal方向のスカウト画像27a上の本撮像位置における被検体領域上の直線の長さをL2aとする。
【００５０】
次に、Sagittal方向の本撮像画像28a上での被検体領域の縦幅と横幅を求める。本撮像による画像は、ノイズにより画像そのものからの割り出しの難しいケースが考えられる。この為、画像データを2次元フーリエ変換して、得られるデータの行方向と列方向の各中央ラインの信号を利用するようにする。より具体的には、画像データを2次元フーリエ変換して得られるデータから、行方向の中央に配置された1行を抽出し、それを1次元フーリエ変換して、基準位相エンコードプロファイル29aとする。また、画像データから2次元フーリエ変換して得られるデータから、列方向の中央に配置された1列を抽出し、それを1次元フーリエ変換して、基準周波数エンコードプロファイル30aとする。なお、このプロファイル生成については、画像再構成前に行っておいても良い。
【００５１】
この基準位相エンコードプロファイル29aと基準周波数エンコードプロファイル30aより、Sagittal方向の本撮像画像28aの被検体領域の横幅、縦幅を算出する。より具体的には、基準位相エンコードプロファイル29aについては、Sagittal方向の本撮像画像28aの列位置を示す横軸の両端からプロファイルをたどり、最初に被検体領域が確認できる点(あらかじめ定めた閾値以上の度数が確認できる点)を点A1’、点B1’とし、この線分A1’B1’の幅を被検体領域幅L1a’とする。また、基準周波数エンコードプロファイル30aについては、Sagittal方向の本撮像画像28aの行位置を示す縦軸の両端からプロファイルをたどり、最初に被検体領域が確認できる点(あらかじめ定めた閾値以上の度数が確認できる点)を点C1’、点D1’とし、この線分C1’D1’の幅を被検体領域幅L2a’とする。
【００５２】
ここで、算出した被検体領域幅L1a’、L2a’を先に算出したL1a、L2aと比較し、予め定めた閾値以上の差があれば位置ずれが生じたこととする。
【００５３】
一方、図2で、Axial方向のスカウト画像26a及び、Coronal方向のスカウト画像27aにおいて規定される被検体領域幅L1a及び被検体領域幅L2a(破線により規定される幅)により、Sagittal方向の本撮像画像28a上に規定される領域(太枠31aで示される領域)内に、実際にどの程度の被検体領域が写っているかを見積もり、その割合が所定の閾値以下であった場合(例えば90％以下)、位置ずれが生じたことにする。なお、被検体領域が写っているかどうかについては、被検体領域が取りうる画素値の閾値以上かどうかで判別する。
【００５４】
以上のようにして、位置決めに使用されたスカウト画像の撮像方向がAxial方向もしくはCoronal方向であり、本撮像の方向がSagittal方向である場合において、位置ずれを検知することができる。ただし、スカウト画像の方向及び本撮像の方向は、スカウト画像がAxial方向とCoronal方向で、本撮像の方向がSagittal方向である場合のみならず、その他の組み合わせであっても良い。例えば、位置決めに使用されたスカウト画像がCoronal方向もしくはSagittal方向で、本撮像の方向がAxial方向であっても良く、位置決めに使用されたスカウト画像がSagittal方向もしくはAxial方向で、本撮像の方向がCoronal方向であっても良い。
【００５５】
次に、本実施例(本撮像の方向がSagittal方向)における位置ずれ検知後の再撮像の方法についてその概略を説明する。
【００５６】
位置決めに使用されたスカウト画像の撮像方向がAxial方向もしくはCoronal方向であり、本撮像の方向がSagittal方向である場合に、位置ずれを補正して再撮像する手法について説明する。この方法は、被検体が直交座標系の3軸のいずれかを中心として回転した場合の位置ずれを補正して再撮像するものである。ここでは、位置ずれ検知後、ユーザより再撮像の要求があった場合、位置ずれの方向および位置ずれ量を割り出し、その位置ずれ量を本撮像位置を示す直線Lに反映させて再撮像を行う。なお、位置ずれの方向は、被検体領域幅L1a,L1a’およびL2a,L2a’とを比較した結果より、Sagittal方向の本撮像画像28aにおいて、被検体領域の縦幅と横幅のどちら方向に長短の差が生じているか、また、その長短差が、各位置ずれのケース毎に予め定めた位置ずれ量の範囲内にあるかどうかで判別する。
【００５７】
次に、位置ずれ方向判別後の再撮像方法について、位置ずれの方向毎に図3〜5を用いて説明する。
【００５８】
先ず、図3は、位置決めに使用されたスカウト画像の撮像方向がAxial方向もしくはCoronal方向であり、本撮像の方向がSagittal方向である場合において、スカウト撮像から本撮像までに被検体領域がZ軸を中心に回転した場合の再撮像を示す図である。ただし、図3において、30bは、Sagittal方向の本撮像画像28bより作成した基準周波数エンコードプロファイル、32aは、本撮像前に撮像されたCoronal方向のスカウト画像を回転させた位置ずれスカウト画像である。この方向の位置ずれにおいては、Coronal方向の位置ずれスカウト画像32aを用いて、位置ずれ量を確認し、本撮像位置を補正して再撮像を行う。以下、位置ずれ量の確認から再撮像までの詳細について、説明する。
【００５９】
先ず、Coronal方向の位置ずれスカウト画像32aでの本撮像位置(直線Lの位置)の被検体領域の幅を算出する。より具体的には、直線Lが横切る画素の画素値を縦軸とするプロファイルを生成する。続いて、そのプロファイルの裾野を所定の画素値の閾値(例えば被検体領域の画素値の1/20)に基づいて削除し、最初に被検体領域が確認できる点を点C2、点D2と定義し、この線分C2D2の長さを被検体領域幅L2bとして長さを求める。
【００６０】
次に、基準周波数エンコードプロファイル30bより、Sagittal方向の本撮像画像28bの被検体領域の縦幅を算出する。より具体的には、基準周波数エンコードプロファイル30bにおいて、Sagittal方向の本撮像画像28bの行位置を示す縦軸の両端からプロファイルをたどり、最初に被検体領域が確認できる点(あらかじめ定めた閾値以上の度数が確認できる点)を点C2’、点D2’とし、この線分C2’D2’の幅を被検体領域幅L2b’とする。ここで、算出した被検体領域幅L2b、L2b’を比較し、Coronal方向の位置ずれスカウト画像32aの中で、その差が最小となるものの回転角度を位置ずれ量とし、その位置ずれ量を本撮像位置を示す直線Lに反映させた上で、再撮像を行う。
【００６１】
次に、図4は、位置決めに使用されたスカウト画像の撮像方向がAxial方向もしくはCoronal方向であり、本撮像の方向がSagittal方向である場合において、スカウト撮像から本撮像までに被検体領域がX軸を中心に回転した場合の再撮像を示す図である。ただし、図4において31bは31aの範囲内を1としたマスク画像であり、32bはこのマスク画像31bを回転させた位置ずれマスク画像である。図4からは、31aの領域外に被検体領域がはみ出ている事がわかり、31aの代わりに32bの位置ずれマスク画像を用い、被検体領域が覆い隠されるかどうか確認することで、位置ずれ量がわかる。よって、この方向の位置ずれにおいては、位置ずれマスク画像32bを用いて、先と同様に、このマスク画像の1の領域内に、本撮像により描出された被検体領域の含まれている割合を確認し、位置ずれマスク画像32bの中で最も割合が高くなるものの回転角度を位置ずれ量とし、その位置ずれ量を本撮像位置を示す直線Lに反映させた上で、再撮像を行う。
【００６２】
続いて、図5は、位置決めに使用されたスカウト画像の撮像方向がAxial方向もしくはCoronal方向であり、本撮像の方向がSagittal方向である場合において、スカウト撮像から本撮像までに被検体領域がY軸を中心に回転した場合の再撮像を示す図である。ただし、図5において、29bは、Sagittal方向の本撮像画像28dより作成した基準位相エンコードプロファイル、32cは、本撮像前に撮像されたAxial方向のスカウト画像を回転させた位置ずれスカウト画像である。
【００６３】
この方向の位置ずれにおいては、Axial方向の位置ずれスカウト画像32cを用いて、位置ずれ量を確認し、本撮像位置を補正して再撮像を行う。以下、位置ずれ量の確認から再撮像までの詳細について、説明する。先ず、Coronal方向の位置ずれスカウト画像32cでの本撮像位置(直線Lの位置)の被検体領域の幅を算出する。
【００６４】
より具体的には、直線Lが横切る画素の画素値を縦軸とするプロファイルを生成する。続いて、そのプロファイルの裾野を所定の画素値の閾値(例えば被検体領域の画素値の1/20)に基づいて削除し、最初に被検体領域が確認できる点を点A2、点B2と定義し、この線分A2B2の長さを被検体領域幅L1bとして長さを求める。
【００６５】
次に、基準位相エンコードプロファイル29bより、Sagittal方向の本撮像画像28dの被検体領域の横幅を算出する。より具体的には、基準位相エンコードプロファイル29bにおいて、Sagittal方向の本撮像画像28dの列位置を示す横軸の両端からプロファイルをたどり、最初に被検体領域が確認できる点(あらかじめ定めた閾値以上の度数が確認できる点)を点A2’、点B2’とし、この線分A2’B2’の幅を被検体領域幅L1b’とする。ここで、算出した被検体領域幅L1b、L1b’を比較し、Axial方向の位置ずれスカウト画像32cの中で、その差が最小となるものの回転角度を位置ずれ量とし、その位置ずれ量を本撮像位置を示す直線Lに反映させた上で、再撮像を行う。
【００６６】
次に、図6(a)〜(c)を用い、本実施例の詳細を表すブロック図を説明する。
【００６７】
図6(a)において、22aは、RAM22上に展開され、CPU8に接続され、MRI装置の撮像を制御する撮像制御部である。22bは、撮像制御部22aに接続され、被検体の位置ずれを検出するための処理を行う位置ずれ検出部である。22cは、撮像制御部22aに接続され、被検体の位置ずれ量を算出し、被検体の動きに合わせて再撮像位置を設定するための再撮像位置設定部である。更に、18aは、磁気ディスク等から成る記憶部、20aは、ディスプレイ等から成る表示部、25aは、操作部である。
【００６８】
更に図6(b)には、位置ずれ検出部22bの内部の構成が示されている。22b1は、被検体のスカウト画像や、本撮像により得られた本撮像画像から所望の直線上の画素値の変化を表すプロファイルを生成するプロファイル生成部である。22b2は、プロファイル生成部により生成されたプロファイルより、被検体の領域の幅を算出する被検体領域幅算出部である。22b3は、スカウト画像より得られた被検体領域の幅と、本撮像画像より得られた被検体領域の幅を比較して、被検体の位置ずれを検出する幅比較部である。22b4は、スカウト画像より求められた、被検体が動いていない場合における被検体領域の幅(被検体領域幅L1a及び被検体領域幅L2a(破線により規定される幅))により、Sagittal方向の本撮像画像28a上に規定される領域(太枠31aで示される領域)内に、実際にどの程度の被検体領域が含まれるかに基づいて、被検体の位置ずれを検出する割合判定部である。
【００６９】
更に、図6(c)は、再撮像位置設定部22cの中身が示されている。22c1は、回転させた位置ずれマスク画像を作成する位置ずれマスク作成部である。
22c2は、位置ずれマスクを用いて、このマスクの中に本撮像により描出された被検体領域の含まれている割合を取得する割合算出部である。22c3は、割合算出部22c2で得られた割合が最も高くなるものの回転角度を位置ずれ量として求める位置ずれ量算出部である。22c4は、位置ずれ量を本撮像位置を表す直線Lに反映させるよう座標変換を行なう座標変換部である。22c5は、回転させたスカウト画像を生成するための回転スカウト画像生成部である。22c6は、回転スカウト画像等より、プロファイルを生成するプロファイル生成部である。22c7は、回転スカウト画像や、本撮像画像より得られたプロファイルを用いて、被検体領域の幅を算出する被検体領域幅算出部である。22c8は、回転スカウト画像より取得したプロファイルより得られる被検体領域幅と、本撮像画像より取得したプロファイルより得られる被検体領域幅を比較する領域幅比較部である。
【００７０】
次に、本実施例(本撮像の方向がSagittal方向)における位置ずれ検知および再撮像の動作について、図7を用いて説明する。
【００７１】
ステップS1：RAM22上に展開される撮像制御部22aは、スカウト画像(Axial,Coronal,Sagittal)を取得し、得られたスカウト画像を記憶部18aに記憶する。
【００７２】
ステップS2：表示部20aに、ステップS1で取得したAxial方向もしくはCoronal方向のスカウト画像を表示させ、表示されたスカウト画像上に、操作者は、本撮像位置を設定する(直線Lをセット)。
【００７３】
ステップS3：プロファイル生成部22b1は、本撮像方向(Sagittal方向)と異なるスカウト画像上(2枚)で、直線Lが横切る画素の画素値を縦軸とするプロファイルを生成する。
【００７４】
ステップS4：被検体領域幅算出部22b2は、ステップS3で生成したプロファイルの裾野を所定の画素値の閾値に基づいて削除し、本撮像位置における被検体領域の幅(L1a,L2a)を算出する。
【００７５】
ステップS5：撮像制御部22aは、本撮像を実施し、本撮像画像を取得する。
【００７６】
ステップS6：プロファイル生成部22b1は、ステップS5で取得した本撮像画像をフーリエ変換し、基準位相エンコードプロファイル29a、基準周波数エンコードプロファイル30aを生成する。
【００７７】
ステップS7：被検体領域幅算出部22b2は、ステップS6で生成したプロファイルより、本撮像画像における被検体領域幅(L1a’,L2a’)を算出する。
【００７８】
ステップS8：幅比較部22b3は、L1aとL1a’、もしくはL2aとL2a’を比較し、その差が既定の閾値以上の差がある場合、ステップS11に進み、無い場合にはステップS9に進む。
【００７９】
ステップS9：割合判定部22b4は、L1aおよびL2aより割り出した31aの領域内に被検体領域(被検体が取り得る画素値の閾値以上の箇所)の含まれる割合が既定の閾値以下である場合、ステップS11に進み、既定の閾値以下でない場合にはステップS10に進む。
【００８０】
ステップS10：本撮像画像を表示する。
【００８１】
ステップS11：位置ずれを検知した本撮像画像枠を赤枠で強調表示し、再撮像するか確認のメッセージを表示して操作者に位置ずれが生じたことを報知する。
【００８２】
ステップS12：再撮像する場合はステップS13に進み、しない場合は終了する。
【００８３】
ステップS13：X軸を中心とした回転の場合(L1a＞L1a’かつL2a＜L2a’であり、既定の位置ずれ量の範囲内である)、ステップS14に進み、そうでない場合はステップS17に進む。
【００８４】
ステップS14：位置ずれマスク生成部22c1は、被検体領域の存在する範囲を示す31aの領域よりマスク画像31bを作成し、マスク画像31bを回転させて位置ずれマスク画像32bも作成する。
【００８５】
ステップS15：割合算出部22c2は、位置ずれマスク画像32bの各対象領域内に被検体領域(被検体が取り得る画素値の閾値以上の箇所)の含まれる割合を算出する。
【００８６】
ステップS16：位置ずれ量算出部22c3は、ステップS15で算出した位置ずれマスク画像32bの内、最も割合が高かったものの位置ずれ量(回転角度)を補正値として設定する。
【００８７】
ステップS17：Y軸を中心とした回転の場合(L1a≠L1a’かつL2a＜L2a’であり、既定の位置ずれ量の範囲内である)、ステップS18に進み、そうでない場合(Z軸を中心とした回転)、ステップS21に進む。
【００８８】
ステップS18：回転スカウト画像生成部22c5は、Axial方向のスカウト画像26dを回転させた位置ずれスカウト画像32cを作成する。
【００８９】
ステップS19：プロファイル生成部22c6は、Axial方向の位置ずれスカウト画像32c上で、直線Lが横切る画素の画素値を縦軸とするプロファイルを作成する。
【００９０】
ステップS20：被検体領域幅算出部22c7は、ステップS19で生成したプロファイルの裾野を所定の画素値の閾値に基づいて削除し、本撮像位置における被検体領域幅L1bを算出する。
【００９１】
ステップS21：領域幅比較部22c8は、本撮像位置の被検体領域幅(L1a’、L1b)を比較する。
【００９２】
ステップS22：位置ずれ量算出部22c3は、ステップS21での比較結果の内、Axial方向の位置ずれスカウト画像32cの中で、その差が最小となるものの位置ずれ量(回転角度)を補正値として設定する。
【００９３】
ステップS23：回転スカウト画像生成部22c5は、Coronal方向のスカウト画像27bを回転させた位置ずれスカウト画像32aを作成する。
【００９４】
ステップS24：プロファイル生成部22c6は、Coronal方向の位置ずれスカウト画像32a上で、直線Lが横切る画素の画素値を縦軸とするプロファイルを作成する。
【００９５】
ステップS25：被検体領域幅算出部22c7は、ステップS24で生成したプロファイルの裾野を所定の画素値の閾値に基づいて削除し、本撮像位置における被検体領域幅L2bを算出する。
【００９６】
ステップS26：領域幅比較部22c8は、本撮像位置(直線Lの位置)の被検体領域幅(L2a’、L2b)を比較する。
【００９７】
ステップS27：位置ずれ量算出部22c3は、ステップS26での比較結果の内、Coronal方向の位置ずれスカウト画像32aの中で、その差が最小となるものの位置ずれ量(回転角度)を補正値として設定する。
【００９８】
ステップS28：座標変換部22c4は、本撮像位置を補正値分ずらして座標変換して再設定し、撮像制御部22aは、再撮像を実施する。
【００９９】
上記本発明の実施例1によれば、スカウト画像が、直交する3断面のいずれかであり、本撮像を行う方向がスカウト画像を撮像した方向のいずれかと平行である場合において、スカウト画像を撮像してから本撮像を撮像するまでに、被検体の向きが第1の角度から第2の角度へ動いたことを検出できる。
【実施例２】
【０１００】
次に、実施例2に係るMRI装置を図8を用い説明する。
【０１０１】
ここでは、本実施例(位置決めに使用されたスカウト画像の撮像方向がCoronal方向であり、本撮像の方向が斜断面方向)における位置ずれの検出についてその概略を説明する。
【０１０２】
先ず、図8は、位置決めに使用されたスカウト画像の撮像方向がCoronal方向であり、本撮像の方向が斜断面方向である場合において、スカウト撮像から本撮像までに被検体1が動かない場合の例を示す図である。ただし、図8において27eは、Coronal方向のスカウト画像であり、この画像上に、本撮像位置を示す直線L’がGUIにより描かれる。また、Coronal方向のスカウト画像27eには、画像上中央の縦方向に、中央線(Center1)が描かれている。このCenter1は、Sagittal方向のスカウト画像33aが取得される位置を示している。なお、Coronal方向のスカウト画像27eにおいて、本撮像位置を示す直線L’と被検体領域の外縁(輪郭)が交わる点の一方を点Eとし、直線L’とCenter1との交点を点Fとし、線分EFの幅をWidth1とする。
【０１０３】
続いて、図8において27fは、Coronal方向のスカウト画像27eを90°右に回転したものである。90°右回転したCoronal方向のスカウト画像27fにも本撮像位置を示す直線L’が描かれており、画像上中央の縦方向に中央線(Center2)が描かれている。このCenter2は、Axial方向のスカウト画像26eが取得される位置を示している。なお、90°右回転したCoronal方向のスカウト画像27fにおいて、本撮像位置を示す直線L’とCenter2との交点を点Gとし、線分EGの幅をWidth2とする。
【０１０４】
また、図8において、33aはSagittal方向のスカウト画像であり、Coronal方向のスカウト画像27eと被検体領域の縦方向の位置が一致するように配置されている。ここで、Coronal方向のスカウト画像27eにおける点FよりSagittal方向のスカウト画像33a上へ水平に破線を描き、この破線が被検体領域上を横切る長さを被検体領域幅L3とする。なお、被検体領域幅L3を求める際には、破線が横切る画素の画素値を縦軸とするプロファイルを生成し、そのプロファイルの裾野を所定の画素値の閾値(例えば被検体領域の画素値の1/20)に基づいて削除し、最初に被検体領域が確認できる点を点H、点Iと定義し、この線分HIの長さを求める。
【０１０５】
続いて、図8において、26eはAxial方向のスカウト画像であり、90°右回転したCoronal方向のスカウト画像27fと被検体領域の縦方向の位置が一致するように配置されている。ここで、90°右回転したCoronal方向のスカウト画像27fにおける点GよりAxial方向のスカウト画像26e上へ水平に破線を描き、この破線が被検体領域上を横切る長さを被検体領域幅L4とする。なお、被検体領域幅L4を求める際には、破線が横切る画素の画素値を縦軸とするプロファイルを生成し、そのプロファイルの裾野を所定の画素値の閾値(例えば被検体領域の画素値の1/20)に基づいて削除し、最初に被検体領域が確認できる点を点J、点Kと定義し、この線分JKの長さを求める。
【０１０６】
更に、図8において、34aは、本撮像により得られた斜断面方向の画像であり、被検体領域が斜線で示されている。ここで、基準周波数エンコードプロファイル30cは、斜断面方向の本撮像画像34aの画像データを、2次元フーリエ変換して得られるデータから、列方向の中央に配置された1列を抽出し、それを1次元フーリエ変換したものである。この基準周波数エンコードプロファイル30cを用い、斜断面方向の本撮像画像34aにおける被検体領域での最上部の行位置をMラインとする。なお、Mラインの行位置は、基準周波数エンコードプロファイル30cの周波数方向の端からプロファイルをたどり、最初に被検体領域が確認できる点(あらかじめ定めた閾値以上の度数が確認できる点)の一方とする。
【０１０７】
これによれば、被検体1がスカウト撮像から本撮像までに動かない場合、斜断面方向の本撮像画像34aのMラインからWidth1の長さぶん下にずれた位置の被検体領域の横幅L3aはL3に等しく、斜断面方向の本撮像画像34aのMラインからWidth2の長さぶん下にずれた位置の被検体領域の横幅L4aはL4に等しいことがわかる。
【０１０８】
次に、位置決めに使用されたスカウト画像の撮像方向がCoronal方向であり、本撮像の方向が斜断面方向である場合について、スカウト撮像から本撮像までに被検体1が動いた場合、本撮像画像上の被検体領域がどのように変化するかを説明する。
【０１０９】
先ず図9を用い、スカウト撮像から本撮像までに被検体1の頭部が左に傾いた場合(以下、Z軸を中心に回転した場合と表現する)、本撮像画像上の被検体領域がどのように変化するかを説明する。ただし、図9において、27e、27f、33b、26f、34bは、図8と同様に、Coronal方向のスカウト画像、90°右回転したCoronal方向のスカウト画像、Sagittal方向のスカウト画像、Axial方向のスカウト画像、本撮像により得られた画像である。図9のL3及びL4は、被検体1がスカウト撮像から本撮像までに動かなかった場合において、斜断面の本撮像画像34bのMラインからWidth1,Width2の長さぶん下にずれた位置の被検体領域の横幅であり、26fおよび33b内には、この場合の被検体1の頭部を表す模式図が薄線で描かれている。
【０１１０】
また、33b、26f、34b上では、被検体1の頭部がZ軸を中心に回転した場合、被検体1の頭部がどのように描出されるかが被検体1の頭部を表す模式図あるいは斜線で示されている。これによれば、被検体1の頭部がZ軸を中心に回転した場合、本撮像位置によって、被検体領域の横幅(L3b及びL4b)は、破線で示した動かなかった場合の被検体領域の横幅(L3及びL4)と比較して、狭くなっていることがわかる(L3＞L3b,L4＞L4b)。
【０１１１】
次に図10を用い、スカウト撮像から本撮像までに被検体1の頭部が顎を引いた場合(以下、X軸を中心に回転した場合と表現する)、本撮像画像上の被検体領域がどのように変化するかを説明する。ただし、図10において、27e、27f、33c、26g、34cは、図8と同様に、Coronal方向のスカウト画像、90°右回転したCoronal方向のスカウト画像、Sagittal方向のスカウト画像、Axial方向のスカウト画像、本撮像により得られた画像である。図10のL3及びL4は、被検体1がスカウト撮像から本撮像までに動かなかった場合において、斜断面の本撮像画像34cのMラインからWidth1の長さぶん下にずれた位置の被検体領域の横幅であり、26gおよび33c内には、この場合の被検体1の頭部を表す模式図が薄線で描かれている。
【０１１２】
また、33c、26g、34c上では、被検体1の頭部がX軸を中心に回転した場合、被検体1の頭部がどのように描出されるかが被検体1の頭部を表す模式図あるいは斜線で示されている。これによれば、被検体1の頭部がX軸を中心に回転した場合、本撮像位置によって、被検体領域の横幅L3cには、破線で示した動かなかった場合の被検体領域の横幅L3と比較して、長短の差ができ(L3≠L3c)、被検体領域の横幅L4cは、破線で示した動かなかった場合の被検体領域の横幅L4と比較して、長くなることがわかる(L4＜L4c)。
【０１１３】
次に図11を用い、スカウト撮像から本撮像までに被検体1の頭部が左を向いた場合(以下、Y軸を中心に回転した場合と表現する)、本撮像画像上の被検体領域がどのように変化するかを説明する。ただし、図11において、27e、27f、33d、26h、34dは、図8と同様に、Coronal方向のスカウト画像、90°右回転したCoronal方向のスカウト画像、Sagittal方向のスカウト画像、Axial方向のスカウト画像、本撮像により得られた画像である。図11のL3及びL4は、被検体1がスカウト撮像から本撮像までに動かなかった場合において、斜断面の本撮像画像34dのMラインからWidth1の長さぶん下にずれた位置の被検体領域の横幅であり、26hおよび33d内には、この場合の被検体1の頭部を表す模式図が薄線で描かれている。また、33d、26h、34d上では、被検体1の頭部がY軸を中心に回転した場合、被検体1の頭部がどのように描出されるかが被検体1の頭部を表す模式図あるいは斜線で示されている。これによれば、被検体1の頭部がY軸を中心に回転した場合、本撮像位置によって、被検体領域の横幅L3dは、破線で示した動かなかった場合の被検体領域の横幅L3と比較して、長くなり(L3＜L3d)、被検体領域の横幅L4bには、破線で示した動かなかった場合の被検体領域の横幅L4と比較して、長短差が生じることがわかる(L4≠L4b)。
【０１１４】
以上より、位置決めに使用されたスカウト画像の撮像方向がCoronal方向であり、本撮像の方向が斜め方向である場合において、位置ずれの特徴を説明したので、この特徴を用いた本実施例における位置ずれの検知の手法について説明する。この方法は、上述した直交座標系の3軸を中心とした回転による位置ずれを検知する方法である。
【０１１５】
また、本撮像位置を設定したスカウト画像と異なる2方向から撮像されたスカウト画像における、本撮像位置の被検体領域幅が、本撮像画像上での被検体領域幅と一致するかにより、位置ずれを検知する。
【０１１６】
以下、位置決めに使用されたスカウト画像の撮像方向がCoronal方向であり、本撮像方向が斜断面方向である場合において、位置ずれの検知方法の詳細を図8を用い説明する。
先ず、Coronal方向のスカウト画像27eと、Sagittal方向のスカウト画像33aの被検体領域の縦方向の位置が一致するように配置し、更にCoronal方向のスカウト画像27e上にSagittal方向のスカウト画像33aが取得される位置を表す中央線(Center1)を描く。続いて、Coronal方向のスカウト画像27e上に描かれた本撮像位置を示す直線L’と中央線(Center1)との交点を点Fとする。
【０１１７】
そして、Coronal方向のスカウト画像27eにおける点FよりSagittal方向のスカウト画像33a上へ水平に破線を描き、この破線が被検体領域上を横切る長さをL3とする。なお、L3を求める際には、破線が横切る画素の画素値を縦軸とするプロファイルを生成し、そのプロファイルの裾野を所定の画素値の閾値(例えば被検体領域の画素値の1/20)に基づいて削除し、それにより被検体領域の破線上の長さを求める。
【０１１８】
また、Coronal方向のスカウト画像27e上で本撮像位置を示す直線L’と被検体領域の外縁(輪郭)が交わる点の一方を点Eとし、線分EFの幅をWidth1とする。なお、点Eは、直線L’の端から確認していき、あらかじめ被検体領域として定めた閾値以上の画素値を最初に確認した点とする。
【０１１９】
次に、90°右回転したCoronal方向のスカウト画像27fと、Axial方向のスカウト画像26eの被検体領域の縦方向の位置が一致するように配置し、更に90°右回転したCoronal方向のスカウト画像27f上にAxial方向のスカウト画像26eの取得される位置を表す中央線(Center2)を描く。続いて、90°右回転したCoronal方向のスカウト画像27f上に描かれた本撮像位置を示す直線L’と中央線(Center2)との交点を点Gとする。そして、90°右回転したCoronal方向のスカウト画像27fにおける点GよりAxial方向のスカウト画像26e上へ水平に破線を描き、この破線が被検体領域上を横切る長さをL4とする。なお、L4を求める際には、破線が横切る画素の画素値を縦軸とするプロファイルを生成し、そのプロファイルの裾野を所定の画素値の閾値(例えば被検体領域の画素値の1/20)に基づいて削除し、それにより被検体領域の破線上の長さを求める。また、90°右回転したCoronal方向のスカウト画像27f上の線分EGの幅をWidth2とする。
【０１２０】
次に、本撮像により得られた画像34aについて、画像データを2次元フーリエ変換して得られるデータから、列方向の中央に配置された1列を抽出し、それを1次元フーリエ変換して、基準周波数エンコードプロファイル30cを作成する。この基準周波数エンコードプロファイル30cを用い、本撮像により得られた画像34aにおける被検体領域での最上部の行位置をMラインとして定義する。なお、このプロファイル生成については、画像再構成前に行っておいても良い。
【０１２１】
更に、Mラインから上述したWidth1,Width2の長さぶん下にずれた位置の被検体領域の横幅L3a,L4aをL3,L4と同様の求め方で算出する。
【０１２２】
ここで、算出した被検体領域幅L3a,L4aを先に算出したL3,L4と比較し、予め定めた閾値以上の差があれば位置ずれが生じたこととする。次に、本実施例(位置決めに使用されたスカウト画像の撮像方向がCoronal方向であり、本撮像の方向が斜断面方向)における位置ずれ検知後の再撮像の方法についてその概略を説明する。この方法は、直交座標系を中心とした回転による位置ずれを補正して再撮像するものであり、ここでは、位置ずれ検知後、ユーザより再撮像の要求があった場合、スカウト画像を再撮像し、元のスカウト画像と比較することで位置ずれ量(回転角)を割り出し、その位置ずれ量を本撮像位置を示す直線L’に反映させて再撮像を行う。
【０１２３】
図12は、本撮像方向が斜断面である場合において、スカウト撮像から本撮像までに被検体領域が直交座標系の3軸を中心に回転した場合の再撮像を示す図である。ただし、図12において27gは、スカウト再撮像前に取得されたCoronal方向のスカウト画像、26iは、スカウト再撮像前に取得されたAxial方向のスカウト画像、33eは、スカウト再撮像前に取得されたSagittal方向のスカウト画像である。
【０１２４】
また、32dは、27gを回転させた位置ずれスカウト画像、32eは、26iを回転させた位置ずれスカウト画像、32fは、33eを回転させた位置ずれスカウト画像である。更に、27hは、再撮像したCoronal方向のスカウト画像であり、26jは、再撮像したAxial方向のスカウト画像であり、33fは、再撮像したSagittal方向のスカウト画像である。この方向の位置ずれにおいては、スカウト画像と本撮像画像のみから位置ずれ量を見積もることは困難な為、再撮像したスカウト画像と、元のスカウト画像を回転させた位置ずれスカウト画像とをテンプレートマッチングすることで位置ずれ量(回転角)を算出し、本撮像位置を示す直線L’を位置ずれ量(回転角)ぶん傾けた上で、再撮像を行う。
【０１２５】
次に、図13(a)〜(c)を用い、本実施例の詳細を表すブロック図を説明する。
【０１２６】
図13(a)において、22aは、RAM22上に展開され、CPU8に接続され、MRI装置の撮像を制御する撮像制御部である。22bは、撮像制御部22aに接続され、被検体の動いたかを検出するための処理を行う位置ずれ検出部である。22cは、撮像制御部22aに接続され、被検体の位置ずれ量を算出し、被検体の動きに合わせて再撮像位置を設定するための再撮像位置設定部である。更に、18aは、磁気ディスク等から成る記憶部、20aは、ディスプレイ等から成る表示部、25aは、操作部等である。
【０１２７】
更に図13(b)には、位置ずれ検出部22bの中身が示されている。22b5は、本実施例における座標調整を行ったりする座標調整部である。22b6は、被検体のスカウト画像や、本撮像により得られた本撮像画像から所望の直線上の画素値の変化を表すプロファイルを生成するプロファイル生成部である。22b7は、スカウト画像や、本撮像により得られる本撮像画像における被検体領域の幅を算出するための被検体領域幅算出部である。22b8は、スカウト画像により得られた幅と本撮像により得られた幅を比較するための幅比較部である。更に、図13(c)は、再撮像位置設定部22cの中身が示されている。22c9は、位置ずれスカウト像を生成する位置ずれスカウト像生成部である。22c10は、プロファイルを生成するプロファイル生成部である。22c11は、相関係数を算出する相関係数算出部である。22c12は、位置ずれ量を算出する位置ずれ量算出部である。22c13は、位置ずれ量を本撮像位置を表す直線Lに反映させるよう座標変換を行なう座標変換部である。
【０１２８】
次に、本実施例(位置決めに使用されたスカウト画像の撮像方向がCoronal方向であり、本撮像の方向が斜断面方向)における位置ずれ検知および再撮像の動作について、図14を用いて説明する。
【０１２９】
ステップT1：RAM22上に展開される撮像制御部22aは、スカウト画像(Axial,Coronal,Sagittal)を取得し、得られたスカウト画像を記憶部18aに記憶する
ステップT2：表示部20aに、ステップT1で取得したCoronal方向のスカウト画像を表示させ、表示されたスカウト画像上に、操作者は、本撮像位置を設定する(直線L’をセット)。
【０１３０】
ステップT3：座標調整部22b5は、本撮像位置が設定されたCoronal方向のスカウト画像27eとSaggital方向のスカウト画像33aで、被検体領域の縦方向の位置が一致するように配置する。
【０１３１】
ステップT4：座標調整部22b5は、Coronal方向のスカウト画像27e上に、中央線(Center1)を描き、直線L’と中央線(Center1)との交点Fを求める。
【０１３２】
ステップT5：プロファイル生成部22b6は、Saggital方向のスカウト画像33a上で点Fの水平ラインが横切る画素の画素値を縦軸とするプロファイルを生成する。
【０１３３】
ステップT6：被検体領域幅算出部22b7は、ステップT5で生成したプロファイルの裾野を所定の画素値の閾値に基づいて削除し、本撮像位置における被検体領域幅L3を算出する。
【０１３４】
ステップT7：座標調整部22b5は、直線L’と被検体領域の外縁(輪郭)が交わる点Eを求め、線分EFの幅(Width1)を算出する。
【０１３５】
ステップT8：座標調整部22b5は、90°右回転したCoronal方向のスカウト画像27fと、Axial方向のスカウト画像26eで、被検体領域の縦方向の位置が一致するように配置する。
【０１３６】
ステップT9：座標調整部22b5は、90°右回転したCoronal方向のスカウト画像27f上に中央線(Center2)を描き、直線L’と中央線(Center2)との交点Gを求める。
【０１３７】
ステップT10：プロファイル生成部22b6は、Axial方向のスカウト画像26e上で点Gの水平ラインが横切る画素の画素値を縦軸とするプロファイルを生成する。
【０１３８】
ステップT11：被検体領域幅算出部22b7は、ステップT10で生成したプロファイルの裾野を所定の画素値の閾値に基づいて削除し、本撮像位置における被検体領域幅L4を算出する。
【０１３９】
ステップT12：座標調整部22b5は、線分EGの幅(Width2)を算出する。
【０１４０】
ステップT13：撮像制御部22aは、本撮像を実施し、本撮像画像34aを取得する。
【０１４１】
ステップT14：座標調整部22b5は、本撮像画像34aをフーリエ変換し、基準周波数エンコードプロファイル30cを生成する。
【０１４２】
ステップT15：座標調整部22b5は、ステップT14で生成したプロファイルより、本撮像画像34aにおける被検体領域での最上部の行位置(Mライン)を求める。
【０１４３】
ステップT16：座標調整部22b5は、本撮像画像34aにおいて、MラインのWidth1の長さぶん下の位置の被検体領域幅L3a、およびMラインのWidth2の長さぶん下の位置の被検体領域幅L4aを算出する。
【０１４４】
ステップT17：幅比較部22b8は、L3とL3aもしくはL4とL4aに既定の閾値以上の差がある場合、ステップT19に進み、既定の閾値以上の差がない場合にはステップT18に進む。
【０１４５】
ステップT18：撮像制御部22aは、本撮像画像を表示する。
【０１４６】
ステップT19：表示部20aは、位置ずれを検知した本撮像画像の枠を赤枠で強調表示し、再撮像するか確認のメッセージを表示し、操作者に位置ずれが生じたことを報知する。
【０１４７】
ステップT20：再撮像する場合はステップT21に進み、しない場合は終了する。
【０１４８】
ステップT21：撮像制御部22aは、スカウト画像(Axial,Coronal,Sagittal)を再取得する(27h,26j,33f)。
【０１４９】
ステップT22：位置ずれスカウト像生成部22c9は、ステップT1で取得したスカウト画像を回転させた位置ずれスカウト画像(32d,32e,32f)を作成する。
【０１５０】
ステップT23：プロファイル生成部22c10は、ステップT21とステップT22の撮像方向が等しいスカウト画像間において、x,y,z方向いずれかにおける全ラインのプロファイルを取得する。
【０１５１】
ステップT24：相関係数算出部22c11は、ステップT21とステップT22の撮像方向が等しいスカウト画像間において、全ラインでの相関係数を算出し、それらの平均値を算出する。
例えば、Coronal方向のスカウト画像を回転させた位置ずれスカウト画像32dと再撮像したCoronal方向のスカウト画像27hの場合、それぞれで得られたx軸方向のプロファイル信号値を、RoImg_xi、ReImg_xi(i=1,2,3,…n、n：ライン数)とし、各プロファイルの相加平均を
2014046123.tif 000003
として、2つのプロファイルの相関係数を以下の(1)式で算出する。
2014046123.tif 000004
ステップT25：位置ずれ量算出部22c12は、各スカウト画像(Axial,Coronal,Sagittal)において、最も相関係数(平均値)の高い位置ずれスカウト画像の位置ずれ量(回転角度)を補正値として設定する。
【０１５２】
ステップT26：座標変換部22c13は、本撮像位置を補正値分ずらして座標変換して再設定し、撮像制御部22aは、再撮像を実施する。
【０１５３】
上記本発明の実施例2によれば、スカウト画像が、直交する3断面のいずれかであり、本撮像を行う方向が斜断面である場合において、スカウト画像を撮像してから本撮像を撮像するまでに、被検体の向きが第1の角度から第2の角度へ動いたことを検知できる。
【産業上の利用可能性】
【０１５４】
本発明は、MRI装置に利用することができる。
【符号の説明】
【０１５５】
26a  Axial方向のスカウト画像(Saggital方向の本撮像前：位置ずれなし)、27a  Coronal方向のスカウト画像(Saggital方向の本撮像前：位置ずれなし)、28a  Sagittal方向の本撮像画像(位置ずれなし)、29a  基準位相エンコードプロファイル(Saggital方向の本撮像時：位置ずれなし)、30a  基準周波数エンコードプロファイル(Saggital方向の本撮像時：位置ずれなし)、31a  被検体領域の存在する範囲
【図１】
2014046123.tif 000005
【図２】
2014046123.tif 000006
【図３】
2014046123.tif 000007
【図４】
2014046123.tif 000008
【図５】
2014046123.tif 000009
【図６】
2014046123.tif 000010
【図７】
2014046123.tif 000011
【図８】
2014046123.tif 000012
【図９】
2014046123.tif 000013
【図１０】
2014046123.tif 000014
【図１１】
2014046123.tif 000015
【図１２】
2014046123.tif 000016
【図１３】
2014046123.tif 000017
【図１４】
2014046123.tif 000018