(57)【要約】【課題】医療用のモダリティーで撮影された動画から、ポイントとなる情報だけを容易に読影できるようにする。【解決手段】放射線画像処理装置１は、複数のフレーム画像からなる放射線動画の動画データを取得する画像取得手段と、画像取得手段が取得した動画データに所定の解析処理を加え、複数の解析フレーム画像からなる解析結果動画を生成する解析手段と、解析手段が生成した解析結果動画を構成する複数の解析フレーム画像の中から複数の特定解析フレーム画像を選択する選択手段と、選択手段によって選択された複数の特定解析フレーム画像における、共通する座標上にそれぞれ位置する各共通座標画素の信号値に基づいて一の演算信号値を導出する演算処理手段と、演算処理手段が座標ごとに生成した複数の演算信号値を用いて前記表示部に演算結果画像を表示する画像表示手段と、を備える。【選択図】図２
2018157968.tif 000002
【特許請求の範囲】
【請求項１】
画像を表示可能な表示部と、
複数のフレーム画像からなる放射線動画の動画データを取得する画像取得手段と、
前記画像取得手段が取得した動画データに所定の解析処理を加え、複数の解析フレーム画像からなる解析結果動画を生成する解析手段と、
前記解析手段が生成した解析結果動画を構成する複数の解析フレーム画像の中から複数の特定解析フレーム画像を選択する選択手段と、
前記選択手段によって選択された複数の特定解析フレーム画像における、共通する座標上にそれぞれ位置する各共通座標画素の信号値に基づいて一の演算信号値を導出する演算処理手段と、
前記演算処理手段が座標ごとに生成した複数の演算信号値を用いて前記表示部に演算結果画像を表示する画像表示手段と、を備えることを特徴とする放射線画像処理装置。
【請求項２】
前記画像表示手段は、
前記表示部の表示領域の少なくとも一部にサマリー表示領域を形成し、
前記サマリー表示領域に前記演算結果画像を表示させることが可能であることを特徴とする請求項１に記載の放射線画像処理装置。
【請求項３】
ユーザーが操作可能な操作部を有し、
前記選択手段は、ユーザーが前記操作部に対して行った操作に基づいて前記特定解析フレーム画像を選択することを特徴とする請求項１又は２に記載の放射線画像処理装置。
【請求項４】
前記選択手段は、予め設定された選択パターンに基づいて前記特定解析フレーム画像を自動的に選択することを特徴とする請求項１又は２に記載の放射線画像処理装置。
【請求項５】
ユーザーが操作可能な操作部を有し、
前記選択手段は、
ユーザーが前記操作部に対して行った操作に基づいて、予め用意された複数の選択パターンの中から何れかの選択パターンを設定し、
設定した選択パターンに基づいて前記特定解析フレーム画像を自動的に選択することを特徴とする請求項４に記載の放射線画像処理装置。
【請求項６】
前記演算処理手段は、前記複数の演算信号値を算出した後、前記演算結果画像の特定領域内に位置することになる複数の領域内画素の信号値に基づいて所定の定量値を導出する定量値算出手段を備え、
前記画像表示手段は、前記定量値算出手段が算出した定量値を前記サマリー表示領域１５Ａに表示させることが可能であることを特徴とする請求項１から５のいずれか一項に記載の放射線画像処理装置。
【請求項７】
前記演算処理手段は、前記複数の特定解析フレーム画像における、前記各共通座標画素の信号値の中から最大値、最小値もしくは最頻値を抽出する、当該共通座標画素の信号値における、選択した特定解析フレーム画像の数に対する最大値と最小値との差の割合、標準偏差、平均値もしくは積算値を算出する、又は前後する２枚の特定解析フレーム画像における各共通座標画素の信号値の差分を算出し、得られた各差分の中の最大値もしくは最小値を抽出することにより、前記一の演算信号値を導出することを特徴とする請求項１から６のいずれか一項に記載の放射線画像処理装置。
【請求項８】
前記定量値算出手段は、前記複数の領域内画素の信号値の中から最大値、最小値もしくは最頻値を抽出する、又は当該複数の領域内画素の信号値における、当該領域内画素の数に対する最大値と最小値との差の割合、標準偏差、平均値もしくは積算値を算出することにより、前記定量値を導出することを特徴とする請求項６又は７に記載の放射線画像処理装置。
【請求項９】
前記特定領域は、肺野内、心臓領域内、肺動脈内、大動脈内又は肺気管支内であることを特徴とする請求項８に記載の放射線画像処理装置。
【請求項１０】
放射線を照射する放射線照射装置と、
前記放射線照射装置から放射線を受けることで放射線画像の画像データを生成する放射線画像撮影装置と、を備え、
前記放射線照射装置による放射線の照射と、放射線画像撮影装置による画像データの読出しを連続的に繰り返すことにより、放射線動画の動画データを生成することが可能な放射線画像撮影システムであって、
前記放射線画像撮影装置と直接又は間接的に通信可能に接続された請求項１から９のいずれか一項に記載の放射線画像処理装置と、を備えることを特徴とする放射線画像撮影システム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、放射線画像処理装置、及びこの放射線画像処理装置を備えた放射線画像撮影システムに関する。
【背景技術】
【０００２】
医療用のモダリティー（例えば、ＣＲ（Computed Radiography）やＣＴ（Computed Tomography）等）を用いて撮影した動画は、表示装置に表示されることで医者の診断に供される。その際、撮影された動画をそのまま表示するだけでなく、様々な解析処理を施してから表示することがある。
従来、動画を構成する複数のフレーム画像の一覧を表示部に表示し、解析に用いるフレーム画像を一覧の中から一括選択できるようにする技術（特許文献１参照）や、画像を即時表示、参照表示あるいは転送する際に、それぞれに最適な画像処理を行う技術（特許文献２参照）が知られている。
【先行技術文献】
【特許文献】
【０００３】
【特許文献１】特開２０１２−１１０３９７号公報
【特許文献２】特開２０１２−１１０４６６号公報
【発明の概要】
【発明が解決しようとする課題】
【０００４】
しかしながら、特許文献１，２に記載されたような従来の技術では、ユーザーが元の動画と解析を加えた動画、あるいは異なる種類の解析を施した動画同士を見比べながら読影することを考えたとき、各動画のそれぞれの関心のある領域を交互に視点を変えながら読影する必要が出てくる。各動画は並行して再生が進むので、一方の動画で詳しく読影したい領域を見つけ、他方の動画の対応する領域を読影するのに視点を変えた時には、本来見たいフレーム画像よりも後のフレーム画像が表示されることになってしまう。このため、一方の動画から他方の動画へ移る際に、フレーム画像を戻して読影しなければならない場合があった。
つまり、動画は情報量が多い上に時間の経過に伴って変化するため、読影（診断）が容易でなく、ユーザーに過度の負担がかかっていた。
【０００５】
本発明は、上記の問題点に鑑みてなされたものであり、医療用のモダリティーで撮影された動画から、ポイントとなる情報だけを容易に読影できるようにすることを課題とする。
【課題を解決するための手段】
【０００６】
前記の問題を解決するために、本発明は、
複数のフレーム画像からなる放射線動画の動画データを取得する画像取得手段と、
前記画像取得手段が取得した動画データに所定の解析処理を加え、複数の解析フレーム画像からなる解析結果動画を生成する解析手段と、
前記解析手段が生成した解析結果動画を構成する複数の解析フレーム画像の中から複数の特定解析フレーム画像を選択する選択手段と、
前記選択手段によって選択された複数の特定解析フレーム画像における、共通する座標上にそれぞれ位置する各共通座標画素の信号値に基づいて一の演算信号値を算出する演算処理手段と、
前記演算処理手段が座標ごとに生成した複数の演算信号値を用いて前記表示部に演算結果画像を表示する画像表示手段と、を備えることを特徴とする。
【発明の効果】
【０００７】
本発明によれば、医療用のモダリティーで撮影された動画から、ポイントとなる情報だけを容易に読影できるようにすることができる。
【図面の簡単な説明】
【０００８】
【図１】本発明の第１〜第３実施形態に係る放射線画像撮影システムの全体構成を示すブロック図である。
【図２】第１実施形態に係る放射線画像処理装置が行う画像処理の概念図である。
【図３】同実施形態に係る放射線画像処理装置が表示する画像処理結果の一例である。
【図４】同実施形態に係る放射線画像処理装置が行う画像処理の概念図である。
【図５】同実施形態に係る放射線画像処理装置が表示する画像処理結果の一例である。
【図６】第２実施形態に係る放射線画像処理装置が行う画像選択の概念図である。
【図７】同実施形態に係る放射線画像処理装置が行う画像選択の概念図である。
【図８】同実施形態に係る放射線画像処理装置が行う画像選択の概念図である。
【図９】同実施形態に係る放射線画像処理装置が行う画像選択の概念図である。
【図１０】本発明の第３実施形態に係る放射線画像処理装置が行う画像選択の概念図である。
【図１１】同実施形態に係る放射線画像処理装置が表示する画像処理結果の一例である。
【図１２】同実施形態に係る放射線画像処理装置が表示する画像処理結果の一例である。
【図１３】第１〜第３実施形態の変形例に係る放射線画像処理装置が行う画像処理の概念図である。
【図１４】同変形例に係る放射線画像処理装置が表示する画像処理結果の一例である。
【図１５】同変形例に係る放射線画像処理装置が行う画像処理の概念図である。
【図１６】同変形例に係る放射線画像処理装置が行う画像処理の概念図である。
【図１７】同変形例に係る放射線画像処理装置が行う画像処理の概念図である。
【図１８】同変形例に係る放射線画像処理装置が表示する画像処理結果の一例である。
【発明を実施するための形態】
【０００９】
以下、図面を参照して本発明の実施の形態を詳細に説明する。ただし、発明の範囲は、図示例に限定されない。
【００１０】
＜第１実施形態＞
〔放射線画像撮影システムの構成〕
まず、本発明の第１実施形態に係る放射線画像撮影システム１００の構成について説明する。図１は本実施形態に係る放射線画像撮影システム１００の機能的構成を表すブロック図である。
【００１１】
本実施形態の放射線画像撮影システム１００は、放射線画像処理装置１や、医療用のモダリティー２等で構成されている。
また、放射線画像撮影システム１００には、必要に応じて、図示しないコンソールやサーバー（例えば、医療用画像管理システム（Picture Archiving and Communication System：ＰＡＣＳ）を構成するもの）等が接続される。
放射線画像撮影システム１００を構成する各装置は、ＤＩＣＯＭ（Digital Image and Communications in Medicine）規格に準じており、各装置間の通信は、ＤＩＣＯＭに則って行われる。
【００１２】
モダリティー２は。放射線照射装置２１と、放射線画像撮影装置２２と、を備えている。
放射線照射装置２１は、図示を省略するが、放射線を生成可能な回転陽極や回転陽極に電子ビームを照射するフィラメント等を有する放射線源や、設定された管電圧や管電流、照射時間（ｍＡｓ値）等に応じた線量の放射線を放射線源から照射させるジェネレーター等を備えている。
そして、ユーザーの操作に基づいて、放射線画像撮影装置２２へ放射線を照射するようになっている。
【００１３】
放射線画像撮影装置２２は、図示を省略するが、放射線を受けることで線量に応じた電荷を蓄積する複数の放射線検出素子が二次元状（マトリクス状）に配列された基板や、各放射線検出素子に蓄積された電荷を画像データとして読み出す読み出し回路、外部装置と通信したり画像データを送信したりするための通信部等を備えている。
そして、放射線画像撮影装置２２は、放射線照射装置２から放射線の照射を受け、画像データを読み出すと、その画像データを、通信部を介して直ちに外部へ送信するようになっている。
【００１４】
このような放射線照射装置２及び放射線画像撮影装置２２を備えたモダリティー２は、放射線照射装置２１による放射線の照射と、放射線画像撮影装置２２による電荷の蓄積・画像データの読み出しと、を連続的に繰り返すことで、複数のフレーム画像からなる放射線動画の動画データを生成することが可能となっている。
なお、放射線照射装置２１と放射線画像撮影装置２２は、ＣＴのように一体化されたものでもよい。
【００１５】
〔放射線画像撮影システムの構成〕
次に、上記放射線画像撮影システム１００を構成する放射線画像処理装置１の構成について説明する。
放射線画像処理装置１は、コンピューター又は専用の制御装置として構成されており、図１に示したように、制御部１１、通信部１２、記憶部１３、操作部１４、表示部１５、を備えて構成され、各部はバス１６により接続されている。
【００１６】
制御部１１は、ＣＰＵ、ＲＡＭ等により構成される。制御部１１のＣＰＵは、操作部１４の操作に応じて、記憶部１３に記憶されているシステムプログラムや、各種処理プログラムを読み出してＲＡＭ内に展開し、展開されたプログラムに従って、後述する画像解析処理を始めとする各種処理を実行したり、表示部１５の表示内容を制御したりする等、放射線画像処理装置１各部の動作を集中制御する。
【００１７】
通信部１２は、ＬＡＮアダプタやモデムやＴＡ等を備え、通信ネットワークＮＴに接続された放射線画像撮影装置２２やコンソール、サーバー等をはじめとする各装置との間のデータ送受信を制御する。
【００１８】
記憶部１３は、不揮発性の半導体メモリやハードディスク等により構成される。記憶部１３は、制御部１１で各種処理（例えば、後述する画像解析処理や画像演算処理）を実行するためのプログラム（画像解析処理プログラムや画像演算処理プログラム等）や、プログラムにより処理の実行に必要なパラメータ、処理結果、あるいは、各種処理を施す前の元動画や静止画、各種処理を施した後の動画や静止画のデータを記憶する。これらの各種プログラムは、読取可能なプログラムコードの形態で格納されている。
【００１９】
操作部１４は、カーソルキー、数字入力キー、及び各種機能キー等を備えたキーボードや、マウス等のポインティングデバイス等によりユーザーが操作可能に構成されており、キーボードに対するキー操作あるいはマウス操作により入力された指示信号を制御部１１に出力する。
なお、操作部１４は、表示部１５の表示画面に備えられたタッチパネルで構成されていても良く、この場合、タッチパネルを介して入力された指示信号を制御部１１に出力する。
【００２０】
表示部１５は、ＬＣＤやＣＲＴ等のモニタにより構成され、制御部１１から入力される表示信号の指示に従って、各種画像や、操作部１４からの入力指示、データ等を表示することが可能となっている。
【００２１】
〔放射線画像処理装置の動作〕
次に、本実施形態に係る放射線画像処理装置１の動作について説明する。図２，４は放射線画像処理装置１における画像処理の概念図であり、図３，５は処理結果の表示例である。
なお、本実施形態に係る放射線画像処理装置１は、静止画像の処理や表示も可能であるが、ここでは、主に動画を処理・表示する場合について説明する。また、各図には、肺野領域の画像を例示したが、本実施形態に係る放射線画像処理装置１は、肺野以外の部位を撮影した動画についても処理・表示することが可能である。
【００２２】
放射線画像処理装置１の制御部１１は、操作部１４を介してユーザーから処理を施す前の元の放射線動画を表示する旨の指示を受けると、複数のフレーム画像からなる放射線動画の元動画データ（ローデータ等）の取得（通信部１２を介してモダリティー２、コンソールもしくはサーバーから受信、又は記憶部１３から読み込み）を行い、準備した元動画データに基づく元動画又は図２（ａ）に示したような元動画を構成する複数の元フレーム画像Ｉ０の一覧（画像群））を表示部１５に表示することが可能となっている。すなわち、制御部１１は、本発明における画像取得手段をなす。なお、図２（ａ）では、元フレーム画像を符号Ｉ０で示したが、以下の説明では、元動画であるか元フレーム画像であるかにかかわらず、元動画データに基づいて表示部１５に表示されるものをＩ０で示すこととする。
【００２３】
また、制御部１１は、操作部１４を介してユーザーから元動画データに動態解析を加える旨の指示を受けると、元動画データ（ローデータ等）の取得を行い、記憶部１３に記憶されている画像解析処理プログラムとの協働により元動画データに各種解析処理（換気解析や血流解析等）を施し、その解析結果である解析動画データを生成することが可能となっている。
そして、制御部１１は、生成した解析動画データに基づく解析結果動画又は図２（ｂ）に示したような解析結果動画を構成する複数の解析フレーム画像Ｉ１の一覧）を、自動で又はユーザーからの指示に基づいて表示部１５に表示することが可能となっている。すなわち、制御部１１は、本発明における解析手段をなす。なお、図２（ｂ）では、解析フレーム画像を符号Ｉ１で示したが、以下の説明では、解析結果動画であるか解析フレーム画像であるかにかかわらず、解析動画データに基づいて表示部１５に表示されるものをＩ１で示すこととする。
【００２４】
また、制御部１１は、操作部１４を介してユーザーから解析動画データに演算処理を加える旨の指示を受けると、解析動画データの取得を行い、記憶部１３に記憶されている画像演算処理プログラムとの協働により解析動画データに各種演算処理を施し、その演算結果である演算画像データを生成することが可能となっている。
なお、元動画データに演算処理を施して演算画像データを生成することも可能である。
【００２５】
演算処理は、表示部１５に表示された、解析結果動画を構成する複数の解析フレーム画像Ｉ１の一覧の中から複数の解析フレーム画像Ｉ１（解析結果動画を構成する全ての解析フレーム画像Ｉ１であってもよいし、その一部であってもよい。）が選択されることにより開始される。以下、選択された解析結果フレーム画像Ｉ１を特定解析フレーム画像Ｉ１と称する。
本実施形態において、特定解析フレーム画像Ｉ１の選択は、ユーザーが操作部１４に対して行った操作に基づいて行われるようになっている。具体的には、例えば図２（ｂ）に示したように、表示部１５に表示された解析結果動画を構成する解析フレーム画像Ｉ１の一覧の中から、ユーザーが処理対象とする特定解析フレーム画像Ｉ１を所定領域（例えば後述するサマリー表示領域）へドラッグ及びドロップしたり、解析フレーム画像Ｉ１の端部に設けられたチェックボックスＢにチェックＣを入れたりすることで選択される。すなわち、制御部１１は、本発明における選択手段をなす。
なお、選択する解析フレーム画像Ｉ１は連続している必要は無く、間を飛ばして任意の解析フレーム画像Ｉ１を選択することも可能である。
【００２６】
演算処理では、選択された複数の特定解析フレーム画像Ｉ１における共通する座標上にそれぞれ位置する各画素の信号値に基づいて一の演算信号値を導出する。そして、この処理を各座標（（１，１）から（ｍ、ｎ）まで）について繰り返し行い、座標ごとに生成した複数の演算信号値を用いて演算結果画像Ｉ２となる演算画像データを生成する。すなわち、制御部１１は、本発明における演算処理手段をなす。以下、各特定解析フレーム画像Ｉ１における共通する座標上にそれぞれ位置する各画素のことを共通座標画素と称する。
【００２７】
演算信号値の導出の仕方としては、例えば、最大値（ＭＩＰ）処理、最小値（ＭｉｎＰ）処理、Ｐｋ−Ｐｋ値処理、標準偏差処理、最頻値処理、積算値処理、平均値処理、最大微分値処理、最小微分値処理等がある。制御部１１は、自動で又はユーザーからの指示に基づいて、上記の処理の中から何れかを設定して演算信号値を導出する。
【００２８】
最大値処理は、複数の特定解析フレーム画像Ｉ１における、各共通座標画素の信号値の中から最大値を抽出して演算信号値とする処理である。この処理により得られた解析画像からは、解析結果の値（換気や血流）の最大値が分かるので、換気や血流の有無を一目で把握することができる。具体的には、無気肺や肺塞栓領域において、換気機能や血流機能がどこまで回復しているのかを把握しやすくなる。また、処理コストを低く抑えることができる。
【００２９】
最小値処理は、複数の特定解析フレーム画像Ｉ１における、各共通座標画素の信号値の中から最小値を抽出して演算信号値とする処理である。この処理により得られた解析画像からは、解析結果の値（換気や血流）の最小値が分かるので、換気機能や血流機能がどこまで低下しているのかを把握しやすくなる。また、処理コストを低く抑えることができる。
【００３０】
最頻値処理は、複数の特定解析フレーム画像Ｉ１における、各共通座標画素の信号値の中から最頻値を抽出して演算信号値とする処理である。この処理により得られた解析画像からは、解析結果の値（換気や血流）の最頻値が分かるので、特異的なノイズまたはアーチファクトの影響を低減させた上で、換気や血流の有無を一目で把握できるようになる。例えば複数呼吸分、複数心拍分の解析結果に対して最頻値処理を実施することで、呼吸や心拍のばらつき影響を低減させることができる。
【００３１】
Ｐｋ−Ｐｋ値処理は、複数の特定解析フレーム画像Ｉ１における、各共通座標画素の信号値の中から最大値と最小値を抽出してその差分または比を算出し、特定解析フレーム画像の数に対する最大値と最小値との差の割合を算出して演算信号値とする処理である。この処理により得られた解析画像からは、解析結果の値（換気や血流）の変動量が分かるので、換気や血流の有無を一目で把握できるようになる。例えば、無気肺や肺塞栓領域の把握を行いやすくなる（信号の変動に強い、ロバストな結果を出すことができる）。
【００３２】
標準偏差処理は、複数の特定解析フレーム画像Ｉ１における、各共通座標画素の信号値の標準偏差を算出して演算信号値とする処理である。この処理により得られた解析画像からは、解析結果の値（換気や血流）の変動量が（複雑な計算アルゴリズム不要で）分かるので、換気や血流の有無を一目で把握できるようになる。例えば、無気肺や肺塞栓領域の把握を行いやすくなる。また、処理コストを低く抑えることができる。
【００３３】
積算値処理は、複数の特定解析フレーム画像Ｉ１における、各共通座標画素の信号値の積算値を算出して演算信号値とする処理である。この処理により得られた解析画像からは、解析結果の値（換気や血流）のある一定時間分の積算結果が分かるので、ノイズまたはアーチファクトの影響を低減させることができ、換気や血流の有無の把握が行いやすくなる。また、複数呼吸分、複数心拍分の解析結果を積算することで、呼吸や心拍のばらつき影響を低減させることができる。また、ある一定時間の間に計測した値を積算した値、または画像を最終結果として出力する他モダリティー(例えばシンチグラフィ)と、類似した物理量を導出できる為、積算値処理を行うことで上記のような他モダリティーの情報と比較が行いやすくなる。その結果、診断の精度が向上する。更に、画素値は線量と関係が有る場合がある為、総照射線量（総被曝量）の推定が行いやすくなる為、被検者の被ばく管理などに有用である。
【００３４】
平均値処理は、複数の特定解析フレーム画像Ｉ１における、各共通座標画素の信号値の平均値を算出して演算信号値とする処理である。この処理により得られた解析画像からは、ノイズまたはアーチファクトの影響を低減させた上で、換気や血流の有無を一目で把握できるようになる。例えば複数呼吸分、複数心拍分の解析結果に対して平均値処理を実施することで、呼吸や心拍のばらつき影響を低減させることができる。
【００３５】
最大微分値処理は、複数の特定解析フレーム画像Ｉ１の前後する２枚のフレーム画像における各共通座標画素の信号値の差分(または比)を算出し、得られた各差分(または各比)の中の最大値を抽出して演算信号値とする処理である。この処理により得られた解析画像からは、手軽に時間変化速度の最大値あるいは、時間変化加速度の最大値を把握することができる。例えば呼気時のフレームを特定解析フレームとするとこで、息の吐きやすさを評価し易くなるため、閉塞性疾患や拘束性疾患など、動きや加速度に変化が生じる疾患領域を検出しやすくなる。
【００３６】
最小微分値処理は、複数の特定解析フレーム画像Ｉ１の前後する２枚のフレーム画像における各共通座標画素の信号値の差分または比を算出し、得られた各差分または各比の中の最小値を抽出して演算信号値とする処理である。この処理により得られた解析画像からは、手軽に時間変化速度の最小値あるいは、時間変化加速度の最小値を把握することができる。例えば呼気時のフレームを特定解析フレームとするとこで、息の吐きやすさを評価し易くなるため、閉塞性疾患や拘束性疾患など、動きや加速度に変化が生じる疾患領域を検出しやすくなる。
【００３７】
なお、制御部１１は、一度に異なる複数種類の解析結果動画に演算処理を加える旨の指示を受ける（サマリー表示領域１５Ａにドラッグ及びドロップされる）と、解析結果動画ごとに演算処理を行うことが可能となっている。
【００３８】
また、制御部１１は、上記演算処理を施すことで生成した演算画像データに基づく演算結果画像Ｉ２を、自動で又はユーザーからの指示に基づいて表示部１５に表示することが可能となっている。すなわち、制御部１１は、本発明における画像表示手段をなす。
【００３９】
また、本実施形態の放射線画像処理装置１は、表示部１５の表示領域の少なくとも一部にサマリー表示領域１５Ａを形成することが可能となっている。サマリー表示領域１５Ａは、主に、上述した演算結果画像Ｉ２を表示するための領域である。
サマリー表示領域１５Ａは、常時表示されるものであってもよいし、必要に応じて表示・非表示が切り替えられるものであってもよい。また、サマリー表示領域１５Ａは、表示部１５の一部に表示されるものであっても、全体に表示されるものであってもよい。
なお、表示部１５とは異なる他の表示装置を備え、そこによりサマリー表示領域１５Ａを表示するようにすることも可能である。
【００４０】
サマリー表示領域１５Ａにおいては、演算結果画像Ｉ２等を様々な表示態様で表示することが可能となっている。本実施形態においては、単独表示態様、並列表示態様、重畳表示態様、切替表示態様の４種類の表示態様がある。
また、演算結果画像Ｉ２を、サマリー表示領域１５Ａだけではなく、表示部１５に表示するようにしてもよい。
【００４１】
単独表示態様は、例えば図２（ｃ），（ｄ）又は（ｅ）に示したように、一のサマリー表示領域１５Ａに一の演算結果画像Ｉ２を表示する態様である。
なお、単独表示態様のサマリー表示領域１５Ａを、同時に複数表示するようにしてもよい。
【００４２】
重畳表示態様は、図２（ｆ）や図４（ｃ）に示したように、演算結果画像Ｉ２を、元動画Ｉ０や、元フレーム画像Ｉ０、解析結果動画Ｉ１、解析フレーム画像Ｉ１、他の演算結果画像Ｉ２に重畳させて表示する表示態様である。因みに、図２（ｆ）は、演算結果画像Ｉ２を元動画又は元フレーム画像Ｉ０に重畳した例、図４（ｃ）は異なる解析処理を施した二種類の解析画像データにそれぞれ同じ演算処理（ＭＩＰ処理）を施して得られた二枚の演算結果画像Ｉ２同士を重畳した例である。
なお、３種類以上の画像を重畳表示することも可能である。
以下、重畳表示態様で表示された画像を重畳画像Ｉ３と称することがある。
【００４３】
並列表示態様は、図３や図４（ｂ），（ｄ）に示したように、演算結果画像Ｉ２を、元動画Ｉ０や、元フレーム画像Ｉ０、解析結果動画Ｉ１、解析フレーム画像Ｉ１、他の演算結果画像Ｉ２、重畳画像Ｉ３と並べて表示する表示態様である。因みに、図３は特定解析フレーム画像Ｉ１の選択枚数だけを変えて同一の演算処理（ＭＩＰ処理）を施した複数の演算結果画像Ｉ２を並べた例、図４（ｂ）は異なる解析処理を施した二種類の解析動画データにそれぞれ同じ演算処理（ＭＩＰ処理）を施して得られた二枚の演算結果画像Ｉ２を並べた例、図４（ｄ）は、異なる複数の演算結果画像Ｉ２、これらの重畳画像Ｉ３及び元動画又は元フレーム画像Ｉ０を並べた例である。
また、上述したように、異なる複数種類の演算処理を行った場合には、各演算結果を並列表示することが可能である。
【００４４】
切替表示態様は、図５に示したように、演算結果画像Ｉ２、元動画Ｉ０、元フレーム画像Ｉ０、解析結果動画Ｉ１、解析フレーム画像Ｉ１、他の演算結果画像Ｉ２及び重畳画像Ｉ３の中から一の画像を所定時間表示した後、当該画像と入れ替えで当該画像とは異なる画像を表示する表示態様である。因みに、図５はある演算結果画像Ｉ２、他の演算結果画像Ｉ２、これらの重畳画像Ｉ３、元動画又は元フレーム画像Ｉ０の順で変更した場合の例である。
【００４５】
なお、並列表示態様と切り替え表示態様を組み合わせることも可能である。すなわち、一のサマリー表示領域１５Ａに複数の画像を切り替えながら表示することも可能である。
また、上記表示態様の説明では、サマリー表示領域１５Ａに表示される画像には、演算結果画像Ｉ２が必ず含まれるものとしたが、含まない形で表示することも可能である。すなわち、元動画又は元フレーム画像Ｉ０と解析結果動画又は解析フレーム画像Ｉ１を重畳表示、並列表示、切替表示することも可能である。
【００４６】
上述したように、本実施形態に係る放射線画像処理装置１は、複数の元フレーム画像からなる元動画（放射線動画）の動画データを取得する画像取得手段と、画像取得手段が取得した動画データに所定の解析処理を加え、複数の解析フレーム画像からなる解析結果動画を生成する解析手段と、解析手段が生成した解析結果動画を構成する複数の解析フレーム画像の中から複数の特定解析フレーム画像を選択する選択手段と、選択手段によって選択された複数の特定解析フレーム画像における、共通する座標上にそれぞれ位置する各共通座標画素の信号値に基づいて一の演算信号値を導出する演算処理手段と、前記演算処理手段が座標ごとに生成した複数の演算信号値を用いて前記表示部に演算結果画像を表示する画像表示手段と、を備えたものとなっている。
【００４７】
サマリー表示領域１５Ａにおいて、上述したような動画又は静止画像の表示を行うことにより、ユーザーにとって関心の高い演算結果画像Ｉ２を直感的に得ることができる。その結果、ユーザー自身が診断に適切な表示方法を容易に探ることができ、読影精度が向上する。
また、表示部１５に表示される動画のどこを注力して読影するべきかユーザーが理解し易くなり、効率的に読影できるようになるとともに、読影洩れを防止することができる。その結果、診断の精度を向上させることができる。
また、任意のフレーム画像を選択可能としているので、アーチファクト（ノイズ）が混入していると推定されるフレーム画像を容易に取り除くことができ、解析結果の精度が向上する。
すなわち、動画から、ポイントとなる情報だけを容易に読影することができ、ユーザーの負担が低減する。
【００４８】
また、図３に示したように、フレーム画像の選択枚数だけを変えて同一の演算処理（ＭＩＰ処理）を施した複数の演算結果画像Ｉ２を並列表示すれば、演算処理の妥当性、すなわち、選択フレームの範囲を広くしすぎたのか、もう少し広げられそうか等をユーザーに直感的にすばやく把握させることができる。
また、図４に示したように、異なる種類の複数の動画（フレーム画像の一覧）からそれぞれ選択したフレーム画像に対して演算処理した結果を並列表示したり、複数の動画からそれぞれ選択したフレーム画像に対して演算処理した結果を重畳表示したりすることで、診断に有用な「Ｖ／Ｑ（換気血流比）画像」をユーザーの手で容易に作成することが可能となる。
【００４９】
＜第２実施形態＞
次に、本発明の第２実施形態に係る放射線画像処理装置について説明する。図６〜９は本実施形態に係る放射線画像処理装置が行う画像選択の概念図である。なお、図６には、肺野領域の画像を例示したが、本実施形態に係る放射線画像処理装置１は、肺野以外にも、心臓領域、肺動脈、大動脈又は肺気管支等、所定周期で同じ動作を繰り返す部位の動画であれば扱うことが可能である。
【００５０】
第１実施形態では、演算処理を施す対象となる特定解析フレーム画像Ｉ１をユーザーが選択するようにしていたが、本実施形態に係る放射線画像処理装置１Ａは、予め設定された選択パターンに基づいて特定解析フレーム画像Ｉ１を自動的に選択するようになっている。
このため、本実施形態に係る放射線画像処理装置１Ａは、実行する制御（記憶部１３Ａの記憶内容）が第１実施形態と異なっている。なお、記憶部１３Ａの記憶内容を除く他の構成は第１実施形態と同様である。
【００５１】
すなわち、放射線画像処理装置１Ａは、記憶部１３Ａに、自動選択の範囲を規定する選択パターンを複数格納しており、ユーザーが操作部１４に対して行った操作に基づいて、予め用意された複数の選択パターンの中から何れかの選択パターンを設定する。そして、ユーザーが特に注意して読影したい関心領域（例えば図６（ａ）で符号Ｒで示したような領域）を選択すると、設定した選択パターンに基づいて前記特定解析フレーム画像を自動的に選択するようになっている。
【００５２】
本実施形態に係る放射線画像処理装置１Ａにおいて設定可能な選択パターンには、大きく分けて（１）極大値（極小値）−極大値（極小値）、（２）極大値（極小値）−極小値（極大値）、（３）極大値（極小値）とその前後所定範囲、の３つがある。
なお、このうちのいずれかのパターンのみが予め設定されている（ユーザーによる設定はできない）ようにしてもよいし、このうちの２パターンの中からユーザーが設定できるものとしてもよい。更に、ここで例示した以外のパターンを含むようにしてもよい。
【００５３】
選択パターン（１）が選択されると、放射線画像処理装置１は、少なくとも肺野、心臓領域、肺動脈、大動脈、あるいは肺気管支の一部が含まれる関心領域Ｒ（Region Of Interest；ROI）内の複数の画素信号値の平均値の時間方向の変化を数値化する（図６（ｂ）に示したようにグラフ化する）。そして、図６（ｂ），（ｃ）に示したように、ある極大値（極小値）から、当該極大値（極小値）と隣り合う極大値（極小値）までの間のフレーム画像（一周期分の動作を表す動画；選択範囲ｒ１）、又はある極大値（極小値）から、当該極大値と隣り合わない極大値（極小値）までの間の解析フレーム画像（複数周期分の動作を表す動画；選択範囲ｒ２）を特定解析フレーム画像として自動選択する。
なお、ｒ１とｒ２のどちらの選択範囲とするかは、ユーザーが選択パターンを設定する際に併せて設定するよう構成してもよいし、制御部１１が自動的に判断するよう構成してもよい。
【００５４】
診断対象が図６（ａ）に示したような肺野である場合、この選択範囲ｒ１，ｒ２は、吸気位から再度吸気位となるまでの動きを表すことになる。一方、診断対象が心臓である場合には、心収縮期から再度心収縮期となるまでの動きを表すことになる。
肺機能を把握する上では、肺の一呼吸分の動作や心臓の一心拍分の動作等、少なくとも一周期分の情報を読影することが重要である。このため、ここで挙げたような解析フレーム画像の選択パターンを設定することで、呼吸周期や心拍周期に関係する肺機能（例えば、肺への空気の出入り量（１回換気量）、肺野内への血液の出入り量（心拍出量、１回拍出量））を正確に把握することができる。
特に、複数周期に亘る選択範囲ｒ２を設定することで、呼吸や心拍のバラつきを低減することが可能である。
【００５５】
選択パターン（２）が選択されると、放射線画像処理装置１は、関心領域Ｒ内の平均値の時間方向の変化を数値化（グラフ化）する。そして、図７（ａ）に示したように、ある極大値から、当該極大値の後（グラフでは右側）に隣り合う極小値までの間の解析フレーム画像（選択範囲ｒ３）、又は図７（ｂ）に示したように、ある極小値から、当該極小値の後（グラフでは右側）に隣り合う極大値までの間の解析フレーム画像（選択範囲ｒ４）を特定解析フレーム画像として自動選択する。
なお、図７に示したように、一のグラフに対して選択範囲ｒ３又はｒ４を複数周期分設定することが可能である。
なお、ｒ３とｒ４のどちらの選択範囲とするかは、ユーザーが解析フレーム画像の選択パターンを設定する際に併せて設定できるようにしてもよいし、制御部１１が自動的に判断するようにしてもよい。
【００５６】
診断対象が肺野である場合には、選択範囲ｒ３は呼気中、選択範囲は吸気中の動きを表すことになる。一方、診断対象が心臓である場合には、選択範囲ｒ３又はｒ４は心臓収縮中又は心臓膨張中の動きを表すことになる。
複数の呼気中（吸気中）あるいは心臓膨張中（心臓収縮中）の解析フレーム画像の結果に対して演算処理を加えることによって、呼吸や心拍のバラつきなどのノイズ影響を低減させることができ、血液が肺野に広がっていく心収縮中の様子、心臓に血液が戻っていく心膨張中の様子、吸気で肺が広がっていく様子、呼気で肺が縮んでいく様子等を精度よく把握することができる。
また、例えば閉塞性疾患患者の息の吐きにくさの評価、拘束性疾患や肥満の方の息の吸い難さの評価、心臓の動き方が不連続になる不整脈などの評価にも有用である。
【００５７】
選択パターン（３）が選択されると、放射線画像処理装置１は、関心領域Ｒ内内の平均値の時間方向の変化を数値化（グラフ化）する。そして、図８（ａ）に示したように、ある極大値となる時間ｔから所定時間ｔ１だけ前に撮影された解析フレーム画像から、当該極大値となる時間から所定時間ｔ２だけ後に撮影された解析フレーム画像まで（ｔ−ｔ１〜ｔ＋ｔ２：選択範囲ｒ５）、又は図８（ｂ）に示したように、ある極小値となる時間ｔから所定時間ｔ１だけ前に撮影された解析フレーム画像から、当該極小値となる時間から所定時間ｔ２だけ後に撮影された解析フレーム画像まで（ｔ−ｔ１〜ｔ＋ｔ２：選択範囲ｒ６）を特定解析フレーム画像として自動選択する。所定時間ｔ１，ｔ２は一周期（呼吸周期、心拍周期）の１／４以下とするのが好ましく、より好ましくは１／８以下とする。
【００５８】
なお、図８には、所定時間ｔ１とｔ２が同じ長さである場合を例示したが、ｔ１とｔ２は異なっていてもよい。
また、時間ではなく、極大からの解析フレーム画像数で範囲を規定するようにしてもよい。
更に、図８に示したように、一のグラフに対して選択範囲ｒ５又はｒ６を複数周期分設定することが可能である。
呼吸周期（心拍周期）は患者ごとに、スパイロメトリやモストグラフ、心電図等で計測した値を用いても良いし、各国、身長や性別、年齢などで定められた平均的な値を用いても良い。
【００５９】
診断対象が肺野である場合には、極大及びその周辺は吸気位辺り、極小及びその周辺は呼気位辺りの動きを表すことになる。一方、診断対象が心臓である場合には、心臓収縮期辺り又は心臓膨張期辺りの動きを表すことになる。
複数の呼気位（吸気位）又は心臓膨張期（心臓収縮期）付近の解析フレーム画像の結果に対して演算処理を加えることによって、血液が肺野に広がり切った心収縮期の様子、心臓に血液が戻り切った心膨張期の様子、呼吸により肺が広がりきった様子、肺が縮み切った様子等を精度よく把握することができる。
【００６０】
心臓から肺に血が広がる速度は速い為、肺野内のどこに血液が流れて行ったかを把握するのは通常の解析結果（動画）を見ているだけでは容易ではない。一方、全ての解析フレーム画像に対して演算処理を加えるとアーチファクトも多く付加されるデメリットがある。このため、血液が肺野に広がり切る前後の解析フレーム画像に対して演算処理を加えることで、アーチファクトの少ない肺血流分布画像を得ることが可能となる。つまり、選択パターン（３）によって得られた特定解析フレーム画像は、肺血流に関する疾患の把握に有用である。例えば、肺塞栓など局所的な肺血流機能欠損領域の把握や、心機能が低下する心不全や心肥大の評価に有用である。一方で、（３）によって得られた特定解析フレーム画像は同様な理由でアーチファクトの少ない換気機能画像となる為、局所的な肺換気機能欠損領域の把握や、横隔膜がどこまで広がっているか等、肺の広がり方の評価にも有用である。
【００６１】
なお、特定解析フレーム画像を複数周期分選択する場合には、サマリー表示領域１５Ａに各選択範囲の演算結果画像を並列表示するようにしてもよい。
こうすることで、どの呼吸（どの心拍）が最も正常/異常かを判断し易くなる。
【００６２】
なお、選択パターン（３）においては、所定時間ｔ１，ｔ２のいずれか一方を０としてもよい。すなわち、図９（ａ）に示したように、ある極大値となる時間ｔから所定時間ｔ２だけ後に撮影された解析フレーム画像まで（極大〜ｔ＋ｔ２：選択範囲ｒ７）、ある極小値となる時間ｔから所定時間ｔ２だけ後に撮影された解析フレーム画像まで（極小〜ｔ＋ｔ２：選択範囲ｒ８）、ある極大値となる時間ｔから所定時間ｔ１だけ前に撮影された解析フレーム画像から当該極大まで（ｔ−ｔ１〜極大：選択範囲ｒ９）、又はある極小値となる時間ｔから所定時間ｔ１だけ前に撮影された解析フレーム画像から当該極小まで（ｔ−ｔ１〜極小：選択範囲ｒ１０）を特定解析フレーム画像として自動選択するようにしてもよい。
【００６３】
なお、本実施形態では、関心領域Ｒ内の画素信号値の平均値の時間方向の変化を数値化（グラフ化）し、そこから特定解析フレーム画像を選択するようにしたが、他の統計値（最大値、最小値、最頻値、積算値、標準偏差等）を用いて特定解析フレーム画像を選択するようにしてもよい。
【００６４】
また、横隔膜移動量の時間変化や、肺野面積の時間変化、肋間距離の変化量の時間変化、肺野領域の水平方向の広がりの時間変化、心臓面積の時間変化、大動脈弓の位置の時間変化等の形態情報の時間方向変化を元に、特定解析フレーム画像を選択するようにしても良い。
撮影対象が、例えば肺野や心臓領域にて脂肪や筋肉などが不均一に分布している体型の患者（肥満の方や女性等に多い）であったり、撮影中に体が上下左右に大きく動いてしまう（体動が大きい）患者であったりする場合、画素値を用いた時間変化情報が正確に呼吸状態や心拍状態を反映していない場合がある。しかし、形態情報の時間方向の変化を用いることで、このような患者に対しても特定解析フレーム画像を精度よく選択することができるようになる。
【００６５】
＜第３実施形態＞
次に、本発明の第３実施形態に係る放射線画像処理装置１Ｂについて説明する。図１０は放射線画像処理装置１Ｂにおける画像処理の概念図であり、図１１，１２は処理結果の表示例である。
なお、図１０には、肺野領域の画像を例示したが、本実施形態に係る放射線画像処理装置１Ｂは、肺野以外の部位を撮影した動画についても処理・表示することが可能である。
【００６６】
第１，第２実施形態では、サマリー表示領域１５Ａに表示する内容を画像のみとしたが、本実施形態に係る放射線画像処理装置１Ｂは、所定の定量値をサマリー表示領域１５Ａに表示するようになっている。
このため、本実施形態に係る放射線画像処理装置１Ｂは、実行する制御（記憶部１３Ｂの記憶内容）が第１，第２実施形態と異なっている。なお、記憶部１３Ｂの記憶内容を除く他の構成は第１，第２実施形態と同様である。
【００６７】
放射線画像処理装置１Ｂの制御部１１は、複数の演算信号値を算出した後、演算結果画像の関心領域Ｒ内に位置することになる複数の画素の信号値に基づいて定量値Ｑを導出するようになっている。すなわち、制御部１１は、本発明における定量値算出手段をなす。なお、関心領域Ｒは、本発明における特定領域である。以下、関心領域Ｒ内に位置することになる画素を領域内画素と称する。なお、関心領域Rは複数設定されていても良く、その場合は、各々の関心領域Rに対応した定量値Ｑを各々個別に導出しても良いし、複数設定された各々の関心領域Rをまとめて一つとみなして、一つの定量値Ｑを導出しても良い。
定量値Ｑの導出の仕方としては、例えば、最大値処理、最小値処理、Ｐｋ−Ｐｋ値処理、標準偏差処理、平均値処理、最頻値処理、積算値処理等がある。
【００６８】
最大（小）値処理は、演算結果画像の関心領域Ｒ内に位置することになる複数の領域内画素の各信号値の中から最大（小）値を抽出して定量値Ｑとする処理である。この処理により得られた定量値Ｑを用いれば、経過観察（換気、血流の悪化・良化等）や、左右肺機能の比較等を行うことができる。特に、低い計算コストで、肺換気機能や肺血流機能、心機能に関する特徴量の最大（小）値を把握できるので、例えば、治療前後などの経過観察（肺活量の変化や血圧の変化を見る）に有用である。
【００６９】
最頻値処理は、複数の領域内画素の各信号値の中から最頻値を抽出して定量値Ｑとする処理である。この処理により得られた定量値Ｑを用いれば、経過観察や左右肺機能の比較等を行うことができる。特に、演算結果画像上に特異的に発生する画欠等のアーチファクトの影響を低減できる為、関心領域Ｒ内の肺換気機能や肺血流機能、心機能を精度よく把握することができる。
【００７０】
Ｐｋ−Ｐｋ値処理は、複数の領域内画素の各信号値の中から最大値と最小値を抽出してその差分を算出し、当該領域内画素の数に対する最大値と最小値との差の割合を算出して定量値Ｑとする処理である。この処理により得られた定量値Ｑを用いれば、経過観察や左右肺機能の比較等を行うことができる。特に、肺換気機能や肺血流機能、心機能に関する特徴量の変化量を把握できるので、関心領域Ｒ内における上記各機能の特徴量のバラつき（局所的に肺換気機能や肺血流機能が異常を示す疾患）を把握しやすくなる。
【００７１】
標準偏差処理は、複数の領域内画素の各信号値の標準偏差を算出して定量値Ｑとする処理である。この処理により得られた定量値Ｑを用いれば、経過観察や左右肺機能の比較等を行うことができる。特に、低い計算コストで、肺換気機能や肺血流機能、心機能に関する特徴量の変化量を把握できるので、関心領域Ｒ内の上記各機能のバラつきを把握しやすくなる。
【００７２】
平均値処理は、複数の領域内画素の各信号値の平均値を算出して定量値Ｑとする処理である。この処理により得られた定量値Ｑを用いれば、経過観察や左右肺機能比較等を行うことができる。特に、低い計算コストで、関心領域Ｒ内の肺換気機能や肺血流機能、心機能を把握することができる。
【００７３】
積算値処理は、複数の領域内画素の各信号値の積算値を算出して定量値Ｑとする処理である。この処理により得られた定量値Ｑを用いれば、経過観察や左右肺機能比較等を行うことができる。特に、関心領域Ｒの取り方にわずかな変化があっても算出される定量値Ｑに変化が少ない為、測定誤差が生じにくくなる。また、関心領域Ｒ内のトータルでの肺機能、心機能を把握しやすくなる。
【００７４】
なお、定量値Ｑは、関心領域Ｒ内の領域内画素の信号値に基づいて導出するのではなく、演算結果画像の全画素のうち、信号値が所定閾値以上、所定閾値以下又は所定範囲内にある画素の個数（面積）としてもよい。
このようにすれば、例えば、肺全体（全体肺面積）に対して、正常又は異常な肺領域が全体の何割あるか等の情報を容易に把握でき、肺の局所的な機能分布(どこに肺機能の低下、欠損があるか)を評価しやすくなる。
【００７５】
本実施形態において、関心領域Ｒとして設定できるのは、全画像領域、肺野内、心臓領域内、肺動脈内、大動脈内、肺気管支内である。
全画像領域を関心領域Ｒとした定量値Ｑを用いれば、計算コストをかけずに、手軽に大局的な肺換気機能、肺血流機能、心機能を把握することができる。
肺野内を関心領域Ｒとした定量値Ｑを用いれば、肺換気機能（呼吸に伴う濃度の変化、横隔膜の動き、肋骨の動き、鎖骨の動き、肩甲骨の動きなど）、肺血流機能（心拍に伴う濃度の変化、血管の動き）を把握することができる。
【００７６】
心臓領域内を関心領域Ｒとした定量値Ｑを用いれば、心機能（心壁の移動量（移動速度、加速度）、心拍に伴う濃度の変化量（濃度変化速度、濃度変化加速度））を把握することができる。
肺動脈内を関心領域Ｒとした定量値Ｑを用いれば、肺血流機能、心機能を把握することができる。
大動脈内を関心領域Ｒとした定量値Ｑを用いれば、体循環の血流機能（心拍に伴う濃度の変化、血管の動き）、心機能を把握することができる。
肺気管支内を関心領域Ｒとした定量値Ｑを用いれば、肺換気機能の把握、気道の肥厚の程度、気道収縮の程度、無気肺/胸水/腫瘤、のう胞/肺気腫の有無などを把握することができる。
【００７７】
なお、定量値Ｑの表示は、図１０（ａ）に示したように単独で行ってもよいし、図１０（ｂ），（ｃ）に示したように、画像と併せて行うようにしてもよい。
画像Ｉ２，Ｉ３と併せて表示する場合には、図１０（ｂ）に示したように並列表示する他、切替表示、重畳表示することも可能である。
また、図１０（ｃ）に示したように、画像Ｉ２，Ｉ３と関心領域Ｒを重畳表示し、更にその中に定量値Ｑを表示するようにしてもよい。このようにすれば、左右の肺機能を直感的に把握することが可能となる。
【００７８】
また、制御部１１は、種類の異なる複数のフレーム画像群からそれぞれフレーム画像を選択（サマリー表示領域１５Ａにドラッグ及びドロップ等）した場合は、各動画についての演算結果画像Ｉ２を生成するだけでなく、各動画についての定量値Ｑも算出する。
このとき、換気解析動画（解析結果動画）について算出した定量値Ｑと、血流解析動画（解析結果動画）について算出した定量値Ｑを用いた特徴量（比や差分など）を算出し、図１１に示したようにＶ／Ｑ比（換気血流比）を示すＶ／Ｑ画像Ｖを表示するようにしても良い。
このようにすれば、ユーザーが正常異常の判断を行いやすくなるだけでなく、経時的な比較が可能となる。
【００７９】
また、図１２（ａ）に示したように、左右の肺ごとに定量値Ｑを算出し、さらに各定量値Ｑの左右肺比率を算出して定量値Ｑとともに表示するようにしても良い。こうすることで、左右肺の肺機能の不均等がわかり易くなり、正常と異常の把握が行いやすくなる。
また、心臓領域を診断対象とする場合には、図１２（ｂ）に示したように、心室、心房、右心系、左心系ごとに定量値Ｑ及び比率を算出し、表示するようにしてもよい。こうすることで、心機能の不均衡がわかり易くなり、正常と異常の把握が行いやすくなる。
【００８０】
また、本実施形態において算出した定量値Ｑについては、演算結果画像Ｉ２を保存する際、画像のヘッダーに保存するようにしても良い。このようにすれば、演算結果画像Ｉ２を再度読み込む際、演算結果画像Ｉ２のファイル名一覧でヘッダー情報の有無を見ることにより、必要な場合のみ定量値Ｑを表示させることができる。
また、画像とは別データ（テキストやＣＳＶファイルなど）として保存しても良い。そうすることで、定量値Ｑのみを再度見直すことが容易になるし、不要に画像を毎回開く手間が減り、経過観察等の作業効率が向上する。
【００８１】
（変形例１）
なお、上記実施形態においては、サマリー表示領域１５Ａに選択した複数の解析フレーム画像に基づく演算結果画像や定量値Ｑを表示するようにしたが、図１３に示したように、単に選択（サマリー表示領域１５Ａにドラッグ及びドロップ）した複数のフレーム画像のみを表示するようにしてもよい。
なお、本変形例においても、上記実施形態と同様、選択するフレーム画像は必ずしも連続のフレーム画像である必要はない。また、サマリー表示領域１５Ａへの表示は、フレーム画像の選択後自動的に行ってもよいし、ユーザーの操作（例えばＯＫボタンをクリックする等）に基づいて行うようにしてもよい。
【００８２】
また、異なる複数種類のフレーム画像の一覧からそれぞれフレーム画像を選択した場合は、図１４（ａ）に示したように重畳画像Ｉ３の表示（ここでは、元フレーム画像と解析フレーム画像の重畳表示）の他、図１４（ｂ）に示した並列表示（ここでは、解析フレーム画像Ｉ２、元フレーム画像Ｉ０及び重畳画像Ｉ３の並列表示）、図１４（ｃ）に示した切替表示（ここでは、解析フレーム画像Ｉ２、元フレーム画像Ｉ０及び重畳画像Ｉ３の切替表示）を行う。
なお、異なる複数種類の動画の中からそれぞれフレーム画像を選択する際、それぞれのフレーム画像数は同一であっても異なっていても良い。
【００８３】
本変形例のようにすれば、ユーザーにとって関心の高いフレーム画像のみがサマリー表示領域１５Ａに表示される為、ユーザーは読影に注力しやすくなり、読影精度が向上する。
特に、１呼吸分など呼吸周期の整数倍に相当するフレーム画像のみをサマリー表示領域１５Ａに表示することで、呼吸が正しく行われているかを評価しやすくなる。
また、１心拍分など、心拍周期の整数倍に相当するフレーム画像のみをサマリー表示領域１５Ａに表示することで、心臓が正常に動いているかを評価しやすくなる。
また、吸気始めあるいは、呼気始めに相当するフレーム画像のみをサマリー表示領域１５Ａに表示することで、呼吸器疾患症例に対しての、息の吸い難さや息の吐きにくさを評価しやすくなる。
また、心膨張あるいは心収縮に相当するフレーム画像のみをサマリー表示領域１５Ａに表示することで、心疾患症例に対しての、血液の拍入のし易さ、血液の拍出のし易さ、不整脈の有無や心肥大の有無を評価しやすくなる。
また、肺機能が正常あるいは異常と判断したフレーム画像群をサマリー表示領域１５Ａに表示させたうえで、表示部１５に元々表示されているフレーム画像の一覧と比較することで、表示部１５に元々表示されているフレーム画像群が総じて正常なのか異常なのかを判断しやすくなる。すなわち、画像表示部に元々表示されているフレーム画像群の読影効率が向上する。
【００８４】
なお、図１５（ａ）に示したように、換気解析結果と、血流解析結果をサマリー表示領域１５Ａに重畳表示することで、容易に換気機能と血流機能を比較する画像（Ｖ／Ｑ画像）と同等の画像を得ることができる。
このようにすれば、ユーザーが選択したフレーム画像のみを用いて換気機能と血流機能の比較が行える為、例えばアーチファクトの影響をあまり受けていないＶ／Ｑ画像を得ることができ、診断精度が向上する。
また、図１５（ｂ）に示したように、換気解析結果と、血流解析結果、これらの重畳画像及び元画像の並列表示や、切替表示も可能である。
【００８５】
（変形例２）
また、上記実施形態では、表示部１５に表示されたフレーム画像の一覧の中から、サマリー表示領域１５Ａに解析結果動画等を表示するためのフレーム画像を選択するようにしていたが、図１６に示したように、表示部１５において動画を再生し、動画から必要な画像を直接ドラックしサマリー表示領域１５Ａにドロップするようにしても良い。
【００８６】
また、図１７に示したように、各フレーム画像にチェックボックスＢを設け、チェックＣを入れたフレーム画像がサマリー表示領域１５Ａに表示されるようにしても良い。このようにすれば、サマリー表示領域１５Ａに表示されたフレーム画像又は動画を見ながら、チェックＣを入れるフレーム画像を選べるため、ユーザーが選択したフレーム画像が妥当かどうか（診断にふさわしいフレーム画像かどうか）判断を素早く行うことが可能となる。
【００８７】
（変形例３）
また、図１８に示したように、選択されサマリー表示領域１５Ａに表示されるフレーム画像については、選択されていないフレーム画像や動画とは異なるルックアップテーブル（Lookup table；ＬＵＴ）を適用した画像ＩＬを表示してもよい。具体的には、濃度、コントラスト、カラーマップの変更し表示する。
このようにすれば、ユーザーが関心の高い画像が、ユーザーの関心の高い色づけでサマリー表示領域１５Ａに表示される為、ユーザーは読影に注力しやすくなり、読影精度が向上する。
【符号の説明】
【００８８】
１，１Ａ，１Ｂ  放射線画像処理装置
１１   制御部（画像取得手段、解析手段、選択手段、演算処理手段、画像表示手段）
１２   通信部
１３，１３Ａ，１３Ｂ  記憶部
１４   操作部
１５   表示部
１５Ａ サマリー表示領域
１６   バス
Ｉ０  元動画、元フレーム画像
Ｉ１  解析結果動画、解析フレーム画像
Ｉ２  演算結果画像
Ｉ３  重畳画像
Ｒ  関心領域
ｒ１〜ｒ１０  （フレーム画像の）選択範囲
ｔ  時間
ｔ１，ｔ２  所定時間
【図１】
2018157968.tif 000003
【図２】
2018157968.tif 000004
【図３】
2018157968.tif 000005
【図４】
2018157968.tif 000006
【図５】
2018157968.tif 000007
【図６】
2018157968.tif 000008
【図７】
2018157968.tif 000009
【図８】
2018157968.tif 000010
【図９】
2018157968.tif 000011
【図１０】
2018157968.tif 000012
【図１１】
2018157968.tif 000013
【図１２】
2018157968.tif 000014
【図１３】
2018157968.tif 000015
【図１４】
2018157968.tif 000016
【図１５】
2018157968.tif 000017
【図１６】
2018157968.tif 000018
【図１７】
2018157968.tif 000019
【図１８】
2018157968.tif 000020