(57)【要約】【課題】診断効率のよい造影ダイナミックスキャン法による撮像が可能な磁気共鳴イメージング装置を提供する。【解決手段】造影ダイナミックスキャンにおける各スキャンが実施された複数の時点の各々を処理対象として、対象の時点における被検体の画像に対応する再構成前のローデータと、被検体の基準画像に対応するローデータとの相違を表す指標値を算出する（ＳＴ３０３）。その指標値がしきい値以上であるかを判定し（ＳＴ３０４）、肯定されるときにその画像を再構成する（ＳＴ３０５）。基準画像は、例えば、非造影による被検体の画像とする。これにより、再構成される不要な画像を低減し、検査時間や診断時間を短縮する。【選択図】図６
2010104459.tif 000002
【特許請求の範囲】
【請求項１】
造影剤が注入された被検体を複数の時点にてＭＲスキャンすることにより前記複数の時点における前記被検体の画像のローデータを取得するスキャン部と、
前記複数の時点の各々を対象として、対象の時点における前記被検体の画像のローデータに基づいて、該画像を再構成すべきであるか否かを判定する判定部と、
前記判定部により再構成すべきであると判定された画像のローデータに基づいて該画像を再構成する再構成部とを備える磁気共鳴イメージング装置。
【請求項２】
前記判定部は、前記対象の時点における前記被検体の画像のローデータと、前記被検体の基準画像のローデータとの相違に基づいて判定する請求項１に記載の磁気共鳴イメージング装置。
【請求項３】
前記判定部は、両ローデータに設定された関心領域における相違に基づいて判定する請求項２に記載の磁気共鳴イメージング装置。
【請求項４】
造影剤が注入された被検体を複数の時点にてＭＲスキャンすることにより前記複数の時点における前記被検体の画像のローデータを取得するスキャン部と、
前記取得されたローデータに基づいて前記複数の時点における前記被検体の画像を再構成する再構成部と、
前記複数の時点の各々を対象として、対象の時点における前記被検体の画像に基づいて該画像を出力すべきであるか否かを判定する判定部とを備える磁気共鳴イメージング装置。
【請求項５】
前記判定部は、前記対象の時点における前記被検体の画像と、前記被検体の基準画像との相違に基づいて判定する請求項４に記載の磁気共鳴イメージング装置。
【請求項６】
前記判定部は、両画像に設定された関心領域における相違に基づいて判定する請求項５に記載の磁気共鳴イメージング装置。
【請求項７】
前記判定部は、
前記相違を表す指標値を算出する指標値算出部と、
前記指標値が所定のしきい値以上であるか否かに基づいて前記判定を行うしきい値判定部とを有する請求項２または請求項３に記載の磁気共鳴イメージング装置。
【請求項８】
前記指標値は、前記相違を形成する一方のローデータを構成するデータ値の総和と、該相違を形成する他方のローデータを構成するデータ値の総和との差分である請求項７に記載の磁気共鳴イメージング装置。
【請求項９】
前記判定部は、
前記相違を表す指標値を算出する指標値算出部と、
前記指標値が所定のしきい値以上であるか否かに基づいて前記判定を行うしきい値判定部とを有する請求項５または請求項６に記載の磁気共鳴イメージング装置。
【請求項１０】
前記指標値は、前記相違を形成する一方の画像データを構成するデータ値の総和と、該相違を形成する他方の画像データを構成するデータ値の総和との差分である請求項９に記載の磁気共鳴イメージング装置。
【請求項１１】
前記基準画像は、非造影による前記被検体の画像である請求項２、請求項３および請求項５から請求項１０のいずれか１項に記載の磁気共鳴イメージング装置。
【請求項１２】
前記基準画像は、前記対象の時点より前の該時点の近傍の時点における前記被検体の画像である請求項２、請求項３および請求項５から請求項１０のいずれか１項に記載の磁気共鳴イメージング装置。
【請求項１３】
前記被検体の画像は、前記被検体に設定された複数のスライスの全部または一部に対応する画像である請求項１から請求項１２のいずれか１項に記載の磁気共鳴イメージング装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、磁気共鳴イメージング（ＭＲＩ：Ｍａｇｎｅｔｉｃ  Ｒｅｓｏｎａｎｃｅ  Ｉｍａｇｉｎｇ）装置に関するものである。
【背景技術】
【０００２】
磁気共鳴イメージング装置は、静磁場空間内の被検体にＲＦパルス（ｐｕｌｓｅ）と呼ばれる電磁波を送信することにより、その被検体内のプロトン（ｐｒｏｔｏｎ）のスピン（ｓｐｉｎ）を核磁気共鳴現象によって励起させ、その励起されたスピンにより発生する磁気共鳴信号を得るスキャン（ｓｃａｎ）を実施する。そして、そのスキャンにより得られる磁気共鳴信号に基づいて、被検体の画像を再構成する装置である。
【０００３】
磁気共鳴イメージング装置を用いた撮像法の１つとして、造影剤が注入された被検体の撮像領域を複数の時相において撮像する造影ダイナミックスキャン（ｄｙｎａｍｉｃ  ｓｃａｎ）法（造影ダイナミック検査法、造影ダイナミック撮像法などとも呼ばれる）が知られている（特許文献１、段落０００７等参照）。具体的には、例えば、造影剤の被検体への注入を開始してから１０秒経過後、２０秒経過後、あるいは３０秒経過後といった各時相において撮像領域をスキャンする。これにより、例えば、造影剤が到達し始めたときの撮像領域の画像、造影剤が充満しているときの撮像領域の画像、あるいは、造影剤が抜け始めたときの撮像領域の画像などを、被検体の診断用画像として得ることができる。
【０００４】
ところが、造影剤が撮像領域に到達するタイミング（ｔｉｍｉｎｇ）や、造影剤が撮像領域から抜け始めるタイミングなどは、被検体や撮像領域によって異なるだけでなく、同一の被検体や撮像領域においても安定しないことがある。また、撮像領域に到達する造影剤の量や濃度も、被検体や撮像領域によって極端に少なくなることがあり、造影剤の到達をモニタ（ｍｏｎｉｔｏｒ）できないことがある。そのため、造影剤が所望の状態にある撮像領域をピンポイント（ｐｉｎ  ｐｏｉｎｔ）でスキャンすることが難しい。
【０００５】
そこで、実際には、次に例示するスキャン方法を用いることが多い。第１のスキャン方法は、撮像領域において造影剤が所望の状態になる時点を推定し、この推定した時点とその前後幾らかの時点についてのみスキャンを実施する方法である。そして、第２のスキャン方法は、スキャンする時点を特に決めず、造影剤を被検体に注入した後、撮像領域において造影剤が所望の状態になる時点を含むと考えられる十分長い期間について略定期的にスキャンを繰り返し実施する方法である。
【特許文献１】特開２００６−３２５９６７号公報
【発明の開示】
【発明が解決しようとする課題】
【０００６】
しかしながら、上記の第１および第２のスキャン方法では、造影剤が到達する前の時点における画像や、造影剤の状態に変化がない画像など、診断に不要な画像も含めて再構成されることになる。
【０００７】
そのため、画像の再構成処理に余分な時間が掛かり、検査時間が延びる原因となる。また、撮像技師は、再検査（再撮像）の必要性を判断するために、再構成された全画像の中に診断に必要な画像が含まれているか否かを確認したり、再構成された全画像の中から診断に必要な画像を分離して画質の良否を判断したりする必要があり、検査時間が延びる原因となる。さらに、読影医は、再構成された全画像の中から診断に必要な画像を検索して読影する必要があり、診断時間が延びる原因となる。このように検査時間や診断時間が延びると、その結果、診断効率が悪くなる。
【０００８】
本発明は、上記事情に鑑み、診断効率のよい造影ダイナミックスキャン法による撮像が可能な磁気共鳴イメージング装置を提供することを目的とする。
【課題を解決するための手段】
【０００９】
第１の観点では、本発明は、造影剤が注入された被検体を複数の時点にてＭＲスキャンすることにより前記複数の時点における前記被検体の画像のローデータ（ｒａｗ  ｄａｔａ）を取得するスキャン部と、前記複数の時点の各々を対象として、対象の時点における前記被検体の画像のローデータに基づいて、該画像を再構成すべきであるか否かを判定する判定部と、前記判定部により再構成すべきであると判定された画像のローデータに基づいて該画像を再構成する再構成部とを備える磁気共鳴イメージング装置を提供する。
【００１０】
第２の観点では、本発明は、前記判定部が、前記対象の時点における前記被検体の画像のローデータと、前記被検体の基準画像のローデータとの相違に基づいて判定する上記第１の観点の磁気共鳴イメージング装置を提供する。
【００１１】
第３の観点では、本発明は、前記判定部が、両ローデータに設定された関心領域における相違に基づいて判定する上記第２の観点の磁気共鳴イメージング装置を提供する。
【００１２】
第４の観点では、本発明は、造影剤が注入された被検体を複数の時点にてＭＲスキャンすることにより前記複数の時点における前記被検体の画像のローデータを取得するスキャン部と、前記取得されたローデータに基づいて前記複数の時点における前記被検体の画像を再構成する再構成部と、前記複数の時点の各々を対象として、対象の時点における前記被検体の画像に基づいて該画像を出力すべきであるか否かを判定する判定部とを備える磁気共鳴イメージング装置を提供する。
【００１３】
第５の観点では、本発明は、前記判定部が、前記対象の時点における前記被検体の画像と、前記被検体の基準画像との相違に基づいて判定する上記第４の観点の磁気共鳴イメージング装置を提供する。
【００１４】
第６の観点では、本発明は、前記判定部が、両画像に設定された関心領域における相違に基づいて判定する上記第５の観点の磁気共鳴イメージング装置を提供する。
【００１５】
第７の観点では、本発明は、前記判定部が、前記相違を表す指標値を算出する指標値算出部と、前記指標値が所定のしきい値以上であるか否かに基づいて前記判定を行うしきい値判定部とを有する上記第２の観点または第３の観点の磁気共鳴イメージング装置を提供する。
【００１６】
第８の観点では、本発明は、前記指標値が、前記相違を形成する一方のデータ（ｄａｔａ）を構成するデータ値の総和と、該相違を形成する他方のデータを構成するデータ値の総和との差分である上記第７の観点の磁気共鳴イメージング装置を提供する。
【００１７】
第９の観点では、本発明は、前記判定部が、前記相違を表す指標値を算出する指標値算出部と、前記指標値が所定のしきい値以上であるか否かに基づいて前記判定を行うしきい値判定部とを有する上記第５の観点または第６の観点の磁気共鳴イメージング装置を提供する。
【００１８】
第１０の観点では、本発明は、前記指標値が、前記相違を形成する一方のデータを構成するデータ値の総和と、該相違を形成する他方のデータを構成するデータ値の総和との差分である上記第９の観点の磁気共鳴イメージング装置を提供する。
【００１９】
第１１の観点では、本発明は、前記基準画像が、非造影による前記被検体の画像である第２の観点、第３の観点および第５の観点から第１０の観点のいずれか１つの観点の磁気共鳴イメージング装置を提供する。
【００２０】
第１２の観点では、本発明は、前記基準画像が、前記対象の時点より前の該時点の近傍の時点における前記被検体の画像である上記第２の観点、第３の観点および第５の観点から第１０の観点のいずれか１つの観点の磁気共鳴イメージング装置を提供する。
【００２１】
第１３の観点では、本発明は、前記被検体の画像が、前記被検体に設定された複数のスライスの全部または一部に対応する画像である上記第１の観点から第１２の観点のいずれか１つの観点の磁気共鳴イメージング装置を提供する。
【００２２】
ここで、「ＭＲスキャン」とは、静磁場空間において被検体にＲＦパルス電磁波を照射することにより磁気共鳴信号または当該信号を表すローデータ（生データ）を取得する処理であり、取得したローデータに基づいて画像を再構成する処理は含まない。
【発明の効果】
【００２３】
本発明によれば、造影ダイナミックスキャンにおける各スキャンが実施された複数の時点の各々について、その時点における被検体の画像に対応する再構成前のローデータに基づいて、その画像を再構成すべきであるか否か判定するので、再構成可能な複数の時点における画像のうち不要な画像を除いて再構成することができる。これにより、不要な画像の再構成に掛かる処理時間や必要な画像の検索に要する時間を低減して、検査時間や診断時間を短縮することができ、診断効率のよい造影ダイナミックスキャン法による撮像が可能となる。
【発明を実施するための最良の形態】
【００２４】
以下より、本発明に係る実施形態の一例について図面を参照して説明する。
【００２５】
（第１の実施形態）
図１は、第１の実施形態による磁気共鳴イメージング装置１の構成を示す構成図である。
【００２６】
図１に示すように、磁気共鳴イメージング装置１は、スキャン部２と、操作コンソール（ｃｏｎｓｏｌｅ）部３とを有している。ここで、スキャン部２は、静磁場マグネット（ｍａｇｎｅｔ）部１２と、勾配コイル（ｃｏｉｌ）部１３と、ＲＦコイル部１４と、撮像テーブル（ｔａｂｌｅ）１５とを有する。そして、操作コンソール部３は、ＲＦ駆動部２２と、勾配駆動部２３と、データ収集部２４と、制御部３０と、記憶部３１と、操作部３２と、データ処理部３３と、表示部３４とを有する。
【００２７】
スキャン部２について説明する。スキャン部２は、図１に示すように、被検体４０における撮像領域が収容される静磁場空間１１を含んでいる。そして、スキャン部２は、操作コンソール部３からの制御信号に基づいて、その静磁場が形成される静磁場空間１１に収容した被検体４０の撮像領域にＲＦパルス電磁波を照射し、その撮像領域から磁気共鳴信号を取得するスキャンを実施する。
【００２８】
スキャン部２の各構成要素について、順次、説明する。
【００２９】
静磁場マグネット部１２は、超伝導磁石（図示なし）を有しており、被検体４０の撮像領域が収容される撮像空間１１に静磁場を形成する。ここでは、静磁場マグネット部１２は、水平磁場型であって、被検体４０の体軸方向（ｚ方向）に沿うように静磁場を形成する。なお、静磁場マグネット部１２は、垂直磁場型であってもよく、また、永久磁石により構成されていてもよい。
【００３０】
勾配コイル部１３は、図１に示すように、被検体４０の撮像領域を囲むように撮像空間１１に設けられている。勾配コイル部１３は、勾配コイル（図示なし）を有しており、勾配コイルに駆動電圧が印加されることにより、静磁場が形成された撮像空間１１に勾配磁場を形成する。これにより、ＲＦコイル部１４が受信する磁気共鳴信号に位置情報を付加する。ここでは、勾配コイル部１３は、３系統の勾配コイルからなり、撮像条件に応じて、周波数エンコード（ｅｎｃｏｄｅ）方向、位相エンコード方向、およびスライス（ｓｌｉｃｅ）選択方向のそれぞれに勾配磁場を形成する。具体的には、勾配コイル部１３は、被検体４０のスライス選択方向に勾配磁場を形成し、ＲＦコイル部１４がＲＦパルス電磁波を送信することによって励起させる被検体４０のスライスを選択する。また、勾配コイル部１３は、被検体４０の位相エンコード方向に勾配磁場を形成し、ＲＦパルス電磁波により励起されたスライスからの磁気共鳴信号を位相エンコードする。そして、勾配コイル部１３は、被検体４０の周波数エンコード方向に勾配磁場を形成し、ＲＦパルス電磁波により励起されたスライスからの磁気共鳴信号を周波数エンコードする。
【００３１】
ＲＦコイル部１４は、図１に示すように、被検体４０の撮像領域を囲むように撮像空間１１に設けられている。ＲＦコイル部１４は、ＲＦコイル（図示なし）を有しており、ＲＦコイルに駆動電圧が印加されることにより、静磁場が形成された撮像空間１１内にＲＦパルス電磁波を送信して高周波磁場を形成する。これにより、撮像空間１１内の被検体４０の撮像領域におけるプロトンのスピンを励起する。そして、ＲＦコイル部１４は、その励起された被検体４０内のプロトンから発生する電磁波を磁気共鳴信号として受信する。
【００３２】
撮像テーブル１５は、被検体４０を載置するクレードル（ｃｒａｄｌｅ）を有する。撮像テーブル１５は、制御部３０からの制御信号に基づいてクレードルを移動させ、クレードルに載置された被検体４０を撮像空間１１の内部と外部との間で移動させる。
【００３３】
操作コンソール部３について説明する。操作コンソール部３は、スキャン部２が被検体４０についてＭＲスキャン（以下、単にスキャンという）を実施するように制御し、そのスキャン部２が実施したスキャンによって得られた磁気共鳴信号に基づいて、被検体４０の画像を再構成すると共に、その再構成した画像を表示する。
【００３４】
操作コンソール部３を構成する各部について、順次、説明する。
【００３５】
ＲＦ駆動部２２は、ＲＦ発振器（図示なし）、ゲート（ｇａｔｅ）変調器（図示なし）、およびＲＦ電力増幅器（図示なし）を有している。ＲＦ駆動部２２は、制御部３０からの制御信号に基づいてこれらを制御することによりＲＦコイル部１４に駆動電圧を印加し、撮像空間１１内に高周波磁場を形成させる。すなわち、ＲＦ駆動部２２は、ゲート変調器を用いてＲＦ発振器からのＲＦ信号を所定のタイミングおよび所定の包絡線の信号に変調する。その後、このゲート変調器により変調された信号をＲＦ電力増幅器により増幅し、この増幅された信号を駆動電圧としてＲＦコイル部１４に印加する。これにより、撮像空間１１内にＲＦパルス電磁波を送信させ、高周波磁場を形成させる。
【００３６】
勾配駆動部２３は、勾配コイル部１３の３系統の勾配コイルに対応して３系統の駆動回路（図示なし）を有している。勾配駆動部２３は、制御部３０からの制御信号に基づいてこれらの駆動回路を制御し、勾配コイルに駆動電圧である勾配パルスを印加する。これにより、撮像空間１１内に勾配磁場を形成させる。
【００３７】
データ収集部２４は、位相検波器（図示なし）と、アナログ（ａｎａｌｏｇ）／デジタル（ｄｉｇｉｔａｌ）変換器（図示なし）とを有する。データ収集部２４は、制御部３０からの制御信号に基づいて位相検波器およびアナログ／デジタル変換器を制御することにより、ＲＦコイル部１４が受信する磁気共鳴信号を収集する。すなわち、データ収集部２４は、ＲＦコイル部１４からの磁気共鳴信号を、ＲＦ駆動部２２のＲＦ発振器の出力を参照信号として位相検波器によって位相検波し、位相検波されたアナログ信号である磁気共鳴信号を、アナログ／デジタル変換器によってローデータに変換する。ローデータは、後述の記憶部３１、データ処理部３３へ出力される。
【００３８】
制御部３０は、コンピュータ（ｃｏｍｐｕｔｅｒ）と、コンピュータを用いて所定のスキャンに対応する動作を各部に実行させるプログラム（ｐｒｏｇｒａｍ）とを有する。そして、制御部３０は、後述する操作部３２に接続されており、操作部３２に入力された操作信号を処理し、撮像テーブル１５とＲＦ駆動部２２と勾配駆動部２３とデータ収集部２４との各部に、制御信号を出力し制御を行う。また、制御部３０は、所望の画像を得るために、操作部３２からの操作信号に基づいてデータ処理部３３、表示部３４を制御する。
【００３９】
本実施形態において、制御部３０は、ＲＦ駆動部２２と勾配駆動部２３とデータ収集部２４とを制御して、造影剤が注入された被検体４０の撮像領域を複数の時点にてスキャンする造影ダイナミックスキャンを実施する。なお、ここでは、造影ダイナミックスキャンを実施する方法として、スキャンする時点を特に決めず、造影剤を被検体に注入した後、被検体４０の撮像領域において造影剤が所望の状態になる時点を含むと考えられる十分長い期間について、略定期的にスキャンを繰り返す方法を採用する。
【００４０】
記憶部３１は、コンピュータとコンピュータに所定のデータ処理を実行させるプログラムとを有する。そして、記憶部３１は、データ収集部２４により収集された磁気共鳴信号であるローデータ（イメージングデータあるいはＫ空間データとも言う）、後述するデータ処理部３３によりローデータをフーリエ変換して得られる画像データ等を記憶する。
【００４１】
操作部３２は、キーボード（ｋｅｙｂｏａｒｄ）やマウス（ｍｏｕｓｅ）などの操作デバイス（ｄｅｖｉｃｅ）により構成されている。操作部３２は、オペレータ（ｏｐｅｒａｔｏｒ）によって操作データ、撮像プロトコル（ｐｒｏｔｏｃｏｌ）、撮像領域やそのスライスの設定に関するデータなどが入力され、これらのデータを制御部３０に出力する。
【００４２】
データ処理部３３は、コンピュータと、そのコンピュータを用いて所定のデータ処理を実行するプログラムを記憶するメモリ（ｍｅｍｏｒｙ）とを有しており、制御部３０からの制御信号に基づいて、データ処理を実行する。
【００４３】
図２は、データ処理部３３の構成を示すブロック（ｂｌｏｃｋ）図である。
【００４４】
図２に示すように、データ処理部３３は、指標値算出部３３１と、しきい値判定部３３２と、画像再構成部３３３と、画像出力部３３４とを有する。
【００４５】
指標値算出部３３１は、造影ダイナミックスキャンにおける個々のスキャンが実施された複数の時点の各々を処理対象として、対象の時点（以下、対象時点という）における被検体４０の画像と、被検体４０の基準画像との相違を表す指標値を算出する。具体的には、対象の時点における撮像領域に設定された所定のスライスに対応する画像のローデータと、上記所定のスライスに対応する基準画像のローデータとの相違を表す指標値を算出する。
【００４６】
所定のスライスは、撮像領域に単一のスライスが設定される場合にはその単一のスライス、撮像領域に複数のスライスが設定される場合にはその複数のスライスの全部または一部を考えることができる。ここでは、その一例として、撮像領域に設定された複数のスライスのうち操作者により指定される１または２以上のスライスを所定のスライスとする。
【００４７】
基準画像は、撮像したい所望の時相に応じて種々の画像を考えることができる。ここでは、その一例として、非造影による被検体４０の画像を基準画像とする。具体的には、その一例として、造影ダイナミックスキャンを開始する前に所定のプレスキャン（ｐｒｅ−ｓｃａｎ）を実施するものとし、そのプレスキャンが実施される時点ｔｐにおける上記所定のスライスに対応する画像を基準画像とする。なお、注入した造影剤が被検体４０の撮像領域に到達するより前に造影ダイナミックスキャンを開始することが前提であれば、造影ダイナミックスキャンにおいて最初のスキャンが実施される時点における所定のスライスに対応する画像等を基準画像としてもよい。
【００４８】
指標値は、上記の相違を表すものであれば、基本的にどのようなものでも適用可能である。例えば、指標値として、上記の相違を形成する一方のデータを構成するデータ値の平均値、中間値、総和、分散等の代表値と、上記相違を形成する他方のデータを構成するデータ値の同代表値との差分や比等を考えることができる。ここでは、その一例として、上記相違を形成する一方のデータを構成するデータ値の総和と、上記相違を形成する他方のデータを構成するデータ値の総和との差分とする。すなわち、対象時点における所定のスライスに対応する画像のローデータを構成するデータ値の総和と、基準画像のローデータを構成するデータ値の総和との差分を指標値とする。
【００４９】
なお、所定のスライスが複数のスライスである場合には、複数のスライスにそれぞれ対応する複数の画像のローデータを構成するすべてのデータ値の総和と、複数のスライスにそれぞれ対応する複数の基準画像のローデータを構成するすべてのデータ値の総和との差分を指標値とする。あるいは、同一のスライスに対応する対象時点における画像と基準画像との組合せごとに、対象時点における所定のスライスに対応する画像のローデータを構成するデータ値の総和と基準画像のローデータを構成するデータ値の総和との差分を順次算出し、算出された差分の中で最大のものを指標値とする。
【００５０】
しきい値判定部３３２は、指標値算出部３３１により算出された指標値が、所定のしきい値以上であるか否かを判定する。そして、指標値がしきい値以上であると判定されたときの対象時点における所定のスライスに対応する画像を再構成すべきであると判定する。このとき、所定のスライスが複数のスライスである場合には、これら複数のスライスに対応する複数の画像すべてを再構成すべきであると判定する。
【００５１】
これにより、被検体４０の撮像領域に造影剤が到達して、撮像領域に対応する画像のローデータのデータ値が通常より変化したと考えられる時相を検出し、このような時相における画像を再構成すべきであると判定する。
【００５２】
なお、しきい値は、所定のスライスとして指定されたスライスの数、関心領域が設定された場合においてはその関心領域のサイズなどに応じて変化させる。例えば、スライスの数が大きくなるほど大きく、関心領域のサイズが大きくなるほど大きくなるようしきい値を変化させる。
【００５３】
画像再構成部３３３は、しきい値判定部３３２により再構成すべきであると判定された画像のみを、そのローデータに基づいて再構成する。
【００５４】
画像出力部３３４は、画像再構成部３３３により再構成された画像を出力する。出力先としては、記憶部３１、表示部３４、磁気共鳴イメージング装置１とネットワーク（ｎｅｔｗｏｒｋ）等により接続された画像データサーバ（ｄａｔａ  ｓｅｒｖｅｒ）（不図示）や、ビューア（ｖｉｅｗｅｒ）端末（不図示）等を考えることができる。
【００５５】
これより、本実施形態による磁気共鳴イメージング装置の動作について説明する。
【００５６】
図３は、本実施形態による磁気共鳴イメージング装置の動作の流れを示すフロー（ｆｌｏｗ）図である。また、図４は、本装置により被検体４０をスキャンする際のシーケンス（ｓｅｑｕｅｎｃｅ）を示す図である。
【００５７】
まず、図３に示すように、プレスキャンを実施する（ＳＴ１０）。
【００５８】
はじめに、撮像領域を設定する。例えば、表示部３４が位置合せ用に撮像した被検体４０のコロナル（ｃｏｒｏｎａｌ）画像を表示する。そして、オペレータが操作部３２を操作することにより、表示されたコロナル画像上において撮像領域を設定する。また、オペレータは操作部３２を操作することにより、設定された撮像領域内に画像化するスライスを設定する。ここでは、図５に示すように、被検体４０の撮像領域に複数のスライスＳ1，Ｓ2，…，Ｓnを設定する。
【００５９】
その後、スキャン部２が、図４に示すように、プレスキャンとしてのスキャンを実施する。すなわち、造影剤の注入を開始する時点ｔ0より前の時点ｔpにおいてスキャンを実施する。スキャンは、例えば、特開２００７−１１１１８８号公報に記載されているようなイメージングシーケンスを実施することにより実施される。なお、１回のスキャンを開始してから終了するまでには所定時間を要するが、便宜上、１の時点に対応して１回のスキャンが実施されるものとして説明する。また、ここでは、プレスキャンとしてスキャンを１回だけ実施しているが、プレスキャンとしてスキャンを複数回、例えば数回程度、繰返し実施してもよい。
【００６０】
次に、図３に示すように、造影ダイナミックスキャンを実施する（ＳＴ２０）。
【００６１】
はじめに、造影ダイナミックスキャンにおける各種の設定を行う。例えば、オペレータが操作部３２を操作することにより、表示された撮像プロトコル上で、スキャンの繰返し周期や繰返し回数などを設定する。なお、このとき、スキャンの繰返し周期、繰返し回数は、造影ダイナミックスキャンが実施される期間内に所望の時相が含まれるよう設定される。
【００６２】
その後、図４に示すように、造影剤の被検体４０への注入を行う。すなわち、オペレータがインジェクタ（ｉｎｊｅｃｔｏｒ）（図示なし）を用いて、時点ｔ0から所定時間が経過するまでの間において造影剤を注入する。
【００６３】
そして、スキャン部２が、図４に示すように、造影ダイナミックスキャンを実施する。すなわち、造影剤の被検体４０への注入が終了した後、複数の時点ｔ1，ｔ2，…，ｔNにおいてスキャンを実施する。これにより、複数の時点ｔ1，ｔ2，…，ｔNの各々におけるスライスＳ1，Ｓ2，…，Ｓnに対応する画像のローデータを取得する。
【００６４】
次に、画像再構成処理を実行する（ＳＴ３０）。なお、この画像再構成処理は、ステップ（ｓｔｅｐ）ＳＴ２０における造影ダイナミックスキャンと略並行して行うことができる。
【００６５】
図６は、ステップＳＴ３０における画像再構成処理のフローを示す図である。
【００６６】
まず、被検体４０の撮像領域に設定された全スライスの中から目的のスライスを指定する（ＳＴ３０１）。例えば、オペレータが、操作部３２を操作することにより、スライスＳ1，Ｓ2，…，Ｓnの中から１または２以上のスライスを指定する（以下、この指定されたスライスを指定スライスという）。
【００６７】
次に、対象時点を選択する（ＳＴ３０２）。すなわち、指標値算出部３３１が、造影ダイナミックスキャンにおける個々のスキャンが実施された複数の時点ｔ1，ｔ2，…，ｔNの中から、処理対象となる時点として未だ選択されていない１つの時点を対象時点ｔqとして選択する。
【００６８】
次に、対象時点における指定スライスに対応する指標値を算出する（ＳＴ３０３）。
【００６９】
図７は、所定のスライスに対応する画像のローデータの概念図である。また、図８はローデータに基づいて指標値を算出する処理の概念図である。
【００７０】
所定のスライスに対応する画像のローデータは、図７に示すように、磁気共鳴信号をデジタル信号に変換して得られるデータであり、時点ｔiにおけるスライスＳjに対応する画像のローデータは、ＲＤ（ｔi，Ｓj，ｄk）（ｋ＝１，２，…，ＮＲ）として表される。なお、図７に示すローデータは、いわゆるＫ空間データであり、本図は、フェーズエンコードを変えながら横軸方向に１ビューずつローデータを収集した例を示しており、波高は各データ値の大きさを表現したものである。
【００７１】
指標値算出部３３１は、対象時点ｔqにおける指定スライスに対応する画像のローデータと、指定スライスに対応する基準画像（ここではプレスキャンが実施された時点ｔpにおける指定スライスに対応する画像）のローデータとの相違を表す指標値を算出する。
【００７２】
具体的には、例えば、指定スライスが単一のスライスＳaであれば、指標値算出部３３１は、図８（ａ）に示すように、対象時点ｔqにおける指定スライスＳaに対応する画像のローデータＲＤ（ｔq，Ｓa，ｄk）を構成するデータ値の絶対値の総和と、指定スライスＳaに対応する基準画像のローデータＲＤ（ｔp，Ｓa，ｄk）を構成するデータ値の絶対値の総和との差分の絶対値を指標値Ｃ（ｔq）として算出する。これを式で表すと次式のように表すことができる。
【数１】
2010104459.tif 000003
また、例えば、指定スライスが撮像領域に設定された全スライスの一部である２以上のスライスＳa，Ｓa+1，…，Ｓa+mであれば、指標値算出部３３１は、図８（ｂ）に示すように、対象時点ｔqにおける指定スライスＳa，Ｓa+1，…，Ｓa+mに対応する画像のローデータＲＤ（ｔq，Ｓa，ｄk），ＲＤ（ｔq，Ｓa+1，ｄk），…，ＲＤ（ｔq，Ｓa+m，ｄk）を構成するデータ値の絶対値の総和と、指定スライスＳa，Ｓa+1，…，Ｓa+mに対応する基準画像のローデータＲＤ（ｔp，Ｓa，ｄk），ＲＤ（ｔp，Ｓa+1，ｄk），…，ＲＤ（ｔp，Ｓa+m，ｄk）を構成するデータ値の絶対値の総和との差分の絶対値を、指標値Ｃ（ｔq）として算出する。これを式で表すと次式のように表すことができる。
【数２】
2010104459.tif 000004
また、例えば、指定スライスが撮像領域に設定された全スライスＳ1，Ｓ2，…，Ｓnであれば、指標値算出部３３１は、図８（ｃ）に示すように、対象時点ｔqにおける指定スライスＳ1，Ｓ2，…，Ｓnに対応する画像のローデータＲＤ（ｔq，Ｓ1，ｄk），ＲＤ（ｔq，Ｓ2，ｄk），…，ＲＤ（ｔq，Ｓn，ｄk）を構成するデータ値の絶対値の総和と、指定スライスＳ1，Ｓ2，…，Ｓnに対応する基準画像のローデータＲＤ（ｔp，Ｓ1，ｄk），ＲＤ（ｔp，Ｓ2，ｄk），…，ＲＤ（ｔp，Ｓn，ｄk）を構成するデータ値の絶対値の総和との差分の絶対値を、指標値Ｃ（ｔq）として算出する。これを式で表すと次式のように表すことができる。
【数３】
2010104459.tif 000005
次に、しきい値判定処理を行う（ＳＴ３０４）。すなわち、しきい値判定部３３２が、ステップＳＴ３０３にて算出された指標値Ｃ（ｔp）が所定のしきい値ｔｈ1以上であるか否かを判定する。ここで、指標値Ｃ（ｔp）がしきい値ｔｈ1以上であると判定された場合には、対象時点における指定スライスに対応する画像のローデータと基準画像のローデータとの間に一定レベル（ｌｅｖｅｌ）以上の相違が存在すると考え、画像再構成部３３３が、この対象時点ｔqにおける指定スライスに対応する画像を、そのローデータに基づいて再構成する（ＳＴ３０５）。一方、指標値Ｃ（ｔp）がしきい値ｔｈ1未満であると判定された場合には、対象時点における指定スライスに対応する画像のローデータと基準画像のローデータとの間に一定レベル以上の相違がないと考え、画像再構成部３３３は、この対象時点ｔqにおける指定スライスに対応する画像は再構成しない。
【００７３】
次に、次期の対象時点として選択すべき時点があるか否かを判定する（ＳＴ３０６）。すなわち、指標値算出部３３１が、複数の時点ｔ1，ｔ2，…，ｔNの中に、対象時点として未だ選択されていない時点が存在するか否かを判定する。ここで、未だ選択されていない時点が存在すると判定された場合には、ステップＳＴ３０２に戻り、指標値算出部３３１が、その未だ選択されていない時点を新たな対象時点として選択し、画像再構成処理を継続する。一方、未だ選択されていない時点が存在しないと判定された場合には、画像再構成処理を終了する。
【００７４】
なお、本実施形態では、対象時点における指定スライスに対応する画像を再構成すべきか否かの判定結果が出る度に、その対象時点における画像を逐次的に再構成しているが、もちろん、先にすべての対象時点における画像について再構成すべきか否かの判定を行い、その後、再構成すべきであると判定された画像だけをまとめて再構成するようにしてもよい。
【００７５】
次に、画像出力処理を実行する（ＳＴ４０）。すなわち、画像出力部３３４が、ステップＳＴ３０において画像再構成部３３３により再構成された画像を記憶部３１、表示部３４、磁気共鳴イメージング装置１とネットワーク等により接続された画像データサーバ（不図示）や、ビューア端末（不図示）等に出力する。
【００７６】
ここで、本実施形態による造影ダイナミックスキャンを実施した場合に得られる結果の一例について説明する。
【００７７】
図９は、本実施形態による造影ダイナミックスキャンの各スキャンが実施された各時点における指標値および画像再構成要否判定結果の一例を示す図である。本実施形態による造影ダイナミックスキャンを実施した場合、例えば、図９に示すような結果が得られる。図９の上段のグラフ（ｇｒａｐｈ）は、横軸を時間ｔ、縦軸を指標値Ｃ（ｔi）として、造影ダイナミックスキャンにおける個々のスキャンが実施された各時点ｔi（ｉ＝１，２，…，５０）に対応する指標値Ｃ（ｔi）を実線の変化曲線Ｃ１で示したものである。なお、破線の変化曲線Ｚ１は、プレスキャンを実施した時点ｔpと各時点ｔiにおける指定スライスＳaに対応する画像のローデータＲＤ（ｔi，Ｓa）を構成するデータ値の絶対値の総和ΣＲＤ（ｔi，Ｓa）の変化を示しており、縦軸は指標値Ｃ（ｔi）と共通である。また、図９の下段のチャートＲ１は、各時点における指定スライスＳaに対応する画像に対する再構成すべきであるか否かの判定結果を示している。
【００７８】
この例では、指定スライスＳaに対応する画像のローデータに関する変化曲線Ｚ１と指標値の変化曲線Ｃ１の形状から、次のような状況が想定される。まず、時点ｔ1からｔ8付近では、造影剤が撮像領域に到達しておらず、時点ｔ9からｔ20付近では、造影剤が撮像領域の動脈に到達してピーク（ｐｅａｋ）を迎える。また、時点ｔ20からｔ31付近では、撮像領域の動脈を通過している造影剤が抜け始めるとともに、撮像領域の静脈に造影剤が到達し始める。時点ｔ32からｔ35付近では、造影剤が撮像領域の静脈に到達してピークを迎える。そして、時点ｔ36から時点ｔ39付近では、撮像領域の静脈を通過していた造影剤も抜け始め、時点ｔ40からｔ50では、造影剤が撮像領域から略抜けている。
【００７９】
指標値Ｃ（ｔi）がしきい値ｔｈ１以上となる時点は、時点ｔ9からｔ39までである。したがって、これらの時点ｔ9からｔ39における指定スライスＳaに対応する画像に対しては、再構成すべきであるとの判定がなされ、一方、時点ｔ1からｔ8までと時点ｔ40からｔ50までの各時点における指定スライスＳaに対応する画像に対しては、再構成すべきでないとの判定がなされる。これにより、造影剤が撮像領域に到達する前の時点と、造影剤が撮像領域から抜けた後の時点における不要な画像は再構成されず、造影剤が撮像領域に到達してから抜けるまでの間の時点における画像のみが再構成される。
【００８０】
（第２の実施形態）
本実施形態は、第１の実施形態と略同様の構成による磁気共鳴イメージング装置であるが、次の点で異なる。
【００８１】
本実施形態では、対象時点における指定スライスに対応する画像と比較する基準画像を、対象時点より前であり当該時点の近傍の時点における指定スライスに対応する画像とする。具体的には、その一例として、造影ダイナミックスキャンにおける個々のスキャンが実施された複数の時点のうち、対象時点より１つ前の時点における指定スライスに対応する画像を基準画像とする。
【００８２】
つまり、指標値算出部３３１は、例えば、指定スライスが単一のスライスＳaであれば、対象時点ｔqにおける指定スライスＳaに対応する画像のローデータＲＤ（ｔq，Ｓa，ｄk）を構成するデータ値の絶対値の総和と、指定スライスＳaに対応する基準画像のローデータＲＤ（ｔq-1，Ｓa，ｄk）を構成するデータ値の絶対値の総和との差分の絶対値を指標値Ｃ（ｔq）として算出する。これを式で表すと次式のように表すことができる。
【数４】
2010104459.tif 000006
また、例えば、指定スライスが撮像領域に設定された全スライスの一部である２以上のスライスＳa，Ｓa+1，…，Ｓa+mであれば、指標値算出部３３１は、対象時点ｔqにおける指定スライスＳa，Ｓa+1，…，Ｓa+mに対応する画像のローデータＲＤ（ｔq，Ｓa，ｄk），ＲＤ（ｔq，Ｓa+1，ｄk），…，ＲＤ（ｔq，Ｓa+m，ｄk）を構成するデータ値の絶対値の総和と、指定スライスＳa，Ｓa+1，…，Ｓa+mに対応する基準画像のローデータＲＤ（ｔq-1，Ｓa，ｄk），ＲＤ（ｔq-1，Ｓa+1，ｄk），…，ＲＤ（ｔq-1，Ｓa+m，ｄk）を構成するデータ値の絶対値の総和との差分の絶対値を、指標値Ｃ（ｔq）として算出する。これを式で表すと次式のように表すことができる。
【数５】
2010104459.tif 000007
また、例えば、指定スライスが撮像領域に設定された全スライスＳ1，Ｓ2，…，Ｓnであれば、指標値算出部３３１は、対象時点ｔqにおける指定スライスＳ1，Ｓ2，…，Ｓnに対応する画像のローデータＲＤ（ｔq，Ｓ1，ｄk），ＲＤ（ｔq，Ｓ2，ｄk），…，ＲＤ（ｔq，Ｓn，ｄk）を構成するデータ値の絶対値の総和と、指定スライスＳ1，Ｓ2，…，Ｓnに対応する基準画像のローデータＲＤ（ｔq-1，Ｓ1，ｄk），ＲＤ（ｔq-1，Ｓ2，ｄk），…，ＲＤ（ｔq-1，Ｓn，ｄk）を構成するデータ値の絶対値の総和との差分の絶対値を、指標値Ｃ（ｔq）として算出する。これを式で表すと次式のように表すことができる。
【数６】
2010104459.tif 000008
また、しきい値判定部３３２は、算出された指標値Ｃ（ｔp）が所定のしきい値ｔｈ2以上であるか否かを判定する。そして、指標値Ｃ（ｔp）がしきい値ｔｈ2以上であると判定された場合には、対象時点における指定スライスＳaに対応する画像のローデータと基準画像のローデータとの間に一定レベル以上の相違が存在すると考えることができ、画像再構成部３３３が、この対象時点ｔqにおける指定スライスに対応する画像を、そのローデータに基づいて再構成する。
【００８３】
これにより、被検体４０の撮像領域における造影剤の濃度等の状態が変化して、撮像領域に対応する画像のローデータのデータ値が所定時間内に一定レベル以上変化したと考えられる時相を検出し、このような時相における画像のみを再構成する。
【００８４】
ここで、本実施形態による造影ダイナミックスキャンを実施した場合に得られる結果の一例について説明する。
【００８５】
図１０は、本実施形態による造影ダイナミックスキャンの各スキャンが実施された各時点における指標値および画像再構成要否判定結果の一例を示す図である。本実施形態による造影ダイナミックスキャンを実施した場合、例えば、図１０に示すような結果が得られる。図１０の上段のグラフは、横軸を時間ｔ、縦軸を指標値Ｃ（ｔi）として、造影ダイナミックスキャンにおける個々のスキャンが実施された各時点ｔi（ｉ＝１，２，…，５０）に対応する指標値Ｃ（ｔi）を実線の変化曲線Ｃ２で示したものである。なお、破線の変化曲線Ｚ２は、プレスキャンを実施した時点ｔpと各時点ｔiにおける指定スライスＳaに対応する画像のローデータＲＤ（ｔi，Ｓa）を構成するデータ値の絶対値の総和ΣＲＤ（ｔi，Ｓa）の変化を示しており、縦軸は指標値Ｃ（ｔi）と共通である。また、図１０の下段のチャートＲ２は、各時点における指定スライスＳaに対応する画像に対する再構成すべきであるか否かの判定結果を示している。
【００８６】
この例では、対象時点における指定スライスＳaに対応する画像のローデータに関する変化曲線Ｚ２と指標値の変化曲線Ｃ２の形状から、次のような状況が想定される。まず、時点ｔ1からｔ8付近では、造影剤が撮像領域に到達しておらず、時点ｔ9からｔ20付近では、造影剤が撮像領域の動脈に到達してピークを迎える。また、時点ｔ20からｔ31付近では、撮像領域の動脈を通過している造影剤が抜け始めるとともに、撮像領域の静脈に造影剤が到達し始める。時点ｔ32からｔ35付近では、造影剤が撮像領域の静脈に到達してピークを迎える。そして、時点ｔ36から時点ｔ39付近では、撮像領域の静脈を通過していた造影剤も抜け始め、時点ｔ40からｔ50では、造影剤が撮像領域から略抜けている。
【００８７】
指標値Ｃ（ｔi）がしきい値ｔｈ2以上となる時点は、時点ｔ9からｔ11までと、時点ｔ18からｔ21までと、時点ｔ31からｔ33までと、時点ｔ35からｔ39までである。したがって、これらの時点における画像に対しては、再構成すべきであるとの判定がなされ、一方、残りの各時点における画像に対しては、再構成すべきでないとの判定がなされる。これにより、撮像領域における造影剤の状態の変化が乏しい時点における不要な画像は再構成されず、撮像領域における造影剤の状態の変化が大きい時点における画像のみが再構成される。
【００８８】
以上のように、上記第１および第２の実施形態によれば、造影ダイナミックスキャンにおける各スキャンが実施された複数の時点の各々について、その時点における被検体の画像に対応する再構成前のローデータに基づいて、その画像を再構成すべきであるか否か判定するので、再構成可能な複数の時点における画像のうち不要な画像を除いて再構成することができる。これにより、不要な画像の再構成に掛かる処理時間や必要な画像の検索に要する時間を低減して、検査時間や診断時間を短縮することができ、診断効率のよい造影ダイナミックスキャン法による撮像が可能となる。
【００８９】
また、上記第１および第２の実施形態によれば、再構成される不要な画像が減少するので、画像を保存するデータサーバの保存スペースを無駄に消費したり、ビューア端末のメモリスペース（ｍｅｍｏｒｙ  ｓｐａｃｅ）を無駄に占有したりすることを防ぐことができ、画像のデータを扱うインフラ（ｉｎｆｒａｓｔｒｕｃｔｕｒｅ）に対する負荷を軽減することができる。
【００９０】
（第３の実施形態）
本実施形態は、第１の実施形態と略同様の構成による磁気共鳴イメージング装置であるが、データ処理部３３の各構成要素の機能が次のように異なる。
【００９１】
画像再構成部３３３は、プレスキャンおよび造影ダイナミックスキャンにおける個々のスキャンを実施して得られたローデータに基づいて、複数の時点の各々における被検体４０の撮像領域に設定された所定のスライスに対応する画像を再構成する。
【００９２】
指標値算出部３３１は、造影ダイナミックスキャンにおける個々のスキャンが実施された複数の時点の各々を処理対象として、対象の時点（以下、対象時点という）における所定のスライスに対応する画像と基準画像との相違を表す指標値を算出する。ここでは、対象時点における所定のスライスに対応する画像の画像データを構成するデータ値の総和と、基準画像の画像データを構成するデータ値の総和との差分を指標値とする。
【００９３】
なお、所定のスライスが複数のスライスである場合には、複数のスライスにそれぞれ対応する複数の画像の画像データを構成するすべてのデータ値の総和と、複数のスライスにそれぞれ対応する複数の基準画像の再構成画像データを構成するすべてのデータ値の総和との差分を指標値とする。あるいは、同一のスライスに対応する対象時点における画像と基準画像との組合せごとに、対象時点における画像の画像データを構成するデータ値の総和と基準画像の画像データを構成するデータ値の総和との差分を順次算出し、算出された各スライスに対応する差分の中で最大のものを指標値とする。
【００９４】
しきい値判定部３３３は、指標値算出部３３１により算出された指標値が、所定のしきい値以上であるか否かを判定する。そして、指標値がしきい値以上であると判定された対象時点の画像を出力すべきであると判定する。このとき、所定のスライスが複数のスライスである場合には、これら複数のスライスに対応する複数の画像すべてを出力すべきであると判定する。
【００９５】
画像出力部３３４は、画像再構成部３３３により再構成された画像のうち、しきい値判定部３３２により出力すべきであると判定された画像を出力する。出力先としては、記憶部３１、表示部３４、磁気共鳴イメージング装置１とネットワーク等により接続された画像データサーバ（不図示）や、ビューア端末（不図示）等を考えることができる。
【００９６】
これより、本実施形態よる磁気共鳴イメージング装置における動作について説明する。なお、本実施形態による磁気共鳴イメージング装置の動作の流れは、図３に示す通りである。また、本装置により被検体４０をスキャンする際のシーケンスは、図４に示す通りである。
【００９７】
まず、図３に示すように、第１の実施形態と同様のプレスキャンを実施する（ＳＴ１０）。
【００９８】
次に、図３に示すように、第１の実施形態と同様の造影ダイナミックスキャンを実施する（ＳＴ２０）。これにより、複数の時点ｔ1，ｔ2，…，ｔNの各々におけるスライスＳ1，Ｓ2，…，Ｓnに対応する画像のローデータを取得する。
【００９９】
次に、図３に示すように、画像再構成処理を実行する（ＳＴ３０）。
【０１００】
ここでは、まず、被検体４０の撮像領域に設定された全スライスの中から目的のスライスを指定する。例えば、オペレータが、操作部３２を操作することにより、スライスＳ1，Ｓ2，…，Ｓnの中から１または２以上のスライスを指定する（以下、この指定されたスライスを指定スライスという）。そして、画像再構成部３３３が、複数の時点ｔ1，ｔ2，…，ｔNの各々における指定スライスに対応する画像を、そのローデータを基に再構成する。
【０１０１】
次に、図３に示すように、画像出力処理を行う（ＳＴ４０）。
【０１０２】
図１１は、ステップＳＴ４０における画像出力処理のフローを示す図である。
【０１０３】
まず、対象時点を選択する（ＳＴ４０１）。すなわち、指標値算出部３３１が、造影ダイナミックスキャンにおける個々のスキャンが実施された複数の時点ｔ1，ｔ2，…，ｔNの中から、対象時点として未だ選択されていない１つの時点を対象時点ｔqとして選択する。
【０１０４】
次に、対象時点における指定スライスに対応する指標値を算出する（ＳＴ４０２）。
【０１０５】
図１２は、所定のスライスに対応する画像の画像データの概念図である。また、図１３は、画像データに基づいて指標値を算出する処理の概念図である。
【０１０６】
所定のスライスに対応する画像の画像データは、図１２に示すように、ローデータをフーリエ変換して得られるデータであり、時点ｔiにおけるスライスＳjに対応する画像の画像データは、ＧＤ（ｔi，Ｓj，ｄk）（ｋ＝１，２，…，ＮＧ）として表される。
【０１０７】
指標値算出部３３１は、対象時点ｔqにおける指定スライスに対応する画像の画像データと、指定スライスに対応する基準画像（ここではプレスキャンが実施された時点ｔpにおける指定スライスに対応する画像）の画像データとの相違を表す指標値を算出する。
【０１０８】
具体的には、例えば、指定スライスが単一のスライスＳaであれば、指標値算出部３３１は、図１３（ａ）に示すように、対象時点ｔqにおける指定スライスＳaに対応する画像の画像データＧＤ（ｔq，Ｓa，ｄk）を構成するデータ値の絶対値の総和と、指定スライスＳaに対応する基準画像の画像データＧＤ（ｔp，Ｓa，ｄk）を構成するデータ値の絶対値の総和との差分の絶対値を指標値Ｃ（ｔq）として算出する。これを式で表すと次式のように表すことができる。
【数７】
2010104459.tif 000009
また、例えば、指定スライスが撮像領域に設定された全スライスの一部である２以上のスライスＳa，Ｓa+1，…，Ｓa+mであれば、指標値算出部３３１は、図１３（ｂ）に示すように、対象時点ｔqにおける指定スライスＳa，Ｓa+1，…，Ｓa+mに対応する画像の画像データＧＤ（ｔq，Ｓa，ｄk），ＧＤ（ｔq，Ｓa+1，ｄk），…，ＧＤ（ｔq，Ｓa+m，ｄk）を構成するデータ値の絶対値の総和と、指定スライスＳa，Ｓa+1，…，Ｓa+mに対応する基準画像の画像データＧＤ（ｔp，Ｓa，ｄk），ＧＤ（ｔp，Ｓa+1，ｄk），…，ＧＤ（ｔp，Ｓa+m，ｄk）を構成するデータ値の絶対値の総和との差分の絶対値を、指標値Ｃ（ｔq）として算出する。これを式で表すと次式のように表すことができる。
【数８】
2010104459.tif 000010
また、例えば、指定スライスが撮像領域に設定された全スライスＳ1，Ｓ2，…，Ｓnであれば、指標値算出部３３１は、図１３（ｃ）に示すように、対象時点ｔqにおける指定スライスＳ1，Ｓ2，…，Ｓnに対応する画像の画像データＧＤ（ｔq，Ｓ1，ｄk），ＧＤ（ｔq，Ｓ2，ｄk），…，ＧＤ（ｔq，Ｓn，ｄk）を構成するデータ値の絶対値の総和と、指定スライスＳ1，Ｓ2，…，Ｓnに対応する基準画像の画像データＧＤ（ｔp，Ｓ1，ｄk），ＧＤ（ｔp，Ｓ2，ｄk），…，ＧＤ（ｔp，Ｓn，ｄk）を構成するデータ値の絶対値の総和との差分の絶対値を、指標値Ｃ（ｔq）として算出する。これを式で表すと次式のように表すことができる。
【数９】
2010104459.tif 000011
次に、しきい値判定処理を行う（ＳＴ４０３）。すなわち、しきい値判定部３３２が、ステップＳＴ４０２にて算出された指標値Ｃ（ｔp）が所定のしきい値ｔｈ2以上であるか否かを判定する。
【０１０９】
ここで、指標値Ｃ（ｔp）がしきい値ｔｈ2以上であると判定された場合には、対象時点における指定スライスに対応する画像と基準画像との間に一定レベル以上の相違が存在すると考え、画像出力部３３４が、この対象時点ｔqにおける指定スライスに対応する画像を出力する（ＳＴ４０４）。
【０１１０】
一方、指標値Ｃ（ｔp）がしきい値ｔｈ2未満であると判定された場合には、対象時点における指定スライスに対応する画像の画像データと基準画像の画像データとの間に一定レベル以上の相違がないと考え、画像出力部３３４は、この対象時点ｔqにおける指定スライスに対応する画像は出力しない。
【０１１１】
なお、再構成画像データの出力先としては、例えば、記憶部３１、表示部３４、磁気共鳴イメージング装置１とネットワーク等により接続された画像データサーバ（不図示）や、ビューア端末（不図示）等を考えることができる。
【０１１２】
また、「画像を出力する」とは、例えば、画像データを出力すること、画像を画面に表示すること、画像を表示可能な状態にしてその画像データを出力すること、画像データのデータ名を表示・検索可能な状態にしてその画像データを出力すること、等を考えることができる。
【０１１３】
次に、次期の対象時点として選択すべき時点があるか否かを判定する（ＳＴ４０５）。すなわち、指標値算出部３３１が、複数の時点ｔ1，ｔ2，…，ｔNの中に、対象時点として未だ選択されていない時点が存在するか否かを判定する。
【０１１４】
ここで、未だ選択されていない時点が存在すると判定された場合には、ステップＳＴ４０１に戻り、指標値算出部３３１が、その未だ選択されていない時点を新たな対象時点として選択し、画像出力処理を継続する。一方、未だ選択されていない時点が存在しないと判定された場合には、画像出力処理を終了する。
【０１１５】
なお、本実施形態では、対象時点における指定スライスに対応する画像を出力すべきであるか否かの判定結果が出る度に、その画像を逐次的に出力しているが、もちろん、先にすべての対象時点における指定スライスに対応する画像について出力すべきであるか否かの判定を行い、その後、出力すべきであると判定された画像だけをまとめて出力するようにしてもよい。
【０１１６】
以上のように、上記第３の実施形態によれば、造影ダイナミックスキャンにおける各スキャンが実施された複数の時点の各々について、その時点における被検体の再構成された画像に基づいて、その画像を出力すべきであるか否かを判定するので、再構成された複数の時点における画像のうち不要な画像を除いて出力することができる。これにより、必要な画像を検索するのに要する時間を低減して、検査時間や診断時間を短縮することができ、診断効率のよい造影ダイナミックスキャン法による撮像が可能となる。
【０１１７】
また、上記第３の実施形態によれば、出力される不要な画像が減少するので、画像を保存するデータサーバの保存スペースを無駄に消費したり、ビューア端末のメモリスペースを無駄に占有したりすることを防ぐことができ、画像のデータを扱うインフラに対する負荷を軽減することができる。
【０１１８】
なお、本発明は、上記の各実施形態に限定されず、種々の変形例を採用することができる。
【０１１９】
上記第１〜第３の実施形態では、造影ダイナミックスキャンを実施する方法として、スキャンする時点を特に決めず、造影剤を被検体に注入した後、被検体４０の撮像領域において造影剤が所望の状態になる時点を含むと考えられる十分長い期間について略定期的にスキャンを繰り返す方法を採用したが、これに限定されない。例えば、被検体４０の撮像領域において造影剤が所望の状態になる時点を推定し、この推定した時点とその前後幾らかの時点についてのみスキャンを実施する方法を採用してもよい。
【０１２０】
また、上記第１および第２の実施形態では、対象時点における画像のローデータと基準画像のローデータとの相違を表す指標値に基づいて、対象時点における画像を再構成すべきであるか否かを判定している。また、上記第３の実施形態では、対象時点における画像と基準画像との相違を表す指標値に基づいて、対象時点における画像を出力すべきであるか否かを判定している。しかし、このような基準画像を用いることなく判定する方法も考えられる。例えば、ローデータや画像データのデータ値を反映する所定の特徴量、例えばローデータのデータ値の総和、あるいは画像データのデータ値の総和が、所定のしきい値以上であるか否かに基づいて、対象時点における画像を再構成すべきであるか否か、あるいは対象時点における画像を出力すべきであるか否かを判定してもよい。
【０１２１】
また、上記第１〜第３の実施形態では、しきい値は、対象時点に関係なく設定しているが、対象時点に応じて変化させてもよい。例えば、対象時点が造影剤を注入してからの経過時間が比較的短い時点である場合には、しきい値を高めに設定し、この経過時間が比較的長い時点である場合には、しきい値を低めに設定する。このようにすれば、造影剤が撮像領域に対応する部位に到達する前に発生した外来ノイズ（ｎｏｉｓｅ）等による誤判定を抑えることができる。
【０１２２】
また、上記第１および第２の実施形態では、対象時点における画像のローデータ全体と基準画像のローデータ全体との相違に基づいて、対象時点における画像を再構成すべきであるか否かを判定している。また、上記第３の実施形態では、対象時点における画像全体と基準画像全体との相違に基づいて、対象時点における画像を出力すべきであるか否かを判定している。しかし、このようにデータ全体を比較するのではなく、データの一部を比較するようにしてもよい。例えば、対象時点における画像のローデータと基準画像のローデータとの間における互いに対応する所定の領域に関心領域ＲＯＩ（Ｒｅｇｉｏｎ  Ｏｆ  Ｉｎｔｅｒｅｓｔ）をそれぞれ設定し、対象時点における画像のローデータの関心領域におけるデータと基準画像のローデータの関心領域におけるデータとの相違に基づいて、対象時点における画像を再構成すべきであるか否かを判定してもよい。あるいは、対象時点における画像と基準画像との間における互いに対応する所定の領域に関心領域ＲＯＩをそれぞれ設定し、対象時点における画像の関心領域におけるデータと基準画像の関心領域におけるデータとの相違に基づいて、対象時点における画像を出力すべきか否かを判定してもよい。このようにすれば、より少ないデータを用いてそれぞれの判定結果を得ることができ、判定に掛かる処理時間の短縮を図ることができる。
【０１２３】
また、上記の各実施形態において、基準画像のローデータを構成するデータ値の総和は、近接する複数の時点の各々における所定のスライスに対応する画像のローデータを合成するなどして求められる所定のローデータを構成するデータ値の総和としてもよい。あるいは、近接する複数の時点の各々における所定のスライスに対応する画像のローデータを構成するデータ値の総和の平均値または中間値としてもよい。例えば、プレスキャンを近接する３つの時点ｔｐ1，ｔｐ2，ｔｐ3について実施する場合、時点ｔｐ１における所定のスライスに対応する画像のローデータを構成するデータ値の総和，時点ｔｐ２における所定のスライスに対応する画像のローデータを構成するデータ値の総和，および時点ｔｐ３における所定のスライスに対応する画像のローデータを構成するデータ値の総和に対する平均値または中間値としてもよい。また例えば、対象時点ｔiの前に近接する３つの時点ｔi-3，ｔi-2，ｔi-1がある場合、時点ｔi-3における所定のスライスに対応する画像のローデータを構成するデータ値の総和，時点ｔi-2における所定のスライスに対応する画像のローデータを構成するデータ値の総和，および時点ｔi-1における所定のスライスに対応する画像のローデータを構成するデータ値の総和に対する平均値または中間値としてもよい。なお、これらは、基準画像の画像データを構成するデータ値の総和に関しても同様のことが言える。
【図面の簡単な説明】
【０１２４】
【図１】図１は、第１の実施形態による磁気共鳴イメージング装置の構成を示す図である。
【図２】図２は、データ処理部の構成を示すブロック図である。
【図３】図３は、第１の実施形態による磁気共鳴イメージング装置の動作の流れを示すフロー図である。
【図４】図４は、磁気共鳴イメージング装置により被検体をスキャンする際のシーケンスを示す図である。
【図５】図５は、被検体の撮像領域に設定された複数のスライスを示す図である。
【図６】図６は、第１の実施形態における画像再構成処理のフローを示す図である。
【図７】図７は、所定のスライスに対応する画像のローデータの概念図である。
【図８】図８はローデータに基づいて指標値を算出する処理の概念図である。
【図９】図９は、第１の実施形態による造影ダイナミックスキャンの各スキャンが実施された各時点における指標値および画像再構成要否判定結果の一例を示す図である。
【図１０】図１０は、第２の実施形態による造影ダイナミックスキャンの各スキャンが実施された各時点における指標値および画像再構成要否判定結果の一例を示す図である。
【図１１】図１１は、第３の実施形態における画像出力処理のフローを示す図である。
【図１２】図１２は、所定のスライスに対応する画像の画像データの概念図である。
【図１３】図１３は、画像データに基づいて指標値を算出する処理の概念図である。
【符号の説明】
【０１２５】
１  磁気共鳴イメージング装置
２  スキャン部
３  操作コンソール部
１１  撮像空間
１２  静磁場マグネット部
１３  勾配コイル部
１４  ＲＦコイル部
１５  撮影テーブル
２２  ＲＦ駆動部
２３  勾配駆動部
２４  データ収集部
３０  制御部
３１  記憶部
３２  操作部
３３  データ処理部
３４  表示部
４０  被検体
３３１  指標値算出部
３３２  しきい値判定部
３３３  画像再構成部
３３４  画像出力部
【図１】
2010104459.tif 000012
【図２】
2010104459.tif 000013
【図３】
2010104459.tif 000014
【図４】
2010104459.tif 000015
【図５】
2010104459.tif 000016
【図６】
2010104459.tif 000017
【図７】
2010104459.tif 000018
【図８】
2010104459.tif 000019
【図９】
2010104459.tif 000020
【図１０】
2010104459.tif 000021
【図１１】
2010104459.tif 000022
【図１２】
2010104459.tif 000023
【図１３】
2010104459.tif 000024