(57)【特許請求の範囲】
【請求項１】
人体または動物の所定部位を時間順次に捉えた動態画像を取得する動態画像取得部と、
前記所定部位の物理的状態の時間変化を検出する検出部と、
前記検出部によって検出された前記所定部位の物理的状態の時間変化に基づいて解析を行い、その解析結果を診断支援情報として生成する診断支援情報生成部と、
前記診断支援情報を前記動態画像と時間的に関連付けて保持する保持部と、
前記動態画像及び前記診断支援情報を表示するための表示用画像を生成する表示画像生成部と、
を備え、
前記診断支援情報は、前記解析に基づく第１の解析結果と第２の解析結果とを含み、
前記表示用画像は、
前記動態画像を表示する動態画像表示部分と、
前記診断支援情報の第１の解析結果と第２の解析結果とを識別可能に時間軸方向に一覧できるように表示する一覧表示部分と、
を含む画像であり、
前記検出部は、
前記所定部位の物理的状態の周期的な時間変化となる所定部位周期を特定する所定部位周期特定部、
を備え、
前記所定部位は肺野であり、
前記第１の解析結果は呼気を示し、
前記第２の解析結果は吸気を示す、ことを特徴とする、
画像生成装置。
【請求項２】
人体または動物の所定部位を時間順次に捉えた動態画像を取得する動態画像取得部と、
前記所定部位の物理的状態の時間変化を検出する検出部と、
前記検出部によって検出された前記所定部位の物理的状態の時間変化に基づいて解析を行い、その解析結果を診断支援情報として生成する診断支援情報生成部と、
前記診断支援情報を前記動態画像と時間的に関連付けて保持する保持部と、
前記動態画像及び前記診断支援情報を表示するための表示用画像を生成する表示画像生成部と、
を備え、
前記診断支援情報は、前記解析に基づく第１の解析結果と第２の解析結果とを含み、
前記表示用画像は、
前記動態画像を表示する動態画像表示部分と、
前記診断支援情報の第１の解析結果と第２の解析結果とを識別可能に時間軸方向に一覧できるように表示する一覧表示部分と、
を含む画像であり、
前記所定部位は複数の所定部位を含み、
前記検出部は、複数の所定部位の物理的状態の時間変化をそれぞれ検出し、
前記診断支援情報生成部は、複数の所定部位の物理的状態の時間変化に基づいてそれぞれ解析を行い、複数の所定部位についての解析結果を複数の診断支援情報としてそれぞれ生成し、
前記一覧表示部分は、前記複数の診断支援情報を表示することを特徴とする、
画像生成装置。
【請求項３】
人体または動物の所定部位を時間順次に捉えた動態画像を取得する動態画像取得部と、
前記所定部位の物理的状態の時間変化を検出する検出部と、
前記検出部によって検出された前記所定部位の物理的状態の時間変化に基づいて解析を行い、その解析結果を診断支援情報として生成する診断支援情報生成部と、
前記診断支援情報を前記動態画像と時間的に関連付けて保持する保持部と、
前記動態画像及び前記診断支援情報を表示するための表示用画像を生成する表示画像生成部と、
を備え、
前記診断支援情報は、前記解析に基づく第１の解析結果と第２の解析結果とを含み、
前記表示用画像は、
前記動態画像を表示する動態画像表示部分と、
前記診断支援情報の第１の解析結果と第２の解析結果とを識別可能に時間軸方向に一覧できるように表示する一覧表示部分と、
を含む画像であり、
前記検出部は、
前記所定部位の物理的状態の時間変化において、特徴点を算出する特徴点算出部、
を備えることを特徴とする、
画像生成装置。
【請求項４】
人体または動物の所定部位を時間順次に捉えた動態画像を取得する動態画像取得部と、
前記所定部位の物理的状態の時間変化を検出する検出部と、
前記検出部によって検出された前記所定部位の物理的状態の時間変化に基づいて解析を行い、その解析結果を診断支援情報として生成する診断支援情報生成部と、
前記診断支援情報を前記動態画像と時間的に関連付けて保持する保持部と、
前記動態画像及び前記診断支援情報を表示するための表示用画像を生成する表示画像生成部と、
を備え、
前記診断支援情報は、前記解析に基づく第１の解析結果と第２の解析結果とを含み、
前記表示用画像は、
前記動態画像を表示する動態画像表示部分と、
前記診断支援情報の第１の解析結果と第２の解析結果とを識別可能に時間軸方向に一覧できるように表示する一覧表示部分と、
を含む画像であり、
前記検出部は、
前記動態画像に基づき、前記所定部位の物理的状態の時間変化を検出することを特徴とする、
画像生成装置。
【請求項５】
人体または動物の所定部位を時間順次に捉えた動態画像を取得する動態画像取得部と、
前記所定部位の物理的状態の時間変化を検出する検出部と、
前記検出部によって検出された前記所定部位の物理的状態の時間変化に基づいて解析を行い、その解析結果を診断支援情報として生成する診断支援情報生成部と、
前記診断支援情報を前記動態画像と時間的に関連付けて保持する保持部と、
前記動態画像及び前記診断支援情報を表示するための表示用画像を生成する表示画像生成部と、
を備え、
前記診断支援情報は、前記解析に基づく第１の解析結果と第２の解析結果とを含み、
前記表示用画像は、
前記動態画像を表示する動態画像表示部分と、
前記診断支援情報の第１の解析結果と第２の解析結果とを識別可能に時間軸方向に一覧できるように表示する一覧表示部分と、
を含む画像であり、
前記診断支援情報生成部による解析結果が所定の条件を満たしたときに、ユーザに前記所定の条件を満たしたことを報知する報知部、
を更に備えることを特徴とする、
画像生成装置。
【請求項６】
請求項２に記載の画像生成装置であって、
前記複数の部位は、左肺野及び右肺野を含むことを特徴とする、
画像生成装置。
【請求項７】
請求項４に記載の画像生成装置であって、
前記検出部は、
前記動態画像に基づき、前記所定部位の血流の有無に関する画像処理結果を出力し、
前記表示画像生成部は、
前記画像処理結果を前記動態画像と同期して表示するように前記表示用画像を生成することを特徴とする、
画像生成装置。
【請求項８】
請求項１ないし請求項７の何れか１項に記載の画像生成装置であって、
前記表示用画像は、
前記一覧表示部分と関連して設けられ、前記一覧表示部分の時間軸方向の特定位置を示す指標を含み、
前記表示画像生成部は、
前記指標によって示された特定位置と対応する時点の動態画像を前記動態画像表示部分に表示するように前記表示用画像を生成する、ことを特徴とする、
画像生成装置。
【請求項９】
請求項１ないし請求項７の何れか１項に記載の画像生成装置であって、
前記一覧表示部分は、前記第１の解析結果と第２の解析結果とを異なる色または濃淡で表示することを特徴とする、
画像生成装置。
【請求項１０】
請求項１ないし請求項７の何れか１項に記載の画像生成装置であって、
前記動態画像表示部分は、矩形であり、
前記一覧表示部分は、前記動態画像表示部分の一辺に沿った長尺状の領域であり、その長手方向と前記診断支援情報の時間軸方向とが対応していることを特徴とする、
画像生成装置。
【請求項１１】
請求項１ないし請求項７の何れか１項に記載の画像生成装置であって、
前記表示用画像は、
前記動態画像表示部分の表示に対応した再生時刻情報を表示する再生時刻表示部分、
を更に含むことを特徴とする、
画像生成装置。
【請求項１２】
請求項１ないし請求項７の何れか１項に記載の画像生成装置であって、
前記動態画像表示部分に表示される動態画像の再生時刻調整用の再生時刻調整インタフェースを更に備える、ことを特徴とする、
画像生成装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、人体または動物の所定部位が撮影された動態画像の画像生成装置に関する。
【背景技術】
【０００２】
医療現場では、Ｘ線等を用いて内臓や骨格等に含まれる患部を撮影することにより、各種検査や診断が行われている。そして、近年では、デジタル技術の適用により、Ｘ線等を用いて患部の動きを捉えた動態画像を比較的容易に取得することが可能となっている。【０００３】
そこでは、ＦＰＤ（flat panel detector）等の半導体イメージセンサを利用し、対象部位を含む被写体領域に対し動態画像を撮影できるため、従来のＸ線撮影による静止画撮影及び診断では実施できなかった対象部位などの動き解析による診断を実施することが可能になってきた。例えば、胸部Ｘ線動態画像から肺野内の換気情報を抽出し、肺野内濃度変化や動きなどから、動態機能の定量的解析により、診断／治療を支援(Ｘ線動画用ＣＡＤ)する検討も実施されている。【０００４】
このように、ＦＰＤによる動態画像撮影では、今まで静止画でしか見られなかった体内状況について、時間変化を見てとれることによる新たな診断効果が期待されている。【０００５】
例えば、特許文献１が開示する技術では、肺のＸ線動態画像を表示するに際して、Ｘ線動態画像を処理して得た各フレームの肺尖部から横隔膜までの距離のグラフを表示し、このグラフに関連して設けられたスクロールバーを操作することで、スクロールバーの位置に対応した動態画像のフレームを表示している。【０００６】
また、特許文献２が開示する技術では、対象部位の動態画像ではないが、生体内で撮影した画像ストリームを表示するに際して、撮影位置（体内組織）や温度等のカラーバーを表示し、カラーバーのストライプを指すことにより、そのストライプに対応するフレームに画像ストリームを進めている。
【先行技術文献】
【特許文献】
【０００７】
【特許文献１】国際公開第２００６／１３７２９４号
【特許文献２】特表２００９−５０８５６７号公報
【発明の概要】
【発明が解決しようとする課題】
【０００８】
しかしながら、上記特許文献１の従来技術では、診断対象である所定部位が時間変化する変化量そのものをグラフ表示するものであり、有効なフレームを探す作業に困難性がある。【０００９】
また、上記特許文献２の従来技術では、撮影位置や温度等の測定結果と画像ストリームとを対応付けているものであり、対象部位の動態画像の解析については何ら触れられていない。【００１０】
本発明は、このような事情に鑑みてなされたものであり、人体または動物の所定部位の動態画像を表示する際に、作業性や視認性の良い画像生成装置を提供することを目的とする。
【課題を解決するための手段】
【００１１】
本発明の１側面による画像生成装置は、人体または動物の所定部位を時間順次に捉えた動態画像を取得する動態画像取得部と、前記所定部位の物理的状態の時間変化を検出する検出部と、前記検出部によって検出された前記所定部位の物理的状態の時間変化に基づいて解析を行い、その解析結果を診断支援情報として生成する診断支援情報生成部と、前記診断支援情報を前記動態画像と時間的に関連付けて保持する保持部と、前記動態画像及び前記診断支援情報を表示するための表示用画像を生成する表示画像生成部とを備え、前記診断支援情報は、前記解析に基づく第１の解析結果と第２の解析結果とを含み、前記表示用画像は、前記動態画像を表示する動態画像表示部分と、前記診断支援情報の第１の解析結果と第２の解析結果とを識別可能に時間軸方向に一覧できるように表示する一覧表示部分と、を含む画像であり、前記検出部は、前記所定部位の物理的状態の周期的な時間変化となる所定部位周期を特定する所定部位周期特定部、を備え、前記所定部位は肺野であり、前記第１の解析結果は呼気を示し、前記第２の解析結果は吸気を示す、ことを特徴とする。
【発明の効果】
【００１２】
本発明によれば、人体または動物の所定部位の動態画像を表示する際に、作業性や視認性を向上することができ、診断効率を上げることができる。
【図面の簡単な説明】
【００１３】
【図１】図１は各実施形態に係る放射線動態画像撮影システムの全体構成を示す図である。
【図２】図２は放射線動態画像撮影によって撮影した動態画像を例示する図である。
【図３】図３は参考例における表示用画像を例示する図である。
【図４】図４は第１実施形態に係る画像生成装置の機能構成を示すブロック図である。
【図５】図５は第１実施形態に係る画像生成装置の機能構成を示すブロック図である。
【図６】図６は第６実施形態に係る画像生成装置の機能構成を示すブロック図である。
【図７】図７は心電計で計測された波形の一部を例示する模式図である。
【図８】図８は画像に対する座標軸の設定状態を例示する模式図である。
【図９】図９は心臓壁の変動を例示する模式図である。
【図１０】図１０は呼吸停止中肺動態画像に係る心臓の横幅の変動サイクルを例示する模式図である。
【図１１】図１１は血流に関する表示を例示する図である。
【図１２】図１２は肺血管領域の血流信号値の時間変化を示す波形を例示する図である。
【図１３】図１３は肺野領域の輪郭抽出を例示する模式図である。
【図１４】図１４は肺野領域の特徴点の位置を例示した模式図である。
【図１５】図１５は呼吸位相の診断支援情報の生成を模式的に例示する図である。
【図１６】図１６は第１実施形態におけるに係る表示用画像を例示する図である。
【図１７】図１７は第２実施形態及び第３実施形態におけるに係る表示用画像を例示する図である。
【図１８】図１８は第１実施形態において実現される画像生成装置の基本動作を説明するフローチャートである。
【図１９】図１９は第４実施形態に係る画像生成装置の機能構成を示すブロック図である。
【図２０】図２０は心拍（血流）の波形データを時系列で示した模式図である。
【図２１】図２１は呼吸情報の波形データを時系列で示した模式図である。
【図２２】図２２は第４実施形態において実現される画像生成装置の基本動作を説明するフローチャートである。
【図２３】図２３は第５実施形態に係る画像生成装置の機能構成を示すブロック図である。
【図２４】図２４は第５実施形態において実現される画像生成装置の基本動作を説明するフローチャートである。
【図２５】図２５は関節屈伸における関節角度を例示する図である。
【図２６】図２６は屈伸位相の診断支援情報の生成を模式的に例示する図である。
【図２７】図２７は診断支援情報として血流情報を表示した場合の表示用画像を例示する図である。
【図２８】図２８は本実施形態に係る画像生成装置の機能構成の一部を例示するブロック図である。
【発明を実施するための形態】
【００１４】
＜１．第１実施形態＞
本発明の第１実施形態に係る放射線動態画像撮影システムについて以下説明する。
【００１５】
＜１−１．放射線動態画像撮影システムの全体構成＞
第１実施形態に係る放射線動態画像撮影システムは、人体または動物の身体を被写体として、被写体の放射線画像の撮影を行い、所望の表示用画像を生成する。
【００１６】
図１は、第１実施形態に係る放射線動態画像撮影システムの全体構成を示す図である。図１に示すように、放射線動態画像撮影システム１００は、撮影装置１と、撮影制御装置２（撮影用コンソール）と、画像生成装置３（診断用コンソール）と、心電計４とを備える。撮影装置１及び心電計４と、撮影制御装置２とが通信ケーブル等により接続され、撮影制御装置２と、画像生成装置３とがＬＡＮ（Local Area Network）等の通信ネットワークＮＴを介して接続されて構成されている。放射線動態画像撮影システム１００を構成する各装置は、ＤＩＣＯＭ（Digital Image and Communications in Medicine）規格に準じており、各装置間の通信は、ＤＩＣＯＭ規格に則って行われる。
【００１７】
＜１−１−１．撮影装置１の構成＞
撮影装置１は、例えば、Ｘ線撮影装置等によって構成され、呼吸に伴う被写体Ｍの胸部の動態を撮影する装置である。動態撮影は、被写体Ｍの胸部に対し、Ｘ線等の放射線を繰り返して照射しつつ、時間順次に複数の画像を取得することにより行う。この連続撮影により得られた一連の画像を動態画像と呼ぶ。また、動態画像を構成する複数の画像のそれぞれをフレーム画像と呼ぶ。
【００１８】
図１に示すように、撮影装置１は、照射部（放射線源）１１と、放射線照射制御装置１２と、撮像部（放射線検出部）１３と、読取制御装置１４と、サイクル検出センサ１５と、サイクル検出装置１６とを備えて構成されている。
【００１９】
照射部１１は、放射線照射制御装置１２の制御に従って、被写体Ｍに対し放射線（Ｘ線）を照射する。図示例は人体用のシステムであり、被写体Ｍは検査対象者に相当する。以下では被写体Ｍを「被検者」とも呼ぶ。
【００２０】
放射線照射制御装置１２は、撮影制御装置２に接続されており、撮影制御装置２から入力された放射線照射条件に基づいて照射部１１を制御して放射線撮影を行う。
【００２１】
撮像部１３は、ＦＰＤ等の半導体イメージセンサにより構成され、照射部１１から照射されて被検者Ｍを透過した放射線を電気信号（画像情報）に変換する。
【００２２】
読取制御装置１４は、撮影制御装置２に接続されている。読取制御装置１４は、撮影制御装置２から入力された画像読取条件に基づいて撮像部１３の各画素のスイッチング部を制御して、当該各画素に蓄積された電気信号の読み取りをスイッチングしていき、撮像部１３に蓄積された電気信号を読み取ることにより、画像データを取得する。そして、読取制御装置１４は、取得した画像データ（フレーム画像）を撮影制御装置２に出力する。画像読取条件は、例えば、フレームレート、フレーム間隔、画素サイズ、画像サイズ（マトリックスサイズ）等である。フレームレートは、１秒あたりに取得するフレーム画像数であり、パルスレートと一致している。フレーム間隔は、連続撮影において、１回のフレーム画像の取得動作開始から次のフレーム画像の取得動作開始までの時間であり、パルス間隔と一致している。
【００２３】
ここで、放射線照射制御装置１２と読取制御装置１４とは互いに接続され、互いに同期信号をやりとりして放射線照射動作と画像の読み取りの動作を同調させるようになっている。
【００２４】
サイクル検出装置１６は、被写体Ｍの呼吸サイクルを検出してサイクル情報を撮影制御装置２の制御部２１に出力する。また、サイクル検出装置１６は、例えばレーザー照射によって被写体Ｍの胸部の動き（被写体Ｍの呼吸サイクル）を検出するサイクル検出センサ１５と、サイクル検出センサ１５により検出された呼吸サイクルの時間を測定し制御部２１に出力する計時部（不図示）とを備える。
【００２５】
＜１−１−２．撮影制御装置２の構成＞
撮影制御装置２は、放射線照射条件や画像読取条件を撮影装置１に出力して撮影装置１による放射線撮影及び放射線画像の読み取り動作を制御するとともに、撮影装置１により取得された動態画像を撮影技師によるポジショニングの確認や診断に適した画像であるか否かの確認用に表示する。
【００２６】
図１に示すように、撮影制御装置２は、制御部２１と、記憶部２２と、操作部２３と、表示部２４と、通信部２５とを備えて構成され、各部はバス２６により接続されている。
【００２７】
制御部２１は、ＣＰＵ（Central Processing Unit）やＲＡＭ（Random Access Memory）等により構成される。制御部２１のＣＰＵは、操作部２３の操作に応じて、記憶部２２に記憶されているシステムプログラムや各種処理プログラムを読み出してＲＡＭ内に展開し、展開されたプログラムに従って後述する撮影制御処理を始めとする各種処理を実行し、撮影制御装置２各部の動作や、撮影装置１の動作を集中制御する。
【００２８】
記憶部２２は、不揮発性の半導体メモリやハードディスク等により構成される。記憶部２２は、制御部２１で実行される各種プログラムやプログラムにより処理の実行に必要なパラメータ、或いは処理結果等のデータを記憶する。
【００２９】
操作部２３は、カーソルキー、数字入力キー、及び各種機能キー等を備えたキーボードと、マウス等のポインティングデバイスとを備えて構成され、キーボードに対するキー操作、マウス操作、あるいは、タッチパネルを介して入力された指示信号を制御部２１に出力する。
【００３０】
表示部２４は、カラーＬＣＤ（Liquid Crystal Display）等のモニタにより構成され、制御部２１から入力される表示信号の指示に従って、操作部２３からの入力指示やデータ等を表示する。
【００３１】
通信部２５は、ＬＡＮアダプタやモデムやＴＡ（Terminal Adapter）等を備え、通信ネットワークＮＴに接続された各装置との間のデータ送受信を制御する。
【００３２】
＜１−１−３．画像生成装置３の構成＞
画像生成装置３は、撮像装置１から送信された動態画像を、撮影制御装置２を介して取得し、医師等が読影診断するための画像を表示する。
【００３３】
図１に示すように、画像生成装置３は、制御部３１と、記憶部３２と、操作部３３と、表示部３４と、通信部３５とを備えて構成され、各部はバス３６により接続されている。
【００３４】
制御部３１は、ＣＰＵ、ＲＡＭ等により構成される。制御部３１のＣＰＵは、操作部３３の操作に応じて、記憶部３２に記憶されているシステムプログラムや、各種処理プログラムを読み出してＲＡＭ内に展開し、展開されたプログラムに従って各種処理を実行し、画像生成装置３各部の動作を集中制御する（詳細は後述する）。
【００３５】
記憶部３２は、不揮発性の半導体メモリやハードディスク等により構成される。記憶部３２は、制御部３１で実行される各種プログラムやプログラムにより処理の実行に必要なパラメータ、或いは処理結果等のデータを記憶する。例えば、記憶部３２は、後述する画像生成処理を実行するための画像生成処理プログラムを記憶している。これらの各種プログラムは、読取可能なプログラムコードの形態で格納され、制御部３１は、当該プログラムコードに従った動作を逐次実行する。
【００３６】
操作部３３は、カーソルキー、数字入力キー、及び各種機能キー等を備えたキーボードと、マウス等のポインティングデバイスを備えて構成され、キーボードに対するキー操作やマウス操作、あるいは、タッチパネルを介して入力された指示信号を制御部３１に出力する。
【００３７】
表示部３４は、カラーＬＣＤ等のモニタにより構成され、制御部３１から入力される表示信号の指示に従って、操作部３３からの入力指示、データ、及び、後述する表示用画像を表示する。
【００３８】
通信部３５は、ＬＡＮアダプタやモデムやＴＡ等を備え、通信ネットワークＮＴに接続された各装置との間のデータ送受信を制御する。
【００３９】
＜１−１−４．心電計４の構成＞
図１では心電計４は被検者Ｍとは離れて示されているが、実際には心電計４の各電極端子は被検者Ｍに装着されており、被検者Ｍの心電波形をデジタル信号として出力する。
【００４０】
図１に示すように、心電計４は、位相検出部４１を備えて構成され、位相検出部４１は、制御部２１のＣＰＵからの制御信号に応答して、撮影装置１による撮影動作を同期させるための基礎情報として、被写体Ｍの心拍（血流）の位相を検出する。なお、本願発明においては、心拍も血流の一種として扱う。また、位相検出部４１は、撮影制御装置２内に設けることも可能である。
【００４１】
＜１−２．参考例の動画診断における問題点＞
この実施形態における画像生成装置３の詳細を説明する前提として、参考例における動画診断の問題点を説明しておく。
【００４２】
図２は、呼吸に伴う被検者Ｍの胸部の動態に対し、放射線動態画像撮影によって撮影した動態画像であり、図３は、参考例において生成された表示用画像を例示する図である。なお、図３の表示用画像ＩＧ０では、フレーム画像ＭＩと、当該フレーム画像ＭＩと同期した、横隔膜の位置を示すグラフＧ１と胸郭の幅を示すグラフＧ２と呼吸情報を示すグラフＧ３と心拍情報（血流情報）を示すグラフＧ４とがフレーム画像ＭＩの片側に合わせて表示されている。
【００４３】
図２で示されるように、フレーム画像Ｍ１〜Ｍ１０は、呼吸サイクルの１周期（詳細は後述する）を一定の撮影タイミングで連続撮影して取得されたものである。具体的には、時刻 ｔ＝ｔ1, ｔ2, ｔ3, …, ｔ10 の撮影タイミングにおいて撮影された画像が、フレーム画像Ｍ１，Ｍ２，Ｍ３，…，Ｍ１０にそれぞれ対応している。また、図２におけるフレーム画像Ｍ７，Ｍ８，Ｍ９が、医者等の専門家であるユーザが探知するための診断において重要なフレーム画像であり、領域ＩＲ（図２及び図３参照）が、診断対象領域に相当するものとし、現時点でのグラフＧ１の位相位置が、ラインＬＴで示されている（図３参照）。
【００４４】
この参考例における動画診断では、図３で示されるように、動態画像の再生の進行に対応したプログレスバーＰＢを操作しながら、所望のフレーム画像ＭＩを探し出し、当該フレーム画像ＭＩを詳細に観察することで診断する。ところが、図３の表示用画像ＩＧ０では、医者は動画に注視したいにもかかわらず、詳細な診断データ（ここでは、グラフＧ１及びラインＬＴ）をも常に注視しなければならないため、視線運動が多くなり、診断効率が極めて悪い。
【００４５】
実際の動画診断では、例えば、１０ Frames Per Second(fps)、３０秒間で肺野動画撮影した場合では、全体で３００枚程度の静止画から診断に関係するフレーム時刻を探す作業と、該当フレーム時刻付近を動画再生して閲覧診断する作業との両方を行う必要があり、非常に手間を要する。特に医者は静止画の読影を日々１００枚（１００人）分診断するというだけでも負荷がある中、動画では１人３００枚の読影をしなければならないとなると、動画読影による診断は困難を極める。
【００４６】
特に、後述する血流動画診断では、血流の一位相（一心拍）分のフレームだけを表示したとしても、十分な情報が得られないことが多い。この理由としては、(i)被爆を抑えるため低線量で撮影されることで画質が悪くなる、(ii)fpsの不足のために肺塞栓疑いを判断できるような最適なタイミングの撮影ができているか不明である、(iii)呼吸や体動の影響で血流の見え方が各時刻で変わるため所望の血流情報がわからない等の要因が挙げられる。したがって、動画全体の血流複数周期の中から該当位相のフレーム（図２の例では、フレーム画像Ｍ７，Ｍ８，Ｍ９に相当）を複数探しだし診断しなければならず、作業負荷が極めて高い現状である。
【００４７】
このような背景の下、簡易に操作者が所望のフレーム画像ＭＩを選択操作することができ、且つ、診断データへの注視を少なくすることで視線運動を少なくし、本来診断すべき動態画像に注力して注視し続けることが望まれている。
【００４８】
そこで、本実施形態では、所望の画像を把握するための診断支援情報が見やすい表示用画像を生成することで、動画読影による診断時間の短縮化を図る。
【００４９】
以下では、第１実施形態における画像生成装置３の具体的構成について説明する。
【００５０】
＜１−３．画像生成装置３の具体的構成＞
本発明の第１実施形態における放射線動態画像撮影システム１００の画像生成装置３は、被検者Ｍの心臓または肺等（所定部位）の周期的な時間変化に基づく状態変化を診断支援情報として時間軸方向にわかりやすく表示した表示用画像を生成することで、診断に関係する所望のフレーム画像ＭＩを探す作業を軽減する。
【００５１】
以下では、画像生成装置３で実現される機能的な構成について説明する。
【００５２】
＜１−３−１．画像生成装置３の機能構成＞
図４及び図５は、放射線動態画像撮影システム１００における画像生成装置３において、ＣＰＵ等が各種プログラムに従って動作することにより制御部３１で実現される機能構成を他の構成とともに示す図である。この実施形態の画像生成装置３は、主として心臓および両肺を含む胸部が撮影された動態画像を使用する。
【００５３】
制御部３１では、主に、動態画像取得部１１０と、検出部１２０と、診断支援情報生成部１５０と、表示画像生成部１６０と、から構成される。
【００５４】
表示部３４では、表示画像生成部１６０において生成された表示用画像を表示するとともに、ユーザが再生時刻を変化させて表示用画像を参照可能な再生時刻調整部３４１を備えて構成される。
【００５５】
以下では、図４及び図５で示されたような制御部３１の機能的な構成が、あらかじめインストールされたプログラムの実行によって、実現されるものとして説明するが、専用のハードウエア構成で実現されても良い。
【００５６】
以降、動態画像取得部１１０、検出部１２０、診断支援情報生成部１５０、及び、表示画像生成部１６０が行う各処理についての具体的内容を、図４及び図５を参照しながら順次説明する。
【００５７】
＜１−３−１−１．動態画像取得部１１０＞
動態画像取得部１１０では、撮像装置１の読取制御装置１４によって撮影された人体または動物の対象部位（所定部位）を時間順次に捉えた動態画像を取得する。ここで、対象部位とは、図４では心臓領域、図５では肺領域であるものとする。
【００５８】
なお、図４及び図５は、撮像装置１と画像生成装置３とが直接接続されて示されているが、実際は、撮影制御装置２が介在し、撮影制御装置２の記憶部２２にて記憶された処理データが通信部２５を介して、画像生成装置３の通信部３５に出力される。
【００５９】
＜１−３−１−２．検出部１２０＞  
検出部１２０は、所定部位周期特定部１３０を備え、対象部位の物理的状態の時間変化を検出する（図４及び図５参照）。ここでいう「物理的状態」という用語は、心臓や肺等の幾何学的形状を指すほか、血流の濃度（血流の有無）などをも包含した意味で用いている。被検者Ｍの心臓または肺の周期的な時間変化、すなわち、血流（心拍も含む）、及び、呼吸の位相情報や周波数（周期）情報を検出する。ここにおける時間変化の検出とは、臓器の外形や血流の濃度などの部位周期物理的状態についての時間的変化の検出を意味する。
【００６０】
所定部位周期特定部１３０では、対象部位の物理的状態の周期的な時間変化となる対象部位周期（所定部位周期）を特定する。
【００６１】
以下では、血流および呼吸による位相情報の算出方法として、本実施形態で用いる、第１の血流情報検出方法と第１及び第２の呼吸情報検出方法とについて説明する。
【００６２】
○  第１の血流情報検出方法：心電計の検出結果
第１の血流情報検出方法として、図４で示されるように検出部１２０（所定部位周期特定部１３０）では、心電計４の位相検出部４１から取得された結果を用いる方法である。図７は、被検者Ｍの心電図波形を例示する図である。なお、図７では、横軸が時刻、縦軸が電気信号の大きさ（電圧）を示しており、いわゆるＰ波、Ｑ波、Ｒ波、Ｓ波、Ｔ波及びＵ波の形状をそれぞれ示す曲線Ｐｐ，Ｑｐ，Ｒｐ，Ｓｐ，Ｔｐ及びＵｐを含む電気信号の変化を示す曲線が示されている。
【００６３】
そこで、所定部位周期特定部１３０では、位相検出部４１から取得された検出結果に基づいて、上記の点（Ｐｐ，Ｑｐ，Ｒｐ，Ｓｐ，Ｔｐ及びＵｐ）を解析することで、心拍（血流）周期を特定する。
【００６４】
なお、位相検出部４１による検出動作は撮像装置１による撮像動作と同期して行われる（図１参照）。また、図４は、心電計４と画像生成装置３とが直接接続されて示されているが、実際は、撮影制御装置２が介在し、撮影制御装置２の記憶部２２にて記憶された検出データが通信部２５を介して、画像生成装置３の通信部３５に出力される。
【００６５】
このように、所定部位周期特定部１３０では、外部より心拍（血流）周期が設定可能に構成されることにより、対象部位の周期的な時間変化を自動的に取得することが可能となる。
【００６６】
例えば、第１の血流情報検出方法では、上述の血流情報を示すグラフＧ４を生成するための情報を得ることが可能となる。
【００６７】
○  第１の呼吸情報検出方法：別機器による計測結果
第１の呼吸情報検出方法として、別機器による計測で実施する。別機器により計測する方法としては、例えば、特許第3793102号に記載されているような装置を用いることができる。また、レーザー光とＣＣＤカメラで構成されたセンサによるモニタリングにより実施する手法（例えば、"FG視覚センサを用いた就寝者の呼吸モニタリングに関する検討",青木  広宙,中島  真人,電子情報通信学会ソサイエティ大会講演論文集 2001年.情報・システムソサイエティ大会講演論文集, 320-321, 2001-08-29.等参照）等を採用することもできる。
【００６８】
なお、本実施形態では、図５で示されるように検出部１２０（所定部位周期特定部１３０）では、サイクル検出装置１６のサイクル検出センサ１５が利用可能である。
【００６９】
また、呼吸サイクルを検出するための別の方法としては、呼吸モニタベルトを用いて被写体の胸部の動きを検出する方法や、気速計により呼吸の気流を検出する方法があり、これらの方法を適用することも可能である。
【００７０】
例えば、第１の呼吸情報検出方法では、上述の呼吸情報を示すグラフＧ３と後述のカラーバーＣ１とを生成するための情報を得ることが可能となる。
【００７１】
このように、所定部位周期特定部１３０では、外部より呼吸周期が設定可能に構成されることにより、対象部位の周期的な時間変化を自動的に取得することが可能となる。
【００７２】
○  第２の呼吸情報検出方法：面積値または特徴点間距離
一方、第２の呼吸情報検出方法として、動態画像取得部１１０によって取得された撮影画像を用いて、肺野部の面積値を算出することで、呼吸情報とするものである。肺野部の面積の求め方は、肺野部の輪郭抽出を行い、輪郭に囲まれた領域の画素数を肺野領域として定義することが可能である。すなわち、横隔膜の位置や胸郭の幅を検出することで、呼吸情報を得ることができる。
【００７３】
例えば、第２の呼吸情報検出方法では、上述の横隔膜の位置を示すグラフＧ１及び胸郭の幅を示すグラフＧ２を生成するための情報を得ることが可能となる。
【００７４】
図１３は、肺野部の輪郭抽出を例示する模式図である。肺野部の抽出は、図１３で示すように、左右ごとに抽出しても、心臓や脊椎の領域を含んだ輪郭ＯＬとして抽出してもよい。抽出方法としては、従来技術（例えば、“Image feature analysis and computer-aided diagnosis: Accurate determination of ribcage boundary in chest radiographs”, Xin-Wei Xu and Kunio Doi, Medical Physics, Volume 22(5), May 1995, pp.617-626.等参照）等を採用することができる。
【００７５】
このように、所定部位周期特定部１３０では、取得された撮影画像を用いて、肺野部の輪郭ＯＬの抽出を実施し、抽出された領域内のピクセル数を特徴量として、該特徴量を肺野領域の面積として検出する。
【００７６】
また、第２の呼吸情報検出方法としては、動態画像取得部１１０によって取得された撮影画像を用いて、肺野領域の特徴点間の距離を算出し、呼吸情報とすることも可能である。すなわち、肺野部の抽出を上記方法と同様に実施し、抽出された領域から、特徴点２点を求め、その２点間の距離を求めることで特徴量を算出する。
【００７７】
図１４は、肺野領域の特徴点の位置を例示した図である。肺領域の上端ＬＴから下端ＢＴまでの長さ（肺野長）の変化を算出する場合、図１４（ａ）では、肺尖部を肺領域の上端ＬＴとし、肺尖部から体軸方向におろした直線と横隔膜との交点を肺領域の下端ＬＢとして抽出した例であり、図１４（ｂ）では、肺尖部を肺領域の上端ＬＴとし、肋横角を肺領域の下端ＬＢとして抽出した例である。
【００７８】
このように、所定部位周期特定部１３０では、取得された撮影画像を用いて、肺野領域の輪郭ＯＬの抽出を実施し、抽出された領域から特徴点間の距離を求めることで、特徴点間距離を検出することで、呼吸周期を設定する。
【００７９】
以上のように、所定部位周期特定部１３０では、動態画像で捉えられた肺野領域の面積値あるいは特徴点間距離の時間的変化（所定部位の形状の変化）に基づき、呼吸周期を検出するため、呼吸周期を自動的に取得することが可能となる。
【００８０】
また、第２の呼吸情報検出方法は、第１の呼吸情報検出方法と比較して間接的に呼吸周期を検出するため、ノイズ成分も含まれやすいと予想される。そこで、所定部位周期特定部１３０では、動態画像で捉えられた肺野領域の面積値あるいは特徴点間距離の時間的変化（所定部位の形状の変化）に基づき、周波数解析等を用いて呼吸周期を検出することが好ましい。これにより、ノイズ成分を取り除いた所望の変動成分を自動的に抽出することができるため、肺野領域の面積値あるいは特徴点間距離の時間的変化（所定部位が時間変化する状態）をより正確に把握することが可能となる。
【００８１】
＜１−３−１−３．診断支援情報生成部１５０＞
診断支援情報生成部１５０では、上述の検出部１２０（所定部位周期特定部１３０）によって検出された心臓または肺などの対象部位の物理的状態の時間変化に基づいて解析を行い、その解析結果を診断支援情報として生成する。
【００８２】
また、診断支援情報とは、解析に基づく第１の解析結果と第２の解析結果とを含んで構成される。ここで、例えば、対象部位が肺野である場合、第１の解析結果は呼気を示し、第２の解析結果は吸気を示す。また、診断支援情報は、第１の解析結果と第２の解析結果とに加えて、第３の解析結果を含んでも構わない。例えば、対象部位が肺野である場合、第３の解析結果としては呼吸停止状態が挙げられる。
【００８３】
なお、解析結果の記憶部（保持部）３２への保持方法に関しては、動態画像にメタデータとして保持される方法、すなわち、診断支援情報生成部１５０において生成された診断支援情報を動態画像（フレーム画像ＭＩ）と時間的に関連付けて、図２８で示されるように、保持部３２に保持する方法でも良いし、別データベースとして管理される方法であっても良い。これにより、診断時に測定や画像特徴量を計算する必要がなくなり望ましい。また、診断支援情報として表示する可能性のある重要な状態変化情報に限定して保持されていてもよい。
【００８４】
以下では、呼吸情報および血流情報における診断支援情報について説明する。
【００８５】
なお、呼吸の周期（呼吸サイクル）Ｂの１周期は、吸気と呼気とから構成され、１回の呼気と１回の吸気とからなる。吸気では、横隔膜が下がって息が吸い込まれるに連れて胸郭中の肺野の領域が大きくなる。息を最大限に吸い込んだとき（吸気と呼気の変換点）が最大吸気時Ｂ１である。呼気では、横隔膜が上がって息が吐き出されるに連れて肺野の領域が小さくなるが、息を最大限に排出したとき（呼気と吸気の変換点）が最大呼気時Ｂ２となる（図２１参照）。
【００８６】
図１５は、所定部位周期特定部１３０において特定された呼吸情報の周期に対して、呼吸位相の診断支援情報の生成を模式的に例示する図である。図１５で示されるように、ここでは、第１の解析結果は呼気を示し、第２の解析結果は吸気を示し、呼吸位相を時間変化に対応させ、呼吸周期において識別可能な診断支援情報を生成する。
【００８７】
なお、本実施形態では呼吸情報における診断支援情報を用いて表示を行うが、血流情報における診断支援情報を用いて表示を行うことも可能である（図２７参照）。この場合、図１１で示される血流情報では、第１の解析結果は血流の「有り」を示し、第２の解析結果は血流の「無し」を示し、血流位相を時間変化に対応させ、血流周期において識別可能な診断支援情報を生成する。
【００８８】
これにより、後述の表示用画像における診断支援情報を参照して、呼吸周期における呼気と吸気との診断情報の変化が識別可能となるため、呼吸周期の肺野領域の状態変化に対して、効率的な医療診断が可能となる。また、解析結果が、心拍（血流）情報（心臓の「収縮期」「拡張期」が第１及び第２の解析結果）である場合においても、同様に心拍（血流）の周期の状態変化に対して、効率的な医療診断が可能となる。
【００８９】
なお、診断支援情報は、診断内容が予め決まっているケースでは、固定しておけば効率的である。たとえば、肺塞栓の疑いがある患者であれば、肺塞栓に有効な肺血流位相を採用すればよい。また、呼吸に異常疑いがある患者であれば、呼吸診断に有効な呼吸位相を採用する他、呼吸器系の異常パターンが幾つか解析可能である場合には、その状態変化が見受けられる複数の診断情報を採用してもよい。
【００９０】
＜１−３−１−４．表示画像生成部１６０＞
表示画像生成部１６０では、フレーム画像ＭＩ（動態画像）及び診断支援情報を表示するための表示用画像を生成する。すなわち、対象部位の位相変化と時間的に対応するフレーム画像ＭＩとを対応付けて、表示用画像が生成される。
【００９１】
また、表示画像生成部１６０は、動態画像を表示する動態画像表示部分１６１と、診断支援情報の第１の解析結果と第２の解析結果とを識別可能に時間軸方向に一覧できるように表示する一覧表示部分１６２と、動態画像表示部分１６１の表示に対応した再生時刻情報を表示する再生時刻表示部分１６３と、が設けられた表示用画像を生成する（図４及び図５参照）。すなわち、表示用画像は、一覧表示部分１６２を含み、一覧表示部分１６２の時間軸方向の特定位置を示す指標を含んで生成される。そして、表示画像生成部１６０は、当該指標によって示された特定位置と対応する時点の動態画像を動態画像表示部分１６１に表示するように表示用画像を生成する。したがって、ある時点での表示用画像において、動態画像表示部分１６１に表示される動態画像（フレーム画像）と、指標によって示される一覧表示部分１６２の時間軸方向の特定位置とは対応関係にある。また、動態画像が動態画像表示部分１６１に動画的な態様で表示される場合においても、一覧表示部分１６２に表示される診断支援情報は静止画的な態様で表示され、一覧表示部分１６２の指標は時間軸方向に移動するように表示される。
【００９２】
また、動態画像表示部分１６１は、矩形であり、一覧表示部分１６２は、動態画像表示部分１６１の一辺に沿った長尺状の領域であり、その長手方向と診断支援情報の時間軸方向とが対応している。
【００９３】
さらに、表示画像生成部１６０では、対象部位が複数の部位であるとき、検出部１２０は、複数の対象部位の物理的状態の時間変化をそれぞれ検出し、診断支援情報生成部１５０は、複数の対象部位の物理的状態の時間変化に基づいてそれぞれ解析を行い、複数の対象部位についての解析結果を複数の診断支援情報としてそれぞれ生成し、一覧表示部分１６２は、該複数の診断支援情報を表示する。
【００９４】
図１６は、表示画像生成部１６０において生成された表示用画像ＩＧが、表示部３４のディスプレイの画面上に表示された例である。この図１６で示されるように、表示用画像ＩＧには被検者Ｍについて撮影されたフレーム画像ＭＩや図３と共通のグラフＧ１〜Ｇ４のほか、グラフィック要素として、
・肺野の呼吸位相に応じたカラーバーＣ１，
・プログレスバーＰＢ，
・再生時刻表示部ＴＭ
が並列的に表示されている。
【００９５】
なお、本実施形態において、呼吸情報を示すグラフＧ３とカラーバーＣ１は、第１の呼吸情報検出方法から得られた情報であり、横隔膜の位置を示すグラフＧ１及び胸郭の幅を示すグラフＧ２は、第２の呼吸情報検出方法から得られた情報であり、血流情報を示すグラフＧ４は、第１の血流情報検出方法から得られた情報である。
【００９６】
なお、図１６に示される表示用画像ＩＧでは、フレーム画像ＭＩを表示する部分が動態画像表示部分１６１に相当し、カラーバーＣ１の部分が一覧表示部分１６２に相当し、再生時刻表示部ＴＭの部分が再生時刻表示部分１６３に相当する。
【００９７】
図１６で示されるように、図３で示された参考例における表示用画像ＩＧ０と異なり、プログレスバーＰＢとカラーバーＣ１とが統合して表示される。
【００９８】
また、一覧表示部分１６２は、呼気及び吸気等の第１の解析結果と第２の解析結果とを異なる色表示（例えば、２色で簡易表示する等）または濃淡（いわゆる、グラデーション）で表示する表示形態であってもよい。色情報として扱う色の選択としては、特に状態変化量を表示する場合は認識しやすくするために状態変化量を正規化するとよい。状態変化の違いの表現方法としては、輝度の違い、色相の違い、彩度の違い等で表現でき、複数の状態変化を、輝度と色相、Ｒ-ＧとＢ-Ｙで表現してもよい。これにより、位相の度合いをより明確に表現することができ、必要に応じて細かい把握も可能となる。したがって、表示用画像ＩＧを参照することにより、対象部位の診断内容が視認しやすくなり、診断効率がより向上する。
【００９９】
また、再生時刻表示部分１６３（すなわち、再生時刻表示部ＴＭ）は、一覧表示部分１６２（すなわち、カラーバーＣ１）に近接させて認識可能にした表示形態を採用することが望ましい。これによって、再生時における診断支援情報と再生時間とを視覚的に同時に把握することが可能となる。
【０１００】
さらに、一覧表示部分１６２は、動態画像表示部分１６１に表示される動態画像の再生時刻調整用の再生時刻調整インタフェース部分（プログレスバーＰＢに相当）が設けられてもよい（図４及び図５，図１６参照）。
【０１０１】
＜１−４．画像生成装置３の基本動作＞
図１８は、本実施形態に係る画像生成装置３において実現される基本動作を説明するフローチャートである。既に各部の個別機能の説明は行ったため（図４及び図５参照）、ここでは全体の流れのみ説明する。
【０１０２】
図１８に示すように、まず、ステップＳ１において、制御部３１の動態画像取得部１１０が、撮像装置１の読取制御装置１４によって撮影された動態画像を、撮影制御装置２を介して取得する。
【０１０３】
ステップＳ２では、検出部１２０が、心臓や肺等の物理的状態の時間変化を検出し、所定部位周期特定部１３０が、血流、呼吸等の周期を特定する。具体的に、血流情報の時間変化に関しては、検出部１２０（所定部位周期特定部１３０）が、心電計４の位相検出部４１から取得された結果（第１の血流情報検出方法）に基づき検出する（図４参照）。また、呼吸情報の時間変化に関しては、検出部１２０（所定部位周期特定部１３０）が、サイクル検出センサ１５から取得された結果（第１の呼吸情報検出方法）に基づき検出するか（図７〜図１４参照）（図５参照）、あるいは、動態画像におけるフレーム画像ＭＩの画像特徴量（第２の呼吸情報検出方法）に基づき検出する。
【０１０４】
ステップＳ３では、診断支援情報生成部１５０が、ステップＳ２において取得された心臓や肺などの物理的状態の時間変化に基づいて解析を行い、その解析結果を、時間変化に対応させた診断支援情報を生成する（図１１，図１５参照）。そして、診断支援情報を動態画像と時間的に関連付けて保持部３２に保持する。
【０１０５】
ステップＳ４では、表示画像生成部１６０が、ステップＳ３において保持されたフレーム画像ＭＩ（動態画像）及び診断支援情報を表示するための表示用画像ＩＧを生成する（図１６参照）。
【０１０６】
最後に、ステップＳ５において、表示画像生成部１６０が、ステップＳ４において生成した表示用画像ＩＧを表示部３４にて出力することで、表示部３４のモニタに表示し、本動作フローが終了される。
【０１０７】
以上、画像生成装置３では、人体または動物における心臓や肺等の対象部位を時間順次に捉えた動態画像を取得し、当該対象部位の物理的状態の時間変化を検出する。そして、当該対象部位の物理的状態の時間変化に基づいて解析を行い、その解析結果を診断支援情報として生成する。続いて、当該診断支援情報を動態画像と時間的に関連付けて保持した後、当該動態画像及び前記診断支援情報を表示するための表示用画像を表示する。ここで、当該診断支援情報は、解析に基づく第１の解析結果と第２の解析結果とを含み、表示用画像は、表示用画像は、動態画像を表示する動態画像表示部分１６１（図１６におけるフレーム画像ＭＩの表示部分）と、当該診断支援情報の第１の解析結果と第２の解析結果とを識別可能に時間軸方向に一覧できるように表示する一覧表示部分１６２（図１６におけるカラーバーＣ１等）と、を含む画像である。したがって、参考例において、時間変化する対象部位の「位置」変化によって専門家により判断する必要があった（図３のラインＬＴを参照）が、本実施形態（以下の第４，第５実施形態も含む）では、診断支援情報の解析結果における複数種の診断情報（例えば、解析結果が呼吸情報である場合は「呼気」「吸気」）を、時間軸方向に視覚的に識別可能に表示するようにしたため（図１６のカラーバーＣ１を参照）、所望の表示用画像の把握において、一覧性を持たせた可視化が実現可能となる。これにより、動画読影による診断時間の短縮化が図れ、ユーザにとっての利便性が向上する。
【０１０８】
＜２．第２実施形態＞
以下では、第２実施形態について説明する。第１実施形態と異なる点は、表示画像生成部１６０にて生成される表示用画像が異なる点である。なお、残余の構成は第１実施形態の画像生成装置３と同様である。
【０１０９】
図１７（ａ）は、第２実施形態における表示用画像ＩＧを示す図であり、表示部３４のディスプレイの画面上に表示された例である。図１７（ａ）で示されるように、表示用画像ＩＧには被検者Ｍについて撮影されたフレーム画像ＭＩや図３と共通のグラフＧ１〜Ｇ４のほか、グラフィック要素として、図１６と同様に、
・肺野の呼吸位相に応じたカラーバーＣ１，
・プログレスバーＰＢ，
・再生時刻表示部ＴＭ，
が表示される一方、第２実施形態における表示用画像ＩＧでは、
・位相変化を示す波形グラフＦ，
・再生時刻調整部３４１
が並列的に表示されている。
【０１１０】
なお、第２実施形態において、カラーバーＣ１と波形グラフＦとは、第１の呼吸情報検出方法から得られた情報であり、残余のグラフィック要素は、第１実施形態と同様である。
【０１１１】
具体的に、第２実施形態における表示用画像ＩＧは、図１７（ａ）で示されるように、対象部位の物理的状態の時間変化を示す波形グラフＦをカラーバーＣ１上に統合して表示される。
【０１１２】
また、一覧表示部分１６２は、動態画像表示部分１６１の再生時刻調整用の再生時刻調整インタフェース部分（プログレスバーＰＢに相当）が設けられるとともに、表示部３４は、ユーザが操作部３５を介して再生時刻調整インタフェース部分を用いて再生時刻を変化させて表示用画像ＩＧを参照可能な再生時刻調整部３４１を備える（図１７（ａ）参照）。これにより、ユーザがプログレスバーＰＢを用いて再生時刻を変化させて表示用画像ＩＧを参照することができるため、表示用画像ＩＧの所望の再生時刻にアクセスすることが可能となる。
【０１１３】
＜３．第３実施形態＞
以下では、第３実施形態について説明する。第１実施形態と異なる点は、表示画像生成部１６０にて生成される表示用画像が異なる点である。また、第３実施形態における表示用画像は第２実施形態における表示用画像とも異なる。なお、残余の構成は第１実施形態の画像生成装置３と同様である。
【０１１４】
図１７（ｂ）は、第３実施形態における表示用画像ＩＧを示す図であり、表示部３４のディスプレイの画面上に表示された例である。図１７（ｂ）で示されるように、表示用画像ＩＧには被検者Ｍについて撮影されたフレーム画像ＭＩや図３と共通のグラフＧ１〜Ｇ４のほか、グラフィック要素として、図１６及び図１７（ａ）と同様に、
・プログレスバーＰＢ，
・再生時刻表示部ＴＭ，
・再生時刻調整部３４１，
が表示される一方、第３実施形態における表示用画像ＩＧでは、
・右肺野の呼吸位相に応じたカラーバーＣ１，
・左肺野の呼吸位相に応じたカラーバーＣ２，
が並列的に表示されている。
【０１１５】
なお、第３実施形態において、カラーバーＣ１，Ｃ２は、第２の呼吸情報検出方法から右肺野及び左肺野を別々に検出することで得られた情報であり、残余のグラフィック要素は、第１実施形態と同様である。
【０１１６】
第３実施形態における表示用画像ＩＧは、図１７（ｂ）で示されるように、図３で示された参考例における表示用画像ＩＧ０と異なり、プログレスバーＰＢとカラーバーＣ１，Ｃ２とが統合して表示される。
【０１１７】
そして、図１７（ｂ）で示されるように、診断支援情報として、例えば、左肺野及び右肺野等のように、複数の解析結果それぞれを時間的に対応して得られる情報として、カラーバーＣ１，Ｃ２を近接させて表示することも可能である。これにより、左右肺の呼吸位相をそれぞれ同時に示すことができ、左右肺の状態変化が異なる箇所が明確に示され、異常点がわかりやすくなるため、動画のフレーム選択操作が容易になる。
【０１１８】
このように、診断支援情報は、複数の部位に対応する複数の解析結果の時間変化を示す情報であるため、表示用画像ＩＧを参照することにより、複数の部位に対応する複数の解析結果を同時に視認することが可能となる。また、複数の部位が、左肺野及び右肺野であることにより、表示用画像ＩＧを参照して、左肺野及び右肺野それぞれの解析結果を同時に視認することが可能となる。
【０１１９】
また、再生時刻表示部分１６３（すなわち、再生時刻表示部ＴＭ）は、一覧表示部分１６２（すなわち、カラーバーＣ１，Ｃ２）に近接させて認識可能にした表示形態を採用することが望ましい。これによって、再生時における診断支援情報と再生時間とを視覚的に同時に把握することが可能となる。
【０１２０】
＜４．第４実施形態＞
ユーザは動画に注視したいため、フレーム選択操作には重要な点のみを表示した簡易表示が効果的である。そこで、第４実施形態では、設定条件下で定められた特徴点を算出し、表示用画像ＩＧに付加される。なお、ここでいう特徴点についての詳細は後述するが、上述の第２の呼吸情報検出方法及び後述の第２の血流情報検出方法における特徴点とは異なることに留意したい。
【０１２１】
図１９は、本発明の第４実施形態として構成された画像生成装置３Ａで用いられる制御部３１Ａの機能構成を示す図である。この制御部３１Ａは、第１実施形態の画像生成装置３における制御部３１（図４参照）の代替としてそれぞれ使用される。第１実施形態と異なる点は、検出部１２０Ａが特徴点算出部１４０をさらに備える点であり、残余の構成は第１実施形態の画像生成装置３と同様である。
【０１２２】
＜４−１．特徴点算出部１４０＞
特徴点算出部１４０は、対象部位の物理的状態の時間変化において、特徴点を算出する。なお、診断支援情報生成部１５０Ａ及び表示画像生成部１６０Ａは、当該特徴点を指示する情報を含んで診断支援情報を生成するものとする。
【０１２３】
図２０は、検出部１２０Ａにおいて検出された血流（心拍）情報の波形データを時系列で示した模式図であり、図２１は、検出部１２０Ａにおいて検出された呼吸情報の波形データを時系列で示した模式図である。なお、図２０は、第１の血流情報検出方法を採用する場合における心電図波形を時間方向にモニタリングした結果であり、図２１は、第２の呼吸情報検出方法を採用する場合における肺野領域の面積値あるいは特徴点間距離の時間的変化を算出した結果を例示したものである。
【０１２４】
設定条件下で定められた特徴点としては、例えば、血流（心拍）情報に関しては、図２０で示されるように、波形中の最大点として点Ｐ１及び点Ｐ４（図７における点Ｒｐに相当する点）、波形中の最小点として点Ｐ２及び点Ｐ５（図７における点Ｓｐに相当する点）、及び、波形中の極大点として点Ｐ３及び点Ｐ６（図７における点Ｔｐに相当する点）等が挙げられる。また、呼吸情報に関しても同様に、図２１で示されるように、最大点の点Ｂ１，点Ｂ３、及び、最小点の点Ｂ２等が挙げられる。
【０１２５】
また、特徴点算出部１４０では、最大点、最小点、極大点、及び、極小点の他にも、例えば、変化点（例えば、１次微分または２次微分における、最大点、最小点、極大点、極小点）等を算出するよう設定条件として設けても良い。
【０１２６】
特徴点算出部１４０により算出された特徴点は、上述のカラーバーＣ１（Ｃ２）上に重畳して示されるように、診断支援情報生成部１５０が表示用画像ＩＧを生成する。すなわち、血流（心拍）情報に関しては、図２０で示されるように、最大点の点Ｐ１及び点Ｐ４はラインＬＰ１，ＬＰ４として、最小点の点Ｐ２及び点Ｐ５はラインＬＰ２，ＬＰ５として、極大点の点Ｐ３及び点Ｐ６はラインＬＰ３，ＬＰ６としてそれぞれ表示される。また、呼吸情報に関しても同様に、図２１で示されるように、最大点の点Ｂ１及び点Ｂ３はラインＬＢ１，ＬＢ３として、最小点の点Ｂ２はラインＬＢ２としてそれぞれ表示される。なお、図２０ではＬＰ１〜ＬＰ６を明確にするため、カラーバーＣ１（Ｃ２）が空白で示されているが、実際は、心臓の「収縮期」「拡張期」（対象部位の「状態」変化）を視覚的に識別可能に表示されているものとする。
【０１２７】
また、当該特徴点を示すラインＬＰ１〜ＬＰ６，ＬＢ１〜ＬＢ３を識別可能に色表示する等、明確にすることにより、診断重要箇所が視認しやすくなり、診断効率がより向上する。また、診断支援情報が複数の部位に対応する複数の解析結果の時間変化を示す場合、例えば、図１７（ｂ）で示されるように、左肺野および右肺野の呼吸位相をそれぞれ示す場合に、上述の特徴点を示すラインを重畳して示されれば、左肺野および右肺野の状態の時間変化に応じた異常点（特徴点）が異なる箇所で現れるなど、状態変化の差異が明確になり有用である。
【０１２８】
また、特徴点が示されることにより、状態変化に異常の疑いが見られ得ると診断されたフレーム画像ＭＩに対しては、色表示等識別可能にカラーバーＣ１（Ｃ２）上に表示することも可能である。さらに、呼気か吸気か等の第１及び第２の解析結果では信頼性が低い状況においても、特徴点を色表示等識別可能にカラーバーＣ１（Ｃ２）上に表示することで、異常が発見しやすくなる。
【０１２９】
＜４−２．画像生成装置３Ａの基本動作＞
続いて、図２２は、第４実施形態に係る画像生成装置３Ａの動作フローを例示した図である。なお、図２２のうち、ステップＳＴ１，ＳＴ２，ＳＴ４，ＳＴ５，ＳＴ６は図１８のステップＳ１〜Ｓ５と同様であるため、その説明は省略する。
【０１３０】
この第４実施形態では、第１実施形態では存在しなかった特徴点算出部１４０が付加されたことで、下記の工程が加わる。
【０１３１】
すなわち、第１実施形態と同様の工程として、ステップＳＴ１〜ステップＳＴ２を経て、図２２で示されるように、ステップＳＴ３にて、検出部１２０Ａにおける特徴点算出部１４０が、ステップＳＴ２において検出された対象部位の時間変化において、設定条件下で定められた特徴点を算出する（図２０及び図２１参照）。
【０１３２】
そして、第１実施形態と同様の工程として、ステップＳＴ４〜ステップＳＴ５を経て、最後に、ステップＳＴ６において、表示画像生成部１６０Ａが、ステップＳＴ５において生成された、当該特徴点を指示する情報を含む表示用画像ＩＧを表示部３４にて出力することで、表示部３４のモニタに表示し、本動作フローが終了される。
【０１３３】
以上のように画像生成装置３Ａでは、診断支援情報は特徴点を指示する情報を含むため、対象部位の時間変化における特徴点が明確になり、診断効率がより向上する。
【０１３４】
＜５．第５実施形態＞
ユーザの所望の条件を満たすフレーム画像ＭＩの再生時刻が到来する、もしくは、到来したタイミングを、ユーザに報知することで、特に経験の浅いユーザには効果的となる。そこで、第５実施形態では、所望の条件を満たすタイミングをユーザに報知する手段を備える。
【０１３５】
図２３は、本発明の第５実施形態として構成された画像生成装置３Ｂで用いられる制御部３１Ｂの機能構成を示す図である。この制御部３１Ｂ及び表示部３４Ｂは、第１（第４）実施形態の画像生成装置３（３Ａ）における制御部３１（３１Ａ）（図４、図１９参照）の代替としてそれぞれ使用される。第１（第４）実施形態と異なる点は、検出部１２０Ｂが報知点算出部１４５をさらに備え、また、画像生成装置３Ｂが報知部３４２をさらに備える点である。なお、残余の構成は第４実施形態の画像生成装置３Ａと同様である。
【０１３６】
なお、図２３における検出部１２０Ｂは、第４実施形態による特徴点算出部１４０を備えて構成されているが、特徴点算出部１４０を備えていなくても良い。
【０１３７】
＜５−１．報知点算出部１４５及び報知部３４２＞
報知点算出部１４５は、対象部位の時間変化において、ユーザが所望する設定条件下で定められた報知するための点（以下「報知点」と称する）を算出し、報知部３４２に出力する。すなわち、設定条件とは、ユーザが指定する条件であり、例えば、対象部位の物理的状態の時間変化が図２１で示される呼吸情報であるときに、ユーザが最大点を報知するよう指定する場合は、報知点算出部１４５が点Ｐ３及び点Ｐ６を報知点として検出する。なお、診断支援情報生成部１５０Ｂ及び表示画像生成部１６０Ｂは、当該報知点を指示する情報を含んで診断支援情報を生成するものとする。
【０１３８】
報知部３４２は、診断支援情報生成部１５０による解析結果が設定条件（所定の条件）を満たしたときに、ユーザに設定条件を満たしたことを報知する。すなわち、報知部３４２は、報知点算出部１４５において検出された報知点を、ユーザに視覚情報、聴覚情報、触覚情報の何れかの手段により報知する。視覚情報として報知するときは、報知部３４２は表示部３４Ｂに視覚情報を表示するよう指示する。具体的に視覚情報としては、現時点から報知点までにかかる時間を視覚的に表現したものであり、インジゲータ、プログレスバー表示、数値による表示、モデル図による表示、周期図による表示などが挙げられ、好ましくは、ユーザが余裕を持って報知点を知得できるように、報知点が実際に到来するよりも前から、その報知点が近づいていることを知ることができるような態様で画面上に表示する。聴覚情報として報知するときは、報知部３４２はブザーやタイミング音、音声などが挙げられる。例えば、報知点まで後何秒かを合成音でアナウンスする方法や、報知点でブザーを鳴らす方法等により報知する。
【０１３９】
また、視覚情報、聴覚情報、及び、触覚情報のいずれを採用した場合においても、プログレスバーＰＢをまったく見ずに操作が可能となる。例えば、経過時間に応じてフレーム画像ＭＩを巻き戻し操作中に、報知情報をユーザが受けるため、動画を注視しながら診断に有用なフレーム画像ＭＩまで到達選択することが可能となる。
【０１４０】
＜５−２．画像生成装置３Ｂの基本動作＞
続いて、図２４は、第５実施形態に係る画像生成装置３Ｂの動作フローを例示した図である。なお、検出部１２０Ｂは、特徴点算出部１４０を備えていない場合を想定するものとする。また、図２４のうち、ステップＳＰ１，ＳＰ２，ＳＰ４，ＳＰ５，ＳＰ６は図１８のステップＳ１〜Ｓ５と同様であるため、その説明は省略する。
【０１４１】
この第５実施形態では、第１実施形態では存在しなかった報知点算出部１４５及び報知部３４２が付加されたことで、下記の工程が加わる。
【０１４２】
すなわち、第１実施形態と同様の工程として、ステップＳＰ１〜ステップＳＰ２を経て、図２４で示されるように、ステップＳＰ３にて、検出部１２０Ａにおける報知点算出部１４５が、ステップＳＰ２において検出された対象部位の時間変化において、設定条件下で定められた報知点を算出する。
【０１４３】
そして、第１実施形態と同様の工程として、ステップＳＰ４〜ステップＳＰ５を経て、最後に、ステップＳＰ６において、表示画像生成部１６０Ｂが、ステップＳＰ５において生成された、当該報知点のタイミングを考慮した表示用画像ＩＧを表示部３４Ｂにて出力することで、表示部３４Ｂのモニタに表示し（ここで、聴覚情報及び触覚情報にて当該報知点のタイミングをユーザに知らせる場合は音声や触覚などにて出力し）、本動作フローが終了される。
【０１４４】
以上のように画像生成装置３Ｂでは、対象部位の解析結果が所望の設定条件を満たしたときに、ユーザに報知することにより、診断経験の浅い医者等に対しても、報知により設定条件を満たした診断内容を認識することが可能となる。
【０１４５】
＜６．第６実施形態＞
図６は、本発明の第６実施形態として構成された放射線動態画像撮影システム１００’における画像生成装置３’で用いられる制御部３１’の機能構成を示す図である。この制御部３１’は、第１実施形態のシステム１００における制御部３１（図４）の代替として使用される。第１実施形態と異なる点は、検出部１２０’における血流情報の検出方法が異なる点である。なお、残余の構成は第１実施形態の画像生成装置３と同様である。また、本実施形態における対象部位とは、心臓領域または肺領域である。以下では、本実施形態で用いられる第２の血流情報検出方法について説明する。
【０１４６】
○  第２の血流情報検出方法：心臓壁の動き量
第２の血流情報検出方法として、図６で示されるように検出部１２０’（所定部位周期特定部１３０’）では、動態画像取得部１１０によって取得された撮影画像を用いて、心臓壁の動き量を算出することで、心拍（血流）情報とするものである。すなわち、呼吸中肺動態画像及び呼吸停止中肺動態画像で、撮影対象の対象部位である肺とともに心臓も捉えられていることが前提条件である。詳細には、呼吸中肺動態画像及び呼吸停止中肺動態画像から心臓壁の変動が検出されることで、各呼吸中フレーム画像及び各呼吸停止中フレーム画像が撮影されたタイミングにおける心臓の拍動の位相が検出される。したがって、心臓壁が心臓の拍動の位相として検出される。
【０１４７】
図８は、各呼吸中フレーム画像及び各呼吸停止中フレーム画像について、所定点（例えば左上の点）を基準点（例えば原点）とし、右方向をＸ軸方向とし、下方向をＹ軸方向とし、１画素毎に座標の値が１つずつ変化するように描かれた画像座標面である。
【０１４８】
図９は、呼吸停止中肺動態画像で捉えられた心臓壁の変動を例示する模式図である。心臓壁の変動の一例として、心臓の横幅の変動を採用する。図９(ａ)〜(ｃ)では、心臓が拡張していく過程で、心臓壁の動きである心臓の横幅がｗ１からｗ３へと大きくなっていく状態が例示されている。
【０１４９】
そこで、所定部位周期特定部１３０では、各呼吸中フレーム画像及び各呼吸停止中フレーム画像から、心臓の横幅を検出することで、心拍（血流）周期を設定する。具体的に、心臓の横幅を検出する手法としては、例えば、心臓の輪郭を検出して行う手法等が挙げられる。そして、この心臓の輪郭を検出する手法としては、種々の公知の手法を採用することができ、例えば、心臓の形状を示すモデル（心臓モデル）を用いて、Ｘ線画像中の特徴点と、心臓モデルの特徴点とを合わせて行くことで、心臓の輪郭を検出する手法（例えば、"Image feature analysis and computer-aided diagnosis in digital radiography: Automated analysis of sizes of heart and lung in chest images", Nobuyuki Nakamori et al., Medical Physics, Volume 17, Issue 3, May, 1990, pp.342-350.等参照）等を採用することができる。
【０１５０】
図１０は、呼吸停止中肺動態画像を構成する複数の呼吸停止中フレーム画像について、撮影された時刻と心臓の横幅との関係を例示する模式図である。図１０では、横軸が時刻、縦軸が心臓の横幅を示し、丸印が検出された心臓の横幅の値を示している。
【０１５１】
ここで、時刻ｔで捉えられた心臓の横幅をＨｗt、時刻ｔ＋１で捉えられた心臓の横幅をＨｗt+1とし、（Ｈｗt+1−Ｈｗt）≧０が成立する場合には、時刻ｔで捉えられた呼吸停止中フレーム画像が心臓の拡張時に分類され、（Ｈｗt+1−Ｈｗt）＜０が成立する場合には、時刻ｔで捉えられた呼吸停止中フレーム画像が心臓の収縮時に分類される。
【０１５２】
以上のように、所定部位周期特定部１３０’では、動態画像で捉えられた心臓壁の動き（所定部位の形状の変化）に基づき、心拍（血流）周期を検出するため、当該周期を自動的に取得することが可能となる。
【０１５３】
また、第２の血流情報検出方法は、第１の血流情報検出方法と比較して間接的に血流周期を検出するため、ノイズ成分も含まれやすいと予想される。そこで、所定部位周期特定部１３０’では、動態画像で捉えられた心臓壁の動き（所定部位の形状の変化）に基づき、周波数解析等を用いて血流周期を検出することが好ましい。これにより、ノイズ成分を取り除いた所望の変動成分を自動的に抽出することができるため、心臓壁の動き量（所定部位が時間変化する状態）をより正確に把握することが可能となる。
【０１５４】
＜７．第７実施形態＞
以下では、第７実施形態について説明する。第１実施形態と異なる点は、検出部１２０における血流情報の検出方法が異なる点である。なお、残余の構成は第１実施形態の画像生成装置３と同様である。また、第７実施形態における血流情報の検出方法は第６実施形態における血流情報の検出方法とは異なるが、図６で示されるように、動態画像に基づき血流情報（対象部位の物理的状態の時間変化）を検出する点においては共通する。以下では、本実施形態で用いられる第３の血流情報検出方法について説明する。
【０１５５】
○  第３の血流情報検出方法：血流位相解析 
第３の血流情報検出方法として、動態画像取得部１１０によって取得された撮影画像を用いて、血流位相解析を行うことで、血流情報とするものである。血流位相とは、血流が流れている位置に応じた血流の有無を示す位相情報である。本発明で用いる血流位相解析処理（血流情報生成処理）に関しては、例えば、本出願人による出願である「特願２０１１−１１５６０１号（出願日：平成２３年５月２４日）」を採用することができる。
【０１５６】
図１１は、肺全体の血流の有無に伴う時空間変化の解析結果を例示する図である。一般的に、血流解析では、心臓の収縮によって右心室から大動脈を介して血液が急激に吐出されることにより肺野血管が拡がるので、この拡がりを、動態画像を解析することにより抽出して肺全体の血流の有無に関する診断支援情報として出力するものである。即ち、肺野において血管が拡張すると、肺血管が拡がった領域の放射線透過量が、肺野（肺胞）領域を透過する放射線透過量よりも比較的大きく減少するので、この領域に対応する放射線検出部１３の出力信号値が低下する。このような心臓の拍動に呼応した肺血管の拡がりは、心臓近傍の動脈から末梢へ伝播する。そこで、動態画像を構成する一連のフレーム画像ＭＩ間の放射線検出部１３の画素（ピクセル）単位、または、複数画素からなる小領域単位（画素ブロック単位）を互いに対応付け、画素単位又は小領域単位毎に、信号値が最も低くなったフレーム画像ＭＩを求め、そのフレーム画像ＭＩの該当する領域を肺血管が血流により拡張したタイミングを示す信号として色付けする。そして、図１１で示されるように、色付け後の一連のフレーム画像ＭＩを表示部３４に順次表示することで、血流の状態を医師等が視認できるようになる。なお、図１１で示される白色箇所は、実際は赤色等で色付けされている。
【０１５７】
各画素（小領域）において、血流により肺血管が拡張したタイミングを示す信号（血流信号という）は、その画素（小領域）の信号値の時間変化を示す波形（出力信号波形という）の極小値を求めることにより取得することができる。この血流信号は心臓の拍動周期と同じ間隔で表れるが、不整脈等の異常個所があると、心臓の拍動周期と異なる間隔で血流に伴う血管の拡張と関係なく極小値が表れる場合がある。そこで、第３の血流情報検出方法では、心臓の拍動を示す拍動信号波形と各小領域の出力信号波形との相関係数を求めることにより、精度良く血流信号を抽出する。
【０１５８】
以下では、肺血管領域における血流情報の検出として、肺血管領域の出力信号波形と血流信号抽出方法との概要を説明する。図１２は、肺血管領域の血流信号値の時間変化を示す波形を例示する図である。図１２（ａ）では、時間順次に取得された一連のフレーム画像ＭＩにおいて、診断対象領域に相当する肺血管領域ＩＲ２の位置を示し、図１２（ｂ）では、横軸を動態画像の撮影開始からの経過時間（フレーム番号）とし、縦軸を肺血管領域ＩＲ２における信号値（代表値）とした座標空間上に各フレーム画像ＭＩの肺血管領域ＩＲ２の信号値（代表値）がプロットされたグラフを示す。
【０１５９】
ここで、図１２（ｂ）のグラフでは、呼吸と血流との位相が混在した状態であるため、次のフィルタリング処理によって、呼吸と血流のそれぞれ互いの影響を排除する。すなわち、当該フィルタリング処理では、呼吸等による低周波数な信号変化を除去し、血流による信号値の時間変化を抽出する。例えば、小領域毎の信号値の時間変化に対して、安静呼吸画像群では低域カットオフ周波数０．７Ｈｚ、深呼吸画像群では低域カットオフ周波数０．５Ｈｚでハイパスフィルタリングを行う。若しくは、さらに高周波数のノイズ成分を除去するために２．５Ｈｚの高域カットオフ周波数で高周波数も遮断するバンドパスフィルタによってフィルタリングを行っても良い。
【０１６０】
ここで、当該カットオフ周波数は、固定値とするよりも、撮影された動態画像毎に最適化することが好ましい。例えば、一連のフレーム画像ＭＩ（図１２（ａ）参照）の心臓領域の信号変化から、心臓の収縮期の時間と拡張期（弛緩期）の時間を算出する（第２の血流情報検出方法を参照）。そして、拡張期の時間の逆数に所定の係数を乗じた値をハイパスフィルタもしくはバンドパスフィルタで低周波数を遮断するカットオフ周波数として設定し、また、バンドパスフィルタの場合は、収縮期の時間の逆数に所定の係数を乗じた値を、高周波数を遮断する高域カットオフ周波数として設定する。更に、低域カットオフ周波数は、呼吸による周波数成分を考慮し、一連のフレーム画像ＭＩから横隔膜の位置の値や肺野領域の面積値あるいは特徴点間距離を解析し（後述の第２の呼吸情報検出方法を参照）、安静換気の場合は安静呼気位及び安静吸気位となるフレーム画像ＭＩを検出し、安静呼気位のフレームとその次の安静吸気位のフレームの間のフレーム数から吸気期の時間を求め、その逆数と、当該拡張期の時間の平均値に所定の係数を乗じた値を低域のカットオフ周波数として設定してもよい。このとき安静換気の場合、自動的に設定されるカットオフ周波数は、低域カットオフ周波数は０．２〜１．０Ｈｚの間に、高域カットオフ周波数は、２．０Ｈｚ以上に制限することが好ましい。また、別途測定された安静時における１分間の呼吸数、脈拍数等（第１の血流情報検出方法、及び、後述の第１の呼吸情報検出方法を参照）のバイタルサインを患者情報として入力しておき、これらの値からカットオフ周波数が算出されるようにしてもよい。例えば、患者情報として入力された１分間の呼吸数を１秒間の呼吸数に変換し、その呼吸数に所定の係数を乗じた値を低域カットオフ周波数としてもよい。そして、入力された１分間の脈拍数を１秒間の脈拍数に変換し、１秒間の呼吸数に所定の係数を乗じた値を高域カットオフ周波数としてもよい。更に、１秒間の呼吸数と１秒間の心拍数の平均値に所定の係数を乗じた値を低域カットオフ周波数として設定するようにしてもよい。
【０１６１】
以上の血流位相解析処理の後、本実施形態の所定部位周期特定部では、動態画像で捉えられた血流位相変化（所定部位の状態の変化）に基づき、血流周期を特定するため、血流周期を自動的に取得することが可能となる。
【０１６２】
＜８．変形例＞
以上、本発明の実施形態について説明してきたが、本発明は、上記実施形態に限定されるものではなく、様々な変形が可能である。
【０１６３】
※  本実施形態では、画像生成装置３，３’，３Ａ，３Ｂを個別に実施されるように各実施形態に分けて記載したが、これらの個別機能は、互いに矛盾しない限り、相互に組み合わせてもよい。
【０１６４】
※  この発明では、身体の撮影対象となる部分のうち、物理的状態が周期的に時間変化する部位が位相検出の対象となるが、それは、心臓や肺だけでなく、蠕動などの不随意運動を行う他の臓器であってもよく、また、筋肉や関節などの随意運動を行う部位であってもよい。後者の場合には、被検者に同一の動作を繰り返して行わせつつ、動態撮影を行う。
【０１６５】
すなわち、本実施形態では、胸部撮影における呼吸情報や血流情報を対象としたが、例えば、関節撮影における関節の屈伸方向情報などを対象としてもよい。
【０１６６】
図２５は、関節屈伸における関節角度を例示する図である。図２５で示されるように、検出部１２０（１２０Ａ，１２０Ｂ）によって、屈伸方向情報は関節角度θの動きから算出される。具体的に、第２の呼吸情報検出方法の輪郭抽出方法や閾値処理などを採用し、骨の輪郭領域を抽出する。そして、その輪郭領域から軸ＡＸ１，ＡＸ２を抽出し、軸ＡＸ１，ＡＸ２の２本が交差する角度θの変動から屈伸方向を算出する。軸ＡＸ１，ＡＸ２は、輪郭領域の中心を通るものであっても、輪郭領域のエッジに沿ったラインであってもよい。
【０１６７】
図２６は、所定部位周期特定部１３０において特定された屈伸方向情報の周期に対して、屈伸位相の診断支援情報の生成を模式的に例示する図である。図２６で示されるように、ここで、第１及び第２の解析結果とは、伸ばす方向と屈する方向とを指し、屈伸位相を時間変化に対応させ、屈伸位相周期においても、図１１，図１５と同様に識別可能な診断支援情報を生成する。
【０１６８】
これにより、表示用画像ＩＧにおける診断支援情報を参照して、屈伸周期における伸ばす方向と屈する方向との変化が識別可能となるため、屈伸周期の関節の状態変化に対して、効率的な医療診断が可能となる。
【０１６９】
※  検出部１２０、特徴点算出部１４０及び報知点算出部１４５の画像解析において、フレーム画像ＭＩの領域うち、画像解析の対象領域を適宜設定することが可能である。これにより、画像解析に掛かる計算時間の短縮化が図れる。
【０１７０】
※  本実施形態では、図１６及び図１７で示されるように、診断支援情報の解析結果として、表示用画像ＩＧにおける一覧表示部分１６２（カラーバーＣ１，Ｃ２）に呼吸情報を表示する場合を示しているが、図１１で示される血流情報を表示しても良い。また、表示画像生成部１６０は、カラーバーＣ１（Ｃ２）上に血流位相を表示するだけでなく、上述の第３の血流情報検出方法（血流位相解析）による画像処理結果ＲＴ（図１１参照）をフレーム画像ＭＩに重畳もしくは近接して表示するように表示用画像ＩＧを生成してもよい。
【０１７１】
図２７は、診断支援情報の解析結果として血流情報を表示する場合の表示用画像ＩＧを例示している。当該表示用画像ＩＧを生成するにあたって、図４で示されるように、検出部１２０は、フレーム画像ＭＩ（動態画像）に基づき、対象部位の血流の有無に関する画像処理結果ＲＴ（図１１参照）を出力し、動態画像と時間的に関連付けて保持部３２に保持する。そして、表示画像生成部１６０は、画像処理結果ＲＴを動態画像と同期して表示するように表示用画像ＩＧを生成する。
【０１７２】
なお、表示用画像ＩＧには、画像処理結果ＲＴを視覚的に認識可能な表示内容として表示する画像処理結果表示部分１６４（図４における破線エリア）が設けられることが望ましい。すなわち、図２７で示されるように、画像処理結果表示部分１６４（すなわち、画像処理結果ＲＴの表示部分）は、動態画像表示部分１６１（すなわち、フレーム画像ＭＩの表示部分）に重畳（不図示）もしくは近接（図２７参照）させて認識可能にした表示形態を採用することが望ましい。これによって、フレーム画像ＭＩと診断支援情報と画像処理結果ＲＴとを視覚的に同時に把握することが可能となり、より効率的な医療診断が可能となる。
【０１７３】
これにより、肺塞栓が疑われる場合には、その疑い箇所を詳細に確認して診断することが可能となる。例えば、肺野内の血流位相状況をカラーバーＣ１上で判断し、画像処理結果ＲＴ（重畳表示した場合はフレーム画像ＭＩ）上で肺塞栓の可能性がある肺血管の血流状況に注視することができる。また、肺塞栓の可能性のある位置が心臓付近の肺血管なのか抹消の肺血管なのか見当がついていれば、ユーザがプログレスバーＰＢを用いて位相情報から診断で見たい再生時刻を調整することができる。
【０１７４】
また、診断支援情報が複数の部位に対応する複数の解析結果の時間変化を示す場合の例として、図１７（ｂ）では右肺野及び左肺野の位相をそれぞれカラーバーＣ１，Ｃ２に示したが、呼吸情報と血流情報との位相を同時に表示してもよい。これは、血流の診断は、呼吸位相が最大吸気の状態が、肺が最も大きくなるために最も視認しやすくなるためである。
【０１７５】
また、上述したように、検出部１２０にて対象領域のみの時間変化を検出し、その注目エリアの位相情報をカラーバーＣ１（Ｃ２）上に表示するようにしてもよい。これは、肺塞栓の疑いがある場所が明確な場合、当該エリアに限定した位相情報の表示が望ましいためである。
【０１７６】
さらに、血流の位相は、指定領域の血流有無による２位相だけでなく、肺野の主血管領域と抹消血管領域とにエリアを分けて、どちらに血流があるかによって３位相で表現してもよい。それにより、主血管と抹消血管とどちらに肺塞栓疑いがあるかに応じた再生時刻の指定選択が容易となる。
【０１７７】
※  表示画像生成部１６０，１６０Ａ，１６０Ｂが生成する表示用画像ＩＧは、本実施形態の例に限られない。すなわち、ユーザーカスタマイズを可能にすることで、あらゆる診断観点に対応して表示用画像ＩＧを生成することができる。例えば、ユーザが、グラフ表示された特定の動き（状態変化）情報をクリックしてフレーム選択操作していた場合、その動き情報をターゲットとして、色表示を切り替え変更できるようにしてもよい。また、ある画素領域を指定することで、その画素領域の時間方向について動きを解析し、分析結果で色表示を切替てもよい。
【０１７８】
※  本実施形態では、プログレスバーＰＢが横方向に動くことでフレーム画像ＭＩの位置を表現する場合は、カラーバーＣ１（Ｃ２）も、プログレスバーＰＢに沿った横方向座標に対応させて色を変化させたが（図１６及び図１７参照）、これに限られない。例えば、プログレスバーＰＢが縦方向や回転方向に動くことでフレーム画像ＭＩの位置を表現する場合には、カラーバーＣ１（Ｃ２）も、その方向に沿って色を変化させるものであってもよい。
【０１７９】
※  被写体は、人体だけでなく、動物の身体であってもよい。
【符号の説明】
【０１８０】
１  撮影装置
２  撮影制御装置
３，３’，３Ａ，３Ｂ  画像生成装置
４  心電計
３１，３１’，３１Ａ，３１Ｂ  制御部
３４，３４Ｂ  表示部
４１  位相検出部
１００，１００’，１００Ａ，１００Ｂ  放射線動態画像撮影システム
１１０  動態画像取得部
１２０  検出部
１３０，１３０’  所定部位周期特定部
１４０  特徴点算出部
１４５  報知点算出部
１５０，１５０Ａ，１５０Ｂ  診断支援情報生成部
１６０，１６０Ａ，１６０Ｂ  表示画像生成部
１６１  動態画像表示部分
１６２  一覧表示部分
１６３  再生時刻表示部分
１６４  画像処理結果表示部分
３４１  再生時刻調整部
３４２  報知部
Ｃ１，Ｃ２  カラーバー
ＰＢ  プログレスバー
Ｍ  被写体（被検者）
ＭＩ  フレーム画像
ＴＭ  再生時刻表示部
ＲＴ  画像処理結果
【要約】
画像生成装置３において、検出部１２０は、対象部位を時間順次に捉えた動態画像（フレーム画像）に対して、対象部位の物理的状態の時間変化状態を検出し、診断支援情報生成部１５０は、対象部位の物理的状態の時間変化に基づいて解析を行い、その解析結果を診断支援情報として生成する。当該解析結果は、例えば、呼気と吸気とからなる２つの診断情報である。表示画像生成部１６０は、動態画像及び前記診断支援情報を表示するための表示用画像を生成する。表示用画像は、動態画像を表示する動態画像表示部分１６１と、診断支援情報の第１の解析結果と第２の解析結果とを識別可能に時間軸方向に一覧できるように表示する一覧表示部分１６２（カラーバー等）とを含む画像である。これにより、動画を再生する際に、作業性や視認性を向上することができ、診断効率を上げることが可能となる。【図１】
0005408399.tif 000002
【図２】
0005408399.tif 000003
【図３】
0005408399.tif 000004
【図４】
0005408399.tif 000005
【図５】
0005408399.tif 000006
【図６】
0005408399.tif 000007
【図７】
0005408399.tif 000008
【図８】
0005408399.tif 000009
【図９】
0005408399.tif 000010
【図１０】
0005408399.tif 000011
【図１２】
0005408399.tif 000012
【図１３】
0005408399.tif 000013
【図１４】
0005408399.tif 000014
【図１５】
0005408399.tif 000015
【図１６】
0005408399.tif 000016
【図１７】
0005408399.tif 000017
【図１８】
0005408399.tif 000018
【図１９】
0005408399.tif 000019
【図２０】
0005408399.tif 000020
【図２１】
0005408399.tif 000021
【図２２】
0005408399.tif 000022
【図２３】
0005408399.tif 000023
【図２４】
0005408399.tif 000024
【図２５】
0005408399.tif 000025
【図２６】
0005408399.tif 000026
【図２８】
0005408399.tif 000027
【図１１】
0005408399.tif 000028
【図２７】
0005408399.tif 000029