(57)【要約】  肺野の換気に係る病態の把握を容易にし、診断に有効となる情報を医師に提供できるようにする。特に、聴診器による実務経験の少ない医師であっても、的確な診断を行える診断情報を提供する。  本発明に係る胸部診断支援システム１００の診断用コンソール３によれば、制御部３１は、胸部の動態を撮影して得られた複数のフレーム画像のそれぞれから肺野領域を抽出し、抽出した肺野領域を複数の小領域に分割し、複数のフレーム画像間において各小領域を対応付ける。次いで、複数のフレーム画像の対応する小領域毎に解析を行って、小領域毎の吸気の特徴量及び呼気の特徴量を算出し、算出された吸気の特徴量と呼気の特徴量との比の値を算出し、当該算出された比の値のヒストグラムを作成する。そして、作成したヒストグラムを表示部３４に表示する。
2012026146.tif 000002
【特許請求の範囲】
【請求項１】
少なくとも１サイクルの呼吸周期を含む胸部の動態を撮影し、連続する複数のフレーム画像を生成する撮影手段と、
前記撮影手段により生成された複数のフレーム画像のそれぞれから肺野領域を抽出する抽出手段と、
前記抽出手段により抽出された肺野領域を複数の小領域に分割し、前記複数のフレーム画像間にわたり小領域を対応付ける領域分割手段と、
前記複数のフレーム画像間において対応付けられた小領域毎に解析を行って、前記小領域毎の吸気の特徴量及び呼気の特徴量を算出し、前記算出された吸気の特徴量と呼気の特徴量との比の値を算出し、当該算出された比の値のヒストグラムを作成する解析手段と、
前記解析手段により作成されたヒストグラムを表示する表示手段と、
を備える胸部診断支援システム。
【請求項２】
前記吸気の特徴量が、吸気期に相当するフレーム画像群の時間的に隣接するフレーム画像間で画素信号値の差分を取ることにより得られるフレーム間差分値の絶対値の最大値であり、前記呼気の特徴量が、呼気期に相当するフレーム画像群の時間的に隣接するフレーム画像間で画素信号値の差分を取ることにより得られるフレーム間差分値の絶対値の最大値である請求項１に記載の胸部診断支援システム。
【請求項３】
前記吸気の特徴量が、呼吸１サイクルにおける吸気期に相当するフレーム画像の数であり、前記呼気の特徴量が、前記呼吸１サイクルにおける呼気期に相当するフレーム画像の数である請求項１に記載の胸部診断支援システム。
【請求項４】
前記吸気の特徴量が、呼吸１サイクル分の複数のフレーム画像における画素信号値の最大値であり、前記呼気の特徴量が、前記呼吸１サイクル分の複数のフレーム画像における画素信号値の最小値である請求項１に記載の胸部診断支援システム。
【請求項５】
前記領域分割手段は、前記複数のフレーム画像のうち一のフレーム画像の肺野領域を複数の小領域に分割し、他のフレーム画像の肺野領域を前記一のフレーム画像の各小領域と同じ画素位置の小領域に分割し、前記複数のフレーム画像間の前記同じ画素位置の各小領域を互いに対応付ける請求項１〜４の何れか一項に記載の胸部診断支援システム。
【請求項６】
前記領域分割手段は、前記複数のフレーム画像のうち一のフレーム画像の肺野領域を複数の小領域に分割し、他のフレーム画像のそれぞれにおいて前記一のフレーム画像における各小領域とマッチング度合いの高い領域を抽出し、前記他のフレーム画像の肺野領域を前記抽出された各領域に分割し、前記複数のフレーム画像間の前記マッチング度合いの高い各小領域を互いに対応付ける請求項１〜４の何れか一項に記載の胸部診断支援システム。
【請求項７】
前記撮影手段は、グリッドを使用せずに動態撮影を行う請求項１〜６の何れか一項に記載の胸部診断支援システム。
【請求項８】
コンピュータを、
少なくとも１サイクルの呼吸周期を含む胸部の動態を示す複数のフレーム画像のそれぞれから肺野領域を抽出する抽出手段、
前記抽出手段により抽出された肺野領域を複数の小領域に分割し、前記複数のフレーム画像間にわたり小領域を対応付ける領域分割手段、
前記複数のフレーム画像間において対応付けられた小領域毎に解析を行って、前記小領域毎の吸気の特徴量及び呼気の特徴量を算出し、前記算出された吸気の特徴量と呼気の特徴量との比の値を算出し、当該算出された比の値のヒストグラムを作成する解析手段、
前記解析手段により作成されたヒストグラムを表示する表示手段、
として機能させるためのプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、胸部診断支援システム及びプログラムに関する。
【背景技術】
【０００２】
従来のフィルム／スクリーンや輝尽性蛍光体プレートを用いた胸部の放射線による静止画撮影及び診断に対し、ＦＰＤ（flat panel detector）等の半導体イメージセンサを利用して胸部の動態画像を撮影し、診断に応用する試みがなされるようになってきている。具体的には、半導体イメージセンサの画像データの読取・消去の応答性の早さを利用し、半導体イメージセンサの読取・消去のタイミングと合わせて放射源からパルス状の放射線を連続照射し、１秒間に複数回の撮影を行って、胸部の動態を撮影する。撮影により取得された一連の複数枚の画像を順次表示することにより、医師は呼吸運動や心臓の拍動等に伴う胸部の一連の動きを観察することが可能となる。【０００３】
胸部の動態画像から診断に有用な情報を抽出するための各種技術も提案されている。例えば、特許文献１には、動態画像を構成する複数のフレーム画像間で画素値の差を示す差分画像を生成し、生成された差分画像から各画素値の最大値（平均値、中間値等）を求めて表示する技術が記載されている。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特開２００４−３１２４３４号公報
【発明の概要】
【発明が解決しようとする課題】
【０００５】
しかしながら、特許文献１に記載のように、フレーム間差分値の最大値（平均値、中間値等）を算出して表示するだけでは、熟練した呼吸器専門医以外の医師が病態（症例名）を把握することは難しい。【０００６】
本発明の課題は、肺野の換気に係る病態の把握を容易にし、診断に有効となる情報を医師に提供できるようにすることである。特に、聴診器による実務経験の少ない医師であっても、的確な診断を行える診断情報を提供することである。
【課題を解決するための手段】
【０００７】
上記課題を解決するため、本発明の第１の側面によると、
胸部診断支援システムは、
少なくとも１サイクルの呼吸周期を含む胸部の動態を撮影し、連続する複数のフレーム画像を生成する撮影手段と、
前記撮影手段により生成された複数のフレーム画像のそれぞれから肺野領域を抽出する抽出手段と、
前記抽出手段により抽出された肺野領域を複数の小領域に分割し、前記複数のフレーム画像間にわたり小領域を対応付ける領域分割手段と、
前記複数のフレーム画像間において対応付けられた小領域毎に解析を行って、前記小領域毎の吸気の特徴量及び呼気の特徴量を算出し、前記算出された吸気の特徴量と呼気の特徴量との比の値を算出し、当該算出された比の値のヒストグラムを作成する解析手段と、
前記解析手段により作成されたヒストグラムを表示する表示手段と、
を備える。【０００８】
前記吸気の特徴量が、吸気期に相当するフレーム画像群の時間的に隣接するフレーム画像間で画素信号値の差分を取ることにより得られるフレーム間差分値の絶対値の最大値であり、前記呼気の特徴量が、呼気期に相当するフレーム画像群の時間的に隣接するフレーム画像間で画素信号値の差分を取ることにより得られるフレーム間差分値の絶対値の最大値であることが好ましい。【０００９】
または、前記吸気の特徴量が、呼吸１サイクルにおける吸気期に相当するフレーム画像の数であり、前記呼気の特徴量が、前記呼吸１サイクルにおける呼気期に相当するフレーム画像の数であることが好ましい。【００１０】
または、前記吸気の特徴量が、呼吸１サイクル分の複数のフレーム画像における画素信号値の最大値であり、前記呼気の特徴量が、前記呼吸１サイクル分の複数のフレーム画像における画素信号値の最小値であることが好ましい。【００１１】
前記領域分割手段は、前記複数のフレーム画像のうち一のフレーム画像の肺野領域を複数の小領域に分割し、他のフレーム画像の肺野領域を前記一のフレーム画像の各小領域と同じ画素位置の小領域に分割し、前記複数のフレーム画像間の前記同じ画素位置の各小領域を互いに対応付けることが好ましい。【００１２】
前記領域分割手段は、前記複数のフレーム画像のうち一のフレーム画像の肺野領域を複数の小領域に分割し、他のフレーム画像のそれぞれにおいて前記一のフレーム画像における各小領域とマッチング度合いの高い領域を抽出し、前記他のフレーム画像の肺野領域を前記抽出された各領域に分割し、前記複数のフレーム画像間の前記マッチング度合いの高い各小領域を互いに対応付けることが好ましい。【００１３】
前記撮影手段は、グリッドを使用せずに動態撮影を行うことが好ましい。【００１４】
本発明の第２の側面によると、
プログラムは、
コンピュータを、
少なくとも１サイクルの呼吸周期を含む胸部の動態を示す複数のフレーム画像のそれぞれから肺野領域を抽出する抽出手段、
前記抽出手段により抽出された肺野領域を複数の小領域に分割し、前記複数のフレーム画像間にわたり小領域を対応付ける領域分割手段、
前記複数のフレーム画像間において対応付けられた小領域毎に解析を行って、前記小領域毎の吸気の特徴量及び呼気の特徴量を算出し、前記算出された吸気の特徴量と呼気の特徴量との比の値を算出し、当該算出された比の値のヒストグラムを作成する解析手段と、
前記解析手段により作成されたヒストグラムを表示する表示手段、
として機能させる。
【発明の効果】
【００１５】
本発明によれば、肺野の換気に係る解析結果（病態情報）を視認可能に提供するので、医師の経験の度合いに依存しない病態の把握が可能となり、病変部の早期発見に基づく早期治療が可能となる。
【図面の簡単な説明】
【００１６】
【図１】本発明の第１の実施の形態における胸部診断支援システムの全体構成を示す図である。
【図２】図１の撮影用コンソールの制御部により実行される撮影制御処理を示すフローチャートである。
【図３】図１の診断用コンソールの制御部により実行される画像解析処理を示すフローチャートである。
【図４】１つの呼吸サイクル（深呼吸時）において撮影された複数の時間位相Ｔ（Ｔ＝ｔ０〜ｔ６）のフレーム画像を示す図である。
【図５】安静呼気位と安静吸気位において肺野の同一部分を描画した領域の位置の変化を示す図である。
【図６】横隔膜の位置の算出方法を説明するための図である。
【図７】最大呼気位と最大吸気位において肺野の同一部分を描画した領域の位置の変化を示す図である。
【図８Ａ】第１の実施の形態において正常者の肺野の動態画像を解析した解析結果を表示した表示画面の一例を示す図である。
【図８Ｂ】第１の実施の形態においてＣＯＰＤ（閉塞性疾患）の肺野の動態画像を解析した解析結果を表示した表示画面の一例を示す図である。
【図８Ｃ】第１の実施の形態において混合性疾患の肺野の動態画像を解析した解析結果を表示した表示画面の一例を示す図である。
【図９】本発明の第２の実施の形態における放射線撮影システムの全体構成を示す図である。
【図１０】図９の保持部４５の近辺を拡大して示す図である。
【図１１】図９の保持部４５の正面を模式的に示す図である。
【図１２Ａ】図９の撮影用コンソールの制御部により実行される撮影制御処理を示すフローチャートである。
【図１２Ｂ】図９の撮影用コンソールの制御部により実行される撮影制御処理を示すフローチャートである。
【図１３】或る正常の肺野をグリッド有りとグリッド無しで撮影した動態画像の動態解析結果の比較を示す図である。
【図１４Ａ】第２の実施の形態において正常者の肺野の動態画像を解析した解析結果を表示した表示画面の一例を示す図である。
【図１４Ｂ】第２の実施の形態において、ＣＯＰＤ（閉塞性疾患）の肺野の動態画像を解析した解析結果を表示した表示画面の一例を示す図である。
【図１４Ｃ】第２の実施の形態において、混合性疾患の肺野の動態画像を解析した解析結果を表示した表示画面の一例を示す図である。
【図１５】「吸気の特徴量／呼気の特徴量」の傾向を示す指標値の表示例を示す図である。
【図１６】「吸気の特徴量／呼気の特徴量」の傾向を示す指標値の他の表示例を示す図である。
【図１７】「吸気の特徴量」又は「呼気の特徴量」の傾向を示す指標値の表示例を示す図である。
【発明を実施するための形態】
【００１７】
＜第１の実施の形態＞
以下、図面を参照して本発明の第１の実施の形態を詳細に説明する。ただし、発明の範囲は、図示例に限定されない。
【００１８】
〔胸部診断支援システム１００の構成〕
まず、構成を説明する。
図１に、本実施の形態における胸部診断支援システム１００の全体構成を示す。
図１に示すように、胸部診断支援システム１００は、撮影装置１と、撮影用コンソール２とが通信ケーブル等により接続され、撮影用コンソール２と、診断用コンソール３とがＬＡＮ（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）等を備え、通信ネットワークＮＴに接続された各装置との間のデータ送受信を制御する。
【００３０】
〔診断用コンソール３の構成〕
診断用コンソール３は、撮影用コンソール２から動態画像を取得し、取得した動態画像や後述するヒストグラム等の診断支援情報を表示して医師が読影診断するための動画像処理装置である。
診断用コンソール３は、図１に示すように、制御部３１、記憶部３２、操作部３３、表示部３４、通信部３５を備えて構成され、各部はバス３６により接続されている。
【００３１】
制御部３１は、ＣＰＵ、ＲＡＭ等により構成される。制御部３１のＣＰＵは、操作部３３の操作に応じて、記憶部３２に記憶されているシステムプログラムや、各種処理プログラムを読み出してＲＡＭ内に展開し、展開されたプログラムに従って、後述する画像解析処理を始めとする各種処理を実行し、診断用コンソール３各部の動作を集中制御する。制御部３１は、後述する画像解析処理を実行することにより、抽出手段、領域分割手段、解析手段としての機能を実現する。
【００３２】
記憶部３２は、不揮発性の半導体メモリやハードディスク等により構成される。記憶部３２は、制御部３１で画像解析処理を実行するための画像解析処理プログラムを始めとする各種プログラムやプログラムにより処理の実行に必要なパラメータ、或いは処理結果等のデータを記憶する。これらの各種プログラムは、読取可能なプログラムコードの形態で格納され、制御部３１は、当該プログラムコードに従った動作を逐次実行する。
【００３３】
操作部３３は、カーソルキー、数字入力キー、及び各種機能キー等を備えたキーボードと、マウス等のポインティングデバイスを備えて構成され、キーボードに対するキー操作やマウス操作により入力された指示信号を制御部３１に出力する。また、操作部３３は、表示部３４の表示画面にタッチパネルを備えても良く、この場合、タッチパネルを介して入力された指示信号を制御部３１に出力する。
【００３４】
表示部３４は、ＬＣＤやＣＲＴ等のモニタにより構成され、制御部３１から入力される表示信号の指示に従って、操作部３３からの入力指示やデータ等を表示する。
【００３５】
通信部３５は、ＬＡＮアダプタやモデムやＴＡ等を備え、通信ネットワークＮＴに接続された各装置との間のデータ送受信を制御する。
【００３６】
〔胸部診断支援システム１００の動作〕
次に、上記胸部診断支援システム１００における動作について説明する。
【００３７】
（撮影装置１、撮影用コンソール２の動作）
まず、撮影装置１、撮影用コンソール２による撮影動作について説明する。
図２に、撮影用コンソール２の制御部２１において実行される撮影制御処理を示す。撮影制御処理は、制御部２１と記憶部２２に記憶されている撮影制御処理プログラムとの協働により実行される。
【００３８】
まず、撮影技師により撮影用コンソール２の操作部２３が操作され、撮影対象（被写体Ｍ）の患者情報（患者の氏名、身長、体重、年齢、性別等）の入力が行われる（ステップＳ１）。
【００３９】
次いで、放射線照射条件が記憶部２２から読み出されて放射線照射制御装置１２に設定されるとともに、画像読取条件が記憶部２２から読み出されて読取制御装置１４に設定される（ステップＳ２）。
【００４０】
次いで、操作部２３の操作による放射線照射の指示が待機される（ステップＳ３）。ここで、撮影実施者は、撮影装置１において被写体Ｍのポジショニング等の撮影準備を行うとともに、安静呼吸の動態を撮影するために被験者（被写体Ｍ）に楽にするように指示し、安静呼吸を促す。撮影準備が整った時点で、操作部２３を操作して放射線照射指示を入力する。
【００４１】
操作部２３により放射線照射指示が入力されると（ステップＳ３；ＹＥＳ）、放射線照射制御装置１２及び読取制御装置１４に撮影開始指示が出力され、動態撮影が開始される（ステップＳ４）。即ち、放射線照射制御装置１２に設定されたパルス間隔で放射線源１１により放射線が照射され、放射線検出部１３によりフレーム画像が取得される。予め定められたフレーム数の撮影が終了すると、制御部２１により放射線照射制御装置１２及び読取制御装置１４に撮影終了の指示が出力され、撮影動作が停止される。撮影されるフレーム数は、少なくとも１呼吸サイクルが撮影できる枚数である。
【００４２】
撮影により取得されたフレーム画像は順次撮影用コンソール２に入力され、各フレーム画像に対して補正処理が行われる（ステップＳ５）。ステップＳ５の補正処理においては、オフセット補正処理、ゲイン補正処理、欠陥画素補正処理の３つの補正処理が行われる。まず最初に、取得された各フレーム画像に対してオフセット補正処理が行われ、取得された各フレーム画像に重畳された暗電流に起因するオフセット値が除去される。オフセット補正処理では、例えば、取得した各フレーム画像の各画素値（濃度値。以下、信号値という）から、予め記憶されたオフセット補正係数を減算する処理が行われる。ここで、オフセット補正係数は、予め放射線非照射時に取得した複数のフレーム画像を平均化した画像である。次いで、ゲイン補正処理が行われ、各フレーム画像の各画素に対応する各検出素子の個体差や読み出しアンプのゲインムラによって生じる画素毎のばらつきが除去される。ゲイン補正処理では、例えば、オフセット補正後の各フレーム画像に、予め記憶されたゲイン補正係数を乗算する処理が行われる。ここで、ゲイン補正係数は、放射線検出部１３に一様に放射線を照射した時に取得した複数のオフセット補正済みフレーム画像を平均化した画像と、このときの放射線照射条件で期待される出力信号値の関係から、補正後の各画素の信号値が一様となるように予め算出され、記憶された係数である。次いで、欠陥画素補正処理が行われ、周囲の画素と比較して感度が非線形な画素や、感度がない欠落画素が除去される。欠陥画素補正処理では、例えば、予め記憶された欠陥画素位置情報マップに従って、欠陥画素位置情報マップに登録された各欠陥画素において、欠陥画素の信号値をその近傍の欠陥でない画素の信号値の平均値で置き換える処理が行われる。ここで、欠陥画素位置情報マップは、放射線検出部１３に一様に放射線を照射した時に取得したオフセット補正、ゲイン補正済みのフレーム画像から、予め複数の欠陥画素が認識され、その欠陥画素の位置が登録されたマップである。上記オフセット補正係数及びゲイン補正係数、欠陥画素位置情報マップは、ビニングやダイナミックレンジ等の収集モードに応じて、予め、それぞれ最適な値が記憶されており、それぞれの収集モードにおいて対応する最適な値が読み出されるようになっている。
【００４３】
次いで、補正処理後の各フレーム画像と、撮影順を示す番号と対応付けて記憶部２２に記憶されるとともに（ステップＳ６）、表示部２４に表示される（ステップＳ７）。ここで、各フレーム画像を記憶する直前に、各フレーム画像の各画素の信号値を真数から対数に変換する対数変換処理を行ってから記憶しても良い。撮影技師は、表示された動態画像によりポジショニング等を確認し、撮影により診断に適した画像が取得された（撮影ＯＫ）か、再撮影が必要（撮影ＮＧ）か、を判断する。そして、操作部２３を操作して、判断結果を入力する。尚、撮影により取得された各フレーム画像は、全撮影の終了後に纏めて入力するようにしても良い。
【００４４】
操作部２３の所定の操作により撮影ＯＫを示す判断結果が入力されると（ステップＳ８；ＹＥＳ）、動態撮影で取得された一連のフレーム画像のそれぞれに、動態画像を識別するための識別ＩＤや、患者情報、検査対象部位、放射線照射条件、画像読取条件、撮影順を示す番号、撮影日時等の情報が付帯され（例えば、ＤＩＣＯＭ形式で画像データのヘッダ領域に書き込まれ）、通信部２５を介して診断用コンソール３に送信される（ステップＳ９）。そして、本処理は終了する。一方、操作部２３の所定の操作により撮影ＮＧを示す判断結果が入力されると（ステップＳ８；ＮＯ）、記憶部２２に記憶された一連のフレーム画像が削除され（ステップＳ１０）、本処理は終了する。尚、このケースに於いては、再撮影が実行されることとなる。
【００４５】
（診断用コンソール３の動作）
次に、診断用コンソール３における動作について説明する。
診断用コンソール３においては、通信部３５を介して撮影用コンソール２から動態画像の一連のフレーム画像が受信されると、制御部３１と記憶部３２に記憶されている画像解析処理プログラムとの協働により図３に示す画像解析処理が実行される。
【００４６】
以下、図３を参照して画像解析処理の流れについて説明する。
まず、各フレーム画像から肺野領域が抽出される（ステップＳ２１）。
肺野領域の抽出方法は何れの方法であってもよい。例えば、一連のフレーム画像中の任意のフレーム画像（ここでは、撮影順が一番（最初）のフレーム画像とする。）の各画素の信号値（濃度値）のヒストグラムから判別分析によって閾値を求め、この閾値より高信号の領域を肺野領域候補として１次抽出する。次いで、１次抽出された肺野領域候補の境界付近でエッジ検出を行い、境界付近の小領域でエッジが最大となる点を境界に沿って抽出すれば肺野領域の境界を抽出することができる。
【００４７】
次いで、各フレーム画像の肺野領域が複数の小領域に分割され、各フレーム画像の小領域が互いに対応付けられる（ステップＳ２２）。各小領域の画素の位置は、制御部３１のＲＡＭに記憶される。
【００４８】
ここで、呼吸サイクルは、呼気期と吸気期により構成される。図４は、１つの呼吸サイクル（深呼吸時）において撮影された複数の時間位相Ｔ（Ｔ＝ｔ０〜ｔ６）のフレーム画像を示す図である。図４に示すように、呼気期は、横隔膜が上がることによって肺から空気が排出され、肺野の領域が小さくなる。最大呼気位では、横隔膜の位置が最も高い状態となる。吸気期は、横隔膜が下がることにより肺に空気が取り込まれ、図４に示すように胸郭中の肺野の領域が大きくなる。最大吸気位では、横隔膜の位置が最も下がった状態となる。即ち、肺野領域の同一部分の位置は呼吸運動に応じて時間とともに変化するため、各フレーム画像間で肺野の同一の部分（特に、下部領域（横隔膜付近））を示す画素位置がずれることになる。
しかし、安静呼吸時に撮影された画像においては、上記の位置ずれは小さく、後述する解析結果が狂ってしまうほどの位置ずれはおきない。図５の画像Ｄ１は、安静呼気位（安静呼吸時において横隔膜の位置が最も高くなったタイミング）のフレーム画像である。図５の画像Ｄ２は、安静吸気位（安静呼吸時において横隔膜の位置が最も低くなったタイミング）のフレーム画像である。即ち、図５の画像Ｄ１とＤ２とは、呼吸１サイクルにおいて最も形状の差の大きいタイミングで撮影された画像である。しかし、図５の画像Ｄ１、Ｄ２間においては、最も位置ずれの大きい肺野領域の下部領域においても位置ずれはわずかであることがわかる（画像Ｄ２のＡ１１は画像Ｄ１のＡ１と同じ画素位置を示し、画像Ｄ２のＡ２は画像Ｄ１のＡ１と肺野における同一部分を描画した領域を示している）。
【００４９】
そこで、ステップＳ２２における具体的な処理としては、まず、一連のフレーム画像の中から一のフレーム画像を基準画像として設定する。次いで、基準画像の抽出された肺野領域を複数の小領域（例えば、２ｍｍ×２ｍｍの矩形領域）に分割する（図５参照）。次いで、他のフレーム画像の肺野領域を基準画像の各小領域と同じ画素位置の小領域（放射線検出部１３の同じ撮像素子から出力される信号値の領域）に分割する。次いで、各フレーム画像間の同じ画素位置の各小領域を互いに対応付ける。この処理では、フレーム画像の小領域への分割及び対応付けを高速に行うことが可能となる。
基準画像としては、安静呼気位のフレーム画像とすることが好ましい。安静呼気位では、安静呼吸時において横隔膜の位置が最も高くなる、即ち、肺野領域の面積が最も小さくなるので、基準画像の小領域を他のフレーム画像に対応付けたときに、小領域が他のフレーム画像の肺野外の領域に対応付けられることがないからである。
安静呼気位の画像は、一連のフレーム画像の中から横隔膜の位置が最も高い位置にある画像を抽出することで取得することができる。横隔膜の位置は、例えば、図６に示す横隔膜の基準位置Ｂ１を横隔膜のカーブＣ（図６に点線で示す）の垂直方向の平均位置と予め定義しておき、肺野領域Ｒから横隔膜のカーブＣ（肺野領域の下端）を抽出し、その垂直方向の平均位置を求め、求めた位置を横隔膜の基準位置Ｂ１として特定する。
【００５０】
なお、撮影された動態画像が深呼吸時の画像である場合は、図７に示すように、肺野の同一の部分を示す画素位置が大幅にずれる。ここで、図７の画像Ｄ３は、深呼吸時の最大呼気位のフレーム画像であり、図７の画像Ｄ４は、深呼吸時の最大吸気位のフレーム画像である（画像Ｄ４のＡ３１は画像Ｄ３のＡ３と同じ画素位置を示し、画像Ｄ４のＡ４は画像Ｄ３のＡ３と肺野における同一部分を描画した領域を示している）。そのため、安静呼吸時と同様に各フレーム画像における、基準画像の各小領域と同じ画素位置の領域をその小領域に対応する領域としてしまうと、後述する解析結果が診断に利用し得ないものとなる。そこで、このような場合には、各フレーム画像間の対応点を抽出する対応点抽出処理（ローカルマッチング処理）及び非線形歪変換処理（ワーピング処理）を行って、肺野領域の同一部分が描画された領域を各フレーム画像間で対応付ける。また、安静呼吸時に於いて、処理速度をあまり重視せず、解析精度をより高めたい場合にも、これらの処理を行うこととしてもよい。
【００５１】
ローカルマッチング処理では、まず、撮影順が１番（最初）のフレーム画像から抽出した肺野領域を、例えば、２ｍｍ角の矩形からなる小領域に分割する。
次いで、撮影順が１番のフレーム画像をＰ１、これと隣接するフレーム画像（撮影順が隣接するフレーム画像（即ち、時間的に隣接するフレーム画像。以下同様。））をＰ２とし、Ｐ２に、Ｐ１の各小領域の探索領域を設定する。ここで、Ｐ２の探索領域は、Ｐ１における各小領域における中心点の座標を（ｘ，ｙ）とすると、同一の中心点（ｘ，ｙ）をもち、Ｐ１の小領域よりも縦横の幅が大きくなるように設定する（例えば、１．５倍）。そして、Ｐ１の各領域毎に、Ｐ２の探索範囲で最もマッチング度合いが高くなる領域を求めることで、Ｐ１の各小領域に対するＰ２上での対応位置を算出する。マッチング度合いとしては、最小二乗法や相互相関係数を指標に用いる。そして、Ｐ２の肺野領域をＰ１の各小領域の対応位置で分割する。
次いで、Ｐ２を、新たにＰ１とみなし、撮影順がＰ２の次のフレーム画像を新たなＰ２とみなして、Ｐ１の各小領域におけるＰ２の対応位置を算出する。以上の処理を繰り返すことで、各フレーム画像の各小領域が隣接するフレーム画像のどの位置に対応するかが求まる。求めた処理結果は、制御部３１のＲＡＭに記憶される。
【００５２】
次いで、ワーピング処理が行われる。具体的には、撮影順が１番のフレーム画像をＰ１、これと撮影順が隣接する（時間的に隣接する）フレーム画像をＰ２とし、上記ローカルマッチング処理で算出された隣接するフレーム画像間の各小領域の対応位置に基づいて、各小領域毎にＰ１からＰ２へのシフトベクトルを算出する。次いで、算出されたシフトベクトルを多項式でフィッティングして、この多項式を用いて各小領域における各画素のシフトベクトルを算出する。そして、算出された各画素のシフトベクトルに基づいて、ワーピング処理を行い、Ｐ２の各小領域内の各画素の位置をＰ１のフレーム画像の対応する画素の位置にシフトする。次いで、ワーピング処理されたＰ２を、新たにＰ１とみなし、撮影順がＰ２の次のフレーム画像を新たなＰ２とみなして、上記処理を行う。以上の処理を撮影順の早い隣接フレーム画像間から順次繰り返すことで、全てのフレーム画像の各小領域の位置を撮影順が１のフレーム画像（この処理での基準画像）に略一致させることが可能となる。各フレーム画像間における小領域の位置の対応関係は、制御部３１のＲＡＭに記憶される。
【００５３】
次いで、一連のフレーム画像のステップＳ２２で対応付けられた各小領域毎に解析が行われ、吸気の特徴量及び呼気の特徴量がそれぞれ算出される（ステップＳ２３）。ここでは、呼気の特徴量及び吸気の特徴量として、例えば、呼気期及び吸気期のそれぞれにおける各小領域毎のフレーム間差分値（微分値）の代表値（絶対値の最大値）を算出する。フレーム間差分値は、そのフレーム画像が撮影されたタイミングにおける信号変化量を示す値である。呼吸により息を吸ったり吐いたりすれば、その息の流れに応じて肺の密度が変化し、これによってＸ線透過量（つまり、画素の出力信号値）が変化する。よって、信号変化量は、そのタイミングにおける気流速度を示す値とみなすことができる。なお、代表値としては、絶対値の最大値に限らず、中央値、平均値、最頻値としてもよい。
【００５４】
具体的には、まず、各フレーム画像の各小領域内の画素の信号値（平均信号値）が算出される。次いで、撮影順が隣接するフレーム画像間で各小領域の信号値の差分を算出するフレーム間差分処理が行われる。ここでは、各小領域毎に、フレーム番号ＮとＮ＋１（Ｎは１、２、３・・・）のフレーム画像について、Ｎ＋１−Ｎの差分値が算出される。そして、呼気期のフレーム間差分値の最大値（絶対値の最大値）が呼気の特徴量として取得され、吸気期のフレーム間差分値の最大値（絶対値の最大値）が吸気の特徴量として取得される。フレーム間差分値の最大値（絶対値の最大値）は、最大微分値に相当する。ここで、各小領域内のフレーム間差分値の符号が正の期間が吸気期、負の期間が呼気期である。
【００５５】
次いで、各小領域毎の吸気の特徴量と呼気の特徴量の比の値（吸気の特徴量／呼気の特徴量）が算出される（ステップＳ２４）。ここでは、「吸気期のフレーム間差分値の最大値／呼気期のフレーム間差分値の最大値」が算出される。
【００５６】
次いで、算出された各小領域毎の「吸気の特徴量／呼気の特徴量」の値のヒストグラムが作成されるとともに、肺野全体における「吸気の特徴量／呼気の特徴量」の傾向を示す指標値（ここでは、平均値、標準偏差）が算出される（ステップＳ２５）。ヒストグラムの縦軸のカウント数は、肺野内の全小領域の数で除算することにより正規化することが好ましい。
【００５７】
次いで、各小領域について求められた「吸気の特徴量／呼気の特徴量」の値を、予め記憶部３２に記憶されている、「吸気の特徴量／呼気の特徴量」の値と表示時のパラメータ値との変換テーブルに基づいて表示用パラメータ値に変換し、変換したパラメータ値にて、基準画像（例えば、安静呼気位のフレーム画像）の各小領域を表示した画像が作成される（ステップＳ２６）。変換テーブルは、例えば、正常／異常（重症度１〜ｎ）の各カテゴリに特徴量を分類するときの各カテゴリの特徴量の大きさの範囲を規定する閾値（各カテゴリの閾値）と色相、明度、輝度、透明度の何れかを一対一に対応付けたテーブルである。ここで、表示時のパラメータ値の変換テーブルとしては、特徴量の大きさに対する識別能を高めるために、各カテゴリの閾値と色相を対応付けることが好ましい。
このとき、例えば、数個（例えば、５〜６個）の色相を上記各カテゴリの閾値と対応させ、さらにその間の特徴量の値に対しては中間の色相を割り当てる（グラデーションさせる）ことで識別能の高い表示とすることができる。
表示用パラメータ値に基づいて色付けされた画像は、基準画像のフレーム画像にオーバーレイ表示してもよい。
【００５８】
そして、表示部３４に、作成されたヒストグラム、及び作成された静止画像等が並列表示され（ステップＳ２７）、画像解析処理は終了する。ヒストグラムの領域は、上述の「吸気の特徴量／呼気の特徴量」の値と表示時のパラメータ値との変換テーブルに基づいて、基準画像の肺野領域の各小領域と同様の基準により色付けして表示される。
【００５９】
図８Ａ〜図８Ｃに、ステップＳ２７において表示部３４に表示される表示画面の例を示す。
図８Ａは、正常者の肺野の動態画像を解析した解析結果を表示した表示画面である。図８Ｂは、ＣＯＰＤ（閉塞性疾患）の肺野の動態画像を解析した解析結果を表示した表示画面である。図８Ｃは、混合性疾患の肺野の動態画像を解析した解析結果を表示した表示画面である。
【００６０】
図８Ａ〜図８Ｃに示すように、ステップＳ２７においては、半切サイズのフレーム画像から抽出された肺野領域内における各小領域（２ｍｍ角の矩形サイズ）の「吸気の特徴量／呼気の特徴量」の値のヒストグラム３４ａと、各小領域の「吸気の特徴量／呼気の特徴量」を一覧表示した静止画像３４ｂと、ヒストグラム３４ａ及び静止画像３４ｂで表示されている色相と「吸気の特徴量／呼気の特徴量」の値との関係を示す表示３４ｃと、肺野全体での「吸気の特徴量／呼気の特徴量」の傾向を示す指標値３４ｄが表示されている。「吸気の特徴量／呼気の特徴量」の分布の特徴を把握することが可能となる。また、図８Ａ〜図８Ｃに示すように、ヒストグラム３４ａは、「吸気の特徴量／呼気の特徴量」の値の大きさに応じて、横軸の領域が５つの色相で色分けして表示される。これにより、医師は、ヒストグラムを一瞥するだけで容易に肺野内における「吸気の特徴量／呼気の特徴量」の分布を把握することができる。また、各小領域の「吸気の特徴量／呼気の特徴量」を示す静止画像３４ｂにおいて、各小領域が「吸気の特徴量／呼気の特徴量」の値に応じてヒストグラムの色分けと同じ基準で色分けされて表示されるので、医師が肺野内の局所的な異常箇所（閉塞性部分、拘束性部分）を容易に把握することが可能となる。更に、肺野全体での「吸気の特徴量／呼気の特徴量」の傾向を示す指標値３４ｄとして、その平均値及び標準偏差を算出して画面上に併せて表示することで、肺野全体での「吸気の特徴量／呼気の特徴量」の傾向を数値として医師に提供することができる。
【００６１】
ここで、呼気期のフレーム間差分値の最大値（絶対値）を呼気の特徴量、吸気期のフレーム間差分値の最大値（絶対値）を吸気の特徴量とした場合、正常者の肺野では、肺野全体の「吸気の特徴量／呼気の特徴量」の平均値は０．９〜１．２、標準偏差は０．１０〜０．２２程度であることが知られている。よって、ステップＳ２７で図８Ａに示す表示画面が表示された場合、医師は、撮影した肺野が正常であることを容易に把握することができる。
一方、ＣＯＰＤ（閉塞性疾患）の肺野では、肺野全体の「吸気の特徴量／呼気の特徴量」の平均値は０．９〜１．２に収まらず（正常者に比べて大きく）、標準偏差も、０．１０〜０．２２に収まらない（正常者に比べて大きい）ことが知られている。よって、ステップＳ２７で図８Ｂに示す表示画面が表示された場合、医師は、撮影した肺野がＣＯＰＤであることを容易に把握することができる。
一方、混合性肺疾患の肺野では、肺野全体の「吸気の特徴量／呼気の特徴量」の値が０．６６以下のデータ数、及び１．５以上のデータ数がともに増加することが知られている。よって、ステップＳ２７で図８Ｃに示す表示画面が表示された場合、医師は、撮影した肺野が混合性疾患であることを容易に把握することができる。
このように、胸部診断支援システム１００においては、ＣＯＰＤ（閉塞性肺疾患）、混合性肺疾患等の病態を特定できるような、有用な診断情報を医師に提供することが可能となる。
【００６２】
呼気の特徴量及び吸気の特徴量としては、上記の例以外の他の特徴量を用いても良い。
例えば、呼気の特徴量を呼吸１サイクルにおける呼気期に相当するフレーム画像数（呼気時間）、吸気の特徴量を呼吸１サイクルにおける吸気期に相当するフレーム画像数（吸気期間）としてもよい。ここで、肺の換気機能が正常である場合、吸気時間と呼気時間はほぼ同じ長さ、もしくは若干呼気時間が長くなる筈である。よって、「呼気期に相当するフレーム画像数／吸気期に相当するフレーム画像数」の値を見れば、医師は肺疾患の疑いがあるか否かを把握することができる。特に、「呼気期のフレーム画像数／吸気期のフレーム画像数」＞１．５の領域は、呼気で換気がしにくく、取り込んだ空気を排出するのが遅れてしまう閉塞性部分であることがわかる。なお、「吸気期のフレーム間差分値の最大値／呼気期のフレーム間差分値の最大値≒呼気時間(呼気のフレーム画像数)／吸気時間(吸気のフレーム画像数)」の関係が成り立つため、医師は、呼気の特徴量を呼気期のフレーム間差分値の最大値、吸気の特徴量を吸気期のフレーム間差分値の最大値とした場合と同様の判断基準により、正常、ＣＯＰＤ(閉塞性肺疾患)、混合性肺疾患の識別を行うことが可能である。
【００６３】
また、１呼吸サイクル分の各フレーム画像において、各小領域の画素の信号値（平均信号値）を算出し、各小領域毎に呼吸１サイクルにおける信号値の最小値及び最大値を求め、求めた最小値をその領域の呼気の特徴量、最大値を吸気の特徴量としてもよい。信号値の最大値と最小値は正常な箇所では両者の値の差が大きくなり、異常のある箇所では両者の差が非常に小さくなると考えられる。異常のある箇所では肺胞の動きが悪くなるため、肺胞の密度変化が小さくなるからである。よって、医師は、「信号値の最大値／信号値の最小値」のヒストグラムを参照し、平均値及び標準偏差を確認することで、肺野が正常であるか疾患であるかの判断材料とすることができる。例えば、肺野全体の「信号値の最大値／信号値の最小値」の平均値が１より大きく、標準偏差が小さい場合、肺の機能は正常であると判断することができる。一方、肺野全体の「信号値の最大値／信号値の最小値」の平均値が１に近い値であり、標準偏差が大きい場合、肺の機能に疾患があると判断することができる。
【００６４】
その他、「吸気の特徴量／呼気の特徴量」の傾向を示す指標値として、平均値、標準偏差の他に、カウント数（ブロック（小領域）数）がヒストグラムのピークとなる「吸気の特徴量／呼気の特徴量」の値、もしくは、ピークのカウント数（ブロック数）、もしくは、「吸気の特徴量／呼気の特徴量」の値が所定以上所定以下となるカウント数の割合を指標値として用いても良い。もしくは、これらの複数の指標値を組み合わせて新たな指標値を作成しても良い。また、「吸気の特徴量／呼気の特徴量」を対数に変換し、対数変換後の値に対して、上記ヒストグラムの作成、指標値の算出、表示時パラメータとの変換テーブルに基づく色付け画像の作成を行っても良い。
【００６５】
以上説明したように、胸部診断支援システム１００の診断用コンソール３によれば、制御部３１は、胸部の動態を撮影して得られた複数のフレーム画像のそれぞれから肺野領域を抽出し、抽出した肺野領域を複数の小領域に分割し、複数のフレーム画像間における各小領域を対応付ける。次いで、複数のフレーム画像の対応する小領域毎に解析を行って、小領域毎の吸気の特徴量及び呼気の特徴量を算出し、算出された吸気の特徴量と呼気の特徴量との比の値を算出し、当該算出された比の値のヒストグラムを作成する。そして、作成したヒストグラムを表示部３４に表示する。
【００６６】
従って、吸気の特徴量と呼気の特徴量の比が肺野内においてどのように分布しているのかを示す、診断に有利となる情報であるヒストグラムを医師に提供することができるので、医師は肺野の換気に係る病態を容易に把握することが可能となる。
【００６７】
吸気の特徴量として、吸気期に相当するフレーム画像群の時間的に隣接するフレーム画像間のフレーム間差分値の絶対値の最大値を算出し、呼気の特徴量として、呼気期に相当するフレーム画像群の時間的に隣接するフレーム画像間のフレーム間差分値の絶対値の最大値を算出し、算出された吸気の特徴量と呼気の特徴量の比の分布をヒストグラムで表示することで、撮影された肺野が正常であるか、閉塞性疾患であるか、混合性疾患であるかを医師が判断する際の有用な情報を提供することが可能となる。
【００６８】
吸気の特徴量として、呼吸１サイクルにおける吸気期に相当するフレーム画像の数を算出し、呼気の特徴量として、呼吸１サイクルにおける呼気期に相当するフレーム画像の数を算出し、算出された吸気の特徴量と呼気の特徴量の比の分布をヒストグラムで表示することで、撮影された肺野が正常であるか、閉塞性疾患であるか、混合性疾患であるかを医師が判断する際の有用な情報を提供することが可能となる。
【００６９】
吸気の特徴量として、呼吸１サイクル分の複数のフレーム画像における画素信号値の最大値を算出し、呼気の特徴量として、呼吸１サイクル分の複数のフレーム画像における画素信号値の最小値を呼気の特徴量として算出し、算出された吸気の特徴量と呼気の特徴量の比の分布をヒストグラムで表示することで、撮影された肺野が正常であるか異常であるかを医師が判断する際の有用な情報を提供することができる。
【００７０】
複数のフレーム画像の肺野領域を小領域に分割する際、複数のフレーム画像のうち一のフレーム画像の肺野領域を複数の小領域に分割し、他のフレーム画像の肺野領域を一のフレーム画像の各小領域と同じ画素位置の小領域に分割し、複数のフレーム画像間の同じ画素位置の各小領域を互いに対応付けることで、フレーム画像の肺野の小領域への分割及びフレーム画像間での小領域の対応付けを高速に行うことができる。その結果、解析処理の高速化を図ることが可能となる。
【００７１】
複数のフレーム画像の肺野領域を小領域に分割する際、複数のフレーム画像のうち一のフレーム画像の肺野領域を複数の小領域に分割し、他のフレーム画像のそれぞれにおいて一のフレーム画像における各小領域とマッチング度合いの高い領域を抽出し、他のフレーム画像の肺野領域を抽出された各領域に分割し、複数のフレーム画像間のマッチング度合いの高い各小領域を互いに対応付けることで、精度良く解析を行うことが可能となる。
【００７２】
＜第２の実施の形態＞
以下、図面を参照して本発明の第２の実施の形態を詳細に説明する。ただし、発明の範囲は、図示例に限定されない。
【００７３】
ここで、従来、胸部静止画像の放射線撮影においては、例えば、特開平１０−９８５８７号公報（公知文献１）に記載のように、散乱放射線を除去するための散乱線除去グリッド（以下、グリッドという）を使用して撮影が行われている。
また、近年、特開２００５−３１２７７５号公報（公知文献２）に記載のように、ＦＰＤ（flat panel detector）等の放射線検出器を利用して胸部の動態撮影を行い、診断に応用するシステムが提案されているが、このような動態撮影用のシステムにおいても、静止画撮影と同様に、グリッドを使用して撮影を行う構成となっている。
【００７４】
公知文献２においては、被写体の被曝線量が増大しがちな動態撮影に於いて、被写体の被曝線量低減のために、静止画撮影１回当たりの照射線量と一連の動態撮影に係るトータルの照射線量とが略一致するように、動態撮影におけるフレーム画像１枚当たりの照射線量を算出している。ここで、従来の手法に従いグリッドを使用すると放射線検出器（検出器）へ到達する線量が減弱されるため、公知文献２に記載のような放射線照射条件で動態撮影を行うと、各フレーム画像に到達する線量が大幅に不足し、Ｓ／Ｎ比が低下する。そのため、一連のフレーム画像を用いて動態解析を行った場合、正確な解析ができず、照射線量を変更して再撮影を行うこととなり、結果として被写体のトータルの被曝線量が増大してしまう恐れもある。
【００７５】
そこで、本件発明者等が鋭意検討を行った結果、検出器の絶対的出力信号（濃度信号）を観察する静止画撮影とは異なり、動態撮影に於いては、各フレーム相互の相対的出力信号差（濃度差）を視認したり、例えばフレーム間差分等の相対的出力信号差に基づく特徴量の解析を行ったりするので、各フレーム画像に略一様に散乱線影響が重畳していても、その目的を達成することが可能であることを見出した。
従って、動態撮影においては、グリッド無しで撮影を行うことができる。グリッド無しでの撮影において、検出器に到達する線量を従来並みに保つ場合、被写体への照射線量はグリッドによる減衰分を予め加算しておく必要が無く、従来手法に比べ、被写体の被曝線量を半減することが可能となる。また、被写体の被曝線量（照射線量）を一定のままとすると、各フレーム画像あたりの検出器到達線量を増加することができるので、より精度を有する解析を行うことができる。
【００７６】
第２の実施の形態においては、動態撮影時にグリッド無しで撮影を行う放射線撮影システム２００について説明する。
【００７７】
〔放射線撮影システム２００の構成〕
まず、構成を説明する。
図９に、本実施の形態における放射線撮影システム２００の全体構成を示す。この放射線撮影システム２００は、胸部診断支援システムとして機能するシステムである。
図９に示すように、放射線撮影システム２００は、撮影装置４と、撮影用コンソール５とが通信ケーブル等により接続され、撮影用コンソール５と、診断用コンソール６とがＬＡＮ（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）等を備え、通信ネットワークＮＴ２に接続された各装置との間のデータ送受信を制御する。
【００９１】
〔診断用コンソール６の構成〕
診断用コンソール６は、撮影用コンソール５から静止画像、又は動態画像の一連のフレーム画像を取得し、取得した画像や後述するヒストグラム等の診断支援情報を表示して医師が読影診断するためのコンピュータ装置である。
診断用コンソール６は、図９に示すように、制御部６１、記憶部６２、操作部６３、表示部６４、通信部６５を備えて構成され、各部はバス６６により接続されている。
【００９２】
制御部６１は、ＣＰＵ、ＲＡＭ等により構成される。制御部６１のＣＰＵは、操作部６３の操作に応じて、記憶部６２に記憶されているシステムプログラムや、各種処理プログラムを読み出してＲＡＭ内に展開し、展開されたプログラムに従って、後述する画像解析処理を始めとする各種処理を実行し、診断用コンソール６各部の動作を集中制御する。
【００９３】
記憶部６２は、不揮発性の半導体メモリやハードディスク等により構成される。記憶部６２は、制御部６１で画像解析処理を実行するための画像解析処理プログラムを始めとする各種プログラムやプログラムにより処理の実行に必要なパラメータ、或いは処理結果等のデータを記憶する。これらの各種プログラムは、読取可能なプログラムコードの形態で格納され、制御部６１は、当該プログラムコードに従った動作を逐次実行する。
【００９４】
操作部６３は、カーソルキー、数字入力キー、及び各種機能キー等を備えたキーボードと、マウス等のポインティングデバイスを備えて構成され、キーボードに対するキー操作やマウス操作により入力された指示信号を制御部６１に出力する。また、操作部６３は、表示部６４の表示画面にタッチパネルを備えても良く、この場合、タッチパネルを介して入力された指示信号を制御部６１に出力する。
【００９５】
表示部６４は、ＬＣＤやＣＲＴ等のモニタにより構成され、制御部６１から入力される表示信号の指示に従って、操作部６３からの入力指示やデータ等を表示する。
【００９６】
通信部６５は、ＬＡＮアダプタやモデムやＴＡ等を備え、通信ネットワークＮＴ２に接続された各装置との間のデータ送受信を制御する。
【００９７】
〔放射線撮影システム２００の動作〕
次に、上記放射線撮影システム２００における動作について説明する。
【００９８】
（撮影装置４、撮影用コンソール５の動作）
まず、撮影装置４、撮影用コンソール５による撮影動作について説明する。
図１２Ａ〜図１２Ｂに、撮影用コンソール５の制御部５１において実行される撮影制御処理を示す。撮影制御処理は、制御部５１と記憶部５２に記憶されている撮影制御処理プログラムとの協働により実行される。
【００９９】
まず、撮影実施者により撮影用コンソール５の操作部５３が操作され、撮影対象（被写体Ｍ）の患者情報（患者の氏名、身長、体重、年齢、性別等）及び撮影種別（動態撮影、静止画撮影の区別）の入力が行われる（ステップＳ３１）。この他、通信ネットワークＮＴ２に接続された図示しない他装置から送信され、通信部５５を経由して受信したデータも、同様の入力情報とすることができる。
【０１００】
次いで、入力された撮影種別が動態撮影であるか静止画撮影であるかが判断される（ステップＳ３２）。入力された撮影種別が動態撮影であると判断されると（ステップＳ３２；ＹＥＳ）、動態撮影用の放射線照射条件が記憶部５２から読み出されて放射線照射制御装置４２に設定されるとともに、動態撮影用の画像読取条件が記憶部５２から読み出されて読取制御装置４４に設定される（ステップＳ３３）。入力された撮影種別が静止画撮影であると判断されると（ステップＳ３２；ＮＯ）、静止画撮影用の放射線照射条件が記憶部５２から読み出されて放射線照射制御装置４２に設定されるとともに、静止画撮影用の画像読取条件が記憶部５２から読み出されて読取制御装置４４に設定される（ステップＳ３７）。
【０１０１】
本実施の形態においては、動態撮影用としてグリッド４６を装着しない状態で（グリッド無しで）の撮影を前提とした放射線照射条件が設定され、静止画撮影用としてグリッド４６を装着した状態で（グリッド有りで）の撮影を前提とした放射線照射条件が設定される。具体的には、動態撮影用の放射線照射条件は、フレーム画像１枚ごとの撮影においてグリッド４６を使用しないで放射線検出部４３に到達する線量が静止画撮影においてグリッド４６を使用して放射線検出部４３に到達する線量と同じになるような放射線照射条件が設定される。即ち、動態撮影では、１枚あたりの照射線量が静止画撮影に比べて低くなるように放射線照射条件が設定される。
【０１０２】
ここで、人体等の被写体Ｍに対して放射線を照射して放射線撮影を行う場合、体内を透過する放射線が体内組織で散乱される。そして、このような散乱放射線が放射線検出部４３に入射すると放射線画像にノイズが発生する。そこで、病変部分の検出や病変部分の観察等に用いられ個々の画素の絶対的出力値（信号値）が重要な静止画の撮影においては、放射線検出部４３の被写体側の面、すなわち被写体を透過した放射線が入射する側の面に、グリッド４６を設けて撮影が行われることが好ましい。グリッド有りで撮影を行う場合、グリッド４６により放射線検出部４３に到達する線量が減弱される（例えば、露出倍数２のグリッドでは約１／２に減弱される）ため、この減弱分を加算した分の放射線を照射する必要がある。従来、動態撮影においても、静止画撮影と同様にグリッドを用いた撮影が行われていた。
【０１０３】
しかしながら、動態撮影の場合は、静止画撮影に比べて撮影するフレーム画像枚数が多いため、１枚当たりに放射線源４１から照射する照射線量を静止画撮影と同等とした場合、被写体Ｍの被曝量が多くなってしまうという問題がある。被写体Ｍの被曝量低減のため１回の静止画撮影と、一連の動態撮影のトータルの照射線量を同じとする技術も開示されているが、この場合は個々のフレーム画像の線量が不足気味となり、Ｓ／Ｎ比が低下してしまう。
【０１０４】
そこで、本件発明者等は、鋭意検討を重ねた結果、動態撮影により得られる胸部動態画像は主として換気機能や血流等の動態解析に使用されるものであり、これらの解析においては、グリッドを配置せずに撮影を行った動態画像を用いても、グリッド有りで撮影を行った動態画像を用いた場合と略同等の結果を得ることができることを見出した。言い換えると、グリッド有無に限らず、放射線検出器に到達する線量が同じであれば、略同等の解析結果を得ることができる（図１３参照。詳細後述）。
【０１０５】
本実施の形態の放射線撮影システム２００においては、この知見に基づき、静止画撮影ではグリッド４６を装着して（グリッド有りで）撮影を行い、動態撮影ではグリッド４６を装着せずに（グリッド無しで）撮影を行う。そして、グリッド４６を使用した静止画撮影時と、グリッド４６無しでの動態撮影の１フレーム当たりの撮影とで放射線検出部４３に到達する線量が略同等となる放射線照射条件で撮影を行うことにより、放射線検出部４３への到達線量は従来と略同等に維持しつつ被写体Ｍの被曝線量の低減する構成となっている。または、１回の静止画撮影と、一連の動態撮影のトータルの照射線量を同じとする放射線照射条件とし、静止画撮影ではグリッド有りで撮影を行い、動態撮影ではグリッド無しで撮影を行うことで、従来方式に比べて各フレーム画像のＳ／Ｎ比を向上させて解析精度の向上を図ることとしてもよい。
【０１０６】
動態撮影用の放射線照射条件及び画像読取条件が設定されると、グリッド装着検知ＭＳ４５３からの出力に基づいて、グリッド４６がグリッド装着部４５２に装着されていない状態であるか否かが判断される（ステップＳ３４）。
【０１０７】
グリッド４６がグリッド装着部４５２に装着されていない状態であると判断されると（ステップＳ３４；ＹＥＳ）、操作部５３の操作による放射線照射の指示が待機される（ステップＳ３５）。ここで、撮影実施者は、安静呼吸の動態を撮影するために被験者（被写体Ｍ）に楽にするように指示し、安静呼吸を促す。撮影準備が整った時点で、操作部５３を操作して放射線照射指示を入力する。
【０１０８】
操作部５３により放射線照射指示が入力されると（ステップＳ３５；ＹＥＳ）、放射線照射制御装置４２及び読取制御装置４４に撮影開始指示が出力され、動態撮影が開始される（ステップＳ３６）。即ち、放射線照射制御装置４２に設定されたパルス間隔で放射線源４１により放射線が照射され、放射線検出部４３によりフレーム画像が取得される。予め定められたフレーム数の撮影が終了すると、制御部５１により放射線照射制御装置４２及び読取制御装置４４に撮影終了の指示が出力され、撮影動作が停止される。撮影されるフレーム数は、少なくとも１呼吸サイクルが撮影できる枚数である。
【０１０９】
一方、静止画撮影用の放射線照射条件及び画像読取条件が設定されると、グリッド装着検知ＭＳ４５３と被写体検知センサ４５４からの出力に基づいて、グリッド４６がグリッド装着部４５２に装着された状態であるか否かが判断される（ステップＳ３８）。
制御部５１は、このステップＳ３８によって、グリッド４６が装着されずに静止画撮影が行われることのないように制御している。
【０１１０】
グリッド４６がグリッド装着部４５２に装着されていると判断されると（ステップＳ３８；ＹＥＳ）、操作部５３の操作による放射線照射の指示が待機される（ステップＳ３９）。ここで、撮影実施者は、被験者に対し、息を吸った後に止めるように指示する。撮影準備が整った時点で、操作部５３を操作して放射線照射指示を入力する。
【０１１１】
操作部５３により放射線照射指示が入力されると（ステップＳ３９；ＹＥＳ）、放射線照射制御装置４２及び読取制御装置４４に撮影開始指示が出力され、静止画撮影が行われる（ステップＳ４０）。
【０１１２】
動態撮影又は静止画撮影が終了すると、撮影により取得された画像（各フレーム画像又は静止画像）は順次撮影用コンソール５に入力され、補正処理が行われる（ステップＳ４１）。ステップＳ４１の補正処理においては、オフセット補正処理、ゲイン補正処理、欠陥画素補正処理の３つの補正処理が必要に応じて行われる。まず最初に、取得された画像に対してオフセット補正処理が行われ、取得された画像に重畳された暗電流に起因するオフセット値が除去される。オフセット補正処理では、例えば、取得した画像の各画素値（濃度値。以下、信号値という）から、予め記憶されたオフセット補正係数を減算する処理が行われる。ここで、オフセット補正係数は、予め放射線非照射時に取得した複数の画像を平均化した画像である。次いで、ゲイン補正処理が行われ、取得された画像の各画素に対応する各検出素子の個体差や読み出しアンプのゲインムラによって生じる画素毎のばらつきが除去される。ゲイン補正処理では、例えば、オフセット補正後の画像に、予め記憶されたゲイン補正係数を乗算する処理が行われる。ここで、ゲイン補正係数は、放射線検出部４３に一様に放射線を照射した時に取得した複数のオフセット補正済み画像を平均化した画像と、このときの放射線照射条件で期待される出力信号値の関係から、補正後の各画素の信号値が一様となるように予め算出され、記憶された係数である。次いで、欠陥画素補正処理が行われ、周囲の画素と比較して感度が非線形な画素や、感度がない欠落画素が除去される。欠陥画素補正処理では、例えば、予め記憶された欠陥画素位置情報マップに従って、欠陥画素位置情報マップに登録された各欠陥画素において、欠陥画素の信号値をその近傍の欠陥でない画素の信号値の平均値で置き換える処理が行われる。ここで、欠陥画素位置情報マップは、放射線検出部４３に一様に放射線を照射した時に取得したオフセット補正、ゲイン補正済みのフレーム画像から、予め複数の欠陥画素が認識され、その欠陥画素の位置が登録されたマップである。上記オフセット補正係数及びゲイン補正係数、欠陥画素位置情報マップは、ビニングやダイナミックレンジ等の収集モードに応じて、予め、それぞれ最適な値が記憶されており、それぞれの収集モードにおいて対応する最適な値が読み出されるようになっている。なお、これらの補正処理は処理時間の短縮化の為、省略することも可能である。また、必要に応じて、取得された画像に対し、例えば間引き処理や階調処理等を施すことが好ましい。
【０１１３】
次いで、補正処理後の画像が記憶部５２に記憶されるとともに（ステップＳ４２）、表示部５４に表示される（ステップＳ４３）。動態撮影を行った場合は、各フレーム画像が撮影順を示す番号と対応付け記憶部５２に記憶される。ここで、取得された画像を記憶する直前に、各画素の信号値を真数から対数に変換する対数変換処理を行ってから記憶しても良い。撮影実施者は、表示された動態画像によりポジショニング等を確認し、撮影により診断に適した画像が取得された（撮影ＯＫ）か、再撮影が必要（撮影ＮＧ）か、を判断する。そして、操作部５３を操作して、判断結果を入力する。尚、撮影により取得された各フレーム画像は、全撮影の終了後に纏めて入力するようにしても良い。
【０１１４】
操作部５３の所定の操作により撮影ＯＫを示す判断結果が入力されると（ステップＳ４４；ＹＥＳ）、静止画撮影で取得された静止画像又は動態撮影で取得された一連のフレーム画像のそれぞれに、画像を識別するための識別ＩＤや、患者情報、検査対象部位、放射線照射条件、画像読取条件、撮影順を示す番号、撮影日時、撮影時のグリッド有無を示す情報（グリッド有無情報）等の情報が付帯され（例えば、ＤＩＣＯＭ形式で画像データのヘッダ領域に書き込まれ）、通信部５５を介して診断用コンソール６に送信される（ステップＳ４５）。そして、本処理は終了する。一方、操作部５３の所定の操作により撮影ＮＧを示す判断結果が入力されると（ステップＳ４４；ＮＯ）、記憶部５２に記憶された一連のフレーム画像が削除され（ステップＳ４６）、本処理は終了する。尚、このケースに於いては、再撮影が実行されることとなる。
【０１１５】
（診断用コンソール６の動作）
次に、診断用コンソール６における動作について説明する。
診断用コンソール６においては、通信部６５を介して撮影用コンソール５から静止画像が受信され、操作部６３により当該画像の表示指示が入力されると、表示部６４に受信された静止画像が表示され、医師の診断に供される。
一方、通信部６５を介して撮影用コンソール５から動態画像の一連のフレーム画像が受信され、操作部６３により動態解析が指示されると、制御部６１と記憶部６２に記憶されている画像解析処理プログラムとの協働により画像解析処理が実行される。
第２の実施の形態における画像解析処理は、第１の実施の形態で図３で示した画像解析処理と大まかな処理ステップは同じであるが詳細が異なるので、以下に図３及び図１３以降の図を参照して説明する。
【０１１６】
まず、各フレーム画像から肺野領域が抽出される（ステップＳ５１）。次いで、各フレーム画像の肺野領域が複数の小領域に分割され、各フレーム画像の小領域が互いに対応付けられる（ステップＳ５２）。ステップＳ５１〜Ｓ５２の具体的な処理については、第１の実施の形態において図３のステップＳ２１〜Ｓ２２で説明したものと同様であるので説明を援用する。
【０１１７】
次いで、一連のフレーム画像のステップＳ５２で対応付けられた各小領域毎に解析が行われ、吸気の特徴量及び呼気の特徴量がそれぞれ算出される（ステップＳ５３）。ここでは、呼気の特徴量及び吸気の特徴量として、例えば、呼気期及び吸気期のそれぞれにおける各小領域毎のフレーム間差分値（微分値）の代表値（絶対値の最大値）を算出する。フレーム間差分値は、そのフレーム画像が撮影されたタイミングにおける信号変化量を示す値である。呼吸により息を吸ったり吐いたりすれば、その息の流れに応じて肺の密度が変化し、これによってＸ線透過量（つまり、画素の出力信号値）が変化する。よって、信号変化量は、そのタイミングにおける気流速度を示す値とみなすことができる。なお、代表値としては、絶対値の最大値に限らず、中央値、平均値、最頻値としてもよい。
【０１１８】
具体的には、まず、各フレーム画像の各小領域内の画素の信号値（平均信号値）が算出される。
次いで、各フレーム画像間で対応づけられた各小領域に対して、時間軸方向のフィルタ処理が施される。このフィルタ処理は、血流等の高周波数な信号変化を除去し、換気による信号値の時間変化を抽出するための処理であり、例えば、小領域毎の信号値の時間変化に対して、安静呼吸画像群ではカットオフ周波数０．７Ｈｚ、深呼吸画像群ではカットオフ周波数０．５Ｈｚでローパスフィルタリングを行う。ここで、ローパスフィルタのカットオフ周波数は、固定値とするよりも、撮影された動態画像毎に最適化することが好ましい。例えば、前述のように一連のフレーム画像の横隔膜の位置を解析し、安静換気の場合は安静呼気位および安静吸気位となるフレームを検出し、安静呼気位のフレームとその次の安静吸気位のフレームの間のフレーム数から吸気期の時間を求め、その逆数に所定の係数を乗じた値をカットオフ周波数としたローパスフィルタを実施する。このとき安静換気の場合、自動的に設定されるカットオフ周波数は、０．２〜１．０Ｈｚの間に制限することが好ましい。また、ステップＳ３１において、別途測定された安静時における１分間の呼吸数、脈拍数等のバイタルサインを患者情報として入力しておき、これらの値からカットオフ周波数が算出されるようにしてもよい。例えば、患者情報として入力された１分間の呼吸数を１秒間の呼吸数に変換し、その呼吸数に所定の係数を乗じた値をカットオフ周波数としてローパスフィルタを実施してもよい。さらには、入力された１分間の脈拍数を１秒間の脈拍数に変換し、１秒間の呼吸数と１秒間の心拍数の平均値をカットオフ周波数としてローパスフィルタを実施するようにしてもよい。
フィルタ処理に次いで、撮影順が隣接するフレーム画像間で各小領域の信号値の差分を算出するフレーム間差分処理が行われる。ここでは、各小領域毎に、フレーム番号ＮとＮ＋１（Ｎは１、２、３・・・）のフレーム画像について、Ｎ＋１−Ｎの差分値が算出される。そして、呼気期のフレーム間差分値の最大値（絶対値の最大値）が呼気の特徴量として取得され、吸気期のフレーム間差分値の最大値（絶対値の最大値）が吸気の特徴量として取得される。フレーム間差分値の最大値（絶対値の最大値）は、最大微分値に相当する。ここで、各小領域内のフレーム間差分値の符号が正の期間が吸気期、負の期間が呼気期である。
【０１１９】
次いで、各小領域毎の吸気の特徴量と呼気の特徴量の比の値（吸気の特徴量／呼気の特徴量）が算出される（ステップＳ５４）。ここでは、「吸気期のフレーム間差分値の最大値／呼気期のフレーム間差分値の最大値」（最大流速比と呼ぶ）が算出される。
【０１２０】
次いで、算出された各小領域毎の「吸気の特徴量／呼気の特徴量」の値のヒストグラムが作成されるとともに、肺野全体における「吸気の特徴量／呼気の特徴量」の傾向を示す指標値（ここでは、平均値、標準偏差）が算出される（ステップＳ５５）。ヒストグラムの縦軸のカウント数は、肺野内の全小領域の数で除算することにより正規化することが好ましい。
【０１２１】
次いで、各小領域について求められた「吸気の特徴量／呼気の特徴量」の値を、予め記憶部６２に記憶されている、「吸気の特徴量／呼気の特徴量」の値と表示時のパラメータ値との変換テーブルに基づいて表示用パラメータ値に変換し、変換したパラメータ値にて、基準画像（例えば、安静呼気位のフレーム画像）の各小領域を表示した画像が作成される（ステップＳ５６）。変換テーブルは、例えば、正常／異常（重症度１〜ｎ）の各カテゴリに特徴量を分類するときの各カテゴリの特徴量の大きさの範囲を規定する閾値（各カテゴリの閾値）と色相、明度、輝度、透明度の何れかを一対一に対応付けたテーブルである。ここで、表示時のパラメータ値の変換テーブルとしては、特徴量の大きさに対する識別能を高めるために、各カテゴリの閾値と色相を対応付けることが好ましい。
このとき、例えば、数個（例えば、５〜６個）の色相を上記各カテゴリの閾値と対応させ、その間の特徴量の値に対しては中間の色相を割り当てる（グラデーションさせる）ことで識別能の高い表示とすることができる。
表示用パラメータ値に基づいて色付けされた画像は、基準画像のフレーム画像にオーバーレイ表示してもよい。
【０１２２】
なお、記憶部６２には、グリッド有りで撮影された動態画像に対応した変換テーブルと、グリッド無しで撮影された動態画像に対応する変換テーブルが記憶されており、ステップＳ５６及び次のステップＳ５７においては、一連のフレーム画像に付帯されているグリッド有無情報に基づいて、グリッド有りで撮影されたかグリッド無しで撮影されたかが判断され、その判断結果に応じた変換テーブルを用いて色付けが行われる。
【０１２３】
そして、表示部６４に、作成されたヒストグラム、及び作成された静止画像等が並列表示され（ステップＳ５７）、画像解析処理は終了する。ヒストグラムの領域は、上述の「吸気の特徴量／呼気の特徴量」の値と表示時のパラメータ値との変換テーブルに基づいて、基準画像の肺野領域の各小領域と同様の基準により色付けして表示される。
【０１２４】
ここで、グリッド有りで撮影したときとグリッド無しで撮影したときの動態解析への影響について説明する。
図１３は、或る正常の肺野をグリッド有りとグリッド無しで撮影した動態画像の動態解析結果の比較を示す図である。図１３においては、グリッド有りとグリッド無しで撮影した動態画像の各小領域を「吸気期のフレーム間差分値の最大値」（最大吸気気流速度）に基づいて色付けした画像、「呼気期のフレーム間差分値の最大値」（最大呼気気流速度）に基づいて色付けした画像、最大流速比に基づいて色付けした画像、及び最大流速比のヒストグラム、を解析結果として示している。
【０１２５】
図１３は、以下の撮影条件で撮影を行った動態画像の解析結果である。
検出器サイズ４０×３０ｃｍ、検出器画素サイズ１９４μｍ、グリッドピッチ８０ｌｉｎｅ／ｃｍ、グリッド比１２：１、管球〜検出器間距離２ｍ、フレーム数７５フレーム（約１０秒間撮影）、トータル被曝線量（検出器への到達線量を一定とした場合の被写体被曝線量）グリッド有り０．２４ｍＧｙ、グリッド無し０．１４ｍＧｙ
また、最大吸気気流速度、最大呼気気流速度、最大流速比の大きさと色（図１３においては濃度で示す）との変換テーブルは、両者を比較するために同一（ここでは、グリッド有り用の変換テーブル）を用いている。
【０１２６】
同一被写体を撮影した場合のグリッド有りの動態画像とグリッド無しの動態画像とでは、最大吸気気流速度、最大呼気気流速度、最大流速比はほぼ同等であるが、図１３に示すように、撮影系の特性等により若干相違が生じることがある。例えば図１３においてはグリッド無しで撮影された動態画像のヒストグラムの形状はグリッド有りのものと比較して幅広な形状となっている。そのため、例えば最大流速比によって肺野内の領域やヒストグラムの領域を正常、異常１〜ｎのカテゴリに分類して色付けする場合、グリッド有りとグリッド無しで同じ閾値（変換テーブル）を用いると、同じ最大流速比であっても異なる色表示、即ち異なる重症度への分類がなされる場合があり、好ましくない。そこで、図１３に示すように、グリッド有りと無しとで診断に影響する差異が生じる場合においては、グリッド有無に応じて特徴量の分類に用いる閾値（変換テーブル）を変更する必要がある。
【０１２７】
グリッド有りで撮影した動態画像を解析した結果とグリッド無しで撮影した動態画像を解析した結果とでどの程度の差異が生じるかは、撮影系の特性や解析内容等によって異なる。よって、撮影系によって、同じ被写体をグリッド有り無しで撮影した複数の画像を解析し、その結果を用いて帰納法的にグリッド有りの動態画像に使用する閾値とグリッド無しの動態画像に使用する閾値を算出しておくことが好ましい。
【０１２８】
なお、本実施の形態においては、撮影装置４は動態画像をグリッド無しで撮影するよう制御しているため、グリッド無しで撮影された動態画像用の閾値が記憶されていれば十分とも考えられる。しかし、グリッド有りで動態画像を撮影可能な撮影系に診断用コンソール６が接続される場合も想定される。この場合、グリッド有無の撮影条件が異なっていると、判断を誤ることになる。そこで、本実施の形態においては、動態画像を構成する各フレーム画像の付帯情報にグリッド有無情報を付帯させておき、制御部６１は、このグリッド有無情報に基づいて、撮影時のグリッド有無に応じた閾値を用いた解析アルゴリズムで解析を行っている。
また、同一患者の過去の画像の解析結果と現在の解析結果とを比較するような場合も想定される。互いにグリッド有無情報が付帯されている場合、一方の解析結果を他方の解析結果に合わせるべく、何れか一方の解析結果の値を補正して使用すること（例えば、両者の変換テーブルの閾値の差分を解析結果としての特徴量に加算又は減算する等）が好ましい。  
【０１２９】
図１４Ａ〜図１４Ｃに、ステップＳ５７において表示部６４に表示される表示画面の例を示す。
図１４Ａは、グリッド有りで撮影された正常者の肺野の動態画像を解析した解析結果を表示した表示画面である。図１４Ｂは、ＣＯＰＤ（閉塞性疾患）の肺野の動態画像を解析した解析結果を表示した表示画面である。図１４Ｃは、混合性疾患の肺野の動態画像を解析した解析結果を表示した表示画面である。
【０１３０】
図１４Ａ〜図１４Ｃに示すように、ステップＳ５７においては、半切サイズのフレーム画像から抽出された肺野領域内における各小領域（２ｍｍ角の矩形サイズ）の「吸気の特徴量／呼気の特徴量」の値のヒストグラム６４ａと、各小領域の「吸気の特徴量／呼気の特徴量」を一覧表示した静止画像６４ｂと、ヒストグラム６４ａ及び静止画像６４ｂで表示されている色相と「吸気の特徴量／呼気の特徴量」の値との関係を示す表示６４ｃと、肺野全体での「吸気の特徴量／呼気の特徴量」の傾向を示す指標値６４ｄが表示されている。また、図１４Ａ〜図１４Ｃに示すように、ヒストグラム６４ａは、「吸気の特徴量／呼気の特徴量」の値の大きさに応じて、横軸の領域が６つの色相で色分けして表示される。これにより、医師は、ヒストグラムを一瞥するだけで容易に肺野内における「吸気の特徴量／呼気の特徴量」の分布を把握することができる。また、各小領域の「吸気の特徴量／呼気の特徴量」を示す静止画像６４ｂにおいて、各小領域が「吸気の特徴量／呼気の特徴量」の値に応じてヒストグラムの色分けと同じ基準で色分けされて表示されるので、医師が肺野内の局所的な異常箇所（閉塞性部分、拘束性部分）を容易に把握することが可能となる。更に、肺野全体での「吸気の特徴量／呼気の特徴量」の傾向を示す指標値６４ｄとして、その平均値及び標準偏差を算出して画面上に併せて表示することで、肺野全体での「吸気の特徴量／呼気の特徴量」の傾向を数値として医師に提供することができる。
尚、本表示例に於いては、縦軸は、２ｍｍ角の小ブロックの出現個数／全ブロック数（（４３０ｍｍ÷２ｍｍ）×（３５４ｍｍ÷２ｍｍ）＝３８０５５個）の割合としているが、図８に示すように、縦軸を２ｍｍ角の小ブロックの出現個数自体としても良い。動態撮影に使用する検出器のサイズが複数存在する場合には、本表示例のような全体の領域個数で正規化した割合表示によって同一症例に対する解析結果表示の統一性が達成され、解析結果に対する医師等の誤認を招き難いので好ましい。
【０１３１】
ここで、呼気期のフレーム間差分値の最大値（絶対値）を呼気の特徴量、吸気期のフレーム間差分値の最大値（絶対値）を吸気の特徴量とした場合、正常者の肺野では、グリッド有りで動態撮影した動態画像で解析を行った場合、肺野全体の「吸気の特徴量／呼気の特徴量」の平均値は０．９〜１．２、標準偏差は０．１０〜０．２２程度であることが知られている。よって、ステップＳ５７で図１４Ａに示す表示画面が表示された場合、医師は、撮影した肺野が正常であることを容易に把握することができる。
一方、ＣＯＰＤ（閉塞性疾患）の肺野では、グリッド有りで動態撮影した動態画像で解析を行った場合、肺野全体の「吸気の特徴量／呼気の特徴量」の平均値は０．９〜１．２に収まらず（正常者に比べて大きく）、標準偏差も、０．１０〜０．２２に収まらない（正常者に比べて大きい）ことが知られている。よって、ステップＳ５７で図１４Ｂに示す表示画面が表示された場合、医師は、撮影した肺野がＣＯＰＤであることを容易に把握することができる。
一方、混合性肺疾患の肺野では、グリッド有りで動態撮影した動態画像で解析を行った場合、肺野全体の「吸気の特徴量／呼気の特徴量」の値が０．６６以下のデータ数、及び１．５以上のデータ数がともに増加することが知られている。よって、ステップＳ５７で図１４Ｃに示す表示画面が表示された場合、医師は、撮影した肺野が混合性疾患であることを容易に把握することができる。
このように、放射線撮影システム２００においては、「吸気の特徴量／呼気の特徴量」の傾向を示す指標値から、ＣＯＰＤ（閉塞性肺疾患）、混合性肺疾患等の換気不均等の病態もしくはその重症度を特定できるような、有用な診断情報を医師に提供することが可能となる。
【０１３２】
なお、「吸気の特徴量／呼気の特徴量」の傾向を示す指標値等を利用して被写体Ｍの正常／異常を判断することとしてもよい。この場合においても、グリッド有無情報に基づいて正常／異常を判断するための閾値を変更することが好ましい。例えば、上述の肺野全体の「吸気の特徴量／呼気の特徴量」の平均値に基づいて正常／異常と判断する場合、グリッド有りでは平均値０．９〜１．２を正常とし、グリッド無しでは平均値０．８〜１．３を正常として判断することが好ましい。
【０１３３】
呼気の特徴量及び吸気の特徴量としては、上記の例以外の他の特徴量を用いても良い。
例えば、呼気の特徴量を呼吸１サイクルにおける呼気期に相当するフレーム画像数（呼気時間）、吸気の特徴量を呼吸１サイクルにおける吸気期に相当するフレーム画像数（吸気期間）としてもよい。ここで、肺の換気機能が正常である場合、吸気時間と呼気時間はほぼ同じ長さ、もしくは若干呼気時間が長くなる筈である。よって、「呼気期に相当するフレーム画像数／吸気期に相当するフレーム画像数」の値を見れば、医師は肺疾患の疑いがあるか否かを把握することができる。特に、「呼気期のフレーム画像数／吸気期のフレーム画像数」＞１．５の領域は、呼気で換気がしにくく、取り込んだ空気を排出するのが遅れてしまう閉塞性部分であることがわかる。なお、「吸気期のフレーム間差分値の最大値／呼気期のフレーム間差分値の最大値≒呼気時間(呼気のフレーム画像数)／吸気時間(吸気のフレーム画像数)」の関係が成り立つため、医師は、呼気の特徴量を呼気期のフレーム間差分値の最大値、吸気の特徴量を吸気期のフレーム間差分値の最大値とした場合と同様の判断基準により、正常、ＣＯＰＤ(閉塞性肺疾患)、混合性肺疾患の識別を行うことが可能である。
【０１３４】
また、１呼吸サイクル分の各フレーム画像において、各小領域の画素の信号値（平均信号値）を算出し、各小領域毎に呼吸１サイクルにおける信号値の最小値及び最大値を求め、求めた最小値をその領域の呼気の特徴量、最大値を吸気の特徴量としてもよい。信号値の最大値と最小値は正常な箇所では両者の値の差が大きくなり、異常のある箇所では両者の差が非常に小さくなると考えられる。異常のある箇所では肺胞の動きが悪くなるため、肺胞の密度変化が小さくなるからである。よって、医師は、「信号値の最大値／信号値の最小値」のヒストグラムを参照し、平均値及び標準偏差を確認することで、肺野が正常であるか疾患であるかの判断材料とすることができる。例えば、肺野全体の「信号値の最大値／信号値の最小値」の平均値が１より大きく、標準偏差が小さい場合、肺の機能は正常であると判断することができる。一方、肺野全体の「信号値の最大値／信号値の最小値」の平均値が１に近い値であり、標準偏差が大きい場合、肺の機能に疾患があると判断することができる。
【０１３５】
その他、「吸気の特徴量／呼気の特徴量」の傾向を示す指標値として、平均値、標準偏差の他に、カウント数（ブロック（小領域）数）がヒストグラムのピークとなる「吸気の特徴量／呼気の特徴量」の値、もしくは、ピークのカウント数（ブロック数）、もしくは、「吸気の特徴量／呼気の特徴量」の値が所定以上所定以下となるカウント数の割合を指標値として用いても良い。もしくは、これらの複数の指標値を組み合わせて新たな指標値を作成しても良い。
【０１３６】
以上説明したように、放射線撮影システム２００によれば、静止画撮影時にはグリッド４６を使用し、動態撮影時にはグリッド４６を使用せずに撮影を行う。これにより、グリッドによる線量の減弱がなくなる分、同じ被写体への照射線量に対して放射線検出部４３への到達する線量が増えるので、高精細な動画解析を行うこと可能となる。また、放射線検出部４３への到達線量を一定とすると、従来のように動態撮影時にグリッドを使用した場合に比べて被写体への照射線量（被写体の被曝線量）の低減を図ることが可能となる。
【０１３７】
また、動態撮影により得られた各フレーム画像には、グリッド有無情報が付帯され、診断用コンソール６の制御部６１は、グリッド有無情報に基づいて解析アルゴリズムを変更する、例えば、解析に用いる閾値を変更するので、動態撮影時のグリッド有無に拘わらず正確な解析結果を提供することが可能となる。
【０１３８】
なお、上述した本実施の形態における記述は、本発明に係る好適な一例であり、これに限定されるものではない。
【０１３９】
例えば、グリッド無しで撮影を行う場合、散乱線影響を緩和して、より良い動態観察や動態解析を行うには、被写体Ｍを放射線検出部４３と離間させて撮影を行うことが好ましい。被写体Ｍを放射線検出部４３と離間させて撮影を行うと、エアギャップにより散乱線の影響を低減することができる。離間により拡大撮影となるが、動態解析結果には影響はない。特許3,861,572号公報に記載のように、位相コントラスト撮影の撮影条件としてもよい。なお、この際には、被写体が離間しているか密着しているかを被写体検知センサ４５４で検知し、グリッド無で密着の場合には動態撮影を禁止するように制御しても良い。
【０１４０】
また、例えば、上記実施の形態においては、「吸気の特徴量／呼気の特徴量」の解析結果として、最大流速比をヒストグラムにより表示することとしたが、これに限定されない。
【０１４１】
例えば、図１５に示すように、「吸気の特徴量／呼気の特徴量」の傾向を示す指標値の一つをＸ軸、もう一つをＹ軸としたグラフに、各々の指標値に対する正常、異常の閾値ＴＨ１を示し、そのグラフ上に、動態画像から算出した肺野全体での「吸気の特徴量／呼気の特徴量」の傾向を示す指標値をプロットした図を解析結果としてもよい。図１５は、Ｘ軸を「吸気期のフレーム間差分値の最大値／呼気期のフレーム間差分値の最大値」の平均値、Ｙ軸をその標準偏差とし、動態画像から算出した肺野全体における「吸気の特徴量／呼気の特徴量」の平均値、標準偏差の指標値をプロットしたグラフである。このようなグラフを用いて「吸気の特徴量／呼気の特徴量」の傾向を示す指標値を表示することにより、プロットされた点から閾値ＴＨ１までの距離によって異常の程度が視覚的に容易に把握することが可能となる。
【０１４２】
また、例えば、２つの指標値（例えば、「吸気の特徴量／呼気の特徴量」の平均値、標準偏差）を線形結合した値を新たな指標値とし、図１６に示すように２つの指標値の一方をＸ軸、他方をＹ軸としたグラフに、新たな指標値（２つの指標値を線形結合した指標値）を重症度により分類するための閾値ｔｈ１〜ｔｈ４を示し、このグラフ上に動態画像から算出した新たな指標値をプロットしてもよい。線形結合の例としては、平均値及び標準偏差に対する多数の測定値データから主成分分析により算出された第１主成分を指標値とすることができる。このようなグラフを用いることにより、異常の程度が視覚的により容易に把握することが可能となる。また、例えば、複数（Ｍ個）の指標値に対する多数（Ｎ個）の測定値データより、Ｎ×Ｍ個のデータに対する共分散行列の最大固有値を算出することで第一主成分を算出し、算出された主成分を指標値として用いてもよい。
【０１４３】
また、「吸気の特徴量／呼気の特徴量」の傾向を示す指標値の他に、吸気の特徴量、もしくは呼気の特徴量のそれぞれに対して傾向を示す指標値を算出してもよい。例えば、図１７に示すように、小領域毎に算出した吸気もしくは呼気の特徴量に対して、左右肺野内をそれぞれ上中下の３つに分割した計６つの領域の各々において変動係数（＝標準偏差／平均値）を算出し、変動係数の大きさに応じた色相（又は輝度又は彩度）によって、６つの領域を色付けし表示した静止画像を作成してもよい。このような表示をすることで、不均等換気の分布を把握することが容易となり、不均等換気が行われている部分が区域的であるかもしくはびまん性であるかが容易に判定可能となる。また、「吸気の特徴量／呼気の特徴量」を対数に変換し、対数変換後の値に対して、上記ヒストグラムの作成、指標値の算出、表示時パラメータとの変換テーブルに基づく色付け画像の作成を行っても良い。
【０１４４】
また、上記実施の形態においては、動態撮影と静止画撮影の双方の撮影が可能な撮影装置を含むシステムについて説明したが、動態撮影のみが可能な撮影装置においても、グリッドを使用せずに撮影を行うことは、被写体の被曝量低減又は画質の向上に有効である。
【０１４５】
また、例えば、上記の説明では、本発明に係るプログラムのコンピュータ読み取り可能な媒体としてハードディスクや半導体の不揮発性メモリ等を使用した例を開示したが、この例に限定されない。その他のコンピュータ読み取り可能な媒体として、ＣＤ-ＲＯＭ等の可搬型記録媒体を適用することが可能である。また、本発明に係るプログラムのデータを通信回線を介して提供する媒体として、キャリアウエーブ(搬送波)も適用される。
【０１４６】
その他、胸部診断支援システム１００、放射線撮影システム２００を構成する各装置の細部構成及び細部動作に関しても、本発明の趣旨を逸脱することのない範囲で適宜変更可能である。
【０１４７】
なお、明細書、請求の範囲、図面及び要約を含む２０１０年８月２７日に出願された日本特許出願Ｎｏ．２０１０−１９０２４１号、及び２０１１年１月２６日に出願された日本特許出願Ｎｏ．２０１１−０１３５３８号の全ての開示は、そのまま本出願の一部に組み込まれる。
【産業上の利用可能性】
【０１４８】
医療の分野における胸部動態等の診断の支援に利用可能性がある。
【符号の説明】
【０１４９】
１００  胸部診断支援システム
２００  放射線撮影システム
１、４  撮影装置
１１、４１  放射線源
１２、４２  放射線照射制御装置
１３、４３  放射線検出部
１４、４４  読取制御装置
２、５  撮影用コンソール
２１、５１  制御部
２２、５２  記憶部
２３、５３  操作部
２４、５４  表示部
２５、５５  通信部
２６、５６  バス
３、６  診断用コンソール
３１、６１  制御部
３２、６２  記憶部
３３、６３  操作部
３４、６４  表示部
３５、６５  通信部
３６、６６  バス
【図１】
2012026146.tif 000003
【図２】
2012026146.tif 000004
【図３】
2012026146.tif 000005
【図４】
2012026146.tif 000006
【図５】
2012026146.tif 000007
【図６】
2012026146.tif 000008
【図７】
2012026146.tif 000009
【図８Ａ】
2012026146.tif 000010
【図８Ｂ】
2012026146.tif 000011
【図８Ｃ】
2012026146.tif 000012
【図９】
2012026146.tif 000013
【図１０】
2012026146.tif 000014
【図１１】
2012026146.tif 000015
【図１２Ａ】
2012026146.tif 000016
【図１２Ｂ】
2012026146.tif 000017
【図１３】
2012026146.tif 000018
【図１４Ａ】
2012026146.tif 000019
【図１４Ｂ】
2012026146.tif 000020
【図１４Ｃ】
2012026146.tif 000021
【図１５】
2012026146.tif 000022
【図１６】
2012026146.tif 000023
【図１７】
2012026146.tif 000024
【国際調査報告】
2012026146.tif 000025
2012026146.tif 000026
2012026146.tif 000027
2012026146.tif 000028