(57)【特許請求の範囲】
【請求項１】
被写体の胸部の動態撮影を行い、連続する複数のフレーム画像を生成する撮影手段と、
前記生成された複数のフレーム画像のそれぞれから被写体領域を抽出し、当該抽出された被写体領域を複数の領域に分割し、前記複数のフレーム画像間において前記分割された領域を互いに対応付けて解析することにより、前記分割された領域の動きを示す特徴量を算出する第１の解析手段と、
前記被写体領域の肋骨、鎖骨、及び肩のうち少なくとも1つの動きに基づいて、前記被写体領域全体の動きを示す特徴量を算出する第２の解析手段と、
前記第１の解析手段及び前記第２の解析手段により算出された特徴量を同時に表示する表示手段と、
を備える診断支援システム。
【請求項２】
前記第２の解析手段は、肋骨の動きに基づいて、胸式呼吸の度合い又は胸郭の動きの制限度合いを示す特徴量を算出する請求項１に記載の診断支援システム。
【請求項３】
前記第２の解析手段は、鎖骨又は肩の動きに基づいて、呼吸補助筋の動きの制限度合いを示す特徴量を算出する請求項１に記載の診断支援システム。
【請求項４】
前記分割された領域の動きを示す特徴量は、当該領域の換気又は血流を示す特徴量である請求項１〜３の何れか一項に記載の診断支援システム。
【請求項５】
前記第１の解析手段は、前記複数のフレーム画像の前記分割された領域内の画素信号値の平均値を算出し、当該算出された平均信号値の時間変化をその領域の動きを示す特徴量として算出し、
前記表示手段は、前記平均信号値の時間変化を波形として表示する請求項１〜４の何れか一項に記載の診断支援システム。
【請求項６】
前記第１の解析手段による解析結果の表示対象となる領域をユーザが選択するための操作手段を有し、
前記第１の解析手段は、前記操作手段により選択された領域の動きを示す特徴量を算出し、
前記表示手段は、前記選択された領域の動きを示す特徴量を表示する請求項１〜５の何れか一項に記載の診断支援システム。
【請求項７】
前記表示手段は、前記操作手段により肺野領域の一方の領域が選択された場合に、前記選択された領域の特徴量とともに、前記選択された領域と体軸対称の他方の肺野の領域の特徴量を同時に表示する請求項６に記載の診断支援システム。
【請求項８】
コンピュータを、
被写体の胸部の動態を示す複数のフレーム画像のそれぞれから被写体領域を抽出し、当該抽出された被写体領域を複数の領域に分割し、前記複数のフレーム画像間において前記分割された領域を互いに対応付けて解析することにより、前記分割された領域の動きを示す特徴量を算出する第１の解析手段、
前記被写体領域の肋骨、鎖骨、及び肩のうち少なくとも1つの動きに基づいて、前記被写体領域全体の動きを示す特徴量を算出する第２の解析手段、
前記第１の解析手段及び前記第２の解析手段により算出された特徴量を同時に表示する表示手段、
として機能させるためのプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、診断支援システム及びプログラムに関する。
【背景技術】
【０００２】
従来、医師は、胸部Ｘ線写真（深呼吸時の肺野面積が最大となったときに撮影した画像）による画像情報及び聴診器の音情報という、目と耳からの両情報により総合的に肺の換気機能についての診断を行っていた。【０００３】
これに対し、近年では、ＦＰＤ（flat panel detector）等の半導体イメージセンサを利用して胸部の動態画像を撮影し、診断に応用する試みがなされるようになってきている。例えば、特許文献１には、動態画像を構成する複数のフレーム画像間で信号値の差を示す差分画像を生成し、生成された差分画像から各信号値の最大値を求めて表示する技術が記載されている。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特開２００４−３１２４３４号公報
【発明の概要】
【発明が解決しようとする課題】
【０００５】
しかしながら、聴診器の音（耳）による診断は、胸部Ｘ線写真のように教示画像に基づくカンファレンス等での診断能の普及、複数の医師による同時確認が困難で、個々の医師の熟練度合いに依存している現状がある。また、聴診器での医師との対面では、患者の緊張により肺が通常時とは異なる挙動を示すことがあり、この影響の判断も含めた熟練が求められる。【０００６】
一方、特許文献１の技術においては、単に動態画像の画素毎のフレーム間差分値の最大値を算出して表示するだけであり、この結果からだけでは、熟練した呼吸器専門医以外の医師が病態を把握することは難しい。【０００７】
本発明の課題は、動態撮影された画像を有効に活用し、診断に使用する情報を目視系に統一し、聴診器の経験の浅い医師であっても精度良い診断を行うことが可能なＧＵＩを含むシステムを提供することである。
【課題を解決するための手段】
【０００８】
上記課題を解決するため、
診断支援システムは、
被写体の胸部の動態撮影を行い、連続する複数のフレーム画像を生成する撮影手段と、
前記生成された複数のフレーム画像のそれぞれから被写体領域を抽出し、当該抽出された被写体領域を複数の領域に分割し、前記複数のフレーム画像間において前記分割された領域を互いに対応付けて解析することにより、前記分割された領域の動きを示す特徴量を算出する第１の解析手段と、
前記被写体領域の肋骨、鎖骨、及び肩のうち少なくとも1つの動きに基づいて、前記被写体領域全体の動きを示す特徴量を算出する第２の解析手段と、
前記第１の解析手段及び前記第２の解析手段により算出された特徴量を同時に表示する表示手段と、
を備える。
【０００９】
前記第２の解析手段は、肋骨の動きに基づいて、胸式呼吸の度合い又は胸郭の動きの制限度合いを示す特徴量を算出することが好ましい。
【００１０】
前記第２の解析手段は、鎖骨又は肩の動きに基づいて、呼吸補助筋の動きの制限度合いを示す特徴量を算出することが好ましい。
【００１１】
前記分割された領域の動きを示す特徴量は、当該領域の換気又は血流を示す特徴量であることが好ましい。
【００１２】
前記第１の解析手段は、前記複数のフレーム画像の前記分割された領域内の画素信号値の平均値を算出し、当該算出された平均信号値の時間変化をその領域の動きを示す特徴量として算出し、
前記表示手段は、前記平均信号値の時間変化を波形として表示することが好ましい。
【００１３】
前記第１の解析手段による解析結果の表示対象となる領域をユーザが選択するための操作手段を有し、
前記第１の解析手段は、前記操作手段により選択された領域の動きを示す特徴量を算出し、
前記表示手段は、前記選択された領域の動きを示す特徴量を表示することが好ましい。
【００１４】
前記表示手段は、前記操作手段により肺野領域の一方の領域が選択された場合に、前記選択された領域の特徴量とともに、前記選択された領域と体軸対称の他方の肺野の領域の特徴量を同時に表示することが好ましい。
【００１８】
本発明の第２の側面によると、
プログラムは、
コンピュータを、
被写体の胸部の動態を示す複数のフレーム画像のそれぞれから被写体領域を抽出し、当該抽出された被写体領域を複数の領域に分割し、前記複数のフレーム画像間において前記分割された領域を互いに対応付けて解析することにより、前記分割された領域の動きを示す特徴量を算出する第１の解析手段、
前記被写体領域の肋骨、鎖骨、及び肩のうち少なくとも1つの動きに基づいて、前記被写体領域全体の動きを示す特徴量を算出する第２の解析手段、
前記第１の解析手段及び前記第２の解析手段により算出された特徴量を同時に表示する表示手段、
として機能させる。
【発明の効果】
【００１９】
本発明によれば、動態撮影された画像を有効に活用し、診断に使用する情報を目視系に統一し、聴診器の経験の浅い医師であっても精度良い診断を行うことが可能なＧＵＩを含むシステムを提供することが可能となる。
【図面の簡単な説明】
【００２０】
【図１】本発明の実施の形態における診断支援システムの全体構成を示す図である。
【図２】図１の放射線発生装置と放射線検出部の詳細構成例を示す図である。
【図３】図１の撮影用コンソールの制御部により実行される撮影制御処理を示すフローチャートである。
【図４】図１の診断用コンソールの制御部により実行される画像解析処理を示すフローチャートである。
【図５】図４のステップＳ１１において実行される前処理を示すフローチャートである。
【図６】図４のステップＳ１２において実行される画像解析・表示処理Ａを示すフローチャートである。
【図７】１つの呼吸サイクル（深呼吸時）において撮影された複数の時間位相Ｔ（Ｔ＝ｔ０〜ｔ６）のフレーム画像を示す図である。
【図８】横隔膜位置の算出を説明するための図である。
【図９】胸郭の幅の算出を説明するための図である。
【図１０】肺野領域の小領域への分割を説明するための図である。
【図１１】安静呼気位と安静吸気位において肺野の同一部分を描画した領域の位置の変化を示す図である。
【図１２】最大呼気位と最大吸気位において肺野の同一部分を描画した領域の位置の変化を示す図である。
【図１３】Ａは、ある肺野領域における正常部分の小領域の換気量の時間変化を示す波形の一例を示す図、Ｂは、拘束性部分の小領域の換気量の時間変化を示す波形の一例を示す図、Ｃは、閉塞性部分の小領域の換気量の時間変化を示す波形の一例を示す図、Ｄは、横隔膜位置の時間変化を示す波形の一例を示す図である。
【図１４】換気の解析結果の表示（標準）が設定されている場合に図６のステップＳ２０９において表示部に表示される表示画面の一例を示す図である。
【図１５Ａ】正常例の波形を選択領域の波形と並列表示させた例を示す図である。
【図１５Ｂ】正常例の波形の周期を患者の波形の周期に一致させて表示させた例を示す図である。
【図１６】換気と血流の解析結果の表示（標準）が設定されている場合に図６のステップＳ２０９において表示部に表示される表示画面の一例を示す図である。
【図１７】追加情報として肺野内の各小領域の解析結果の静止画上への一覧表示が設定されている場合に、図６のステップＳ２０９において表示部に表示される表示画面の一例を示す図である。
【図１８】追加情報として波形の解析結果の動画表示が設定されている場合に、図６のステップＳ２０９において表示部に表示される表示画面の一例を示す図である。
【図１９】解析結果の表示対象領域を予め定められた領域（自動）とすることが選択されている場合に表示部に表示される表示画面の一例を示す図である。
【図２０】第２の実施の形態において図１の診断用コンソールの制御部により実行される画像解析・表示処理Ｂを示すフローチャートである。
【図２１】図２０のステップＳ３０９において表示部に表示される表示画面の一例を示す図である。
【発明を実施するための形態】
【００２１】
以下、図面を参照して本発明に係る実施の形態を詳細に説明する。ただし、発明の範囲は、図示例に限定されない。
【００２２】
〈第１の実施の形態〉 
〔診断支援システム１００の構成〕
まず、構成を説明する。
図１に、本実施の形態における診断支援システム１００の全体構成を示す。
図１に示すように、診断支援システム１００は、撮影装置１と、撮影用コンソール２とが通信ケーブル等により接続され、撮影用コンソール２と、診断用コンソール３とがＬＡＮ（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）等を備え、通信ネットワークＮＴに接続された各装置との間のデータ送受信を制御する。
【００３５】
〔診断用コンソール３の構成〕
診断用コンソール３は、撮影用コンソール２から動態画像を取得し、取得した動態画像を表示して医師が読影診断するための動態画像処理装置である。
診断用コンソール３は、図１に示すように、制御部３１、記憶部３２、操作部３３、表示部３４、通信部３５を備えて構成され、各部はバス３６により接続されている。
【００３６】
制御部３１は、ＣＰＵ、ＲＡＭ等により構成される。制御部３１のＣＰＵは、操作部３３の操作に応じて、記憶部３２に記憶されているシステムプログラムや、各種処理プログラムを読み出してＲＡＭ内に展開し、展開されたプログラムに従って、後述する画像解析処理を始めとする各種処理を実行し、診断用コンソール３各部の動作を集中制御する。
【００３７】
記憶部３２は、不揮発性の半導体メモリやハードディスク等により構成される。記憶部３２は、制御部３１で画像解析処理を実行するための画像解析処理プログラムを始めとする各種プログラムやプログラムにより処理の実行に必要なパラメータ、或いは処理結果等のデータを記憶する。これらの各種プログラムは、読取可能なプログラムコードの形態で格納され、制御部３１は、当該プログラムコードに従った動作を逐次実行する。例えば、記憶部３２は、操作部３３により入力された表示態様の設定情報（動態画像の解析結果の表示態様の設定情報）を記憶している。
【００３８】
操作部３３は、カーソルキー、数字入力キー、及び各種機能キー等を備えたキーボードと、マウス等のポインティングデバイスを備えて構成され、キーボードに対するキー操作やマウス操作により入力された指示信号を制御部３１に出力する。また、操作部３３は、表示部３４の表示画面にタッチパネルを備えても良く、この場合、タッチパネルを介して入力された指示信号を制御部３１に出力する。
【００３９】
本実施の形態において、操作部３３は、後述する画像解析処理における解析結果の表示態様を設定可能に構成されている。表示態様としては、例えば、（ａ）解析結果の表示対象領域をユーザが選択するか、予め定められた領域（自動）とするか（ｂ）換気についての解析結果のみを表示するか、換気とともに血流についての解析結果を表示するか、（ｃ）追加情報として、肺野内の各小領域の解析結果を併せて静止画上に一覧表示するか否か、（ｄ）追加情報として、波形の解析結果の動画表示を行うか否か、（ｄ）複数の特徴量の中から何れの特徴量を解析対象とするか、等を設定可能に構成されている。
【００４０】
表示部３４は、ＬＣＤやＣＲＴ等のモニタにより構成され、制御部３１から入力される表示信号の指示に従って、操作部３３からの入力指示やデータ等を表示する。
【００４１】
通信部３５は、ＬＡＮアダプタやモデムやＴＡ等を備え、通信ネットワークＮＴに接続された各装置との間のデータ送受信を制御する。
【００４２】
〔診断支援システム１００の動作〕
次に、上記診断支援システム１００における動作について説明する。
【００４３】
（撮影装置１、撮影用コンソール２の動作）
まず、撮影装置１、撮影用コンソール２による撮影動作について説明する。
図３に、撮影用コンソール２の制御部２１において実行される撮影制御処理を示す。撮影制御処理は、制御部２１と記憶部２２に記憶されている撮影制御処理プログラムとの協働により実行される。
【００４４】
まず、撮影技師により撮影用コンソール２の操作部２３が操作され、撮影対象（被写体Ｍ）の患者情報（患者の氏名、身長、体重、年齢、性別等）の入力が行われる（ステップＳ１）。      
【００４５】
次いで、放射線照射条件が記憶部２２から読み出されて放射線照射制御装置１２に設定されるとともに、画像読取条件が記憶部２２から読み出されて読取制御装置１４に設定される（ステップＳ２）。ここで、フレームレート（パルスレート）としては、人間の心拍周期を考慮して７．５フレーム／秒以上とすることが好ましい。
【００４６】
次いで、操作部２３の操作による放射線照射の指示が待機される。
撮影技師は、撮影装置１において患者のポジショニング等の撮影準備を行う。具体的には、被写体Ｍ（患者）の背丈に応じて、図示しないフットスイッチにより、放射線検出器１３１が装着された検出器保持部１３２の高さを調整する。また、被写体Ｍに体動補正用のＸ線非透過マーカを２箇所以上（ここでは、マーカＭ１、マーカＭ２の２箇所）貼付する。また、安静呼吸時の動態画像を得るため、患者に楽にするよう指示する。
【００４７】
撮影装置１においては、撮影技師の操作により検出器保持部１３２の高さが調整されると、測距センサＳＥ１により床〜放射線検出器１３１中央までの距離が取得され、読取制御装置１４に出力される。読取制御装置１４においては、測距センサＳＥ１の出力値が高さ情報として放射線照射制御装置１２に出力される。放射線照射制御装置１２においては、測距センサＳＥ２から出力される床〜放射線源１１１の焦点位置までの距離の値が読取制御装置１４から出力された値と同じになるように図示しない駆動機構が駆動され、放射線源保持部１１２の高さが調整される。
被写体Ｍのポジショニングが終了すると、撮影技師は、撮影用コンソール２の操作部２３により放射線照射指示を入力する。
【００４８】
操作部２３により放射線照射指示が入力されると（ステップＳ３；ＹＥＳ）、放射線照射制御装置１２及び読取制御装置１４に撮影開始指示が出力され、動態撮影が開始される（ステップＳ４）。即ち、放射線照射制御装置１２に設定されたパルス間隔で放射線源１１１により放射線が照射され、放射線検出器１３１によりフレーム画像が取得される。動態撮影開始から予め定められた時間が経過すると、制御部２１により放射線照射制御装置１２及び読取制御装置１４に撮影終了の指示が出力され、撮影動作が停止される。なお、動態撮影開始から撮影動作停止までの時間は、複数の呼吸サイクルの動態を十分に撮影可能な時間である。
【００４９】
撮影により取得されたフレーム画像は順次撮影用コンソール２に入力され、各フレーム画像に対して補正処理が行われる（ステップＳ５）。ステップＳ５の補正処理においては、オフセット補正処理、ゲイン補正処理、欠陥画素補正処理の３つの補正処理が行われる。まず最初に、取得された各フレーム画像に対してオフセット補正処理が行われ、取得された各フレーム画像に重畳された暗電流に起因するオフセット値が除去される。オフセット補正処理では、例えば、取得した各フレーム画像の各画素値（濃度値。以下、信号値という）から、予め記憶されたオフセット補正係数を減算する処理が行われる。ここで、オフセット補正係数は、予め放射線非照射時に取得した複数のフレーム画像を平均化した画像である。次いで、ゲイン補正処理が行われ、各フレーム画像の各画素に対応する各検出素子の個体差や読み出しアンプのゲインムラによって生じる画素毎のばらつきが除去される。ゲイン補正処理では、例えば、オフセット補正後の各フレーム画像に、予め記憶されたゲイン補正係数を乗算する処理が行われる。ここで、ゲイン補正係数は、放射線検出器１３１に一様に放射線を照射した時に取得した複数のオフセット補正済みフレーム画像を平均化した画像と、このときの放射線照射条件で期待される出力信号値の関係から、補正後の各画素の信号値が一様となるように予め算出され、記憶された係数である。次いで、欠陥画素補正処理が行われ、周囲の画素と比較して感度が非線形な画素や、感度がない欠落画素が除去される。欠陥画素補正処理では、例えば、予め記憶された欠陥画素位置情報マップに従って、欠陥画素位置情報マップに登録された各欠陥画素において、欠陥画素の信号値をその近傍の欠陥でない画素の信号値の平均値で置き換える処理が行われる。ここで、欠陥画素位置情報マップは、放射線検出器１３１に一様に放射線を照射した時に取得したオフセット補正、ゲイン補正済みのフレーム画像から、予め複数の欠陥画素が認識され、その欠陥画素の位置が登録されたマップである。上記オフセット補正係数及びゲイン補正係数、欠陥画素位置情報マップは、ビニングやダイナミックレンジ等の収集モードに応じて、予め、それぞれ最適な値が記憶されており、それぞれの収集モードにおいて対応する最適な値が読み出されるようになっている。
【００５０】
次いで、補正処理後の各フレーム画像と、撮影順を示す番号と対応付けて記憶部２２に記憶されるとともに（ステップＳ６）、表示部２４に表示される（ステップＳ７）。ここで、各フレーム画像を記憶する直前に、各フレーム画像の各画素の信号値を真数から対数に変換する対数変換処理を行ってから記憶しても良い。撮影技師は、表示された動態画像によりポジショニング等を確認し、撮影により診断に適した画像が取得された（撮影ＯＫ）か、再撮影が必要（撮影ＮＧ）か、を判断する。そして、操作部２３を操作して、判断結果を入力する。尚、撮影により取得された各フレーム画像は、全撮影の終了後に纏めて入力するようにしても良い。
【００５１】
操作部２３の所定の操作により撮影ＯＫを示す判断結果が入力されると（ステップＳ８；ＹＥＳ）、動態撮影で取得された一連のフレーム画像のそれぞれに、動態画像を識別するための識別ＩＤや、患者情報、検査対象部位、放射線照射条件、画像読取条件、撮影順を示す番号、撮影日時等の情報が付帯され（例えば、ＤＩＣＯＭ形式で画像データのヘッダ領域に書き込まれ）、通信部２５を介して診断用コンソール３に送信される（ステップＳ９）。そして、本処理は終了する。一方、操作部２３の所定の操作により撮影ＮＧを示す判断結果が入力されると（ステップＳ８；ＮＯ）、記憶部２２に記憶された一連のフレーム画像が削除され（ステップＳ１０）、本処理は終了する。尚、このケースに於いては、再撮影が実行されることとなる。
【００５２】
（診断用コンソール３の動作）
次に、診断用コンソール３における動作について説明する。
診断用コンソール３においては、通信部３５を介して撮影用コンソール２から動態画像の一連のフレーム画像が受信されると、制御部３１と記憶部３２に記憶されている画像解析処理プログラムとの協働により図４に示す画像解析処理が実行される。
【００５３】
画像解析処理においては、まず、前処理が行われる（ステップＳ１１）。
【００５４】
図５に、ステップＳ１１において実行される前処理のフローチャートを示す。前処理は、制御部３１と記憶部３２に記憶されているプログラムとの協働により実行される。
前処理においては、まず、対数変換処理が行われ、動態画像の各フレーム画像の各画素の信号値が真数から対数に変換される（ステップＳ１０１）。尚、対数変換済の信号が放射線検出器１３１や撮影用コンソール２から出力される場合には、当該ステップは省略される。
【００５５】
次いで、時間変化信号（濃度）補正処理（トレンド補正処理）が行われ、各フレーム画像の直接Ｘ線領域の信号値が同じ値となるように補正される（ステップＳ１０２）。
ステップＳ１０２においては、まず、撮影用コンソール２から入力された一連のフレーム画像の中から任意の基準フレーム画像を選択し、以下の（式１）により各フレーム画像の補正値を算出する。
（式１）各フレーム画像の補正値＝（各フレーム画像の直接Ｘ線領域の平均信号値）−（基準フレーム画像の直接Ｘ線領域の平均信号値）
次いで、各フレーム画像において、画素毎に、信号値から（式１）により算出された補正値が減算される。
【００５６】
次いで、各フレーム画像に対し、グリッド目除去処理が行われる（ステップＳ１０３）。グリッド目除去処理は、散乱放射線を除去するために被写体Ｍ（患者）と放射線検出器１３１の間に設けられた散乱線除去グリッドのグリッド配列に起因する縞模様を除去する処理である。グリッド目除去処理は、公知の技術を用いて行うことができる。例えば、各フレーム画像に離散的フーリエ変換等の周波数変換処理を施した後、ローパスフィルタ処理を行ってグリッド像の周波数を含む高周波数領域を除去し、逆フーリエ変換処理を施すことにより行うことができる（石田隆行著「医用画像処理入門」３．４Ｘ線画像のグリッドによる縦縞影の除去参照）。尚、散乱線除去グリッドが移動グリッド(Ｘ線照射時にグリッドをグリッド面と水平な方向に移動させる機構を持つ)の場合には、当該ステップは省略しても良い。
【００５７】
次いで、各フレーム画像に対し、低周波数ムラ補正処理が実行される（ステップＳ１０４）。ここで、記憶部３２には、ムラ補正データが予め記憶されている。ムラ補正データは、撮影装置１において被写体なしでＸ線を照射することにより取得された撮影画像に基づいて作成された、各画素の信号値の補正値が記憶されたデータである。低周波数ムラ補正処理では、記憶部３２からムラ補正データが読み出され、各フレーム画像の各画素の信号値からその画素の補正値が差し引かれることにより、ムラ補正が行われる。
【００５８】
次いで、各フレーム画像に対し、体動補正処理が行われる（ステップＳ１０５）。体動補正処理では、全フレーム画像の２箇所のＸ線非透過マーカＭ１、Ｍ２の線分が一致するように、各フレーム画像が回転、平行移動され、位置合わせが行われる。
【００５９】
前処理が終了すると、図４に戻り、画像解析・表示処理が実行される（ステップＳ１２）。
図６は、図４のステップＳ１２において実行される画像解析・表示処理を示すフローチャートである。画像解析・表示処理は、制御部３１と記憶部３２に記憶されているプログラムとの協働により実行される。なお、第２の実施の形態との区別のため、第１の実施の形態における画像解析・表示処理を画像解析・表示処理Ａとする。
【００６０】
まず、各フレーム画像から肺野領域の抽出が行われる（ステップＳ２０１）。肺野領域の抽出方法は何れの方法であってもよい。例えば、一連のフレーム画像中の任意のフレーム画像（ここでは、撮影順が一番（最初）のフレーム画像とする。）の各画素の信号値のヒストグラムから判別分析によって閾値を求め、各.フレーム画像においてこの閾値より高信号の領域を肺野領域候補として１次抽出する。次いで、１次抽出された肺野領域候補の境界付近でエッジ検出を行い、境界付近の小領域でエッジが最大となる点を境界に沿って抽出すれば肺野領域の境界を抽出することができる。抽出された肺野領域の情報は、制御部３１のＲＡＭに記憶される。
【００６１】
次いで、各フレーム画像において、横隔膜の位置の検出が行われる（ステップＳ２０２）。
横隔膜は、その上下運動によって肺の呼吸運動を促すものである。図７に、１つの呼吸サイクル（深呼吸時）において撮影された複数の時間位相Ｔ（Ｔ＝ｔ０〜ｔ６）のフレーム画像を示す。図７に示すように、呼吸サイクルは、呼気期と吸気期により構成される。呼気期は、横隔膜が上がることによって肺から空気が排出され、図７に示すように肺野の領域(肺野領域の面積)が小さくなる。最大呼気位では、横隔膜の位置が最も高い状態となる。吸気期は、横隔膜が下がることにより肺に空気が取り込まれ、図７に示すように胸郭中の肺野の領域(肺野領域の面積)が大きくなる。最大吸気位では、横隔膜の位置が最も下がった状態となる。安静呼吸時においては肺野の面積の変化は図７に示す深呼吸時より小さいが、同様の動きとなる。安静呼気位(安静呼吸時における呼気から吸気への変換点)にて、横隔膜の位置が最も高い状態となり、安静吸気位(吸気から呼気への変換点)にて、横隔膜の位置が最も低い状態となる。
【００６２】
ここで、図７からもわかるように、肺尖の上下位置は呼吸運動の影響をほとんど受けず、その位置がほとんど変わらないため、肺尖と横隔膜の垂直方向の距離は、横隔膜の上下方向の位置を表しているといえる。よって、動態画像の各フレーム画像の肺尖と横隔膜の垂直方向の距離を横隔膜の位置を示す情報として取得することができる。
【００６３】
例えば、ステップＳ２０２においては、各フレーム画像について以下の処理が行われる。
左右それぞれの肺野領域（図８のＲ）から肺尖及び横隔膜の基準位置Ｐ１、Ｐ２が特定される。例えば、肺尖の基準位置Ｐ１を肺野領域Ｒの最上端の位置と予め定義しておき、肺野領域Ｒにおける垂直方向の最も上方にある位置を抽出することによって肺尖の基準位置Ｐ１を特定する。また、横隔膜の基準位置Ｐ２を横隔膜のカーブＣ（図８に点線で示す）の垂直方向の平均位置と予め定義しておき、肺野領域Ｒから横隔膜のカーブＣを抽出し、その垂直方向の平均位置を求め、求めた位置を横隔膜の基準位置Ｐ２として特定する。そして、特定された肺尖の基準位置Ｐ１と横隔膜の基準位置Ｐ２の垂直方向の位置（Ｙ座標）の間の距離Ｄを算出する。算出された距離Ｄが横隔膜の位置を示す情報として取得され、制御部３１のＲＡＭに記憶される。
【００６４】
次いで、胸郭の幅の算出が行われる（ステップＳ２０３）。
胸郭は、上部胸郭（第２〜第６肋骨（上肋骨））と下部胸郭（第７〜第１０肋骨（下肋骨））により構成される。図９に示すように、上部胸郭は、吸気期に肋骨と胸骨を上方に引き上げ、胸郭の前後径を増す動きをする。下部胸郭は、吸気期に肋骨を外上方に引き上げ、胸郭の左右系を増す動きをする。呼気期は、これと逆の動きをする。ステップＳ２０３においては、各フレーム画像の上部胸郭と下部胸郭のそれぞれの幅を算出する。
【００６５】
ステップＳ２０３における胸郭の幅の算出方法としては、例えば、肺尖から所定の高さ(垂直方向の距離)における左右肺野領域の外側端(胸郭端)の距離を、各フレーム画像における胸郭の幅（上部胸郭の幅、下部胸郭の幅）として算出する。肺尖からの所定の高さとしては、例えば、まず、安静吸気位もしくは安静呼気位のフレーム画像に対して、胸郭の幅(左右肺野領域の外側端間の距離)の最大値を算出し、この最大値の関数として予め記憶部３２に記憶された、上肋骨、下肋骨のそれぞれの肺尖からの距離を、算出された最大値をもとに記憶部から読み出し、読み出した肺尖からの距離に対して、胸郭の幅を算出しても良い。
また、上肋骨、下肋骨の肺尖からの距離は、性別、年齢、身長、体重等の患者情報にも依存するので、患者情報(性別、年齢、身長、体重の何れか又は組み合わせ)と、上肋骨、下肋骨の肺尖からの距離とを対応付けて予め記憶部にテーブルとして記憶しておき、患者情報をもとに記憶部３２から上肋骨、下肋骨の肺尖からの距離を読み出し、読み出した肺尖からの距離に対して、胸郭の幅を算出しても良い。
また、肋骨を認識し、所定の肋骨 (例えば、上肋骨として第４肋骨、下肋骨として第８肋骨)の位置における胸郭の幅を算出することとしても良い。所定の肋骨を認識する方法としては、例えば、以下の公知文献に記載の手法により、安静吸気位もしくは安静呼気位のフレーム画像に対して、まずRobinsonオペレータ等のエッジ抽出フィルタを用いてエッジ画像を生成し、次にエッジ画像から円弧形状を検出するハフ変換などを用いて肋骨らしい円弧形状線成分を見つけることで肋骨形状を抽出し、この抽出した肋骨を上から数えていくことにより、所定の肋骨を認識することができる。（「ハフ変換および線形状性を利用した胸部X線像の主要陰影の辺縁抽出」電子情報通信学会論文誌 D-II Vol.J77-D-II No.7 pp1375-1381）。
算出された胸郭の幅の情報は、制御部３１のＲＡＭに記憶される。
【００６６】
次いで、記憶部３２に記憶されている表示態様の設定情報が読み出され、解析結果の表示対象となる領域をユーザが選択することが設定されているか又は予め定められた領域とすることが設定されているかが判断される（ステップＳ２０４）。表示対象領域をユーザが選択することが設定されていると判断されると（ステップＳ２０４；ＹＥＳ）、処理はステップＳ２０５に移行する。表示対象領域を予め定められた領域とすることが設定されていると判断されると（ステップＳ２０４；ＮＯ）、処理はステップＳ２１０に移行する。
【００６７】
ステップＳ２０５においては、各フレーム画像の肺野領域が複数の小領域に分割され、各フレーム画像の小領域が互いに対応付けられる（ステップＳ２０５）。各フレーム画像間における小領域の位置の対応関係は、制御部３１のＲＡＭに記憶される。
【００６８】
ここで、図７に示すように、呼吸運動によって肺野領域の各部分の位置は時間とともに変化する(下肺野ほど上下方向の位置の変化量が大きく、また、各部分の位置の変化量は呼吸が深いほど大きくなる)。すなわち、各フレーム画像間で肺野の同一部分を示す画素の位置はずれることになる。しかし、安静呼吸時に撮影された画像においては、図７に示す深呼吸時と比べて肺野の動きは小さく、後述する解析結果に大きな誤差を生じるほど肺野の同一部分に対する位置ズレはおきない。
【００６９】
図１１の画像Ｄ１は、安静呼気位（安静呼吸時において横隔膜の位置が最も高くなったタイミング）のフレーム画像である。図１１の画像Ｄ２は、安静吸気位（安静呼吸時において横隔膜の位置が最も低くなったタイミング）のフレーム画像である。即ち、図１１の画像Ｄ１とＤ２とは、呼吸１サイクルにおいて最も形状の差の大きいタイミングで撮影された画像である。しかし、図１１の画像Ｄ１、Ｄ２間においては、最も位置ずれの大きい肺野領域の下部領域においても位置ずれはわずかであることがわかる（画像Ｄ２のＡ１１は画像Ｄ１のＡ１と同じ画素位置を示し、画像Ｄ２のＡ２は画像Ｄ１のＡ１と肺野における同一部分を描画した領域を示している）。
【００７０】
そこで、ステップＳ２０５における具体的な処理としては、まず、一連のフレーム画像の中から一のフレーム画像を基準画像として設定する。次いで、基準画像の抽出された肺野領域を複数の小領域（例えば、０．２〜４ｃｍの矩形領域）に分割する（図１０参照）。次いで、他のフレーム画像の肺野領域を基準画像の各小領域と同じ画素位置の小領域（放射線検出部１３１の同じ撮像素子から出力される信号値の領域）に分割する。次いで、各フレーム画像間の同じ画素位置の各小領域を互いに対応付ける。この処理では、フレーム画像の小領域への分割及び対応付けを高速に行うことが可能となる。
基準画像としては、安静呼気位のフレーム画像とすることが好ましい。安静呼気位では、安静呼吸時において横隔膜の位置が最も高くなる、即ち、肺野領域の面積が最も小さくなるので、基準画像の小領域を他のフレーム画像に対応付けたときに、小領域が他のフレーム画像の肺野外の領域に対応付けられることがないためである。
安静呼気位の画像は、一連のフレーム画像の中から横隔膜の位置が最も高い位置にある画像を抽出することで取得することができる。
【００７１】
上記実施例は、一連のフレーム画像のうちの一つのフレーム画像を基準画像と設定し、当該基準画像を用いて抽出された肺野領域に対応する放射線検出器１３１の画素群が、他のフレーム画像に於いても、肺野領域に係る放射線検出器１３１の画素群と見做し、放射線検出器１３１の画素単位に以後の演算処理を行う、いわゆる見做し方式である。
これに対し、処理時間を要するものの解析精度をより高める方式として、いわゆるローカルマッチング処理及びワーピング処理を用い、一連のフレーム画像にわたる実際の肺野領域の精確な対応付けをはかることも可能である。
【００７２】
撮影された動態画像が深呼吸時の画像である場合は、図１２に示すように、肺野の同一の部分を示す画素位置が大幅にずれる。ここで、図１２の画像Ｄ３は、深呼吸時の最大呼気位のフレーム画像であり、図１２の画像Ｄ４は、深呼吸時の最大吸気位のフレーム画像である（画像Ｄ４のＡ３１は画像Ｄ３のＡ３と同じ画素位置を示し、画像Ｄ４のＡ４は画像Ｄ３のＡ３と肺野における同一部分を描画した領域を示している）。そのため、安静呼吸時と同様に各フレーム画像における、基準画像の各小領域と同じ画素位置の領域をその小領域に対応する領域としてしまうと、後述する解析結果に位置ずれの信号値変化に起因する大きな誤差を生じ、診断に利用し得ないものとなる。そこで、このような場合においても、各フレーム画像間の対応点を抽出する対応点抽出処理（ローカルマッチング処理）及び非線形歪変換処理（ワーピング処理）を行って、肺野領域の同一部分が描画された領域を各フレーム画像間で対応付けることが好ましい。
【００７３】
ローカルマッチング処理では、まず、撮影順が１番（最初）のフレーム画像から抽出した肺野領域を、例えば、例えば、０．２〜４ｃｍ角の矩形からなる小領域に分割する。
次いで、撮影順が１番のフレーム画像をＦ１、これと隣接するフレーム画像（撮影順が隣接するフレーム画像（即ち、時間的に隣接するフレーム画像。以下同様。））をＦ２とし、Ｆ２に、Ｆ１の各小領域の探索領域を設定する。ここで、Ｆ２の探索領域は、Ｆ１における各小領域における中心点の座標を（ｘ，ｙ）とすると、同一の中心点（ｘ，ｙ）をもち、Ｆ１の小領域よりも縦横の幅が大きくなるように設定する（例えば、１．５倍）。そして、Ｆ１の各領域毎に、Ｆ２の探索範囲で最もマッチング度合いが高くなる領域を求めることで、Ｆ１の各小領域に対するＦ２上での対応位置を算出する。マッチング度合いとしては、最小二乗法や相互相関係数を指標に用いる。そして、Ｆ２の肺野領域をＦ１の各小領域の対応位置で分割する。
次いで、Ｆ２を、新たにＦ１とみなし、撮影順がＦ２の次のフレーム画像を新たなＦ２とみなして、Ｆ１の各小領域におけるＦ２の対応位置を算出する。以上の処理を繰り返すことで、各フレーム画像の各小領域が隣接するフレーム画像のどの位置に対応するかが求まる。求めた処理結果は、制御部３１のＲＡＭに記憶される。
【００７４】
次いで、ワーピング処理が行われる。具体的には、撮影順が１番のフレーム画像をＦ１、これと撮影順が隣接する（時間的に隣接する）フレーム画像をＦ２とし、上記ローカルマッチング処理で算出された隣接するフレーム画像間の各小領域の対応位置に基づいて、各小領域毎にＦ１からＦ２へのシフトベクトルを算出する。次いで、算出されたシフトベクトルを多項式でフィッティングして、この多項式を用いて各小領域における各画素のシフトベクトルを算出する。そして、算出された各画素のシフトベクトルに基づいて、ワーピング処理を行い、Ｆ２の各小領域内の各画素の位置をＦ１のフレーム画像の対応する画素の位置にシフトする。次いで、ワーピング処理されたＦ２を、新たにＦ１とみなし、撮影順がＦ２の次のフレーム画像を新たなＦ２とみなして、上記処理を行う。以上の処理を撮影順の早い隣接フレーム画像間から順次繰り返すことで、全てのフレーム画像の各小領域の位置を撮影順が１のフレーム画像（ローカルマッチング及びワーピング処理における基準画像）に略一致させることが可能となる。各フレーム画像間における小領域の位置の対応関係は、制御部３１のＲＡＭに記憶される。
【００７５】
次いで、各フレーム画像において、各小領域内の画素の平均信号値が算出される（ステップＳ２０６）。平均化された画素信号値は、制御部３１のＲＡＭに記憶される。
【００７６】
次いで、横隔膜の位置、胸郭の幅、各小領域の平均信号値のそれぞれについて、時間変化を示す波形が作成される（ステップＳ２０７）。具体的には、横軸を撮影開始からの経過時間、縦軸を各値（横隔膜の位置を示す値、胸郭の幅の値、又は画素の平均信号値）とした座標平面を作成し、各フレーム画像の撮影開始からの経過時間及び算出された各値の点がプロットされることにより、各値の時間変化を示す波形が取得される。更に、各小領域の平均信号値の時間変化を時間方向のローパスフィルタ（例えば、カットオフ周波数０．５Ｈｚ）でフィルタリングすることで、換気量を示す信号値の時間変化（換気量の時間変化）を示す波形を取得することができる。また、各小領域の平均信号値の時間変化を時間方向のハイパスフィルタ（例えば、カットオフ周波数０．７Ｈｚ）でフィルタリングすることで、血流量を示す信号値の時間変化（血流量の時間変化）を示す波形を取得することができる。横隔膜の位置、胸郭の幅の時間変化は、肺野全体の動きを示す特徴量である。各小領域の平均信号値の時間変化をローパスフィルタ処理したものは、各小領域毎の動き（換気）を示す特徴量である。各小領域の平均信号値の時間変化をハイパスフィルタ処理したものは、各小領域毎の動き（血流）を示す特徴量である。
【００７７】
そして、取得された波形が解析され、肺野全体の動き、及び各小領域の動きを示す特徴量（時間変化以外の特徴量）が更に算出される（ステップＳ２０８）。各小領域について算出される特徴量は、操作部３３を介してユーザにより予め選択された特徴量が算出される。即ち、表示態様の設定情報に基づいて算出される。例えば、以下の（１）〜（５）のうち、ユーザにより選択された特徴量が算出される。
なお、呼吸１サイクルは、ここでは縦軸の値が極小値となるタイミングから次の極小値となるタイミングまでをさす。吸気期（吸気時）は、縦軸の値が極小から次に極大になるまでの期間をさす。呼気期（呼気時）は、縦軸の値が極大から次に極小になるまでの期間をさす。波形に複数の呼吸サイクルが含まれている場合は、何れか一つについての特徴量を算出してもよいし、複数サイクルのそれぞれにおいて特徴量を算出して平均化してもよい。
【００７８】
ここで、図１３は、一人の患者に正常系及び異常系が混在すると想定される場合に於ける各波形を模式的に示している。図１３のＡは、ある肺野領域における正常部分の小領域の換気量を示す信号値の時間変化を示す波形の一例である。図１３のＢは、拘束性部分の小領域の換気量を示す信号値の時間変化を示す波形の一例である。図１３のＣは、閉塞性部分の小領域の換気量を示す信号値の時間変化を示す波形の一例である。図１３のＤは、横隔膜位置の時間変化を示す波形の一例である。
【００７９】
（１）換気量の振幅
換気量の振幅（図１３のＡ〜ＣにＡｎで示す）は、各小領域の換気量を示す指標である。
換気量の振幅は、各小領域の上記換気量を示す信号値の時間変化を示す波形の、呼吸１サイクルにおける最大信号値−最小信号値を算出することにより求めることができる。
各小領域の換気量の振幅を算出することによって、局所的に換気機能が低下し、換気量が小さい（換気が十分に行われていない）部分を特定することが可能となる。例えば、図１３のＢ、Ｃに示す波形を有する小領域のように、換気量の振幅が予め定められた閾値より小さい小領域は、局所的に換気機能が低下していると判断することができる。
また、肺野全体において、各小領域の換気量の振幅のばらつきが大きい場合、局所的な換気機能低下部分が散在していると考えられる。この場合、例えば、ブラ、気腫化、気胸、無気肺、肺炎、肺水腫等の局所的換気障害の疑いがあると判断することができる。
一方、肺野全体において、各小領域の換気量の振幅のばらつきが小さい場合であって、全体的に換気量の振幅が小さい場合には、一様に換気機能が低下していると考えることができる。この場合、間質性肺炎等の拘束性肺疾患の疑いがあると判断することができる。
【００８０】
（２）吸気時間に対する呼気時間の延長度合い
吸気時間に対する呼気時間の延長度合い（図１３のＣのｔ１／ｔ２）は、各小領域の呼気時の動きを示す指標である。
吸気時間に対する呼気時間の延長度合いは、「呼気時間／吸気時間」によって算出することができる。呼気時間は、上記換気量を示す信号値の時間変化を示す波形の、呼吸１サイクルにおいて信号値が最大値となった時点から最小値となった時点までの時間ｔ２を算出することにより求めることができる。吸気時間は、上記換気量を示す信号値の時間変化を示す波形の、呼吸１サイクルにおいて信号値が最小値となった時点から最大値となった時点までの時間ｔ１を算出することにより求めることができる。または、吸気時間に対する呼気時間の延長度合いは、「吸気の最大微分値／呼気の最大微分値」によっても算出することができる。吸気、呼気の最大微分値は、呼吸１サイクルにおいて、隣接するフレーム画像間（撮影順が（時間的に）隣接するフレーム画像間。以下同様）の対応する小領域毎に換気量を示す信号値の差分値を求め、吸気期、呼気期のそれぞれにおいてフレーム間差分値の絶対値が最大となる値を求めることにより得られる。
閉塞性肺疾患では、特に呼気における末梢気道抵抗（粘性抵抗）が増加し、吸気時間に対して呼気時間が増大する。そのため、肺野全体的に、吸気時間に対する呼気時間の延長度合いの値が大きくなる。
一方、肺野全体において、吸気時間に対する呼気時間の延長度合いのばらつきが大きい場合は、混合性肺疾患の疑いがあると考えられる。
【００８１】
（３）吸気開始遅延時間
吸気開始遅延時間（図１３のＡ〜ＣにＴで示す）は、小領域毎の肺の硬さを示す指標である。
吸気開始遅延時間は、次のように求めることができる。まず、横隔膜の位置（肺尖から横隔膜までの距離）の時間変化を示す波形から、肺尖から横隔膜までの距離が極小値となる安静呼気位のタイミング（時刻）を取得する。次いで、上記各小領域の換気量を示す信号値の時間変化を示す波形において、上記取得したタイミングと略同タイミングに得られる極小値を基準として、信号値とその基準との差分が予め定められた閾値ｔｈ以上となったタイミング（時刻）を取得する。そして、上記取得した２つのタイミング間の時間を、吸気開始遅延時間Ｔとして算出する。閾値ｔｈは、小領域での平均信号値に対するノイズレベルをもとに一意に定められた信号値であり、ノイズによる信号変化が吸気開始時間と誤認識されない程度の大きさをもった信号値である。この閾値ｔｈは、信号値ではなく微分値（傾き）で定義してもよい。このとき、各小領域に対して、撮影順が隣接するフレーム画像間でのフレーム間差分値を微分値として算出し、肺尖から横隔膜までの距離が極小値となるタイミングから対象とする小領域の平均信号値の微分値がこの閾値を超えたタイミングまでの時間を吸気遅延時間として算出する。
肺が硬いと、吸気遅延時間Ｔが増大するので、吸気の開始の遅延時間Ｔを求めることによって弾性抵抗が増加し、局所的に肺が硬くなっている部分を特定することができる。吸気開始遅延時間が全体的に大きいと、間質性肺炎などの拘束性肺疾患の疑いがあると判断することができる。吸気開始遅延時間が全体的に小さいと、肺が柔らかくなる閉塞性肺疾患の疑いがあると判断することができる。
なお、呼気遅延算出時間を算出することとしてもよい。
【００８２】
（４）吸気気流速度
吸気気流速度は、各小領域の肺のやわらかさ（肺コンプライアンス）を示す指標である。
吸気気流速度は、各小領域の上記換気量を示す信号値の時間変化を示す波形から吸気時の微分値の代表値を算出することにより求めることができる。具体的には、各小領域において、隣接するフレーム画像間の差分値を求め、吸気時における差分値の代表値を吸気気流速度として算出する。代表値としては、吸気時の微分値の最大値、又は、吸気時の微分値の平均値とすることができる。または、吸気開始直後の気流速度に焦点をあてて、肺尖から横隔膜までの距離が極小値となる安静呼気位のタイミングから所定時間経過後の微分値を吸気気流速度としてもよい。
肺が硬いと吸気時の気流速度が低下するので、吸気時の気流速度を求めることによって、弾性抵抗が増加し局所的に肺が硬くなっている部分を特定することができる。吸気時の気流速度が小さい部分は肺が硬い（拘束性）部分である。吸気時の気流速度が大きい部分は、肺がやわらかい（閉塞性）部分である。全体的に吸気時気流速度が小さいと、拘束性肺疾患の疑いがあると判断することができる。全体的に吸気時気流速度が大きいと、閉塞性肺疾患の疑いがあると判断することができる。
【００８３】
（５）血流の振幅
血流の振幅は、各小領域の血流量を示す指標である。
血流量の振幅は、各小領域の上記血流量を示す信号値の時間変化を示す波形の、心拍１サイクルにおける最大信号値−最小信号値を算出することにより求めることができる。
各小領域の血流量の振幅を算出することによって、肺血流量が小さい部分を特定することが可能となる。
【００８４】
肺野全体の特徴量としては、以下の（６）〜（８）が算出される。
（６）肋骨位置（上肋骨、下肋骨）・横隔膜位置の移動量
肋骨位置・横隔膜位置の移動量は、胸式呼吸の度合い、胸郭の動きの制限度合いを示す指標である。
上肋骨の移動量は、胸郭の幅の時間変化を示す波形から、呼吸１サイクルにおける上肋骨の幅の最大値及び最小値を取得し、「最大値−最小値」（差分）を算出することにより求めることができる。同様に、下肋骨の移動量は、胸郭の幅の時間変化を示す波形から、呼吸１サイクルにおける下肋骨の幅の最大値及び最小値を取得し、「最大値−最小値」（差分）を算出することにより求めることができる。横隔膜の移動量は、上記肺尖から横隔膜までの距離の時間変化を示す波形から、呼吸１サイクルにおける最大距離と最小距離を取得し、「最大距離−最小距離」を算出することにより求めることができる。
「上肋骨の移動量／横隔膜の移動量」が大きい場合、即ち、横隔膜位置の移動量が少なく、上部胸式呼吸が優位である場合、横隔膜を使わずに呼吸をしている、即ち、重症な慢性換気障害の疑いがあると判断することができる。下肋骨の移動量が少ない場合には、慢性閉塞性肺疾患（ＣＯＰＤ）の疑いがあると判断することができる。
【００８５】
（７）呼吸数
呼吸周期は、換気障害があるか否かを示す指標である。
呼吸周期は、上記肺尖から横隔膜までの距離の時間変化を示す波形から呼吸周期を求め、呼吸周期の逆数から単位時間当たりの呼吸数を算出することができる。具体的には、上記肺尖から横隔膜までの距離の時間変化を示す波形における極小値から次の極小値までの時間間隔を呼吸周期として算出することができる。または、肺野全体の信号値を平均化し、肺野全体の信号値の時間変化を示す波形を求め、求めた波形における極小値から次の極小値までの時間間隔を呼吸周期として算出することができる。
呼吸運動が正常であれば、呼吸数は１５回前後（１２〜１８回）／分である。この範囲を超えている場合、換気障害があると判断することができる。
【００８６】
（８）鎖骨の移動量
鎖骨の移動量は、呼吸補助筋の動きを示す指標である。
鎖骨の移動量は、各フレーム画像から鎖骨領域を認識し、認識した鎖骨領域に基準点を設定し、その基準点の位置を追跡することで、鎖骨の移動量を算出することができる。鎖骨領域の認識は、例えば、まず、各フレーム画像に対して、Robinsonオペレータ等のエッジ抽出フィルタを用いてエッジ画像を生成し、エッジ画像から円弧形状を検出するハフ変換などを用いて肋骨らしい円弧形状線をみつけて肋骨形状を検出する。また、直線を検出するハフ変換などを用いて鎖骨らしい直線部分を見つけて鎖骨形状を検出する。この検出した鎖骨形状をもとに鎖骨領域を特定する（「ハフ変換および線形状性を利用した胸部Ｘ線像の主要陰影の辺縁抽出」電子情報通信学会論文誌D-II Vol.J77-D-II No.7 pp1375-1381参照）。
同様に、肩の移動量についても、呼吸補助筋の動きを示す指標として用いることができる。
横隔膜の運動が十分でない場合、呼吸補助筋を使ってその動きを補おうとする。呼吸補助筋を使った呼吸の場合、鎖骨位置、肩位置の移動量が大きくなる。したがって、安静呼吸時に鎖骨位置、肩位置の移動量が大きい場合、呼吸補助筋を使って安静呼吸をしていると判断でき、重症な慢性換気障害の疑いがある判断することができる。
【００８７】
特徴量の算出が終了すると、設定されている表示態様に基づいて、波形及び波形の解析結果が表示部３４の画面上に表示される（ステップＳ２０９）。
【００８８】
図１４に、換気の解析結果の表示（標準）が設定されている場合にステップＳ２０９において表示部３４に表示される表示画面３４１の一例を示す。ステップＳ２０９においては、まず、撮影された複数のフレーム画像のうちの一つのフレーム画像（基準画像）３４１ａと、肺野全体の動きを示す指標となる、横隔膜位置の時間変化を示す波形３４１ｂと、胸郭の幅の時間変化を示す波形３４１ｃと、が表示される。横隔膜位置の時間変化を示す波形３４１ｂ及び胸郭の幅の時間変化を示す波形３４１ｃを観察することにより、ユーザである医師は、例えば、胸式呼吸の度合い、胸郭の動きの制限度合いを把握し、肺野全体にわたる換気障害があるか否かを判定することができる。
表示画面３４１において、基準画像３４１ａの肺野領域内の一又は複数の小領域が操作部３３のマウス等により選択されると、選択された領域における換気量を示す信号値の時間変化を示す波形３４１ｄ、及び予め設定された特徴量の数値３４１ｅが表示される。なお、図１４におけるＸＸは、何らかの数値が入っていることを示す（図１６、図１７、図１８についても同様）。波形３４１ｄ、特徴量の数値３４１ｅは、選択された領域についてのみ表示してもよいし、図１４に示すように、選択された領域及びその領域と体軸対称の領域（右が選択されたら左、左が選択されたら右）の双方について表示してもよい。選択された領域についての波形３４１ｄを観察することにより、医師は局所的な換気機能の低下等があるか否かを把握することができる。また、この波形３４１ｄから算出された特徴量の数値３４１ｅを確認することで、波形から判定した換気障害についての確認を行うことが可能となる。また、選択された領域の波形を左右対称の領域の波形と比較して診断を行うことが可能となる。
なお、医師による波形観察での換気障害の判定を容易にするために、記憶部３２に正常な換気量の変化を示す波形（正常例）を予め記憶しておき、図１５Ａに示すように、正常例の波形を波形３４１ｄに並列表示させてもよい。また、正常例の波形との比較を行いやすくするために、図１５Ｂに示すように正常例の波形の周期を患者の波形の周期に一致するよう時間軸方向に伸縮させて表示してもよい（図１６、図１７、図１８、図１９、図２１についても同様）。また、上記の正常例の波形の代わりに、同一患者が過去に撮影した画像における、同じ領域に対する信号値の時間変化の波形を表示してもよい。
【００８９】
図１６に、換気と血流の解析結果の表示（標準）が設定されている場合にステップＳ２０９において表示部３４に表示される表示画面３４２の一例を示す。ステップＳ２０９においては、まず、撮影されたフレーム画像（基準画像）３４２ａと、肺野全体の動きを示す指標となる、横隔膜位置の時間変化を示す波形３４２ｂと、胸郭の幅の時間変化を示す波形３４２ｃと、が表示される。基準画像３４２ａの肺野領域内の一又は複数の小領域が操作部３３のマウス等により選択されると、選択された領域についての、換気量を示す信号値の時間変化を示す波形３４２ｄ、血流量を示す信号値の時間変化を示す波形３４２ｅ、及び予め設定された特徴量の数値３４２ｆが表示される。表示画面３４２は、上記表示画面３４１に血流量の信号値の時間変化を示す波形３４２ｅ及びその特徴量が追加されたものであるので、上記表示画面３４１と同様の効果を奏することができる。更に、換気と血流の情報を同時に表示するので、医師が肺の病態をより精度良く把握することが可能となる。例えば、ＣＯＰＤでは気腫化が進行した部分やブラなどの局所的な換気異常部分にて、局所的な血流の変化(血流量の減少)のほうがより顕著に現れる(より良く対応する)場合もある。このような場合においても異常を見落とすことなく診断を行うことができる。また、血流のみに異常が現れる病変についても把握することが可能となる。例えば、血流が小さく換気が正常な場合は肺塞栓症の可能性あり、といった診断も可能となる。
【００９０】
図１７に、追加情報として肺野内の各小領域の解析結果の静止画上への一覧表示が設定されている場合に、ステップＳ２０９において表示部３４に表示される表示画面３４３の一例を示す。ステップＳ２０９においては、まず、撮影されたフレーム画像（基準画像）３４３ａと、肺野全体の動きを示す指標となる、横隔膜位置の時間変化を示す波形３４３ｂと、胸郭の幅の時間変化を示す波形３４３ｃと、各小領域についての換気量（換気量及び血流量）の波形の解析結果を示す静止画３４３ｄ〜３４３ｈが表示される。図１７においては、波形の解析結果として、換気量、呼気時間の延長度合、吸気開始遅延時間、吸気時気流速度、血流量の特徴量の静止画表示が設定されている場合の表示の一例を示している。静止画３４３ｄ〜３４３ｈには、それぞれ各小領域の呼吸１サイクル毎の特徴量（項目に示す特徴量）の値を複数周期に渡って算出したものを平均化し、平均化された数値を予め記憶部３２に記憶されている特徴量と表示時のパラメータ値との変換テーブル(例えば、特徴量の大きさと色相、彩度、明度、輝度、透明度のいずれかを一対一に対応付けたテーブル)に基づいた表示用パラメータ値）に変換し、変換したパラメータ値、例えば、色で各小領域を色づけした画像を上記基準画像にオーバーレイした静止画像が表示される。静止画３４３ｄ〜３４３ｈを表示することにより、医師は、各小領域毎の特徴量の値の大小の全体像を色により把握することができる。そして、例えば、換気の少ない箇所、血流の少ない箇所、呼気時間の遅延度合いの大きい箇所等、局所的な異常のある小領域を特定することができる。
基準画像３４３ａの肺野領域内の一又は複数の小領域が操作部３３のマウス等により選択されると、選択された小領域についての、換気量を示す信号値の時間変化を示す波形３４３ｉ、血流量を示す信号値の時間変化を示す波形３４３ｊ、及び特徴量の数値３４３ｋが表示される。医師は、上述の静止画３４３ｄ〜３４３ｈを観察し、局所的に異常のある小領域が発見された場合、その小領域を基準画像３４３ａにおいて選択することで、局所的に異常のある小領域の波形や特徴量の数値を確認することが可能となる。
【００９１】
図１８に、追加情報として波形の解析結果の動画表示が設定されている場合に、ステップＳ２０９において表示部３４に表示される表示画面３４４の一例を示す。ステップＳ２０９においては、まず、撮影されたフレーム画像（動画像）３４４ａと、肺野全体の動きを示す指標となる、横隔膜位置の時間変化を示す波形３４４ｂと、胸郭の幅の時間変化を示す波形３４４ｃと、換気量を示す信号値の波形の解析結果を示す動画３４４ｄと、血流量を示す信号値の波形の解析結果を示す動画３４４ｅと、動画再生経過時間を示すインジケータ３４４ｆと、動画の再生、コマ送り、一時停止、停止等を指示するための操作ボタン３４４ｇと、が表示される。動画３４４ｄ、３４４ｅとしては、隣接するフレーム画像間の対応する小領域毎に換気量（血流量）を示す信号値の差分値を求め、求めたフレーム間差分値を予め記憶部３２に記憶されている特徴量の大きさと色との変換テーブルに基づいて対応する色に変換し、変換した色で各小領域を色づけした画像を、原画像（動画）にオーバーレイして表示する。この表示により、呼吸サイクルと換気量及び血流量の関係（呼吸サイクルにおける換気及び血流の変化）を医師が観察することが可能となる。その結果、医師は生理学的挙動をもとに、病態（異常か否か、重症度）を把握することが可能となる。
動画像３４４ａの肺野領域内の一又は複数の小領域が操作部３３のマウス等により選択されると、選択された小領域についての、換気量を示す信号値の時間変化を示す波形３４４ｈ、血流量を示す信号値の時間変化を示す波形３４４ｉ、及び特徴量の数値３４４ｊが表示される。医師は、上述の動画像３４４ｂ、３４４ｃ等を観察し、異常と思われる挙動をする領域が発見された場合、その小領域を動画像３４４ａにおいて選択することで、局所的に異常な挙動の見られる小領域の波形や特徴量の数値を確認し、診断の参考にすることが可能となる。
【００９２】
一方、ステップＳ２０４において、表示対象領域を予め定められた領域とすることが設定されていると判断されると（ステップＳ２０４；ＮＯ）、各フレーム画像の左右の肺野領域のそれぞれが上肺野、中肺野、下肺野の３つの小領域に分割される（ステップＳ２１０）。例えば、まず、基準画像の肺尖から横隔膜までの垂直方向の距離に基づいて、基準画像の肺野が３つの小領域に分割される。次いで、各フレーム画像において、基準画像における各領域の境界と同じ位置に境界が設定され、領域分割が行われる。
【００９３】
次いで、各フレーム画像において、各小領域内の画素の平均信号値が算出される（ステップＳ２１１）。平均信号値は、制御部３１のＲＡＭに記憶される。
【００９４】
次いで、横隔膜の位置、胸郭の幅、左右の上肺野、中肺野、下肺野の平均信号値、のそれぞれについて、時間変化を示す波形が作成される（ステップＳ２１２）。具体的には、横軸を撮影開始からの経過時間、縦軸をそれぞれの値（横隔膜の位置、胸郭の幅、画素の平均信号値）とした座標平面を作成し、各フレーム画像の撮影開始からの経過時間及び算出された各値の点をプロットすることにより、各値の時間変化を示す波形を取得する。更に、上肺野、中肺野、下肺野のそれぞれの平均信号値の時間変化を時間方向のローパスフィルタ（例えば、カットオフ周波数０．５Ｈｚ）でフィルタリングすることで、換気量を示す信号値の時間変化を示す波形を取得する。
【００９５】
次いで、上記作成された波形を表示する表示画面が表示部３４に表示される（ステップＳ２１３）。
図１９に、ステップＳ２１３において表示部３４に表示される表示画面３４５の一例を示す。図１９に示すように、表示画面３４５には、上肺野、中肺野、下肺野の各領域の位置を示すフレーム画像（基準画像）３４５ａと、肺野全体の動きを示す指標となる、横隔膜位置の時間変化を示す波形３４５ｂと、胸郭の幅の時間変化を示す波形３４５ｃと、上肺野、中肺野、下肺野のそれぞれの換気量を示す信号値の時間変化を示す波形３４５ｄ〜３４５ｆと、が表示される。波形３４５ｄ〜３４５ｆについては、それぞれ左右の肺野の波形が同じ座標空間上に表示される。横隔膜位置の時間変化を示す波形３４５ｂ及び胸郭の幅の時間変化を示す波形３４５ｃを観察することにより、医師は胸式呼吸の度合い、胸郭の動きの制限度合いを把握し、例えば、胸式呼吸優位、ＣＯＰＤ等の肺野全体にわたる換気障害があるか否かを診断することができる。また、左右の上肺野、中肺野、下肺野のそれぞれについて、換気量を示す信号値の時間変化を示す波形３４５ｄ〜３４５ｆを観察することにより、局所的な換気機能の低下等があるか否かを把握することができる。
【００９６】
なお、上肺野、中肺野、下肺野の各領域についての換気量を示す信号値の時間変化を示す波形とともに、片肺全体の平均信号値の時間変化を示す波形を表示することとしてもよい。これにより、医師が片肺全体の換気機能と、各領域の換気機能とを比較観察することができる。また、血流量を示す信号値の時間変化を算出し、この波形も併せて表示することとしてもよい。
【００９７】
以上、画像解析・表示処理について説明したが、各小領域の特徴量の値の分布や平均値から、自動的に被写体Ｍの肺全体としての病態を判定し、その結果を上記表示画面３４１〜３４５に併せて表示することとしてもよい。
【００９８】
例えば、肺野内小領域から取得された振幅の分散を算出し、分散が所定閾値より大きい場合、局所的な換気機能の低下部分が散在していると判定する。また、肺野内小領域から取得された振幅の平均値を算出し、算出した平均値が所定閾値より小さい場合、間質性肺炎等の拘束性肺疾患であると判定する。
また、肺野内小領域の吸気遅延時間の平均値が所定閾値より大きい場合、間質性肺炎等の拘束性疾患であると判定する。
肺野内小領域の吸気気流速度の平均値が所定閾値より小さい場合、間質性肺炎等の拘束性疾患であると判定する。
上記の各所定の閾値は、上記の各所定閾値は、性別、年齢、身長、体重 に依存するため、性別、年齢、身長、体重 に応じた閾値を予めテーブルとして記憶部３２に記憶しておき、図２のステップＳ１で入力された患者情報(性別、年齢、身長、体重)をもとに、対応する閾値を読み出し、病態判定に用いる。
【００９９】
〈第２の実施の形態〉
次に、本発明の第２の実施の形態について説明する。第２の実施の形態においては、図４のステップＳ１２で実行される画像解析・表示処理の内容が第１の実施の形態と異なるので、以下、第２の実施の形態における画像解析・表示処理（画像解析・表示処理Ｂ）について説明する。
【０１００】
図２０は、図４のステップＳ１２において制御部３１により実行される画像解析・表示処理Ｂのフローチャートである。画像解析・表示処理Ｂは、制御部３１と記憶部３２に記憶されているプログラムとの協働により実行される。
【０１０１】
まず、各フレーム画像から肺野領域の抽出（ステップＳ３０１）、横隔膜の位置の検出（ステップＳ３０２）、胸郭の幅の算出が行われる（ステップＳ３０３）。ステップＳ３０１〜ステップＳ３０３の処理は、図６のステップＳ２０１〜ステップＳ２０３で説明したものと同様であるので説明を援用する。
【０１０２】
次いで、各フレーム画像の肺野領域が複数の小領域に分割される（ステップＳ３０４）。ステップＳ３０４の処理は、ステップＳ２０６で説明したものと同様であるので説明を援用する。
【０１０３】
次いで、各フレーム画像において、各小領域内の画素の平均信号値が算出される（ステップＳ３０５）。
【０１０４】
次いで、横隔膜の位置、胸郭の幅、各小領域の平均信号値のそれぞれについて、時間変化を示す波形が作成される（ステップＳ３０６）。具体的には、横軸を撮影開始からの経過時間、縦軸をそれぞれの値（横隔膜の位置、胸郭の幅、画素の平均信号値）とした座標平面を作成し、各フレーム画像の撮影開始からの経過時間及び算出された値の点がプロットされることにより、各値の時間変化を示す波形が取得される。更に、各小領域の平均信号値の時間変化を時間方向のローパスフィルタ（例えば、カットオフ周波数０．５Ｈｚ）でフィルタリングすることで、換気量を示す信号値の時間変化を示す波形を取得することができる。
【０１０５】
次いで、取得された波形が解析され、各小領域の動きを示す特徴量が算出される（ステップＳ３０７）。各小領域について算出される特徴量としては、ここでは、第１の実施の形態において説明した換気量の振幅、呼気時間の延長度合い、吸気開始遅延時間が算出される。
【０１０６】
次いで、ステップＳ３０７で算出された特徴量に基づいて、各小領域が異常の可能性があるか否かが判断される（ステップＳ３０８）。例えば、各小領域の換気量の振幅、呼気時間の延長度合い、吸気開始遅延時間の何れかの値が、（１）体軸対称の領域（左右対称領域）で乖離が大きい（差分が予め定められた閾値以上である）場合、若しくは（２）予め定められた閾値以上を越えた（又は下回った）場合に、その小領域は異常の可能性があると判断される。
【０１０７】
そして、表示部３４に、異常の可能性の判定結果が表示される（ステップＳ３０９）。図２１に、ステップＳ３０９において表示される表示画面３４６の一例を示す。図２１に示すように、表示画面３４６には、異常の可能性があると判定された小領域を色付きの枠等で囲ったアノテーションが表示された基準画像３４６ａと、肺野全体の動きを示す指標となる、横隔膜位置の時間変化を示す波形３４６ｂと、胸郭の幅の時間変化を示す波形３４６ｃと、異常と判断された小領域における換気量を示す信号値の時間変化を示す波形３４６ｄ、３４６ｆと、が表示される。なお、図２１においては、異常の可能性があると判定された小領域が２つ検出された場合を例としているが、異常の可能性があると判定された各小領域の数に応じた波形が表示される。このように、異常の可能性がある領域を自動的に判定し、その領域の位置及び換気量を示す信号値の時間変化を示す波形を表示するので、医師は、効率よく、異常の可能性がある領域の波形を確認することができる。同時に、医師の疾患の見落としを低減することが可能となる。
【０１０８】
なお、波形は換気量を示す信号値の時間変化の波形だけでなく、血流量を示す信号値の時間変化の波形（図２１の波形３４６ｅ、波形３４６ｇ）を併せて表示することとしてもよい。この場合、ステップＳ３０６においては、各小領域の信号値の時間変化を時間方向のローパスフィルタ（例えば、カットオフ周波数０．５Ｈｚ）でフィルタリングするとともに、ハイパスフィルタ（例えば、カットオフ周波数０．７Ｈｚ）でフィルタリングすることで、血流量を示す信号値の時間変化を示す波形を取得する。ステップＳ３０７においては、換気量の振幅、呼気時間の延長度合い、吸気開始遅延時間及び血流の振幅が算出される。ステップＳ３０８においては、各小領域の換気量の振幅、呼気時間の延長度合い、吸気開始遅延時間、血流の振幅の何れかの値が、（１）体軸対称の領域（左右対称領域）で乖離が大きい（差分が予め定められた閾値以上である）場合、若しくは（２）予め定められた閾値以上を越えた（又は下回った）場合に、その小領域は異常の可能性があると判断される。そして、ステップＳ３０９においては、図２１に示すように、異常の可能性があると判定された領域の換気量を示す波形とともに、血流量を示す信号値の波形を併せて表示する。
【０１０９】
このように、換気と血流の情報を同時に表示することで、医師が肺の病態をより精度良く把握することが可能となる。例えば、ＣＯＰＤでは気腫化が進行した部分やブラなどの局所的な換気異常部分にて、局所的な血流の変化(血流量の減少)のほうがより顕著に現れる(より良く対応する)場合もある。このような場合においても異常を見落とすことなく診断を行うことができる。また、血流のみに異常が現れる病変についても把握することが可能となる。例えば、血流が小さく換気が正常な場合は肺塞栓症の可能性あり、といった診断も可能となる。
【０１１０】
また、算出した各特徴量を複合的に判断して異常の可能性があると判断した領域（異常領域）の病態（症例名）を推定し、推定される病態を表示画面３４６等に併せて表示することとしてもよい。また、病態と色との対応関係を予め記憶部３２に記憶しておき、推定された病態に対応する色で基準画像３４６ａ上の異常領域を表示することとしてもよい。このようにすれば、医師に病態についてのより詳しい診断支援情報を提供することが可能となる。病態の推定は、例えば、上記ステップＳ３０７で算出された特徴量のそれぞれが予め定められた閾値を越えたか否かに基づいて、各特徴量を大、中、小、に分類し、各特徴量の分類の組み合わせに基づいて推定することができる。以下の（例１）〜（例２）に、病態の推定基準の例を示す。
【０１１１】
（例１）
換気量（振幅）：小、呼気時間の延長度合い：小、吸気開始遅延時間：大（吸気時気流速度：小）、血流量：小  →拘束性部分
（例２）
換気量（振幅）：小、呼気時間の延長度合い：大、吸気開始遅延時間：小（吸気時気流速度：大）、血流量：小  →閉塞性部分
【０１１２】
各小領域に対して異常判定を行う際の各特徴量の閾値としては、例えば、正常値(標準値)をもとに算出した値を用いる。
各小領域における振幅(換気量)、遅延時間、(吸気)気流速度、血流量は、正常な場合、肺底から上肺野に向かってトレンドがあり、上肺野ほど、振幅は小さく、遅延時間は大きく、(吸気)気流速度は小さく、血流量は小さくなる。
そこで、上記各特徴量に対して、肺底からの垂直位置ｔｒ(安静吸気位の肺野の長さに対する、対象とする局所領域の肺底を基準とした垂直位置 の割合[%]として定義)に応じた各特徴量の正常値(標準値)を予め記憶しておく。また、ここでの正常値は、肺底からの垂直位置が略同じ複数の小領域の特徴量を水平方向に平均化した、平均値に対する正常値とする。さらに、各特徴量は、性別、年齢、身長、体重 に依存するため、性別、年齢、身長、体重に応じて、肺底からの垂直位置に対する各特徴量の正常値を、予めテーブルとして記憶部３２に記憶しておく。そして、図２のステップＳ１で入力された患者情報(性別、年齢、身長、体重)と、対象とする小領域の肺底からの垂直位置をもとに、対応する正常値を読み出し、例えば、対応する正常値±２０％(下限と比較する場合は−２０%、上限と比較する場合は＋２０％)の値を、異常判定の閾値として用いる。
この性別、年齢、身長、体重に応じた正常値は、例えば、性別、年齢、身長、体重が異なる各グループに対して、複数の健常者のデータから、肺底からの垂直位置に対する各特徴量の平均値を算出し、この平均値を正常値として用いることができる。
【０１１３】
以上説明したように、本発明に係る診断支援システム１００によれば、診断用コンソール３の制御部３１は、撮影用コンソール２から送信された胸部の動態を示す複数のフレーム画像のそれぞれのフレーム画像から肺野領域を抽出し、当該抽出された肺野領域を複数の小領域に分割し、複数のフレーム画像間において分割された小領域を互いに対応付けて解析することにより、分割された小領域の動きを示す予め設定された特徴量を算出する。そして、操作部３３により解析結果の表示対象となる領域が選択されると、選択された小領域についての特徴量を表示する。
【０１１４】
従って、医師であるユーザは、聴診器をあてるように肺野内の小領域を操作部３３により選択していくことで、選択した小領域の動きを示す特徴量を取得することができる。このように、診断支援システム１００によれば、診断に使用する情報を目視系に統一し、聴診器の経験の浅い医師であっても精度良い診断を行うことが可能なＧＵＩを含むシステムを提供することができる。
【０１１５】
また、制御部３１は、更に、肺野領域全体についての予め定められた特徴量を算出し、操作部３３により選択された領域の動きを示す特徴量と、肺野領域全体の動きを示す特徴量を同時に表示部３４に表示するので、肺野全体と選択した局所的な領域の双方についての診断に係る情報を一度に把握することが可能となる。
【０１１６】
また、制御部３１は、各小領域の動きを示す一又は複数の特徴量を算出するので、医師は、選択した小領域についての一又は複数の特徴量を取得することが可能となる。
【０１１７】
また、表示部３４に表示された静止画のフレーム画像上から操作部３３により解析結果の表示対象となる一又は複数の小領域を選択可能な構成とすることで、ユーザは、聴診器をあてるのと同じように容易な操作で一又は複数の領域を選択してその領域の解析結果を表示させることが可能となる。
【０１１８】
また、表示部３４に表示された動画から操作部３３により解析結果の表示対象となる一又は複数の領域を選択可能な構成とすることで、ユーザは、動画を見て異常と思われる挙動をする一又は複数の小領域を容易に選択してその領域の解析結果を表示させることが可能となる。
【０１１９】
また、表示部３４に、更に複数のフレーム画像のうちの一つを静止画表示し、そのフレーム画像の各小領域を算出された特徴量の値に応じた色で表示することで、ユーザは、各小領域における特徴量を一瞥して把握することが可能となり、肺野全体の中での局所的な異常個所を容易に把握することが可能となる。
【０１２０】
また、操作部３３により選択された小領域の動きを示す特徴量として、換気量（及び血流量）の時間変化を波形として表示することで、選択された領域についての換気及び血流の時間変化をユーザが容易に把握することが可能となる。
【０１２１】
また、操作部３３により肺野領域の一方の領域が選択された場合に、選択された領域の特徴量とともに、選択された領域と体軸対称の他方の肺野の領域の特徴量を同時に表示することで、左右の対称領域での特徴量の比較を容易に行うことができ、診断に有用な情報をユーザに提供することができる。
【０１２２】
また、算出される特徴量は、操作部３３により複数のなかから予め選択可能に構成することで、ユーザが所望する特徴量を表示させることが可能となる。
【０１２３】
なお、上述した本実施の形態における記述は、本発明に係る好適な診断支援システムの一例であり、これに限定されるものではない。
例えば、上記実施の形態において、各小領域毎の換気量を示す波形として、縦軸を換気量を示す信号値、横軸を時間としてプロットしたグラフを表示したが、波形の表示方法はこれに限定されない。例えば、縦軸を気流速度に相当する信号値の微分値（信号値／Ｓ）、横軸を換気量を示す信号値としてプロットしたグラフを表示してもよい。これは、スパイロのフローボリューム曲線に相当するグラフであるため、医師はスパイロと同様な感覚で診断情報を読み取ることが可能となる。
【０１２４】
また、上記実施の形態においては、予め各小領域について解析を行って特徴量を算出しておき、操作部３３により選択された小領域についての解析結果（波形、特徴量の値）を表示することとしたが、操作部３３により選択された領域について解析を行って特徴量を算出することとしてもよい。
【０１２５】
また、上記実施の形態においては、横隔膜の位置の時間変化を示す波形及び胸郭の幅の時間変化を示す波形を肺野全体の動きを示す特徴量として表示することとしたが、これに限定されない。例えば、肺野全体の平均信号値の時間変化を示す波形、若しくは、左右それぞれの肺に対する平均信号値の時間変化を示す波形を肺野全体の動きを示す特徴量として表示してもよい。また、波形で表示する特徴量を操作部３３を介してユーザが設定できることとしてもよい。
【０１２６】
また、上記実施の形態においては、各小領域について同一の種類の特徴量を算出することとして説明したが、小領域毎に異なる特徴量を算出することとしてもよい。
【０１２７】
また、例えば、上記の説明では、本発明に係るプログラムのコンピュータ読み取り可能な媒体としてハードディスクや半導体の不揮発性メモリ等を使用した例を開示したが、この例に限定されない。その他のコンピュータ読み取り可能な媒体として、ＣＤ-ＲＯＭ等の可搬型記録媒体を適用することが可能である。また、本発明に係るプログラムのデータを通信回線を介して提供する媒体として、キャリアウエーブ(搬送波)も適用される。
【０１２８】
その他、診断支援システム１００を構成する各装置の細部構成及び細部動作に関しても、本発明の趣旨を逸脱することのない範囲で適宜変更可能である。
【０１２９】
なお、明細書、請求の範囲、図面及び要約を含む２０１０年８月２７日に出願された日本特許出願Ｎｏ．２０１０−１９０２４０号の全ての開示は、そのまま本出願の一部に組み込まれる。
【産業上の利用可能性】
【０１３０】
医療の分野における胸部動態等の診断に利用可能性がある。
【符号の説明】
【０１３１】
１００  診断支援システム
１  撮影装置
１１  放射線発生装置
１１１  放射線源
１１２  放射線源保持部
１１３  支持基軸
１２  放射線照射制御装置
１３  放射線検出部
１３１  放射線検出器
１３２  検出器保持部
１３３  支持基軸
１４  読取制御装置
２  撮影用コンソール
２１  制御部
２２  記憶部
２３  操作部
２４  表示部
２５  通信部
２６  バス
３  診断用コンソール
３１  制御部
３２  記憶部
３３  操作部
３４  表示部
３５  通信部
３６  バス
【図１】
0005874636.tif 000002
【図２】
0005874636.tif 000003
【図３】
0005874636.tif 000004
【図４】
0005874636.tif 000005
【図５】
0005874636.tif 000006
【図６】
0005874636.tif 000007
【図７】
0005874636.tif 000008
【図８】
0005874636.tif 000009
【図９】
0005874636.tif 000010
【図１０】
0005874636.tif 000011
【図１１】
0005874636.tif 000012
【図１２】
0005874636.tif 000013
【図１３】
0005874636.tif 000014
【図１４】
0005874636.tif 000015
【図１５Ａ】
0005874636.tif 000016
【図１５Ｂ】
0005874636.tif 000017
【図１６】
0005874636.tif 000018
【図１７】
0005874636.tif 000019
【図１８】
0005874636.tif 000020
【図１９】
0005874636.tif 000021
【図２０】
0005874636.tif 000022
【図２１】
0005874636.tif 000023