JP5177367B2 - 磁気共鳴イメージング装置及びナビゲータデータ解析方法 - Google Patents

磁気共鳴イメージング装置及びナビゲータデータ解析方法 Download PDF

Info

Publication number
JP5177367B2
JP5177367B2 JP2007221489A JP2007221489A JP5177367B2 JP 5177367 B2 JP5177367 B2 JP 5177367B2 JP 2007221489 A JP2007221489 A JP 2007221489A JP 2007221489 A JP2007221489 A JP 2007221489A JP 5177367 B2 JP5177367 B2 JP 5177367B2
Authority
JP
Japan
Prior art keywords
navigator
subject
phase
unit
tissue
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
JP2007221489A
Other languages
English (en)
Other versions
JP2009050544A (ja
Inventor
健一 神田
雄治 岩舘
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
GE Medical Systems Global Technology Co LLC
Original Assignee
GE Medical Systems Global Technology Co LLC
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by GE Medical Systems Global Technology Co LLC filed Critical GE Medical Systems Global Technology Co LLC
Priority to JP2007221489A priority Critical patent/JP5177367B2/ja
Priority to US12/200,818 priority patent/US7782053B2/en
Publication of JP2009050544A publication Critical patent/JP2009050544A/ja
Application granted granted Critical
Publication of JP5177367B2 publication Critical patent/JP5177367B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00—Arrangements or instruments for measuring magnetic variables
    • G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48—NMR imaging systems
    • G01R33/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
    • G01R33/56—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
    • G01R33/567—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution gated by physiological signals, i.e. synchronization of acquired MR data with periodical motion of an object of interest, e.g. monitoring or triggering system for cardiac or respiratory gating
    • G01R33/5676—Gating or triggering based on an MR signal, e.g. involving one or more navigator echoes for motion monitoring and correction
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00—Arrangements or instruments for measuring magnetic variables
    • G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48—NMR imaging systems
    • G01R33/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
    • G01R33/56—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
    • G01R33/565—Correction of image distortions, e.g. due to magnetic field inhomogeneities
    • G01R33/56509—Correction of image distortions, e.g. due to magnetic field inhomogeneities due to motion, displacement or flow, e.g. gradient moment nulling

Landscapes

  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Cardiology (AREA)
  • Radiology & Medical Imaging (AREA)
  • Physiology (AREA)
  • Power Engineering (AREA)
  • Pulmonology (AREA)
  • General Health & Medical Sciences (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Biophysics (AREA)
  • Signal Processing (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Condensed Matter Physics & Semiconductors (AREA)
  • General Physics & Mathematics (AREA)
  • Magnetic Resonance Imaging Apparatus (AREA)
  • Image Processing (AREA)
  • Image Analysis (AREA)

Description

本発明は、磁気共鳴イメージング装置（ＭＲＩ：Ｍａｇｎｅｔｉｃ Ｒｅｓｏｎａｎｃｅ Ｉｍａｇｉｎｇ）及びナビゲータデータ（ｎａｖｉｇａｔｏｒ ｄａｔａ）解析方法に関し、特に静磁場空間内において被検体にＲＦパルス（ｐｕｌｓｅ）を送信し、被検体から磁気共鳴信号を収集するスキャン（ｓｃａｎ）を実施することによって、被検体についての画像を生成する磁気共鳴イメージング装置、及びスキャンを実施することによって得られる、被検体における組織の位置を検出するためのナビゲータデータを解析するナビゲータデータ解析方法に関するものである。
磁気共鳴イメージング装置は、静磁場空間内の被検体に電磁波を照射することにより、その被検体内のプロトン（ｐｒｏｔｏｎ）のスピン（ｓｐｉｎ）を核磁気共鳴現象によって励起させ、その励起されたスピンにより発生する磁気共鳴信号を得るスキャンを実施する。そして、そのスキャンにより得られる磁気共鳴信号に基づいて、被検体の断層面についてスライス画像を生成する装置である。
この磁気共鳴イメージング装置を用いて被検体を撮影する際に被検体に体動が生じた場合には、生成したスライス画像に体動アーチファクト（ａｒｔｉｆａｃｔ）が発生する場合がある。たとえば、被検体の心臓や腹部を撮影する場合においては、呼吸運動や心拍運動などの体動によって、体動アーチファクトが発生し、画像の質が低下する。
そこで、この体動アーチファクトによる画像の低下の問題を解決する手法が提案されている。その一つの方法は、例えば、通常呼吸下での撮像において、横隔膜の位置の変化に合わせて、被検体の励起断面をリアルタイムで補正し、常に同一の断面から磁気共鳴信号を計測することにより体動アーチファクトによる画像の低下を防止する。また、取得したナビゲータエコーを利用して、イメージングシーケンスを変更したり、イメージングデータを選択したりすることにより体動アーチファクトによる画像の質の低下を防止する（例えば、特許文献１、２参照）。
これらの技術に用いる横隔膜の位置検出の手法として、ナビゲータエコーを利用して横隔膜の動きを追跡し、得たれたナビゲータデータを利用して呼吸同期、ゲーティング（ｇａｔｉｎｇ）を行う方法が知られている（例えば、特許文献３参照）。
そして、ナビゲータデータの解析方法として、例えばＡｈｎ法、ＬＳＱ法、エッジ（ｅｄｇｅ）検出法などが知られている（例えば、非特許文献１参照）。
特開２００７−１１１１８８号公報（例えば、段落００４７から段落００５５参照。） 特開２００７−９８０２６号公報（例えば、段落００５７から段落００６２参照。） 特開平１０−２７７０１０号公報（例えば、段落０００２から段落０００６参照。） Yiping P. Du, Manojkumar Saranathan, and Thomas K. F. Foo. "An Accurate, Robust, and Computationally Efficient Navigator Algorithm for Measuring Diaphragm Positions",JOURNAL OF CARDIOVASCULAR MAGNETIC RESONANCE Vol.6, No.2, pp. 483-490, 2004
しかし、図１８のコロナル（ｃｏｒｏｎａｌ）画像に示すように、イメージングスキャン（ｉｍａｇｉｎｇ ｓｃａｎ）を実施してイメージングデータ（ｉｍａｇｉｎｇ ｄａｔａ）を取得するイメージング領域ＩＡがナビゲータデータの取得位置であるナビゲータ領域ＮＡと重なった結果、取得したナビゲータデータにスライス干渉による信号の乱れが発生し、図１９の破線部に示すように、取得したナビゲータデータの信号強度Ｉとナビゲータ領域の位置Ｌとの関係をプロットした信号強度プロファイル（ｐｒｏｆｉｌｅ）にノイズ（ｎｏｉｓｅ）が発生する。ここで、図１９における破線部は、図１８におけるイメージング領域ＩＡとナビゲータ領域ＮＡとの重複部分における信号強度プロファイルを示す。そうすると、上記に示す従来のナビゲータデータの解析法で安定した解析結果を得ることは困難であった。
したがって、本発明は、ナビゲータデータの安定した解析結果を得ることにより、画像の低下を防止できる磁気イメージング装置、及びナビゲータデータの安定した解析結果を得ることができるナビゲータデータ解析方法を提供する。
さらに、本発明は、イメージングスキャンの撮影スライスとナビゲータ領域とが重なった場合であっても、ナビゲータデータの安定した解析結果を得ることにより、画像の低下を防止できる磁気イメージング装置、及びナビゲータデータの安定した解析結果を得ることができるナビゲータデータ解析方法を提供する。
本発明は、被検体において体動する組織を含むナビゲータ領域から磁気共鳴信号をナビゲータデータとして得るナビゲータシーケンス（ｎａｖｉｇａｔｏｒ ｓｅｑｕｅｎｃｅ）を実施すると共に、前記被検体においてイメージング領域から磁気共鳴信号をイメージングデータとして得るイメージングシーケンス（ｉｍａｇｉｎｇ ｓｅｑｕｅｎｃｅ）を実施するスキャンを実施することによって、前記イメージング領域について画像を生成する磁気共鳴イメージング装置であって、前記ナビゲータデータの位相と前記ナビゲータ領域の位置との関係を示すように位相プロファイルを生成する位相プロファイル生成部と、前記位相プロファイル生成部によって生成された位相プロファイルについて折り返りを補正する位相補正部と、前記位相補正部によって補正された位相プロファイルに基づいて、前記ナビゲータ領域において体動する組織の位置を検出する位置検出部とを有する。
好適には、前記位相補正部により位相の補正が実施された前記磁気共鳴信号の位相プロファイルにおいて、前記被検体における組織の位置を検出する範囲を設定する範囲設定部を有する。
好適には、前記位置検出部は、エッジ検出法により前記被検体における組織の位置を検出する。
好適には、前記位置検出部は、ＬＳＱ法により前記被検体における組織の位置を検出する。
好適には、前記位置検出部は、相関係数法により前記被検体における組織の位置を検出する。
好適には、前記位置検出部は、相互情報量法により前記被検体における組織の位置を検出する。
好適には、前記組織が横隔膜である。
好適には、前記位相補正部により位相の補正が実施された位相プロファイルを平滑化させる位相プロファイル平滑部を有する。
本発明のナビゲータデータ解析方法は、静磁場空間内において被検体にＲＦパルスを送信し、当該被検体から磁気共鳴信号を収集するスキャンを実施することにより、前記被検体のスライス画像を生成するためのイメージングデータと前記被検体における組織の位置を検出するためのナビゲータデータを取得し、取得した当該ナビゲータデータを解析するナビゲータデータ解析方法であって、前記ナビゲータデータの位相と前記ナビゲータ領域の位置との関係を示すように位相プロファイルを生成する第１ステップと、前記第１ステップにおいて生成された前記位相プロファイルの折り返りを補正する第２ステップと、前記第２ステップにおいて折り返りの補正がされた位相プロファイルに基づいて、前記ナビゲータ領域において体動する組織の位置を検出する第３ステップとを含む。
好適には、前記第２ステップの後、位相の補正が実施された前記位相プロファイルにおいて、前記被検体における組織の位置を検出する範囲を設定するステップを有する。
好適には、前記第３ステップは、エッジ検出法により前記被検体における組織の位置を検出する。
好適には、前記第３ステップは、ＬＳＱ法により前記被検体における組織の位置を検出する。
好適には、前記第３ステップは、相関係数法により前記被検体における組織の位置を検出する。
好適には、前記第３ステップは、相互情報量法により前記被検体における組織の位置を検出する。
好適には、前記組織が横隔膜である。
好適には、前記第２ステップの後、前記第２ステップにおいて位相の補正が実施された前記位相プロファイルを平滑化するステップを有する。
本発明によれば、イメージングスキャンの撮影スライスとナビゲータ領域とが重なった場合であっても、ナビゲータデータの安定した解析結果を得ることができる磁気イメージング装置及びナビゲータデータ解析方法を提供することができる。
以下より、本発明に係る一実施形態について図面を参照して説明する。
＜第１の実施形態＞
（装置構成）
図１は、本発明に係る一実施形態におけるＲＦコイル部により構成される磁気共鳴イメージング装置の構成を示す構成図である。本装置は、本発明の実施形態の一例である。
図１に示すように、磁気共鳴イメージング装置１は、スキャン部２と、操作コンソール（ｃｏｎｓｏｌｅ）部３とを有している。ここで、スキャン部２は、静磁場マグネット（ｍａｇｎｅｔ）部１２と、勾配コイル（ｃｏｉｌ）部１３と、ＲＦコイル部１４と、クレードル（ｃｒａｄｌｅ）１５とを有する。そして、操作コンソール部３は、ＲＦ駆動部２２と、勾配駆動部２３と、データ収集部２４と、制御部３０と、記憶部３１と、操作部３２と、データ処理部３３と、表示部３４とを有する。
スキャン部２について説明する。
スキャン部２は、図１に示すように、被検体４０における撮影スライス領域が収容される静磁場空間１１を含んでいる。そして、スキャン部２は、操作コンソール部３からの制御信号に基づいて、その静磁場が形成される静磁場空間１１に収容した被検体４０の撮影領域にＲＦパルスを照射し、その撮影領域から磁気共鳴信号を取得するスキャンを実施する。
本実施形態において、スキャン部２は、被検体４０のイメージング領域ＩＡにおいて発生する磁気共鳴信号をイメージングデータとして得るイメージングシーケンスＩＳと、被検体４０のナビゲータ領域ＮＡにおいて発生する磁気共鳴信号をナビゲータデータとして取得するナビゲータシーケンスＮＳとを繰り返し実施する。
スキャン部２の各構成要素について、順次、説明する。
静磁場マグネット部１２は、被検体４０が収容される静磁場空間１１に静磁場を形成するために設けられている。静磁場マグネット部１２は、水平磁場型であって、被検体４０が収容される静磁場空間１１において載置される被検体４０の体軸方向（ｚ方向）に沿うように、超伝導磁石（図示なし）が静磁場を形成する。なお、静磁場マグネット部１２は、水平磁場型の他に、垂直磁場型であってもよく、永久磁石により構成されていてもよい。
勾配コイル部１３は、ＲＦコイル部１４が受信する磁気共鳴信号に３次元の位置情報を持たせるために、静磁場空間１１に勾配磁場を形成する。勾配コイル部１３は、スライス選択勾配磁場、読み取り勾配磁場、位相エンコード勾配磁場の３種類の勾配磁場を形成するために勾配コイルを３系統有する。
ＲＦコイル部１４は、たとえば、被検体４０を囲むように配置される。ＲＦコイル部１４は、静磁場マグネット部１２によって静磁場が形成される静磁場空間１１内において、制御部３０からの制御信号に基づいて、電磁波であるＲＦパルスを被検体４０に送信して高周波磁場を形成する。これにより、被検体４０の撮影スライス領域におけるプロトンのスピンを励起する。そして、ＲＦコイル部１４は、その励起された被検体４０の撮影スライス領域におけるプロトンのスピンが元の磁化ベクトル（ｖｅｃｔｏｒ）へ戻る際に生ずる電磁波を、磁気共鳴信号として受信する。ＲＦコイル部１４は、同一のＲＦコイルによりＲＦパルスの送受信を行ってもよい。
本実施形態において、ＲＦコイル部１４は、ナビゲータシーケンスＮＳとイメージングシーケンスＩＳとにおいて、ＲＦパルスを送信する。
クレードル１５は、被検体４０を載置するテーブル（ｔａｂｌｅ）を有する。クレードル１５は、制御部３０からの制御信号に基づいて、静磁場空間１１の内部と外部との間でテーブルに載置された被検体４０を移動する。
操作コンソール部３について説明する。
操作コンソール部３は、スキャン部２が被検体４０についてスキャンを実施するように制御し、そのスキャン部２が実施したスキャンによって得られた磁気共鳴信号に基づいて、被検体４０の画像を生成すると共に、その生成した画像を表示する。
操作コンソール部３を構成する各部について、順次、説明する。
ＲＦ駆動部２２は、ＲＦコイル部１４を駆動させて静磁場空間１１内に高周波磁場を形成するために、ゲート（ｇａｔｅ）変調器（図示なし）とＲＦ電力増幅器（図示なし）とＲＦ発振器（図示なし）とを有する。ＲＦ駆動部２２は、制御部３０からの制御信号に基づいて、ＲＦ発振器からのＲＦ信号を、ゲート変調器を用いて所定のタイミング（ｔｉｍｉｎｇ）および所定の包絡線の信号に変調する。そして、ゲート変調器により変調されたＲＦ信号を、ＲＦ電力増幅器により増幅した後、ＲＦコイル部１４に出力する。
勾配駆動部２３は、制御部３０の制御信号に基づいて勾配コイル部１３を駆動させて、静磁場空間１１内に勾配磁場を発生させる。勾配駆動部２３は、勾配コイル部１３の３系統の勾配コイルに対応して３系統の駆動回路（図示なし）を有する。
データ収集部２４は、ＲＦコイル部１４が受信する磁気共鳴信号を収集するために、位相検波器（図示なし）とアナログ／デジタル変換器（図示なし）とを有する。データ収集部２４は、ＲＦコイル部１４からの磁気共鳴信号を、ＲＦ駆動部２２のＲＦ発振器の出力を参照信号として、位相検波器によって位相検波し、アナログ／デジタル変換器に出力する。そして、位相検波器により位相検波されたアナログ信号である磁気共鳴信号を、アナログ／デジタル変換器によってデジタル信号に変換して、データ処理部３３に出力する。
制御部３０は、コンピュータ（ｃｏｍｐｕｔｅｒ）と、コンピュータを用いて所定のスキャンに対応する動作を各部に実行させるプログラム（ｐｒｏｇｒａｍ）とを有する。そして、制御部３０は、後述する操作部３２に接続されており、操作部３２に入力された操作信号を処理し、クレードル１５とＲＦ駆動部２２と勾配駆動部２３とデータ収集部２４との各部に、制御信号を出力し制御を行う。また、制御部３０は、所望の画像を得るために、操作部３２からの操作信号に基づいてデータ処理部３３、表示部３４を制御する。
本実施形態においては、制御部３０は、ＲＦ駆動部２２と勾配駆動部２３とを制御して、スキャン部２にナビゲータシーケンスＮＳ及びイメージングシーケンスＩＳを実施させる。制御部３０は、ナビゲータシーケンスＮＳを実施して得られたナビゲータデータに基づいて、スキャン部２が実施するイメージングシーケンスＩＳの制御を実施し、またデータ処理部３３を制御することにより、イメージングデータの補正をする。
記憶部３１は、コンピュータとコンピュータに所定のデータ処理を実行させるプログラムとを有する。そして、記憶部３１は、データ収集部２４に収集された画像再構成処理前の磁気共鳴信号、後述するデータ処理部３３で画像再構成処理された画像データ等を記憶する。
操作部３２は、キーボード（ｋｅｙｂｏａｒｄ）やマウス（ｍｏｕｔｈ）などの操作デバイス（ｄｅｖｉｃｅ）により構成されている。操作部３２は、オペレータ（ｏｐｅｒａｔｏｒ）によって操作データ、撮像プロトコル（ｐｒｏｔｏｃｏｌ）などが入力され、またイメージングシーケンスＩＳを実施する領域、ナビゲータシーケンスＮＳを実施する領域の設定がされ、その操作データ、撮像プロトコル、設定領域に関するデータを制御部３０に出力する。
データ処理部３３は、コンピュータと、そのコンピュータを用いて所定のデータ処理を実行するプログラムを記憶するメモリとを有しており、制御部３０からの制御信号に基づいて、データ処理を実施する。
ここで、図２は、本発明に係る一実施形態におけるデータ処理部３３の構成を示すブロック（ｂｌｏｃｋ）図である。
図２に示すように、データ処理部３３は、画像再構成部３３１と、位相プロファイル生成部３３２と、位相補正部３３３と、位置検出部３３４とを有する。
データ処理部３３の各部について順次説明する。
画像再構成部３３１は、コンピュータとコンピュータに所定のデータ処理を実行させるプログラムを記憶するメモリとを有する。そして、画像再構成部３３１は、被検体４０の撮影領域についてのイメージングシーケンスＩＳの実施によってイメージングデータとして得られた磁気共鳴信号を、ローデータとして用いることによって、その被検体４０の撮影領域についての画像を画像再構成する。すなわち、本実施形態においては、スキャン部２によるイメージングシーケンスＩＳの実施において、ＲＦコイル部１４が撮影領域へＲＦパルスを送信し、その撮影領域において発生する磁気共鳴信号を、ＲＦコイル部１４が受信することによって収集したイメージングデータに基づいて、スライス画像を画像再構成する。そして、画像再構成されたスライス画像を表示部３４に出力する。
位相プロファイル生成部３３２は、例えば、被検体４０のイメージング領域ＩＡについてのイメージングシーケンスＩＳの実施前に実施されたナビゲータシーケンスＮＳによってナビゲータデータとして得られた磁気共鳴信号を、ローデータとして用いることによって、その被検体４０のナビゲータ領域ＮＡについての位相プロファイルを生成する。すなわち、本実施形態においては、スキャン部２によるナビゲータシーケンスＮＳの実施において、ＲＦコイル部１４がナビゲータ領域ＮＡへＲＦパルスを送信し、その撮影領域において発生する磁気共鳴信号を、ＲＦコイル部１４が受信することによってデータ収集部２４に収集されたナビゲータデータを一次元フーリエ（Ｆｏｕｒｉｅｒ）変換し、一次元フーリエ変換したデータの位相Ｐと位置Ｌとの関係をプロット（ｐｌｏｔ）することにより位相プロファイルを生成する。
位相補正部３３３は、位相プロファイル生成部３３２が生成した位相プロファイルの補正を実施する。本実施形態において、位相補正部３３３は、折り返された位相プロファイルを補正する。例えば、位相プロファイルが折り返されて、０°から３６０°の範囲内にある位相を補正し、実際の位相とする。
位置検出部３３４は、コンピュータとコンピュータに所定のデータ処理を実行させるプログラムを記憶するメモリとを有しており、制御部３０からの制御信号に基づいて、被検体４０における組織の位置を検出する。本実施形態において、位置検出部３３４は、位相プロファイル生成部３３２により１次元フーリエ変換が実施されたナビゲータデータの位相プロファイルにおいて、被検体４０における組織の境界ｂの位置Ｌを検出する。例えば、予め閾値Ｐｔｈを設定しておき、位相プロファイルにおける位相Ｐが閾値Ｐｔｈとなる位置Ｌｔｈを境界ｂとして検出する。被検体４０における組織の境界ｂの位置Ｌとして、例えば、被検体４０における肺臓と肝臓の境界にある横隔膜の位置を検出する。
表示部３４は、ディスプレイ（ｄｉｓｐｌａｙ）などの表示デバイスにより構成されており、制御部３０からの制御信号に基づいて、表示画面に画像を表示する。表示部３４は、例えば、オペレータによって操作部３２に操作データが入力される入力項目についての画像を表示画面に表示する。また、表示部３４は、データ処理部３３が生成する被検体４０のスライス画像を表示する。
（動作）
以下より、本実施形態の磁気共鳴イメージング装置１を用いて、被検体４０を撮像する際の動作について説明する。
図３は、本発明に係る一実施形態において、被検体を撮像する際の動作を示すフロー（ｆｌｏｗ）図である。また、図４は、本発明に係る一実施形態において、被検体４０をスキャンする際のシーケンスを示すシーケンス図である。横軸が時間軸ｔである。また、図５は、本発明に係る一実施形態において、ナビゲータ領域ＮＡ及びイメージング領域ＩＡを示すコロナル画像を示す図である。図５におけるｚ方向は、後述する図８における横軸の位置Ｌの方向である。なお、ここでは、黒色領域は肺臓を示す。また、灰色領域は肝臓を示し、肺臓と肝臓の間に横隔膜が位置する。また、横隔膜にほぼ直交するｚ軸方向に長辺が配置された矩形領域は、ナビゲータシーケンスＮＳを実施するナビゲータ領域ＮＡである。また、コロナル画像における肝臓において、横隔膜にほぼ平行なｘ軸方向に平行となるように配置された領域は、イメージングシーケンスＩＳを実施するイメージング領域ＩＡである。図６は、本発明に係る一実施形態におけるナビゲータシーケンスを示すパルスシーケンスを示す図である。Ｇｓはスライス選択傾斜磁場、Ｇφは位相エンコーディング（ｅｎｃｏｄｉｎｇ）傾斜磁場、Ｇｆは周波数エンコーディング傾斜磁場を示し、ＲＦは高周波パルスを示す。なお、ここでは、縦軸が強度を示し、横軸が時間を示す。図７は、本発明に係る一実施形態におけるシリンドリカルグラジエントエコー型ナビゲータシーケンスであるパルスシーケンスを示す図である。ＲＦは高周波パルスを示す。ｘｇｒａｄはｘ軸方向に印加する傾斜磁場、ｙｇｒａｄはｙ軸方向に印加する傾斜磁場、ｚｇｒａｄはｚ軸方向に印加する傾斜磁場を示し、ＲＦは高周波パルスを示す。なお、ここでは、縦軸が強度を示し、横軸が時間を示す。
まず、図３に示すように、ナビゲータシーケンスＮＳを実施する（ＳＴ１０）。
ここでは、被検体４０における組織の境界ｂの位置Ｌを検出するために、スキャン部２が被検体４０のナビゲータシーケンスＮＳを実施することによりナビゲータデータを取得する。
例えば、ナビゲータシーケンスＮＳは、図４に示すように、ナビゲータシーケンスＮＳの実施が開始される時点ｔ１から、所定時間Ｄ１が経過した時点ｔ２までの間において実施される。具体的には、ナビゲータシーケンスＮＳを実施する前に、図５に示すように、コロナル画像５０を表示部３４が表示する。そして、表示部３４が表示するコロナル画像５０上にナビゲータシーケンスを実施する領域をナビゲータ領域ＮＡとして、例えば、オペレータが操作部３２により設定する。設定されるナビゲータ領域ＮＡはイメージングシーケンスＩＳを実施するイメージング領域ＩＡと重なっていてもよい。図５に示すように、本実施形態の場合、ナビゲータ領域ＮＡは、例えば肺臓６０と肝臓６１との間に位置する横隔膜６２と交差し、横隔膜６２とほぼ垂直であるｚ軸方向に平行な長辺を有する矩形領域として設定する。
そして、ナビゲータ領域ＮＡにおいて、図６に示すように、フリップアングル（ｆｌｉｐ ａｎｇｌｅ）９０°の高周波パルスＲＦ１と共に第１のスライス選択傾斜磁場Ｇｓ１を印加することによって、被検体４０において横隔膜を含む第１のスライス面を選択的に９０°励起する。その後、第２のスライス選択傾斜磁場Ｇｓ２を被検体４０に印加することによって位相を戻した後、フリップアングル１８０°の高周波パルスＲＦ２と共に第３のスライス選択傾斜磁場Ｇｓ３と第１の位相エンコーディング傾斜磁場Ｇφ１とを印加することによって、横隔膜を含む領域において第１のスライス面と交差する第２のスライス面を１８０°励起する。
そして、周波数エンコードするように第１と第２との周波数エンコーディング傾斜磁場Ｇｆ１，Ｇｆ２を印加し、被検体４０において第１のスライス面と第２のスライス面とが交差する領域からの磁気共鳴信号Ｅを、スキャン部２がナビゲータデータとして取得する。
また、ナビゲータ領域ＮＡにおいて、図７に示すように、例えば、ＲＦパルスの送信と同時に、ｙ方向およびｚ方向に極性が交互に連続的に変化するように傾斜磁場を印加し線状に被検体４０のスライスを選択し、そしてｘ方向に読み取り傾斜磁場を印加するシリンドリカルグラジエントエコー型ナビゲータシーケンスを実施してもよい。
ここで、本実施形態において図５に示すように、ナビゲータデータは、横隔膜とほぼ垂直な方向（ｚ軸方向）の位置情報を必要とし、横隔膜とほぼ平行な方向（ｘ軸方向）の位置情報は必要としない。従って、ナビゲータシーケンスＮＳは、ｘ軸方向の空間情報を取得する位相エンコード傾斜磁場を印加しないナビゲータデータを取得するシーケンスである。
つぎに、図３に示すように、ナビゲータデータの位相プロファイルを生成する（ＳＴ２０）。
ここでは、ステップＳＴ１０においてナビゲータシーケンスＮＳを実施することにより取得したナビゲータデータに対して、位相プロファイル生成部３３２が１次元フーリエ変換を実施する。そして、一次元フーリエ変換を実施したナビゲータデータにおける位相と、図５に示すナビゲータ領域ＮＡにおけるｚ軸方向の位置との関係をプロットすることにより、位相プロファイルを生成する。ここで、図８は、本発明に係る一実施形態において、ナビゲータデータを１次元フーリエ変換し、位相をプロットした位相プロファイルを示す図である。縦軸は位相Ｐ、横軸は位置Ｌを示す。ここで、位置Ｌは、図５に示すナビゲータ領域ＮＡにおけるｚ軸方向の位置である。また、位相Ｐが小さい領域Ｓは図５における肺臓６０の位相を示し、位相Ｐが大きい領域Ｂは図５における肝臓６１の位相を示す。以下の図においても同様である。
具体的には、位相プロファイル生成部３３２がスキャン部２によって取得したナビゲータデータに対し、図５に示すｚ軸方向に１次元フーリエ変換を実施する。
そして、図８に示すように、１次元フーリエ変換を実施したナビゲータデータの位相と位置との関係をプロットし、ナビゲータデータの位相プロファイルを生成する。
つぎに、図３に示すように、位相プロファイルの補正を実施する（ＳＴ３０）。
ここでは、ステップＳＴ２０において位相をプロットすることにより生成した位相プロファイルの補正を、位相補正部３３３が実施する。ここで、図９は、本発明に係る一実施形態において、位相の補正を実施した位相プロファイルを示す図である。縦軸は位相Ｐ、横軸は位置Ｌを示す。
具体的には、図８に示すように、位相プロファイルにおいて、破線で囲んだ位相Ｐが急激に低下している部分１００、すなわち位相Ｐが離散的に変化している部分が存在する。これは、位相Ｐが３６０°になると０°に折り返され、位相Ｐが３６０°を超えると３６０°を超えた分が０°から起算した値になるためである。従って、０°に折り返された位相を実際の位相に補正するために、折り返された位相プロファイルの補正を実施する。例えば、０°位相は３６０°位相に、１０°位相は３７０°位相とする補正を実施する。補正を実施した後の位相プロファイルは、図９に示すようなプロファイルを奏する。より具体的には、例えば、隣り合うピクセルの位相差が１８０°を超える場合、その差が１８０°以下になるように３６０°×ｎ（ｎは正負の整数）を片方のピクセルにおける位相に加え、全領域において隣り合うピクセルの位相差が１８０°以下になるように、位相の補正を実施して、位相プロファイルにおける折り返しを除去する。
つぎに、図３に示すように、位相プロファイルにおいて、被検体の組織の境界を検出する（ＳＴ４０）。
ここでは、ステップＳＴ３０において位相の補正が実施された位相プロファイルにおいて、被検体４０の組織の境界ｂの位置Ｌを、例えばエッジ検出法、ＬＳＱ法、相関係数法、相互情報量法などを用いて、位置検出部３３４が検出する。
ここで、図１０は、本発明の一実施形態において、組織の境界ｂの位置Ｌが検出された位相プロファイルを示す図である。縦軸は位相Ｐ、横軸は位置Ｌを示す。
具体的には、図１０に示すように、例えば、エッジ検出法により、位置検出部３３４が、位相の補正を実施した位相プロファイルにおける位相Ｐの最大位相Ｐｍａｘと最小位相Ｐｍｉｎの平均値（Ｐｍａｘ＋Ｐｍｉｎ）／２を算出し、この平均値を閾値Ｐｔｈとする。そして、位相プロファイルにおいて、位相が閾値Ｐｔｈ以上となる位置Ｌが所定の長さ、例えば１０ピクセル分以上継続した場合、位相Ｐが閾値Ｐｔｈとなる位置Ｌｔｈを被検体４０における組織の境界ｂの位置Ｌとして、位置検出部３３４が検出する。
また、例えば、ＬＳＱ法（最小二乗法）により、ナビゲータデータの空間プロファイルｆ（ｘ）と基準ナビゲータデータプロファイルｆｒ（ｘ）を用い、以下の数式（１）で示される最小二乗誤差ｌ(ｄ)が最小となる変位ｄを位置検出部３３４が算出し、算出した変位ｄを被検体４０における組織の境界ｂの位置Ｌとして、位置検出部３３４が検出する。
ｘ１、ｘ２：処理を行う区間
また、例えば、相関値から被検体４０における組織の境界ｂの位置Ｌを位置検出部３３４が検出する相関係数法、相互情報量から被検体４０における組織の境界ｂの位置Ｌを位置検出部３３４が検出する相互情報量法などを用いてもよい。
ステップＳＴ３０において補正を実施していない図８に示す位相プロファイルにおいて、本ステップの被検体の組織の境界を検出した場合、破線で囲んだ位相Ｐが急激に低下している部分１００にも位相が閾値Ｐｔｈ以上となる位置Ｌが所定の長さ、例えば、１０ピクセル分以上継続する場合がある。その場合、実際には部分１００には被検体４０における組織の境界ｂは存在しないが、部分１００における位相Ｐが閾値Ｐｔｈとなる位置Ｌｔｈを被検体４０における組織の境界ｂの位置Ｌとして検出してしまう可能性がある。従って、ステップＳＴ３０において、位相プロファイルの補正を実施する。
以上のように検出できる被検体４０における組織の境界として、例えば、被検体４０における肺臓６０と肝臓６１の境界にある横隔膜６２などが挙げられる。
つぎに、図３に示すように、組織の境界の位置が許容範囲内であるかの判断を行う（ＳＴ５０）。
ここでは、ステップＳＴ４０において位置検出部３３４が検出した位相プロファイルにおいて検出された被検体４０における組織の境界ｂの位置Ｌが許容範囲内であるかの判断を制御部３０が比較処理する。
具体的には、位置検出部３３４により検出された被検体４０における組織の境界ｂの位置Ｌと、予め定められた許容範囲の上下限値とを制御部３０が比較処理することにより、許容範囲内であるか否かを判断する。
境界ｂの位置と予め定められた許容範囲の上下限値とを比較した結果、境界ｂの位置が上限値以下であって下限値以上である場合、許容範囲内であると判断し、境界ｂの位置が上限値を超える、あるいは下限値未満である場合、許容範囲でないと判断する。
そして、図３に示すように、境界ｂの位置が許容範囲内でないと判断した場合（ＮＯ）においては、再びナビゲータシーケンスを実施（ＳＴ１０）、ナビゲータデータの位相プロファイルを生成（ＳＴ２０）、位相プロファイルの補正を実施（ＳＴ３０）、位相プロファイルにおいて、被検体の組織の境界を検出（ＳＴ４０）する。
つぎに、図３に示すように、イメージングシーケンスを実施する（ＳＴ６０）。
ここでは、ステップＳＴ５０において境界ｂの位置が許容範囲内であると判断した場合（ＹＥＳ）に、スキャン部２が被検体４０のイメージングシーケンスＩＳを実施することによりイメージングデータを取得する。
具体的には、図４、５に示すように、スキャン部２が被検体４０におけるイメージング領域ＩＡにおいて、イメージングシーケンスＩＳを実施して、イメージングデータとしての磁気共鳴信号を取得する。例えば、被検体４０の肝臓６１におけるイメージング領域ＩＡにおいて、イメージングシーケンスＩＳを実施してイメージングデータを取得する。このイメージングシーケンスＩＳは、図４に示すように、ナビゲータシーケンスが終了した時点ｔ２から所定の時間Ｄ２の経過後であって、イメージングシーケンスＩＳが開始する時点ｔ３から所定の時間Ｄ３を経過後の時点ｔ４までの間において実施される。
そして、ｋ空間において全ての位相エンコードステップに対応するように、イメージング領域ＩＡにおいて複数回イメージングシーケンスＩＳを実施して、イメージングデータを取得する。
つぎに、図３に示すように、スライス画像の画像再構成を実施する（ＳＴ７０）。
ここでは、ステップＳＴ６０において取得されたイメージングデータを用いて、画像再構成部３３１が被検体４０のスライス画像を画像再構成する。例えば、画像再構成部３３１は被検体４０における肝臓６１のスライス画像を画像再構成する。そして、画像再構成部３３１は、画像再構成したスライス画像を表示部３４に出力する。
以上のように、本発明の一実施形態は、ナビゲータ領域ＮＡを設定し、設定したナビゲータ領域ＮＡにおいてナビゲータシーケンスＮＳを実施する。そして、取得したナビゲータデータに対して、１次元フーリエ変換を実施する。そして、一次元フーリエ変換を実施したナビゲータデータにおける位相をプロットした位相プロファイルを生成する。そして、生成した位相プロファイルの補正を実施する。そしてエッジ検出法により、補正を実施した位相プロファイルの最大位相と最小位相の平均値を算出し、その値を閾値Ｐｔｈと設定する。そして、位相プロファイルにおいて、位相が閾値Ｐｔｈ以上となる位置Ｌが所定の長さ、例えば、１０ピクセル分以上継続した場合、位相Ｐが閾値Ｐｔｈとなる位置Ｌｔｈを被検体４０における組織の境界ｂの位置Ｌとして検出する。また、ＬＳＱ法の場合、ナビゲータデータの空間プロファイルｆ（ｘ）と基準ナビゲータデータプロファイルｆｒ（ｘ）との最小二乗誤差ｌ（ｄ）が最小となるｄを算出し、変位ｄを被検体４０における組織の境界ｂの位置Ｌとする。そして、境界ｂの位置Ｌが許容範囲内である場合、イメージングシーケンスＩＳを実施してイメージングデータを取得する。そして、取得されたイメージングデータを用いて、画像再構成部が被検体のスライス画像を画像再構成する。
このように、ナビゲータ領域ＮＡとイメージング領域ＩＡとが重なっていても、位相をプロットした位相プロファイルは、イメージングスキャンの干渉による信号の乱れが発生しないため、補正を実施した位相プロファイルにおいて、被検体における組織の境界の位置を正確に検出することができる。そして、組織の境界の位置を正確に検出することができるため、画像再構成されたスライス画像において体動アーチファクトの発生を防止することができ、画像品質の良いスライス画像を得ることができる。
＜第２の実施形態＞
（装置構成）
図１１は、本発明に係る一実施形態におけるデータ処理部３３の構成を示すブロック図である。本データ処理部は、本発明の実施形態の一例である。
第２の実施形態は、磁気共鳴イメージング装置１の構成である範囲設定部３３５以外の部分は第１の実施形態と同じである。そのため、重複する箇所については、記載を省略する。
図１１に示すように、データ処理部３３は、画像再構成部３３１と、位相プロファイル生成部３３２と、位相補正部３３３と、位置検出部３３４と、範囲設定部３３５を有する。
範囲設定部３３５は、例えば、位相補正部３３３により補正が実施された位相プロファイルにおいて、被検体４０における組織の境界ｂの位置Ｌを検出する範囲を設定する。本実施形態において、範囲設定部３３５は、オペレータが操作部３２により、位相プロファイルに設定した被検体４０における組織の境界ｂの位置Ｌを検出する範囲を認識し、決定する。
（動作）
以下より、本実施形態の磁気共鳴イメージング装置１を用いて、被検体４０を撮像する際の動作について図１２により説明する。
動作フローにおけるステップＳＴ３５である被検体の組織の境界を検出する範囲の設定及び、ステップＳＴ３６である被検体の組織の境界を検出以外の部分は第１の実施形態と同じである。そのため、重複する箇所については、記載を省略する。
図１２は、本発明に係る一実施形態において、被検体を撮像する際の動作を示すフロー図である。
図１２に示すように、位相プロファイルにおける組織の境界を検出する範囲を設定する（ＳＴ３５）。
ここでは、ステップＳＴ３０において、補正を実施した位相プロファイルにおいて、被検体４０における組織の境界ｂの位置Ｌを検出する範囲２００の設定を、オペレータが操作部３２により実施する。ここで、図１３は、本発明に係る一実施形態における組織の境界ｂの位置Ｌを検出する範囲を設定した位相プロファイルを示す図である。破線部が、組織の境界ｂの位置を検出する範囲２００である。
具体的には、図１３に示すように、ステップＳＴ３０において、オペレータが操作部３２により、位相の補正を実施した位相プロファイルに被検体４０における組織の境界ｂの位置Ｌを検出する矩形状の領域を設定する。例えば、この矩形状の領域は、ナビゲータの分解能が２５６の場合、位相プロファイルにおける位置Ｌの中央を中心に１２８ピクセル（ｐｉｘｅｌ）分の位置Ｌの領域を、境界ｂを検出する範囲２００として設定する。
つぎに、図１２に示すように、位相プロファイルにおいて、被検体の組織の境界を検出する（ＳＴ３６）。
ここでは、ステップＳＴ３５において位相プロファイルにおける組織の境界ｂの位置Ｌを検出する範囲２００が設定された位相プロファイルにおいて、被検体４０における組織の境界ｂの位置Ｌを、例えばエッジ検出法、ＬＳＱ法、相関係数法、相互情報量法などを用いて、位置検出部３３４が検出する。
ここで、図１４は、本発明の一実施形態において、組織の境界ｂの位置Ｌが検出された位相プロファイルを示す図である。縦軸は位相Ｐ、横軸は位置Ｌを示す。
具体的には、図１４に示すように、例えば、エッジ検出法により、境界ｂを検出する範囲２００を設定した位相プロファイルにおいて、設定した範囲２００内における位相Ｐの最大値Ｐｍａｘと最小値Ｐｍｉｎの平均値（Ｐｍａｘ＋Ｐｍｉｎ）／２を、位置検出部３３４が算出し、この平均値を閾値Ｐｔｈとする。そして、位相プロファイルにおいて、位相が閾値Ｐｔｈ以上となる位置Ｌが所定の長さ、例えば１０ピクセル分以上継続した場合、位相Ｐが閾値Ｐｔｈとなる位置Ｌｔｈを被検体４０における組織の境界ｂの位置Ｌとして、位置検出部３３４が検出する。
また、ＬＳＱ法により、ナビゲータデータの空間プロファイルｆ（ｘ）と基準ナビゲータデータプロファイルｆｒ（ｘ）を用い、以下の数式（１）で示される最小二乗誤差ｌ(ｄ)が最小となる変位ｄを位置検出部３３４が算出し、算出した値を被検体４０における組織の境界ｂの位置Ｌとして、位置検出部３３４が検出する。
また、例えば、相関値から被検体４０における組織の境界ｂの位置Ｌを位置検出部３３４が検出する相関係数法、相互情報量から被検体４０における組織の境界ｂの位置Ｌを位置検出部３３４が検出する相互情報量法などを用いてもよい。
以上のように、本発明の第２の実施形態は、位相の補正を実施した位相プロファイルにおいて、オペレータが被検体４０における組織の境界ｂの位置Ｌを検出する範囲２００を設定する。そして、位相プロファイルに設定した範囲２００において、エッジ検出法により補正を実施した位相プロファイルの最大位相と最小位相の平均値を算出し、その値を閾値Ｐｔｈと設定し、位相プロファイルにおいて、位相が閾値Ｐｔｈ以上となる位置Ｌが所定の長さ、例えば１０ピクセル分以上継続した場合、位相Ｐが閾値Ｐｔｈとなる位置Ｌｔｈを被検体４０における組織の境界ｂの位置Ｌとして検出する。また、ＬＳＱ法の場合、位相プロファイルに設定した範囲２００において、ナビゲータデータの空間プロファイルｆ（ｘ）と基準ナビゲータデータプロファイルｆｒ（ｘ）との最小二乗誤差ｌ（ｄ）が最小となるｄを算出し、変位ｄを被検体４０における組織の境界ｂの位置Ｌとする。
このように、位相プロファイルに被検体における組織の境界の位置を検出する範囲を設定することにより、位相プロファイルにおいて被検体における組織の境界の位置をより正確に検出することができる。そして、組織の境界の位置をより正確に検出することができるため、画像再構成されたスライス画像において体動アーチファクトの発生を防止することができ、画像品質の良いスライス画像を得ることができる。
＜第３の実施形態＞
（装置構成）
図１５は、本発明に係る一実施形態におけるデータ処理部３３の構成を示すブロック図である。本データ処理部は、本発明の実施形態の一例である。
第３の実施形態は、磁気共鳴イメージング装置１の構成である位相プロファイル平滑部３３６以外の部分は第２の実施形態と同じである。そのため、重複する箇所については、記載を省略する。
図１５に示すように、データ処理部３３は、画像再構成部３３１と、位相プロファイル生成部３３２と、位相補正部３３３と、位置検出部３３４と、範囲設定部３３５と、位相プロファイル平滑部３３６とを有する。
位相プロファイル平滑部３３６は、例えば、位相補正部３３３により補正が実施された位相プロファイルが不連続の場合、位相プロファイルの平滑化を行う。本実施形態において、位相プロファイル平滑部３３６は、例えば、入力される位相プロファイルのデータを近傍平均やメディアンフィルタなどのフィルタによりフィルタリングし、位相プロファイルの平滑化を行う。
（動作）
以下より、本実施形態の磁気共鳴イメージング装置１を用いて、被検体４０を撮像する際の動作について図１６により説明する。
動作フローにおけるステップＳＴ３１である位相プロファイルの平滑化以外の部分は第２の実施形態と同じである。そのため、重複する箇所については、記載を省略する。
図１６は、本発明に係る一実施形態において、被検体を撮像する際の動作を示すフロー図である。
図１６に示すように、位相プロファイルの平滑化を実施する（ＳＴ３１）。
ここでは、ステップＳＴ３０において、位相の補正を実施した位相プロファイルの平滑化を位相プロファイル平滑部３３６が実施する。ここで、図１７は、本発明に係る一実施形態における位相プロファイルを示す図である。図１７（ａ）は、位相Ｐの不連続を平滑化する前の位相プロファイルを示す図であり、図１７（ｂ）は、位相Ｐの不連続を平滑化した後の位相プロファイルを示す図である。
具体的には、図１７（ａ）に示すように、ステップＳＴ３０において位相の補正を実施した位相プロファイルにおいて、位相の不連続が生じる場合がある。その場合、例えば、近傍平均やメディアンフィルタ(ｍｅｄｉａｎ ｆｉｌｔｅｒ)などの平滑化フィルタをかけて位相プロファイルの平滑化を行う。図１７（ａ）に示す位相プロファイルの平滑化を行うことにより、図１７（ｂ）に示す位相プロファイルとなる。
ここで、近傍平均は、例えば、位相プロファイルにおけるある位置Ｌにおいて、当該位置Ｌにおける位相Ｐと隣り合う位置Ｌにおける位相Ｐとの３点の位相Ｐの平均値を当該位置Ｌにおける位相Ｐとするフィルタリング（ｆｉｌｔｅｒｉｎｇ）である。メディアンフィルタは、例えば、位相プロファイルにおけるある位置Ｌにおいて、当該位置Ｌにおける位相Ｐと隣り合う位置Ｌにおける位相Ｐとの３点の位相Ｐを位相の大きさ順に並べたときの中央に位置する位相を当該位置Ｌにおける位相Ｐとするフィルタリングである。
以上のように、本発明の第３の実施形態は、補正を実施した位相プロファイルの平滑化を位相プロファイル平滑部３３６が実施する。そして、平滑化として、例えば、近傍平均やメディアンフィルタなどの平滑化フィルタをかけて、位相プロファイルの平滑化を実施する。
このように、位相プロファイルを平滑化することにより、被検体における組織の境界の位置をより正確に検出することができる。そして、組織の境界の位置をより正確に検出することができるため、画像再構成されたスライス画像において体動アーチファクトの発生を防止することができ、画像品質の良いスライス画像を得ることができる。
なお、上記の本実施形態における位置検出部３３４は、本発明の境界検出部に相当する。
なお、本発明の実施に際しては、上記した実施形態に限定されるものではなく、種々の変形形態を採用することができる。
本発明の第３の実施形態において、位相プロファイルの平滑化を実施した後に、位相プロファイルにおける組織の境界の位置を検出する範囲を設定しているが、これには限定されず、位相プロファイルにおける組織の境界の位置を検出する範囲を設定した後に、位相プロファイルの平滑化を実施してもよい。
また、本発明の実施形態において、イメージングシーケンスは、被検体における組織の境界の位置を検出した後に実施しているが、これに限定されず、被検体における組織の境界の位置を検出する前に実施してもよい。また、ナビゲータシーケンスを実施する前であってもよい。
また、本発明の実施形態において、取得したナビゲータデータを用いて生成した位相プロファイルにおいて、被検体における組織の境界の位置を検出し、その境界の位置によりイメージングシーケンスを実施するか判断しているが、これには限定されず、境界の位置によりイメージングデータの選択をしてもよい。また、境界の位置によりイメージングデータの補正を行ってもよい。
また、本発明の実施形態において、イメージングシーケンスは、位置検出部により検出された被検体における境界の位置と、予め定められた許容範囲の上下限値とを比較し、許容範囲内の場合、イメージングシーケンスを実施しているが、これには限定されず、許容範囲外の場合、取得したイメージングデータを除去し、許容範囲内の場合、取得したイメージングデータを選択してもよい。また、位置検出部により検出された被検体における境界の位置に基づいて、取得したイメージングデータを補正してもよい。
図１は、本発明に係る一実施形態のＲＦコイル部により構成される磁気共鳴イメージング装置の構成を示す構成図である。 図２は、本発明に係る一実施形態におけるデータ処理部の構成を示すブロック図である。 図３は、本発明に係る一実施形態において、被検体を撮像する際の動作を示すフロー図である。 図４は、本発明に係る一実施形態において、被検体をスキャンする際のシーケンスを示すシーケンス図である。 図５は、本発明に係る一実施形態において、ナビゲータ領域ＮＡ及びイメージング領域ＩＡを示すコロナル画像を示す図である。 図６は、本発明に係る一実施形態におけるナビゲータシーケンスを示すパルスシーケンスを示す図である。 図７は、本発明に係る一実施形態におけるシリンドリカルナビゲータシーケンスを示すパルスシーケンスを示す図である。 図８は、本発明に係る一実施形態において、ナビゲータデータを１次元フーリエ変換し、位相をプロットした位相プロファイルを示す図である。 図９は、本発明に係る一実施形態において、位相の補正を実施した位相プロファイルを示す図である。 。 図１０は、本発明に係る一実施形態において、組織の境界の位置が検出された位相プロファイルを示す図である。 図１１は、本発明に係る一実施形態におけるデータ処理部の構成を示すブロック図である。 図１２は、本発明に係る一実施形態において、被検体を撮像する際の動作を示すフロー図である。 図１３は、本発明に係る一実施形態における組織の境界の位置を検出する範囲を設定した位相プロファイルを示す図である。 図１４は、本発明に係る一実施形態において、組織の境界の位置が検出された位相プロファイルを示す図である。 図１５は、本発明に係る一実施形態におけるデータ処理部の構成を示すブロック図である。 図１６は、本発明に係る一実施形態において、被検体を撮像する際の動作を示すフロー図である。 図１７は、本発明に係る一実施形態における位相プロファイルを示す図である。 図１８は、従来技術を説明するためのナビゲータ領域ＮＡ及びイメージング領域ＩＡを示すコロナル画像を示す図である。 図１９は、従来技術を説明するための信号強度プロファイルを示す図である。
符号の説明
１：磁気共鳴イメージング装置、２:スキャン部、３：操作コンソール部、１１：静磁場空間、１２：静磁場マグネット部、１３：勾配コイル部、１４：ＲＦコイル部、１５：クレードル、２２：ＲＦ駆動部、２３：勾配駆動部、２４：データ収集部、３０：制御部、３１：記憶部、３２：操作部、３３：データ処理部、３４：表示部、４０：被検体、６０：肺臓、６１：肝臓、６２：横隔膜、３３１：画像再構成部、３３２：位相プロファイル生成部、３３３：位相補正部、３３４：位置検出部、３３５：範囲設定部、３３６：位相プロファイル平滑部、ＮＡ：ナビゲータ領域、ＩＡ：イメージング領域

Claims (18)

  1. 被検体において体動する組織を含むナビゲータ領域から磁気共鳴信号をナビゲータデータとして得るナビゲータシーケンスを実施すると共に、前記被検体においてイメージング領域から磁気共鳴信号をイメージングデータとして得るイメージングシーケンスを実施するスキャンを実施することによって、前記イメージング領域について画像を生成する磁気共鳴イメージング装置であって、
    前記ナビゲータデータの位相と前記ナビゲータ領域の位置との関係を示すように位相プロファイルを生成する位相プロファイル生成部と、
    前記位相プロファイル生成部によって生成された位相プロファイルについて折り返りを補正する位相補正部と、
    前記位相補正部によって補正された位相プロファイルに基づいて、前記ナビゲータ領域において体動する組織の位置を検出する位置検出部と
    を有する、
    磁気共鳴イメージング装置。
  2. 前記スキャンとして、重複する前記イメージング領域と前記ナビゲータ領域とのスキャンを実施する、
    請求項１に記載の磁気共鳴イメージング装置。
  3. 前記位相補正部により位相の補正が実施された前記磁気共鳴信号の位相プロファイルにおいて、前記被検体における組織の位置を検出する範囲を設定する範囲設定部を有する、
    請求項１又は２に記載の磁気共鳴イメージング装置。
  4. 前記位置検出部は、エッジ検出法により前記被検体における組織の位置を検出する、
    請求項１から３のいずれか１項に記載の磁気共鳴イメージング装置。
  5. 前記位置検出部は、ＬＳＱ法により前記被検体における組織の位置を検出する、
    請求項１から３のいずれか１項に記載の磁気共鳴イメージング装置。
  6. 前記位置検出部は、相関係数法により前記被検体における組織の位置を検出する、
    請求項１から３のいずれか１項に記載の磁気共鳴イメージング装置。
  7. 前記位置検出部は、相互情報量法により前記被検体における組織の位置を検出する、
    請求項１から３のいずれか１項に記載の磁気共鳴イメージング装置。
  8. 前記組織が横隔膜である、
    請求項４から７のいずれか１項に記載の磁気共鳴イメージング装置。
  9. 前記位相補正部により位相の補正が実施された位相プロファイルを平滑化させる位相プロファイル平滑部を有する、
    請求項１から８のいずれか１項に記載の磁気共鳴イメージング装置。
  10. 静磁場空間内において被検体にＲＦパルスを送信し、当該被検体から磁気共鳴信号を収集するスキャンを実施することにより、前記被検体のスライス画像を生成するためのイメージングデータと前記被検体における組織の位置を検出するためのナビゲータデータを取得し、取得した当該ナビゲータデータを解析するナビゲータデータ解析方法であって、
    前記ナビゲータデータの位相と前記ナビゲータ領域の位置との関係を示すように位相プロファイルを生成する第１ステップと、
    前記第１ステップにおいて生成された前記位相プロファイルの折り返りを補正する第２ステップと、
    前記第２ステップにおいて折り返りの補正がされた位相プロファイルに基づいて、前記ナビゲータ領域において体動する組織の位置を検出する第３ステップと
    を含む、
    ナビゲータデータ解析方法。
  11. 前記スキャンとして、重複する前記ナビゲータデータを取得するナビゲータ領域と前記イメージングデータを取得するイメージング領域とのスキャンを実施する、
    請求項１０に記載のナビゲータデータ解析方法。
  12. 前記第２ステップの後、位相の補正が実施された前記位相プロファイルにおいて、前記被検体における組織の位置を検出する範囲を設定するステップを有する、
    請求項１０又は１１に記載のナビゲータデータ解析方法。
  13. 前記第３ステップは、エッジ検出法により前記被検体における組織の位置を検出する、
    請求項１０から１２のいずれか１項に記載のナビゲータデータ解析方法。
  14. 前記第３ステップは、ＬＳＱ法により前記被検体における組織の位置を検出する、
    請求項１０から１２のいずれか１項に記載のナビゲータデータ解析方法。
  15. 前記第３ステップは、相関係数法により前記被検体における組織の位置を検出する、
    請求項１０から１２のいずれか１項に記載のナビゲータデータ解析方法。
  16. 前記第３ステップは、相互情報量法により前記被検体における組織の位置を検出する、
    請求項１０から１２のいずれか１項に記載のナビゲータデータ解析方法。
  17. 前記組織が横隔膜である、
    請求項１３から１６のいずれか１項に記載のナビゲータデータ解析方法。
  18. 前記第２ステップの後、前記第２ステップにおいて位相の補正が実施された前記位相プロファイルを平滑化するステップを有する、
    請求項１０から１７のいずれか１項に記載のナビゲータデータ解析方法。
JP2007221489A 2007-08-28 2007-08-28 磁気共鳴イメージング装置及びナビゲータデータ解析方法 Active JP5177367B2 (ja)

Priority Applications (2)

Application Number Priority Date Filing Date Title
JP2007221489A JP5177367B2 (ja) 2007-08-28 2007-08-28 磁気共鳴イメージング装置及びナビゲータデータ解析方法
US12/200,818 US7782053B2 (en) 2007-08-28 2008-08-28 Magnetic resonance imaging apparatus and navigator data analyzing method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2007221489A JP5177367B2 (ja) 2007-08-28 2007-08-28 磁気共鳴イメージング装置及びナビゲータデータ解析方法

Publications (2)

Publication Number Publication Date
JP2009050544A JP2009050544A (ja) 2009-03-12
JP5177367B2 true JP5177367B2 (ja) 2013-04-03

Family

ID=40431176

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2007221489A Active JP5177367B2 (ja) 2007-08-28 2007-08-28 磁気共鳴イメージング装置及びナビゲータデータ解析方法

Country Status (2)

Country Link
US (1) US7782053B2 (ja)
JP (1) JP5177367B2 (ja)

Families Citing this family (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP5227244B2 (ja) * 2008-08-18 2013-07-03 日立アロカメディカル株式会社 画像処理装置
JP5305925B2 (ja) * 2009-01-05 2013-10-02 ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー 磁気共鳴イメージング装置
JP2011131000A (ja) * 2009-12-25 2011-07-07 Ge Medical Systems Global Technology Co Llc 磁気共鳴イメージング装置
JP5461255B2 (ja) * 2010-03-17 2014-04-02 ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー 磁気共鳴イメージング装置
JP5704876B2 (ja) * 2010-09-29 2015-04-22 ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー 磁気共鳴イメージング装置
JP5977158B2 (ja) * 2012-11-30 2016-08-24 ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー 検出装置、磁気共鳴装置、検出方法、およびプログラム
US9648531B2 (en) * 2013-01-20 2017-05-09 Eci Telecom Ltd. Communication services to a moving platform
JP6109598B2 (ja) * 2013-02-26 2017-04-05 東芝メディカルシステムズ株式会社 磁気共鳴イメージング装置
JP6489744B2 (ja) * 2014-02-28 2019-03-27 ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー 磁気共鳴装置およびプログラム
DE102015222835B4 (de) * 2015-11-19 2019-06-13 Siemens Healthcare Gmbh Magnetresonanzbildgebungsverfahren mit simultaner Bildaufnahme von mehreren Teilvolumen mit einer synchronen Bildaufnahme von Navigatoren
CN110226099B (zh) * 2017-01-23 2022-04-29 皇家飞利浦有限公司 在对象运动期间采集四维磁共振数据
US20190120923A1 (en) * 2017-10-24 2019-04-25 Bruker Biospin Mri Gmbh Method for generating positron-emission tomography (pet) images
US11163029B2 (en) * 2019-08-14 2021-11-02 GE Precision Healthcare LLC MRI system with improved navigator

Family Cites Families (19)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6690963B2 (en) 1995-01-24 2004-02-10 Biosense, Inc. System for determining the location and orientation of an invasive medical instrument
CA2187964C (en) * 1996-10-16 2005-02-01 Kecheng Liu Sliding interleaved motsa for magnetic resonance imaging
JPH10277010A (ja) 1997-04-10 1998-10-20 Toshiba Corp Ｍｒｉ装置
US6275720B1 (en) 1999-03-16 2001-08-14 General Electric Company Iterative method for reducing errors in navigator echo MR imaging
JP3454760B2 (ja) * 1999-10-22 2003-10-06 ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー 位相分布測定方法および装置、位相補正方法および装置、並びに、磁気共鳴撮像装置
US6535754B2 (en) 2001-03-30 2003-03-18 The United States Of America As Represented By The Department Of Health And Human Services Lever coil sensor for respiratory and cardiac motion
US6718004B2 (en) 2001-06-28 2004-04-06 General Electric Company Methods and apparatus for coronary-specific imaging reconstruction
US6426990B1 (en) 2001-06-28 2002-07-30 General Electric Company Methods and apparatus for coronary-specific imaging reconstruction
US6771997B2 (en) 2001-09-11 2004-08-03 The Board Of Trustees Of The Leland Stanford Junior University Respiratory compensation in MRI coronary imaging using diminishing variance
US7012603B2 (en) 2001-11-21 2006-03-14 Viatronix Incorporated Motion artifact detection and correction
US7245786B2 (en) * 2002-05-10 2007-07-17 976076 Alberta Inc. Filtering artifact from fMRI data using the stockwell transform
US20040092809A1 (en) 2002-07-26 2004-05-13 Neurion Inc. Methods for measurement and analysis of brain activity
US6894494B2 (en) 2002-11-25 2005-05-17 Her Majesty The Queen In Right Of Canada, As Represented By The Minister Of National Defence Of Her Majesty's Canadian Government Method and device for correcting organ motion artifacts in MRI systems
US7012428B1 (en) 2003-12-17 2006-03-14 General Electric Company Method and apparatus of reducing artifacts in phase-cycled steady-state free precession imaging
US7945305B2 (en) 2005-04-14 2011-05-17 The Board Of Trustees Of The University Of Illinois Adaptive acquisition and reconstruction of dynamic MR images
US8155419B2 (en) * 2005-05-04 2012-04-10 Mayo Foundation For Medical Education And Research MRI acquisition using sense and highly undersampled fourier space sampling
JP2007098026A (ja) 2005-10-07 2007-04-19 Ge Medical Systems Global Technology Co Llc 磁気共鳴イメージング装置
JP4639136B2 (ja) 2005-10-19 2011-02-23 ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー 磁気共鳴イメージング装置
JP4717608B2 (ja) 2005-11-24 2011-07-06 ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー 磁気共鳴イメージング装置

Also Published As

Publication number Publication date
US20090066329A1 (en) 2009-03-12
US7782053B2 (en) 2010-08-24
JP2009050544A (ja) 2009-03-12

Similar Documents

Publication Publication Date Title
JP5177367B2 (ja) 磁気共鳴イメージング装置及びナビゲータデータ解析方法
US9797974B2 (en) Nonrigid motion correction in 3D using autofocusing with localized linear translations
JP4639136B2 (ja) 磁気共鳴イメージング装置
JP3976684B2 (ja) 画像における動きの影響を低減する方法および装置
US8290566B2 (en) Magnetic resonance imaging apparatus and image generating method
JP5171373B2 (ja) 磁気共鳴イメージング装置及びｒｆパルス印加方法
US20080205730A1 (en) Independent Motion Correction In Respective Signal Channels Of A Magnetic Resonance Imaging System
US9301706B2 (en) Magnetic resonance imaging system for non-contrast MRA and magnetic resonance signal acquisition method employed by the same
WO2005082245A1 (ja) 磁気共鳴イメージング装置および磁気共鳴イメージング装置のデータ処理方法
US7663365B2 (en) Magnetic resonance imaging apparatus and analysis method for fat suppression effect in magnetic resonance imaging apparatus
JP5818631B2 (ja) 磁気共鳴イメージング装置
KR20070038929A (ko) 자기공명촬상장치
JP6008839B2 (ja) 磁気共鳴イメージング装置および磁気共鳴イメージング方法
US8890526B2 (en) Method and apparatus for making distinction in a magnetic resonance imaging water-fat image
JP2010075573A (ja) 磁気共鳴イメージング装置
JP3993799B2 (ja) 磁気共鳴イメージング装置
US7560924B2 (en) Magnetic resonance imaging apparatus and magnetic resonance imaging method
JP4072879B2 (ja) 核磁気共鳴撮影装置
US8055038B2 (en) Image processing apparatus, image processing method, and magnetic resonance imaging apparatus
WO2003032830A1 (fr) Appareil et procede d'imagerie par resonance magnetique
JP2008055023A (ja) 磁気共鳴イメージング装置
JP2001292976A (ja) 核磁気共鳴イメージング方法及び装置
JP4675936B2 (ja) 核磁気共鳴撮影装置
JP4282152B2 (ja) Ｍｒｉ装置
JP2001149342A (ja) Ｍｒｉ装置及びｍｒイメージング方法

Legal Events

Date Code Title Description
A625 Written request for application examination (by other person)

Free format text: JAPANESE INTERMEDIATE CODE: A625

Effective date: 20100630

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20120724

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20121022

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20121211

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20121225

R150 Certificate of patent or registration of utility model

Ref document number: 5177367

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250