"少なくとも1サイクルの呼吸周期を含む胸部の動態を撮影し、連続する複数のフレーム画像を生成する撮影手段と、前記撮影手段により生成された複数のフレーム画像のそれぞれから肺野領域を抽出する抽出手段と、前記抽出手段により抽出された肺野領域を複数の小領域に分割し、前記複数のフレーム画像間にわたり小領域を対応付ける領域分割手段と、前記複数のフレーム画像間において対応付けられた小領域毎に解析を行って、前記小領域毎の吸気の特徴量及び呼気の特徴量を算出し、前記算出された吸気の特徴量と呼気の特徴量との比の値を算出し、当該算出された比の値のヒストグラムを作成する解析手段と、前記解析手段により作成されたヒストグラムを表示する表示手段と、を備える胸部診断支援システム。前記吸気の特徴量が、吸気期に相当するフレーム画像群の時間的に隣接するフレーム画像間で画素信号値の差分を取ることにより得られるフレーム間差分値の絶対値の最大値であり、前記呼気の特徴量が、呼気期に相当するフレーム画像群の時間的に隣接するフレーム画像間で画素信号値の差分を取ることにより得られるフレーム間差分値の絶対値の最大値である請求項1に記載の胸部診断支援システム。前記吸気の特徴量が、呼吸1サイクルにおける吸気期に相当するフレーム画像の数であり、前記呼気の特徴量が、前記呼吸1サイクルにおける呼気期に相当するフレーム画像の数である請求項1に記載の胸部診断支援システム。前記吸気の特徴量が、呼吸1サイクル分の複数のフレーム画像における画素信号値の最大値であり、前記呼気の特徴量が、前記呼吸1サイクル分の複数のフレーム画像における画素信号値の最小値である請求項1に記載の胸部診断支援システム。前記領域分割手段は、前記複数のフレーム画像のうち一のフレーム画像の肺野領域を複数の小領域に分割し、他のフレーム画像の肺野領域を前記一のフレーム画像の各小領域と同じ画素位置の小領域に分割し、前記複数のフレーム画像間の前記同じ画素位置の各小領域を互いに対応付ける請求項1~4の何れか一項に記載の胸部診断支援システム。前記領域分割手段は、前記複数のフレーム画像のうち一のフレーム画像の肺野領域を複数の小領域に分割し、他のフレーム画像のそれぞれにおいて前記一のフレーム画像における各小領域とマッチング度合いの高い領域を抽出し、前記他のフレーム画像の肺野領域を前記抽出された各領域に分割し、前記複数のフレーム画像間の前記マッチング度合いの高い各小領域を互いに対応付ける請求項1~4の何れか一項に記載の胸部診断支援システム。前記撮影手段は、グリッドを使用せずに動態撮影を行う請求項1~6の何れか一項に記載の胸部診断支援システム。コンピュータを、少なくとも1サイクルの呼吸周期を含む胸部の動態を示す複数のフレーム画像のそれぞれから肺野領域を抽出する抽出手段、前記抽出手段により抽出された肺野領域を複数の小領域に分割し、前記複数のフレーム画像間にわたり小領域を対応付ける領域分割手段、前記複数のフレーム画像間において対応付けられた小領域毎に解析を行って、前記小領域毎の吸気の特徴量及び呼気の特徴量を算出し、前記算出された吸気の特徴量と呼気の特徴量との比の値を算出し、当該算出された比の値のヒストグラムを作成する解析手段、前記解析手段により作成されたヒストグラムを表示する表示手段、として機能させるためのプログラム。しかしながら、特許文献1に記載のように、フレーム間差分値の最大値平均値、中間値等を算出して表示するだけでは、熟練した呼吸器専門医以外の医師が病態症例名を把握することは難しい。本発明の課題は、肺野の換気に係る病態の把握を容易にし、診断に有効となる情報を医師に提供できるようにすることである。特に、聴診器による実務経験の少ない医師であっても、的確な診断を行える診断情報を提供することである。上記課題を解決するため、本発明の第1の側面によると、胸部診断支援システムは、少なくとも1サイクルの呼吸周期を含む胸部の動態を撮影し、連続する複数のフレーム画像を生成する撮影手段と、前記撮影手段により生成された複数のフレーム画像のそれぞれから肺野領域を抽出する抽出手段と、前記抽出手段により抽出された肺野領域を複数の小領域に分割し、前記複数のフレーム画像間にわたり小領域を対応付ける領域分割手段と、前記複数のフレーム画像間において対応付けられた小領域毎に解析を行って、前記小領域毎の吸気の特徴量及び呼気の特徴量を算出し、前記算出された吸気の特徴量と呼気の特徴量との比の値を算出し、当該算出された比の値のヒストグラムを作成する解析手段と、前記解析手段により作成されたヒストグラムを表示する表示手段と、を備える。前記吸気の特徴量が、吸気期に相当するフレーム画像群の時間的に隣接するフレーム画像間で画素信号値の差分を取ることにより得られるフレーム間差分値の絶対値の最大値であり、前記呼気の特徴量が、呼気期に相当するフレーム画像群の時間的に隣接するフレーム画像間で画素信号値の差分を取ることにより得られるフレーム間差分値の絶対値の最大値であることが好ましい。または、前記吸気の特徴量が、呼吸1サイクルにおける吸気期に相当するフレーム画像の数であり、前記呼気の特徴量が、前記呼吸1サイクルにおける呼気期に相当するフレーム画像の数であることが好ましい。または、前記吸気の特徴量が、呼吸1サイクル分の複数のフレーム画像における画素信号値の最大値であり、前記呼気の特徴量が、前記呼吸1サイクル分の複数のフレーム画像における画素信号値の最小値であることが好ましい。前記領域分割手段は、前記複数のフレーム画像のうち一のフレーム画像の肺野領域を複数の小領域に分割し、他のフレーム画像の肺野領域を前記一のフレーム画像の各小領域と同じ画素位置の小領域に分割し、前記複数のフレーム画像間の前記同じ画素位置の各小領域を互いに対応付けることが好ましい。前記領域分割手段は、前記複数のフレーム画像のうち一のフレーム画像の肺野領域を複数の小領域に分割し、他のフレーム画像のそれぞれにおいて前記一のフレーム画像における各小領域とマッチング度合いの高い領域を抽出し、前記他のフレーム画像の肺野領域を前記抽出された各領域に分割し、前記複数のフレーム画像間の前記マッチング度合いの高い各小領域を互いに対応付けることが好ましい。前記撮影手段は、グリッドを使用せずに動態撮影を行うことが好ましい。本発明の第2の側面によると、プログラムは、コンピュータを、少なくとも1サイクルの呼吸周期を含む胸部の動態を示す複数のフレーム画像のそれぞれから肺野領域を抽出する抽出手段、前記抽出手段により抽出された肺野領域を複数の小領域に分割し、前記複数のフレーム画像間にわたり小領域を対応付ける領域分割手段、前記複数のフレーム画像間において対応付けられた小領域毎に解析を行って、前記小領域毎の吸気の特徴量及び呼気の特徴量を算出し、前記算出された吸気の特徴量と呼気の特徴量との比の値を算出し、当該算出された比の値のヒストグラムを作成する解析手段と、前記解析手段により作成されたヒストグラムを表示する表示手段、として機能させる。本発明によれば、肺野の換気に係る解析結果病態情報を視認可能に提供するので、医師の経験の度合いに依存しない病態の把握が可能となり、病変部の早期発見に基づく早期治療が可能となる。<第1の実施の形態>以下、図面を参照して本発明の第1の実施の形態を詳細に説明する。ただし、発明の範囲は、図示例に限定されない。胸部診断支援システム100の構成まず、構成を説明する。図1に、本実施の形態における胸部診断支援システム100の全体構成を示す。図1に示すように、胸部診断支援システム100は、撮影装置1と、撮影用コンソール2とが通信ケーブル等により接続され、撮影用コンソール2と、診断用コンソール3とがLANLocalAreaNetwork等の通信ネットワークNTを介して接続されて構成されている。胸部診断支援システム100を構成する各装置は、DICOMDigitalImageandCommunicationsinMedicine規格に準じており、各装置間の通信は、DICOMに則って行われる。撮影装置1の構成撮影装置1は、例えば、呼吸運動に伴う肺の膨張及び収縮の形態変化、心臓の拍動等の、周期性サイクルを持つ胸部の動態を撮影する装置である。動態撮影は、人体の胸部に対し、X線等の放射線を連続照射して複数の画像を取得即ち、連続撮影することにより行う。この連続撮影により得られた一連の画像を動態画像と呼ぶ。また、動態画像を構成する複数の画像のそれぞれをフレーム画像と呼ぶ。撮影装置1は、図1に示すように、放射線源11、放射線照射制御装置12、放射線検出部13、読取制御装置14等を備えて構成されている。放射線源11は、被写体Mを挟んで放射線検出部13と対向する位置に配置され、放射線照射制御装置12の制御に従って、被写体Mに対し放射線X線を照射する。放射線照射制御装置12は、撮影用コンソール2に接続されており、撮影用コンソール2から入力された放射線照射条件に基づいて放射線源11を制御して放射線撮影を行う。撮影用コンソール2から入力される放射線照射条件は、例えば、連続照射時のパルスレート、パルス幅、パルス間隔、1撮影あたりの撮影フレーム数、X線管電流の値、X線管電圧の値、フィルタ種等である。パルスレートは、1秒あたりの放射線照射回数であり、後述するフレームレートと一致している。パルス幅は、放射線照射1回当たりの放射線照射時間である。パルス間隔は、連続撮影において、1回の放射線照射開始から次の放射線照射開始までの時間であり、後述するフレーム間隔と一致している。放射線検出部13は、FPD等の半導体イメージセンサにより構成される。FPDは、例えば、ガラス基板等を有しており、基板上の所定位置に、放射線源11から照射されて少なくとも被写体Mを透過した放射線をその強度に応じて検出し、検出した放射線を電気信号に変換して蓄積する複数の画素がマトリックス状に配列されている。各画素は、例えばTFTThinFilmTransistor等のスイッチング部により構成されている。読取制御装置14は、撮影用コンソール2に接続されている。読取制御装置14は、撮影用コンソール2から入力された画像読取条件に基づいて放射線検出部13の各画素のスイッチング部を制御して、当該各画素に蓄積された電気信号の読み取りをスイッチングしていき、放射線検出部13に蓄積された電気信号を読み取ることにより、画像データを取得する。この画像データがフレーム画像である。そして、読取制御装置14は、取得したフレーム画像を撮影用コンソール2に出力する。画像読取条件は、例えば、フレームレート、フレーム間隔、画素サイズ、画像サイズマトリックスサイズ等である。フレームレートは、1秒あたりに取得するフレーム画像数であり、パルスレートと一致している。フレーム間隔は、連続撮影において、1回のフレーム画像の取得動作開始から次のフレーム画像の取得動作開始までの時間であり、パルス間隔と一致している。ここで、放射線照射制御装置12と読取制御装置14は互いに接続され、互いに同期信号をやりとりして放射線照射動作と画像の読み取りの動作を同調させるようになっている。この他、後述するオフセット補正に用いるオフセット補正係数を算出するための複数の暗画像を取得するキャリブレーション時は、放射線照射動作と同期せず、放射線が照射されない状態で、リセット~蓄積~データ読取~リセットの一連の画像の読み取り動作を行うが、一連の動態撮影前、一連の動態撮影後のいずれかのタイミングで行うようにしてもよい。撮影用コンソール2の構成撮影用コンソール2は、放射線照射条件や画像読取条件を撮影装置1に出力して撮影装置1による放射線撮影及び放射線画像の読み取り動作を制御するとともに、撮影装置1により取得された動態画像を撮影技師によるポジショニングの確認や診断に適した画像であるか否かの確認用に表示する。撮影用コンソール2は、図1に示すように、制御部21、記憶部22、操作部23、表示部24、通信部25を備えて構成され、各部はバス26により接続されている。制御部21は、CPUCentralProcessingUnit、RAMRandomAccessMemory等により構成される。制御部21のCPUは、操作部23の操作に応じて、記憶部22に記憶されているシステムプログラムや各種処理プログラムを読み出してRAM内に展開し、展開されたプログラムに従って後述する撮影制御処理を始めとする各種処理を実行し、撮影用コンソール2各部の動作や、撮影装置1の放射線照射動作及び読み取り動作を集中制御する。記憶部22は、不揮発性の半導体メモリやハードディスク等により構成される。記憶部22は、制御部21で実行される各種プログラムやプログラムにより処理の実行に必要なパラメータ、或いは処理結果等のデータを記憶する。例えば、記憶部22は、図2に示す撮影制御処理を実行するための撮影制御処理プログラムを記憶している。また、記憶部22は、検査対象部位に対応付けて放射線照射条件及び画像読取条件を記憶している。各種プログラムは、読取可能なプログラムコードの形態で格納され、制御部21は、当該プログラムコードに従った動作を逐次実行する。操作部23は、カーソルキー、数字入力キー、及び各種機能キー等を備えたキーボードと、マウス等のポインティングデバイスを備えて構成され、キーボードに対するキー操作やマウス操作により入力された指示信号を制御部21に出力する。また、操作部23は、表示部24の表示画面にタッチパネルを備えても良く、この場合、タッチパネルを介して入力された指示信号を制御部21に出力する。表示部24は、LCDLiquidCrystalDisplayやCRTCathodeRayTube等のモニタにより構成され、制御部21から入力される表示信号の指示に従って、操作部23からの入力指示やデータ等を表示する。通信部25は、LANアダプタやモデムやTATerminalAdapter等を備え、通信ネットワークNTに接続された各装置との間のデータ送受信を制御する。診断用コンソール3の構成診断用コンソール3は、撮影用コンソール2から動態画像を取得し、取得した動態画像や後述するヒストグラム等の診断支援情報を表示して医師が読影診断するための動画像処理装置である。診断用コンソール3は、図1に示すように、制御部31、記憶部32、操作部33、表示部34、通信部35を備えて構成され、各部はバス36により接続されている。制御部31は、CPU、RAM等により構成される。制御部31のCPUは、操作部33の操作に応じて、記憶部32に記憶されているシステムプログラムや、各種処理プログラムを読み出してRAM内に展開し、展開されたプログラムに従って、後述する画像解析処理を始めとする各種処理を実行し、診断用コンソール3各部の動作を集中制御する。制御部31は、後述する画像解析処理を実行することにより、抽出手段、領域分割手段、解析手段としての機能を実現する。記憶部32は、不揮発性の半導体メモリやハードディスク等により構成される。記憶部32は、制御部31で画像解析処理を実行するための画像解析処理プログラムを始めとする各種プログラムやプログラムにより処理の実行に必要なパラメータ、或いは処理結果等のデータを記憶する。これらの各種プログラムは、読取可能なプログラムコードの形態で格納され、制御部31は、当該プログラムコードに従った動作を逐次実行する。操作部33は、カーソルキー、数字入力キー、及び各種機能キー等を備えたキーボードと、マウス等のポインティングデバイスを備えて構成され、キーボードに対するキー操作やマウス操作により入力された指示信号を制御部31に出力する。また、操作部33は、表示部34の表示画面にタッチパネルを備えても良く、この場合、タッチパネルを介して入力された指示信号を制御部31に出力する。表示部34は、LCDやCRT等のモニタにより構成され、制御部31から入力される表示信号の指示に従って、操作部33からの入力指示やデータ等を表示する。通信部35は、LANアダプタやモデムやTA等を備え、通信ネットワークNTに接続された各装置との間のデータ送受信を制御する。胸部診断支援システム100の動作次に、上記胸部診断支援システム100における動作について説明する。撮影装置1、撮影用コンソール2の動作まず、撮影装置1、撮影用コンソール2による撮影動作について説明する。図2に、撮影用コンソール2の制御部21において実行される撮影制御処理を示す。撮影制御処理は、制御部21と記憶部22に記憶されている撮影制御処理プログラムとの協働により実行される。まず、撮影技師により撮影用コンソール2の操作部23が操作され、撮影対象被写体Mの患者情報患者の氏名、身長、体重、年齢、性別等の入力が行われるステップS1。次いで、放射線照射条件が記憶部22から読み出されて放射線照射制御装置12に設定されるとともに、画像読取条件が記憶部22から読み出されて読取制御装置14に設定されるステップS2。次いで、操作部23の操作による放射線照射の指示が待機されるステップS3。ここで、撮影実施者は、撮影装置1において被写体Mのポジショニング等の撮影準備を行うとともに、安静呼吸の動態を撮影するために被験者被写体Mに楽にするように指示し、安静呼吸を促す。撮影準備が整った時点で、操作部23を操作して放射線照射指示を入力する。操作部23により放射線照射指示が入力されるとステップS3;YES、放射線照射制御装置12及び読取制御装置14に撮影開始指示が出力され、動態撮影が開始されるステップS4。即ち、放射線照射制御装置12に設定されたパルス間隔で放射線源11により放射線が照射され、放射線検出部13によりフレーム画像が取得される。予め定められたフレーム数の撮影が終了すると、制御部21により放射線照射制御装置12及び読取制御装置14に撮影終了の指示が出力され、撮影動作が停止される。撮影されるフレーム数は、少なくとも1呼吸サイクルが撮影できる枚数である。撮影により取得されたフレーム画像は順次撮影用コンソール2に入力され、各フレーム画像に対して補正処理が行われるステップS5。ステップS5の補正処理においては、オフセット補正処理、ゲイン補正処理、欠陥画素補正処理の3つの補正処理が行われる。まず最初に、取得された各フレーム画像に対してオフセット補正処理が行われ、取得された各フレーム画像に重畳された暗電流に起因するオフセット値が除去される。オフセット補正処理では、例えば、取得した各フレーム画像の各画素値濃度値。以下、信号値というから、予め記憶されたオフセット補正係数を減算する処理が行われる。ここで、オフセット補正係数は、予め放射線非照射時に取得した複数のフレーム画像を平均化した画像である。次いで、ゲイン補正処理が行われ、各フレーム画像の各画素に対応する各検出素子の個体差や読み出しアンプのゲインムラによって生じる画素毎のばらつきが除去される。ゲイン補正処理では、例えば、オフセット補正後の各フレーム画像に、予め記憶されたゲイン補正係数を乗算する処理が行われる。ここで、ゲイン補正係数は、放射線検出部13に一様に放射線を照射した時に取得した複数のオフセット補正済みフレーム画像を平均化した画像と、このときの放射線照射条件で期待される出力信号値の関係から、補正後の各画素の信号値が一様となるように予め算出され、記憶された係数である。次いで、欠陥画素補正処理が行われ、周囲の画素と比較して感度が非線形な画素や、感度がない欠落画素が除去される。欠陥画素補正処理では、例えば、予め記憶された欠陥画素位置情報マップに従って、欠陥画素位置情報マップに登録された各欠陥画素において、欠陥画素の信号値をその近傍の欠陥でない画素の信号値の平均値で置き換える処理が行われる。ここで、欠陥画素位置情報マップは、放射線検出部13に一様に放射線を照射した時に取得したオフセット補正、ゲイン補正済みのフレーム画像から、予め複数の欠陥画素が認識され、その欠陥画素の位置が登録されたマップである。上記オフセット補正係数及びゲイン補正係数、欠陥画素位置情報マップは、ビニングやダイナミックレンジ等の収集モードに応じて、予め、それぞれ最適な値が記憶されており、それぞれの収集モードにおいて対応する最適な値が読み出されるようになっている。次いで、補正処理後の各フレーム画像と、撮影順を示す番号と対応付けて記憶部22に記憶されるとともにステップS6、表示部24に表示されるステップS7。ここで、各フレーム画像を記憶する直前に、各フレーム画像の各画素の信号値を真数から対数に変換する対数変換処理を行ってから記憶しても良い。撮影技師は、表示された動態画像によりポジショニング等を確認し、撮影により診断に適した画像が取得された撮影OKか、再撮影が必要撮影NGか、を判断する。そして、操作部23を操作して、判断結果を入力する。尚、撮影により取得された各フレーム画像は、全撮影の終了後に纏めて入力するようにしても良い。操作部23の所定の操作により撮影OKを示す判断結果が入力されるとステップS8;YES、動態撮影で取得された一連のフレーム画像のそれぞれに、動態画像を識別するための識別IDや、患者情報、検査対象部位、放射線照射条件、画像読取条件、撮影順を示す番号、撮影日時等の情報が付帯され例えば、DICOM形式で画像データのヘッダ領域に書き込まれ、通信部25を介して診断用コンソール3に送信されるステップS9。そして、本処理は終了する。一方、操作部23の所定の操作により撮影NGを示す判断結果が入力されるとステップS8;NO、記憶部22に記憶された一連のフレーム画像が削除されステップS10、本処理は終了する。尚、このケースに於いては、再撮影が実行されることとなる。診断用コンソール3の動作次に、診断用コンソール3における動作について説明する。診断用コンソール3においては、通信部35を介して撮影用コンソール2から動態画像の一連のフレーム画像が受信されると、制御部31と記憶部32に記憶されている画像解析処理プログラムとの協働により図3に示す画像解析処理が実行される。以下、図3を参照して画像解析処理の流れについて説明する。まず、各フレーム画像から肺野領域が抽出されるステップS21。肺野領域の抽出方法は何れの方法であってもよい。例えば、一連のフレーム画像中の任意のフレーム画像ここでは、撮影順が一番最初のフレーム画像とする。の各画素の信号値濃度値のヒストグラムから判別分析によって閾値を求め、この閾値より高信号の領域を肺野領域候補として1次抽出する。次いで、1次抽出された肺野領域候補の境界付近でエッジ検出を行い、境界付近の小領域でエッジが最大となる点を境界に沿って抽出すれば肺野領域の境界を抽出することができる。次いで、各フレーム画像の肺野領域が複数の小領域に分割され、各フレーム画像の小領域が互いに対応付けられるステップS22。各小領域の画素の位置は、制御部31のRAMに記憶される。ここで、呼吸サイクルは、呼気期と吸気期により構成される。図4は、1つの呼吸サイクル深呼吸時において撮影された複数の時間位相TT=t0~t6のフレーム画像を示す図である。図4に示すように、呼気期は、横隔膜が上がることによって肺から空気が排出され、肺野の領域が小さくなる。最大呼気位では、横隔膜の位置が最も高い状態となる。吸気期は、横隔膜が下がることにより肺に空気が取り込まれ、図4に示すように胸郭中の肺野の領域が大きくなる。最大吸気位では、横隔膜の位置が最も下がった状態となる。即ち、肺野領域の同一部分の位置は呼吸運動に応じて時間とともに変化するため、各フレーム画像間で肺野の同一の部分特に、下部領域横隔膜付近を示す画素位置がずれることになる。しかし、安静呼吸時に撮影された画像においては、上記の位置ずれは小さく、後述する解析結果が狂ってしまうほどの位置ずれはおきない。図5の画像D1は、安静呼気位安静呼吸時において横隔膜の位置が最も高くなったタイミングのフレーム画像である。図5の画像D2は、安静吸気位安静呼吸時において横隔膜の位置が最も低くなったタイミングのフレーム画像である。即ち、図5の画像D1とD2とは、呼吸1サイクルにおいて最も形状の差の大きいタイミングで撮影された画像である。しかし、図5の画像D1、D2間においては、最も位置ずれの大きい肺野領域の下部領域においても位置ずれはわずかであることがわかる画像D2のA11は画像D1のA1と同じ画素位置を示し、画像D2のA2は画像D1のA1と肺野における同一部分を描画した領域を示している。そこで、ステップS22における具体的な処理としては、まず、一連のフレーム画像の中から一のフレーム画像を基準画像として設定する。次いで、基準画像の抽出された肺野領域を複数の小領域例えば、2mm×2mmの矩形領域に分割する図5参照。次いで、他のフレーム画像の肺野領域を基準画像の各小領域と同じ画素位置の小領域放射線検出部13の同じ撮像素子から出力される信号値の領域に分割する。次いで、各フレーム画像間の同じ画素位置の各小領域を互いに対応付ける。この処理では、フレーム画像の小領域への分割及び対応付けを高速に行うことが可能となる。基準画像としては、安静呼気位のフレーム画像とすることが好ましい。安静呼気位では、安静呼吸時において横隔膜の位置が最も高くなる、即ち、肺野領域の面積が最も小さくなるので、基準画像の小領域を他のフレーム画像に対応付けたときに、小領域が他のフレーム画像の肺野外の領域に対応付けられることがないからである。安静呼気位の画像は、一連のフレーム画像の中から横隔膜の位置が最も高い位置にある画像を抽出することで取得することができる。横隔膜の位置は、例えば、図6に示す横隔膜の基準位置B1を横隔膜のカーブC図6に点線で示すの垂直方向の平均位置と予め定義しておき、肺野領域Rから横隔膜のカーブC肺野領域の下端を抽出し、その垂直方向の平均位置を求め、求めた位置を横隔膜の基準位置B1として特定する。なお、撮影された動態画像が深呼吸時の画像である場合は、図7に示すように、肺野の同一の部分を示す画素位置が大幅にずれる。ここで、図7の画像D3は、深呼吸時の最大呼気位のフレーム画像であり、図7の画像D4は、深呼吸時の最大吸気位のフレーム画像である画像D4のA31は画像D3のA3と同じ画素位置を示し、画像D4のA4は画像D3のA3と肺野における同一部分を描画した領域を示している。そのため、安静呼吸時と同様に各フレーム画像における、基準画像の各小領域と同じ画素位置の領域をその小領域に対応する領域としてしまうと、後述する解析結果が診断に利用し得ないものとなる。そこで、このような場合には、各フレーム画像間の対応点を抽出する対応点抽出処理ローカルマッチング処理及び非線形歪変換処理ワーピング処理を行って、肺野領域の同一部分が描画された領域を各フレーム画像間で対応付ける。また、安静呼吸時に於いて、処理速度をあまり重視せず、解析精度をより高めたい場合にも、これらの処理を行うこととしてもよい。ローカルマッチング処理では、まず、撮影順が1番最初のフレーム画像から抽出した肺野領域を、例えば、2mm角の矩形からなる小領域に分割する。次いで、撮影順が1番のフレーム画像をP1、これと隣接するフレーム画像撮影順が隣接するフレーム画像即ち、時間的に隣接するフレーム画像。以下同様。をP2とし、P2に、P1の各小領域の探索領域を設定する。ここで、P2の探索領域は、P1における各小領域における中心点の座標をx,yとすると、同一の中心点x,yをもち、P1の小領域よりも縦横の幅が大きくなるように設定する例えば、1.5倍。そして、P1の各領域毎に、P2の探索範囲で最もマッチング度合いが高くなる領域を求めることで、P1の各小領域に対するP2上での対応位置を算出する。マッチング度合いとしては、最小二乗法や相互相関係数を指標に用いる。そして、P2の肺野領域をP1の各小領域の対応位置で分割する。次いで、P2を、新たにP1とみなし、撮影順がP2の次のフレーム画像を新たなP2とみなして、P1の各小領域におけるP2の対応位置を算出する。以上の処理を繰り返すことで、各フレーム画像の各小領域が隣接するフレーム画像のどの位置に対応するかが求まる。求めた処理結果は、制御部31のRAMに記憶される。次いで、ワーピング処理が行われる。具体的には、撮影順が1番のフレーム画像をP1、これと撮影順が隣接する時間的に隣接するフレーム画像をP2とし、上記ローカルマッチング処理で算出された隣接するフレーム画像間の各小領域の対応位置に基づいて、各小領域毎にP1からP2へのシフトベクトルを算出する。次いで、算出されたシフトベクトルを多項式でフィッティングして、この多項式を用いて各小領域における各画素のシフトベクトルを算出する。そして、算出された各画素のシフトベクトルに基づいて、ワーピング処理を行い、P2の各小領域内の各画素の位置をP1のフレーム画像の対応する画素の位置にシフトする。次いで、ワーピング処理されたP2を、新たにP1とみなし、撮影順がP2の次のフレーム画像を新たなP2とみなして、上記処理を行う。以上の処理を撮影順の早い隣接フレーム画像間から順次繰り返すことで、全てのフレーム画像の各小領域の位置を撮影順が1のフレーム画像この処理での基準画像に略一致させることが可能となる。各フレーム画像間における小領域の位置の対応関係は、制御部31のRAMに記憶される。次いで、一連のフレーム画像のステップS22で対応付けられた各小領域毎に解析が行われ、吸気の特徴量及び呼気の特徴量がそれぞれ算出されるステップS23。ここでは、呼気の特徴量及び吸気の特徴量として、例えば、呼気期及び吸気期のそれぞれにおける各小領域毎のフレーム間差分値微分値の代表値絶対値の最大値を算出する。フレーム間差分値は、そのフレーム画像が撮影されたタイミングにおける信号変化量を示す値である。呼吸により息を吸ったり吐いたりすれば、その息の流れに応じて肺の密度が変化し、これによってX線透過量つまり、画素の出力信号値が変化する。よって、信号変化量は、そのタイミングにおける気流速度を示す値とみなすことができる。なお、代表値としては、絶対値の最大値に限らず、中央値、平均値、最頻値としてもよい。具体的には、まず、各フレーム画像の各小領域内の画素の信号値平均信号値が算出される。次いで、撮影順が隣接するフレーム画像間で各小領域の信号値の差分を算出するフレーム間差分処理が行われる。ここでは、各小領域毎に、フレーム番号NとN1Nは1、2、3のフレーム画像について、N1-Nの差分値が算出される。そして、呼気期のフレーム間差分値の最大値絶対値の最大値が呼気の特徴量として取得され、吸気期のフレーム間差分値の最大値絶対値の最大値が吸気の特徴量として取得される。フレーム間差分値の最大値絶対値の最大値は、最大微分値に相当する。ここで、各小領域内のフレーム間差分値の符号が正の期間が吸気期、負の期間が呼気期である。次いで、各小領域毎の吸気の特徴量と呼気の特徴量の比の値吸気の特徴量/呼気の特徴量が算出されるステップS24。ここでは、吸気期のフレーム間差分値の最大値/呼気期のフレーム間差分値の最大値が算出される。次いで、算出された各小領域毎の吸気の特徴量/呼気の特徴量の値のヒストグラムが作成されるとともに、肺野全体における吸気の特徴量/呼気の特徴量の傾向を示す指標値ここでは、平均値、標準偏差が算出されるステップS25。ヒストグラムの縦軸のカウント数は、肺野内の全小領域の数で除算することにより正規化することが好ましい。次いで、各小領域について求められた吸気の特徴量/呼気の特徴量の値を、予め記憶部32に記憶されている、吸気の特徴量/呼気の特徴量の値と表示時のパラメータ値との変換テーブルに基づいて表示用パラメータ値に変換し、変換したパラメータ値にて、基準画像例えば、安静呼気位のフレーム画像の各小領域を表示した画像が作成されるステップS26。変換テーブルは、例えば、正常/異常重症度1~nの各カテゴリに特徴量を分類するときの各カテゴリの特徴量の大きさの範囲を規定する閾値各カテゴリの閾値と色相、明度、輝度、透明度の何れかを一対一に対応付けたテーブルである。ここで、表示時のパラメータ値の変換テーブルとしては、特徴量の大きさに対する識別能を高めるために、各カテゴリの閾値と色相を対応付けることが好ましい。このとき、例えば、数個例えば、5~6個の色相を上記各カテゴリの閾値と対応させ、さらにその間の特徴量の値に対しては中間の色相を割り当てるグラデーションさせることで識別能の高い表示とすることができる。表示用パラメータ値に基づいて色付けされた画像は、基準画像のフレーム画像にオーバーレイ表示してもよい。そして、表示部34に、作成されたヒストグラム、及び作成された静止画像等が並列表示されステップS27、画像解析処理は終了する。ヒストグラムの領域は、上述の吸気の特徴量/呼気の特徴量の値と表示時のパラメータ値との変換テーブルに基づいて、基準画像の肺野領域の各小領域と同様の基準により色付けして表示される。図8A~図8Cに、ステップS27において表示部34に表示される表示画面の例を示す。図8Aは、正常者の肺野の動態画像を解析した解析結果を表示した表示画面である。図8Bは、COPD閉塞性疾患の肺野の動態画像を解析した解析結果を表示した表示画面である。図8Cは、混合性疾患の肺野の動態画像を解析した解析結果を表示した表示画面である。図8A~図8Cに示すように、ステップS27においては、半切サイズのフレーム画像から抽出された肺野領域内における各小領域2mm角の矩形サイズの吸気の特徴量/呼気の特徴量の値のヒストグラム34aと、各小領域の吸気の特徴量/呼気の特徴量を一覧表示した静止画像34bと、ヒストグラム34a及び静止画像34bで表示されている色相と吸気の特徴量/呼気の特徴量の値との関係を示す表示34cと、肺野全体での吸気の特徴量/呼気の特徴量の傾向を示す指標値34dが表示されている。吸気の特徴量/呼気の特徴量の分布の特徴を把握することが可能となる。また、図8A~図8Cに示すように、ヒストグラム34aは、吸気の特徴量/呼気の特徴量の値の大きさに応じて、横軸の領域が5つの色相で色分けして表示される。これにより、医師は、ヒストグラムを一瞥するだけで容易に肺野内における吸気の特徴量/呼気の特徴量の分布を把握することができる。また、各小領域の吸気の特徴量/呼気の特徴量を示す静止画像34bにおいて、各小領域が吸気の特徴量/呼気の特徴量の値に応じてヒストグラムの色分けと同じ基準で色分けされて表示されるので、医師が肺野内の局所的な異常箇所閉塞性部分、拘束性部分を容易に把握することが可能となる。更に、肺野全体での吸気の特徴量/呼気の特徴量の傾向を示す指標値34dとして、その平均値及び標準偏差を算出して画面上に併せて表示することで、肺野全体での吸気の特徴量/呼気の特徴量の傾向を数値として医師に提供することができる。ここで、呼気期のフレーム間差分値の最大値絶対値を呼気の特徴量、吸気期のフレーム間差分値の最大値絶対値を吸気の特徴量とした場合、正常者の肺野では、肺野全体の吸気の特徴量/呼気の特徴量の平均値は0.9~1.2、標準偏差は0.10~0.22程度であることが知られている。よって、ステップS27で図8Aに示す表示画面が表示された場合、医師は、撮影した肺野が正常であることを容易に把握することができる。一方、COPD閉塞性疾患の肺野では、肺野全体の吸気の特徴量/呼気の特徴量の平均値は0.9~1.2に収まらず正常者に比べて大きく、標準偏差も、0.10~0.22に収まらない正常者に比べて大きいことが知られている。よって、ステップS27で図8Bに示す表示画面が表示された場合、医師は、撮影した肺野がCOPDであることを容易に把握することができる。一方、混合性肺疾患の肺野では、肺野全体の吸気の特徴量/呼気の特徴量の値が0.66以下のデータ数、及び1.5以上のデータ数がともに増加することが知られている。よって、ステップS27で図8Cに示す表示画面が表示された場合、医師は、撮影した肺野が混合性疾患であることを容易に把握することができる。このように、胸部診断支援システム100においては、COPD閉塞性肺疾患、混合性肺疾患等の病態を特定できるような、有用な診断情報を医師に提供することが可能となる。呼気の特徴量及び吸気の特徴量としては、上記の例以外の他の特徴量を用いても良い。例えば、呼気の特徴量を呼吸1サイクルにおける呼気期に相当するフレーム画像数呼気時間、吸気の特徴量を呼吸1サイクルにおける吸気期に相当するフレーム画像数吸気期間としてもよい。ここで、肺の換気機能が正常である場合、吸気時間と呼気時間はほぼ同じ長さ、もしくは若干呼気時間が長くなる筈である。よって、呼気期に相当するフレーム画像数/吸気期に相当するフレーム画像数の値を見れば、医師は肺疾患の疑いがあるか否かを把握することができる。特に、呼気期のフレーム画像数/吸気期のフレーム画像数>1.5の領域は、呼気で換気がしにくく、取り込んだ空気を排出するのが遅れてしまう閉塞性部分であることがわかる。なお、吸気期のフレーム間差分値の最大値/呼気期のフレーム間差分値の最大値≒呼気時間/吸気時間の関係が成り立つため、医師は、呼気の特徴量を呼気期のフレーム間差分値の最大値、吸気の特徴量を吸気期のフレーム間差分値の最大値とした場合と同様の判断基準により、正常、COPD、混合性肺疾患の識別を行うことが可能である。また、1呼吸サイクル分の各フレーム画像において、各小領域の画素の信号値平均信号値を算出し、各小領域毎に呼吸1サイクルにおける信号値の最小値及び最大値を求め、求めた最小値をその領域の呼気の特徴量、最大値を吸気の特徴量としてもよい。信号値の最大値と最小値は正常な箇所では両者の値の差が大きくなり、異常のある箇所では両者の差が非常に小さくなると考えられる。異常のある箇所では肺胞の動きが悪くなるため、肺胞の密度変化が小さくなるからである。よって、医師は、信号値の最大値/信号値の最小値のヒストグラムを参照し、平均値及び標準偏差を確認することで、肺野が正常であるか疾患であるかの判断材料とすることができる。例えば、肺野全体の信号値の最大値/信号値の最小値の平均値が1より大きく、標準偏差が小さい場合、肺の機能は正常であると判断することができる。一方、肺野全体の信号値の最大値/信号値の最小値の平均値が1に近い値であり、標準偏差が大きい場合、肺の機能に疾患があると判断することができる。その他、吸気の特徴量/呼気の特徴量の傾向を示す指標値として、平均値、標準偏差の他に、カウント数ブロック小領域数がヒストグラムのピークとなる吸気の特徴量/呼気の特徴量の値、もしくは、ピークのカウント数ブロック数、もしくは、吸気の特徴量/呼気の特徴量の値が所定以上所定以下となるカウント数の割合を指標値として用いても良い。もしくは、これらの複数の指標値を組み合わせて新たな指標値を作成しても良い。また、吸気の特徴量/呼気の特徴量を対数に変換し、対数変換後の値に対して、上記ヒストグラムの作成、指標値の算出、表示時パラメータとの変換テーブルに基づく色付け画像の作成を行っても良い。以上説明したように、胸部診断支援システム100の診断用コンソール3によれば、制御部31は、胸部の動態を撮影して得られた複数のフレーム画像のそれぞれから肺野領域を抽出し、抽出した肺野領域を複数の小領域に分割し、複数のフレーム画像間における各小領域を対応付ける。次いで、複数のフレーム画像の対応する小領域毎に解析を行って、小領域毎の吸気の特徴量及び呼気の特徴量を算出し、算出された吸気の特徴量と呼気の特徴量との比の値を算出し、当該算出された比の値のヒストグラムを作成する。そして、作成したヒストグラムを表示部34に表示する。従って、吸気の特徴量と呼気の特徴量の比が肺野内においてどのように分布しているのかを示す、診断に有利となる情報であるヒストグラムを医師に提供することができるので、医師は肺野の換気に係る病態を容易に把握することが可能となる。吸気の特徴量として、吸気期に相当するフレーム画像群の時間的に隣接するフレーム画像間のフレーム間差分値の絶対値の最大値を算出し、呼気の特徴量として、呼気期に相当するフレーム画像群の時間的に隣接するフレーム画像間のフレーム間差分値の絶対値の最大値を算出し、算出された吸気の特徴量と呼気の特徴量の比の分布をヒストグラムで表示することで、撮影された肺野が正常であるか、閉塞性疾患であるか、混合性疾患であるかを医師が判断する際の有用な情報を提供することが可能となる。吸気の特徴量として、呼吸1サイクルにおける吸気期に相当するフレーム画像の数を算出し、呼気の特徴量として、呼吸1サイクルにおける呼気期に相当するフレーム画像の数を算出し、算出された吸気の特徴量と呼気の特徴量の比の分布をヒストグラムで表示することで、撮影された肺野が正常であるか、閉塞性疾患であるか、混合性疾患であるかを医師が判断する際の有用な情報を提供することが可能となる。吸気の特徴量として、呼吸1サイクル分の複数のフレーム画像における画素信号値の最大値を算出し、呼気の特徴量として、呼吸1サイクル分の複数のフレーム画像における画素信号値の最小値を呼気の特徴量として算出し、算出された吸気の特徴量と呼気の特徴量の比の分布をヒストグラムで表示することで、撮影された肺野が正常であるか異常であるかを医師が判断する際の有用な情報を提供することができる。複数のフレーム画像の肺野領域を小領域に分割する際、複数のフレーム画像のうち一のフレーム画像の肺野領域を複数の小領域に分割し、他のフレーム画像の肺野領域を一のフレーム画像の各小領域と同じ画素位置の小領域に分割し、複数のフレーム画像間の同じ画素位置の各小領域を互いに対応付けることで、フレーム画像の肺野の小領域への分割及びフレーム画像間での小領域の対応付けを高速に行うことができる。その結果、解析処理の高速化を図ることが可能となる。複数のフレーム画像の肺野領域を小領域に分割する際、複数のフレーム画像のうち一のフレーム画像の肺野領域を複数の小領域に分割し、他のフレーム画像のそれぞれにおいて一のフレーム画像における各小領域とマッチング度合いの高い領域を抽出し、他のフレーム画像の肺野領域を抽出された各領域に分割し、複数のフレーム画像間のマッチング度合いの高い各小領域を互いに対応付けることで、精度良く解析を行うことが可能となる。<第2の実施の形態>以下、図面を参照して本発明の第2の実施の形態を詳細に説明する。ただし、発明の範囲は、図示例に限定されない。ここで、従来、胸部静止画像の放射線撮影においては、例えば、特開平10-98587号公報公知文献1に記載のように、散乱放射線を除去するための散乱線除去グリッド以下、グリッドというを使用して撮影が行われている。また、近年、特開2005-312775号公報公知文献2に記載のように、FPDflatpaneldetector等の放射線検出器を利用して胸部の動態撮影を行い、診断に応用するシステムが提案されているが、このような動態撮影用のシステムにおいても、静止画撮影と同様に、グリッドを使用して撮影を行う構成となっている。公知文献2においては、被写体の被曝線量が増大しがちな動態撮影に於いて、被写体の被曝線量低減のために、静止画撮影1回当たりの照射線量と一連の動態撮影に係るトータルの照射線量とが略一致するように、動態撮影におけるフレーム画像1枚当たりの照射線量を算出している。ここで、従来の手法に従いグリッドを使用すると放射線検出器検出器へ到達する線量が減弱されるため、公知文献2に記載のような放射線照射条件で動態撮影を行うと、各フレーム画像に到達する線量が大幅に不足し、S/N比が低下する。そのため、一連のフレーム画像を用いて動態解析を行った場合、正確な解析ができず、照射線量を変更して再撮影を行うこととなり、結果として被写体のトータルの被曝線量が増大してしまう恐れもある。そこで、本件発明者等が鋭意検討を行った結果、検出器の絶対的出力信号濃度信号を観察する静止画撮影とは異なり、動態撮影に於いては、各フレーム相互の相対的出力信号差濃度差を視認したり、例えばフレーム間差分等の相対的出力信号差に基づく特徴量の解析を行ったりするので、各フレーム画像に略一様に散乱線影響が重畳していても、その目的を達成することが可能であることを見出した。従って、動態撮影においては、グリッド無しで撮影を行うことができる。グリッド無しでの撮影において、検出器に到達する線量を従来並みに保つ場合、被写体への照射線量はグリッドによる減衰分を予め加算しておく必要が無く、従来手法に比べ、被写体の被曝線量を半減することが可能となる。また、被写体の被曝線量照射線量を一定のままとすると、各フレーム画像あたりの検出器到達線量を増加することができるので、より精度を有する解析を行うことができる。第2の実施の形態においては、動態撮影時にグリッド無しで撮影を行う放射線撮影システム200について説明する。放射線撮影システム200の構成まず、構成を説明する。図9に、本実施の形態における放射線撮影システム200の全体構成を示す。この放射線撮影システム200は、胸部診断支援システムとして機能するシステムである。図9に示すように、放射線撮影システム200は、撮影装置4と、撮影用コンソール5とが通信ケーブル等により接続され、撮影用コンソール5と、診断用コンソール6とがLANLocalAreaNetwork等の通信ネットワークNT2を介して接続されて構成されている。放射線撮影システム200を構成する各装置は、DICOMDigitalImageandCommunicationsinMedicine規格に準じており、各装置間の通信は、DICOMに則って行われる。撮影装置4の構成撮影装置4は、被写体M人体の胸部に放射線を照射して被写体Mの動態撮影又は静止画撮影を行う装置である。動態撮影とは、被写体Mに対し、X線等の放射線をパルス的に連続照射して複数の画像を取得即ち、連続撮影することをいう。動態撮影では、例えば、呼吸運動に伴う肺の膨張及び収縮の形態変化、心臓の拍動等の、周期性サイクルを持つ被写体Mの動態を撮影する。この連続撮影により得られた一連の画像を動態画像と呼ぶ。また、動態画像を構成する複数の画像のそれぞれをフレーム画像と呼ぶ。静止画撮影とは、従来のフィルム方式やCR方式と同様に撮影部位の濃度分解能に基づく診断に使用されるもので、被写体Mに対し、X線等の放射線を1回照射して一枚の静止画像を取得することをいう。撮影装置4は、図9に示すように、放射線源41、放射線照射制御装置42、放射線検出部43、読取制御装置44、保持部45、グリッド46等を備えて構成されている。放射線源41は、単射及び連射パルス照射が可能な放射線発生装置である。即ち、静止画撮影と動態撮影の双方に対応した放射線発生装置である。放射線源41は、被写体Mを挟んで放射線検出部43と対向する位置に配置され、放射線照射制御装置42の制御に従って、被写体Mに対し放射線X線を照射する。放射線照射制御装置42は、撮影用コンソール5に接続されており、撮影用コンソール5から入力された放射線照射条件に基づいて放射線源41を制御して放射線撮影を行う。撮影用コンソール5から入力される放射線照射条件は、例えば、連続照射時のパルスレート、パルス幅、パルス間隔、1撮影あたりの撮影フレーム数、X線管電流の値、X線管電圧の値、フィルタ種等である。パルスレートは、1秒あたりの放射線照射回数であり、後述するフレームレートと一致している。パルス幅は、放射線照射1回当たりの放射線照射時間である。パルス間隔は、連続撮影において、1回の放射線照射開始から次の放射線照射開始までの時間であり、後述するフレーム間隔と一致している。放射線検出部43は、動態撮影及び静止画撮影対応のFPD等により構成される。FPDは、例えば、ガラス基板等を有しており、基板上の所定位置に、放射線源41から照射されて少なくとも被写体Mを透過した放射線をその強度に応じて検出し、検出した放射線を電気信号に変換して蓄積する複数の画素がマトリックス状に配列されている。各画素は、例えばTFTThinFilmTransistor等のスイッチング部により構成されている。放射線検出部43は、図10に示すように、被写体Mを挟んで放射線源41と対向するように保持部45によって保持されている。読取制御装置44は、撮影用コンソール5に接続されている。読取制御装置44は、撮影用コンソール5から入力された画像読取条件に基づいて放射線検出部43の各画素のスイッチング部を制御して、当該各画素に蓄積された電気信号の読み取りをスイッチングしていき、放射線検出部43に蓄積された電気信号を読み取ることにより、画像データ静止画像又はフレーム画像を取得する。そして、読取制御装置44は、取得した画像データを撮影用コンソール5に出力する。画像読取条件は、例えば、フレームレート、フレーム間隔、画素サイズビニングサイズ、画像サイズマトリックスサイズ等である。フレームレートは、1秒あたりに取得するフレーム画像数であり、パルスレートと一致している。フレーム間隔は、連続撮影において、1回のフレーム画像の取得動作開始から次のフレーム画像の取得動作開始までの時間であり、パルス間隔と一致している。ここで、放射線照射制御装置42と読取制御装置44は互いに接続され、互いに同期信号をやりとりして放射線照射動作と画像の読み取りの動作を同調させるようになっている。この他、後述するオフセット補正に用いるオフセット補正係数を算出するための少なくとも1つの暗画像を取得するダーク読取時は、放射線照射動作と同期せず、放射線が照射されない状態で、リセット~蓄積~データ読取~リセットの一連の画像の読み取り動作を行うが、一連の動態撮影前、一連の動態撮影後のいずれかのタイミングで行うようにしてもよい。保持部45は、図10に示すように、検出器保持部451を有し、撮影時に放射線検出部43を放射線源41及び被写体Mに対向するように保持する。また、保持部45は、放射線検出部43よりも被写体側放射線源41側に散乱放射線を除去するためのグリッド46を装着するためのグリッド装着部452を有している。即ち、保持部45は、グリッド46が着脱可能に構成されている。グリッド装着部452には、図11に示すようにグリッド46が装着されているか否かを検知するためのグリッド装着検知MSマイクロスイッチ453が設けられており、保持部45は、グリッド装着検知MS453の検知信号を読取制御装置44に出力する。また、保持部45には、図10に示すように被写体Mが所定距離離間して存在するか否かを検知するための被写体検知センサ454が設けられており、保持部45は、被写体検知センサ454の検知信号を読取制御装置44を介して撮影用コンソール5に出力する。撮影用コンソール5の構成撮影用コンソール5は、放射線照射条件や画像読取条件を撮影装置4に出力して撮影装置4による放射線撮影及び放射線画像の読み取り動作を制御するとともに、撮影装置4により取得された静止画像又は動態画像に基づく画像、例えば間引きビニング処理されたプレビュ画像や諧調処理等を施した処理済画像等を適宜作成し、撮影実施者によるポジショニングの確認や診断に適した画像であるか否かの確認用に表示する。撮影用コンソール5は、図9に示すように、制御部51、記憶部52、操作部53、表示部54、通信部55を備えて構成され、各部はバス56により接続されている。制御部51は、CPUCentralProcessingUnit、RAMRandomAccessMemory等により構成される。制御部51のCPUは、操作部53の操作に応じて、記憶部52に記憶されているシステムプログラムや各種処理プログラムを読み出してRAM内に展開し、展開されたプログラムに従って後述する撮影制御処理を始めとする各種処理を実行し、撮影用コンソール5各部の動作や、撮影装置4の放射線照射動作及び読み取り動作を集中制御する。記憶部52は、不揮発性の半導体メモリやハードディスク等により構成される。記憶部52は、制御部51で実行される各種プログラムやプログラムにより処理の実行に必要なパラメータ、或いは処理結果等のデータを記憶する。例えば、記憶部52は、図12A~図12Bに示す撮影制御処理を実行するための撮影制御処理プログラムを記憶している。各種プログラムは、読取可能なプログラムコードの形態で格納され、制御部51は、当該プログラムコードに従った動作を逐次実行する。また、記憶部52は、動態撮影用、静止画撮影用のそれぞれの放射線照射条件及び画像読取条件を記憶している。操作部53は、カーソルキー、数字入力キー、及び各種機能キー等を備えたキーボードと、マウス等のポインティングデバイスを備えて構成され、キーボードに対するキー操作やマウス操作により入力された指示信号を制御部51に出力する。また、操作部53は、表示部54の表示画面にタッチパネルを備えても良く、この場合、タッチパネルを介して入力された指示信号を制御部51に出力する。表示部54は、LCDLiquidCrystalDisplayやCRTCathodeRayTube等のモニタにより構成され、制御部51から入力される表示信号の指示に従って、操作部53からの入力指示やデータ等を表示する。通信部55は、LANアダプタやモデムやTATerminalAdapter等を備え、通信ネットワークNT2に接続された各装置との間のデータ送受信を制御する。診断用コンソール6の構成診断用コンソール6は、撮影用コンソール5から静止画像、又は動態画像の一連のフレーム画像を取得し、取得した画像や後述するヒストグラム等の診断支援情報を表示して医師が読影診断するためのコンピュータ装置である。診断用コンソール6は、図9に示すように、制御部61、記憶部62、操作部63、表示部64、通信部65を備えて構成され、各部はバス66により接続されている。制御部61は、CPU、RAM等により構成される。制御部61のCPUは、操作部63の操作に応じて、記憶部62に記憶されているシステムプログラムや、各種処理プログラムを読み出してRAM内に展開し、展開されたプログラムに従って、後述する画像解析処理を始めとする各種処理を実行し、診断用コンソール6各部の動作を集中制御する。記憶部62は、不揮発性の半導体メモリやハードディスク等により構成される。記憶部62は、制御部61で画像解析処理を実行するための画像解析処理プログラムを始めとする各種プログラムやプログラムにより処理の実行に必要なパラメータ、或いは処理結果等のデータを記憶する。これらの各種プログラムは、読取可能なプログラムコードの形態で格納され、制御部61は、当該プログラムコードに従った動作を逐次実行する。操作部63は、カーソルキー、数字入力キー、及び各種機能キー等を備えたキーボードと、マウス等のポインティングデバイスを備えて構成され、キーボードに対するキー操作やマウス操作により入力された指示信号を制御部61に出力する。また、操作部63は、表示部64の表示画面にタッチパネルを備えても良く、この場合、タッチパネルを介して入力された指示信号を制御部61に出力する。表示部64は、LCDやCRT等のモニタにより構成され、制御部61から入力される表示信号の指示に従って、操作部63からの入力指示やデータ等を表示する。通信部65は、LANアダプタやモデムやTA等を備え、通信ネットワークNT2に接続された各装置との間のデータ送受信を制御する。放射線撮影システム200の動作次に、上記放射線撮影システム200における動作について説明する。撮影装置4、撮影用コンソール5の動作まず、撮影装置4、撮影用コンソール5による撮影動作について説明する。図12A~図12Bに、撮影用コンソール5の制御部51において実行される撮影制御処理を示す。撮影制御処理は、制御部51と記憶部52に記憶されている撮影制御処理プログラムとの協働により実行される。まず、撮影実施者により撮影用コンソール5の操作部53が操作され、撮影対象被写体Mの患者情報患者の氏名、身長、体重、年齢、性別等及び撮影種別動態撮影、静止画撮影の区別の入力が行われるステップS31。この他、通信ネットワークNT2に接続された図示しない他装置から送信され、通信部55を経由して受信したデータも、同様の入力情報とすることができる。次いで、入力された撮影種別が動態撮影であるか静止画撮影であるかが判断されるステップS32。入力された撮影種別が動態撮影であると判断されるとステップS32;YES、動態撮影用の放射線照射条件が記憶部52から読み出されて放射線照射制御装置42に設定されるとともに、動態撮影用の画像読取条件が記憶部52から読み出されて読取制御装置44に設定されるステップS33。入力された撮影種別が静止画撮影であると判断されるとステップS32;NO、静止画撮影用の放射線照射条件が記憶部52から読み出されて放射線照射制御装置42に設定されるとともに、静止画撮影用の画像読取条件が記憶部52から読み出されて読取制御装置44に設定されるステップS37。本実施の形態においては、動態撮影用としてグリッド46を装着しない状態でグリッド無しでの撮影を前提とした放射線照射条件が設定され、静止画撮影用としてグリッド46を装着した状態でグリッド有りでの撮影を前提とした放射線照射条件が設定される。具体的には、動態撮影用の放射線照射条件は、フレーム画像1枚ごとの撮影においてグリッド46を使用しないで放射線検出部43に到達する線量が静止画撮影においてグリッド46を使用して放射線検出部43に到達する線量と同じになるような放射線照射条件が設定される。即ち、動態撮影では、1枚あたりの照射線量が静止画撮影に比べて低くなるように放射線照射条件が設定される。ここで、人体等の被写体Mに対して放射線を照射して放射線撮影を行う場合、体内を透過する放射線が体内組織で散乱される。そして、このような散乱放射線が放射線検出部43に入射すると放射線画像にノイズが発生する。そこで、病変部分の検出や病変部分の観察等に用いられ個々の画素の絶対的出力値信号値が重要な静止画の撮影においては、放射線検出部43の被写体側の面、すなわち被写体を透過した放射線が入射する側の面に、グリッド46を設けて撮影が行われることが好ましい。グリッド有りで撮影を行う場合、グリッド46により放射線検出部43に到達する線量が減弱される例えば、露出倍数2のグリッドでは約1/2に減弱されるため、この減弱分を加算した分の放射線を照射する必要がある。従来、動態撮影においても、静止画撮影と同様にグリッドを用いた撮影が行われていた。しかしながら、動態撮影の場合は、静止画撮影に比べて撮影するフレーム画像枚数が多いため、1枚当たりに放射線源41から照射する照射線量を静止画撮影と同等とした場合、被写体Mの被曝量が多くなってしまうという問題がある。被写体Mの被曝量低減のため1回の静止画撮影と、一連の動態撮影のトータルの照射線量を同じとする技術も開示されているが、この場合は個々のフレーム画像の線量が不足気味となり、S/N比が低下してしまう。そこで、本件発明者等は、鋭意検討を重ねた結果、動態撮影により得られる胸部動態画像は主として換気機能や血流等の動態解析に使用されるものであり、これらの解析においては、グリッドを配置せずに撮影を行った動態画像を用いても、グリッド有りで撮影を行った動態画像を用いた場合と略同等の結果を得ることができることを見出した。言い換えると、グリッド有無に限らず、放射線検出器に到達する線量が同じであれば、略同等の解析結果を得ることができる図13参照。詳細後述。本実施の形態の放射線撮影システム200においては、この知見に基づき、静止画撮影ではグリッド46を装着してグリッド有りで撮影を行い、動態撮影ではグリッド46を装着せずにグリッド無しで撮影を行う。そして、グリッド46を使用した静止画撮影時と、グリッド46無しでの動態撮影の1フレーム当たりの撮影とで放射線検出部43に到達する線量が略同等となる放射線照射条件で撮影を行うことにより、放射線検出部43への到達線量は従来と略同等に維持しつつ被写体Mの被曝線量の低減する構成となっている。または、1回の静止画撮影と、一連の動態撮影のトータルの照射線量を同じとする放射線照射条件とし、静止画撮影ではグリッド有りで撮影を行い、動態撮影ではグリッド無しで撮影を行うことで、従来方式に比べて各フレーム画像のS/N比を向上させて解析精度の向上を図ることとしてもよい。動態撮影用の放射線照射条件及び画像読取条件が設定されると、グリッド装着検知MS453からの出力に基づいて、グリッド46がグリッド装着部452に装着されていない状態であるか否かが判断されるステップS34。グリッド46がグリッド装着部452に装着されていない状態であると判断されるとステップS34;YES、操作部53の操作による放射線照射の指示が待機されるステップS35。ここで、撮影実施者は、安静呼吸の動態を撮影するために被験者被写体Mに楽にするように指示し、安静呼吸を促す。撮影準備が整った時点で、操作部53を操作して放射線照射指示を入力する。操作部53により放射線照射指示が入力されるとステップS35;YES、放射線照射制御装置42及び読取制御装置44に撮影開始指示が出力され、動態撮影が開始されるステップS36。即ち、放射線照射制御装置42に設定されたパルス間隔で放射線源41により放射線が照射され、放射線検出部43によりフレーム画像が取得される。予め定められたフレーム数の撮影が終了すると、制御部51により放射線照射制御装置42及び読取制御装置44に撮影終了の指示が出力され、撮影動作が停止される。撮影されるフレーム数は、少なくとも1呼吸サイクルが撮影できる枚数である。一方、静止画撮影用の放射線照射条件及び画像読取条件が設定されると、グリッド装着検知MS453と被写体検知センサ454からの出力に基づいて、グリッド46がグリッド装着部452に装着された状態であるか否かが判断されるステップS38。制御部51は、このステップS38によって、グリッド46が装着されずに静止画撮影が行われることのないように制御している。グリッド46がグリッド装着部452に装着されていると判断されるとステップS38;YES、操作部53の操作による放射線照射の指示が待機されるステップS39。ここで、撮影実施者は、被験者に対し、息を吸った後に止めるように指示する。撮影準備が整った時点で、操作部53を操作して放射線照射指示を入力する。操作部53により放射線照射指示が入力されるとステップS39;YES、放射線照射制御装置42及び読取制御装置44に撮影開始指示が出力され、静止画撮影が行われるステップS40。動態撮影又は静止画撮影が終了すると、撮影により取得された画像各フレーム画像又は静止画像は順次撮影用コンソール5に入力され、補正処理が行われるステップS41。ステップS41の補正処理においては、オフセット補正処理、ゲイン補正処理、欠陥画素補正処理の3つの補正処理が必要に応じて行われる。まず最初に、取得された画像に対してオフセット補正処理が行われ、取得された画像に重畳された暗電流に起因するオフセット値が除去される。オフセット補正処理では、例えば、取得した画像の各画素値濃度値。以下、信号値というから、予め記憶されたオフセット補正係数を減算する処理が行われる。ここで、オフセット補正係数は、予め放射線非照射時に取得した複数の画像を平均化した画像である。次いで、ゲイン補正処理が行われ、取得された画像の各画素に対応する各検出素子の個体差や読み出しアンプのゲインムラによって生じる画素毎のばらつきが除去される。ゲイン補正処理では、例えば、オフセット補正後の画像に、予め記憶されたゲイン補正係数を乗算する処理が行われる。ここで、ゲイン補正係数は、放射線検出部43に一様に放射線を照射した時に取得した複数のオフセット補正済み画像を平均化した画像と、このときの放射線照射条件で期待される出力信号値の関係から、補正後の各画素の信号値が一様となるように予め算出され、記憶された係数である。次いで、欠陥画素補正処理が行われ、周囲の画素と比較して感度が非線形な画素や、感度がない欠落画素が除去される。欠陥画素補正処理では、例えば、予め記憶された欠陥画素位置情報マップに従って、欠陥画素位置情報マップに登録された各欠陥画素において、欠陥画素の信号値をその近傍の欠陥でない画素の信号値の平均値で置き換える処理が行われる。ここで、欠陥画素位置情報マップは、放射線検出部43に一様に放射線を照射した時に取得したオフセット補正、ゲイン補正済みのフレーム画像から、予め複数の欠陥画素が認識され、その欠陥画素の位置が登録されたマップである。上記オフセット補正係数及びゲイン補正係数、欠陥画素位置情報マップは、ビニングやダイナミックレンジ等の収集モードに応じて、予め、それぞれ最適な値が記憶されており、それぞれの収集モードにおいて対応する最適な値が読み出されるようになっている。なお、これらの補正処理は処理時間の短縮化の為、省略することも可能である。また、必要に応じて、取得された画像に対し、例えば間引き処理や階調処理等を施すことが好ましい。次いで、補正処理後の画像が記憶部52に記憶されるとともにステップS42、表示部54に表示されるステップS43。動態撮影を行った場合は、各フレーム画像が撮影順を示す番号と対応付け記憶部52に記憶される。ここで、取得された画像を記憶する直前に、各画素の信号値を真数から対数に変換する対数変換処理を行ってから記憶しても良い。撮影実施者は、表示された動態画像によりポジショニング等を確認し、撮影により診断に適した画像が取得された撮影OKか、再撮影が必要撮影NGか、を判断する。そして、操作部53を操作して、判断結果を入力する。尚、撮影により取得された各フレーム画像は、全撮影の終了後に纏めて入力するようにしても良い。操作部53の所定の操作により撮影OKを示す判断結果が入力されるとステップS44;YES、静止画撮影で取得された静止画像又は動態撮影で取得された一連のフレーム画像のそれぞれに、画像を識別するための識別IDや、患者情報、検査対象部位、放射線照射条件、画像読取条件、撮影順を示す番号、撮影日時、撮影時のグリッド有無を示す情報グリッド有無情報等の情報が付帯され例えば、DICOM形式で画像データのヘッダ領域に書き込まれ、通信部55を介して診断用コンソール6に送信されるステップS45。そして、本処理は終了する。一方、操作部53の所定の操作により撮影NGを示す判断結果が入力されるとステップS44;NO、記憶部52に記憶された一連のフレーム画像が削除されステップS46、本処理は終了する。尚、このケースに於いては、再撮影が実行されることとなる。診断用コンソール6の動作次に、診断用コンソール6における動作について説明する。診断用コンソール6においては、通信部65を介して撮影用コンソール5から静止画像が受信され、操作部63により当該画像の表示指示が入力されると、表示部64に受信された静止画像が表示され、医師の診断に供される。一方、通信部65を介して撮影用コンソール5から動態画像の一連のフレーム画像が受信され、操作部63により動態解析が指示されると、制御部61と記憶部62に記憶されている画像解析処理プログラムとの協働により画像解析処理が実行される。第2の実施の形態における画像解析処理は、第1の実施の形態で図3で示した画像解析処理と大まかな処理ステップは同じであるが詳細が異なるので、以下に図3及び図13以降の図を参照して説明する。まず、各フレーム画像から肺野領域が抽出されるステップS51。次いで、各フレーム画像の肺野領域が複数の小領域に分割され、各フレーム画像の小領域が互いに対応付けられるステップS52。ステップS51~S52の具体的な処理については、第1の実施の形態において図3のステップS21~S22で説明したものと同様であるので説明を援用する。次いで、一連のフレーム画像のステップS52で対応付けられた各小領域毎に解析が行われ、吸気の特徴量及び呼気の特徴量がそれぞれ算出されるステップS53。ここでは、呼気の特徴量及び吸気の特徴量として、例えば、呼気期及び吸気期のそれぞれにおける各小領域毎のフレーム間差分値微分値の代表値絶対値の最大値を算出する。フレーム間差分値は、そのフレーム画像が撮影されたタイミングにおける信号変化量を示す値である。呼吸により息を吸ったり吐いたりすれば、その息の流れに応じて肺の密度が変化し、これによってX線透過量つまり、画素の出力信号値が変化する。よって、信号変化量は、そのタイミングにおける気流速度を示す値とみなすことができる。なお、代表値としては、絶対値の最大値に限らず、中央値、平均値、最頻値としてもよい。具体的には、まず、各フレーム画像の各小領域内の画素の信号値平均信号値が算出される。次いで、各フレーム画像間で対応づけられた各小領域に対して、時間軸方向のフィルタ処理が施される。このフィルタ処理は、血流等の高周波数な信号変化を除去し、換気による信号値の時間変化を抽出するための処理であり、例えば、小領域毎の信号値の時間変化に対して、安静呼吸画像群ではカットオフ周波数0.7Hz、深呼吸画像群ではカットオフ周波数0.5Hzでローパスフィルタリングを行う。ここで、ローパスフィルタのカットオフ周波数は、固定値とするよりも、撮影された動態画像毎に最適化することが好ましい。例えば、前述のように一連のフレーム画像の横隔膜の位置を解析し、安静換気の場合は安静呼気位および安静吸気位となるフレームを検出し、安静呼気位のフレームとその次の安静吸気位のフレームの間のフレーム数から吸気期の時間を求め、その逆数に所定の係数を乗じた値をカットオフ周波数としたローパスフィルタを実施する。このとき安静換気の場合、自動的に設定されるカットオフ周波数は、0.2~1.0Hzの間に制限することが好ましい。また、ステップS31において、別途測定された安静時における1分間の呼吸数、脈拍数等のバイタルサインを患者情報として入力しておき、これらの値からカットオフ周波数が算出されるようにしてもよい。例えば、患者情報として入力された1分間の呼吸数を1秒間の呼吸数に変換し、その呼吸数に所定の係数を乗じた値をカットオフ周波数としてローパスフィルタを実施してもよい。さらには、入力された1分間の脈拍数を1秒間の脈拍数に変換し、1秒間の呼吸数と1秒間の心拍数の平均値をカットオフ周波数としてローパスフィルタを実施するようにしてもよい。フィルタ処理に次いで、撮影順が隣接するフレーム画像間で各小領域の信号値の差分を算出するフレーム間差分処理が行われる。ここでは、各小領域毎に、フレーム番号NとN1Nは1、2、3のフレーム画像について、N1-Nの差分値が算出される。そして、呼気期のフレーム間差分値の最大値絶対値の最大値が呼気の特徴量として取得され、吸気期のフレーム間差分値の最大値絶対値の最大値が吸気の特徴量として取得される。フレーム間差分値の最大値絶対値の最大値は、最大微分値に相当する。ここで、各小領域内のフレーム間差分値の符号が正の期間が吸気期、負の期間が呼気期である。次いで、各小領域毎の吸気の特徴量と呼気の特徴量の比の値吸気の特徴量/呼気の特徴量が算出されるステップS54。ここでは、吸気期のフレーム間差分値の最大値/呼気期のフレーム間差分値の最大値最大流速比と呼ぶが算出される。次いで、算出された各小領域毎の吸気の特徴量/呼気の特徴量の値のヒストグラムが作成されるとともに、肺野全体における吸気の特徴量/呼気の特徴量の傾向を示す指標値ここでは、平均値、標準偏差が算出されるステップS55。ヒストグラムの縦軸のカウント数は、肺野内の全小領域の数で除算することにより正規化することが好ましい。次いで、各小領域について求められた吸気の特徴量/呼気の特徴量の値を、予め記憶部62に記憶されている、吸気の特徴量/呼気の特徴量の値と表示時のパラメータ値との変換テーブルに基づいて表示用パラメータ値に変換し、変換したパラメータ値にて、基準画像例えば、安静呼気位のフレーム画像の各小領域を表示した画像が作成されるステップS56。変換テーブルは、例えば、正常/異常重症度1~nの各カテゴリに特徴量を分類するときの各カテゴリの特徴量の大きさの範囲を規定する閾値各カテゴリの閾値と色相、明度、輝度、透明度の何れかを一対一に対応付けたテーブルである。ここで、表示時のパラメータ値の変換テーブルとしては、特徴量の大きさに対する識別能を高めるために、各カテゴリの閾値と色相を対応付けることが好ましい。このとき、例えば、数個例えば、5~6個の色相を上記各カテゴリの閾値と対応させ、その間の特徴量の値に対しては中間の色相を割り当てるグラデーションさせることで識別能の高い表示とすることができる。表示用パラメータ値に基づいて色付けされた画像は、基準画像のフレーム画像にオーバーレイ表示してもよい。なお、記憶部62には、グリッド有りで撮影された動態画像に対応した変換テーブルと、グリッド無しで撮影された動態画像に対応する変換テーブルが記憶されており、ステップS56及び次のステップS57においては、一連のフレーム画像に付帯されているグリッド有無情報に基づいて、グリッド有りで撮影されたかグリッド無しで撮影されたかが判断され、その判断結果に応じた変換テーブルを用いて色付けが行われる。そして、表示部64に、作成されたヒストグラム、及び作成された静止画像等が並列表示されステップS57、画像解析処理は終了する。ヒストグラムの領域は、上述の吸気の特徴量/呼気の特徴量の値と表示時のパラメータ値との変換テーブルに基づいて、基準画像の肺野領域の各小領域と同様の基準により色付けして表示される。ここで、グリッド有りで撮影したときとグリッド無しで撮影したときの動態解析への影響について説明する。図13は、或る正常の肺野をグリッド有りとグリッド無しで撮影した動態画像の動態解析結果の比較を示す図である。図13においては、グリッド有りとグリッド無しで撮影した動態画像の各小領域を吸気期のフレーム間差分値の最大値最大吸気気流速度に基づいて色付けした画像、呼気期のフレーム間差分値の最大値最大呼気気流速度に基づいて色付けした画像、最大流速比に基づいて色付けした画像、及び最大流速比のヒストグラム、を解析結果として示している。図13は、以下の撮影条件で撮影を行った動態画像の解析結果である。検出器サイズ40×30cm、検出器画素サイズ194μm、グリッドピッチ80line/cm、グリッド比12:1、管球~検出器間距離2m、フレーム数75フレーム約10秒間撮影、トータル被曝線量検出器への到達線量を一定とした場合の被写体被曝線量グリッド有り0.24mGy、グリッド無し0.14mGyまた、最大吸気気流速度、最大呼気気流速度、最大流速比の大きさと色図13においては濃度で示すとの変換テーブルは、両者を比較するために同一ここでは、グリッド有り用の変換テーブルを用いている。同一被写体を撮影した場合のグリッド有りの動態画像とグリッド無しの動態画像とでは、最大吸気気流速度、最大呼気気流速度、最大流速比はほぼ同等であるが、図13に示すように、撮影系の特性等により若干相違が生じることがある。例えば図13においてはグリッド無しで撮影された動態画像のヒストグラムの形状はグリッド有りのものと比較して幅広な形状となっている。そのため、例えば最大流速比によって肺野内の領域やヒストグラムの領域を正常、異常1~nのカテゴリに分類して色付けする場合、グリッド有りとグリッド無しで同じ閾値変換テーブルを用いると、同じ最大流速比であっても異なる色表示、即ち異なる重症度への分類がなされる場合があり、好ましくない。そこで、図13に示すように、グリッド有りと無しとで診断に影響する差異が生じる場合においては、グリッド有無に応じて特徴量の分類に用いる閾値変換テーブルを変更する必要がある。グリッド有りで撮影した動態画像を解析した結果とグリッド無しで撮影した動態画像を解析した結果とでどの程度の差異が生じるかは、撮影系の特性や解析内容等によって異なる。よって、撮影系によって、同じ被写体をグリッド有り無しで撮影した複数の画像を解析し、その結果を用いて帰納法的にグリッド有りの動態画像に使用する閾値とグリッド無しの動態画像に使用する閾値を算出しておくことが好ましい。なお、本実施の形態においては、撮影装置4は動態画像をグリッド無しで撮影するよう制御しているため、グリッド無しで撮影された動態画像用の閾値が記憶されていれば十分とも考えられる。しかし、グリッド有りで動態画像を撮影可能な撮影系に診断用コンソール6が接続される場合も想定される。この場合、グリッド有無の撮影条件が異なっていると、判断を誤ることになる。そこで、本実施の形態においては、動態画像を構成する各フレーム画像の付帯情報にグリッド有無情報を付帯させておき、制御部61は、このグリッド有無情報に基づいて、撮影時のグリッド有無に応じた閾値を用いた解析アルゴリズムで解析を行っている。また、同一患者の過去の画像の解析結果と現在の解析結果とを比較するような場合も想定される。互いにグリッド有無情報が付帯されている場合、一方の解析結果を他方の解析結果に合わせるべく、何れか一方の解析結果の値を補正して使用すること例えば、両者の変換テーブルの閾値の差分を解析結果としての特徴量に加算又は減算する等が好ましい。図14A~図14Cに、ステップS57において表示部64に表示される表示画面の例を示す。図14Aは、グリッド有りで撮影された正常者の肺野の動態画像を解析した解析結果を表示した表示画面である。図14Bは、COPD閉塞性疾患の肺野の動態画像を解析した解析結果を表示した表示画面である。図14Cは、混合性疾患の肺野の動態画像を解析した解析結果を表示した表示画面である。図14A~図14Cに示すように、ステップS57においては、半切サイズのフレーム画像から抽出された肺野領域内における各小領域2mm角の矩形サイズの吸気の特徴量/呼気の特徴量の値のヒストグラム64aと、各小領域の吸気の特徴量/呼気の特徴量を一覧表示した静止画像64bと、ヒストグラム64a及び静止画像64bで表示されている色相と吸気の特徴量/呼気の特徴量の値との関係を示す表示64cと、肺野全体での吸気の特徴量/呼気の特徴量の傾向を示す指標値64dが表示されている。また、図14A~図14Cに示すように、ヒストグラム64aは、吸気の特徴量/呼気の特徴量の値の大きさに応じて、横軸の領域が6つの色相で色分けして表示される。これにより、医師は、ヒストグラムを一瞥するだけで容易に肺野内における吸気の特徴量/呼気の特徴量の分布を把握することができる。また、各小領域の吸気の特徴量/呼気の特徴量を示す静止画像64bにおいて、各小領域が吸気の特徴量/呼気の特徴量の値に応じてヒストグラムの色分けと同じ基準で色分けされて表示されるので、医師が肺野内の局所的な異常箇所閉塞性部分、拘束性部分を容易に把握することが可能となる。更に、肺野全体での吸気の特徴量/呼気の特徴量の傾向を示す指標値64dとして、その平均値及び標準偏差を算出して画面上に併せて表示することで、肺野全体での吸気の特徴量/呼気の特徴量の傾向を数値として医師に提供することができる。尚、本表示例に於いては、縦軸は、2mm角の小ブロックの出現個数/全ブロック数430mm÷2mm×354mm÷2mm=38055個の割合としているが、図8に示すように、縦軸を2mm角の小ブロックの出現個数自体としても良い。動態撮影に使用する検出器のサイズが複数存在する場合には、本表示例のような全体の領域個数で正規化した割合表示によって同一症例に対する解析結果表示の統一性が達成され、解析結果に対する医師等の誤認を招き難いので好ましい。ここで、呼気期のフレーム間差分値の最大値絶対値を呼気の特徴量、吸気期のフレーム間差分値の最大値絶対値を吸気の特徴量とした場合、正常者の肺野では、グリッド有りで動態撮影した動態画像で解析を行った場合、肺野全体の吸気の特徴量/呼気の特徴量の平均値は0.9~1.2、標準偏差は0.10~0.22程度であることが知られている。よって、ステップS57で図14Aに示す表示画面が表示された場合、医師は、撮影した肺野が正常であることを容易に把握することができる。一方、COPD閉塞性疾患の肺野では、グリッド有りで動態撮影した動態画像で解析を行った場合、肺野全体の吸気の特徴量/呼気の特徴量の平均値は0.9~1.2に収まらず正常者に比べて大きく、標準偏差も、0.10~0.22に収まらない正常者に比べて大きいことが知られている。よって、ステップS57で図14Bに示す表示画面が表示された場合、医師は、撮影した肺野がCOPDであることを容易に把握することができる。一方、混合性肺疾患の肺野では、グリッド有りで動態撮影した動態画像で解析を行った場合、肺野全体の吸気の特徴量/呼気の特徴量の値が0.66以下のデータ数、及び1.5以上のデータ数がともに増加することが知られている。よって、ステップS57で図14Cに示す表示画面が表示された場合、医師は、撮影した肺野が混合性疾患であることを容易に把握することができる。このように、放射線撮影システム200においては、吸気の特徴量/呼気の特徴量の傾向を示す指標値から、COPD閉塞性肺疾患、混合性肺疾患等の換気不均等の病態もしくはその重症度を特定できるような、有用な診断情報を医師に提供することが可能となる。なお、吸気の特徴量/呼気の特徴量の傾向を示す指標値等を利用して被写体Mの正常/異常を判断することとしてもよい。この場合においても、グリッド有無情報に基づいて正常/異常を判断するための閾値を変更することが好ましい。例えば、上述の肺野全体の吸気の特徴量/呼気の特徴量の平均値に基づいて正常/異常と判断する場合、グリッド有りでは平均値0.9~1.2を正常とし、グリッド無しでは平均値0.8~1.3を正常として判断することが好ましい。呼気の特徴量及び吸気の特徴量としては、上記の例以外の他の特徴量を用いても良い。例えば、呼気の特徴量を呼吸1サイクルにおける呼気期に相当するフレーム画像数呼気時間、吸気の特徴量を呼吸1サイクルにおける吸気期に相当するフレーム画像数吸気期間としてもよい。ここで、肺の換気機能が正常である場合、吸気時間と呼気時間はほぼ同じ長さ、もしくは若干呼気時間が長くなる筈である。よって、呼気期に相当するフレーム画像数/吸気期に相当するフレーム画像数の値を見れば、医師は肺疾患の疑いがあるか否かを把握することができる。特に、呼気期のフレーム画像数/吸気期のフレーム画像数>1.5の領域は、呼気で換気がしにくく、取り込んだ空気を排出するのが遅れてしまう閉塞性部分であることがわかる。なお、吸気期のフレーム間差分値の最大値/呼気期のフレーム間差分値の最大値≒呼気時間/吸気時間の関係が成り立つため、医師は、呼気の特徴量を呼気期のフレーム間差分値の最大値、吸気の特徴量を吸気期のフレーム間差分値の最大値とした場合と同様の判断基準により、正常、COPD、混合性肺疾患の識別を行うことが可能である。また、1呼吸サイクル分の各フレーム画像において、各小領域の画素の信号値平均信号値を算出し、各小領域毎に呼吸1サイクルにおける信号値の最小値及び最大値を求め、求めた最小値をその領域の呼気の特徴量、最大値を吸気の特徴量としてもよい。信号値の最大値と最小値は正常な箇所では両者の値の差が大きくなり、異常のある箇所では両者の差が非常に小さくなると考えられる。異常のある箇所では肺胞の動きが悪くなるため、肺胞の密度変化が小さくなるからである。よって、医師は、信号値の最大値/信号値の最小値のヒストグラムを参照し、平均値及び標準偏差を確認することで、肺野が正常であるか疾患であるかの判断材料とすることができる。例えば、肺野全体の信号値の最大値/信号値の最小値の平均値が1より大きく、標準偏差が小さい場合、肺の機能は正常であると判断することができる。一方、肺野全体の信号値の最大値/信号値の最小値の平均値が1に近い値であり、標準偏差が大きい場合、肺の機能に疾患があると判断することができる。その他、吸気の特徴量/呼気の特徴量の傾向を示す指標値として、平均値、標準偏差の他に、カウント数ブロック小領域数がヒストグラムのピークとなる吸気の特徴量/呼気の特徴量の値、もしくは、ピークのカウント数ブロック数、もしくは、吸気の特徴量/呼気の特徴量の値が所定以上所定以下となるカウント数の割合を指標値として用いても良い。もしくは、これらの複数の指標値を組み合わせて新たな指標値を作成しても良い。以上説明したように、放射線撮影システム200によれば、静止画撮影時にはグリッド46を使用し、動態撮影時にはグリッド46を使用せずに撮影を行う。これにより、グリッドによる線量の減弱がなくなる分、同じ被写体への照射線量に対して放射線検出部43への到達する線量が増えるので、高精細な動画解析を行うこと可能となる。また、放射線検出部43への到達線量を一定とすると、従来のように動態撮影時にグリッドを使用した場合に比べて被写体への照射線量被写体の被曝線量の低減を図ることが可能となる。また、動態撮影により得られた各フレーム画像には、グリッド有無情報が付帯され、診断用コンソール6の制御部61は、グリッド有無情報に基づいて解析アルゴリズムを変更する、例えば、解析に用いる閾値を変更するので、動態撮影時のグリッド有無に拘わらず正確な解析結果を提供することが可能となる。なお、上述した本実施の形態における記述は、本発明に係る好適な一例であり、これに限定されるものではない。例えば、グリッド無しで撮影を行う場合、散乱線影響を緩和して、より良い動態観察や動態解析を行うには、被写体Mを放射線検出部43と離間させて撮影を行うことが好ましい。被写体Mを放射線検出部43と離間させて撮影を行うと、エアギャップにより散乱線の影響を低減することができる。離間により拡大撮影となるが、動態解析結果には影響はない。特許3861572号公報に記載のように、位相コントラスト撮影の撮影条件としてもよい。なお、この際には、被写体が離間しているか密着しているかを被写体検知センサ454で検知し、グリッド無で密着の場合には動態撮影を禁止するように制御しても良い。また、例えば、上記実施の形態においては、吸気の特徴量/呼気の特徴量の解析結果として、最大流速比をヒストグラムにより表示することとしたが、これに限定されない。例えば、図15に示すように、吸気の特徴量/呼気の特徴量の傾向を示す指標値の一つをX軸、もう一つをY軸としたグラフに、各々の指標値に対する正常、異常の閾値TH1を示し、そのグラフ上に、動態画像から算出した肺野全体での吸気の特徴量/呼気の特徴量の傾向を示す指標値をプロットした図を解析結果としてもよい。図15は、X軸を吸気期のフレーム間差分値の最大値/呼気期のフレーム間差分値の最大値の平均値、Y軸をその標準偏差とし、動態画像から算出した肺野全体における吸気の特徴量/呼気の特徴量の平均値、標準偏差の指標値をプロットしたグラフである。このようなグラフを用いて吸気の特徴量/呼気の特徴量の傾向を示す指標値を表示することにより、プロットされた点から閾値TH1までの距離によって異常の程度が視覚的に容易に把握することが可能となる。また、例えば、2つの指標値例えば、吸気の特徴量/呼気の特徴量の平均値、標準偏差を線形結合した値を新たな指標値とし、図16に示すように2つの指標値の一方をX軸、他方をY軸としたグラフに、新たな指標値2つの指標値を線形結合した指標値を重症度により分類するための閾値th1~th4を示し、このグラフ上に動態画像から算出した新たな指標値をプロットしてもよい。線形結合の例としては、平均値及び標準偏差に対する多数の測定値データから主成分分析により算出された第1主成分を指標値とすることができる。このようなグラフを用いることにより、異常の程度が視覚的により容易に把握することが可能となる。また、例えば、複数M個の指標値に対する多数N個の測定値データより、N×M個のデータに対する共分散行列の最大固有値を算出することで第一主成分を算出し、算出された主成分を指標値として用いてもよい。また、吸気の特徴量/呼気の特徴量の傾向を示す指標値の他に、吸気の特徴量、もしくは呼気の特徴量のそれぞれに対して傾向を示す指標値を算出してもよい。例えば、図17に示すように、小領域毎に算出した吸気もしくは呼気の特徴量に対して、左右肺野内をそれぞれ上中下の3つに分割した計6つの領域の各々において変動係数=標準偏差/平均値を算出し、変動係数の大きさに応じた色相又は輝度又は彩度によって、6つの領域を色付けし表示した静止画像を作成してもよい。このような表示をすることで、不均等換気の分布を把握することが容易となり、不均等換気が行われている部分が区域的であるかもしくはびまん性であるかが容易に判定可能となる。また、吸気の特徴量/呼気の特徴量を対数に変換し、対数変換後の値に対して、上記ヒストグラムの作成、指標値の算出、表示時パラメータとの変換テーブルに基づく色付け画像の作成を行っても良い。また、上記実施の形態においては、動態撮影と静止画撮影の双方の撮影が可能な撮影装置を含むシステムについて説明したが、動態撮影のみが可能な撮影装置においても、グリッドを使用せずに撮影を行うことは、被写体の被曝量低減又は画質の向上に有効である。また、例えば、上記の説明では、本発明に係るプログラムのコンピュータ読み取り可能な媒体としてハードディスクや半導体の不揮発性メモリ等を使用した例を開示したが、この例に限定されない。その他のコンピュータ読み取り可能な媒体として、CDROM等の可搬型記録媒体を適用することが可能である。また、本発明に係るプログラムのデータを通信回線を介して提供する媒体として、キャリアウエーブも適用される。その他、胸部診断支援システム100、放射線撮影システム200を構成する各装置の細部構成及び細部動作に関しても、本発明の趣旨を逸脱することのない範囲で適宜変更可能である。なお、明細書、請求の範囲、図面及び要約を含む2010年8月27日に出願された日本特許出願No.2010-190241号、及び2011年1月26日に出願された日本特許出願No.2011-013538号の全ての開示は、そのまま本出願の一部に組み込まれる。"