(57)【要約】【課題】動態画像またはその解析結果画像において、呼吸や心拍に伴う構造物等の移動によって生じるアーチファクトをより短時間の処理で精度良く減弱する。【解決手段】診断用コンソール３の制御部３１によれば、胸部の動態画像に注目領域を設定し、注目領域の移動に関する移動情報を取得し、取得した注目領域の移動情報に基づいて、注目領域に対応する注目画素の画素信号値の強調度を決定する。そして、決定した強調度に基づいて、動態画像を解析することにより生成されたフレーム間差分画像等の解析結果画像の注目画素の画素信号値を補正する。【選択図】図４
2018007949.tif 000002
【特許請求の範囲】
【請求項１】
検査対象部位に放射線を照射して動態撮影することにより得られた動態画像に注目領域を設定し、前記注目領域の移動に関する移動情報を取得する取得手段と、
前記取得手段により取得された前記注目領域の移動情報に基づいて、前記注目領域内の画素である注目画素の画素信号値の強調度を決定する決定手段と、
前記決定手段により決定された強調度に基づいて、前記動態画像又は前記動態画像を解析することにより生成された解析結果画像の前記注目画素の画素信号値を補正する補正手段と、
を備える動態解析装置。
【請求項２】
前記移動情報は、移動量を含み、
前記決定手段は、前記注目領域の移動量が大きいほど前記強調度を低い強調度に決定する請求項１に記載の動態解析装置。
【請求項３】
前記動態画像に参照領域を設定し、前記参照領域に含まれる構造物の移動に関する移動情報を取得する第２の取得手段を備え、
前記決定手段は、前記注目領域の移動情報及び前記参照領域に含まれる構造物の移動情報に基づいて、前記強調度を決定する請求項１又は２に記載の動態解析装置。
【請求項４】
前記注目領域の前記移動情報と、前記参照領域に含まれる構造物の移動に関する移動情報は、移動方向を含み、
前記決定手段は、前記注目領域の移動方向と前記参照領域に含まれる構造物の移動方向の角度差に基づいて、前記強調度を決定する請求項３に記載の動態解析装置。
【請求項５】
前記動態画像は、胸部の動態画像であり、
前記第２の取得手段は、前記動態画像の横隔膜を含む領域に参照領域を設定して横隔膜の移動方向を取得し、
前記決定手段は、前記注目領域の移動方向と前記横隔膜の移動方向の角度差が大きいほど前記強調度を低い強調度に決定する請求項４に記載の動態解析装置。
【請求項６】
前記動態画像は、胸部の動態画像であり、
前記第２の取得手段は、前記動態画像の肩を含む領域に参照領域を設定して肩の移動方向を取得し、
前記決定手段は、前記注目領域の移動方向と前記肩の移動方向の角度差が小さいほど前記強調度を低い強調度に決定する請求項４に記載の動態解析装置。
【請求項７】
前記移動情報は、移動方向を含み、
前記動態画像に参照領域を設定し、前記参照領域の画素信号値の変化に基づいて所定の構造物の移動方向を特定する特定手段を備え、
前記決定手段は、前記注目領域の移動方向と前記所定の構造物の移動方向の角度差に基づいて前記強調度を決定する請求項１又は２に記載の動態解析装置。
【請求項８】
前記動態画像は、胸部の動態画像であり、
前記特定手段は、前記動態画像の肋間領域に参照領域を設定し、前記参照領域の画素信号値の変化に基づいて肋骨の移動方向を特定し、
前記決定手段は、前記注目領域の移動方向と前記肋骨の移動方向の角度差が小さいほど前記強調度を低い強調度に決定する請求項７に記載の動態解析装置。
【請求項９】
前記決定手段は、前記注目領域の移動情報に対する強調度の対応テーブル又はグラフによって前記強調度を決定する請求項１又は２に記載の動態解析装置。
【請求項１０】
前記決定手段は、前記注目領域の移動情報及び前記角度差に対する強調度の対応テーブル又はグラフによって前記強調度を決定する請求項４〜８の何れか一項に記載の動態解析装置。
【請求項１１】
前記決定手段は、前記注目領域の移動情報及びπから前記角度差を引いた角度に対する強調度の対応テーブル又はグラフによって前記強調度を決定する請求項４〜８の何れか一項に記載の動態解析装置。
【請求項１２】
前記解析結果画像は、前記動態画像のフレーム画像間で対応する画素同士の差分値を算出することにより生成されるフレーム間差分画像である請求項１〜１１の何れか一項に記載の動態解析装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、動態解析装置に関する。
【背景技術】
【０００２】
従来、胸部の動態を撮影し、得られた動態画像を解析して診断に提供する技術が知られている。しかし、胸部の動態画像には呼吸や心拍に伴う骨格等の動きが含まれており、これらの動きによるアーチファクトは、例えば、フレーム間差分等の差分処理を基本とした動態解析処理全般において強調され、注目すべき信号の観察の邪魔になり、解析精度及び診断精度を低下させていた。
【０００３】
そこで、上記の解決策として、例えば、非特許文献１には、胸部動態画像（オリジナル画像）に肋骨陰影低減処理を実行して軟組織動画像を生成し、生成した軟組織動画像から肺換気マッピング画像（フレーム間差分画像）を生成することで、オリジナル画像から作成された肺換気マッピング画像上で見られた肋骨アーチファクトが低減されたことが記載されている。また、特許文献１においても、Ｘ線動画像から骨陰影を除去して軟組織に関するＸ線動画像を生成し、生成したＸ線動画像に基づいて解析を行うことで、軟組織を高精度に解析できることが記載されている。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特開２０１５−４３８９４号公報
【非特許文献】
【０００５】
【非特許文献１】田中 利恵、真田茂著「低コスト・低被ばくポータブルＸ線肺機能イメージング (診る聴診器) の開発 」医用画像情報学会雑誌、 31 ( 2 ) ，2014
【発明の概要】
【発明が解決しようとする課題】
【０００６】
しかしながら、肋骨陰影低減処理は肋骨陰影を確実に低減させる処理ではなく、ある領域では肋骨が残ることがあり、その影響を受けアーチファクトが除去しきれない。その結果、場所によってアーチファクトの出方が変わるためユーザは処理の挙動を理解しづらい。特に限られた時間の中で大量の読影を行わなければならない医師にとって、いちいちアーチファクトの出方の強弱を勘案しながら読影することは煩雑であり、かえって診断効率を下げる結果となる。
【０００７】
また、非特許文献１には、肋骨陰影低減処理の処理時間が１枚あたり１５秒かかることが報告されており、動態画像の全てのフレーム画像へ適用した場合、撮影完了後速やかに解析結果を参照しながら診断することはできない。
【０００８】
本発明の課題は、動態画像またはその解析結果画像において、呼吸や心拍に伴う構造物等の移動によって生じるアーチファクトをより短時間の処理で精度良く減弱することである。
【課題を解決するための手段】
【０００９】
上記課題を解決するため、請求項１に記載の発明の動態解析装置は、
検査対象部位に放射線を照射して動態撮影することにより得られた動態画像に注目領域を設定し、前記注目領域の移動に関する移動情報を取得する取得手段と、
前記取得手段により取得された前記注目領域の移動情報に基づいて、前記注目領域内の画素である注目画素の画素信号値の強調度を決定する決定手段と、
前記決定手段により決定された強調度に基づいて、前記動態画像又は前記動態画像を解析することにより生成された解析結果画像の前記注目画素の画素信号値を補正する補正手段と、
を備える。
【００１０】
請求項２に記載の発明は、請求項１に記載の発明において、
前記移動情報は、移動量を含み、
前記決定手段は、前記注目領域の移動量が大きいほど前記強調度を低い強調度に決定する。
【００１１】
請求項３に記載の発明は、請求項１又は２に記載の発明において、
前記動態画像に参照領域を設定し、前記参照領域に含まれる構造物の移動に関する移動情報を取得する第２の取得手段を備え、
前記決定手段は、前記注目領域の移動情報及び前記参照領域に含まれる構造物の移動情報に基づいて、前記強調度を決定する。
【００１２】
請求項４に記載の発明は、請求項３に記載の発明において、
前記注目領域の前記移動情報と、前記参照領域に含まれる構造物の移動に関する移動情報は、移動方向を含み、
前記決定手段は、前記注目領域の移動方向と前記参照領域に含まれる構造物の移動方向の角度差に基づいて、前記強調度を決定する。
【００１３】
請求項５に記載の発明は、請求項４に記載の発明において、
前記動態画像は、胸部の動態画像であり、
前記第２の取得手段は、前記動態画像の横隔膜を含む領域に参照領域を設定して横隔膜の移動方向を取得し、
前記決定手段は、前記注目領域の移動方向と前記横隔膜の移動方向の角度差が大きいほど前記強調度を低い強調度に決定する。
【００１４】
請求項６に記載の発明は、請求項４に記載の発明において、
前記動態画像は、胸部の動態画像であり、
前記第２の取得手段は、前記動態画像の肩を含む領域に参照領域を設定して肩の移動方向を取得し、
前記決定手段は、前記注目領域の移動方向と前記肩の移動方向の角度差が小さいほど前記強調度を低い強調度に決定する。
【００１５】
請求項７に記載の発明は、請求項１又は２に記載の発明において、
前記移動情報は、移動方向を含み、
前記動態画像に参照領域を設定し、前記参照領域の画素信号値の変化に基づいて所定の構造物の移動方向を特定する特定手段を備え、
前記決定手段は、前記注目領域の移動方向と前記所定の構造物の移動方向の角度差に基づいて前記強調度を決定する。
【００１６】
請求項８に記載の発明は、請求項７に記載の発明において、
前記動態画像は、胸部の動態画像であり、
前記特定手段は、前記動態画像の肋間領域に参照領域を設定し、前記参照領域の画素信号値の変化に基づいて肋骨の移動方向を特定し、
前記決定手段は、前記注目領域の移動方向と前記肋骨の移動方向の角度差が小さいほど前記強調度を低い強調度に決定する。
【００１７】
請求項９に記載の発明は、請求項１又は２に記載の発明において、
前記決定手段は、前記注目領域の移動情報に対する強調度の対応テーブル又はグラフによって前記強調度を決定する。
【００１８】
請求項１０に記載の発明は、請求項４〜８の何れか一項に記載の発明において、
前記決定手段は、前記注目領域の移動情報及び前記角度差に対する強調度の対応テーブル又はグラフによって前記強調度を決定する。
【００１９】
請求項１１に記載の発明は、請求項４〜８の何れか一項に記載の発明において、
前記決定手段は、前記注目領域の移動情報及びπから前記角度差を引いた角度に対する強調度の対応テーブル又はグラフによって前記強調度を決定する。
【００２０】
請求項１２に記載の発明は、請求項１〜１１の何れか一項に記載の発明において、
前記解析結果画像は、前記動態画像のフレーム画像間で対応する画素同士の差分値を算出することにより生成されるフレーム間差分画像である。
【発明の効果】
【００２１】
本発明によれば、動態画像またはその解析結果画像において、呼吸や心拍に伴う構造物等の移動によって生じるアーチファクトをより短時間の処理で精度良く減弱することが可能となる。
【図面の簡単な説明】
【００２２】
【図１】本発明の実施形態における動態解析システムの全体構成を示す図である。
【図２】図１の撮影用コンソールの制御部により実行される撮影制御処理を示すフローチャートである。
【図３】図１の診断用コンソールの制御部により実行される画像解析処理を示すフローチャートである。
【図４】第１の実施形態において図３のステップＳ１２において実行される画素信号値補正処理Ａを示すフローチャートである。
【図５】強調度テーブルの一例を示す図である。
【図６】強調度テーブルの一例を示す図である。
【図７】強調度テーブルの一例を示す図である。
【図８】第２の実施形態において図３のステップＳ１２において実行される画素信号値補正処理Ｂを示すフローチャートである。
【図９】横隔膜への参照領域の設定例を示す図である。
【図１０】肩への参照領域の設定例を示す図である。
【図１１】移動方向の他の軸方向への投影例を示す図である。
【図１２】異方度を説明するための図である。
【図１３】同方向度を説明するための図である。
【図１４】第３の実施形態において図３のステップＳ１２で実行される画素信号値補正処理Ｃを示すフローチャートである。
【図１５】肋間への参照領域の設定例を示す図である。
【発明を実施するための形態】
【００２３】
以下、図面を参照して本発明の実施形態を詳細に説明する。ただし、発明の範囲は、図示例に限定されない。
【００２４】
＜第１の実施形態＞
〔動態解析システム１００の構成〕
まず、構成を説明する。
図１に、本発明の第１の実施形態における動態解析システム１００の全体構成を示す。
図１に示すように、動態解析システム１００は、撮影装置１と、撮影用コンソール２とが通信ケーブル等により接続され、撮影用コンソール２と、診断用コンソール３とがＬＡＮ（Local Area Network）等の通信ネットワークＮＴを介して接続されて構成されている。なお、本実施形態においては、撮影用コンソール２及び診断用コンソール３が通信ネットワークＮＴに有線により接続されている場合を例にとり説明するが、無線により接続されていることとしてもよい。
【００２５】
〔撮影装置１の構成〕
撮影装置１は、例えば、呼吸運動に伴う肺の膨張及び収縮の形態変化、心臓の拍動等の、周期性（サイクル）を持つ人体の動態を撮影する撮影手段である。動態撮影とは、被写体Ｍに対し、Ｘ線等の放射線をパルス状にして所定時間間隔で繰り返し照射するか(パルス照射)、もしくは、低線量率にして途切れなく継続して照射する(連続照射)ことで、動態を示す複数の画像を取得することをいう。動態撮影により得られた一連の画像を動態画像と呼ぶ。また、動態画像を構成する複数の画像のそれぞれをフレーム画像と呼ぶ。なお、以下の実施形態では、パルス照射により動態撮影を行う場合を例にとり説明する。
【００２６】
放射線源１１は、被写体Ｍを挟んで放射線検出部１３と対向する位置に配置され、放射線照射制御装置１２の制御に従って、被写体Ｍに対し放射線（Ｘ線）を照射する。
放射線照射制御装置１２は、撮影用コンソール２に接続されており、撮影用コンソール２から入力された放射線照射条件に基づいて放射線源１１を制御して放射線撮影を行う。撮影用コンソール２から入力される放射線照射条件は、例えば、パルスレート、パルス幅、パルス間隔、１撮影あたりの撮影フレーム数、Ｘ線管電流の値、Ｘ線管電圧の値、グリッドの種類、付加フィルタ種等である。パルスレートは、１秒あたりの放射線照射回数であり、後述するフレームレートと一致している。パルス幅は、放射線照射１回当たりの放射線照射時間である。パルス間隔は、１回の放射線照射開始から次の放射線照射開始までの時間であり、後述するフレーム間隔と一致している。
【００２７】
放射線検出部１３は、ＦＰＤ等の半導体イメージセンサーにより構成される。ＦＰＤは、例えば、ガラス基板等を有しており、基板上の所定位置に、放射線源１１から照射されて少なくとも被写体Ｍを透過した放射線をその強度に応じて検出し、検出した放射線を電気信号に変換して蓄積する複数の検出素子（画素）がマトリックス状に配列されている。各画素は、例えばＴＦＴ（Thin Film Transistor）等のスイッチング部を備えて構成されている。ＦＰＤにはＸ線をシンチレーターを介して光電変換素子により電気信号に変換する間接変換型、Ｘ線を直接的に電気信号に変換する直接変換型があるが、何れを用いてもよい。放射線検出部１３は、被写体Ｍを挟んで放射線源１１と対向するように設けられている。
【００２８】
ここで、放射線検出部１３において取得される画像の各画素の信号値（画素信号値）は、放射線検出部１３に到達した放射線強度を電気信号に変換した値、即ち、放射線検出部１３に到達した放射線強度に相関する値であり、到達した放射線強度が高いほど画素信号値は大きくなる。一方、本実施形態において、撮影用コンソール２及び診断用コンソール３では、各画素信号値を吸収線量を表すものとして取り扱う。即ち、撮影用コンソール２は、通信部２５により受信した一連のフレーム画像の各画素信号値を吸収線量を表す値に変換し、診断用コンソール３に送信する。即ち、撮影用コンソール２及び診断用コンソールでは、動態画像の画素信号値は、吸収線量が多いほど高くなり、画像上では白く（低濃度で）描画される。
【００２９】
読取制御装置１４は、撮影用コンソール２に接続されている。読取制御装置１４は、撮影用コンソール２から入力された画像読取条件に基づいて放射線検出部１３の各画素のスイッチング部を制御して、当該各画素に蓄積された電気信号の読み取りをスイッチングしていき、放射線検出部１３に蓄積された電気信号を読み取ることにより、画像データを取得する。この画像データがフレーム画像である。そして、読取制御装置１４は、取得したフレーム画像を撮影用コンソール２に出力する。画像読取条件は、例えば、フレームレート、フレーム間隔、画素サイズ、画像サイズ（マトリックスサイズ）等である。フレームレートは、１秒あたりに取得するフレーム画像数であり、パルスレートと一致している。フレーム間隔は、１回のフレーム画像の取得動作開始から次のフレーム画像の取得動作開始までの時間であり、パルス間隔と一致している。
【００３０】
ここで、放射線照射制御装置１２と読取制御装置１４は互いに接続され、互いに同期信号をやりとりして放射線照射動作と画像の読み取りの動作を同調させるようになっている。
【００３１】
〔撮影用コンソール２の構成〕
撮影用コンソール２は、放射線照射条件や画像読取条件を撮影装置１に出力して撮影装置１による放射線撮影及び放射線画像の読み取り動作を制御するとともに、撮影装置１により取得された動態画像を撮影技師等の撮影実施者によるポジショニングの確認や診断に適した画像であるか否かの確認用に表示する。
撮影用コンソール２は、図１に示すように、制御部２１、記憶部２２、操作部２３、表示部２４、通信部２５を備えて構成され、各部はバス２６により接続されている。
【００３２】
制御部２１は、ＣＰＵ（Central Processing Unit）、ＲＡＭ（Random Access Memory
）等により構成される。制御部２１のＣＰＵは、操作部２３の操作に応じて、記憶部２２に記憶されているシステムプログラムや各種処理プログラムを読み出してＲＡＭ内に展開し、展開されたプログラムに従って後述する撮影制御処理を始めとする各種処理を実行し、撮影用コンソール２各部の動作や、撮影装置１の放射線照射動作及び読み取り動作を集中制御する。
【００３３】
記憶部２２は、不揮発性の半導体メモリーやハードディスク等により構成される。記憶部２２は、制御部２１で実行される各種プログラムやプログラムにより処理の実行に必要なパラメーター、或いは処理結果等のデータを記憶する。例えば、記憶部２２は、図２に示す撮影制御処理を実行するためのプログラムを記憶している。また、記憶部２２は、検査対象部位（ここでは胸部）及び診断対象の種類（例えば、換気、血流）に対応付けて放射線照射条件及び画像読取条件を記憶している。各種プログラムは、読取可能なプログラムコードの形態で格納され、制御部２１は、当該プログラムコードに従った動作を逐次実行する。
【００３４】
操作部２３は、カーソルキー、数字入力キー、及び各種機能キー等を備えたキーボードと、マウス等のポインティングデバイスを備えて構成され、キーボードに対するキー操作やマウス操作により入力された指示信号を制御部２１に出力する。また、操作部２３は、表示部２４の表示画面にタッチパネルを備えても良く、この場合、タッチパネルを介して入力された指示信号を制御部２１に出力する。
【００３５】
表示部２４は、ＬＣＤ（Liquid Crystal Display）やＣＲＴ（Cathode Ray Tube）等のモニターにより構成され、制御部２１から入力される表示信号の指示に従って、操作部２３からの入力指示やデータ等を表示する。
【００３６】
通信部２５は、ＬＡＮアダプターやモデムやＴＡ（Terminal Adapter）等を備え、通信ネットワークＮＴに接続された各装置との間のデータ送受信を制御する。
【００３７】
〔診断用コンソール３の構成〕
診断用コンソール３は、撮影用コンソール２から動態画像を取得し、取得した動態画像を解析し、その解析結果画像を表示して医師の診断を支援するための動態解析装置である。
診断用コンソール３は、図１に示すように、制御部３１、記憶部３２、操作部３３、表示部３４、通信部３５を備えて構成され、各部はバス３６により接続されている。
【００３８】
制御部３１は、ＣＰＵ、ＲＡＭ等により構成される。制御部３１のＣＰＵは、操作部３３の操作に応じて、記憶部３２に記憶されているシステムプログラムや、各種処理プログラムを読み出してＲＡＭ内に展開し、展開されたプログラムに従って、後述する画像解析処理を始めとする各種処理を実行し、診断用コンソール３各部の動作を集中制御する。制御部３１は、取得手段、決定手段、補正手段として機能する。
【００３９】
記憶部３２は、不揮発性の半導体メモリーやハードディスク等により構成される。記憶部３２は、制御部３１で画像解析処理を実行するためのプログラムを始めとする各種プログラムやプログラムにより処理の実行に必要なパラメーター、或いは処理結果等のデータを記憶する。これらの各種プログラムは、読取可能なプログラムコードの形態で格納され、制御部３１は、当該プログラムコードに従った動作を逐次実行する。
【００４０】
また、記憶部３２は、注目領域の移動量Ｍと画素信号値の強調度合い（強調度）の関係を示す強調度テーブル３２１（図５参照）を記憶している。
【００４１】
操作部３３は、カーソルキー、数字入力キー、及び各種機能キー等を備えたキーボードと、マウス等のポインティングデバイスを備えて構成され、キーボードに対するキー操作やマウス操作により入力された指示信号を制御部３１に出力する。また、操作部３３は、表示部３４の表示画面にタッチパネルを備えても良く、この場合、タッチパネルを介して入力された指示信号を制御部３１に出力する。
【００４２】
表示部３４は、ＬＣＤやＣＲＴ等のモニターにより構成され、制御部３１から入力される表示信号の指示に従って、各種表示を行う。
【００４３】
通信部３５は、ＬＡＮアダプターやモデムやＴＡ等を備え、通信ネットワークＮＴに接続された各装置との間のデータ送受信を制御する。
【００４４】
〔動態解析システム１００の動作〕
次に、上記動態解析システム１００における動作について説明する。
【００４５】
（撮影装置１、撮影用コンソール２の動作）
まず、撮影装置１、撮影用コンソール２による撮影動作について説明する。
図２に、撮影用コンソール２の制御部２１において実行される撮影制御処理を示す。撮影制御処理は、制御部２１と記憶部２２に記憶されているプログラムとの協働により実行される。
【００４６】
まず、撮影実施者により撮影用コンソール２の操作部２３が操作され、被写体Ｍとなる被検者の患者情報（例えば、患者ＩＤ、氏名、身長、体重、性別等）、検査情報（検査対象部位（ここでは、胸部）、診断対象の種類（例えば、換気、血流）等）の入力が行われる（ステップＳ１）。なお、患者情報や検査情報は、通信部２５を介してＲＩＳ（Radiology Information System）から受信する構成としてもよい。
【００４７】
次いで、入力された検査情報に基づいて、放射線照射条件が記憶部２２から読み出されて放射線照射制御装置１２に設定されるとともに、画像読取条件が記憶部２２から読み出されて読取制御装置１４に設定される（ステップＳ２）。なお、放射線照射条件及び画像読取条件は、操作部２３により撮影実施者が調整し、設定し直すことができる。
【００４８】
次いで、操作部２３の操作による放射線照射の指示が待機される（ステップＳ３）。ここで、撮影実施者は、被写体Ｍを放射線源１１と放射線検出部１３の間に配置してポジショニングを行う。また、本実施形態においては、呼吸状態下で撮影を行うため、被検者（被写体Ｍ）に楽にするように指示し、安静呼吸を促す。或いは、「吸って、吐いて」等の深呼吸の誘導を行うこととしてもよい。撮影準備が整った時点で、操作部２３を操作して放射線照射指示を入力する。
【００４９】
操作部２３により放射線照射指示が入力されると（ステップＳ３；ＹＥＳ）、放射線照射制御装置１２及び読取制御装置１４に撮影開始指示が出力され、動態撮影が開始される（ステップＳ４）。即ち、放射線照射制御装置１２に設定されたパルス間隔で放射線源１１により放射線が照射され、放射線検出部１３によりフレーム画像が取得される。
【００５０】
予め定められたフレーム数の撮影が終了すると、制御部２１により放射線照射制御装置１２及び読取制御装置１４に撮影終了の指示が出力され、撮影動作が停止される。撮影されるフレーム数は、少なくとも１呼吸サイクルが撮影できる枚数である。
【００５１】
撮影により取得されたフレーム画像は順次撮影用コンソール２に入力され、撮影順を示す番号（フレーム番号）と対応付けて記憶部２２に記憶されるとともに（ステップＳ５）、表示部２４に表示される（ステップＳ６）。撮影実施者は、表示された動態画像によりポジショニング等を確認し、撮影により診断に適した画像が取得された（撮影ＯＫ）か、再撮影が必要（撮影ＮＧ）か、を判断する。そして、操作部２３を操作して、判断結果を入力する。
【００５２】
操作部２３の所定の操作により撮影ＯＫを示す判断結果が入力されると（ステップＳ７；ＹＥＳ）、動態撮影で取得された一連のフレーム画像のそれぞれに、動態画像を識別するための識別ＩＤや、患者情報、検査情報、放射線照射条件、画像読取条件、撮影順を示す番号（フレーム番号）等の情報が付帯され（例えば、画像データのヘッダ領域に書き込まれ）、通信部２５を介して診断用コンソール３に送信される（ステップＳ８）。そして、本処理は終了する。一方、操作部２３の所定の操作により撮影ＮＧを示す判断結果が入力されると（ステップＳ７；ＮＯ）、記憶部２２に記憶された一連のフレーム画像が削除され（ステップＳ９）、本処理は終了する。この場合、再撮影が必要となる。
【００５３】
（診断用コンソール３の動作）
次に、診断用コンソール３における動作について説明する。
診断用コンソール３においては、通信部３５を介して撮影用コンソール２から動態画像の一連のフレーム画像が受信されると、制御部３１と記憶部３２に記憶されているプログラムとの協働により図３に示す画像解析処理が実行される。
【００５４】
以下、図３を参照して画像解析処理の流れについて説明する。
まず、受信した動態画像に前処理が施される（ステップＳ１０）。
前処理では、診断対象の種類（換気、血流）に応じた信号成分（換気信号成分、血流信号成分）を抽出する処理が施される。
まず、動態画像の各フレーム画像が診断対象の種類に応じた所定サイズの小領域に分割され、分割された小領域内において、画素信号値の代表値（例えば、平均値、中央値等）が算出され、小領域内の画素信号値が算出した代表値に置き換えられる。これにより、診断対象に応じた空間周波数の信号成分を抽出することができる。次いで、一連のフレーム画像間で各小領域の対応付けが行われ、小領域毎に、領域内の画素信号値の時間変化が診断対象の種類に応じたカットオフ周波数の時間方向のハイパスフィルタ、ローパスフィルタ、またはバンドパスフィルタでフィルタリングされる。これにより、診断対象の種類に応じた時間周波数の信号成分を抽出することができる。
【００５５】
次いで、前処理後の動態画像に解析が行われ、解析結果画像が生成される（ステップＳ１１）。
例えば、動態画像の隣接するフレーム画像間で対応する画素毎に画素信号値の差分をとる（例えば、ｎ番目のフレーム画像の各画素の画素信号値からｎ−１番目（或いは、所定フレーム前）のフレーム画像の対応する各画素の画素信号値を引く）ことにより解析結果画像（フレーム間差分画像）が生成される。或いは、動態画像の各フレーム画像と基準のフレーム画像との間で対応する画素毎に画素信号値の差分をとることにより解析結果画像を生成することとしてもよい。なお、解析結果画像の種類は、特に限定されない。また、ステップＳ１０やＳ１１を除いてステップＳ１２に進んでも良い。
【００５６】
次いで、画素信号値補正処理が実行される（ステップＳ１２）。
【００５７】
図４に、第１の実施形態において図３のステップＳ１２で実行される画素信号値補正処理Ａのフローチャートを示す。
画素信号値補正処理Ａにおいては、まず、前処理後の動態画像において、各画素に対応する注目領域が設定され、設定された注目領域の移動情報（移動量Ｍ）が取得される（ステップＳ１２０１）。
【００５８】
ステップＳ１２０１においては、例えば、一のフレーム画像（例えば、１番目のフレーム画像）の左上の画素から順に、その画素（注目画素）を中心とする注目領域（例えば、ｍ×ｍ画素（ｍは正の整数））が設定され、設定された注目領域が他の各フレーム画像のどこに対応するかがテンプレートマッチングにより特定される。そして、例えば、ｎ−１番目の注目領域の中心からｎ番目のフレーム画像の注目領域の中心までの移動量がｎ番目のフレーム画像のその注目領域の移動情報（移動量Ｍ）として取得される。なお、ここでは隣り合うフレーム画像間で移動量Ｍを求める例について説明したが、所定間隔離れたフレーム画像との移動量Ｍを求めることとしてもよい。
【００５９】
なお、本実施形態においては、動態画像の各画素に対応する注目領域を設定することとしているが、注目領域は、少なくとも診断対象となる検査対象部位の領域内の各画素に対応して設定すればよい。例えば、本実施形態では、検査対象部位が肺野であるため、肺野領域を抽出し、肺野領域内にのみ注目領域を設定しても構わない。
【００６０】
次いで、取得された各注目領域の移動情報に基づいて、各画素の画素信号値の強調度αが決定される（ステップＳ１２０２）。
例えば、記憶部３２に記憶されている強調度テーブル３２１が参照され、注目領域の移動量Ｍに応じた強調度αがその注目領域に対応する注目画素の画素信号値の強調度αとして決定される。図５に、強調度テーブル３２１の一例を示す。例えば、移動量Ｍが図５に矢印で示す値である場合には、０．５が注目画素の画素信号値の強調度αとして決定される。なお、図５の強調度テーブル３２１は、移動量Ｍに対する強調度αをグラフとして示しているが、移動量Ｍの値に対する強調度αの値を対応付けた一覧表（テーブル）としてもよい。
【００６１】
次いで、取得された各画素の強調度αに基づいて、ステップＳ１１で生成された解析結果画像の画素信号値が補正され（ステップＳ１２０３）、画像解析処理は終了する。
ステップＳ１２０３においては、画素（ｘ１、ｙ１）の画素信号値ｉｎ（ｘ１、ｙ１）は、以下の（式１）により画素信号値ｏｕｔ（ｘ１、ｙ１）に補正される。
ｏｕｔ（ｘ１、ｙ１）＝α×ｉｎ（ｘ１、ｙ１）・・・（式１）
【００６２】
ここで、胸部の動態画像には呼吸や心拍に伴う骨格等の動きが含まれており、これらの動きによるアーチファクトは、例えば、フレーム間差分等の解析によって強調され、注目すべき信号の観察の邪魔になり、解析精度及び診断精度を低下させている。そこで、図５に示すように、移動量Ｍが大きいほど低い強調度αを画素信号値に乗算することで、従来のように処理時間のかかる肋骨陰影低減処理等を施すことなく、短時間で、精度良く、解析結果画像から、病変の観察の邪魔になる肋骨等の動きの大きい構造物を減弱することが可能となる。
【００６３】
画像解析処理において生成された補正済みの解析結果画像は、表示部３４に表示される。ここで表示される解析結果画像は、病変の観察の邪魔になる肋骨等の動きの大きい構造物が減弱された画像であるので、医師は容易に解析結果画像を観察することが可能となり、診断性能を向上させることができる。
【００６４】
＜第２の実施形態＞
次に、本発明の第２の実施形態について説明する。
呼吸運動の呼気位相では、横隔膜や肺組織は上昇し、肋骨、鎖骨、肩等の骨格は下降する。吸気位相では、横隔膜や肺組織は下降し、肋骨、鎖骨、肩等の骨格は上昇する。そこで、第２の実施形態では、横隔膜と逆の方向に移動をする領域ほど画素信号値の強調度αを低くし、或いは、肩と同じ方向に移動する領域ほど画素信号値の強調度αを低くすることで、肋骨等の診断に邪魔な構造物を減弱する。
【００６５】
第２の実施形態において、診断用コンソール３の記憶部３２には、注目領域の移動量Ｍ及び異方度φと画素信号値の強調度αとの関係を示す強調度テーブル３２２（図６参照）、又は、注目領域の移動量Ｍ及び同方向度θと画素信号値の強調度αとの関係を示す強調度テーブル３２３（図７参照）が記憶されている。動態解析システム１００のその他の構成及び撮影動作は、第１の実施形態で説明したものと同様であるので説明を援用し、第２の実施形態の診断用コンソール３の動作について説明する。
【００６６】
第２の実施形態において、診断用コンソール３の制御部３１は、第１の実施形態で説明した図３に示す画像解析処理を実行するが、ステップＳ１２において実行される処理が第１の実施形態と異なるので、以下ステップＳ１２において実行される画素信号値補正処理Ｂについて説明する。
【００６７】
図８に、第２の実施形態において図３のステップＳ１２で実行される画素信号値補正処理Ｂのフローチャートを示す。画素信号値補正処理Ｂは、制御部３１と記憶部３２に記憶されているプログラムとの協働により実行される。制御部３１は、画素信号値補正処理Ｂを実行することにより、取得手段、決定手段、補正手段、第２の取得手段として機能する。
【００６８】
まず、前処理後の動態画像に参照領域が設定される（ステップＳ１２１１）。
参照領域は、例えば、横隔膜、肩（体表）等に設定される。例えば、表示部３４に動態画像の一のフレーム画像（例えば、１フレーム目の画像）が表示され、操作部３３によるユーザ操作により指定された直線又は矩形が参照領域として設定される。他のフレーム画像については、１フレーム目のフレーム画像で指定された領域と同じ画素位置の領域が参照領域として設定される。
【００６９】
横隔膜領域に参照領域を設定する場合、ユーザは、図９にＲで示すように、横隔膜付近で体軸と略平行になる直線又は矩形を参照領域として指定する。肩領域に参照領域を設定する場合、ユーザは、図１０にＲで示すように、肩の体表ラインに略垂直になるように直線又は矩形を指定する。
なお、参照領域は、ユーザ操作ではなく、制御部３１により自動的に設定することとしてもよい。
【００７０】
次いで、前処理後の動態画像において、参照領域に含まれる構造物の移動情報が取得される（ステップＳ１２１２）。
例えば、動態画像の各フレーム画像の参照領域のそれぞれにおいてエッジ検出により最も強いエッジを検出し、最も強いエッジのフレーム画像間での移動方向を参照領域に含まれる構造物の移動情報として取得することができる。ここで、肩の移動方向は斜め上下方向であるが、図１１に示すように、別の軸（例えば、ｙ軸）に投影したデータを用いてもよい。これにより、斜め方向の移動を上下方向の移動に変換することができるため、後段の処理で同方向度θを算出する際の基準を上下方向とすることができ、呼吸により上下方向に移動する肋骨等の骨部の画素信号値を減弱しやすくなるためである。  
【００７１】
或いは、動態画像の各フレーム画像において、参照領域内の画素信号値の代表値（例えば、平均値、中央値等）を算出し、ｎ番目のフレーム画像から算出された代表値からｎ−１番目のフレーム画像から算出された代表値を引いた差分値に基づいて参照領域に含まれる構造物の移動方向を特定し、特定した移動方向を参照領域に含まれる構造物の移動情報として取得することができる。
【００７２】
横隔膜が上昇した場合、これに連動して肺野領域が上昇するため、横隔膜に設定した参照領域においては、画素信号値の低い肺野領域が占める割合が減少し、画素信号値の高い横隔膜下の白い領域の割合が増加して、画素信号値の代表値は上昇する。即ち、参照領域を横隔膜に設定し、参照領域の代表値の差分値が＋である場合、横隔膜の移動方向は上であると特定することができる。一方、横隔膜が下降した場合、これに連動して肺野領域が下降するため、横隔膜に設定した参照領域においては、画素信号値の低い肺野領域の占める割合が増加し、画素信号値の高い横隔膜下の白い領域の割合が減少して、画素信号値の代表値は減少する。即ち、参照領域を横隔膜に設定し、参照領域の代表値の差分値が−である場合、横隔膜の移動方向は下であると特定することができる。
【００７３】
肩が上昇した場合、肩に設定した参照領域においては、素抜け領域（放射線が直接到達した領域）の占める割合が減少し、画素信号値の高い人体領域（白い領域）の占める割合が増加するため、画素信号値の代表値は上昇する。即ち、参照領域を肩に設定し、参照領域の代表値の差分値が＋である場合、肩の移動方向は上方向であると特定することができる。一方、肩が下降した場合、肩に設定した参照領域においては、素抜け領域の割合が増加し、画素信号値の高い人体領域の占める割合が減少するため、画素信号値の代表値は減少する。即ち、参照領域を肩に設定し、参照領域の代表値の差分値が−である場合、肩の移動方向は下方向であると特定することができる。
【００７４】
次いで、前処理後の動態画像において、各画素に対応する注目領域が設定され、設定された注目領域の移動情報（移動量Ｍ及び移動方向）が取得される（ステップＳ１２１３）。
【００７５】
ステップＳ１２１３においては、例えば、基準画像（例えば、１番目のフレーム画像）の対象領域の左上の画素から順に、その画素（注目画素）を中心とする注目領域（例えば、ｍ×ｍ画素（ｍは正の整数））が設定され、設定された注目領域が他の各フレーム画像のどこに対応するかがテンプレートマッチングにより特定される。そして、例えば、ｎ−１番目の注目領域の中心位置からｎ番目のフレーム画像の注目領域の中心位置までの移動量Ｍ及び移動方向をｎ番目のフレーム画像のその注目領域の移動情報として求める。なお、ここでは隣合うフレーム画像間で移動量Ｍ及び移動方向を求める例について説明したが、所定間隔離れたフレーム画像との移動量を求めることとしてもよい。
【００７６】
次いで、取得された注目領域の移動情報及び参照領域に含まれる構造物の移動情報に基づいて、各画素の画素信号値の強調度αが決定される（ステップＳ１２１４）。
【００７７】
例えば、参照領域が横隔膜に設定された場合、まず、図１２に示すように、参照領域に含まれる構造物の移動方向と注目画素（ｘ１、ｙ１）を中心とする注目領域の移動方向の角度差である異方度φ（ｘ１、ｙ１）が算出される。次いで、記憶部３２に記憶されている強調度テーブル３２２（図６参照）が参照され、注目領域の移動量Ｍ（ｘ１、ｙ１）及び異方度φ（ｘ１、ｙ１）に応じた強調度α（ｘ１、ｙ１）がその注目領域に対応する注目画素の画素信号値の強調度αとして決定される。
【００７８】
例えば、参照領域が肩に設定された場合、まず、図１３に示すように、参照領域に含まれる構造物の移動方向と注目画素（ｘ１、ｙ１）を中心とする注目領域の移動方向の角度差φをπから引いて同方向度θ（ｘ１、ｙ１）が算出される。次いで、記憶部３２に記憶されている強調度テーブル３２３（図７参照）が参照され、注目領域の移動量Ｍ（ｘ１、ｙ１）及び同方向度θ（ｘ１、ｙ１）に応じた強調度α（ｘ１、ｙ１）がその注目領域に対応する注目画素の強調度αとして決定される。
なお、参照領域が肩に設定された場合も、横隔膜と同様に、異方度φを算出し、強調度テーブル３２２を用いて強調度αを決定してもよい。ただし、強調度テーブル３２２のグラフ形状は図６に示したものとは異なる。
また、図６の強調度テーブル３２２及び図７の強調度テーブル３２３は、移動量Ｍ及び異方度φ（同方向度θ）に対する強調度αをグラフとして示しているが、移動量Ｍ及び異方度φ（同方向度θ）の値に対する強調度αの値を対応付けた一覧表（テーブル）としてもよい。
【００７９】
次いで、取得された各画素の強調度αに基づいて、ステップＳ１１で生成された解析結果画像の画素信号値が補正され（ステップＳ１２１５）、画像解析処理は終了する。
ステップＳ１２１５においては、注目画素（ｘ１、ｙ１）の画素信号値ｉｎ（ｘ１、ｙ１）は、上述の（式１）により画素信号値ｏｕｔ（ｘ１、ｙ１）に補正される。
【００８０】
ここで、胸部の動態画像には呼吸や心拍に伴う骨格等の動きが含まれており、これらの動きによるアーチファクトは、例えば、フレーム間差分等の解析によって強調され、注目すべき信号の観察の邪魔になり、解析精度及び診断精度を低下させている。肋骨等の骨格は、上述のように、呼吸により、横隔膜と異なる方向に動き、肩と同方向に動く。そこで、図６に示すように、移動量Ｍが大きく、かつ横隔膜との異方度φが大きい（即ち、横隔膜との角度差が大きい）注目領域の注目画素ほど、解析結果画像の画素信号値に低い強調度αを乗算する。或いは、図７に示すように、移動量Ｍが大きく、かつ肩との同方向度θが大きい（即ち、肩との角度差が小さい）注目領域の注目画素ほど、解析結果画像の画素信号値に低い強調度αを乗算する。このように、注目領域の移動量Ｍだけでなく、注目領域の移動方向も加味して強調度αを決定することで、従来のように処理時間のかかる肋骨陰影低減処理等を施すことなく、短時間で、より精度良く病変の観察の邪魔になる肋骨等の構造物を減弱することが可能となる。
【００８１】
画像解析処理において生成された補正済みの解析結果画像は、表示部３４に表示される。ここで表示される解析結果画像は、病変の観察の邪魔になる肋骨等の構造物が減弱された画像であるので、医師は容易に解析結果画像を観察することが可能となり、診断性能を向上させることができる。
【００８２】
＜第３の実施形態＞
次に、本発明の第３の実施形態について説明する。
呼吸運動の呼気位相では、肋骨、鎖骨等の骨部は下降する。また、肺胞の空気量が減少するため肺胞の密度は増加し、肋間の画素信号値は上昇する。一方、呼吸運動の吸気位相では、肋骨、鎖骨等の骨格は上昇する。また、肺胞の空気量が増加するため肺胞の密度は減少し、肋間の画素信号値は減少する。第３の実施形態では、この原理を用いて、肋間の画素信号値に基づいて肋骨の移動方向を特定し、肋骨の移動方向と同方向に移動する領域の画素信号値の強調度αを低くすることで、肋骨等の診断に邪魔な構造物を精度良く減弱する。
【００８３】
第３の実施形態において、診断用コンソール３の記憶部３２には、注目領域の移動量Ｍ及び同方向度θと、画素信号値の強調度αとの関係を示す強調度テーブル３２３が記憶されている。動態解析システム１００のその他の構成及び撮影動作は、第１の実施形態で説明したものと同様であるので説明を援用し、第３の実施形態の診断用コンソール３の動作について説明する。
【００８４】
第３の実施形態において、診断用コンソール３の制御部３１は、第１の実施形態で説明した図３に示す画像解析処理を実行するが、図３のステップＳ１２において実行される処理が第１の実施形態と異なるので、以下ステップＳ１２において実行される画素信号値補正処理Ｃについて説明する。
【００８５】
図１４に、第３の実施形態において図３のステップＳ１２で実行される画素信号値補正処理Ｃのフローチャートを示す。画素信号値補正処理Ｃは、制御部３１と記憶部３２に記憶されているプログラムとの協働により実行される。制御部３１は、画素信号値補正処理Ｃを実行することにより、取得手段、決定手段、補正手段、特定手段として機能する。
【００８６】
まず、前処理後の動態画像に参照領域が設定される（ステップＳ１２２１）。
ここでは、肋間に参照領域が設定される。例えば、表示部３４に動態画像の一のフレーム画像（例えば、１フレーム目の画像）が表示され、図１５にＲで示すように、操作部３３によるユーザ操作により指定された肋骨と肋骨の間の領域（例えば、矩形、円等）が参照領域として設定される。ここで、ユーザは、全フレーム画像において肋骨が進入しない領域を参照領域として指定する必要がある。他のフレーム画像については、１フレーム目のフレーム画像で指定された領域と同じ画素位置の領域が参照領域として設定される。
なお、参照領域は、ユーザ操作ではなく、制御部３１により自動的に設定することとしてもよい。
【００８７】
次いで、前処理後の動態画像の参照領域の画素信号値の変化に基づいて、肋骨の移動方向が特定される（ステップＳ１２２２）。
ステップＳ１２２２においては、まず、前処理後の動態画像の各フレーム画像において、参照領域内の画素信号値の代表値（例えば、平均値、中央値等）が算出され、ｎ番目のフレーム画像から算出された代表値からｎ−１番目のフレーム画像から算出された代表値を引いた差分値が算出される。差分値が＋の場合、上述のように呼気位相に相当し、この場合、肋骨の移動方向は下であると特定される。差分値が−の場合、上述のように吸気位相に相当し、この場合、肋骨の移動方向は上であると特定される。
【００８８】
次いで、前処理後の動態画像において、各画素に対応する注目領域が設定され、設定された注目領域の移動情報（移動量Ｍ及び移動方向）が取得される（ステップＳ１２２３）。
ステップＳ１２２３の処理は、ステップＳ１２１３の処理と同様であるので説明を援用する。
【００８９】
次いで、取得された注目領域の移動情報及び肋骨の移動情報に基づいて、各画素の画素信号値の強調度αが決定される（ステップＳ１２２４）。
まず、肋骨の移動方向と注目画素（ｘ１、ｙ１）を中心とする注目領域の移動方向の角度差φをπから引いて同方向度θ（ｘ１、ｙ１）が算出される。次いで、記憶部３２に記憶されている強調度テーブル３２３が参照され、注目領域の移動量Ｍ（ｘ１、ｙ１）及び同方向度θ（ｘ１、ｙ１）に応じた強調度α（ｘ１、ｙ１）がその注目領域における画素信号値の強調度αとして決定される。
なお、角度差φを異方度φとして、強調度テーブル３２２を用いて強調度αを決定してもよい。ただし、強調度テーブル３２２のグラフ形状は図６に示したものとは異なる。
【００９０】
次いで、取得された各画素の強調度αに基づいて、ステップＳ１１で生成された解析結果画像の各画素の画素信号値が補正され（ステップＳ１２２５）、画像解析処理は終了する。
ステップＳ１２２５においては、注目領域（ｘ１、ｙ１）の画素信号値ｉｎ（ｘ１、ｙ１）は、上述の（式１）により画素信号値ｏｕｔ（ｘ１、ｙ１）に補正される。
【００９１】
ここで、第３の実施形態では、肋骨の影響を受けずに肺の密度変化を示す画素信号値を得ることのできる肋間領域の画素信号値の変化に基づいて肋骨の移動方向を特定し、移動量Ｍが大きく、かつ肋骨との同方向度θが大きい（即ち、肋骨との角度差が小さい）注目領域の注目画素ほど、解析結果画像の画素信号値に低い強調度αを乗算する。このように、注目領域の移動量Ｍだけでなく、移動方向も加味して強調度αを導出することで、従来のように処理時間のかかる肋骨陰影低減処理等を施すことなく、短時間で、より精度良く病変の観察の邪魔になる肋骨等の構造物を減弱することが可能となる。
【００９２】
画像解析処理において生成された補正済みの解析結果画像は、表示部３４に表示される。ここで表示される解析結果画像は、病変の観察の邪魔になる肋骨等の構造物が減弱された画像であるので、医師は容易に解析結果画像を観察することが可能となり、診断性能を向上させることができる。
【００９３】
以上説明したように、診断用コンソール３の制御部３１によれば、胸部の動態画像に注目領域を設定し、注目領域の移動に関する移動情報を取得し、取得した注目領域の移動情報に基づいて、注目領域に対応する注目画素の画素信号値の強調度を決定する。そして、決定した強調度に基づいて、動態画像を解析することにより生成されたフレーム間差分画像等の解析結果画像の注目画素の画素信号値を補正する。
従って、従来のように処理時間のかかる肋骨陰影低減処理等を施すことなく、短時間で、より精度良く病変の観察の邪魔になる肋骨等の構造物を減弱することが可能となる。
【００９４】
例えば、制御部３１は、注目領域の移動量が大きいほど注目画素の画素信号値の強調度を低い強調度に決定する。従って、病変の観察の邪魔になる肋骨等の動きの大きい構造物を減弱することが可能となる。
【００９５】
また、例えば、制御部３１は、動態画像に参照領域を設定し、参照領域に含まれる構造物の移動に関する移動情報を取得し、注目領域の移動情報及び参照領域に含まれる構造物の移動情報に基づいて、注目領域に対応する注目画素の画素信号値の強調度を決定する。従って、病変の観察の邪魔になる肋骨等の構造物をより精度良く減弱することが可能となる。
【００９６】
例えば、制御部３１は、動態画像の横隔膜を含む領域に参照領域を設定して横隔膜の移動方向を取得し、注目領域の移動方向と横隔膜の移動方向の角度差が大きいほど注目画素の画素信号値の強調度を低い強調度に決定する。従って、横隔膜との移動方向の角度差の大きい肋骨等の構造物をより精度良く減弱することができる。
また、例えば、制御部３１は、動態画像の肩を含む領域に参照領域を設定して肩の移動方向を取得し、注目領域の移動方向と肩の移動方向の角度差が小さいほど注目画素の画素信号値の強調度を低い強調度に決定する。従って、肩との移動方向の角度差が小さい肋骨等の構造物をより精度良く減弱することができる。
【００９７】
また、例えば、制御部３１は、動態画像に参照領域を設定し、参照領域の画素信号値の変化に基づいて所定の構造物の移動方向を特定し、注目領域の移動方向と所定の構造物の移動方向の角度差に基づいて注目領域に対応する注目画素の強調度を決定する。例えば、動態画像の肋間領域に参照領域を設定し、参照領域の画素信号値の変化に基づいて肋骨の移動方向を特定し、注目領域の移動方向と肋骨の移動方向の角度差が小さいほど注目画素の画素信号値の強調度を低い強調度に決定する。従って、肋骨や、肋骨との移動方向の角度差が小さい骨格等の構造物をより精度良く減弱することができる。
【００９８】
なお、本実施形態における記述は、本発明に係る好適な動態解析システムの一例であり、これに限定されるものではない。
【００９９】
例えば、上記実施形態においては、動態画像が胸部の動態画像である場合を例にとり説明したが、これに限定されず、他の部位を撮影した動態画像としてもよい。
【０１００】
また、上記実施形態においては、注目領域の移動情報、又は注目領域及び参照領域の移動情報に基づいて、解析結果画像の各画素の画素信号値を補正する場合を例にとり説明したが、これに限定されず、検査対象部位に放射線を照射して動態撮影することにより得られた動態画像の各画素の画素信号値を補正することとしてもよい。これにより、動態画像における病変の観察の邪魔になる肋骨等の動きの大きい構造物を減弱することが可能となる。
また、動態画像の各画素の画素信号値を本発明により補正してから解析を行って解析結果画像を取得することとしてもよい。
なお、本願の請求項における「検査対象部位に放射線を照射して動態撮影することにより得られた動態画像」には、動態撮影により放射線検出部１３から取得された動態画像（所謂ＲＡＷデータ）だけでなく、これに加工（例えば、画素信号値を吸収線量に変換する処理や上述の前処理をはじめとする画像処理）を施した動態画像も含まれる。
【０１０１】
また、上記実施形態においては、前処理を行う場合を例にとり説明したが、これは、精度良く換気や血流の解析結果を得るために行うものであり、本発明の実施においては省略してもかまわない。
【０１０２】
また、第１の実施形態、第２の実施形態、第３の実施形態を組み合わせて用いてもよい。
【０１０３】
また、例えば、上記の説明では、本発明に係るプログラムのコンピュータ読み取り可能な媒体としてハードディスクや半導体の不揮発性メモリー等を使用した例を開示したが、この例に限定されない。その他のコンピュータ読み取り可能な媒体として、ＣＤ-ＲＯＭ
等の可搬型記録媒体を適用することが可能である。また、本発明に係るプログラムのデータを通信回線を介して提供する媒体として、キャリアウエーブ(搬送波)も適用される。
【０１０４】
その他、動態解析システム１００を構成する各装置の細部構成及び細部動作に関しても、本発明の趣旨を逸脱することのない範囲で適宜変更可能である。
【符号の説明】
【０１０５】
１００  画像表示システム
１  撮影装置
１１  放射線源
１２  放射線照射制御装置
１３  放射線検出部
１４  読取制御装置
２  撮影用コンソール
２１  制御部
２２  記憶部
２３  操作部
２４  表示部
２５  通信部
２６  バス
３  診断用コンソール
３１  制御部
３２  記憶部
３３  操作部
３４  表示部
３５  通信部
３６  バス
【図１】
2018007949.tif 000003
【図２】
2018007949.tif 000004
【図３】
2018007949.tif 000005
【図４】
2018007949.tif 000006
【図５】
2018007949.tif 000007
【図６】
2018007949.tif 000008
【図７】
2018007949.tif 000009
【図８】
2018007949.tif 000010
【図９】
2018007949.tif 000011
【図１０】
2018007949.tif 000012
【図１１】
2018007949.tif 000013
【図１２】
2018007949.tif 000014
【図１３】
2018007949.tif 000015
【図１４】
2018007949.tif 000016
【図１５】
2018007949.tif 000017