(57)【要約】  被写体の被曝線量を増やすことなく診断に要求される画質の動態画像を取得できるようにする。  本発明に係る動態画像撮影システムによれば、撮影用コンソールの制御部は、動態撮影の初期に得られた少なくとも一のフレーム画像を解析して画質を示す指標値を算出し、当該算出された画質を示す指標値に基づいて、動態撮影により得られるフレーム画像の画質を示す指標値が予め定められた基準値以下となる上限フレームレートｆｓｕｐを算出する。そして、この算出した上限フレームレートｆｓｕｐと撮影部位の動態を診断するのに必要な最低限の下限フレームレートｆｉｎｆとに基づいて、動態撮影に用いるフレームレートを決定し、撮影装置１に決定したフレームレートで撮影を行わせる。
2009139206.tif 000002
【特許請求の範囲】
【請求項１】
撮影部位に放射線を複数回連続照射することにより動態撮影を行って、前記撮影部位の動態を示す複数のフレーム画像を取得する撮影装置に接続された動態画像撮影制御装置であって、
撮影対象となる撮影部位を指定するための操作手段と、
前記操作手段により指定された撮影部位に応じて、動態撮影に用いる照射放射線量及びフレームレートを含む撮影条件を前記撮影装置に設定して動態撮影を行わせ、前記撮影部位の動態を示すフレーム画像を取得する撮影制御手段と、
動態撮影の本撮影前に前記撮影装置において取得された少なくとも一のフレーム画像を解析して画質を示す指標値を算出する指標値算出手段と、
前記算出された画質を示す指標値に基づいて、前記動態撮影の本撮影により取得されるフレーム画像の画質を示す指標値を予め定められた基準値以下にするフレームレートの上限を算出する上限フレームレート算出手段と、
を備え、
前記撮影制御手段は、前記算出された上限フレームレートに基づいて前記動態撮影の本撮影に用いるフレームレートを決定し、当該決定したフレームレートで前記撮影装置に動態撮影の本撮影を行わせる動態画像撮影制御装置。
【請求項２】
前記画質を示す指標値は、フレーム画像のノイズレベルであり、
前記指標値算出手段は、前記撮影装置により取得されたフレーム画像に関心領域を設定し、当該関心領域における非エッジ領域を抽出し、この抽出された非エッジ領域における高周波数帯域の画像信号成分の局所分散値を算出することにより、前記フレーム画像のノイズレベルを算出する請求項１に記載の動態画像撮影制御装置。
【請求項３】
撮影部位毎に、当該撮影部位の動態を診断するのに必要な下限フレームレートを記憶する記憶手段と、
前記記憶手段から前記操作手段により指定された撮影部位に応じた下限フレームレートを読み出し、前記算出された上限フレームレートと前記読み出された下限フレームレートとを比較し、比較の結果、前記上限フレームレートが前記下限フレームレートを下回る場合に、その旨を通知する通知手段と、
を備える請求項１又は２に記載の動態画像撮影制御装置。
【請求項４】
前記撮影部位の動態は、周期性をもつ動態であり、
前記撮影制御手段は、前記上限フレームレートが前記下限フレームレートを下回る場合には、前記撮影装置に前記下限フレームレートで動態撮影の本撮影を行わせ、
前記本撮影により取得された前記撮影部位の動態の複数サイクルのフレーム画像間で画素信号値の平均化を行うことによりノイズ低減処理を行うノイズ低減処理手段を備える請求項１〜３の何れか一項に記載の動態画像撮影制御装置。
【請求項５】
前記下限フレームレートで動態撮影を行うことにより得られるフレーム画像の画質を示す指標値を算出し、当該算出された画質を示す指標値に基づいて、前記ノイズ低減処理に用いる動態サイクル数を決定するサイクル数決定手段を備え、
前記ノイズ低減処理手段は、前記本撮影により取得された複数のフレーム画像のうち前記決定された動態サイクル数のフレーム画像間で画素信号値の平均化を行うことによりノイズ低減処理を行う請求項４に記載の動態画像撮影制御装置。
【請求項６】
撮影対象となる撮影部位を指定するための操作手段と、
前記指定された撮影部位に放射線を複数回連続照射して動態撮影を行い、前記撮影部位の動態を示す複数のフレーム画像を取得する撮影手段と、  前記操作手段により指定された撮影部位に応じて、動態撮影に用いる照射放射線量及びフレームレートを含む撮影条件を前記撮影手段に設定して動態撮影を行わせる撮影制御手段と、
動態撮影の本撮影前に前記撮影手段において取得された少なくとも一のフレーム画像を解析して画質を示す指標値を算出する指標値算出手段と、
前記算出された画質を示す指標値に基づいて、前記動態撮影の本撮影により取得されるフレーム画像の画質を示す指標値を予め定められた基準値以下にするフレームレートの上限を算出する上限フレームレート算出手段と、
を備え、
前記撮影制御手段は、前記算出された上限フレームレートに基づいて前記動態撮影の本撮影に用いるフレームレートを決定し、当該決定したフレームレートで前記撮影手段に動態撮影の本撮影を行わせる動態画像撮影システム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、動態画像撮影制御装置及び動態画像撮影システムに関する。
【背景技術】
【０００２】
医療の分野では、近年、従来のフィルム／スクリーンや輝尽性蛍光体プレートを用いた放射線の静止画撮影及び診断に対し、ＦＰＤ（flat panel detector）等の半導体イメージセンサを利用して被写体の動態画像を撮影し、診断に応用する試みがなされるようになってきている。具体的には、半導体イメージセンサの画像データの読取・消去の応答性の早さを利用し、半導体イメージセンサの読取・消去のタイミングと合わせて放射線源からパルス状の放射線を連続照射し、１秒間に複数回の撮影を行って被写体の動態を撮影する。撮影により取得された一連の複数枚の画像を順次表示することにより、医師は撮影部位の一連の動きを認識することが可能となる。【０００３】
動態撮影においては複数回の撮影を行うため、１回当たりの撮影において静止画撮影と同等の放射線を照射してしまうと被写体の被曝線量は相当量になってしまう。そこで、例えば、特許文献１には、動態撮影の総撮影線量を静止画撮影と同等若しくはそれ以下に抑えるように動態撮影を行う技術が記載されている。
【特許文献１】特開２００５−３１２７７５号公報
【発明の開示】
【発明が解決しようとする課題】
【０００４】
しかしながら、特許文献１においては、撮影するフレーム数を高く設定してしまうと画像１枚当たりの照射放射線量が低下してしまうことから各画像の画質低下が生じ、画像解析や診断に必要な画質が得られない場合がある。一方、動態撮影トータルでの照射放射線量を上げるようにすれば、画質低下は防ぐことはできるが、被写体の被曝線量が増加してしまう。【０００５】
本発明の課題は、被写体の被曝線量を増やすことなく診断に要求される画質の動態画像を取得できるようにすることである。
【課題を解決するための手段】
【０００６】
上記課題を解決するために、本発明の第１の側面によると、
撮影部位に放射線を複数回連続照射することにより動態撮影を行って、前記撮影部位の動態を示す複数のフレーム画像を取得する撮影装置に接続された動態画像撮影制御装置であって、
撮影対象となる撮影部位を指定するための操作手段と、
前記操作手段により指定された撮影部位に応じて、動態撮影に用いる照射放射線量及びフレームレートを含む撮影条件を前記撮影装置に設定して動態撮影を行わせ、前記撮影部位の動態を示すフレーム画像を取得する撮影制御手段と、
動態撮影の本撮影前に前記撮影装置において取得された少なくとも一のフレーム画像を解析して画質を示す指標値を算出する指標値算出手段と、
前記算出された画質を示す指標値に基づいて、前記動態撮影の本撮影により取得されるフレーム画像の画質を示す指標値を予め定められた基準値以下にするフレームレートの上限を算出する上限フレームレート算出手段と、
を備え、
前記撮影制御手段は、前記算出された上限フレームレートに基づいて前記動態撮影の本撮影に用いるフレームレートを決定し、当該決定したフレームレートで前記撮影装置に動態撮影の本撮影を行わせる。【０００７】
前記画質を示す指標値は、フレーム画像のノイズレベルであり、
前記指標値算出手段は、前記撮影装置により取得されたフレーム画像に関心領域を設定し、当該関心領域における非エッジ領域を抽出し、この抽出された非エッジ領域における高周波数帯域の画像信号成分の局所分散値を算出することにより、前記フレーム画像のノイズレベルを算出することが好ましい。【０００８】
前記動態画像撮影制御装置は、
撮影部位毎に、当該撮影部位の動態を診断するのに必要な下限フレームレートを記憶する記憶手段と、
前記記憶手段から前記操作手段により指定された撮影部位に応じた下限フレームレートを読み出し、前記算出された上限フレームレートと前記読み出された下限フレームレートとを比較し、比較の結果、前記上限フレームレートが前記下限フレームレートを下回る場合に、その旨を通知する通知手段と、
を備えることが好ましい。【０００９】
前記撮影部位の動態は、周期性をもつ動態であり、
前記撮影制御手段は、前記上限フレームレートが前記下限フレームレートを下回る場合には、前記撮影装置に前記下限フレームレートで動態撮影の本撮影を行わせ、
前記動態画像撮影制御装置は、更に、前記本撮影により取得された前記撮影部位の動態の複数サイクルのフレーム画像間で画素信号値の平均化を行うことによりノイズ低減処理を行うノイズ低減処理手段を備えることが好ましい。【００１０】
前記動態画像撮影制御装置は、更に、
前記下限フレームレートで動態撮影を行うことにより得られるフレーム画像の画質を示す指標値を算出し、当該算出された画質を示す指標値に基づいて、前記ノイズ低減処理に用いる動態サイクル数を決定するサイクル数決定手段を備え、
前記ノイズ低減処理手段は、前記本撮影により取得された複数のフレーム画像のうち前記決定された動態サイクル数のフレーム画像間で画素信号値の平均化を行うことによりノイズ低減処理を行うことが好ましい。【００１１】
本発明の第２の側面によると、動態画像撮影システムは、
撮影対象となる撮影部位を指定するための操作手段と、
前記指定された撮影部位に放射線を複数回連続照射して動態撮影を行い、前記撮影部位の動態を示す複数のフレーム画像を取得する撮影手段と、
前記操作手段により指定された撮影部位に応じて、動態撮影に用いる照射放射線量及びフレームレートを含む撮影条件を前記撮影手段に設定して動態撮影を行わせる撮影制御手段と、
動態撮影の本撮影前に前記撮影手段において取得された少なくとも一のフレーム画像を解析して画質を示す指標値を算出する指標値算出手段と、
前記算出された画質を示す指標値に基づいて、前記動態撮影の本撮影により取得されるフレーム画像の画質を示す指標値を予め定められた基準値以下にするフレームレートの上限を算出する上限フレームレート算出手段と、
を備え、
前記撮影制御手段は、前記算出された上限フレームレートに基づいて前記動態撮影の本撮影に用いるフレームレートを決定し、当該決定したフレームレートで前記撮影手段に動態撮影の本撮影を行わせる。 
【発明の効果】
【００１２】
本発明によれば、被写体の被曝線量を増やすことなく診断に要求される画質の動態画像を取得することが可能となる。
【図面の簡単な説明】
【００１３】
【図１】本発明の実施の形態における動態画像撮影システムの全体構成例を示す図である。
【図２】図１の撮影用コンソールの制御部により実行される撮影制御処理を示すフローチャートである。
【図３】図２のステップＳ５において実行されるノイズレベル算出処理を示すフローチャートである。
【図４】図３のステップＳ１０２において実行される多重解像度分解処理を説明するための図である。
【図５Ａ】多重解像度分解処理を施す前の関心領域における各画素（座標（ｉ）で画素位置を示す）の信号値を示す図である。
【図５Ｂ】図５Ａを多重解像度分解処理した結果得られる複数周波数帯域の画像信号成分のうち高周波数帯域の画像信号成分を示す図である。
【図５Ｃ】図５Ａを多重解像度分解処理した結果得られる複数周波数帯域の画像信号成分のうち低周波数帯域の画像信号成分を示す図である。
【図６】局所領域を説明するための図である。
【図７】規格化ノイズレベルとフレームレートの関係を示す図である。
【図８】規格化ノイズレベルと平均化に用いるサイクル数の関係を示す図である。
【図９】平均化に用いる動態サイクル数３である場合のサイクル間平均化によるノイズ低減処理を説明するための図である。
【図１０】濃度ヒストグラムのシフト量を説明するための図である。
【発明を実施するための最良の形態】
【００１４】
（動態画像撮影システム１００の構成）
まず、構成を説明する。
図１に、第１の実施の形態における動態画像撮影システム１００の全体構成を示す。
図１に示すように、動態画像撮影システム１００は、撮影装置１と、撮影用コンソール２とが通信ケーブル等により接続され、撮影用コンソール２と、診断用コンソール３とがＬＡＮ（Local Area Network）等の通信ネットワークＮを介して接続されて構成されている。【００１５】
（撮影装置１の構成）
撮影装置１は、例えば、呼吸に伴う肺の膨張及び収縮の形態変化、心臓の拍動等の、周期性（サイクル）をもつ人体の動態を撮影する撮影手段である。動態撮影は、被写体となる撮影部位に対し、放射線を連続照射して複数の画像を取得（即ち、連続撮影）することにより行う。この連続撮影により得られた一連の画像を動態画像と呼ぶ。また、動態画像を構成する複数の画像のそれぞれをフレーム画像と呼ぶ。
撮影装置１は、図１に示すように、放射線源１１、放射線照射制御装置１２、放射線検出部１３、読取制御装置１４、サイクル検出センサ１５、サイクル検出装置１６等を備えて構成されている。【００１６】
放射線源１１は、放射線照射制御装置１２の制御に従って、患者Ｐに対し放射線（Ｘ線）を照射する。
放射線照射制御装置１２は、撮影用コンソール２に接続されており、撮影用コンソール２から入力された放射線照射条件に基づいて放射線源１１を制御して放射線撮影を行う。撮影用コンソール２から入力される放射線照射条件は、例えば、連続照射時のパルスレート、パルス幅、パルス間隔、撮影開始／終了タイミング、撮影サイクル数、Ｘ線管電流の値、Ｘ線管電圧の値、フィルタ種等である。パルスレートは、単位時間あたり（ここでは１秒あたり）の放射線照射回数であり、後述するフレームレートと一致している。パルス幅は、放射線照射１回当たりの放射線照射時間である。パルス間隔は、連続撮影において、１回の放射線照射開始から次の放射線照射開始までの時間であり、後述するフレーム間隔と一致している。【００１７】
放射線検出部１３は、ＦＰＤ等の半導体イメージセンサにより構成される。ＦＰＤは、例えば、ガラス基板等を有しており、基板上の所定位置に、放射線源１１から照射されて少なくとも患者Ｐの撮影部位を透過した放射線をその強度に応じて検出し、検出した放射線を電気信号に変換して蓄積する複数の画素がマトリックス状に配列されている。各画素は、例えばＴＦＴ（Thin Film Transistor）等のスイッチング部により構成されている。
読取制御装置１４は、撮影用コンソール２に接続されており、撮影用コンソール２から入力された画像読取条件に基づいて放射線検出部１３の各画素のスイッチング部を制御して、当該各画素に蓄積された電気信号の読み取りをスイッチングしていき、放射線検出部１３に蓄積された電気信号を読み取ることにより、画像データを取得する。そして、取得した画像データを撮影用コンソール２に出力する。画像読取条件は、例えば、フレームレート、フレーム間隔、画素サイズ、画像サイズ（マトリックスサイズ）等である。フレームレートは、単位時間あたり（ここでは１秒あたり）に取得するフレーム画像数であり、パルスレートと一致している。フレーム間隔は、連続撮影において、１回のフレーム画像の取得動作開始から次のフレーム画像の取得動作開始までの時間であり、パルス間隔と一致している。
なお、パルス間隔及びフレーム間隔は、パルスレート、フレームレートから求められる。【００１８】
ここで、放射線照射制御装置１２と読取制御装置１４は互いに接続され、互いに同期信号をやりとりして放射線照射動作と画像の読み取りの動作を同調させるようになっている。【００１９】
サイクル検出センサ１５は、患者Ｐの撮影部位の状態を検出して検出情報をサイクル検出装置１６に出力する。サイクル検出センサ１５としては、例えば、撮影部位が肺（換気）の場合、呼吸モニタベルト、ＣＣＤ（Charge Coupled Device）カメラ、光学カメラ、スパイロメータ等を適用することができる。また、撮影部位が心臓（血流）である場合、心電計等を適用することができる。【００２０】
サイクル検出装置１６は、サイクル検出センサ１５により入力された検出情報に基づいて、撮影部位の動態のサイクル数（単位時間当たりのサイクル数。例えば、撮影部位が肺（換気）の場合は呼吸数（回／秒）、心臓の場合は心拍数（回／秒））、及び現在の撮影部位の状態が１サイクル中のどの状態であるか等を検出し、検出結果（サイクル情報）を撮影用コンソール２の制御部２１に出力する。
例えば、サイクル検出装置１６は、撮影部位が肺（換気）の場合は、サイクル検出センサ１５（呼吸モニタベルト、ＣＣＤカメラ、光学カメラ、スパイロメータ等）により肺の状態が吸気から呼気への変換点であることを示す検出情報が入力されたタイミングを１サイクルの基点とし、次にこの状態が検出されるタイミングまでの間を１サイクルとして認識する。
また、サイクル検出装置１６は、撮影部位が心臓（血流を含む）の場合は、サイクル検出センサ１５（心電計等）によりＲ波が入力されたタイミングを基点とし、次にＲ波が検出されるタイミングまでの間を１サイクルとして認識する。
そして、１秒当たりに認識したサイクル数をサイクル数として検出する。【００２１】
（撮影用コンソール２の構成）
撮影用コンソール２は、放射線照射条件や画像読取条件を撮影装置１に設定して撮影装置１による放射線撮影及び放射線画像の読み取り動作を制御するとともに、撮影装置１により取得された動態画像を撮影技師による確認のために表示する動態画像撮影制御装置である。
撮影用コンソール２は、図１に示すように、制御部２１、記憶部２２、操作部２３、表示部２４、通信部２５等を備えて構成され、各部はバス２６により接続されている。【００２２】
制御部２１は、ＣＰＵ（Central Processing Unit）、ＲＡＭ（Random Access Memory）等により構成される。制御部２１のＣＰＵは、操作部２３の操作に応じて、記憶部２２に記憶されているシステムプログラムや、各種処理プログラムを読み出してＲＡＭ内に展開し、展開されたプログラムに従って撮影用コンソール２各部の動作や、撮影装置１の放射線照射動作及び読み取り動作を集中制御する。
具体的に、制御部２１は、記憶部２２に記憶されている撮影制御処理プログラムを読み出して撮影制御処理を実行することにより、撮影制御手段、上限フレームレート算出手段、ノイズ低減処理手段、サイクル数決定手段としての機能を実現する。また、制御部２１は、記憶部２２に記憶されているノイズレベル算出処理プログラムを読み出してノイズレベル算出処理を実行することにより、指標値算出手段としての機能を実現する。【００２３】
記憶部２２は、不揮発性の半導体メモリやハードディスク等により構成される。記憶部２２は、制御部２１で実行される各種プログラムやプログラムにより処理の実行に必要なパラメータ、或いは処理結果等のデータを記憶する。各種プログラムは、読取可能なプログラムコードの形態で格納され、制御部２１は、当該プログラムコードに従った動作を逐次実行する。【００２４】
例えば、記憶部２２は、後述する撮影制御処理を実行するための撮影制御処理プログラムを記憶している。また、記憶部２２は、撮影部位毎に、その撮影部位に応じた撮影条件（放射線照射条件及び画像読取条件）を記憶している。記憶部２２に撮影部位毎に記憶されている撮影条件としては、例えば、フレームレート（パルスレート）、撮影開始／終了タイミング、撮影サイクル数、Ｘ線管電流の値、Ｘ線管電圧の値、総照射放射線量、画素サイズ、画像サイズ等が挙げられる。これらの撮影部位毎の撮影条件は、実験的経験的に得られた、各撮影部位の一般的な動態サイクル時間（１サイクル当たりに要する時間）、診断に必要とされる１サイクル当たりのフレーム画像数、許容される被爆線量等に基づき設定されたものである。【００２５】
また、記憶部２２は、撮影部位毎に、当該撮影部位に応じて予め定められた下限フレームレートｆｉｎｆを記憶している。下限フレームレートｆｉｎｆは、撮影部位の動態を診断するのに最低限必要として設定されているフレームレートであり、実験的経験的に得られた、撮影部位の単位時間当たりのサイクル数、診断に必要な１サイクルあたりのフレーム画像数等に基づき算出された値である。【００２６】
操作部２３は、カーソルキー、数字入力キー、及び各種機能キー等を備えたキーボードと、マウス等のポインティングデバイスを備えて構成され、キーボードに対するキー操作やマウス操作により入力された指示信号を制御部２１に出力する。また、操作部２３は、表示部２４の表示画面にタッチパネルを備えても良く、この場合、タッチパネルを介して入力された指示信号を制御部２１に出力する。【００２７】
表示部２４は、ＬＣＤ（Liquid Crystal Display）やＣＲＴ（Cathode Ray Tube）等のモニタにより構成され、制御部２１から入力される表示信号の指示に従って、操作部２３からの入力指示やデータ等を表示する。【００２８】
通信部２５は、ＬＡＮアダプタやルータやＴＡ（Terminal Adapter）等を備え、通信ネットワークＮに接続された各装置との間のデータ送受信を制御する。【００２９】
（診断用コンソール３の構成）
診断用コンソール３は、撮影用コンソール２から送信された動態画像を表示して医師が読影診断するための装置である。
診断用コンソール３は、図１に示すように、制御部３１、記憶部３２、操作部３３、表示部３４、通信部３５等を備えて構成され、各部はバス３６により接続されている。【００３０】
制御部３１は、ＣＰＵ、ＲＡＭ等により構成される。制御部３１のＣＰＵは、操作部３３の操作に応じて、記憶部３２に記憶されているシステムプログラムや、各種処理プログラムを読み出してＲＡＭ内に展開し、展開されたプログラムに従って診断用コンソール３各部の動作を集中制御する。【００３１】
記憶部３２は、不揮発性の半導体メモリやハードディスク等により構成される。記憶部３２は、制御部３１で実行される各種プログラムやプログラムにより処理の実行に必要なパラメータ、或いは処理結果等のデータを記憶する。これらの各種プログラムは、読取可能なプログラムコードの形態で格納され、制御部３１は、当該プログラムコードに従った動作を逐次実行する。【００３２】
操作部３３は、カーソルキー、数字入力キー、及び各種機能キー等を備えたキーボードと、マウス等のポインティングデバイスを備えて構成され、キーボードに対するキー操作やマウス操作により入力された指示信号を制御部３１に出力する。また、操作部３３は、表示部３４の表示画面にタッチパネルを備えても良く、この場合、タッチパネルを介して入力された指示信号を制御部３１に出力する。【００３３】
表示部３４は、ＬＣＤやＣＲＴ等のモニタにより構成され、制御部３１から入力される表示信号の指示に従って、操作部３３からの入力指示やデータ等を表示する。【００３４】
通信部３５は、ＬＡＮアダプタやルータやＴＡ等を備え、通信ネットワークＮに接続された各装置との間のデータ送受信を制御する。【００３５】
（動態画像撮影システム１００の動作）
次に、上記動態画像撮影システム１００における動作について説明する。
図２に、撮影用コンソール２の制御部２１により実行される撮影制御処理の流れを示す。【００３６】
まず、撮影技師により操作部２３による患者Ｐの患者情報（氏名、身長、体重、年齢、性別等）の入力、撮影部位の指定入力等が受け付けられる（ステップＳ１）。【００３７】
次いで、操作部２３により指定された撮影部位に応じた放射線照射条件が放射線照射制御装置１２に設定されるとともに、指定された撮影部位に応じた画像読取条件が読取制御装置１４に設定される（ステップＳ２）。ステップＳ２においては、まず、記憶部２２から撮影条件（放射線照射条件及び画像読取条件）が読み出される。次いで、読み出された撮影条件と、サイクル検出装置１６から取得された撮影部位のサイクル情報に基づいて、更に必要な撮影条件が算出される。例えば、記憶部２２から読み出されたフレームレート（パルスレート）、総照射放射線量、及びサイクル検出装置１６から検出された撮影部位のサイクル情報に基づいて、動態撮影トータルの照射放射線量が予め設定されている総照射放射線量以下となるように、他の撮影条件、例えば、パルス幅等が算出される。そして、読み出された撮影条件及び算出された撮影条件のうち、放射線照射条件が放射線照射制御装置１２に設定され、画像読取条件が読取制御装置１４に設定される。【００３８】
次いで、操作部２３の操作による撮影開始の指示の入力が待機され、操作部２３により撮影開始指示が入力されると（ステップＳ３；ＹＥＳ）、放射線照射制御装置１２及び読取制御装置１４に撮影指示が出力され、設定された撮影条件（フレームレートを含む）に基づいて１又は数フレーム分のプレ撮影が行われる（ステップＳ４）。プレ撮影は、動態撮影の本撮影前に行われる撮影である。本実施の形態では、本撮影時のフレームレートの調整を行うための画像解析用のフレーム画像を取得するためにプレ撮影が実施される。そして、撮影により取得されたフレーム画像に対しノイズレベル算出処理が実行される（ステップＳ５）。【００３９】
図３に、図２のステップＳ５において制御部２１により実行されるノイズレベル算出処理の流れを示す。
まず、撮影により取得された何れか１のフレーム画像（フレーム画像Ｆとする）に対して関心領域が設定される（ステップＳ１０１）。関心領域は、撮影部位に応じて予め定められており、撮影部位に応じた位置に関心領域が設定される。なお、フレーム画像Ｆを画像解析して関心領域を抽出するようにしてもよい。【００４０】
次いで、関心領域の画像に対して多重解像度分解処理が行われる（ステップＳ１０２）。
ここで、図４を参照して多重解像度分解処理について説明する。なお、本実施の形態においては、ラプラシアンピラミッド法により入力画像（ここでは、関心領域の画像）を多重解像分解して互いに異なる複数の周波数帯域の画像信号成分に分解する例について説明する。このラプラシアンピラミッド法によれば迅速な処理が可能となるが、ウェーブレット変換等、他の手法を用いることも可能である。【００４１】
まず、関心領域の画像信号（原画像信号Ｓとする）が、ローパスフィルタを構成するフィルタリング手段１２１ａによりフィルタリングされる。【００４２】
フィルタリング手段１２１ａによりフィルタリングされた原画像信号Ｓは１画素おきにサンプリング（ダウンサンプリング）されることで、低解像度近似画像ｇ1が生成される。この低解像度近似画像ｇ1は、原画像の１／４の大きさになっている。【００４３】
次いで、補間手段１２２ａにおいて、この低解像度近似画像ｇ１のサンプリングされた間隔に値が０の画素が補間される。この補間は、低解像度近似画像ｇ１の列毎および１行ごとに値が０の行および列を挿入することにより行なう。なお、このように補間された低解像度近似画像は、１画素おきに値が０の画素が挿入されているため、信号値の変化が滑らかではない状態になっている。そして、このような補間が行われた後に、補間手段１２２ａに含まれるローパスフィルタにおいて、再度フィルタリングを施し、低解像度近似画像ｇ１′を得る。この低解像度近似画像ｇ１′は、前記した補間直後の低解像度近似画像に比べると信号値の変化が滑らかな状態になっている。【００４４】
この低解像度近似画像ｇ１′は、画像を１／４にした後に１画素おきに０の補間とフィルタリングとをすることにより、原画像の空間周波数の半分より高い周波数が消えた状態になっている。【００４５】
そして、減算器１２３ａにより、原画像から低解像度近似画像ｇ１′の減算を行って、細部画像ｂ０を得る。この減算は、原画像と低解像度近似画像ｇ１′との対応する画素の間で実行される。これにより、細部画像ｂ０は、原画像の空間周波数の半分より高い周波数帯域のみを示す画像になる。すなわち、細部画像ｂ０は、ナイキスト周波数をＮｙとすると、Ｎｙ／２〜Ｎｙの周波数帯域の画像信号成分である。【００４６】
さらに、前述したローパスフィルタ１２１ａからの低解像度近似画像ｇ1が、ローパスフィルタ１２１ｂによってフィルタリング処理される。これにより低解像度近似画像ｇ１は、さらに１画素おきにサンプリングされて１／４（もとの１／１６）の低解像度近似画像ｇ２に変換される。そして、この低解像度近似画像ｇ２にも、補間手段１２２ｂと減算器１２３ｂとによって同様な処理が施されて、低解像度近似画像ｇ２′から細部画像ｂ１が生成される。すなわち、細部画像ｂ１は、ナイキスト周波数をＮｙとすると、Ｎｙ／４〜Ｎｙ／２の周波数帯域の画像信号成分である。【００４７】
このような処理を順次繰り返すことで、ローパスフィルタ１２１により生成された低解像度近似画像ｇｋ （ここで、ｋ＝１〜Ｌ）から、細部画像ｂｋ-１、および、残留画像ｇＬを得る。【００４８】
ここで、細部画像ｂｋ-１は、ｂ０から順に解像度が低くなる、すなわち、周波数帯域が低くなるものであり、細部画像ｂｋ-１は（Ｎｙ／２＾(（ｋ-１）+１)）〜（Ｎｙ／２＾（ｋ-１））の周波数帯域の画像信号成分になっている。
このように、多重解像度分解処理により、関心領域の画像信号が複数の周波数帯域の画像信号成分に分解される。【００４９】
図５Ａ〜図５Ｃに、上述した多重解像度分解処理を施す前の関心領域における各画素の信号値及びこれに多重解像度分解処理を施して得られる複数周波数帯域の画像信号成分の信号値を模式的に示す。
図５Ａは、多重解像度分解処理を施す前の関心領域における各画素（座標（ｉ）で画素位置を示す）の信号値を示す図である。図５Ｂは、図５Ａを多重解像度分解処理した結果得られる複数周波数帯域の画像信号成分のうち高周波数帯域の画像信号成分（例えば、細部画像ｂ０）を示す図である。図５Ｃは、図５Ａを多重解像度分解処理した結果得られる複数周波数帯域の画像信号成分のうち低周波数帯域の画像信号成分（例えば、細部画像ｂＬ−１）を示す図である。
なお、説明を簡略化するために、以降ステップＳ１０１で設定された関心領域が図５Ａに示す領域Ａのように水平方向に延びる１次元の領域であることとして説明するが、２次元の領域であってもよい。  【００５０】
多重解像度分解処理が終了すると、この処理結果を用いて、ステップＳ１０３、Ｓ１０４の処理によりフレーム画像ＦのノイズレベルσＮが算出される。
ここで、構造物の輪郭の存在しない画像領域（非エッジ領域）においては、低周波帯域の画像信号成分の分散値及び高周波数帯域の画像信号成分の分散値はともに低い値を示す筈である。しかし、その領域にノイズがある場合、高周波数帯域の画像信号成分の分散値は高くなる。そこで、本実施の形態においては、関心領域内の非エッジ領域における高周波数帯域の画像信号成分の平均分散値をフレーム画像ＦのノイズレベルσＮとして算出する。ノイズレベルσＮは画質を示す指標値である。ノイズレベルσＮが高いほど画質は悪いことを示す。
以下、ノイズレベルσＮの算出について説明する。【００５１】
まず、低周波数帯域の画像信号成分を用いて非エッジ領域が抽出される（ステップＳ１０３）。
非エッジ領域の抽出は、以下の（１）〜（２）により行うことができる。
（１）関心領域がｎ画素毎（ｎは予め定められた正の数）の局所領域に分割され、局所領域毎に分散値σが算出される。
例えば、図６に示す座標ｉ＝ｘ＋１〜ｘ＋ｎの局所領域での分散値（局所分散値）σは、〔数１〕に示す式で求めることができる。
【数１】
2009139206.tif 000003
ただし、Ａｖｅとは、当該局所領域の各画素の信号値の平均値を指す。【００５２】
（２）次いで、各局所領域における分散値σと予め設定された閾値σｍｉｄとが比較され、σ＜σｍｉｄとなる局所領域が非エッジ領域として抽出される。
例えば、図５Ｃでは、領域１、領域３が非エッジ領域として抽出される。なお、σｍｉｄは、例えば、下記の〔数２〕で算出される値を用いることができるが、これに限定されるものではない。
【数２】
2009139206.tif 000004
ただし、〔数２〕においてσは低周波数帯域における関心領域全体の分散値、ｎは画素数。【００５３】
非エッジ領域が抽出されると、この抽出された非エッジ領域に対応する高周波数帯域の領域（例えば、図５Ｂに示す領域１及び領域３）における画像信号成分の平均分散値がノイズレベルσＮとして算出される（ステップＳ１０４）。
具体的には、まず、抽出された各非エッジ領域毎の分散値σｓｕｐｔ（ｔ＝１〜Ｎ（Ｎは非エッジ領域の数）が算出される。分散値σｓｕｐｔは〔数１〕で示したものと同様の式で算出できる。次いで、Ｎ個の非エッジ領域の平均分散値（ノイズレベルσＮ）が算出される。ノイズレベルσＮを算出する式は、〔数３〕に示すとおりである。
【数３】
2009139206.tif 000005
ノイズレベルσＮが算出されると、処理は図２のステップＳ６に移行する。  【００５４】
図２のステップＳ６においては、算出されたノイズレベルσＮと、図７及び〔数４〕に示す規格化ノイズレベルとフレームレートの関係式とに基づいて、動態撮影により得られるフレーム画像のノイズレベルが予め定められた基準値σｌｉｍｉｔ以下となる上限フレームレートｆｓｕｐが算出され、表示部２４に表示される等により撮影技師に通知される（ステップＳ６）。
【数４】
2009139206.tif 000006
【００５５】
ここで、１回の動態撮影の総照射放射線量は撮影部位毎に予め定められており、制御部２１は、この総照射放射線量に基づいて、１フレーム画像当たりの照射放射線量を制御している。即ち、フレームレートを高くすると１フレーム画像当たりの照射放射線量が少なくなる。被写体としての撮影部位の動態撮影全体での被曝量を抑えるためである。従って、フレームレートを高くするほど１フレーム画像毎のノイズレベルは高くなる。予め定められた基準値σｌｉｍｉｔは、診断に許容されるノイズレベルとして実験的経験的に得られた値である。ステップＳ６においては、ノイズレベルがこのσｌｉｍｉｔ以下となる上限のフレームレートｆｓｕｐが算出される。σｌｉｍｉｔは、例えば上述の〔数２〕に示す式で算出される値を用いることができるが、これに限定されるものではない。【００５６】
図７及び〔数４〕に示す式より、フレーム画像Ｆの撮影に用いられたフレームレートがｆ１のときのノイズレベルがσＮであった場合、ｆｓｕｐは〔数５〕に示す式により求めることができる。  
【数５】
2009139206.tif 000007
【００５７】
次いで、撮影部位に応じて予め定められた下限フレームレートｆｉｎｆが記憶部２２から読み出され、下限フレームレートｆｉｎｆが取得される（ステップＳ７）。下限フレームレートｆｉｎｆは、上述のように、撮影部位の動態を診断するのに最低限必要として設定されているフレームレートである。例えば、撮影部位が心臓で、診断には心拍サイクルを１０％以上の時間分解能で撮影した動態画像が必要となるとすると、心拍数が６０回／分のときｆｉｎｆ＞１０フレーム／秒である。
なお、本実施の形態においては、下限フレームレートｆｉｎｆが予め記憶部２２に記憶されていることとするが、サイクル検出装置１６から取得した撮影部位の動態１サイクルに要する時間と、診断に必要な単位時間当たりのフレーム画像数とに基づき制御部２１において下限フレームレートｆｉｎｆを算出するようにしてもよい。【００５８】
次いで、算出された上限フレームレートｆｓｕｐと下限フレームレートｆｉｎｆが比較され、ｆｓｕｐ≧ｆｉｎｆである場合（ステップＳ８；ＹＥＳ）、フレームレートがｆｓｕｐに決定されて放射線照射制御装置１２に設定され、フレームレートｆｓｕｐで動態撮影の本撮影が行われる（ステップＳ９）。ステップＳ９においては、フレームレートｆｓｕｐに基づいて、本撮影におけるトータルの照射放射線量が予め設定されている総照射放射線量となるように、他の撮影条件、例えば、パルス幅等が併せて変更される。そして、変更された撮影条件が放射線照射制御装置１２及び読取制御装置１４に設定され、変更された撮影条件で本撮影が行われるよう制御が行われる。なお、本実施の形態においては、以降の撮影におけるフレームレートをフレームレートｆｓｕｐに決定とすることとしたが、ｆｓｕｐ以下でｆｉｎｆ以上の何れのフレームレートに決定してもよい。フレームレートをｆｓｕｐとすることで、ノイズレベルが許容される範囲内で撮影部位の動きを最も精細に再現する動態画像を取得することができる。【００５９】
撮影により取得されたフレーム画像は順次撮影用コンソール２に入力される。サイクル検出装置１６により、フレームレートｆｓｕｐで撮影を開始してから予め定められた撮影サイクル数の動態サイクルが検出されると、放射線照射制御装置１２及び読取制御装置１４に撮影終了の指示が出力され、撮影動作が停止される。撮影終了後、処理はステップＳ１５に移行する。【００６０】
一方、ｆｓｕｐ＜ｆｉｎｆである場合（ステップＳ８；ＮＯ）、通知手段としての表示部２４にエラーメッセージが表示され、ｆｓｕｐ＜ｆｉｎｆであること等がユーザである撮影技師に通知される（ステップＳ１０）。
ここで、撮影に用いるフレームレートが診断に必要な下限フレームレートｆｉｎｆに満たない場合、撮影で得られる動態画像では撮影部位の動きを十分に表しておらず診断に使用できない。そのため、ｆｓｕｐ＜ｆｉｎｆの場合、撮影を続行するとすればｆｉｎｆで本撮影をすることとなる。しかし、ｆｉｎｆで撮影すると、本撮影により得られる動態画像のノイズレベルが許容されるノイズレベルσｌｉｍｉｔより高くなってしまう。そこで、表示部２４には、ｆｓｕｐ＜ｆｉｎｆであること、撮影を続行するのであればフレームレートｆｉｎｆで本撮影を行うこと、本撮影により取得される画像はノイズレベルが許容されるレベルに満たないのでノイズ低減処理を行うこと等が表示される。また、撮影を続行する指示をユーザが入力するための「撮影続行」ボタン、撮影を中止する指示をユーザが入力するための「撮影中止」ボタンが表示される。ユーザは、操作部２３から何れかのボタンを操作する。【００６１】
操作部２３により「撮影中止」ボタンが押下されると（ステップＳ１１；ＮＯ）、本処理は終了する。【００６２】
一方、操作部２３により「撮影続行」ボタンが押下されると（ステップＳ１１；ＹＥＳ）、
サイクル間平均化によるノイズ低減処理に用いる動態画像の動態サイクル数Ｍが算出される（ステップＳ１２）。【００６３】
ここで、サイクル間平均化によるノイズ低減処理について説明する。
周期性のある動態の複数サイクルをフレームレートを一定にして撮影すると、各サイクルにおいてｎ番目に撮影される撮影部位の状態（サイズ、位置等）は、略同じ状態となる。また、Ｍサイクル分のフレーム画像における互いに対応する位置にある画素の信号値を平均化すると、図８に示すようにノイズレベルは略１／√Ｍに低減されることが知られている。そこで、本実施の形態においては、下限フレームレートｆｉｎｆで撮影した場合のノイズレベルがσｌｉｍｉｔとなる動態サイクル数Ｍを算出し、Ｍサイクルの動態撮影を行う。そして、Ｍサイクルの各サイクルにおいてｎ番目に撮影されたフレーム画像同士の互いに対応する画素の信号値を平均化した画像を作成することで、ノイズ低減処理を行う。【００６４】
具体的に、ステップＳ１２においては、図２のステップＳ５においてフレーム画像Ｆから算出されたノイズレベルと、〔数４〕及び図７に示す規格化ノイズレベルとフレームレートとの関係式に基づいて、フレームレートｆｉｎｆにおけるノイズレベルσＮが算出される。そして、算出されたノイズレベルσＮ、図８に示す規格化されたノイズレベルと平均化に用いるサイクル数との関係式に基づいて、ノイズレベルが上記σｌｉｍｉｔ以下になる動態サイクル数Ｍが算出される。動態サイクル数Ｍは、以下に示す〔数６〕により算出することができる。
【数６】
2009139206.tif 000008
【００６５】
また、動態撮影のフレームレートがｆｉｎｆに決定されて放射線照射制御装置１２に設定され、フレームレートｆｉｎｆで本撮影が開始される（ステップＳ１３）。なお、ステップＳ１３においては、フレームレートｆｉｎｆに基づいて、本撮影におけるトータルの照射放射線量が予め設定されている総照射放射線量となるように、他の撮影条件、例えば、パルス幅等が併せて変更される。そして、変更された撮影条件が放射線照射制御装置１２及び読取制御装置１４に設定され、変更された撮影条件での撮影が行われるよう制御が行われる。【００６６】
撮影により取得された各フレーム画像の画像データは順次撮影用コンソール２に入力される。  
サイクル検出装置１６により、フレームレートｆｉｎｆで撮影を開始してから予め定められた撮影サイクル数の動態サイクルが検出されると算出された撮影すべき撮影サイクル数の動態が検出されると、放射線照射制御装置１２及び読取制御装置１４に撮影終了の指示が出力され、撮影動作が停止される。【００６７】
撮影が終了すると、サイクル間平均化によるノイズ低減処理が実行される（ステップＳ１４）。サイクル間平均化によるノイズ低減処理では、Ｍサイクルの各サイクルにおいてｎ番目に撮影されたフレーム画像間の互いに対応する画素の信号値を平均化する。図９に、動態サイクル数Ｍが３である場合のノイズ低減処理の例を示す。図９に示すように、例えばｉサイクル目のｊフレーム目のフレーム画像Ｉｏｕｔは、ｉ−１サイクル、ｉサイクル、ｉ＋１サイクルのｊフレーム目の対応する画素毎に信号値を加算して撮影サイクル数３で除算することにより取得される。なお、図９においては、単純平均の場合について説明しているが、平均化するサイクル間に重みを付けて加重平均を行っても良い。ノイズ低減処理が終了すると、処理はステップＳ１５に移行する。【００６８】
ステップＳ１５においては、本撮影により取得されたフレーム画像に濃度補正等の画像処理が施され、撮影順を示す番号と対応付けて記憶部２２に保存されるとともに、表示部２４に表示される（ステップＳ１５）。撮影技師は、表示された動態画像によりポジショニングや画質等を確認し、本撮影により診断に適した画像が取得された（撮影ＯＫ）か、再撮影が必要（撮影ＮＧ）か、を判断する。そして、操作部２３を操作して、判断結果を入力する。【００６９】
操作部２３の所定の操作により撮影ＯＫを示す判断結果が入力されると（ステップＳ１６；ＹＥＳ）、動態撮影の本撮影で取得された一連のフレーム画像のそれぞれに、動態画像を識別するための識別ＩＤや、患者情報、撮影部位、放射線照射条件、画像読取条件、撮影順を示す番号、サイクル情報、ノイズ低減処理の有無等の情報が付帯され（例えば、ＤＩＣＯＭ形式で画像データのヘッダ領域に書き込まれ）、通信部２５を介して診断用コンソール３に送信される（ステップＳ１７）。そして、本処理は終了する。一方、操作部２３の所定の操作により撮影ＮＧを示す判断結果が入力されると（ステップＳ１６；ＮＯ）、記憶部２２に記憶された一連のフレーム画像が削除され（ステップＳ１８）、本処理は終了する。【００７０】
診断用コンソール３においては、撮影用コンソール２から動態画像の一連のフレーム画像が通信部３５により受信されると、制御部３１により、記憶部３２に受信されたフレーム画像が記憶される。【００７１】
診断用コンソール３において、操作部３３により識別ＩＤ等が入力され、診断対象の動態画像が指定されると、制御部３１により、記憶部３２から指定された動態画像の一連のフレーム画像が読み出される。そして、表示部３４に順次表示されることにより、医師の読影診断に供される。【００７２】
以上説明したように、動態画像撮影システム１００によれば、撮影用コンソール２の制御部２１は、プレ撮影により得られた少なくとも一のフレーム画像を解析して画質を示す指標値を算出し、当該算出された画質を示す指標値に基づいて、動態撮影の本撮影により得られるフレーム画像の画質を示す指標値を予め定められた基準値以下にする上限フレームレートｆｓｕｐを算出する。そして、この算出した上限フレームレートｆｓｕｐと撮影部位の動態を診断するのに必要な最低限の下限フレームレートｆｉｎｆとに基づいて、動態撮影に用いるフレームレートを決定し、撮影装置１に決定したフレームレートで撮影を行わせる。【００７３】
従って、被写体の被曝線量を増やすことなく診断に要求される画質の動態画像を取得することが可能となる。【００７４】
制御部２１が算出する画質を示す指標値はノイズレベルであり、制御部２１は、ノイズレベルを予め定められた基準値以下にする上限フレームレートｆｓｕｐを算出し、算出した上限フレームレートｆｓｕｐに基づいてフレームレートを決定して撮影装置１に動態撮影を行わせるので、ノイズレベルが診断に許容されるレベル以下の画質の動態画像を取得することが可能となる。【００７５】
また、制御部２１は、上限フレームレートｆｓｕｐと撮影部位の動態を診断するのに必要な最低限の下限フレームレートｆｉｎｆとを比較し、ｆｓｕｐ＜ｆｉｎｆである場合にはその旨を表示部２４に表示させるので、動態撮影により得られる動態画像がそのままでは診断に適したものではない場合にその旨を撮影技師が把握することが可能となる。【００７６】
また、ｆｓｕｐ＜ｆｉｎｆであっても撮影を続行する指示が操作部２３により入力された場合には、下限フレームレートｆｉｎｆで本撮影を行った場合にフレーム画像をサイクル間で平均化することによりノイズレベルが予め定められた基準値以下となる動態サイクル数Ｍを決定し、下限フレームレートｆｉｎｆで本撮影されたサイクル数Ｍのフレーム画像間で画素信号値の平均化を行うことによりノイズ低減処理を行う。従って、ノイズ低減処理により診断に許容されるノイズレベル以下の画質の動態画像を取得することが可能となる。【００７７】
なお、上述した本実施の形態における記述は、本発明に係る好適な動態画像撮影システムの一例であり、これに限定されるものではない。
例えば、上記実施の形態においては、画質を示す指標値として、分散値から算出したノイズレベルを使用することとして説明したが、これに限定されず、例えば、エントロピー値を用いても良い。【００７８】
エントロピー値について説明すると、ある領域において、中周波数帯域でエントロピー値が所定の値より小さく、最高周波数帯域でエントロピー値が大きい場合、当該領域ではノイズ成分が支配的であるといえる。
図６に示すような座標ｉ＝ｘ＋１〜ｘ＋ｎの局所領域でのエントロピー値Ｓｅは、下記の〔数７〕の式により求めることができる。
【数７】
2009139206.tif 000009
ただし、ｍとは階調数を、Ｐ（ｚ）とは、当該領域内の画素が画素値（信号値）“ｚ”をとる確率を表す。【００７９】
エントロピー値を用いた場合の撮影制御処理は、図２に示す処理と略同様である。まず、撮影された１のフレーム画像を多重解像度分解する。次いで、中周波数帯域でエントロピー値が所定の値より小さい領域を抽出し、最高周波数帯域において、上記抽出された領域のエントロピー値（平均）を算出する。そして、算出されたエントロピー値、予め定められている、規格化されたエントロピー値とフレームレートの関係式に基づいて、エントロピー値を所定値以下とする上限フレームレートをｆｓｕｐとして決定する。平均化によるノイズ低減処理に用いる動態サイクル数Ｍを算出する際には、予め定められている、規格化されたエントロピー値と平均化に使用する動態サイクル数Ｍとの関係式に基づいて、撮影サイクル数を決定する。【００８０】
また、撮影により取得された画像は、濃度ヒストグラムが作成され、所定の濃度（例えば、最高濃度等）が予め定められた濃度値となるように濃度ヒストグラムをシフトさせる濃度補正が行われる。そこで、初期に撮影した１フレームの画像から濃度ヒストグラムのシフト量ｓを求め（図１０参照）、このシフト量ｓを画質を示す指標値として用いても良い。シフト量ｓが大きいほど基準の濃度値とかけ離れていることから、シフト量ｓは大きいほど画質が悪いことを示す。
シフト量ｓを用いた場合の撮影制御処理では、まず、撮影された１のフレーム画像から濃度ヒストグラムを作成し、シフト量ｓを算出する。次いで、算出されたシフト量ｓ、予め定められている、規格化された濃度ヒストグラムのシフト量とフレームレートの関係式に基づいて、シフト量が所定値以下となる上限フレームレートをｆｓｕｐとして決定する。平均化によるノイズ低減処理に用いる動態サイクル数Ｍを算出する際には、予め定められている、規格化された濃度ヒストグラムのシフト量と平均化に使用する撮影サイクル数との関係式に基づいて、動態サイクル数Ｍを決定する。【００８１】
また、上記の説明では、本発明に係るプログラムのコンピュータ読み取り可能な媒体としてハードディスクや半導体の不揮発性メモリ等を使用した例を開示したが、この例に限定されない。その他のコンピュータ読み取り可能な媒体として、ＣＤ-ＲＯＭ等の可搬型記録媒体を適用することが可能である。また、本発明に係るプログラムのデータを通信回線を介して提供する媒体として、キャリアウエーブ(搬送波)も適用される。【００８２】
その他、動態画像撮影システム１００を構成する各装置の細部構成及び細部動作に関しても、本発明の趣旨を逸脱することのない範囲で適宜変更可能である。【００８３】
なお、明細書、請求の範囲、図面及び要約を含む２００８年５月１４日に出願された日本特許出願Ｎｏ．２００８−１２６６３４号の全ての開示は、そのまま本出願の一部に組み込まれる。
【産業上の利用可能性】
【００８４】
医療の分野において、患者の動態画像を撮影する動態画像撮影システムに利用可能性がある。
【符号の説明】
【００８５】
１００  動態画像撮影システム
１  撮影装置
１１  放射線源
１２  放射線照射制御装置
１３  放射線検出部
１４  読取制御装置
１５  サイクル検出センサ
１６  サイクル検出装置
２  撮影用コンソール
２１  制御部
２２  記憶部
２３  操作部
２４  表示部
２５  通信部
２６  バス
３  診断用コンソール
３１  制御部
３２  記憶部
３３  操作部
３４  表示部
３５  通信部
３６  バス
【図１】
2009139206.tif 000010
【図２】
2009139206.tif 000011
【図３】
2009139206.tif 000012
【図４】
2009139206.tif 000013
【図５Ａ】
2009139206.tif 000014
【図５Ｂ】
2009139206.tif 000015
【図５Ｃ】
2009139206.tif 000016
【図６】
2009139206.tif 000017
【図７】
2009139206.tif 000018
【図８】
2009139206.tif 000019
【図９】
2009139206.tif 000020
【図１０】
2009139206.tif 000021
【国際調査報告】
2009139206.tif 000022
2009139206.tif 000023
2009139206.tif 000024
2009139206.tif 000025