(57)【要約】【課題】被検者、撮像対象部位、PLD値に依存せず、適切なBGSパルス印加タイミングで撮像することができるASL Perfusion撮像可能な磁気共鳴イメージング装置を提供する。【解決手段】血液プロトンをラベルしたラベル画像とラベルしないコントロール画像との差分により血行動態を描出する灌流画像の撮像において、プリスキャンを行って撮像対象部位のT1マップを取得し、背景信号として抑制する対象として設定された組織からのエコー信号を抑制するBGSパルスの印加タイミングを、取得されたT1値を用いて算出し、パルスシーケンスに反映させる。【選択図】図３
2019076182.tif 000002
【特許請求の範囲】
【請求項１】
被検体が配置される空間に静磁場を発生させる静磁場発生源と、
前記静磁場に重畳する傾斜磁場を発生させる傾斜磁場コイルと、
前記被検体に高周波磁場パルスを照射する送信RFコイルと、
前記送信RFコイルによる高周波磁場パルスの印加に応じて前記被検体より放出されるエコー信号を検出する受信RFコイルと、
所定のパルスシーケンスに基づいて、前記傾斜磁場コイルによる傾斜磁場の発生、前記送信RFコイルによる高周波磁場パルスの照射及び前記受信RFコイルによるエコー信号の検出を制御する計測制御部と、
前記計測制御部を制御する全体制御部とを有し、
前記全体制御部は、血液プロトンをラベルしたラベル画像とラベルしないコントロール画像との差分により血行動態を描出する灌流画像の撮像において、プリスキャンを行って撮像対象部位のT1マップを取得し、背景信号として抑制する対象として設定された組織からのエコー信号を抑制するBGSパルスの印加タイミングを、取得されたT1値を用いて算出し、前記計測制御部のパルスシーケンスに反映させる磁気共鳴イメージング装置。
【請求項２】
請求項１において、
前記BGSパルスは、プリサチレーションパルスと１以上のBGS IRパルスを含み、
前記BGS IRパルスの数は、前記背景信号として抑制する対象として設定された組織の数に応じて設定される磁気共鳴イメージング装置。
【請求項３】
請求項２において、
制御情報を入力するための制御情報入力画面が表示される表示部と、
前記制御情報を入力する操作部とを有し、
前記制御情報入力画面には、取得した前記撮像対象部位のT1マップが表示され、表示された前記撮像対象部位のT1マップに対して前記操作部によりROIが設定され、
前記全体制御部は、設定された前記ROIの数およびそのT1値を用いて前記BGSパルスの印加タイミングを算出する磁気共鳴イメージング装置。
【請求項４】
請求項２において、
制御情報を入力するための制御情報入力画面が表示される表示部と、
前記制御情報を入力する操作部とを有し、
前記制御情報入力画面には、取得した前記撮像対象部位のT1マップのT1値分布が表示され、表示された前記T1値分布に対して前記操作部により１以上のT1値が指定され、
前記全体制御部は、指定された前記T1値の数およびその値を用いて前記BGSパルスの印加タイミングを算出する磁気共鳴イメージング装置。
【請求項５】
請求項２において、
制御情報を入力するための制御情報入力画面が表示される表示部と、
前記制御情報を入力する操作部とを有し、
前記制御情報入力画面には、取得した前記撮像対象部位のT1マップのT1値分布が表示され、表示された前記T1値分布に対して前記操作部により１以上のT1値帯が指定され、
前記全体制御部は、指定された前記T1値帯の数およびその代表値を用いて前記BGSパルスの印加タイミングを算出する磁気共鳴イメージング装置。
【請求項６】
請求項１〜５のいずれか１項において、
前記撮像対象部位は、頭部である磁気共鳴イメージング装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は磁気共鳴イメージング装置に関する。
【背景技術】
【０００２】
磁気共鳴イメージング装置（MRI（Magnetic Resonance Imaging）装置）は、核磁気共鳴（NMR：Nuclear Magnetic Resonance）現象を利用して被検体の画像を得る計測装置である。被検体に高周波磁場パルス（以下、RFパルスという）を照射し、その応答として、被検体の組織を構成する原子核スピンが発生するNMR信号を計測する。そして計測したNMR信号に基づいて、被検体の頭部、腹部、四肢等の形態や機能を２次元的にあるいは３次元的に画像化する。撮像の際には、傾斜磁場によって、NMR信号は異なる位相エンコードやスライスエンコードが付与されると共に、周波数エンコードが付与されて、時系列データとして計測される。計測されたNMR信号は、２次元または３次元フーリエ変換されることにより画像に再構成される。
【０００３】
MRI装置において、血液プロトンをラベルした画像（ラベル画像）とラベルしない画像（コントロール画像）の差分により血行動態を描出する撮像が行われている。その方法の代表例として、特許文献１や特許文献２のような方法がある。これらは非造影で灌流画像を得るものである。この灌流画像は、一般的にArterial Spin Labeling Perfusion（ASL Perfusion）画像と呼ばれている。なお、ラベル／コントロールした血液プロトンが灌流領域に到達する以前のタイミングで撮像した場合は血管画像を得ることができる。ラベル／コントロールしてからリードアウトシーケンス開始までの時間を一般的にPost Labeling Delay（PLD）という。
【０００４】
ASL Perfusion撮像は１％程度の極めて小さい信号差を画像化する撮像法であるため、ラベル画像およびコントロール画像（元画像）の背景信号を抑制した状態でデータを取得する。一般的に背景信号は、特許文献３のようにプリサチレーションパルス（BGS Presat）やIRパルス（BGS IRパルス）の組み合わせ（BGSパルス）で抑制する。BGSパルスの印加タイミングは、頭部であれば、例えば非特許文献１や非特許文献２の白質や灰白質、脳脊髄液（CSF）などのT1値（縦緩和時間）、T2値（横緩和時間）、M0値（初期縦磁化）を使用して求めている。
【０００５】
一方、非特許文献３や非特許文献４のような手法でT1値を計測し、T1マップを取得する方法が知られている。これらの手法は一般的に心臓領域を対象とすることが多いが、ASL Perfusion撮像においてT1値の分布を知ることはASL Perfusion画像の信頼性および定量性を向上するために重要である。
【先行技術文献】
【特許文献】
【０００６】
【特許文献１】米国特許第５，８４６，１９７号
【特許文献２】米国特許第７，５４５，１４２号
【特許文献３】米国特許第５，５８８，４３１号
【非特許文献】
【０００７】
【非特許文献１】Magnetic Resonance in Medicine 54:507-512 (2005)
【非特許文献２】Magnetic Resonance in Medicine 60:1155-1168 (2008)
【非特許文献３】Magnetic Resonance in Medicine 52:141-146 (2004)
【非特許文献４】Magnetic Resonance in Medicine 71:2082-2095 (2014)
【発明の概要】
【発明が解決しようとする課題】
【０００８】
各組織のT1値は被検者によっては文献値と数10ms〜数100msの誤差があり、従来技術ではその誤差により元画像が背景信号抑制不足となる可能性がある。元画像が背景信号抑制不足となった結果、ASL Perfusion画像の安定性およびその画像から求める定量値マップの信頼性が低下する可能性がある。また、撮像対象部位毎にBGSパルスの印加タイミングを最適化する必要がある（対象部位によりT1値が異なる）ため、例えば頭部のBGSパルス印加タイミングを他の撮像部位に展開することはできない。
【課題を解決するための手段】
【０００９】
本発明の一実施の形態である磁気共鳴イメージング装置においては、被検体が配置される空間に静磁場を発生させる静磁場発生源と、静磁場に重畳する傾斜磁場を発生させる傾斜磁場コイルと、被検体に高周波磁場パルスを照射する送信RFコイルと、送信RFコイルによる高周波磁場パルスの印加に応じて被検体より放出されるエコー信号を検出するRFコイルと、所定のパルスシーケンスに基づいて、傾斜磁場コイルによる傾斜磁場の発生、送信RFコイルによる高周波磁場パルスの照射及び受信RFコイルによるエコー信号の検出を制御する計測制御部と、計測制御部を制御する全体制御部とを有し、全体制御部は、血液プロトンをラベルしたラベル画像とラベルしないコントロール画像との差分により血行動態を描出する灌流画像の撮像において、プリスキャンを行って撮像対象部位のT1マップを取得し、背景信号として抑制する対象として設定された組織からのエコー信号を抑制するBGSパルスの印加タイミングを、取得されたT1値を用いて算出し、計測制御部のパルスシーケンスに反映させる。
【発明の効果】
【００１０】
被検者、撮像対象部位、PLD値に依存せず、ASL Perfusion撮像シーケンスを適切化することができ、背景信号抑制不足がなくなり、ASL Perfusion画像の安定性およびその画像から求める定量値マップの信頼性が向上する。
【図面の簡単な説明】
【００１１】
【図１】MRI装置の全体基本構成を示す図である。
【図２】一般的なASL Perfusion撮像シーケンスのシーケンス図である。
【図３】実施例１におけるASL Perfusion検査のフローチャートである。
【図４】全体制御部内の機能ブロックを示す図である。
【図５】ASL Perfusion撮像シーケンスにおける各時間間隔を示す図である。
【図６】ROI設定組織選択画面の一例である。
【図７】実施例２におけるユーザー選択画面の一例である。
【図８】実施例２におけるフローチャートである。
【図９】実施例３におけるユーザー選択画面の一例である。
【図１０】実施例３におけるフローチャートである。
【発明を実施するための形態】
【００１２】
以下、添付図面に従って本発明のMRI装置の好ましい実施形態について詳説する。なお、発明の実施形態を説明するための全図において、同一機能を有するものは同一符号を付け、その繰り返しの説明は省略する。
【００１３】
図１は、MRI装置の一実施例の全体構成を示すブロック図である。MRI装置は、NMR現象を利用して被検体101の断層画像を得るもので、静磁場発生磁石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の３軸方向に巻かれたコイルであり、それぞれの傾斜磁場コイルは、それを駆動する傾斜磁場電源109に接続され、電流が供給される。具体的には、各傾斜磁場コイルの傾斜磁場電源109は、それぞれ後述の計測制御部111からの命令に従って駆動されて、それぞれの傾斜磁場コイルに電流を供給する。これにより、X, Y, Zの３軸方向に傾斜磁場Gx，Gy，Gzが発生する。
【００１６】
２次元スライス面の撮像時には、スライス面（撮像断面）に直交する方向にスライス傾斜磁場パルス（Gs）が印加されて被検体101に対するスライス面が設定され、そのスライス面に直交して且つ互いに直交する残りの２つの方向に位相エンコード傾斜磁場パルス（Gp）と周波数エンコード（リードアウト）傾斜磁場パルス（Gf）が印加されて、エコー信号にそれぞれの方向の位置情報がエンコードされる。
【００１７】
送信RFコイル104は、被検体101にRFパルスを照射するコイルであり、RF送信部110に接続され高周波パルス電流が供給される。これにより、被検体101の生体組織を構成する原子の原子核スピンにNMR現象が誘起される。具体的には、RF送信部110が、後述の計測制御部111からの命令に従って駆動されて、高周波パルスが振幅変調され、増幅された後に被検体101に近接して配置された送信RFコイル104に供給されることにより、RFパルスが被検体101に照射される。
【００１８】
受信RFコイル105は、被検体101の生体組織を構成する原子核スピンのNMR現象により放出されるNMR信号（エコー信号）を受信するコイルであり、信号検出部106に接続され、受信RFコイル105で受信されたエコー信号は信号検出部106に送られる。信号検出部106は、受信RFコイル105で受信されたエコー信号の検出処理を行う。具体的には、送信RFコイル104から照射されたRFパルスによって誘起された被検体101の応答のエコー信号が被検体101に近接して配置された受信RFコイル105で受信され、後述の計測制御部111からの命令に従って、信号検出部106が、受信されたエコー信号を増幅し、直交位相検波により直交する２系統の信号に分割し、それぞれを所定数（例えば128、256、512等）サンプリングし、各サンプリング信号をＡ／Ｄ変換してディジタル量に変換し、後述の信号処理部107に送る。 従って、エコー信号は所定数のサンプリングデータからなる時系列のデジタルデータ（以下、エコーデータという）として得られる。信号処理部107は、エコーデータに対して各種処理を行い、処理されたエコーデータを計測制御部111に送る。
【００１９】
計測制御部111は、被検体101の断層画像の再構成に必要なデータ収集のための種々の命令を、主に、傾斜磁場電源109と、RF送信部110と、信号検出部106に送信してこれらを制御する制御部である。具体的には、計測制御部111は、後述する全体制御部108の制御で動作し、ある所定のパルスシーケンスに基づいて、傾斜磁場電源109、RF送信部110及び信号検出部106を制御して、被検体101へのRFパルスと傾斜磁場パルスの印加及び被検体101からのエコー信号の検出を繰り返し実行し、被検体101の撮像領域についての画像の再構成に必要なエコーデータを収集する。
【００２０】
全体制御部108は、計測制御部111の制御、及び、各種データ処理と処理結果の表示及び保存等の制御を行うものであって、CPU及びメモリを内部に有する演算処理部114と、光ディスク、磁気ディスク等の記憶部115とを有する。具体的には、計測制御部111を制御してエコーデータの収集を実行させ、計測制御部111からのエコーデータが入力されると、演算処理部114がそのエコーデータに印加されたエンコード情報に基づいて、メモリのｋ空間に相当する領域に記憶させる。メモリのｋ空間に相当する領域に記憶されたエコーデータ群をｋ空間データともいう。そして演算処理部114はこのｋ空間データに対して信号処理やフーリエ変換による画像再構成等の処理を実行し、その結果である被検体101の画像を、後述の表示・操作部113に表示させると共に記憶部115に記録する。
【００２１】
表示・操作部113は、再構成された被検体101の画像等を表示する表示部と、MRI装置の各種制御情報や全体制御部108で行う処理の制御情報を入力するトラックボール又はマウス及びキーボード等の操作部とを有する。この操作部は表示部に近接して配置され、表示部には制御情報を入力するための制御情報入力画面が表示され、操作者は表示部を見ながら操作部を通してインタラクティブにMRI装置の各種処理を制御する。
【００２２】
また、本実施の形態に係るMRI装置は、被検体101の体動・心電情報を検出する体動・心電情報検出機構を備える。この体動・心電情報検出機構は、被検体101に装着されて被検体の体動・心電情報を検出するセンサー部116と、センサー部116からの信号を処理して、その処理した体動・心電情報を計測制御部111に送る体動・心電情報処理部117とを有している。体動・心電情報検出機構が被検体101の呼吸波形を検出するものであれば、センサー部116は呼吸波形を検出するセンサーであり、体動・心電情報処理部117は呼吸センサーからのアナログ信号を処理する。計測制御部111は、体動・心電情報処理部117で検出された被検体101の体動・心電情報に同期させて、パルスシーケンスの実行による撮像を行なう同期撮像を制御する。
【００２３】
なお、図１において、送信側の送信RFコイル104と傾斜磁場コイル103は、被検体101が挿入される静磁場発生磁石102の静磁場空間内に、垂直磁場方式であれば被検体101に対向して、水平磁場方式であれば被検体101を取り囲むようにして設置されている。また、受信側の受信RFコイル105は、被検体101に対向して、あるいは取り囲むように設置されている。
【００２４】
現在のMRI装置の撮像対象核種は、臨床で普及しているものとしては、被検体の主たる構成物質である水素原子核（プロトン）である。プロトン密度の空間分布や、励起状態の緩和時間の空間分布に関する情報を画像化することで、人体頭部、腹部、四肢等の形態または、機能を２次元もしくは３次元的に撮像する。
【００２５】
ASL Perfusion撮像シーケンスは図２に示すようにLabel RFパルス201を照射してデータを取得するラベルスキャン、Control RFパルス202を照射してデータを取得するコントロールスキャンが入れ子になった形である。Label RFパルス201は血液プロトンをラベル（静磁場方向に対して反転）するためのパルス、Control RFパルス202は実質的には血液プロトンをラベルせず、ＭＴ（磁化移動：Magnetization Transfer）効果をLabel RFパルス201と揃えるためのパルスである。
【００２６】
ASL Perfusion撮像は１％程度の極めて小さい信号差を画像化する撮像法であるため、ラベル画像およびコントロール画像（元画像）の背景信号を抑制した状態でデータを取得する。背景信号はプリサチレーションパルス（BGS （Background Suppressing）Presat）やIRパルス（BGS IRパルス）の組み合わせ（BGSパルス）で抑制する。図２のシーケンスでは、BGS Presat 203の１発、BGS IRパルス204〜206の３発でBGSパルスを構成している。BGSパルスの印加タイミングは、従来、例えば頭部であれば白質や灰白質、CSFなどのT1値の文献値を使用して求めていた。
【００２７】
データ取得207（リードアウトシーケンスの開始）は血液プロトンをラベルまたはコントロールしてから待ち時間PLDを待ってから行うが、これはラベルまたはコントロールされた血液プロトンが撮像スライス（スラブ）に流入するための待ち時間である。データ取得は２Ｄまたは３Ｄで行う。ラベルスキャンから得られたラベル画像とコントロールスキャンから得られたコントロール画像の差分により、ASL Perfusion画像を得る。
【実施例１】
【００２８】
実施例１におけるASL Perfusion検査のフローチャートを図３に示す。ASL Perfusion検査を開始（301）し、プリスキャンとしてT1マップを取得（302）する。ステップ302で取得したT1マップに対し、抑制したい組織上（撮像対象部位が頭部であれば白質／灰白質／CSFなど）にユーザーがROI（Region of Interest）を設定し、設定した組織を選択（303）する。ステップ303でユーザーが設定したROIの数およびそのT1値に対し、BGS IRパルス印加数およびBGSパルス印加タイミングを最適化（304）する。BGS IRパルス印加数はROI設定数と同一とする。抑制するROIの数以上のBGS IRパルスを印加する必要のあるところ、同数とするものである。BGS IRパルス印加数および最適化したBGSパルス印加タイミングをシーケンスに反映（305）し、ASL Perfusion撮像を開始（306）する。撮像終了後、ASL Perfusion画像を表示（307）し、ASL Perfusion検査が終了となる。
【００２９】
図４に全体制御部108内部の機能ブロックを示す。表示・操作部113でユーザーがT1マップ上に設定したROIに対し、設定したROIの数をBGS IRパルス印加数判定部1081で判定する。判定したBGS IRパルス印加数およびROIのT1値に対し、BGSパルス印加タイミング算出部1082でBGSパルスの印加タイミングを計算する。計算したBGSパルス印加タイミングをシーケンス調整部1083でシーケンスに反映し、計測制御部111に送る。
【００３０】
対象部位が頭部、ユーザーが白質、灰白質、CSFにROIを設定した場合のシーケンス調整の例について、図５を用いて詳細に説明する。まず、ユーザーのROI設定数が３個なのでBGS IRパルスは３発となる。次に、設定したROIのT1値に対するBGSパルス印加タイミングの最適化方法について説明する。Bloch方程式を用いて複数IRパルスに対する縦磁化は一般的に（数１）および（数２）で表現できる。
【００３１】
【数１】
2019076182.tif 000003
【００３２】
【数２】
2019076182.tif 000004
【００３３】
ここで、全IRパルスのFA（Flip Angle）を180deg.とし、BGS Presat 203で全ての組織の縦磁化がゼロになると仮定すると、（数１）および（数２）はそれぞれ（数３）および（数４）に書き換えることができる。
【００３４】
【数３】
2019076182.tif 000005
【００３５】
【数４】
2019076182.tif 000006
【００３６】
（数３）および（数４）より、３発目のIRパルス印加後の縦磁化Ｍｚ(3)は（数５）で表現できる。
【００３７】
【数５】
2019076182.tif 000007
【００３８】
Imagingの励起RFパルス207の照射直前の縦磁化ＭｚImagingは、Ｍｚ(3)においてBGS IRパルス206の照射より待ち時間B4が経過した状態であるから、（数５）にＴｄ(3)＝B4を代入して（数６）となる。
【００３９】
【数６】
2019076182.tif 000008
【００４０】
ｋ空間中心に配置されるデータ信号の最大値および最小値はＭｚImagingを用いて（数７）で表現できる。なお、TEはエコー時間である。
【００４１】
【数７】
2019076182.tif 000009
【００４２】
（数７）を用いて、評価関数ｆ１は（数８）として定義される。
【００４３】
【数８】
2019076182.tif 000010
【００４４】
評価関数ｆ１は制約付非線形多変数関数であり、所定の制約条件のもとで評価関数ｆ１を最小化するB1、B2、B3およびB4を探索することでBGSパルス印加タイミングを最適化する。BGSパルスは背景となる組織（この場合は白質、灰白質、CSF）からの信号を抑制することを目的としている。（数８）に示した評価関数ｆ１を最小化することは、これらの組織からの信号値をできるだけコントラストのつかない、フラットな状態にすることになるため、BGSパルス印加タイミングが最適化されるものである。
【００４５】
この場合の制約条件は（数９）となる。なお、Label DurationはLabel RFパルス201のパルス幅でありシステムに予め設定された値、Imaging DurationはImagingの励起RFパルス207のパルス幅であり、ユーザーが入力したパラメータにより自動計算される値、PLDはユーザーが入力するパラメータ値である。
【００４６】
【数９】
2019076182.tif 000011
【００４７】
また、Ｓmaxの上限やＳminの下限も予めシステムに設定しておく。Ｓmaxの上限を0.1、Ｓminの下限を0.02とした場合、ROIを設定した組織の信号値を元の信号値の２〜１０％の範囲に均一に抑制することができる。これは制約条件としては、（数１０）のように書ける。
【００４８】
【数１０】
2019076182.tif 000012
【００４９】
以上より、（数９）、（数１０）の制約条件を満たしつつ、評価関数ｆ１を最小化するB1、B2、B3およびB4を探索して求めることにより、白質、灰白質、CSFに対してBGSパルス印加タイミングを最適化する。
【００５０】
また、以上の数式で用いるT2値、M0値は非特許文献１、非特許文献２に記載の値を用いればよい。本実施例においては、BGS IRパルスによる信号の抑制の程度は、スピンが平衡状態に戻るためにRFパルスから獲得したエネルギーを周囲の格子に与えるT1緩和が支配的であるため、T1値については実測するが、T2値、M0値については文献値を用いる。ユーザーは選択画面（図６にその一例を示す）において、ROIを設定した組織に対応する組織のチェックボックスにチェックを入れる（図３のステップ303を参照）。非特許文献１、非特許文献２に記載されているT2値、M0値を予めデータベースに保持しておき、チェックが入れられた組織に対するT2値、M0値を読み出して、BGSパルス印加タイミングを最適化する。ユーザーがチェックボックスにチェックを入れる手順を省略した場合は、M0＝１およびT2＝∞（すなわち、exp(-TE/T2)＝１）としてBGSパルス印加タイミングを最適化してもよい。
【実施例２】
【００５１】
実施例２では、T1マップのT1値分布をバーで表示し、バーをユーザーがクリックして抑制対象とするT1値を決定する。実施例１ではT1マップ上で抑制したい組織を指定するのに対して、実施例２は、抑制したい組織によってそのT1値がおよそ定まっていることを利用して、取得したT1値に基づき、抑制するT1値を直接指定するものである。ユーザーがバーをクリックした回数が抑制するT1値の個数となる。図７に実施例２におけるユーザー選択画面の一例を示す。組織におけるT1値は400ms〜2000msの分布をもつことがバー700によって示され、これに対してユーザーがT1値を３つ選択した状態を示す。クリックしたT1値およびその個数に対応した値が表示（701〜703）され、表示されたT1値に対してBGSパルス印加タイミングを最適化する。
【００５２】
実施例２におけるフローチャートを図８に示す。ASL Perfusion検査を開始（801）し、プリスキャンとしてT1マップを取得（802）する。ステップ802で取得したT1マップのT1値分布を表したバーに対し、ユーザーが抑制したいT1値およびその個数を選択し、抑制したい組織も選択（803）する。ステップ803でユーザーが選択したT1値および選択数に対し、BGS IRパルス印加数およびBGSパルス印加タイミングを最適化（804）する。BGS IRパルス印加数は実施例１と同様に選択数と同一とする。BGS IRパルス印加数および最適化したBGSパルス印加タイミングをシーケンスに反映（805）し、ASL Perfusion撮像を開始（806）する。撮像終了後、ASL Perfusion画像を表示（807）し、ASL Perfusion検査が終了となる。これにより、様々なT1値を持つ組織を均一に抑制した元画像から灌流画像を取得することができる。
【実施例３】
【００５３】
実施例３では、T1マップのT1値分布を最小値から最大値まで例えば50msステップで分割し、ユーザーが抑制したいT1値帯を選択する。実施例２が抑制するT1値をポイントで指定するのに対し、実施例３では抑制するT1値を範囲で指定できる。分割するためのステップ数は予めシステムに登録しておいてもよいし、ユーザーが入力できるようにしてもよい。図９に実施例３におけるユーザー選択画面の一例を示す。バー700に示されるように400ms〜2000msのT1値分布を持つ組織に対し、ユーザーが抑制対象のT1値帯として、1900〜1950ms（902）、1000〜1050ms（903）、450〜500ms（904）を選択した例を示す。選択したT1値帯および選択数に対してBGSパルス印加タイミングを最適化する。BGSパルス印加タイミングの算出にはT1値帯の代表値としてその中間値（T1値帯902であれば1925ms）を用いる。
【００５４】
実施例３におけるフローチャートを図１０に示す。ASL Perfusion検査を開始（1001）し、プリスキャンとしてT1マップを取得（1002）する。ステップ1002で取得したT1マップのT1値分布に対し、ユーザーが抑制したいT1値帯を選択肢の中から選択し、抑制したい組織も選択（1003）する。ステップ1003でユーザーが選択したT1値帯および選択数に対し、BGS IRパルス印加数およびBGSパルス印加タイミングを最適化（1004）する。BGS IRパルス印加数は実施例１と同様に選択数と同一とする。BGS IRパルス印加数および最適化したBGSパルス印加タイミングをシーケンスに反映（1005）し、ASL Perfusion撮像を開始（1006）する。撮像終了後、ASL Perfusion画像を表示（1007）し、ASL Perfusion検査が終了となる。これにより、様々なT1値を持つ組織を均一に抑制した元画像から灌流画像を取得することができる。
【符号の説明】
【００５５】
101…被検体、102…静磁場発生磁石、103…傾斜磁場コイル、104…送信RFコイル、105…受信RFコイル、106…信号検出部106、107…信号処理部、108…全体制御部、109…傾斜磁場電源、110…RF送信部、111…計測制御部、112…ベッド、113…表示・操作部、114…演算処理部、116…センサー部、117…体動・心電情報処理部。
【図１】
2019076182.tif 000013
【図２】
2019076182.tif 000014
【図３】
2019076182.tif 000015
【図４】
2019076182.tif 000016
【図５】
2019076182.tif 000017
【図６】
2019076182.tif 000018
【図７】
2019076182.tif 000019
【図８】
2019076182.tif 000020
【図９】
2019076182.tif 000021
【図１０】
2019076182.tif 000022