JP5437610B2 - 磁気共鳴イメージング装置 - Google Patents

磁気共鳴イメージング装置 Download PDF

Info

Publication number
JP5437610B2
JP5437610B2 JP2008223151A JP2008223151A JP5437610B2 JP 5437610 B2 JP5437610 B2 JP 5437610B2 JP 2008223151 A JP2008223151 A JP 2008223151A JP 2008223151 A JP2008223151 A JP 2008223151A JP 5437610 B2 JP5437610 B2 JP 5437610B2
Authority
JP
Japan
Prior art keywords
image
subject
slice position
images
data
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
JP2008223151A
Other languages
English (en)
Other versions
JP2010057532A (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.)
Hitachi Healthcare Manufacturing Ltd
Original Assignee
Hitachi Medical Corp
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 Hitachi Medical Corp filed Critical Hitachi Medical Corp
Priority to JP2008223151A priority Critical patent/JP5437610B2/ja
Publication of JP2010057532A publication Critical patent/JP2010057532A/ja
Application granted granted Critical
Publication of JP5437610B2 publication Critical patent/JP5437610B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Magnetic Resonance Imaging Apparatus (AREA)

Description

本発明は、磁気共鳴イメージング(以下、「MRI」という)装置において、被検体が動いた場合の撮像位置の設定に関する。
MRI装置は、被検体、特に人体の組織を構成する原子核スピンが発生する核磁気共鳴(以下、NMRという)信号を計測し、その頭部、腹部、四肢等の形態や機能を2次元的に或いは3次元的に画像化する装置である。撮像においては、NMR信号は、傾斜磁場によって異なる位相エンコードが付与されるとともに周波数エンコードされて、時系列データとして計測される。計測されたNMR信号は、2次元又は3次元フーリエ変換されることにより画像に再構成される。
上記MRIを用いた検査では、1回の検査において、同一のスライス位置にて、複数のコントラスト画像を撮像し、それらを比較することで診断が行われる。例えば、脳ドックの場合、スライス位置決めの後、T1強調画像、T2強調画像、FLAIR(もしくは、プロトン密度画像)、T2*強調画像が同一スライス位置にて撮像される(例えば非特許文献1)。
脳ドックのガイドライン2003 p.5、日本脳ドック学会発行
上述したように、MRI装置における検査では、同一のスライス位置にて、複数のコントラスト画像が撮像される。この場合に問題となるのが、検査がある程度進行した状態で、被検体が動いた場合である。被検体が検査中に動くと、撮像目的部位が設定済みのスライス位置からずれてしまう。そのため、各コントラスト画像に写っている部位が異なり、診断に支障を来たすことになる。このような場合、従来では、操作者は、スカウト画像を撮り直して、被検体が動く前のスライス位置と同じになるようにスライス位置を合わせ直す必要がある。実際には、操作者の手動で、スライス位置を完全に同位置に合わせることは困難であるため、スカウト画像も含め、全コントラスト画像を撮り直すことも有り得る。
そこで、本発明の目的は、MRI装置において、検査中に被検体が動いて、撮像目的部位が設定済みのスライス位置からずれた場合でも、容易に位置修正を可能とすることである。
上記目的を達成するために、本発明のMRI装置は、予め取得された3次元画像データと、被検体が動いた後に撮像された画像の少なくとも一部のデータと、から、被検体が動く前のスライス位置を求め、その求めた被検体が動く前のスライス位置を撮像することを特徴とする。
本発明のMRI装置によれば、検査中に被検体が動いて、撮像目的部位が設定済みのスライス位置からずれた場合でも、容易に位置修正が可能となる。そのため、従来必要であった、スカウト画像の撮り直しや、スライス位置の再設定を行うことなく、被検体が動く前のスライス位置で撮像できることになる。
以下、添付図面に従って本発明のMRI装置の好ましい実施形態について詳説する。なお、発明の実施形態を説明するための全図において、同一機能を有するものは同一符号を付け、その繰り返しの説明は省略する。
最初に、本発明に係るMRI装置の一例の全体概要を図1に基づいて説明する。図1は、本発明に係るMRI装置の一実施例の全体構成を示すブロック図である。このMRI装置は、NMR現象を利用して被検体101の断層画像を得るもので、図1に示すように、静磁場発生磁石102と、傾斜磁場コイル103及び傾斜磁場電源109と、送信RFコイル104及びRF送信部110と、受信RFコイル105及び信号検出部106と、信号処理部107と、計測制御部111と、全体制御部108と、表示・操作部113と、被検体101を搭載してその被検体101を静磁場発生磁石102の内部に出し入れするベッド112と、を備えて構成される。
静磁場発生磁石102は、垂直磁場方式であれば被検体101の体軸と直交する方向に、水平磁場方式であれば体軸方向に、それぞれ均一な静磁場を発生させるもので、被検体101の周りに永久磁石方式、常電導方式あるいは超電導方式の静磁場発生源が配置されている。
傾斜磁場コイル103は、MRI装置の座標系(静止座標系)であるX,Y,Zの3軸方向に巻かれたコイルであり、それぞれの傾斜磁場コイルは、それを駆動する傾斜磁場電源109に接続され電流が供給される。具体的には、各傾斜磁場コイルの傾斜磁場電源109は、それぞれ後述の計測制御部111からの命令に従って駆動されて、それぞれの傾斜磁場コイルに電流を供給する。これにより、X,Y,Zの3軸方向に傾斜磁場Gx,Gy,Gzが発生する。撮像時には、スライス面(撮像断面)に直交する方向にスライス傾斜磁場パルス(Gs)が印加されて被検体101に対するスライス面が設定され、そのスライス面に直交して且つ互いに直交する残りの2つの方向に位相エンコード傾斜磁場パルス(Gp)と周波数エンコード傾斜磁場パルス(Gf)が印加されて、エコー信号にそれぞれの方向の位置情報がエンコードされる。
送信RFコイル104は、被検体101に高周波磁場(以下、RFという)パルスを照射するコイルであり、RF送信部110に接続され高周波パルス電流が供給される。これにより、被検体101の生体組織を構成する原子の原子核スピンに核磁気共鳴が誘起される。具体的には、RF送信部110が、後述の計測制御部111からの命令に従って駆動されて、高周波パルスを振幅変調し、増幅した後に被検体101に近接して配置された送信RFコイル104に供給することにより、RFパルスが被検体101に照射される。
受信RFコイル105は、被検体101の生体組織を構成する原子核スピンのNMR現象により放出されるエコー信号(NMR信号)を受信するコイルであり、信号検出部106に接続されて受信したエコー信号を信号検出部106に送る。信号検出部106は、受信RFコイル105で受信したエコー信号の検出処理を行う。具体的には、RF送信コイル104から照射されたRFパルスによって誘起された被検体101の応答のエコー信号が被検体101に近接して配置された受信RFコイル105で受信され、後述の計測制御部111からの命令に従って、信号検出部106が、受信したエコー信号を増幅し、直交位相検波により直交する二系統の信号に分割し、それぞれを所定数(例えば128,256,512等)サンプリングし、各サンプリング信号をA/D変換してディジタル量に変換し、後述の信号処理部107に送る。 従って、エコー信号は所定数のサンプリングデータからなる時系列のデジタルデータ(以下、エコーデータという)として得られる。
計測制御部111は、被検体101の断層画像の再構成に必要なデータ収集のための種々の命令を、主に、傾斜磁場電源109と、RF送信部110と、信号検出部106に送信してこれらを制御する制御部である。具体的には、計測制御部111は、後述する全体制御部108の制御で動作し、ある所定のパルスシーケンスに基づいて、傾斜磁場電源109、RF送信部110及び信号検出部106を制御して、被検体101へのRFパルスと傾斜磁場パルスの印加及び被検体101からのエコー信号の検出を繰り返し実行し、被検体101の断層画像の再構成に必要なエコーデータを収集する。
全体制御部108は、計測制御部111の制御、及び、各種データ処理と処理結果の表示及び保存等の制御を行うものであってCPU及びメモリを内部に有する演算処理部と、光ディスク、磁気ディスク等の記憶部とを有して成る。具体的には、計測制御部111を制御してエコーデータの収集を実行させ、信号処理部107からのエコーデータが入力されると、演算処理部が信号処理、フーリエ変換による画像再構成等の処理を実行し、その結果である被検体101の断層画像を、後述の表示・操作部108に表示させると共に記憶部に記録する。
表示・操作部113は、被検体101の断層画像を表示するディスプレイと、MRI装置の各種制御情報や上記全体制御部108で行う処理の制御情報を入力するトラックボール又はマウス及びキーボード等の操作部と、から成る。この操作部はディスプレイに近接して配置され、操作者がディスプレイを見ながら操作部を通してインタラクティブにMRI装置の各種処理を制御する。
なお、図1において、送信側のRF送信コイル104と傾斜磁場コイル103は、被検体101が挿入される静磁場発生磁石102の静磁場空間内に、垂直磁場方式であれば被検体101に対向して、水平磁場方式であれば被検体101を取り囲むようにして設置されている。また、受信側の受信RFコイル105は、被検体101に対向して、或いは取り囲むように設置されている。
現在のMRI装置の撮像対象核種は、臨床で普及しているものとしては、被検体の主たる構成物質である水素原子核(プロトン)である。プロトン密度の空間分布や、励起状態の緩和時間の空間分布に関する情報を画像化することで、人体頭部、腹部、四肢等の形態または、機能を2次元もしくは3次元的に撮像する。
上記のようなMRI装置における典型的な撮像は計測制御部111の制御の下に次の様に行われる。即ち、所定のパルスシーケンスに基づいて、被検体101に、スライス傾斜磁場パルス(Gs)が印加された状態でRFパルスが照射されることにより所望のスライス領域にNMR現象が誘起される。次に、発生したエコー信号に位相エンコード傾斜射場(Gp)と周波数エンコード傾斜磁場(Gf)が印加されて、RF受信コイル105で受信され、信号検出部106で受信されたエコー信号が検出されてエコーデータとされ、そのエコーデータが演算処理部のメモリに記憶される。この操作が、位相エンコード傾斜磁場の印加量(傾斜磁場波形と時間軸との囲む面積)を変えて、所定回数(例えば、128,256,512等)繰り返される。これにより2次元断層像の画像再構成に必要な所定数(例えば、128X128, 128X256, 256X256,256X512, 512X512等)のデータが取得される。
(第1の実施形態)
次に、本発明のMRI装置の第1の実施形態を説明する。本実施形態は、被検体が動いた後に撮像された画像上で設定されたROI内データと3次元画像データとのマッチングを行い、被検体が動く前のスライス位置を求め、その求めたスライス位置を撮像する。以下、図2,3に基づいて本実施形態を説明する。
図2は、3つのスカウト画像(スカウトAX201,スカウトSAG202,スカウトCOL203)を用いて位置決めを行い、診断用としてT1強調画像(211)、T2強調画像(212)、FLAIR(213)の順で撮像する場合を示す。そして、T1強調画像とT2強調画像の撮像完了した時点で被検体が動き、最後のFLAIR画像213は、スカウト画像上で設定したスライス位置で撮像できなかった場合を示している。
このような場合であっても本実施形態では、FLAIR画像213のスライス位置を、T1強調画像211、T2強調画像212のスライス位置に合わせるために、スカウト画像の撮り直し、および、再位置決めは行わない。その代わりに、被検体が動いた後に取得された画像(この例では、FLAIR画像213)と、予め取得された3次元画像データ(以下、3Dデータ)301を用いて、被検体が動いた量を算出し、被検体の動き分だけスライス位置を調整し、調整されたスライス位置の画像(この例ではFLAIR画像303)を撮像する。
以下に、スライス位置調整方法の具体的処理内容を図3に示すフローチャートに基づいて説明する。本処理はプログラムとして記憶部に記憶されており、必要に応じて演算処理部のCPUがメモリに読み出して実行することにより実施される。
ステップ311で、操作者は、被検体が動いたと判断される画像213を選択する。操作者は、撮像された全ての画像(211〜213)を見比べて、被検体が動いて撮像された画像213(この例ではFLAIR画像)を選択する。
ステップ312で、操作者は、表示・操作部113のトラックボール又はマウスを操作して、被検体が動いたと判断した画像213上でその特徴部位(例えば、頭部であれば脳室など)を内部に含む様にROI302を設定する。
ステップ313で、CPUは、ステップ312で設定されたROI302内の画像パターンを抽出し、この抽出した画像パターンを、被検体が動く前に取得された3Dデータ301内で検索し、3Dデータ301内で一致するスライス位置を検索する。なお、ここで使用する3Dデータ301は、RF受信コイルの感度補正やパラレルイメージングの展開処理に使用する等の目的で、あらかじめ取得されたものを用いる。そのため、スライス位置調整のために新たな3Dデータを取得する必要は無い。
ステップ314で、ステップ313における3Dデータ301内でROI302内の画像パターンの検索が行われた結果、一致するスライス位置が見つかれば、CPUは、被検体が動いたことによる、最初に設定されたスライス位置の位置ずれ量、つまり被検体が動いたことによる、一連の診断画像の撮像前に設定されたスライス位置の移動量及び回転量を算出する。
ステップ315で、CPUは、ステップ314で取得したスライス位置ずれ量に基づいて、被検体が動く前の位置を再撮像するように、最初に設定されたスライス位置を調整する。具体的には、被検体が動いたことによる最初に設定されたスライス位置の移動量に基づいて、調整後のスライス位置の静磁場中心の法線ベクトル成分、法線距離及びFOVのOff-Center量を求め、被検体が動いたことによる最初に設定されたスライス位置の回転量から調整後のスライス位置の法線の(静止座標系における)方向ベクトルを求める。
ステップ316で、CPUは、計測制御部111に、スライス位置調整した位置のスライスを、被検体が動いた画像と同じコントラスト種となる撮像条件での撮像を指示し、計測制御部111は指示された撮像条件でスライス位置調整されたスライスを撮像する。具体的には、CPUは、調整後のスライス位置の静磁場中心の法線距離に応じて照射パルスの周波数を設定し、調整後のスライス位置の静磁場中心からの法線距離に応じてスライス位置のオブリーク角度を求め、このオリーブ角度に応じて、各方向の傾斜磁場(Gx,Gy,Gz)を求め、調整後のスライス位置のFOVのOff-Center量に応じて、エコー信号の検波周波数及び位相エンコード方向の位相回転量を求め、スライス面の方向ベクトル(Row,Col)からスライスの面内回転量を求める。そして、CPUは、求めたこれらの値を計測制御部111に指示して、調整されたスライス位置を撮像させる。これにより、被検体が動く前のT1強調画像211、T2強調画像212と同位置のFLAIR画像303が得られる。
以上までが、スライス位置調整方法の具体的処理フローの概要である。以上の処理により、調整後のスライス位置で取得されたFLAIR画像303は、被検体が動く前に取得されたT1強調画像211、T2強調画像212と、撮像対象のスライス位置及び表示された画像の位置関係とが一致することになる。
次に、ROI302内の画像パターンを3Dデータ301内で検索するステップ313の詳細を図4に示すフローチャートに基づいて以下説明する。
ステップ421で、前述のステップ312で、被検体が動いた後のFLAIR画像213上にROI302が設定されると、CPUは、被検体が動く前に取得された3Dデータ301内で、ROI302と同じ位置情報を持つ位置にROI412を設定する。ここで、FLAIR画像213のスライス位置は、スライス面内の1点(Pos1)と、2方向の単位ベクトル(Row1,Col1)で定義できる。Pos1、Row1、Col1は、撮像時に既に指定済みである。また、ROI302はFLAIR画像213のスライス面内にあるため、ROI302内の1点(Pos1')とROI302のサイズ(FOVr1,FOVr2)を指定すれば、FLAIR画像213と同じ2つの単位ベクトル(Row1,Col1)を用いてROI302を定義できる。したがって、これらのPos1',FOVr1,FOVr2、及びRow1,Col1とを用いて、FLAIR画像213上で設定されたROI302と同じROI412及びそのROI412を含むスライス位置411を3Dデータ303内に容易に設定可能となる。Pos1'、FOVr1、FOVr2は、操作者がディスプレイに表示されたUI上でマウスやキーボード入力によって設定する。この状態では、ROI302とROI412の位置情報は一致しているが、被検体の位置が異なるため、両ROI内の画像は異なっている。
ステップ422で、CPUは、被検体が動く前の画像であるROI412を3Dデータ301内で平行移動と回転移動させて移動させながら、移動後のROI412内画像パターンと被検体が動いた後の画像データであるROI302内画像パターンとのマッチングを繰り返し行い、一致するスライス位置を検索する。ROI内画像パターンの一致判定には、相関計算や残差計算(いずれも詳細は後述する)などが用いられる。一致判定の計算は、ROI412の3Dデータ301内での平行移動と回転移動の度に実行され、一致度合いが高いと判定される方向に移動方向を絞り込みながらマッチングが繰り返される。
ステップ423で、CPUは、相関計算や残差計算により、ROI302とROI412の画像パターンが一致するROI413が見つかったら、そのROI413を含むスライス位置403を調整後のスライス位置とする。
以上までがROI302内の画像パターンを3Dデータ301内で検索するステップ313の詳細である。
次に、前述のステップ422における相関計算の処理フローの一例を図5に示すフローチャートに基づいて説明する。
ステップ520で、CPUは、図4に示したROI302と3Dデータ301内のROI412それぞれで、x方向、もしくは、y方向の信号値プロファイルを取得する。図5の501,511は、それぞれROI302とROI412のy方向に信号値プロファイルを取得する位置を示す。また、501,511の信号値プロファイルをそれぞれ502,512とする。
ステップ521で、CPUは、得られた2つの信号値プロファイルを、(y1,y2)={(y1i,y2i)}(i=1,2,・・・n、nはライン数)とすると、2つの信号値プロファイルの相関係数を(1)式のように計算する。

Figure 0005437610

Figure 0005437610

ステップ522で、CPUは、隣接するラインでも同様の処理を行い、ROI内の全ラインで相関係数を算出する。図5に示す例では、ROI内の全x方向ラインについて、ステップ520〜521の処理が繰り返される。
ステップ523で、CPUは、全ライン分(nライン分)の相関係数が算出されたら、それらの平均値を求める。
ステップ524で、CPUは、ステップ523で求めた平均値と、検索用画像(3Dデータからの切り出し画像)のスライス位置情報と、を共に記憶部に保存する。
ステップ525で、CPUは、ステップ520〜524を、ステップ527で検索用画像のスライス位置を変更して、繰り返し実行し、3Dデータ301内の全ての位置を探索する。
ステップ526で、CPUは、全検索が終了したら、ステップ524にて保存した相関計数の平均値の中で、最大値となるスライス位置を、両ROI内の画像パターンが、最も一致したスライス位置と判定する。
以上までが、相関計算の処理フローの一例の概要である。なお、本例は、y方向の信号値プロファイル間の相関係数演算をx方向に繰り返す例をであったが、x方向の信号値プロファイル間の相関係数演算をy方向に繰り返してもよい。
次に、前述のステップ422における残差計算の処理フローの一例を図6に示すフローチャートに基づいて説明する。
ステップ601で、CPUは、図4に示したROI302と、3Dデータ301内のROI412の画素値を、それぞれ(x1i,y1j)、(x2i,y2j) (i=1,2,・・・,n:nはx方向のROI内画素数、j=1,2,・・・,m:mはy方向のROI内画素数)とすると、(x1i-x2i,y1j-y2j)で表すことができる差分画像を作成する。
ステップ602で、CPUは、ステップ601で作成した差分画像の全画素の平均画素値を算出する。
ステップ603で、CPUは、ステップ602で算出した平均画素値を、検索用画像のスライス位置と共に記憶部に保存する。
ステップ604で、CPUは、ステップ601〜603の処理を、ステップ606で検索用画像のスライス位置を変更しながら繰り返し実行し、3Dデータ301内の全ての位置を探索する。
ステップ605で、全検索が終了したら、CPUは、ステップ603にて保存した平均画素値の中で、最小値となるスライス位置を、両ROI内の画像パターンが、最も一致したスライス位置と判定する。
以上までが、残差計算の処理フローの一例の概要である。
なお、前述の相関計算と残差計算は、ROI内画像の画素値をそのまま用いて計算を行う例であったが、比較する画像が、必ずしも同じコントラスト画像とは限らない。コントラストの違いが計算結果に影響を及ぼさないよう、相関係数や残差を計算する前に、画像データに二値化などの前処理を入れることも可能である。二値化処理は、例えば一般的な固定閾値による二値化などを用いることができる。閾値(t)には、画像内の全画素の平均画素値やヒストグラムの谷となる値などを用いることができる。二値化の一例として、i軸、j軸からなる空間内で、ある画素の画素値をf(i,j)とした場合、二値化後の画素値f'(i,j)は、(2)式で表される。

Figure 0005437610

以上までが本発明のMRI装置の第1の実施形態の説明である。本実施形態のMRI装置によれば、ROIで設定された画像の特徴的部分で3Dデータ内を検索するので、検索精度を向上させることができる。そのため、スライス位置調整の精度を向上させることができ、検査中の被検体の動きによらずに、同じ撮像部位についての複数のコントラスト画像を撮像することが可能になる。その結果、従来必要であった、スカウト画像の撮り直しや、スライス位置の再設定を行うことなく、被検体が動く前のスライス位置で撮像できることになる。
(第2の実施形態)
次に、本発明のMRI装置の第2の実施形態を説明する。前述の第1の実施形態は、被検体が動いた後に撮像された画像に設定されたROIの画像パターンを用いて、被検体が動く前のスライス位置検索を行う例であったが、本実施形態は、ROIの設定を行うことなく、被検体が動いた後に撮像された画像を用いて被検体が動く前のスライス位置検索を行う。以下、前述の第1の実施形態と異なる点のみを図7に示す本実施形態の処理フローを用いて説明し、同一の点についての説明は省略する。(a)は残差計算を用いる場合の処理フローを、(b)は相関計算を用いる場合の処理フローを、それぞれ示す。
本実施形態でも、前述の第1の実施形態と同様に、スライス位置の検索は、図4に示す被検体が動いた後の画像213と3Dデータ301からの切り出し画像411とで行う。また、本実施形態では、操作者によるROI設定が無い。つまり、画像パターンのマッチング処理を行いやすい特徴的な領域としてのROIの、操作者による意図的な選択を行わない。以下、図7の処理フローを各処理ステップ毎に説明する。
ステップ700で、CPUは、必要に応じて、後述するマッチング処理の精度を向上させるために、例えば、予め画像213と411に輪郭強調処理710を行い、両画像内の被検体形状を認識しやすい状態にする。輪郭強調の必要がなければ、このステップを省略しても良い。
輪郭強調は、一般的なフィルタ処理によって行うことができる。一例として、2×2マスクの微分フィルタの場合を示す。i軸、j軸からなる空間内で、ある画素の画素値をf(i,j)とした場合、i、j両方向のフィルタ処理後の画素値gi(i,j)、gj(i,j)は(3)式で表すことができる。

Figure 0005437610

gi(i,j)、gj(i,j)を用いて、フィルタ処理後の新たな画素値は(4)式で定義する。
Figure 0005437610
ステップ701で、CPUは、ステップ700で輪郭強調処理した画像213と411を用いて、前述の第1の実施形態と同様に、(a)図6の残差計算によるマッチング、又は、(b)図5の相関計算によるマッチングなどを行い、3Dデータ301内で被検体が動いた後に撮像されたFLAIR画像と同じスライス位置を検索する。
ステップ702で、CPUは、図3に示す前述の第1の実施形態のステップ314以降の処理を行い、スライス位置調整した位置のスライス位置の再撮像を行う。
以上までが本実施形態の処理フローの概要である。本実施形態のMRI装置によれば、前述の第1の実施形態と比較して、ROIを設定する処理が無いので、スライス位置調整が簡略化される。そのため、操作者の負担を軽減しつつ、検査中の被検体の動きによらずに、同じ撮像部位についての複数のコントラスト画像を撮像することが可能になる。
(第3の実施形態)
次に、本発明のMRI装置の第3の実施形態を説明する。前述の第1,2の実施形態は、被検体が動いた後の診断用画像(これまでの例ではFLAIR画像)を見て、被検体が動いたことを認識し、スライス位置調整、および、再撮像を行う場合であった。しかし、診断用画像を撮像する前に、インルームモニタ等で被検者が動いたことを認知できた場合、診断用画像を撮像する必要は無い。そこで本実施形態は、撮像時間の長い診断用画像の代わりに、数秒で撮像可能なスカウト画像を撮像し、このスカウト画像と予め取得しておいた3Dデータとのマッチングを行う。以下、前述の第1、2の実施形態と異なる点のみを図8,9を用いて説明し、同一の点についての説明は省略する。
図8は、3つのスカウト画像(スカウトAx801、スカウトSAG802、スカウトCOL803)を用いてスライス位置調整し、T1強調画像811、T2強調画像812まで撮像が完了した状態し、この時点で、インルームモニタ等で被検体が動いたことを、操作者が認識したとする。そのまま撮像すると、T2強調画像の次の画像(例えばFLAIR画像)は、異なる部位が写った画像となる。そこで本実施形態は、被検体が動いた後には診断画像を撮像せずにスカウト画像を撮像して、このスカウト画像を用いてスライス位置調整を行い、調整されたスライス位置の撮像を行う。以下、図9に示すフローチャートを用いて本実施形態の処理フローにおける各処理ステップを詳細に説明する。
ステップ911で、操作者は、インルームモニタ等で被検体が動いたことを認識する。
ステップ912で、操作者は一連の診断画像の撮像を中止し、スカウト画像813の撮像を行うため、表示・操作部113のキーボード等を介してスカウト画像813の撮像を起動する。そこで、CPUは、計測制御部111に指示して、スカウト画像813を撮像させる。第1,2の実施形態の説明で示したように、スライス位置検索には、相関計算や残差計算を用いるため、比較する画像は、同じコントラスト画像であることが望ましい。そこで、スライス位置調整用のスカウト画像813は、比較対照である予め取得された3Dデータ301と同じコントラスト画像となるような撮像条件で撮像される。
ステップ913で、CPUは、スライス位置調整用のスカウト画像813が取得されたら、3Dデータ301内でスライス位置検索を行う。このとき、前述の第1の実施形態で説明したように、特徴部位にROI902を設定して、ROI内画像パターンを用いてスライス位置検索することも可能である。また、前述の第2の実施形態で説明したように、ROIを設定せずに、スカウト画像を用いたスライス位置検索も可能である。
ステップ914で、CPUは、ステップ912のスライス位置検索の結果に基づいて、被検体の位置ずれ量を算出する。
ステップ915で、CPUは、ステップ914で算出された被検体の位置ずれ量に基づいて、図3に示す前述の第1の実施形態の処理ステップ315及び316の各処理を行う。これにより、被検体が動く前のT1強調811、T2強調画像812と同位置のFLAIR画像903が得られる。
以上までが本実施形態の処理フローの概要である。
以上説明したように、本実施形態のMRI装置によれば、前述の第1、2の実施形態と比較して、被検体が動いた後の診断画像を撮像することなく、診断画像の撮像時間より短い時間でスカウト画像を撮像し、このスカウト画像を用いてスライス位置調整を行うので、検査中の被検体の動きによらずに、同じ撮像部位についての複数のコントラスト画像を、より短時間で、撮像することが可能になる。
(第4の実施形態)
次に、本発明のMRI装置の第4の実施形態を説明する。前述の各実施形態は、自動的に被検体が動く前のスライス位置に合わせて撮像を実行したが、本実施形態は、スライス位置検索課程で、3Dデータ301から切り出された複数のスライス位置の画像をディスプレイ上にそれぞれ表示し、それらの画像の内から選択された画像と同じスライス位置を撮像する。以下、前述の各実施形態と異なる点のみを図10を用いて説明し、同一の点についての説明は省略する。
検査用画像の撮像が終了すると、取得された各コントラスト画像(T1強調画像(T1W)1011,T2強調画像(T2W)1012,FLAIR画像1013)がディスプレイ1001上の、結果画像表示エリア1002に表示される。これらの画像を操作者が確認し、スライス位置を調整して撮り直したい画像を選択する(本例では、1013のFLAIR画像が選択される)。操作者が画像を選択すると、CPUは、前述の第1,2実施形態と同様に、3Dデータ301内で一致するスライスを検索する。その際、CPUは、3Dデータ301から切り出されて、マッチング処理に用いられた各画像をディスプレイ上のマッチング画像表示エリア1003に表示させる。画像表示エリア1003内には、マッチング処理の結果、一致している確率が高い画像を中心にして各画像が表示される。操作者は、マッチング画像表示エリア1003内の画像から、FLAIR画像1013と最も一致している切り出し画像1004を選択し、再撮像スタートボタン1005を押す。これにより、CPUは、図3に示す前述の第1の実施形態におけるステップ314以降の処理を行うことにより、計測制御部111に調整されたスライス位置を撮像させる。撮像が完了すると、結果表示エリア1006に撮り直した画像(この場合はFLAIR画像)が表示される。
以上までが本実施形態の処理フローの概要である。本実施形態のMRI装置によれば、操作者により、3Dデータからの切り出し画像の内から、被検体が動いた後に撮像された画像と一致する画像が選択されるので、被検体が動く前のスライス位置検索を確実に行うことが可能になる。その結果、検査中の被検体の動きによらずに、より高い精度で、同じ撮像部位についての複数のコントラスト画像を撮像することが可能になる。
以上までが、本発明のMRI装置における、検査中の被検体の動きによらずに、同じ撮像部位についての複数のコントラスト画像を撮像する各実施形態の説明である。しかし、本発明のMRI装置は、上記各実施形態の説明で開示された内容にとどまらず、本発明の趣旨を踏まえた上で他の形態を取り得る。
例えば、被検体の動きの検出は、インルームモニタに限らず、周知のナビゲータエコーや、体動検出センサを被検体の腹部等に配置してもよい。
本発明に係るMRI装置の一例の全体概要を示すブロック図。 本発明の第1の実施形態に係り、T1強調画像とT2強調画像の撮像完了した時点で被検体が動き、最後のFLAIR画像は、スカウト画像上で設定したスライス位置で撮像できなかった場合を示す図。 本発明の第1の実施形態に係り、スライス位置調整方法の具体的処理内容を示すフルーチャート。 ROI302内の画像パターンを3Dデータ301内で検索するステップ313の詳細を示すフローチャート。 相関計算の処理フローの一例を示すフローチャート。 残差計算の処理フローの一例を示すフローチャート。 本発明の第2の実施形態に係る処理内容を示すフルーチャート。 本発明の第3の実施形態に係り、T1強調画像811、T2強調画像812まで撮像が完了した状態した時点で、被検体が動いたことを操作者が認識し、スカウト画像を撮像する場合を示す図。 本発明の第3の実施形態に係る処理内容を示すフルーチャート。 本発明の第4の実施形態に係るディスプレイ表示の一例を示す図。
符号の説明
101 被検体、102 静磁場発生磁石、103 傾斜磁場コイル、104 送信RFコイル、105 受信RFコイル、106 信号検出部106、107 信号処理部、108 全体制御部、109 傾斜磁場電源、110 RF送信部、111 計測制御部、112 ベッド、113 表示・操作部

Claims (4)

  1. 被検体を撮像して画像の再構成に必要なエコーデータを取得する計測制御部と、
    前記エコーデータを用いて前記被検体の画像を再構成する演算処理部と、
    前記画像を表示する表示部と、
    を備え
    前記演算処理部は、予め取得された3次元画像ら複数の画像を切り出し、該複数の画像の中から前記被検体が動いた後に撮像された画像とのマッチング処理をそれぞれ行い、
    前記表示部は、一致度が相対的に高い複数の画像をそれぞれ表示し、
    前記計測制御部は、前記一致度が相対的に高い複数の画像の内から選択された画像と同じスライス位置の撮像を行うことを特徴とする磁気共鳴イメージング装置。
  2. 被検体を撮像して画像の再構成に必要なエコーデータを取得する計測制御部と、
    前記エコーデータを用いて前記被検体の画像を再構成する演算処理部と、
    前記画像を表示する表示部と、
    を備え、
    前記演算処理部は、予め取得された3次元画像データと、前記被検体が動いた後に撮像された画像上でその特徴部位を含むように設定されたROI内のデータと、のマッチング処理を行い、前記被検体が動く前のスライス位置を求め、
    前記計測制御部は、前記被検体が動く前のスライス位置を撮像することを特徴とする磁気共鳴イメージング装置。
  3. 請求項2記載の磁気共鳴イメージング装置において、
    前記被検体の動きを検出する動き検出部を備え、
    前記被検体が動いた後に撮像される画像は、前記3次元画像と略同一コントラストのスカウト画像であり、
    前記計測制御部は、前記動き検出部により前記被検体の動きが検出されると、前記スカウト画像の撮像を行うことを特徴とする磁気共鳴イメージング装置。
  4. 請求項1乃至3のいずれか一項に記載の磁気共鳴イメージング装置において、
    前記演算処理部は、前記被検体が動いた後に撮像された画像をフィルタ処理し、該フィルタ処理された画像のデータと前記3次元画像のデータとのマッチング処理を行い、前記被検体が動く前のスライス位置を求めることを特徴とする磁気共鳴イメージング装置。
JP2008223151A 2008-09-01 2008-09-01 磁気共鳴イメージング装置 Active JP5437610B2 (ja)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2008223151A JP5437610B2 (ja) 2008-09-01 2008-09-01 磁気共鳴イメージング装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2008223151A JP5437610B2 (ja) 2008-09-01 2008-09-01 磁気共鳴イメージング装置

Publications (2)

Publication Number Publication Date
JP2010057532A JP2010057532A (ja) 2010-03-18
JP5437610B2 true JP5437610B2 (ja) 2014-03-12

Family

ID=42185046

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2008223151A Active JP5437610B2 (ja) 2008-09-01 2008-09-01 磁気共鳴イメージング装置

Country Status (1)

Country Link
JP (1) JP5437610B2 (ja)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP5451010B2 (ja) * 2008-08-29 2014-03-26 キヤノン株式会社 X線検出装置、情報処理方法および記録媒体
JP2012115399A (ja) * 2010-11-30 2012-06-21 Ge Medical Systems Global Technology Co Llc 医用装置およびプログラム
JP5745444B2 (ja) 2012-03-05 2015-07-08 富士フイルム株式会社 医用画像表示装置および医用画像表示方法、並びに、医用画像表示プログラム
US20180045800A1 (en) * 2015-02-24 2018-02-15 Koninklijke Philips N.V. Scan geometry planning method for mri or ct

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP3498980B2 (ja) * 1992-09-16 2004-02-23 東芝医用システムエンジニアリング株式会社 磁気共鳴イメージング用スキャン位置決め方法及び磁気共鳴イメージング装置
JPH07299061A (ja) * 1994-04-30 1995-11-14 Shimadzu Corp 断層撮影装置
JP2001157667A (ja) * 1999-12-02 2001-06-12 Fuji Photo Film Co Ltd 画像表示方法および画像表示装置
US7206460B2 (en) * 2001-11-01 2007-04-17 General Electric Company Method for contrast matching of multiple images of the same object or scene to a common reference image

Also Published As

Publication number Publication date
JP2010057532A (ja) 2010-03-18

Similar Documents

Publication Publication Date Title
US9254112B2 (en) Respiratory interval-based correlation and processing of dynamic imaging data
US8502532B2 (en) Magnetic resonance data acquisition system and method with recursively adapted object-specific measurement parameter adjustment during patient movement through the MRI system
JP6141883B2 (ja) ディクソン法を使って取得されるナビゲーターを使う動き補正のあるmri
CN106233154B (zh) 使用预脉冲和导航器的具有运动校正的磁共振成像
US9626777B2 (en) Method and apparatus to generate image data
US9766316B2 (en) Magnetic resonance imaging device and quantitative susceptibility mapping method
KR101652387B1 (ko) 영상 데이터를 생성하는 방법
US9952303B2 (en) Method for automatic calibration of motion detection techniques in medical imaging systems
US8781552B2 (en) Localization of aorta and left atrium from magnetic resonance imaging
JP6549612B2 (ja) 改善された多相動的造影磁気共鳴撮像方法
US7457657B2 (en) Magnetic resonance method and apparatus for determining the position and/or orientation of the image plane of slice image exposures of a vessel region in a contrast agent bolus examination
EP3595525B1 (en) Systems and methods for real-time motion prediction in dynamic imaging
US10018699B2 (en) Method and magnetic resonance apparatus for acquiring magnetic resonance data with a prospective motion correction
US10641854B2 (en) Systems and methods for automatic voxel positioning in magnetic resonance spectroscopy
JP5437610B2 (ja) 磁気共鳴イメージング装置
JP2009005789A (ja) 核磁気共鳴装置
US10330764B2 (en) Magnetic resonance measuring apparatus and image processing method
US8995738B2 (en) System and method for magnetic resonance imaging parametric mapping using confidence maps
US10459055B2 (en) System and method for reduced field of view MR fingerprinting for parametric mapping
WO2016021440A1 (ja) 磁気共鳴イメージング装置
US11971467B2 (en) Free-breathing abdominal magnetic resonance fingerprinting using a pilot tone navigator
JP2004305454A (ja) 磁気共鳴イメージング装置
US10859653B2 (en) Blind source separation in magnetic resonance fingerprinting

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20110826

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20130624

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20130820

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: 20131125

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20131212

R150 Certificate of patent or registration of utility model

Ref document number: 5437610

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

Free format text: JAPANESE INTERMEDIATE CODE: R150

S111 Request for change of ownership or part of ownership

Free format text: JAPANESE INTERMEDIATE CODE: R313111

S533 Written request for registration of change of name

Free format text: JAPANESE INTERMEDIATE CODE: R313533

R350 Written notification of registration of transfer

Free format text: JAPANESE INTERMEDIATE CODE: R350

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

S111 Request for change of ownership or part of ownership

Free format text: JAPANESE INTERMEDIATE CODE: R313111

R350 Written notification of registration of transfer

Free format text: JAPANESE INTERMEDIATE CODE: R350

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

S111 Request for change of ownership or part of ownership

Free format text: JAPANESE INTERMEDIATE CODE: R313111

R350 Written notification of registration of transfer

Free format text: JAPANESE INTERMEDIATE CODE: R350