(57)【要約】【課題】動態画像の解析結果の画像の視認性を向上させることができる動態画像解析装置及びプログラムを提供する。【解決手段】診断用コンソール３によれば、制御部３１は、放射線撮影により得られた動態画像から所定の特徴量Ｆを解析し当該所定の特徴量Ｆを表示部３４に表示する。また、制御部３１は、前記動態画像の信頼性を評価する。次いで、制御部３１は、所定の特徴量Ｆが表示部３４に表示される際に、当該所定の特徴量Ｆと対応付けて動態画像の信頼性を示す信頼性情報Ｄ２を表示する。【選択図】図３
2020000807.tif 000002
【特許請求の範囲】
【請求項１】
放射線撮影により得られた動態画像から所定の特徴量を解析し当該所定の特徴量を表示手段に表示する動態画像解析装置であって、
前記動態画像の信頼性を評価する評価手段と、
前記所定の特徴量が前記表示手段に表示される際に、当該所定の特徴量と対応付けて前記評価手段によって評価された前記動態画像の信頼性を示す信頼性情報を表示する表示制御手段と、
を備えることを特徴とする動態画像解析装置。
【請求項２】
前記動態画像を構成する一連のフレーム画像から被験者の体動を検出する検出手段と、
前記検出手段によって検出された前記被験者の体動の量を算出する第１の算出手段と、を備え、
前記評価手段は、前記第１の算出手段によって算出された前記被験者の体動の量に基づき、前記動態画像の信頼性を評価することを特徴とする請求項１に記載の動態画像解析装置。
【請求項３】
前記評価手段は、前記第１の算出手段によって算出された前記被験者の体動の量が当該体動の量に関する閾値を超えている場合、前記動態画像の信頼性が低い旨の評価を行い、
前記表示制御手段は、前記評価手段によって前記動態画像の信頼性が低い旨の評価が行われた場合、当該動態画像の信頼性が低い旨を示す前記信頼性情報を前記表示手段に表示するとともに、前記被験者の体動の量が前記閾値を超えたときの当該体動の発生箇所を前記表示手段に表示することを特徴とする請求項２に記載の動態画像解析装置。
【請求項４】
前記表示制御手段は、前記評価手段によって前記動態画像の信頼性が低い旨の評価が行われた場合、更に、前記閾値を超えた体動の量を前記表示手段に表示することを特徴とする請求項３に記載の動態画像解析装置。
【請求項５】
前記動態画像の画素信号ごとに当該画素信号のＳＮ比を算出する第２の算出手段を備え、
前記評価手段は、前記第２の算出手段によって算出された前記ＳＮ比に基づき、前記動態画像の信頼性を評価することを特徴とする請求項１〜４のいずれか一項に記載の動態画像解析装置。
【請求項６】
前記評価手段は、前記第２の算出手段によって算出された前記ＳＮ比が当該ＳＮ比に関する閾値未満の場合、前記動態画像の信頼性が低い旨の評価を行い、
前記表示制御手段は、前記評価手段によって前記動態画像の信頼性が低い旨の評価が行われた場合、当該動態画像の信頼性が低い旨を示す前記信頼性情報を前記表示手段に表示するとともに、前記ＳＮ比が前記閾値未満である前記画素信号の画素位置を前記表示手段に表示することを特徴とする請求項５に記載の動態画像解析装置。
【請求項７】
前記動態画像の画素信号ごとに当該画素信号から抽出される前記所定の特徴量の解析対象となる周波数成分の信号と所定の理想の信号との相関係数を算出する第３の算出手段を備え、
前記評価手段は、前記第３の算出手段によって算出された前記相関係数に基づき、前記動態画像の信頼性を評価することを特徴とする請求項１〜６のいずれか一項に記載の動態画像解析装置。
【請求項８】
前記評価手段は、前記第３の算出手段によって算出された前記相関係数が当該相関係数に関する閾値未満の場合、前記動態画像の信頼性が低い旨の評価を行い、
前記表示制御手段は、前記評価手段によって前記動態画像の信頼性が低い旨の評価が行われた場合、当該動態画像の信頼性が低い旨を示す前記信頼性情報を前記表示手段に表示するとともに、前記相関係数が前記閾値未満である前記画素信号の画素位置を前記表示手段に表示することを特徴とする請求項７に記載の動態画像解析装置。
【請求項９】
放射線撮影により得られた動態画像から所定の特徴量を解析し当該所定の特徴量を表示手段に表示する動態画像解析装置のコンピュータを、
前記動態画像の信頼性を評価する評価手段、
前記所定の特徴量が前記表示手段に表示される際に、当該所定の特徴量と対応付けて前記評価手段によって評価された前記動態画像の信頼性を示す信頼性情報を表示する表示制御手段、
として機能させるためのプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、動態画像解析装置及びプログラムに関する。
【背景技術】
【０００２】
従来、被写体の周期性を持つ動態を放射線撮影することにより得られた動態画像が診断に用いられている。動態画像では、静止画では捉えられなかった被写体の動態を表示・解析することができる。
【０００３】
例えば、特許文献１には、時刻情報を有する複数のボリュームデータ（動態画像）を生成し、それぞれの画質情報を生成するとともに、当該複数のボリュームデータのそれぞれの動き解析情報を生成し、当該動き解析情報及び当該画質情報に基づき、複数のボリュームデータからそれぞれ画像（解析結果の画像）を生成するときの画質を向上させることができる医療用画像処理装置が記載されている。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特開２０１４−５０７５３号公報
【発明の概要】
【発明が解決しようとする課題】
【０００５】
上記特許文献１に記載されている医療用画像処理装置では、複数のボリュームデータから生成された画像（解析結果の画像）は、当該ボリュームデータのボクセルの信頼性の高さに応じて変わるパラメータのフィルタを用いて補正されるようになっているが、当該画像を見ただけでは当該画像のどの領域がどれだけ補正されたのかが分からず、また、当該画像がユーザの意図した精度を担保しているか分からないため、却って視認性が低下してしまうおそれがある。
【０００６】
本発明は、上記課題に鑑みてなされたもので、動態画像の解析結果の画像の視認性を向上させることができる動態画像解析装置及びプログラムを提供することを目的とする。
【課題を解決するための手段】
【０００７】
上記課題を解決するため、請求項１に記載の発明の動態画像解析装置は、
放射線撮影により得られた動態画像から所定の特徴量を解析し当該所定の特徴量を表示手段に表示する動態画像解析装置であって、
前記動態画像の信頼性を評価する評価手段と、
前記所定の特徴量が前記表示手段に表示される際に、当該所定の特徴量と対応付けて前記評価手段によって評価された前記動態画像の信頼性を示す信頼性情報を表示する表示制御手段と、
を備える。
【０００８】
請求項２に記載の発明は、請求項１に記載の発明において、
前記動態画像を構成する一連のフレーム画像から被験者の体動を検出する検出手段と、
前記検出手段によって検出された前記被験者の体動の量を算出する第１の算出手段と、を備え、
前記評価手段は、前記第１の算出手段によって算出された前記被験者の体動の量に基づき、前記動態画像の信頼性を評価する。
【０００９】
請求項３に記載の発明は、請求項２に記載の発明において、
前記評価手段は、前記第１の算出手段によって算出された前記被験者の体動の量が当該体動の量に関する閾値を超えている場合、前記動態画像の信頼性が低い旨の評価を行い、
前記表示制御手段は、前記評価手段によって前記動態画像の信頼性が低い旨の評価が行われた場合、当該動態画像の信頼性が低い旨を示す前記信頼性情報を前記表示手段に表示するとともに、前記被験者の体動の量が前記閾値を超えたときの当該体動の発生箇所を前記表示手段に表示する。
【００１０】
請求項４に記載の発明は、請求項３に記載の発明において、
前記表示制御手段は、前記評価手段によって前記動態画像の信頼性が低い旨の評価が行われた場合、更に、前記閾値を超えた体動の量を前記表示手段に表示する。
【００１１】
請求項５に記載の発明は、請求項１〜４のいずれか一項に記載の発明において、
前記動態画像の画素信号ごとに当該画素信号のＳＮ比を算出する第２の算出手段を備え、
前記評価手段は、前記第２の算出手段によって算出された前記ＳＮ比に基づき、前記動態画像の信頼性を評価する。
【００１２】
請求項６に記載の発明は、請求項５に記載の発明において、
前記評価手段は、前記第２の算出手段によって算出された前記ＳＮ比が当該ＳＮ比に関する閾値未満の場合、前記動態画像の信頼性が低い旨の評価を行い、
前記表示制御手段は、前記評価手段によって前記動態画像の信頼性が低い旨の評価が行われた場合、当該動態画像の信頼性が低い旨を示す前記信頼性情報を前記表示手段に表示するとともに、前記ＳＮ比が前記閾値未満である前記画素信号の画素位置を前記表示手段に表示する。
【００１３】
請求項７に記載の発明は、請求項１〜６のいずれか一項に記載の発明において、
前記動態画像の画素信号ごとに当該画素信号から抽出される前記所定の特徴量の解析対象となる周波数成分の信号と所定の理想の信号との相関係数を算出する第３の算出手段を備え、
前記評価手段は、前記第３の算出手段によって算出された前記相関係数に基づき、前記動態画像の信頼性を評価する。
【００１４】
請求項８に記載の発明は、請求項７に記載の発明において、
前記評価手段は、前記第３の算出手段によって算出された前記相関係数が当該相関係数に関する閾値未満の場合、前記動態画像の信頼性が低い旨の評価を行い、
前記表示制御手段は、前記評価手段によって前記動態画像の信頼性が低い旨の評価が行われた場合、当該動態画像の信頼性が低い旨を示す前記信頼性情報を前記表示手段に表示するとともに、前記相関係数が前記閾値未満である前記画素信号の画素位置を前記表示手段に表示する。
【００１５】
請求項９に記載の発明のプログラムは、
放射線撮影により得られた動態画像から所定の特徴量を解析し当該所定の特徴量を表示手段に表示する動態画像解析装置のコンピュータを、
前記動態画像の信頼性を評価する評価手段、
前記所定の特徴量が前記表示手段に表示される際に、当該所定の特徴量と対応付けて前記評価手段によって評価された前記動態画像の信頼性を示す信頼性情報を表示する表示制御手段、
として機能させる。
【発明の効果】
【００１６】
本発明によれば、動態画像の解析結果の画像の視認性を向上させることができる。
【図面の簡単な説明】
【００１７】
【図１】本発明の実施形態における動態画像解析システムの全体構成を示す図である。
【図２】図１の撮影用コンソールの制御部により実行される撮影制御処理を示すフローチャートである。
【図３】第１の実施形態において図１の診断用コンソールの制御部により実行される画像解析処理を示すフローチャートである。
【図４】第１の実施形態において図１の診断用コンソールの表示部に表示される画像解析結果の一例を示す図である。
【図５】第２の実施形態において図１の診断用コンソールの制御部により実行される画像解析処理を示すフローチャートである。
【図６】画素信号のＳＮ比の算出過程を説明するための図である。
【図７】第２の実施形態において図１の診断用コンソールの表示部に表示される画像解析結果の一例を示す図である。
【図８】第３の実施形態において図１の診断用コンソールの制御部により実行される画像解析処理を示すフローチャートである。
【図９】第３の実施形態において図１の診断用コンソールの表示部に表示される画像解析結果の一例を示す図である。
【発明を実施するための形態】
【００１８】
以下、図面を参照して本発明の実施形態について説明する。ただし、発明の範囲は、図示例に限定されない。
【００１９】
＜第１の実施形態＞
〔動態画像解析システム１００の構成〕
まず、第１の実施形態の構成を説明する。
図１に、第１の実施形態における動態画像解析システム１００の全体構成を示す。
図１に示すように、動態画像解析システム１００は、撮影装置１と、撮影用コンソール２とが通信ケーブル等により接続され、撮影用コンソール２と、診断用コンソール３とがＬＡＮ（Local Area Network）等の通信ネットワークＮＴを介して接続されて構成されている。動態画像解析システム１００を構成する各装置は、ＤＩＣＯＭ（Digital Image and Communications in Medicine）規格に準じており、各装置間の通信は、ＤＩＣＯＭに則って行われる。
【００２０】
〔撮影装置１の構成〕
撮影装置１は、例えば、呼吸運動に伴う肺の膨張及び収縮の形態変化、心臓の拍動等の、周期性（サイクル）を持つ動態を撮影する撮影手段である。動態撮影とは、被写体に対し、Ｘ線等の放射線をパルス状にして所定時間間隔で繰り返し照射するか(パルス照射)、もしくは、低線量率にして途切れなく継続して照射する(連続照射)ことで、被写体の動態を示す複数の画像を取得することをいう。動態撮影により得られた一連の画像を動態画像と呼ぶ。また、動態画像を構成する複数の画像のそれぞれをフレーム画像と呼ぶ。なお、以下の実施形態では、パルス照射により胸部の動態撮影を行う場合を例にとり説明する。
【００２１】
放射線源１１は、被写体Ｍを挟んで放射線検出部１３と対向する位置に配置され、放射線照射制御装置１２の制御に従って、被写体Ｍに対し放射線（Ｘ線）を照射する。
放射線照射制御装置１２は、撮影用コンソール２に接続されており、撮影用コンソール２から入力された放射線照射条件に基づいて放射線源１１を制御して放射線撮影を行う。撮影用コンソール２から入力される放射線照射条件は、例えば、パルスレート、パルス幅、パルス間隔、１撮影あたりの撮影フレーム数、Ｘ線管電流の値、Ｘ線管電圧の値、付加フィルタ種等である。パルスレートは、１秒あたりの放射線照射回数であり、後述するフレームレートと一致している。パルス幅は、放射線照射１回当たりの放射線照射時間である。パルス間隔は、１回の放射線照射開始から次の放射線照射開始までの時間であり、後述するフレーム間隔と一致している。
【００２２】
放射線検出部１３は、ＦＰＤ等の半導体イメージセンサーにより構成される。ＦＰＤは、例えば、ガラス基板等を有しており、基板上の所定位置に、放射線源１１から照射されて少なくとも被写体Ｍを透過した放射線をその強度に応じて検出し、検出した放射線を電気信号に変換して蓄積する複数の検出素子（画素）がマトリックス状に配列されている。各画素は、例えばＴＦＴ（Thin Film Transistor）等のスイッチング部を備えて構成されている。ＦＰＤにはＸ線をシンチレーターを介して光電変換素子により電気信号に変換する間接変換型、Ｘ線を直接的に電気信号に変換する直接変換型があるが、何れを用いてもよい。
放射線検出部１３は、被写体Ｍを挟んで放射線源１１と対向するように設けられている。
【００２３】
読取制御装置１４は、撮影用コンソール２に接続されている。読取制御装置１４は、撮影用コンソール２から入力された画像読取条件に基づいて放射線検出部１３の各画素のスイッチング部を制御して、当該各画素に蓄積された電気信号の読み取りをスイッチングしていき、放射線検出部１３に蓄積された電気信号を読み取ることにより、画像データを取得する。この画像データがフレーム画像である。そして、読取制御装置１４は、取得したフレーム画像を撮影用コンソール２に出力する。画像読取条件は、例えば、フレームレート、フレーム間隔、画素サイズ、画像サイズ（マトリックスサイズ）等である。フレームレートは、１秒あたりに取得するフレーム画像数であり、パルスレートと一致している。フレーム間隔は、１回のフレーム画像の取得動作開始から次のフレーム画像の取得動作開始までの時間であり、パルス間隔と一致している。
【００２４】
ここで、放射線照射制御装置１２と読取制御装置１４は互いに接続され、互いに同期信号をやりとりして放射線照射動作と画像の読み取りの動作を同調させるようになっている。
【００２５】
〔撮影用コンソール２の構成〕
撮影用コンソール２は、放射線照射条件や画像読取条件を撮影装置１に出力して撮影装置１による放射線撮影及び放射線画像の読み取り動作を制御するとともに、撮影装置１により取得された動態画像を撮影技師等の撮影実施者によるポジショニングの確認や診断に適した画像であるか否かの確認用に表示する。
撮影用コンソール２は、図１に示すように、制御部２１、記憶部２２、操作部２３、表示部２４、通信部２５を備えて構成され、各部はバス２６により接続されている。
【００２６】
制御部２１は、ＣＰＵ（Central Processing Unit）、ＲＡＭ（Random Access Memory
）等により構成される。制御部２１のＣＰＵは、操作部２３の操作に応じて、記憶部２２に記憶されているシステムプログラムや各種処理プログラムを読み出してＲＡＭ内に展開し、展開されたプログラムに従って後述する撮影制御処理を始めとする各種処理を実行し、撮影用コンソール２各部の動作や、撮影装置１の放射線照射動作及び読み取り動作を集中制御する。
【００２７】
記憶部２２は、不揮発性の半導体メモリやハードディスク等により構成される。記憶部２２は、制御部２１で実行される各種プログラムやプログラムにより処理の実行に必要なパラメータ、或いは処理結果等のデータを記憶する。例えば、記憶部２２は、図２に示す撮影制御処理を実行するためのプログラムを記憶している。また、記憶部２２は、被写体部位（ここでは、胸部とする）に対応付けて放射線照射条件及び画像読取条件を記憶している。各種プログラムは、読取可能なプログラムコードの形態で格納され、制御部２１は、当該プログラムコードに従った動作を逐次実行する。
【００２８】
操作部２３は、カーソルキー、数字入力キー、及び各種機能キー等を備えたキーボードと、マウス等のポインティングデバイスを備えて構成され、キーボードに対するキー操作やマウス操作により入力された指示信号を制御部２１に出力する。また、操作部２３は、表示部２４の表示画面にタッチパネルを備えても良く、この場合、タッチパネルを介して入力された指示信号を制御部２１に出力する。
【００２９】
表示部２４は、ＬＣＤ（Liquid Crystal Display）やＣＲＴ（Cathode Ray Tube）等のモニターにより構成され、制御部２１から入力される表示信号の指示に従って、操作部２３からの入力指示やデータ等を表示する。
【００３０】
通信部２５は、ＬＡＮアダプターやモデムやＴＡ（Terminal Adapter）等を備え、通信ネットワークＮＴに接続された各装置との間のデータ送受信を制御する。
【００３１】
〔診断用コンソール３の構成〕
診断用コンソール３は、撮影用コンソール２から動態画像を取得し、取得した動態画像や動態画像の解析結果を表示して医師の診断を支援するための動態画像解析装置である。
診断用コンソール３は、図１に示すように、制御部３１、記憶部３２、操作部３３、表示部３４、通信部３５を備えて構成され、各部はバス３６により接続されている。
【００３２】
制御部３１は、ＣＰＵ、ＲＡＭ等により構成される。制御部３１のＣＰＵは、操作部３３の操作に応じて、記憶部３２に記憶されているシステムプログラムや、各種処理プログラムを読み出してＲＡＭ内に展開し、展開されたプログラムに従って、後述する画像解析処理Ａを始めとする各種処理を実行し、診断用コンソール３の各部の動作を集中制御する。制御部３１は、評価手段、表示制御手段、検出手段、第１の算出手段として機能する。
【００３３】
記憶部３２は、不揮発性の半導体メモリやハードディスク等により構成される。記憶部３２は、制御部３１で画像解析処理Ａを実行するためのプログラムを始めとする各種プログラムやプログラムによる処理の実行に必要なパラメータ、或いは処理結果等のデータを記憶する。これらの各種プログラムは、読取可能なプログラムコードの形態で格納され、制御部３１は、当該プログラムコードに従った動作を逐次実行する。
【００３４】
また、記憶部３２には、過去に撮影された動態画像が患者情報（例えば、患者ＩＤ、患者の氏名、身長、体重、年齢、性別等）、検査情報（例えば、検査ＩＤ、検査日、被写体部位（ここでは、胸部）、診断対象の動態の種類（例えば、安静呼吸、深呼吸、心拍等）に対応付けて記憶されている。また、動態画像に対応する電子カルテ情報を図示しない電子カルテ装置から取得して動態画像に対応付けて記憶しておくこととしてもよい。
【００３５】
操作部３３は、カーソルキー、数字入力キー、及び各種機能キー等を備えたキーボードと、マウス等のポインティングデバイスを備えて構成され、ユーザによるキーボードに対するキー操作やマウス操作により入力された指示信号を制御部３１に出力する。また、操作部３３は、表示部３４の表示画面にタッチパネルを備えても良く、この場合、タッチパネルを介して入力された指示信号を制御部３１に出力する。
【００３６】
表示部３４は、ＬＣＤやＣＲＴ等のモニターにより構成され、制御部３１から入力される表示信号の指示に従って、各種表示を行う。
【００３７】
通信部３５は、ＬＡＮアダプターやモデムやＴＡ等を備え、通信ネットワークＮＴに接続された各装置との間のデータ送受信を制御する。
【００３８】
〔動態画像解析システム１００の動作〕
次に、本実施形態における上記動態画像解析システム１００の動作について説明する。
【００３９】
（撮影装置１、撮影用コンソール２の動作）
まず、撮影装置１、撮影用コンソール２による撮影動作について説明する。
図２に、撮影用コンソール２の制御部２１において実行される撮影制御処理を示す。撮影制御処理は、制御部２１と記憶部２２に記憶されているプログラムとの協働により実行される。
【００４０】
まず、撮影実施者により撮影用コンソール２の操作部２３が操作され、被検者の患者情報、検査情報の入力が行われる（ステップＳ１）。
【００４１】
次いで、放射線照射条件が記憶部２２から読み出されて放射線照射制御装置１２に設定されるとともに、画像読取条件が記憶部２２から読み出されて読取制御装置１４に設定される（ステップＳ２）。
【００４２】
次いで、操作部２３の操作による放射線照射の指示が待機される（ステップＳ３）。ここで、撮影実施者は、被写体Ｍを放射線源１１と放射線検出部１３の間に配置してポジショニングを行う。また、被検者に対し、診断対象の動態の種類に応じた呼吸状態を指示する。撮影準備が整った時点で、操作部２３を操作して放射線照射指示を入力する。
【００４３】
操作部２３により放射線照射指示が入力されると（ステップＳ３；ＹＥＳ）、放射線照射制御装置１２及び読取制御装置１４に撮影開始指示が出力され、動態撮影が開始される（ステップＳ４）。即ち、放射線照射制御装置１２に設定されたパルス間隔で放射線源１１により放射線が照射され、放射線検出部１３によりフレーム画像が取得される。
【００４４】
予め定められたフレーム数の撮影が終了すると、制御部２１により放射線照射制御装置１２及び読取制御装置１４に撮影終了の指示が出力され、撮影動作が停止される。撮影されるフレーム数は、少なくとも１呼吸サイクルが撮影できる枚数である。
【００４５】
撮影により取得されたフレーム画像は順次撮影用コンソール２に入力され、撮影順を示す番号（フレーム番号）と対応付けて記憶部２２に記憶されるとともに（ステップＳ５）、表示部２４に表示される（ステップＳ６）。撮影実施者は、表示された動態画像によりポジショニング等を確認し、撮影により診断に適した画像が取得された（撮影ＯＫ）か、再撮影が必要（撮影ＮＧ）か、を判断する。そして、操作部２３を操作して、判断結果を入力する。
【００４６】
操作部２３の所定の操作により撮影ＯＫを示す判断結果が入力されると（ステップＳ７；ＹＥＳ）、動態撮影で取得された一連のフレーム画像のそれぞれに、動態画像を識別するための識別ＩＤや、患者情報、検査情報、放射線照射条件、画像読取条件、撮影順を示す番号（フレーム番号）等の情報が付帯され（例えば、ＤＩＣＯＭ形式で画像データのヘッダ領域に書き込まれ）、通信部２５を介して診断用コンソール３に送信される（ステップＳ８）。そして、本処理は終了する。一方、操作部２３の所定の操作により撮影ＮＧを示す判断結果が入力されると（ステップＳ７；ＮＯ）、記憶部２２に記憶された一連のフレーム画像が削除され（ステップＳ９）、本処理は終了する。この場合、再撮影が必要となる。
【００４７】
（診断用コンソール３の動作）
次に、診断用コンソール３における動作について説明する。
診断用コンソール３においては、通信部３５を介して撮影用コンソール２から動態画像の一連のフレーム画像が受信されると、制御部３１と記憶部３２に記憶されているプログラムとの協働により図３に示す画像解析処理Ａが実行される。
【００４８】
画像解析処理Ａにおいては、まず、受信した動態画像から所定の特徴量が解析される（ステップＳ１１）。ここで、所定の特徴量としては、換気機能に係る特徴量や血流機能に係る特徴量が挙げられる。
【００４９】
次いで、解析対象の動態画像の一連のフレーム画像の各々が比較され（ステップＳ１２）、例えば、公知のフレーム間差分法を用いて被験者の体動領域の取得を試みることにより、被験者の体動が検出されたか否かが判定される（ステップＳ１３）。
ここで、体動としては、例えば、被写体部位（胸部）全体のずれ（横ずれ、前後ずれ等）や、呼吸による肋骨、外胸郭、乳房の移動などが挙げられるが、これらに限定されるものではない。
【００５０】
ステップＳ１３において、被験者の体動が検出されないと判定された場合（ステップＳ１３；ＮＯ）、ステップＳ１４〜ステップＳ１６をスキップして、ステップＳ１７の処理へ移行する。
一方、ステップＳ１３において、被験者の体動が検出されたと判定された場合（ステップＳ１３；ＹＥＳ）、検出された体動の量が算出される（ステップＳ１４）。
【００５１】
次いで、ステップＳ１４で算出された体動の量が当該体動の量に係る閾値を超えているか否かが判定される（ステップＳ１５）。
例えば、ステップＳ１３で検出された被験者の体動が被写体部位（胸部）全体のずれである場合には、当該被写体部位全体のずれの量に係る閾値（例えば、５ｃｍ）を用いてステップＳ１５の判定処理が行われる。また、ステップＳ１３で検出された被験者の体動が肋骨の移動である場合には、当該肋骨の移動の量に係る閾値（例えば、３ｃｍ）を用いてステップＳ１５の判定処理が行われる。
【００５２】
ステップＳ１５において、ステップＳ１４で算出された体動の量が当該体動の量に係る閾値を超えていないと判定された場合（ステップＳ１５；ＮＯ）、ステップＳ１６をスキップして、ステップＳ１７の処理へ移行する。
一方、ステップＳ１５において、ステップＳ１４で算出された体動の量が当該体動の量に係る閾値を超えていると判定された場合（ステップＳ１５；ＹＥＳ）、解析対象の動態画像について信頼性が低いとの評価がなされる（ステップＳ１６）。
【００５３】
次いで、動態画像の各フレーム画像の比較が終了したか否かが判定される（ステップＳ１７）。
【００５４】
ステップＳ１７において、動態画像の各フレーム画像の比較が終了していないと判定された場合（ステップＳ１７；ＮＯ）、ステップＳ１２の処理へ戻り、それ以降の処理が繰り返し行われる。
一方、ステップＳ１７において、動態画像の各フレーム画像の比較が終了したと判定された場合（ステップＳ１７；ＹＥＳ）、解析対象の動態画像について信頼性が低いとの評価がなされたか否かが判定される（ステップＳ１８）。
【００５５】
ステップＳ１８において、解析対象の動態画像について信頼性が低いとの評価がなされていると判定された場合（ステップＳ１８；ＹＥＳ）、ステップＳ１１で解析された特徴量Ｆを示す画像ＩＭ１と、ステップＳ１５で閾値を超えたと判定された体動の検出領域や当該体動の量などを示す体動情報Ｄ１、及び、信頼性が低いとする評価結果を示す信頼性情報Ｄ２が重畳された画像ＩＭ２と、が横に並べられた状態で表示部３４に表示され（ステップＳ１９）、画像解析処理Ａは終了する。ここで、画像ＩＭ２は、例えば、解析対象の動態画像のうちの一のフレーム画像であって、上記閾値を超える体動が検出されたときのフレーム画像であるものとする。
【００５６】
例えば、ステップＳ１３で検出された被験者の体動が肋骨の移動であり、当該肋骨の移動量が５ｃｍである場合、図４（ａ）に示すように、閾値（３ｃｍ）を超えたと判定された肋骨の検出領域や当該肋骨の移動量を示す体動情報Ｄ１、及び、信頼性が低いとする評価結果を示す信頼性情報Ｄ２が重畳された画像ＩＭ２が、ステップＳ１１で解析された特徴量Ｆを示す画像ＩＭ１と対応付けられて表示部３４に表示されることとなる。
また、上述の肋骨の移動に加えて、例えば、閾値を超える横ずれが検出されている場合、図４（ｂ）に示すように、横ずれの発生を示す体動情報Ｄ３が更に重畳された画像ＩＭ３が、上記の画像ＩＭ１と対応付けられて表示部３４に表示されることとなる。なお、横ずれの発生を示す体動情報Ｄ３には、当該横ずれのずれ量が含まれるようにし、当該体動情報Ｄ３が表示部３４に表示される際は、当該横ずれのずれ量も表示されるようにしてもよい。
【００５７】
また、ステップＳ１８において、解析対象の動態画像について信頼性が低いとの評価がなされていないと判定された場合（ステップＳ１８；ＮＯ）、ステップＳ１１で解析された特徴量Ｆを示す画像ＩＭ１のみが表示部３４に表示され（ステップＳ１９）、画像解析処理Ａは終了する。
【００５８】
このように、第１の実施形態においては、放射線撮影により得られた動態画像の信頼性が評価され、当該動態画像から解析された所定の特徴量Ｆを示す画像ＩＭ１が表示部３４に表示される際に、当該所定の特徴量Ｆと対応付けて当該動態画像の信頼性を示す信頼性情報Ｄ２が重畳された画像ＩＭ２が表示されるので、当該画像ＩＭ１に表示されている所定の特徴量Ｆの信頼性を明確にすることができ、当該所定の特徴量Ｆについて慎重な判断をすべきか示唆を与えることができる。したがって、第１の実施形態においては、所定の特徴量Ｆについて慎重な判断をすべきか示唆を与えることができるようになるので、当該所定の特徴量Ｆの視認性を向上させることができる。
【００５９】
また、第１の実施形態においては、動態画像を構成する一連のフレーム画像から被験者の体動が検出され、検出された被験者の体動の量が算出され、算出された被験者の体動の量に基づき、動態画像の信頼性が評価されるので、動態画像の信頼性を客観的に評価することができる。したがって、第１の実施形態においては、動態画像の信頼性を客観的に評価することができるので、当該動態画像の信頼性について精度の高い評価を行うことができる。
【００６０】
また、第１の実施形態においては、算出された被験者の体動の量が当該体動の量に関する閾値を超えている場合、動態画像の信頼性が低い旨の評価が行われ、動態画像の信頼性が低い旨の評価が行われた場合、当該動態画像の信頼性が低い旨を示す信頼性情報Ｄ２が表示部３４に表示されるとともに、被験者の体動の量が上記閾値を超えたときの当該体動の発生箇所（検出領域）を示す体動情報Ｄ１が表示部３４に表示されるので、画像ＩＭ１に表示されている所定の特徴量Ｆのうちどの領域の特徴量の信頼性が低いかをわかりやすくすることができる。したがって、第１の実施形態においては、所定の特徴量Ｆについて慎重な判断を要する箇所を示すことができるようになるので、当該所定の特徴量Ｆの視認性をより向上させることができる。
【００６１】
また、第１の実施形態においては、動態画像の信頼性が低い旨の評価が行われた場合、更に、上記体動情報Ｄ１として上記閾値を超えた体動の量が表示部３４に表示されるので、画像ＩＭ１に表示されている所定の特徴量Ｆのうちどの領域の特徴量の信頼性が低いかをわかりやすくすることができるとともに、体動の影響をどれだけ受けているかを明確にすることができる。したがって、第１の実施形態においては、所定の特徴量Ｆの慎重な判断を要する箇所について、更にどの程度慎重に判断すべきか示唆を与えることができるようになるので、当該所定の特徴量Ｆの視認性をより一層向上させることができる。
【００６２】
＜第２の実施形態＞
以下、本発明の第２の実施形態について説明する。
第２の実施形態における構成は、診断用コンソール３の記憶部３２に画像解析処理Ｂを実行するためのプログラムが記憶されている他は、第１の実施形態で説明したものと同様であるので説明を省略し、以下第２の実施形態の動作について説明する。
【００６３】
図５は、第２の実施形態において、診断用コンソール３により実行される画像解析処理Ｂを示すフローチャートである。画像解析処理Ｂは、制御部３１と記憶部３２に記憶されているプログラムとの協働により実行される。制御部３１は、評価手段、表示制御手段、第２の算出手段として機能する。
【００６４】
まず、画像解析処理Ａと同様、撮影用コンソール２より受信した動態画像から所定の特徴量Ｆを解析する（ステップＳ２１）。
【００６５】
次いで、解析対象の動態画像の画素信号ごとに当該画素信号のＳＮ比が算出される（ステップＳ２２）。
画素信号のＳＮ比の算出方法としては、例えば、図６に示すように、動態画像の各画素Ｐに対応する画素信号について周波数分析を行う。そして、当該周波数分析により得られた解析対象の周波数成分の信号強度とノイズ成分の信号強度とから画素信号のＳＮ比を算出する。なお、ＳＮ比が高ければノイズの影響が小さく、ＳＮ比が小さければノイズの影響が大きい。
【００６６】
次いで、ステップＳ２２で算出されたＳＮ比が当該ＳＮ比に係る閾値未満であるか否かが判定される（ステップＳ２３）。
【００６７】
ステップＳ２３において、ステップＳ２２で算出されたＳＮ比が当該ＳＮ比に係る閾値未満ではない、すなわちノイズの影響が小さいと判定された場合（ステップＳ２３；ＮＯ）、ステップＳ２４をスキップして、ステップＳ２５の処理へ移行する。
一方、ステップＳ２３において、ステップＳ２２で算出されたＳＮ比が当該ＳＮ比に係る閾値未満である、すなわちノイズの影響が大きいと判定された場合（ステップＳ２３；ＹＥＳ）、解析対象の動態画像について信頼性が低いとの評価がなされる（ステップＳ２４）。
【００６８】
次いで、全ての画素信号のＳＮ比が算出されたか否かが判定される（ステップＳ２５）。
【００６９】
ステップＳ２５において、全ての画素信号のＳＮ比が算出されていないと判定された場合（ステップＳ２５；ＮＯ）、ステップＳ２２の処理へ戻り、それ以降の処理が繰り返し行われる。
一方、ステップＳ２５において、全ての画素信号のＳＮ比が算出されたと判定された場合（ステップＳ２５；ＹＥＳ）、解析対象の動態画像について信頼性が低いとの評価がなされたか否かが判定される（ステップＳ２６）。
【００７０】
ステップＳ２６において、解析対象の動態画像について信頼性が低いとの評価がなされていると判定された場合（ステップＳ２６；ＹＥＳ）、ステップＳ２１で解析された特徴量Ｆを示す画像ＩＭ１と、ステップＳ２３で閾値未満であると判定された画素信号の領域情報Ｒ１、並びに、信頼性が低いとする評価結果及び当該信頼性が低いと評価された要因を示す信頼性情報Ｄ４が重畳された画像ＩＭ４と、が横に並べられた状態で表示部３４に表示され（ステップＳ９）、画像解析処理Ｂは終了する。ここで、画像ＩＭ４は、解析対象の動態画像のうちの一のフレーム画像であるものとする。
【００７１】
例えば、画素信号のＳＮ比が閾値未満であると判定された画素が右肺の下部付近に集中している場合、図７に示すように、画素信号のＳＮ比が閾値未満であると判定された画素信号の領域情報Ｒ１、並びに、信頼性が低いとする評価結果及び当該信頼性が低いと評価された要因（ノイズ）を示す信頼性情報Ｄ４が重畳された画像ＩＭ４が、ステップＳ２１で解析された特徴量Ｆを示す画像ＩＭ１と対応付けられて表示部３４に表示されることとなる。
【００７２】
また、ステップＳ２６において、解析対象の動態画像について信頼性が低いとの評価がなされていないと判定された場合（ステップＳ２６；ＮＯ）、ステップＳ２１で解析された特徴量を示す画像ＩＭ１のみが表示部３４に表示され（ステップＳ２８）、画像解析処理Ｂは終了する。
【００７３】
このように、第２の実施形態においては、動態画像の画素信号ごとに当該画素信号のＳＮ比が算出され、算出されたＳＮ比に基づき、動態画像の信頼性を評価するので、動態画像の信頼性を客観的に評価することができる。したがって、第２の実施形態においては、動態画像の信頼性を客観的に評価することができるので、当該動態画像の信頼性について精度の高い評価を行うことができる。
【００７４】
また、第２の実施形態においては、算出されたＳＮ比が当該ＳＮ比に関する閾値未満の場合、動態画像の信頼性が低い旨の評価が行われ、動態画像の信頼性が低い旨の評価が行われた場合、当該動態画像の信頼性が低い旨を示す信頼性情報Ｄ４が表示部３４に表示されるとともに、ＳＮ比が上記閾値未満である画素信号の領域情報Ｒ１が表示部３４に表示されるので、画像ＩＭ１に表示されている所定の特徴量Ｆのうちどの領域の特徴量の信頼性が低いかをわかりやすくすることができる。したがって、第２の実施形態においては、所定の特徴量Ｆについて慎重な判断を要する箇所を示すことができるようになるので、当該所定の特徴量Ｆの視認性をより向上させることができる。
【００７５】
＜第３の実施形態＞
以下、本発明の第３の実施形態について説明する。
第３の実施形態における構成は、診断用コンソール３の記憶部３２に画像解析処理Ｃを実行するためのプログラムが記憶されている他は、第１の実施形態で説明したものと同様であるので説明を省略し、以下第３の実施形態の動作について説明する。
【００７６】
図８は、第３の実施形態において、診断用コンソール３により実行される画像解析処理Ｃを示すフローチャートである。画像解析処理Ｃは、制御部３１と記憶部３２に記憶されているプログラムとの協働により実行される。制御部３１は、評価手段、表示制御手段、第３の算出手段として機能する。
【００７７】
まず、画像解析処理Ａと同様、撮影用コンソール２より受信した動態画像から所定の特徴量を解析する（ステップＳ３１）。
【００７８】
次いで、解析対象の動態画像の画素信号ごとに、当該画素信号のうちの解析対象の周波数成分の信号と理想の信号との相関係数が算出される（ステップＳ３２）。
相関係数の算出方法としては、例えば、図６に示すように、動態画像の各画素に対応する画素信号について周波数分析を行う。そして、当該周波数分析により得られた解析対象の周波数成分の信号を抽出し、抽出された信号と予め設定がなされている理想の信号との相関係数を算出する。
【００７９】
次いで、ステップＳ３２で算出された相関係数が当該相関係数に係る閾値未満であるか否かが判定される（ステップＳ３３）。
【００８０】
ステップＳ３３において、ステップＳ３２で算出された相関係数が当該相関係数に係る閾値未満ではない、すなわち理想の信号から逸脱していないと判定された場合（ステップＳ３３；ＮＯ）、ステップＳ３４をスキップして、ステップＳ３５の処理へ移行する。
一方、ステップＳ３３において、ステップＳ３２で算出された相関係数が当該相関係数に係る閾値未満である、すなわち理想の信号から逸脱していると判定された場合（ステップＳ３３；ＹＥＳ）、解析対象の動態画像について信頼性が低いとの評価がなされる（ステップＳ３４）。
【００８１】
次いで、全ての画素信号について相関係数が算出されたか否かが判定される（ステップＳ３５）。
【００８２】
ステップＳ３５において、全ての画素信号について相関係数が算出されていないと判定された場合（ステップＳ３５；ＮＯ）、ステップＳ３２の処理へ戻り、それ以降の処理が繰り返し行われる。
一方、ステップＳ３５において、全ての画素信号について相関係数が算出されたと判定された場合（ステップＳ３５；ＹＥＳ）、解析対象の動態画像について信頼性が低いとの評価がなされたか否かが判定される（ステップＳ３６）。
【００８３】
ステップＳ３６において、解析対象の動態画像について信頼性が低いとの評価がなされていると判定された場合（ステップＳ３６；ＹＥＳ）、ステップＳ３１で解析された特徴量を示す画像ＩＭ１と、ステップＳ３３で閾値未満であると判定された画素信号の領域情報Ｒ２、並びに、信頼性が低いとする評価結果及び当該信頼性が低いと評価された要因を示す信頼性情報Ｄ５が重畳された画像ＩＭ５と、が横に並べられた状態で表示部３４に表示され（ステップＳ３７）、画像解析処理Ｃは終了する。ここで、画像ＩＭ５は、解析対象の動態画像のうちの一のフレーム画像であるものとする。
【００８４】
例えば、相関係数が閾値未満であると判定された画素が右肺の下部付近に集中している場合、図９に示すように、相関係数が閾値未満であると判定された画素信号の領域情報Ｒ２、並びに、信頼性が低いとする評価結果及び当該信頼性が低いと評価された要因（理想の信号から逸脱）を示す信頼性情報Ｄ５が重畳された画像ＩＭ５が、ステップＳ３１で解析された特徴量Ｆを示す画像ＩＭ１と対応付けられて表示部３４に表示されることとなる。
【００８５】
また、ステップＳ３６において、解析対象の動態画像について信頼性が低いとの評価がなされていないと判定された場合（ステップＳ３６；ＮＯ）、ステップＳ３１で解析された特徴量を示す画像ＩＭ１のみが表示部３４に表示され（ステップＳ３８）、画像解析処理Ｃは終了する。
【００８６】
このように、第３の実施形態においては、動態画像の画素信号ごとに当該画素信号から抽出される所定の特徴量Ｆの解析対象となる周波数成分の信号と所定の理想の信号との相関係数が算出され、算出された相関係数に基づき、動態画像の信頼性を評価するので、動態画像の信頼性を客観的に評価することができる。したがって、第３の実施形態においては、動態画像の信頼性を客観的に評価することができるので、当該動態画像の信頼性について精度の高い評価を行うことができる。
【００８７】
また、第３の実施形態においては、算出された相関係数が当該相関係数に関する閾値未満の場合、動態画像の信頼性が低い旨の評価が行われ、動態画像の信頼性が低い旨の評価が行われた場合、当該動態画像の信頼性が低い旨を示す信頼性情報Ｄ５が表示部３４に表示されるとともに、相関係数が上記閾値未満である画素信号の領域情報Ｒ２が表示部３４に表示されるので、画像ＩＭ１に表示されている所定の特徴量Ｆのうちどの領域の特徴量の信頼性が低いかをわかりやすくすることができる。したがって、第３の実施形態においては、所定の特徴量Ｆについて慎重な判断を要する箇所を示すことができるようになるので、当該所定の特徴量Ｆの視認性をより向上させることができる。
【００８８】
以上、本発明の第１〜第３の実施形態について説明したが、実施形態における記述内容は、本発明の好適な一例であり、これに限定されるものではない。
【００８９】
例えば、上記実施形態においては、被写体部位が胸部である場合を例にとり説明したが、他の部位を撮影した動態画像について解析処理を行う場合においても本発明を適用することができる。
【００９０】
また、上記実施形態１においては、診断用コンソール３の制御部３１が画像解析処理Ａを実行する場合を例にとり説明したが、例えば、制御部３１が画像解析処理Ａに加えて、第２の実施形態で説明した画像解析処理Ｂと第３の実施形態で説明した画像解析処理Ｃのうちの少なくとも何れか一方の処理を合わせて実行するようにしてもよい。
【００９１】
また、上記実施形態２においては、診断用コンソール３の制御部３１が画像解析処理Ｂを実行する場合を例にとり説明したが、例えば、制御部３１が画像解析処理Ｂに加えて、第３の実施形態で説明した画像解析処理Ｃを合わせて実行するようにしてもよい。
【００９２】
また、上記実施形態１においては、診断用コンソール３の制御部３１で算出された被験者の体動の量が当該体動の量に関する閾値を超えている場合、動態画像の信頼性が低い旨の評価が行われるようにしているが、例えば、各フレーム画像の比較が終了した段階で、被験者の体動の量が当該体動の量に関する閾値を超えたとする判定が一度もなされなかった場合には、動態画像の信頼性が高い旨の評価が行われるようにしてもよい。かかる場合、例えば、少なくとも信頼性が高いとする評価結果を示す信頼性情報と、診断用コンソール３の制御部３１で解析された所定の特徴量Ｆを示す画像ＩＭ１とが横に並べられた状態で表示部３４に表示されるようにする。
【００９３】
また、上記実施形態２においては、診断用コンソール３の制御部３１で算出されたＳＮ比が当該ＳＮ比に関する閾値未満の場合、動態画像の信頼性が低い旨の評価が行われるようにしているが、例えば、診断用コンソール３の制御部３１で算出されたＳＮ比が当該ＳＮ比に関する閾値以上であると判定された画素信号が存在する場合には、少なくとも当該画素信号の領域情報と信頼性が高いとする評価結果を示す信頼性情報とを含む画像と、診断用コンソール３で解析された所定の特徴量Ｆを示す画像ＩＭ１とが横に並べられた状態で表示部３４に表示されるようにする。
【００９４】
また、上記実施形態３においては、診断用コンソール３の制御部３１で算出された相関係数が当該相関係数に関する閾値未満の場合、動態画像の信頼性が低い旨の評価が行われるようにしているが、例えば、診断用コンソール３の制御部３１で算出された相関係数が当該相関係数に関する閾値以上であると判定された画素信号が存在する場合には、少なくとも当該画素信号の領域情報と信頼性が高いとする評価結果を示す信頼性情報とを含む画像と、診断用コンソール３で解析された特徴量Ｆを示す画像ＩＭ１とが横に並べられた状態で表示部３４に表示されるようにする。
【００９５】
また、例えば、上記の説明では、本発明に係るプログラムのコンピュータ読み取り可能な媒体としてハードディスクや半導体の不揮発性メモリ等を使用した例を開示したが、この例に限定されない。その他のコンピュータ読み取り可能な媒体として、ＣＤ−ＲＯＭ等の可搬型記録媒体を適用することが可能である。また、本発明に係るプログラムのデータを通信回線を介して提供する媒体として、キャリアウエーブ（搬送波）も適用される。
【００９６】
その他、動態画像解析システムを構成する各装置の細部構成及び細部動作に関しても、本発明の趣旨を逸脱することのない範囲で適宜変更可能である。
【符号の説明】
【００９７】
１００  動態画像解析システム
１  撮影装置
１１  放射線源
１２  放射線照射制御装置
１３  放射線検出部
１４  読取制御装置
２  撮影用コンソール
２１  制御部
２２  記憶部
２３  操作部
２４  表示部
２５  通信部
２６  バス
３  診断用コンソール
３１  制御部
３２  記憶部
３３  操作部
３４  表示部
３５  通信部
３６  バス
【図１】
2020000807.tif 000003
【図２】
2020000807.tif 000004
【図３】
2020000807.tif 000005
【図４】
2020000807.tif 000006
【図５】
2020000807.tif 000007
【図６】
2020000807.tif 000008
【図７】
2020000807.tif 000009
【図８】
2020000807.tif 000010
【図９】
2020000807.tif 000011