(57)【特許請求の範囲】
【請求項１】
撮影領域の位置決め及び前記撮影領域内へのラベル領域の位置決めのために事前に取得された参照画像にラベル領域を表すマークを表示し、
第１のパルスシーケンスのもとで、スライス選択用傾斜磁場印加と共にインバージョンパルスを印加することによって、前記撮影領域内の前記ラベル領域中の移動体を含む撮影対象をラベルし、前記インバージョンパルスから待ち時間の後に、前記撮影領域から信号を収集するための走査を実行することによって前記撮影対象がラベルされた第１の画像を取得し、
第２のパルスシーケンスのもとで、前記インバージョンパルスを印加しないで前記撮影領域から信号を収集するための走査を実行することによって前記撮影対象がラベルされていない第２の画像を取得し、
前記第１、第２の画像を差分することにより、前記ラベルされた移動体を選択的に抽出した画像を生成することを特徴とする磁気共鳴映像化装置。
【請求項２】
前記第１、第２のパルスシーケンスの走査は交互に実行されることを特徴とする請求項１記載の磁気共鳴映像化装置。
【請求項３】
前記第１のパルスシーケンスは、前記スライス選択用傾斜磁場と共に印加される前記移動体をラベルするための前記インバージョンパルスの他に、前記ラベル領域以外の領域を低信号化するためにスライス選択用傾斜磁場が印加されない状態のもとで印加される高周波パルスを含むことを特徴とする請求項１記載の磁気共鳴映像化装置。
【請求項４】
前記ラベル領域を表すマークは操作部の操作に従って回転されることを特徴とする請求項１記載の磁気共鳴映像化装置。
【請求項５】
前記ラベル領域を表すマークは操作部の操作に従って形状変更されることを特徴とする請求項１記載の磁気共鳴映像化装置。
【請求項６】
前記ラベルされた移動体の移動経路に基づいて前記移動体の平均速度が計算されることを特徴とする請求項１記載の磁気共鳴映像化装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、ＭＲアンギオグラフィ（ＭＲＡ）機能を備えた磁気共鳴映像化装置に関する。
【背景技術】
【０００２】
磁気共鳴映像化は、静磁場中に置かれた被検体の原子核スピンをラーモア周波数の高周波信号で励起し、この励起に伴って発生するＦＩＤ（自由誘導減衰）信号やエコー信号から画像を得る手法である。
【０００３】
この磁気共鳴映像化の一つのカテゴリーとして、血液の動態を画像化するＭＲアンギオグラフィがある。さらにＭＲアンギオグラフィの一形態として、インバージョンパルス（高周波反転パルス）により局所的にスピンに標識（ラベル）を付けて、信号抜けさせる方法がある。
【０００４】
例えば、インバージョンパルスによりプリサチュレートさせる場所（ラベル領域）を、インバージョンパルスから本スキャンの励起パルスまでの時間間隔とともに設定し、プリサチュレーションの影響により信号変化（低信号化）する部分の移動軌跡を画像上で観察することにより、血液などの動きを観察するフローイメージング法（ボーラストラッキング法、Edelman, et.al., Radiology,Vol.171,551-556,1989）がある。
【０００５】
図９にこの方法によるパルスシーケンスを、図１０に従来のラベル領域の設定画面を示す。まず、インバージョンパルス１が印加される。このとき、関心のある領域ＲＯＩへの流入元と想定される領域（ラベル領域）が選択的に反転励起されるように、スライス選択傾斜磁場２が印加されるとともに、そのラベル領域の位置に対応するオフセット周波数３がインバージョンパルス１に付加されている。なお、ラベル領域は、例えば図９に点線で示した方向に印加することにより、、任意の方向に設定できる。これらの設定法については、Matt Bernstein et.al., JMRI, Vol.4, 105-108,1994等に示されている。
【０００６】
インバージョンパルス１から、待ち時間（反転時間）ＴＩ[s] 後に、励起パルス４により撮影断面を励起する。位相エンコード傾斜磁場５を印加したのち、エコー時間ＴＥ後にエコー信号６を収集する。この一連の操作を画像の再構成に必要なエンコードステップ数だけ繰り返すことによって１回の撮影が完了する。
【０００７】
エコー時間ＴＥが反転時間ＴＩにくらべ十分に短いとすると、血流が平均速度ｖ[m/s] で一定方向に移動するとき、ラベル領域７内でプレサチュレートされた血流は、ほぼｖ×ＴＩ[m] だけ移動した位置でとらえらえる。これにより、撮影対象の撮影中の動きを観察することができる。
【０００８】
この方法を診断に利用する場合、ラベル領域の場所の設定と、反転時間ＴＩの設定とが、最も重要な撮影パラメータであり、しかしそれは操作者にとって直感的でなく使いにくいという問題があった。
【０００９】
例えば、関心領域ＲＯＩ付近での血流の動態を観察したいとき、プリサチュレートされた血流が励起パルス４を印加する時点で、ちょうど当該関心領域ＲＯＩに到達しているように、当該関心領域に対して適切な距離だけ上流の位置にラベル領域を設定する必要がある。この距離は、血流速度ｖと、反転時間ＴＩとに依存して決まる。このため、目的の画像を収集するためには、反転時間ＴＩを様々に変えながら撮影したり、ラベル領域を移動させながら撮影を繰り返す必要があり、非常に不便である。特に、撮影対象のＴ１値等の事情により、反転時間ＴＩとして設定できる範囲は限られており、設定は困難な場合がある。
【発明の開示】
【発明が解決しようとする課題】
【００１０】
本発明の目的は、磁気共鳴映像化装置において、撮影対象として例えば血流に対して動きのない部分からの信号を抑制して、動きのある部分をより明瞭に描出することにある。
【課題を解決するための手段】
【００１１】
本発明は、撮影領域の位置決め及び前記撮影領域内へのラベル領域の位置決めのために事前に取得された参照画像にラベル領域を表すマークを表示し、第１のパルスシーケンスのもとで、スライス選択用傾斜磁場印加と共にインバージョンパルスを印加することによって、前記撮影領域内の前記ラベル領域中の移動体を含む撮影対象をラベルし、前記インバージョンパルスから待ち時間の後に、前記撮影領域から信号を収集するための走査を実行することによって前記撮影対象がラベルされた第１の画像を取得し、第２のパルスシーケンスのもとで、前記インバージョンパルスを印加しないで前記撮影領域から信号を収集するための走査を実行することによって前記撮影対象がラベルされていない第２の画像を取得し、前記第１、第２の画像を差分することにより、前記ラベルされた移動体を選択的に抽出した画像を生成する。
【発明の効果】
【００１２】
本発明によれば、撮影対象として例えば血流に対して動きのない部分からの信号を抑制して、動きのある部分をより明瞭に描出することができる。
【発明を実施するための最良の形態】
【００１３】
以下、図面を参照して本発明による磁気共鳴映像化装置（ＭＲＩ装置）を実施形態により説明する。このＭＲＩ装置は、特徴的には、造影剤を使用することなく、被検体内の動く対象物（移動体）としての例えば血流の動態を表す画像を取得する非造影ＭＲアンギオグラフィ（ＭＲＡ）を実行する機能を有している。この非造影ＭＲＡを行なうパルスシーケンスには、選択励起の高調波反転（インバージョン）パルスを用いる。
【００１４】
図１に本実施形態によるＭＲＩ装置の概略構成を示す。このＭＲＩ装置は、被検体としての患者Ｐを載せる寝台部と、静磁場を発生させる静磁場発生部と、静磁場に位置情報を付加するための傾斜磁場発生部と、高周波信号を送受信する送受信部と、システム全体のコントロール及び画像再構成を担う制御・演算部と、被検体Ｐの心時相を表す信号としてのＥＣＧ信号を計測する心電計測部とを備えている。
【００１５】
静磁場発生部は、例えば超電導方式の磁石１０１と、この磁石１０１に電流を供給する静磁場電源１０２とを備え、被検体Ｐが遊挿される円筒状の開口部（診断用空間）の軸方向（Ｚ軸方向）に静磁場Ｈ０ を発生させる。なお、この磁石部にはシムコイル１１４が設けられている。このシムコイル１１４には、ホスト計算機１０６の制御下で、シムコイル電源１１５から静磁場均一化のための電流が供給される。寝台部は、被検体Ｐを載せた天板を磁石１０１の開口部に退避可能に挿入できる。
【００１６】
傾斜磁場発生部は、磁石１０１に組み込まれた傾斜磁場コイルユニット１０３を備える。この傾斜磁場コイルユニット１０３は、互いに直交するＸ、Ｙ及びＺ軸方向の傾斜磁場を発生させるための３組（種類）のｘ，ｙ，ｚコイル１０３ｘ、１０３ｙ、１０３ｚを備える。傾斜磁場部はまた、ｘ，ｙ，ｚコイル１０３ｘ、１０３ｙ、１０３ｚに電流を供給する傾斜磁場電源４を備える。この傾斜磁場電源４は、後述するシーケンサ５の制御のもと、ｘ，ｙ，ｚコイル１０３ｘ、１０３ｙ、１０３ｚに傾斜磁場を発生させるためのパルス電流を供給する。
【００１７】
傾斜磁場電源１０４からｘ，ｙ，ｚコイル１０３ｘ、１０３ｙ、１０３ｚに供給されるパルス電流を制御することにより、物理軸である３軸Ｘ，Ｙ，Ｚ方向の傾斜磁場を合成して、互いに直交するスライス方向傾斜磁場Ｇｓ、位相エンコード方向傾斜磁場Ｇｅ、および読出し方向（周波数エンコード方向）傾斜磁場Ｇｒの各論理軸方向を任意に設定・変更することができる。スライス方向、位相エンコード方向、および読出し方向の各傾斜磁場は、静磁場Ｈ０に重畳される。
【００１８】
送受信部は、磁石１０１内の撮影空間にて被検体Ｐの近傍に配設されるＲＦコイル１０７と、このコイル１０７に接続された送信器１０８Ｔ及び受信器１０８Ｒとを備える。この送信器１０８Ｔ及び受信器１０８Ｒは、シーケンサ１０５の制御のもとで動作する。送信器１０８Ｔは、核磁気共鳴（ＮＭＲ）を起こさせるためのラーモア周波数のＲＦ電流パルスをＲＦコイル１０７に供給する。受信器１０８Ｒは、ＲＦコイル１０７が受信したエコー信号（高周波信号）を取り込み、これに前置増幅、中間周波変換、位相検波、低周波増幅、フィルタリングなどの各種の信号処理を施した後、Ａ／Ｄ変換してエコー信号に応じたデジタル量のエコーデータ（原データ）を生成する。
【００１９】
制御・演算部は、シーケンサ（シーケンスコントローラとも呼ばれる）１０５、ホスト計算機１０６、演算ユニット１１０、記憶ユニット１１１、表示器１１２、入力器１１３、音声発生器１１６、さらにＭＲＡエキスパートシステム１２０を備える。この内、ホスト計算機１０６は、記憶したソフトウエア手順により、シーケンサ１０５にパルスシーケンス情報を指令するとともに、装置全体の動作を統括する機能を有する。
【００２０】
ホスト計算機１０６は、位置決め用スキャンなどの準備作業に引き続いて、パルスシーケンスに基づいてイメージングスキャンを実施する。このイメージングスキャンは、画像再構成に必要なエコーデータの組を収集するスキャンであり、ここでは２次元スキャンに設定されている。イメージングスキャンは、ＥＣＧ信号に依るＥＣＧゲート法を併用して行われる。なお、このＥＣＧゲート法は場合によっては併用しなくてもよい。
【００２１】
このパルスシーケンスとしては、３次元（３Ｄ）スキャンまたは２次元（２Ｄ）スキャン）である。そのパルス列の形態としては、ＳＥ（スピンエコー）法、ＦＳＥ（高速ＳＥ）法、ＦＡＳＥ（高速 Ａｓｙｍｍｅｔｒｉｃ ＳＥ）法（すなわち、高速ＳＥ法にハーフフーリエ法を組み合わせたイメージング法）、ＥＰＩ（エコープラナーイメージング）法、などが用いられる。
【００２２】
シーケンサ１０５は、ＣＰＵおよびメモリを備えており、ホスト計算機１０６から送られてきたパルスシーケンス情報を記憶し、この情報にしたがって傾斜磁場電源１０４、送信器１０８Ｔ、受信器１０８Ｒの動作を制御するとともに、受信器１０８Ｒが出力したエコーデータを一旦入力し、これを演算ユニット１０１０に転送するように構成されている。ここで、パルスシーケンス情報とは、一連のパルスシーケンスにしたがって傾斜磁場電源１０４、送信器１０８Ｔおよび受信器１０８Ｒを動作させるために必要な全ての情報であり、例えばｘ，ｙ，ｚコイル１０３ｘ、１０３ｙ、１０３ｚに印加するパルス電流の強度、印加時間、印加タイミングなどに関する情報を含む。
【００２３】
また、演算ユニット１１０は、受信器１０８Ｒが出力したエコーデータ（原データ又は生データ）をシーケンサ１０５を通して入力し、その内部メモリ上のフーリエ空間（ｋ空間または周波数空間とも呼ばれる）にエコーデータを配置し、このエコーデータを各組毎に２次元又は３次元のフーリエ変換に付して実空間の画像データに再構成する。また演算ユニットは、必要に応じて、画像に関するデータの合成処理、差分演算処理などを行うことができる。
【００２４】
この合成処理には、２次元の複数フレームの画像データを対応する画素毎に加算する加算処理、３次元データに対して視線方向の最大値又は最小値を選択する最大値投影（ＭＩＰ）又は最小値（ＭＩＰ）投影処理などが含まれる。また、合成処理の別の例として、フーリエ空間上で複数フレームの軸の整合をとってエコーデータのまま１フレームのエコーデータに合成するようにしてもよい。なお、加算処理には、単純加算処理、加算平均処理、重み付け加算処理などが含まれる。
【００２５】
記憶ユニット１１１は、再構成された画像データのみならず、上述の合成処理や差分処理が施された画像データを保管することができる。表示器１１２は画像を表示する。また入力器１１３を介して、術者が希望する撮影条件、パルスシーケンス、画像合成や差分演算に関する情報をホスト計算機１０６に入力できる。
【００２６】
さらに、心電計測部は、被検体の体表に付着させてＥＣＧ信号を電気信号として検出するＥＣＧセンサ１１７と、このセンサ信号にデジタル化処理を含む各種の処理を施してホスト計算機１０６およびシーケンサ１０５に出力するＥＣＧユニット１１８とを備える。この心電計測部による計測信号は、イメージングスキャンを実行するときにシーケンサ１０５により用いられる。これにより、ＥＣＧゲート法（心電同期法）による同期タイミングを適切に設定でき、この同期タイミングに基づくＥＣＧゲート法のイメージングスキャンを行ってデータ収集できるようになっている。
【００２７】
ＭＲＡエキスパートシステム１２０は、ラベル領域の位置と、インバージョンパルスから励起パルスまでの待ち時間（反転時間）ＴＩとの設定を支援する機能を備えている。この支援機能としては、撮影断面の位置決め及びラベル領域のポジショニングのために予め取得された参照画像を、ラベル領域を表すマークと、インバージョンパルスでラベルされた血流がラベル領域から反転時間ＴＩを経て励起パルス印加時点に到達する位置を表す指標とともに表示器１１２に表示させるというもので、その支援モードとして、本実施形態では、２種類容易する。操作者は、第１と第２の２種類の支援モードを選択的に用いることができる。
【００２８】
第１の支援モードでは、ある反転時間ＴＩを想定し、その想定した反転時間ＴＩのもとで様々な血流速度ｖに対応する複数の到達位置が提示され、一方、第２の支援モードは、ある血流速度ｖを想定し、その想定した血流速度ｖのもとで様々な反転時間ＴＩに対応する複数の到達位置が提示される。
【００２９】
図２（ａ）に第１の支援モードにおけるラベル領域の設定支援画面の一例を示している。参照画像９は、本スキャン前に、スピンエコー又はエコープラナー法等を使って比較的低解像度仕様で取得された画像である。ラベル領域の設定シーンでは、ＭＲＡエキスパートシステム１２０により、この参照画像９上に、ラベル領域を表す矩形のマーク１０が、血流がラベル領域から、想定した反転時間ＴＩを経て、様々な血流速度ｖのもとで到達する複数の想定位置を点線で表す複数のラインマーク１１と、それぞれ対応する想定流速Ｖと、想定した反転時間ＴＩとともに表示される。図２では、３本のラインマーク１１が、それぞれ対応する流速１０ｃｍ／ｓ、２０ｃｍ／ｓ、３０ｃｍ／ｓの付箋とともに表示される。ラインマーク１１とそれぞれ対応する流速との関係は、任意に設定可能であり、操作者が数値設定するようにしてもよいし、検査部位に応じてデフォルト値として設定されるようにしてもよい。
【００３０】
各ラインマーク１１は、それぞれ対応する流速をｖとして、それらに、想定した反転時間ＴＩをかけた距離の画像縮尺換算距離を、ラベル領域マーク１０の端（又は中央）から、離間させた位置であって、ラベル領域マーク１０の長軸と平行に配置される。
【００３１】
操作者は、参照画像９で関心領域ＲＯＩを確認し、そのＲＯＩの位置に対してラインマーク１１を参考にして、関心領域ＲＯＩに対して好適な位置にラベル領域をそのマーク１０により設定でき、またラベル領域を固定した状態で、想定した反転時間ＴＩを変更することで、反転時間ＴＩを好適な時間に設定することができる。具体的には、操作者は、入力器１１３の例えばマウス等のポインタを操作して、参照画像９上でラベル領域のマーク１０を移動し、また必要に応じてその形状を変化させる。また、操作者は、入力器１１３の例えばキーボードをたたいて想定反転時間ＴＩを変更する。
【００３２】
ラベル領域のマーク１０が移動されたとき、ＭＲＡエキスパートシステム１２０は、参照画像９上のラインマーク１１の位置を再計算し、その位置にラインマーク１１を再表示させる。それによりラベル領域のマーク１０の移動に追従して、ラインマーク１１の表示位置が移動する。
【００３３】
想定反転時間ＴＩが変更されたとき、ＭＲＡエキスパートシステム１２０は、その変更された反転時間ＴＩに従って、参照画像９上のラインマーク１１の位置を再計算し、その位置にラインマーク１１を再表示させる（図２（ｂ）参照）。
【００３４】
さらに、操作者は、ラベル領域のマーク１０を回転操作することが可能である。図３には、回転されたラベル領域のマーク１３を示している。ラベル領域のマーク１３の長軸との平行を維持するように、ＭＲＡエキスパートシステム１２０は、このラベル領域のマーク１３の回転に追従して、ラインマーク１１を傾斜させる。撮影対象はラベル領域からどの方向にも移動する可能性があるので、ラベル領域のマーク１３の長軸と平行に、上下どちらもラインマーク１４を点線で表示する。
【００３５】
また、操作者は、ラベル領域のマーク１０の形状を、円形に変更することが可能である。図４には、円形に変更されたラベル領域のマーク１９を示している。ラベル領域のマーク１３と同心円を描くように、ＭＲＡエキスパートシステム１２０は、ラベル血流の到達位置を表すマーク２０を円形に変更する。ラベル領域のマーク１３の中心又は辺縁から、各マーク２０までの距離は、矩形の場合と同様に、ＴＩ×ｖに相当する距離に設定される。なお、このような矩形でない領域を反転励起する方法に関しては、P.A.Bottomley, et.al., J.of Appl. Phys. Vol62, 4284, 1987 あるいは、C.J. Hardy ,et.al., J.of Magn. Reso.,Vol77, 233-250, 1988に示されている。
【００３６】
次に、第２の支援モードについて説明する。上述の第１の支援モードでは、反転時間ＴＩを想定し、それのもとで様々な血流速度ｖに対応する複数の到達位置が提示されていたが、この第２の支援モードは、ある血流速度ｖを想定し、それのもとでの様々な反転時間ＴＩに対応する複数の到達位置を提示する。
【００３７】
図５（ａ）に第２の支援モードにおけるラベル領域の設定支援画面の一例を示している。ＭＲＡエキスパートシステム１２０により、参照画像９上に、ラベル領域を表す矩形のマーク１５が、血流がラベル領域から、想定した流速で流れ出て、様々な反転時間ＴＩを経て到達する複数の想定位置を点線で表す複数のラインマーク１６と、それぞれ対応する反転時間ＴＩと、想定した流速Ｖともに表示される。
【００３８】
図５（ａ）では、３本のラインマーク１６が、それぞれ対応する反転時間ＴＩ２００ｍｓ、４００ｍｓ、６００ｍｓの付箋とともに表示される。ラインマーク１６とそれぞれ対応する反転時間ＴＩとの関係は、任意に設定可能であり、操作者が数値設定するようにしてもよいし、検査部位に応じてデフォルト値として設定されるようにしてもよい。
【００３９】
各ラインマーク１６は、それぞれ対応する反転時間ＴＩに、想定した流速ｖをかけた距離に相当する距離を、ラベル領域マーク１５の端（又は中央）から、画像の縮尺で離間させた位置であって、ラベル領域マーク１５の長軸と平行に配置される。
【００４０】
操作者により、入力器１１３の例えばマウス等のポインタを操作して、参照画像上でラベル領域のマーク１５が移動されたとき、ＭＲＡエキスパートシステム１２０は、ラベル領域のマーク１５の移動に従って、参照画像上のラインマーク１６の位置を再計算し、その位置にラインマーク１６を再表示させる。それによりラベル領域のマーク１５の移動に追従して、ラインマーク１６の表示位置が移動する。また、操作者により想定流速ｖが変更されたとき、ＭＲＡエキスパートシステム１２０は、その変更された想定流速ｖに従って、参照画像上のラインマーク１８の位置を再計算し、その位置にラインマーク１８を再表示させる（図５（ｂ）参照）。さらに、第１の支援モードと同様に、ラベル領域のマーク１５を回転操作することが可能であり、そのときラベル領域のマーク１５の長軸との平行を維持するようにラインマーク１６は傾斜される。また、ラベル領域のマーク１５の形状が円形に変更されたとき、そのラベル領域のマーク１５と同心円を描くように、ラベル血流の到達位置を表すマーク１６が円形に変更される。
【００４１】
このように本実施形態では、反転時間ＴＩを想定して様々な流速ｖでのそれぞれの到達位置を、ラベル領域の位置設定又は反転時間ＴＩの調整のための指標として表示し、また逆に、流速ｖを想定して様々な反転時間ＴＩでのそれぞれの到達位置を、ラベル領域の位置設定の指標として表示するようにしたことで、好適な位置にラベル領域を設定し、又は好適な反転時間ＴＩを設定することができる。
【００４２】
以上のように設定されたラベル領域の位置、反転時間ＴＩに従って、ＭＲＡエキスパートシステム１２０は、標準パルスシーケンスを更新し、そのデータをシーケンサ１０５に供給する。シーケンサ１０５は、供給されたパルスシーケンスのデータに従って、傾斜磁場電源１０４、送信器１０８Ｔ及び受信器１０８Ｒを制御することにより、当該パルスシーケンスを実行する。
【００４３】
次に、パルスシーケンスの実行により収集されたＭＲデータにより演算ユニット１１０で画像データが生成され、表示器１１２で表示される。ＭＲＡエキスパートシステム１２０は、この画像データの表示に際して、診断支援情報を提供する機能を備える。
【００４４】
図６は、画像の表示画面の一例を示している。２４は生成された画像であり、この画像２４上に、ラベル領域を示すマーク２１と、血流速度スケール２２とが重ねられる。なお、２３は、ラベルされた血流像を示している。速度スケール２２には、撮影のために使われたパルスシーケンスの反転時間ＴＩに対して、様々な血流速度ｖとかけた距離に相当する画像上の距離で目盛られている。このスケール２２により、ラベルされた血流２３の流動経路を推定することができる。
【００４５】
さらに、ＭＲＡエキスパートシステム１２０は、血流の平均速度を計算する機能を備えている。図７は、血流の平均速度計算時の画面を示している。２５はラベル領域マークである。画像上には、血管像３１とともに、ラベルされた血流３０が例えば低信号として表示される。操作者は、ポインタ２９を操作して、ラベル領域から、測定対象としてのラベルされた血流３０までの経路を、血管像３１に従って複数の直線２７を関節（頂点）２６で近似的になぞっていく。つまり、操作者は、ラベル領域の下の端から、撮影対象３０が移動したと思われる経路をポインタ２９でなぞりながら順次マウスクリックにより関節２６を指定していく。
【００４６】
ＭＲＡエキスパートシステム１２０は、折れ線２７の距離Ｌ[m] を実際の距離に換算すると共に、その距離Ｌを、反転時間ＴＩ[m/s] でわり算することにより、平均速度を計算し、それを数値２８で表示する。これにより、撮影対象の速度が正確にしかも簡便に測定できるようになり、操作者にとって便利である。必要に応じて、エコー時間ＴＥ[s] の間の撮影対象の移動による速度の補正をしてもよい。
【００４７】
図８は、ＭＲＡエキスパートシステム１２０により提供されるパルスシーケンスを示している。Ａで示す部分は、図９で示した従来例と同じである。Ｂで示す部分は、Ａ部分から、インバージョンパルス１，スライス選択用傾斜磁場パルス２，ラベル領域のロケーションのための周波数シフトパルス３を含むプリサチュレーション部分を削除して、信号発生部分４０だけの部分である。機能的には、インバージョンパルス１だけを削除してもよい。Ａ部分とＢ部分の撮影を交互に行い、それぞれ別に画像を再構成し、２組の画像データを得る。得られた画像データの差画像を作成することにより、撮影対象（ここでは血流）の動きのない部分からの信号を抑制できる。これにより、動きのある部分をより明瞭に描出することができる。図６の速度スケールは２２は図８の差画像にも表示される。また図８の差画像上で、図７の平均速度計算のための設定を行うようにしてもよい。
【００４８】
図８のパルスシーケンスは、プレサチュレーション部分（インバージョンパルス１，スライス選択用傾斜磁場パルス２，ラベル領域のロケーションのための周波数シフトパルス３）の後であって、信号発生部分（実際のイメージングパルス列部分）４０の前に、高周波反転パルス１を、スライス選択用傾斜磁場パルス２及びラベル領域のロケーションのための周波数シフトパルス３を印加しないで、もう一度印加することによって、プレサチュレート領域を、Ａ部分の反転励起領域以外の領域に設定するようにしてもよい。つまり、上記方法では、ラベル領域を、それ以外の領域に対して、低信号化することでラベル化を実現していたが、逆に、ラベル領域以外の領域を低信号化することでラベル領域を抽出できるようにしてもよい。
【００４９】
また、図８で４０で示される実際のイメージングの部分は、ＦＥ法の場合を示したが、その他のパルスシーケンス例えば、スピンエコー(spin echo) 法、エコープラナー(echo planar) 法、レア(RARE)法（fast spin echo法ともいう）、trueFISP法などでも可能である。また、撮影対象の動きは、心時相、呼吸による体動に関係したものも多いため、高周波反転パルスの印加に関しては、心電同期、脈波同期、呼吸同期などの波形に合わせて同期させて収集してもよい。呼吸同期波形のトリガ後、始めての心電波形に同期させて収集してもよい。
【００５０】
（変形例）
本発明は、上述した実施形態に限定されるものではなく、実施段階ではその要旨を逸脱しない範囲で種々変形して実施することが可能である。さらに、上記実施形態には種々の段階が含まれており、開示される複数の構成要件における適宜な組み合わせにより種々の発明が抽出され得る。例えば、実施形態に示される全構成要件から幾つかの構成要件が削除されてもよい。
【図面の簡単な説明】
【００５１】
【図１】本発明による磁気共鳴映像化の実施形態の構成を示すブロック図。
【図２】図１のＭＲＡエキスパートシステムによる流速指標を含むラベル領域設定画面の一例を示す図。
【図３】図１のＭＲＡエキスパートシステムによる流速指標を含むラベル領域設定画面の他の例を示す図。
【図４】図１のＭＲＡエキスパートシステムによる流速指標を含むラベル領域設定画面のさらに他の例を示す図。
【図５】図１のＭＲＡエキスパートシステムによるＴＩ指標を含むラベル領域設定画面の一例を示す図。
【図６】図１のＭＲＡエキスパートシステムにより再構成画像とともに表示されるラベル領域マーク及び流速スケールの一例を示す図。
【図７】図１のＭＲＡエキスパートシステムによる平均速度計算のための設定画面の一例を示す図。
【図８】図１のＭＲＡエキスパートシステムによるＭＲＡサブトラクション用のパルスシーケンスの一例を示す図。
【図９】一般的なフローイメージング法のパルスシーケンスの一例を示す図。
【図１０】従来のラベル領域設定画面を示す図。
【符号の説明】
【００５２】
１０１…磁石
１０２…静磁場電源
１０３…傾斜磁場コイルユニット
１０４…傾斜磁場電源
１０５…シーケンサ
１０６…ホスト計算機
１０７…ＲＦコイル
１０８Ｔ…送信器
１０８Ｒ…受信器
１１０…演算ユニット
１１１…記憶ユニット
１１２…表示器
１１３…入力器
１２０…ＭＲＡエキスパートシステム。
【図１】
0004519827.tif 000002
【図２】
0004519827.tif 000003
【図３】
0004519827.tif 000004
【図４】
0004519827.tif 000005
【図５】
0004519827.tif 000006
【図６】
0004519827.tif 000007
【図７】
0004519827.tif 000008
【図８】
0004519827.tif 000009
【図９】
0004519827.tif 000010
【図１０】
0004519827.tif 000011