特許 6805918 放射線画像処理装置及び放射線画像撮影システム

(19)【発行国】日本国特許庁(JP)

(12)【公報種別】特許公報(B2)

(11)【特許番号】6805918

(24)【登録日】2020年12月8日

(45)【発行日】2020年12月23日

(54)【発明の名称】放射線画像処理装置及び放射線画像撮影システム

(51)【国際特許分類】

   A61B 6/00 20060101AFI20201214BHJP

【ＦＩ】

   A61B6/00 330A

   A61B6/00 350A

【請求項の数】11

【全頁数】23

(21)【出願番号】特願2017-56991(P2017-56991)

(22)【出願日】2017年3月23日

(65)【公開番号】特開2018-157968(P2018-157968A)

(43)【公開日】2018年10月11日

【審査請求日】2019年6月20日

(73)【特許権者】

【識別番号】000001270

【氏名又は名称】コニカミノルタ株式会社

(74)【代理人】

【識別番号】110001254

【氏名又は名称】特許業務法人光陽国際特許事務所

(72)【発明者】

【氏名】金子 志行

(72)【発明者】

【氏名】板屋 敬子

(72)【発明者】

【氏名】勝原 慎介

【審査官】 亀澤 智博

(56)【参考文献】

【文献】 国際公開第２０１４／１９２５０４（ＷＯ，Ａ１）

【文献】 特開２０１６−２１４７２５（ＪＰ，Ａ）

【文献】 特開２０１７−０１８６８１（ＪＰ，Ａ）

【文献】 特開２０１５−０４３８８０（ＪＰ，Ａ）

(58)【調査した分野】（Int.Cl.，ＤＢ名）

Ａ６１Ｂ ６／００ − ６／１４

Ｇ０６Ｔ １／００

(57)【特許請求の範囲】

【請求項1】

画像を表示可能な表示部と、
複数のフレーム画像からなる放射線動画の動画データを取得する画像取得手段と、
前記画像取得手段が取得した動画データに所定の解析処理を加え、複数の解析フレーム画像からなる解析結果動画を生成する解析手段と、
前記解析手段が生成した解析結果動画を構成する複数の解析フレーム画像の中から複数の特定解析フレーム画像を選択する選択手段と、
前記選択手段によって選択された複数の特定解析フレーム画像における、共通する座標上にそれぞれ位置する各共通座標画素の信号値に基づいて一の演算信号値を導出する演算処理手段と、
前記演算処理手段が座標ごとに生成した複数の演算信号値を用いて前記表示部に演算結果画像を表示する画像表示手段と、を備え、
前記演算処理手段は、前記複数の演算信号値を算出した後、前記演算結果画像の特定領域内に位置することになる複数の領域内画素の信号値に基づいて所定の定量値を導出する定量値算出手段を備え、
前記画像表示手段は、前記定量値算出手段が算出した定量値を前記表示部に表示させることが可能であることを特徴とする放射線画像処理装置。

【請求項2】

前記画像表示手段は、
前記表示部の表示領域の少なくとも一部にサマリー表示領域を形成し、
前記サマリー表示領域に前記演算結果画像を表示させることが可能であることを特徴とする請求項１に記載の放射線画像処理装置。

【請求項3】

ユーザーが操作可能な操作部を有し、
前記選択手段は、ユーザーが前記操作部に対して行った操作に基づいて前記特定解析フレーム画像を選択することを特徴とする請求項１又は請求項２に記載の放射線画像処理装置。

【請求項4】

前記選択手段は、予め設定された選択パターンに基づいて前記特定解析フレーム画像を自動的に選択することを特徴とする請求項１又は請求項２に記載の放射線画像処理装置。

【請求項5】

ユーザーが操作可能な操作部を有し、
前記選択手段は、
ユーザーが前記操作部に対して行った操作に基づいて、予め用意された複数の選択パターンの中から何れかの選択パターンを設定し、
設定した選択パターンに基づいて前記特定解析フレーム画像を自動的に選択することを特徴とする請求項４に記載の放射線画像処理装置。

【請求項6】

前記演算処理手段は、前記複数の特定解析フレーム画像における、前記各共通座標画素の信号値の中から最大値、最小値もしくは最頻値を抽出する、当該共通座標画素の信号値における、選択した特定解析フレーム画像の数に対する最大値と最小値との差の割合、標準偏差、平均値もしくは積算値を算出する、又は前後する２枚の特定解析フレーム画像における各共通座標画素の信号値の差分を算出し、得られた各差分の中の最大値もしくは最小値を抽出することにより、前記一の演算信号値を導出することを特徴とする請求項１から請求項５のいずれか一項に記載の放射線画像処理装置。

【請求項7】

前記定量値算出手段は、前記複数の領域内画素の信号値の中から最大値、最小値もしくは最頻値を抽出する、又は当該複数の領域内画素の信号値における、当該領域内画素の数に対する最大値と最小値との差の割合、標準偏差、平均値もしくは積算値を算出することにより、前記定量値を導出することを特徴とする請求項１から請求項６のいずれか一項に記載の放射線画像処理装置。

【請求項8】

前記特定領域は、肺野内、心臓領域内、肺動脈内、大動脈内又は肺気管支内であることを特徴とする請求項７に記載の放射線画像処理装置。

【請求項9】

画像を表示可能な表示部と、
ユーザーが操作可能な操作部と、
複数のフレーム画像からなる放射線動画の動画データを取得する画像取得手段と、
前記画像取得手段が取得した動画データに所定の解析処理を加え、複数の解析フレーム画像からなる解析結果動画を生成する解析手段と、
前記解析手段が生成した解析結果動画を構成する複数の解析フレーム画像の中から複数の特定解析フレーム画像を選択する選択手段と、
前記選択手段によって選択された複数の特定解析フレーム画像における、共通する座標上にそれぞれ位置する各共通座標画素の信号値に基づいて一の演算信号値を導出する演算処理手段と、
前記演算処理手段が座標ごとに生成した複数の演算信号値を用いて前記表示部に演算結果画像を表示する画像表示手段と、を備え、
前記選択手段は、ユーザーにより前記操作部になされた操作に応じた解析フレーム画像を前記特定解析フレーム画像として選択することを特徴とする放射線画像処理装置。

【請求項10】

前記画像表示手段は、
前記表示部の表示領域の少なくとも一部にサマリー表示領域を形成し、
前記サマリー表示領域に前記演算結果画像を表示させることが可能であることを特徴とする請求項９に記載の放射線画像処理装置。

【請求項11】

放射線を照射する放射線照射装置と、
前記放射線照射装置から放射線を受けることで放射線画像の画像データを生成する放射線画像撮影装置と、を備え、
前記放射線照射装置による放射線の照射と、放射線画像撮影装置による画像データの読出しを連続的に繰り返すことにより、放射線動画の動画データを生成することが可能な放射線画像撮影システムであって、
前記放射線画像撮影装置と直接又は間接的に通信可能に接続された請求項１から請求項１０のいずれか一項に記載の放射線画像処理装置と、を備えることを特徴とする放射線画像撮影システム。

【発明の詳細な説明】

【技術分野】

【0001】

本発明は、放射線画像処理装置、及びこの放射線画像処理装置を備えた放射線画像撮影システムに関する。

【背景技術】

【0002】

医療用のモダリティー（例えば、ＣＲ（Computed Radiography）やＣＴ（Computed Tomography）等）を用いて撮影した動画は、表示装置に表示されることで医者の診断に供される。その際、撮影された動画をそのまま表示するだけでなく、様々な解析処理を施してから表示することがある。
従来、動画を構成する複数のフレーム画像の一覧を表示部に表示し、解析に用いるフレーム画像を一覧の中から一括選択できるようにする技術（特許文献１参照）や、画像を即時表示、参照表示あるいは転送する際に、それぞれに最適な画像処理を行う技術（特許文献２参照）が知られている。

【先行技術文献】

【特許文献】

【0003】

【特許文献1】特開２０１２−１１０３９７号公報

【特許文献2】特開２０１２−１１０４６６号公報

【発明の概要】

【発明が解決しようとする課題】

【0004】

しかしながら、特許文献１，２に記載されたような従来の技術では、ユーザーが元の動画と解析を加えた動画、あるいは異なる種類の解析を施した動画同士を見比べながら読影することを考えたとき、各動画のそれぞれの関心のある領域を交互に視点を変えながら読影する必要が出てくる。各動画は並行して再生が進むので、一方の動画で詳しく読影したい領域を見つけ、他方の動画の対応する領域を読影するのに視点を変えた時には、本来見たいフレーム画像よりも後のフレーム画像が表示されることになってしまう。このため、一方の動画から他方の動画へ移る際に、フレーム画像を戻して読影しなければならない場合があった。
つまり、動画は情報量が多い上に時間の経過に伴って変化するため、読影（診断）が容易でなく、ユーザーに過度の負担がかかっていた。

【0005】

本発明は、上記の問題点に鑑みてなされたものであり、医療用のモダリティーで撮影された動画から、ポイントとなる情報だけを容易に読影できるようにすることを課題とする。

【課題を解決するための手段】

【0006】

前記の問題を解決するために、本発明は、
画像を表示可能な表示部と、
複数のフレーム画像からなる放射線動画の動画データを取得する画像取得手段と、
前記画像取得手段が取得した動画データに所定の解析処理を加え、複数の解析フレーム画像からなる解析結果動画を生成する解析手段と、
前記解析手段が生成した解析結果動画を構成する複数の解析フレーム画像の中から複数の特定解析フレーム画像を選択する選択手段と、
前記選択手段によって選択された複数の特定解析フレーム画像における、共通する座標上にそれぞれ位置する各共通座標画素の信号値に基づいて一の演算信号値を導出する演算処理手段と、
前記演算処理手段が座標ごとに生成した複数の演算信号値を用いて前記表示部に演算結果画像を表示する画像表示手段と、を備え、
前記演算処理手段は、前記複数の演算信号値を算出した後、前記演算結果画像の特定領域内に位置することになる複数の領域内画素の信号値に基づいて所定の定量値を導出する定量値算出手段を備え、
前記画像表示手段は、前記定量値算出手段が算出した定量値を前記表示部に表示させることが可能であることを特徴とする。

【発明の効果】

【0007】

本発明によれば、医療用のモダリティーで撮影された動画から、ポイントとなる情報だけを容易に読影できるようにすることができる。

【図面の簡単な説明】

【0008】

【図1】本発明の第１〜第３実施形態に係る放射線画像撮影システムの全体構成を示すブロック図である。

【図2】第１実施形態に係る放射線画像処理装置が行う画像処理の概念図である。

【図3】同実施形態に係る放射線画像処理装置が表示する画像処理結果の一例である。

【図4】同実施形態に係る放射線画像処理装置が行う画像処理の概念図である。

【図5】同実施形態に係る放射線画像処理装置が表示する画像処理結果の一例である。

【図6】第２実施形態に係る放射線画像処理装置が行う画像選択の概念図である。

【図7】同実施形態に係る放射線画像処理装置が行う画像選択の概念図である。

【図8】同実施形態に係る放射線画像処理装置が行う画像選択の概念図である。

【図9】同実施形態に係る放射線画像処理装置が行う画像選択の概念図である。

【図10】本発明の第３実施形態に係る放射線画像処理装置が行う画像選択の概念図である。

【図11】同実施形態に係る放射線画像処理装置が表示する画像処理結果の一例である。

【図12】同実施形態に係る放射線画像処理装置が表示する画像処理結果の一例である。

【図13】第１〜第３実施形態の変形例に係る放射線画像処理装置が行う画像処理の概念図である。

【図14】同変形例に係る放射線画像処理装置が表示する画像処理結果の一例である。

【図15】同変形例に係る放射線画像処理装置が行う画像処理の概念図である。

【図16】同変形例に係る放射線画像処理装置が行う画像処理の概念図である。

【図17】同変形例に係る放射線画像処理装置が行う画像処理の概念図である。

【図18】同変形例に係る放射線画像処理装置が表示する画像処理結果の一例である。

【発明を実施するための形態】

【0009】

以下、図面を参照して本発明の実施の形態を詳細に説明する。ただし、発明の範囲は、図示例に限定されない。

【0010】

＜第１実施形態＞
〔放射線画像撮影システムの構成〕
まず、本発明の第１実施形態に係る放射線画像撮影システム１００の構成について説明する。図１は本実施形態に係る放射線画像撮影システム１００の機能的構成を表すブロック図である。

【0011】

本実施形態の放射線画像撮影システム１００は、放射線画像処理装置１や、医療用のモダリティー２等で構成されている。
また、放射線画像撮影システム１００には、必要に応じて、図示しないコンソールやサーバー（例えば、医療用画像管理システム（Picture Archiving and Communication System：ＰＡＣＳ）を構成するもの）等が接続される。
放射線画像撮影システム１００を構成する各装置は、ＤＩＣＯＭ（Digital Image and Communications in Medicine）規格に準じており、各装置間の通信は、ＤＩＣＯＭに則って行われる。

【0012】

モダリティー２は。放射線照射装置２１と、放射線画像撮影装置２２と、を備えている。
放射線照射装置２１は、図示を省略するが、放射線を生成可能な回転陽極や回転陽極に電子ビームを照射するフィラメント等を有する放射線源や、設定された管電圧や管電流、照射時間（ｍＡｓ値）等に応じた線量の放射線を放射線源から照射させるジェネレーター等を備えている。
そして、ユーザーの操作に基づいて、放射線画像撮影装置２２へ放射線を照射するようになっている。

【0013】

放射線画像撮影装置２２は、図示を省略するが、放射線を受けることで線量に応じた電荷を蓄積する複数の放射線検出素子が二次元状（マトリクス状）に配列された基板や、各放射線検出素子に蓄積された電荷を画像データとして読み出す読み出し回路、外部装置と通信したり画像データを送信したりするための通信部等を備えている。
そして、放射線画像撮影装置２２は、放射線照射装置２から放射線の照射を受け、画像データを読み出すと、その画像データを、通信部を介して直ちに外部へ送信するようになっている。

【0014】

このような放射線照射装置２及び放射線画像撮影装置２２を備えたモダリティー２は、放射線照射装置２１による放射線の照射と、放射線画像撮影装置２２による電荷の蓄積・画像データの読み出しと、を連続的に繰り返すことで、複数のフレーム画像からなる放射線動画の動画データを生成することが可能となっている。
なお、放射線照射装置２１と放射線画像撮影装置２２は、ＣＴのように一体化されたものでもよい。

【0015】

〔放射線画像撮影システムの構成〕
次に、上記放射線画像撮影システム１００を構成する放射線画像処理装置１の構成について説明する。
放射線画像処理装置１は、コンピューター又は専用の制御装置として構成されており、図１に示したように、制御部１１、通信部１２、記憶部１３、操作部１４、表示部１５、を備えて構成され、各部はバス１６により接続されている。

【0016】

制御部１１は、ＣＰＵ、ＲＡＭ等により構成される。制御部１１のＣＰＵは、操作部１４の操作に応じて、記憶部１３に記憶されているシステムプログラムや、各種処理プログラムを読み出してＲＡＭ内に展開し、展開されたプログラムに従って、後述する画像解析処理を始めとする各種処理を実行したり、表示部１５の表示内容を制御したりする等、放射線画像処理装置１各部の動作を集中制御する。

【0017】

通信部１２は、ＬＡＮアダプタやモデムやＴＡ等を備え、通信ネットワークＮＴに接続された放射線画像撮影装置２２やコンソール、サーバー等をはじめとする各装置との間のデータ送受信を制御する。

【0018】

記憶部１３は、不揮発性の半導体メモリやハードディスク等により構成される。記憶部１３は、制御部１１で各種処理（例えば、後述する画像解析処理や画像演算処理）を実行するためのプログラム（画像解析処理プログラムや画像演算処理プログラム等）や、プログラムにより処理の実行に必要なパラメータ、処理結果、あるいは、各種処理を施す前の元動画や静止画、各種処理を施した後の動画や静止画のデータを記憶する。これらの各種プログラムは、読取可能なプログラムコードの形態で格納されている。

【0019】

操作部１４は、カーソルキー、数字入力キー、及び各種機能キー等を備えたキーボードや、マウス等のポインティングデバイス等によりユーザーが操作可能に構成されており、キーボードに対するキー操作あるいはマウス操作により入力された指示信号を制御部１１に出力する。
なお、操作部１４は、表示部１５の表示画面に備えられたタッチパネルで構成されていても良く、この場合、タッチパネルを介して入力された指示信号を制御部１１に出力する。

【0020】

表示部１５は、ＬＣＤやＣＲＴ等のモニタにより構成され、制御部１１から入力される表示信号の指示に従って、各種画像や、操作部１４からの入力指示、データ等を表示することが可能となっている。

【0021】

〔放射線画像処理装置の動作〕
次に、本実施形態に係る放射線画像処理装置１の動作について説明する。図２，４は放射線画像処理装置１における画像処理の概念図であり、図３，５は処理結果の表示例である。
なお、本実施形態に係る放射線画像処理装置１は、静止画像の処理や表示も可能であるが、ここでは、主に動画を処理・表示する場合について説明する。また、各図には、肺野領域の画像を例示したが、本実施形態に係る放射線画像処理装置１は、肺野以外の部位を撮影した動画についても処理・表示することが可能である。

【0022】

放射線画像処理装置１の制御部１１は、操作部１４を介してユーザーから処理を施す前の元の放射線動画を表示する旨の指示を受けると、複数のフレーム画像からなる放射線動画の元動画データ（ローデータ等）の取得（通信部１２を介してモダリティー２、コンソールもしくはサーバーから受信、又は記憶部１３から読み込み）を行い、準備した元動画データに基づく元動画又は図２（ａ）に示したような元動画を構成する複数の元フレーム画像Ｉ_０の一覧（画像群））を表示部１５に表示することが可能となっている。すなわち、制御部１１は、本発明における画像取得手段をなす。なお、図２（ａ）では、元フレーム画像を符号Ｉ_０で示したが、以下の説明では、元動画であるか元フレーム画像であるかにかかわらず、元動画データに基づいて表示部１５に表示されるものをＩ_０で示すこととする。

【0023】

また、制御部１１は、操作部１４を介してユーザーから元動画データに動態解析を加える旨の指示を受けると、元動画データ（ローデータ等）の取得を行い、記憶部１３に記憶されている画像解析処理プログラムとの協働により元動画データに各種解析処理（換気解析や血流解析等）を施し、その解析結果である解析動画データを生成することが可能となっている。
そして、制御部１１は、生成した解析動画データに基づく解析結果動画又は図２（ｂ）に示したような解析結果動画を構成する複数の解析フレーム画像Ｉ_１の一覧）を、自動で又はユーザーからの指示に基づいて表示部１５に表示することが可能となっている。すなわち、制御部１１は、本発明における解析手段をなす。なお、図２（ｂ）では、解析フレーム画像を符号Ｉ_１で示したが、以下の説明では、解析結果動画であるか解析フレーム画像であるかにかかわらず、解析動画データに基づいて表示部１５に表示されるものをＩ_１で示すこととする。

【0024】

また、制御部１１は、操作部１４を介してユーザーから解析動画データに演算処理を加える旨の指示を受けると、解析動画データの取得を行い、記憶部１３に記憶されている画像演算処理プログラムとの協働により解析動画データに各種演算処理を施し、その演算結果である演算画像データを生成することが可能となっている。
なお、元動画データに演算処理を施して演算画像データを生成することも可能である。

【0025】

演算処理は、表示部１５に表示された、解析結果動画を構成する複数の解析フレーム画像Ｉ_１の一覧の中から複数の解析フレーム画像Ｉ_１（解析結果動画を構成する全ての解析フレーム画像Ｉ_１であってもよいし、その一部であってもよい。）が選択されることにより開始される。以下、選択された解析結果フレーム画像Ｉ_１を特定解析フレーム画像Ｉ_１と称する。
本実施形態において、特定解析フレーム画像Ｉ_１の選択は、ユーザーが操作部１４に対して行った操作に基づいて行われるようになっている。具体的には、例えば図２（ｂ）に示したように、表示部１５に表示された解析結果動画を構成する解析フレーム画像Ｉ_１の一覧の中から、ユーザーが処理対象とする特定解析フレーム画像Ｉ_１を所定領域（例えば後述するサマリー表示領域）へドラッグ及びドロップしたり、解析フレーム画像Ｉ_１の端部に設けられたチェックボックスＢにチェックＣを入れたりすることで選択される。すなわち、制御部１１は、本発明における選択手段をなす。
なお、選択する解析フレーム画像Ｉ_１は連続している必要は無く、間を飛ばして任意の解析フレーム画像Ｉ_１を選択することも可能である。

【0026】

演算処理では、選択された複数の特定解析フレーム画像Ｉ_１における共通する座標上にそれぞれ位置する各画素の信号値に基づいて一の演算信号値を導出する。そして、この処理を各座標（（１，１）から（ｍ、ｎ）まで）について繰り返し行い、座標ごとに生成した複数の演算信号値を用いて演算結果画像Ｉ_２となる演算画像データを生成する。すなわち、制御部１１は、本発明における演算処理手段をなす。以下、各特定解析フレーム画像Ｉ_１における共通する座標上にそれぞれ位置する各画素のことを共通座標画素と称する。

【0027】

演算信号値の導出の仕方としては、例えば、最大値（ＭＩＰ）処理、最小値（ＭｉｎＰ）処理、Ｐｋ−Ｐｋ値処理、標準偏差処理、最頻値処理、積算値処理、平均値処理、最大微分値処理、最小微分値処理等がある。制御部１１は、自動で又はユーザーからの指示に基づいて、上記の処理の中から何れかを設定して演算信号値を導出する。

【0028】

最大値処理は、複数の特定解析フレーム画像Ｉ_１における、各共通座標画素の信号値の中から最大値を抽出して演算信号値とする処理である。この処理により得られた解析画像からは、解析結果の値（換気や血流）の最大値が分かるので、換気や血流の有無を一目で把握することができる。具体的には、無気肺や肺塞栓領域において、換気機能や血流機能がどこまで回復しているのかを把握しやすくなる。また、処理コストを低く抑えることができる。

【0029】

最小値処理は、複数の特定解析フレーム画像Ｉ_１における、各共通座標画素の信号値の中から最小値を抽出して演算信号値とする処理である。この処理により得られた解析画像からは、解析結果の値（換気や血流）の最小値が分かるので、換気機能や血流機能がどこまで低下しているのかを把握しやすくなる。また、処理コストを低く抑えることができる。

【0030】

最頻値処理は、複数の特定解析フレーム画像Ｉ_１における、各共通座標画素の信号値の中から最頻値を抽出して演算信号値とする処理である。この処理により得られた解析画像からは、解析結果の値（換気や血流）の最頻値が分かるので、特異的なノイズまたはアーチファクトの影響を低減させた上で、換気や血流の有無を一目で把握できるようになる。例えば複数呼吸分、複数心拍分の解析結果に対して最頻値処理を実施することで、呼吸や心拍のばらつき影響を低減させることができる。

【0031】

Ｐｋ−Ｐｋ値処理は、複数の特定解析フレーム画像Ｉ_１における、各共通座標画素の信号値の中から最大値と最小値を抽出してその差分または比を算出し、特定解析フレーム画像の数に対する最大値と最小値との差の割合を算出して演算信号値とする処理である。この処理により得られた解析画像からは、解析結果の値（換気や血流）の変動量が分かるので、換気や血流の有無を一目で把握できるようになる。例えば、無気肺や肺塞栓領域の把握を行いやすくなる（信号の変動に強い、ロバストな結果を出すことができる）。

【0032】

標準偏差処理は、複数の特定解析フレーム画像Ｉ_１における、各共通座標画素の信号値の標準偏差を算出して演算信号値とする処理である。この処理により得られた解析画像からは、解析結果の値（換気や血流）の変動量が（複雑な計算アルゴリズム不要で）分かるので、換気や血流の有無を一目で把握できるようになる。例えば、無気肺や肺塞栓領域の把握を行いやすくなる。また、処理コストを低く抑えることができる。

【0033】

積算値処理は、複数の特定解析フレーム画像Ｉ_１における、各共通座標画素の信号値の積算値を算出して演算信号値とする処理である。この処理により得られた解析画像からは、解析結果の値（換気や血流）のある一定時間分の積算結果が分かるので、ノイズまたはアーチファクトの影響を低減させることができ、換気や血流の有無の把握が行いやすくなる。また、複数呼吸分、複数心拍分の解析結果を積算することで、呼吸や心拍のばらつき影響を低減させることができる。また、ある一定時間の間に計測した値を積算した値、または画像を最終結果として出力する他モダリティー(例えばシンチグラフィ)と、類似した物理量を導出できる為、積算値処理を行うことで上記のような他モダリティーの情報と比較が行いやすくなる。その結果、診断の精度が向上する。更に、画素値は線量と関係が有る場合がある為、総照射線量（総被曝量）の推定が行いやすくなる為、被検者の被ばく管理などに有用である。

【0034】

平均値処理は、複数の特定解析フレーム画像Ｉ_１における、各共通座標画素の信号値の平均値を算出して演算信号値とする処理である。この処理により得られた解析画像からは、ノイズまたはアーチファクトの影響を低減させた上で、換気や血流の有無を一目で把握できるようになる。例えば複数呼吸分、複数心拍分の解析結果に対して平均値処理を実施することで、呼吸や心拍のばらつき影響を低減させることができる。

【0035】

最大微分値処理は、複数の特定解析フレーム画像Ｉ_１の前後する２枚のフレーム画像における各共通座標画素の信号値の差分(または比)を算出し、得られた各差分(または各比)の中の最大値を抽出して演算信号値とする処理である。この処理により得られた解析画像からは、手軽に時間変化速度の最大値あるいは、時間変化加速度の最大値を把握することができる。例えば呼気時のフレームを特定解析フレームとするとこで、息の吐きやすさを評価し易くなるため、閉塞性疾患や拘束性疾患など、動きや加速度に変化が生じる疾患領域を検出しやすくなる。

【0036】

最小微分値処理は、複数の特定解析フレーム画像Ｉ_１の前後する２枚のフレーム画像における各共通座標画素の信号値の差分または比を算出し、得られた各差分または各比の中の最小値を抽出して演算信号値とする処理である。この処理により得られた解析画像からは、手軽に時間変化速度の最小値あるいは、時間変化加速度の最小値を把握することができる。例えば呼気時のフレームを特定解析フレームとするとこで、息の吐きやすさを評価し易くなるため、閉塞性疾患や拘束性疾患など、動きや加速度に変化が生じる疾患領域を検出しやすくなる。

【0037】

なお、制御部１１は、一度に異なる複数種類の解析結果動画に演算処理を加える旨の指示を受ける（サマリー表示領域１５Ａにドラッグ及びドロップされる）と、解析結果動画ごとに演算処理を行うことが可能となっている。

【0038】

また、制御部１１は、上記演算処理を施すことで生成した演算画像データに基づく演算結果画像Ｉ_２を、自動で又はユーザーからの指示に基づいて表示部１５に表示することが可能となっている。すなわち、制御部１１は、本発明における画像表示手段をなす。

【0039】

また、本実施形態の放射線画像処理装置１は、表示部１５の表示領域の少なくとも一部にサマリー表示領域１５Ａを形成することが可能となっている。サマリー表示領域１５Ａは、主に、上述した演算結果画像Ｉ_２を表示するための領域である。
サマリー表示領域１５Ａは、常時表示されるものであってもよいし、必要に応じて表示・非表示が切り替えられるものであってもよい。また、サマリー表示領域１５Ａは、表示部１５の一部に表示されるものであっても、全体に表示されるものであってもよい。
なお、表示部１５とは異なる他の表示装置を備え、そこによりサマリー表示領域１５Ａを表示するようにすることも可能である。

【0040】

サマリー表示領域１５Ａにおいては、演算結果画像Ｉ_２等を様々な表示態様で表示することが可能となっている。本実施形態においては、単独表示態様、並列表示態様、重畳表示態様、切替表示態様の４種類の表示態様がある。
また、演算結果画像Ｉ_２を、サマリー表示領域１５Ａだけではなく、表示部１５に表示するようにしてもよい。

【0041】

単独表示態様は、例えば図２（ｃ），（ｄ）又は（ｅ）に示したように、一のサマリー表示領域１５Ａに一の演算結果画像Ｉ_２を表示する態様である。
なお、単独表示態様のサマリー表示領域１５Ａを、同時に複数表示するようにしてもよい。

【0042】

重畳表示態様は、図２（ｆ）や図４（ｃ）に示したように、演算結果画像Ｉ_２を、元動画Ｉ_０や、元フレーム画像Ｉ_０、解析結果動画Ｉ_１、解析フレーム画像Ｉ_１、他の演算結果画像Ｉ_２に重畳させて表示する表示態様である。因みに、図２（ｆ）は、演算結果画像Ｉ_２を元動画又は元フレーム画像Ｉ_０に重畳した例、図４（ｃ）は異なる解析処理を施した二種類の解析画像データにそれぞれ同じ演算処理（ＭＩＰ処理）を施して得られた二枚の演算結果画像Ｉ_２同士を重畳した例である。
なお、３種類以上の画像を重畳表示することも可能である。
以下、重畳表示態様で表示された画像を重畳画像Ｉ_３と称することがある。

【0043】

並列表示態様は、図３や図４（ｂ），（ｄ）に示したように、演算結果画像Ｉ_２を、元動画Ｉ_０や、元フレーム画像Ｉ_０、解析結果動画Ｉ_１、解析フレーム画像Ｉ_１、他の演算結果画像Ｉ_２、重畳画像Ｉ_３と並べて表示する表示態様である。因みに、図３は特定解析フレーム画像Ｉ_１の選択枚数だけを変えて同一の演算処理（ＭＩＰ処理）を施した複数の演算結果画像Ｉ_２を並べた例、図４（ｂ）は異なる解析処理を施した二種類の解析動画データにそれぞれ同じ演算処理（ＭＩＰ処理）を施して得られた二枚の演算結果画像Ｉ_２を並べた例、図４（ｄ）は、異なる複数の演算結果画像Ｉ_２、これらの重畳画像Ｉ_３及び元動画又は元フレーム画像Ｉ_０を並べた例である。
また、上述したように、異なる複数種類の演算処理を行った場合には、各演算結果を並列表示することが可能である。

【0044】

切替表示態様は、図５に示したように、演算結果画像Ｉ_２、元動画Ｉ_０、元フレーム画像Ｉ_０、解析結果動画Ｉ_１、解析フレーム画像Ｉ_１、他の演算結果画像Ｉ_２及び重畳画像Ｉ_３の中から一の画像を所定時間表示した後、当該画像と入れ替えで当該画像とは異なる画像を表示する表示態様である。因みに、図５はある演算結果画像Ｉ_２、他の演算結果画像Ｉ_２、これらの重畳画像Ｉ_３、元動画又は元フレーム画像Ｉ_０の順で変更した場合の例である。

【0045】

なお、並列表示態様と切り替え表示態様を組み合わせることも可能である。すなわち、一のサマリー表示領域１５Ａに複数の画像を切り替えながら表示することも可能である。
また、上記表示態様の説明では、サマリー表示領域１５Ａに表示される画像には、演算結果画像Ｉ_２が必ず含まれるものとしたが、含まない形で表示することも可能である。すなわち、元動画又は元フレーム画像Ｉ_０と解析結果動画又は解析フレーム画像Ｉ_１を重畳表示、並列表示、切替表示することも可能である。

【0046】

上述したように、本実施形態に係る放射線画像処理装置１は、複数の元フレーム画像からなる元動画（放射線動画）の動画データを取得する画像取得手段と、画像取得手段が取得した動画データに所定の解析処理を加え、複数の解析フレーム画像からなる解析結果動画を生成する解析手段と、解析手段が生成した解析結果動画を構成する複数の解析フレーム画像の中から複数の特定解析フレーム画像を選択する選択手段と、選択手段によって選択された複数の特定解析フレーム画像における、共通する座標上にそれぞれ位置する各共通座標画素の信号値に基づいて一の演算信号値を導出する演算処理手段と、前記演算処理手段が座標ごとに生成した複数の演算信号値を用いて前記表示部に演算結果画像を表示する画像表示手段と、を備えたものとなっている。

【0047】

サマリー表示領域１５Ａにおいて、上述したような動画又は静止画像の表示を行うことにより、ユーザーにとって関心の高い演算結果画像Ｉ_２を直感的に得ることができる。その結果、ユーザー自身が診断に適切な表示方法を容易に探ることができ、読影精度が向上する。
また、表示部１５に表示される動画のどこを注力して読影するべきかユーザーが理解し易くなり、効率的に読影できるようになるとともに、読影洩れを防止することができる。その結果、診断の精度を向上させることができる。
また、任意のフレーム画像を選択可能としているので、アーチファクト（ノイズ）が混入していると推定されるフレーム画像を容易に取り除くことができ、解析結果の精度が向上する。
すなわち、動画から、ポイントとなる情報だけを容易に読影することができ、ユーザーの負担が低減する。

【0048】

また、図３に示したように、フレーム画像の選択枚数だけを変えて同一の演算処理（ＭＩＰ処理）を施した複数の演算結果画像Ｉ_２を並列表示すれば、演算処理の妥当性、すなわち、選択フレームの範囲を広くしすぎたのか、もう少し広げられそうか等をユーザーに直感的にすばやく把握させることができる。
また、図４に示したように、異なる種類の複数の動画（フレーム画像の一覧）からそれぞれ選択したフレーム画像に対して演算処理した結果を並列表示したり、複数の動画からそれぞれ選択したフレーム画像に対して演算処理した結果を重畳表示したりすることで、診断に有用な「Ｖ／Ｑ（換気血流比）画像」をユーザーの手で容易に作成することが可能となる。

【0049】

＜第２実施形態＞
次に、本発明の第２実施形態に係る放射線画像処理装置について説明する。図６〜９は本実施形態に係る放射線画像処理装置が行う画像選択の概念図である。なお、図６には、肺野領域の画像を例示したが、本実施形態に係る放射線画像処理装置１は、肺野以外にも、心臓領域、肺動脈、大動脈又は肺気管支等、所定周期で同じ動作を繰り返す部位の動画であれば扱うことが可能である。

【0050】

第１実施形態では、演算処理を施す対象となる特定解析フレーム画像Ｉ_１をユーザーが選択するようにしていたが、本実施形態に係る放射線画像処理装置１Ａは、予め設定された選択パターンに基づいて特定解析フレーム画像Ｉ_１を自動的に選択するようになっている。
このため、本実施形態に係る放射線画像処理装置１Ａは、実行する制御（記憶部１３Ａの記憶内容）が第１実施形態と異なっている。なお、記憶部１３Ａの記憶内容を除く他の構成は第１実施形態と同様である。

【0051】

すなわち、放射線画像処理装置１Ａは、記憶部１３Ａに、自動選択の範囲を規定する選択パターンを複数格納しており、ユーザーが操作部１４に対して行った操作に基づいて、予め用意された複数の選択パターンの中から何れかの選択パターンを設定する。そして、ユーザーが特に注意して読影したい関心領域（例えば図６（ａ）で符号Ｒで示したような領域）を選択すると、設定した選択パターンに基づいて前記特定解析フレーム画像を自動的に選択するようになっている。

【0052】

本実施形態に係る放射線画像処理装置１Ａにおいて設定可能な選択パターンには、大きく分けて（１）極大値（極小値）−極大値（極小値）、（２）極大値（極小値）−極小値（極大値）、（３）極大値（極小値）とその前後所定範囲、の３つがある。
なお、このうちのいずれかのパターンのみが予め設定されている（ユーザーによる設定はできない）ようにしてもよいし、このうちの２パターンの中からユーザーが設定できるものとしてもよい。更に、ここで例示した以外のパターンを含むようにしてもよい。

【0053】

選択パターン（１）が選択されると、放射線画像処理装置１は、少なくとも肺野、心臓領域、肺動脈、大動脈、あるいは肺気管支の一部が含まれる関心領域Ｒ（Region Of Interest；ROI）内の複数の画素信号値の平均値の時間方向の変化を数値化する（図６（ｂ）に示したようにグラフ化する）。そして、図６（ｂ），（ｃ）に示したように、ある極大値（極小値）から、当該極大値（極小値）と隣り合う極大値（極小値）までの間のフレーム画像（一周期分の動作を表す動画；選択範囲ｒ_１）、又はある極大値（極小値）から、当該極大値と隣り合わない極大値（極小値）までの間の解析フレーム画像（複数周期分の動作を表す動画；選択範囲ｒ_２）を特定解析フレーム画像として自動選択する。
なお、ｒ_１とｒ_２のどちらの選択範囲とするかは、ユーザーが選択パターンを設定する際に併せて設定するよう構成してもよいし、制御部１１が自動的に判断するよう構成してもよい。

【0054】

診断対象が図６（ａ）に示したような肺野である場合、この選択範囲ｒ_１，ｒ_２は、吸気位から再度吸気位となるまでの動きを表すことになる。一方、診断対象が心臓である場合には、心収縮期から再度心収縮期となるまでの動きを表すことになる。
肺機能を把握する上では、肺の一呼吸分の動作や心臓の一心拍分の動作等、少なくとも一周期分の情報を読影することが重要である。このため、ここで挙げたような解析フレーム画像の選択パターンを設定することで、呼吸周期や心拍周期に関係する肺機能（例えば、肺への空気の出入り量（１回換気量）、肺野内への血液の出入り量（心拍出量、１回拍出量））を正確に把握することができる。
特に、複数周期に亘る選択範囲ｒ_２を設定することで、呼吸や心拍のバラつきを低減することが可能である。

【0055】

選択パターン（２）が選択されると、放射線画像処理装置１は、関心領域Ｒ内の平均値の時間方向の変化を数値化（グラフ化）する。そして、図７（ａ）に示したように、ある極大値から、当該極大値の後（グラフでは右側）に隣り合う極小値までの間の解析フレーム画像（選択範囲ｒ_３）、又は図７（ｂ）に示したように、ある極小値から、当該極小値の後（グラフでは右側）に隣り合う極大値までの間の解析フレーム画像（選択範囲ｒ_４）を特定解析フレーム画像として自動選択する。
なお、図７に示したように、一のグラフに対して選択範囲ｒ_３又はｒ_４を複数周期分設定することが可能である。
なお、ｒ_３とｒ_４のどちらの選択範囲とするかは、ユーザーが解析フレーム画像の選択パターンを設定する際に併せて設定できるようにしてもよいし、制御部１１が自動的に判断するようにしてもよい。

【0056】

診断対象が肺野である場合には、選択範囲ｒ_３は呼気中、選択範囲は吸気中の動きを表すことになる。一方、診断対象が心臓である場合には、選択範囲ｒ_３又はｒ_４は心臓収縮中又は心臓膨張中の動きを表すことになる。
複数の呼気中（吸気中）あるいは心臓膨張中（心臓収縮中）の解析フレーム画像の結果に対して演算処理を加えることによって、呼吸や心拍のバラつきなどのノイズ影響を低減させることができ、血液が肺野に広がっていく心収縮中の様子、心臓に血液が戻っていく心膨張中の様子、吸気で肺が広がっていく様子、呼気で肺が縮んでいく様子等を精度よく把握することができる。
また、例えば閉塞性疾患患者の息の吐きにくさの評価、拘束性疾患や肥満の方の息の吸い難さの評価、心臓の動き方が不連続になる不整脈などの評価にも有用である。

【0057】

選択パターン（３）が選択されると、放射線画像処理装置１は、関心領域Ｒ内内の平均値の時間方向の変化を数値化（グラフ化）する。そして、図８（ａ）に示したように、ある極大値となる時間ｔから所定時間ｔ_１だけ前に撮影された解析フレーム画像から、当該極大値となる時間から所定時間ｔ_２だけ後に撮影された解析フレーム画像まで（ｔ−ｔ_１〜ｔ＋ｔ_２：選択範囲ｒ_５）、又は図８（ｂ）に示したように、ある極小値となる時間ｔから所定時間ｔ_１だけ前に撮影された解析フレーム画像から、当該極小値となる時間から所定時間ｔ_２だけ後に撮影された解析フレーム画像まで（ｔ−ｔ_１〜ｔ＋ｔ_２：選択範囲ｒ_６）を特定解析フレーム画像として自動選択する。所定時間ｔ_１，ｔ_２は一周期（呼吸周期、心拍周期）の１／４以下とするのが好ましく、より好ましくは１／８以下とする。

【0058】

なお、図８には、所定時間ｔ_１とｔ_２が同じ長さである場合を例示したが、ｔ_１とｔ_２は異なっていてもよい。
また、時間ではなく、極大からの解析フレーム画像数で範囲を規定するようにしてもよい。
更に、図８に示したように、一のグラフに対して選択範囲ｒ_５又はｒ_６を複数周期分設定することが可能である。
呼吸周期（心拍周期）は患者ごとに、スパイロメトリやモストグラフ、心電図等で計測した値を用いても良いし、各国、身長や性別、年齢などで定められた平均的な値を用いても良い。

【0059】

診断対象が肺野である場合には、極大及びその周辺は吸気位辺り、極小及びその周辺は呼気位辺りの動きを表すことになる。一方、診断対象が心臓である場合には、心臓収縮期辺り又は心臓膨張期辺りの動きを表すことになる。
複数の呼気位（吸気位）又は心臓膨張期（心臓収縮期）付近の解析フレーム画像の結果に対して演算処理を加えることによって、血液が肺野に広がり切った心収縮期の様子、心臓に血液が戻り切った心膨張期の様子、呼吸により肺が広がりきった様子、肺が縮み切った様子等を精度よく把握することができる。

【0060】

心臓から肺に血が広がる速度は速い為、肺野内のどこに血液が流れて行ったかを把握するのは通常の解析結果（動画）を見ているだけでは容易ではない。一方、全ての解析フレーム画像に対して演算処理を加えるとアーチファクトも多く付加されるデメリットがある。このため、血液が肺野に広がり切る前後の解析フレーム画像に対して演算処理を加えることで、アーチファクトの少ない肺血流分布画像を得ることが可能となる。つまり、選択パターン（３）によって得られた特定解析フレーム画像は、肺血流に関する疾患の把握に有用である。例えば、肺塞栓など局所的な肺血流機能欠損領域の把握や、心機能が低下する心不全や心肥大の評価に有用である。一方で、（３）によって得られた特定解析フレーム画像は同様な理由でアーチファクトの少ない換気機能画像となる為、局所的な肺換気機能欠損領域の把握や、横隔膜がどこまで広がっているか等、肺の広がり方の評価にも有用である。

【0061】

なお、特定解析フレーム画像を複数周期分選択する場合には、サマリー表示領域１５Ａに各選択範囲の演算結果画像を並列表示するようにしてもよい。
こうすることで、どの呼吸（どの心拍）が最も正常/異常かを判断し易くなる。

【0062】

なお、選択パターン（３）においては、所定時間ｔ_１，ｔ_２のいずれか一方を０としてもよい。すなわち、図９（ａ）に示したように、ある極大値となる時間ｔから所定時間ｔ_２だけ後に撮影された解析フレーム画像まで（極大〜ｔ＋ｔ_２：選択範囲ｒ_７）、ある極小値となる時間ｔから所定時間ｔ_２だけ後に撮影された解析フレーム画像まで（極小〜ｔ＋ｔ_２：選択範囲ｒ_８）、ある極大値となる時間ｔから所定時間ｔ_１だけ前に撮影された解析フレーム画像から当該極大まで（ｔ−ｔ_１〜極大：選択範囲ｒ_９）、又はある極小値となる時間ｔから所定時間ｔ_１だけ前に撮影された解析フレーム画像から当該極小まで（ｔ−ｔ_１〜極小：選択範囲ｒ_１０）を特定解析フレーム画像として自動選択するようにしてもよい。

【0063】

なお、本実施形態では、関心領域Ｒ内の画素信号値の平均値の時間方向の変化を数値化（グラフ化）し、そこから特定解析フレーム画像を選択するようにしたが、他の統計値（最大値、最小値、最頻値、積算値、標準偏差等）を用いて特定解析フレーム画像を選択するようにしてもよい。

【0064】

また、横隔膜移動量の時間変化や、肺野面積の時間変化、肋間距離の変化量の時間変化、肺野領域の水平方向の広がりの時間変化、心臓面積の時間変化、大動脈弓の位置の時間変化等の形態情報の時間方向変化を元に、特定解析フレーム画像を選択するようにしても良い。
撮影対象が、例えば肺野や心臓領域にて脂肪や筋肉などが不均一に分布している体型の患者（肥満の方や女性等に多い）であったり、撮影中に体が上下左右に大きく動いてしまう（体動が大きい）患者であったりする場合、画素値を用いた時間変化情報が正確に呼吸状態や心拍状態を反映していない場合がある。しかし、形態情報の時間方向の変化を用いることで、このような患者に対しても特定解析フレーム画像を精度よく選択することができるようになる。

【0065】

＜第３実施形態＞
次に、本発明の第３実施形態に係る放射線画像処理装置１Ｂについて説明する。図１０は放射線画像処理装置１Ｂにおける画像処理の概念図であり、図１１，１２は処理結果の表示例である。
なお、図１０には、肺野領域の画像を例示したが、本実施形態に係る放射線画像処理装置１Ｂは、肺野以外の部位を撮影した動画についても処理・表示することが可能である。

【0066】

第１，第２実施形態では、サマリー表示領域１５Ａに表示する内容を画像のみとしたが、本実施形態に係る放射線画像処理装置１Ｂは、所定の定量値をサマリー表示領域１５Ａに表示するようになっている。
このため、本実施形態に係る放射線画像処理装置１Ｂは、実行する制御（記憶部１３Ｂの記憶内容）が第１，第２実施形態と異なっている。なお、記憶部１３Ｂの記憶内容を除く他の構成は第１，第２実施形態と同様である。

【0067】

放射線画像処理装置１Ｂの制御部１１は、複数の演算信号値を算出した後、演算結果画像の関心領域Ｒ内に位置することになる複数の画素の信号値に基づいて定量値Ｑを導出するようになっている。すなわち、制御部１１は、本発明における定量値算出手段をなす。なお、関心領域Ｒは、本発明における特定領域である。以下、関心領域Ｒ内に位置することになる画素を領域内画素と称する。なお、関心領域Rは複数設定されていても良く、その場合は、各々の関心領域Rに対応した定量値Ｑを各々個別に導出しても良いし、複数設定された各々の関心領域Rをまとめて一つとみなして、一つの定量値Ｑを導出しても良い。
定量値Ｑの導出の仕方としては、例えば、最大値処理、最小値処理、Ｐｋ−Ｐｋ値処理、標準偏差処理、平均値処理、最頻値処理、積算値処理等がある。

【0068】

最大（小）値処理は、演算結果画像の関心領域Ｒ内に位置することになる複数の領域内画素の各信号値の中から最大（小）値を抽出して定量値Ｑとする処理である。この処理により得られた定量値Ｑを用いれば、経過観察（換気、血流の悪化・良化等）や、左右肺機能の比較等を行うことができる。特に、低い計算コストで、肺換気機能や肺血流機能、心機能に関する特徴量の最大（小）値を把握できるので、例えば、治療前後などの経過観察（肺活量の変化や血圧の変化を見る）に有用である。

【0069】

最頻値処理は、複数の領域内画素の各信号値の中から最頻値を抽出して定量値Ｑとする処理である。この処理により得られた定量値Ｑを用いれば、経過観察や左右肺機能の比較等を行うことができる。特に、演算結果画像上に特異的に発生する画欠等のアーチファクトの影響を低減できる為、関心領域Ｒ内の肺換気機能や肺血流機能、心機能を精度よく把握することができる。

【0070】

Ｐｋ−Ｐｋ値処理は、複数の領域内画素の各信号値の中から最大値と最小値を抽出してその差分を算出し、当該領域内画素の数に対する最大値と最小値との差の割合を算出して定量値Ｑとする処理である。この処理により得られた定量値Ｑを用いれば、経過観察や左右肺機能の比較等を行うことができる。特に、肺換気機能や肺血流機能、心機能に関する特徴量の変化量を把握できるので、関心領域Ｒ内における上記各機能の特徴量のバラつき（局所的に肺換気機能や肺血流機能が異常を示す疾患）を把握しやすくなる。

【0071】

標準偏差処理は、複数の領域内画素の各信号値の標準偏差を算出して定量値Ｑとする処理である。この処理により得られた定量値Ｑを用いれば、経過観察や左右肺機能の比較等を行うことができる。特に、低い計算コストで、肺換気機能や肺血流機能、心機能に関する特徴量の変化量を把握できるので、関心領域Ｒ内の上記各機能のバラつきを把握しやすくなる。

【0072】

平均値処理は、複数の領域内画素の各信号値の平均値を算出して定量値Ｑとする処理である。この処理により得られた定量値Ｑを用いれば、経過観察や左右肺機能比較等を行うことができる。特に、低い計算コストで、関心領域Ｒ内の肺換気機能や肺血流機能、心機能を把握することができる。

【0073】

積算値処理は、複数の領域内画素の各信号値の積算値を算出して定量値Ｑとする処理である。この処理により得られた定量値Ｑを用いれば、経過観察や左右肺機能比較等を行うことができる。特に、関心領域Ｒの取り方にわずかな変化があっても算出される定量値Ｑに変化が少ない為、測定誤差が生じにくくなる。また、関心領域Ｒ内のトータルでの肺機能、心機能を把握しやすくなる。

【0074】

なお、定量値Ｑは、関心領域Ｒ内の領域内画素の信号値に基づいて導出するのではなく、演算結果画像の全画素のうち、信号値が所定閾値以上、所定閾値以下又は所定範囲内にある画素の個数（面積）としてもよい。
このようにすれば、例えば、肺全体（全体肺面積）に対して、正常又は異常な肺領域が全体の何割あるか等の情報を容易に把握でき、肺の局所的な機能分布(どこに肺機能の低下、欠損があるか)を評価しやすくなる。

【0075】

本実施形態において、関心領域Ｒとして設定できるのは、全画像領域、肺野内、心臓領域内、肺動脈内、大動脈内、肺気管支内である。
全画像領域を関心領域Ｒとした定量値Ｑを用いれば、計算コストをかけずに、手軽に大局的な肺換気機能、肺血流機能、心機能を把握することができる。
肺野内を関心領域Ｒとした定量値Ｑを用いれば、肺換気機能（呼吸に伴う濃度の変化、横隔膜の動き、肋骨の動き、鎖骨の動き、肩甲骨の動きなど）、肺血流機能（心拍に伴う濃度の変化、血管の動き）を把握することができる。

【0076】

心臓領域内を関心領域Ｒとした定量値Ｑを用いれば、心機能（心壁の移動量（移動速度、加速度）、心拍に伴う濃度の変化量（濃度変化速度、濃度変化加速度））を把握することができる。
肺動脈内を関心領域Ｒとした定量値Ｑを用いれば、肺血流機能、心機能を把握することができる。
大動脈内を関心領域Ｒとした定量値Ｑを用いれば、体循環の血流機能（心拍に伴う濃度の変化、血管の動き）、心機能を把握することができる。
肺気管支内を関心領域Ｒとした定量値Ｑを用いれば、肺換気機能の把握、気道の肥厚の程度、気道収縮の程度、無気肺/胸水/腫瘤、のう胞/肺気腫の有無などを把握することができる。

【0077】

なお、定量値Ｑの表示は、図１０（ａ）に示したように単独で行ってもよいし、図１０（ｂ），（ｃ）に示したように、画像と併せて行うようにしてもよい。
画像Ｉ_２，Ｉ_３と併せて表示する場合には、図１０（ｂ）に示したように並列表示する他、切替表示、重畳表示することも可能である。
また、図１０（ｃ）に示したように、画像Ｉ_２，Ｉ_３と関心領域Ｒを重畳表示し、更にその中に定量値Ｑを表示するようにしてもよい。このようにすれば、左右の肺機能を直感的に把握することが可能となる。

【0078】

また、制御部１１は、種類の異なる複数のフレーム画像群からそれぞれフレーム画像を選択（サマリー表示領域１５Ａにドラッグ及びドロップ等）した場合は、各動画についての演算結果画像Ｉ_２を生成するだけでなく、各動画についての定量値Ｑも算出する。
このとき、換気解析動画（解析結果動画）について算出した定量値Ｑと、血流解析動画（解析結果動画）について算出した定量値Ｑを用いた特徴量（比や差分など）を算出し、図１１に示したようにＶ／Ｑ比（換気血流比）を示すＶ／Ｑ画像Ｖを表示するようにしても良い。
このようにすれば、ユーザーが正常異常の判断を行いやすくなるだけでなく、経時的な比較が可能となる。

【0079】

また、図１２（ａ）に示したように、左右の肺ごとに定量値Ｑを算出し、さらに各定量値Ｑの左右肺比率を算出して定量値Ｑとともに表示するようにしても良い。こうすることで、左右肺の肺機能の不均等がわかり易くなり、正常と異常の把握が行いやすくなる。
また、心臓領域を診断対象とする場合には、図１２（ｂ）に示したように、心室、心房、右心系、左心系ごとに定量値Ｑ及び比率を算出し、表示するようにしてもよい。こうすることで、心機能の不均衡がわかり易くなり、正常と異常の把握が行いやすくなる。

【0080】

また、本実施形態において算出した定量値Ｑについては、演算結果画像Ｉ_２を保存する際、画像のヘッダーに保存するようにしても良い。このようにすれば、演算結果画像Ｉ_２を再度読み込む際、演算結果画像Ｉ_２のファイル名一覧でヘッダー情報の有無を見ることにより、必要な場合のみ定量値Ｑを表示させることができる。
また、画像とは別データ（テキストやＣＳＶファイルなど）として保存しても良い。そうすることで、定量値Ｑのみを再度見直すことが容易になるし、不要に画像を毎回開く手間が減り、経過観察等の作業効率が向上する。

【0081】

（変形例１）
なお、上記実施形態においては、サマリー表示領域１５Ａに選択した複数の解析フレーム画像に基づく演算結果画像や定量値Ｑを表示するようにしたが、図１３に示したように、単に選択（サマリー表示領域１５Ａにドラッグ及びドロップ）した複数のフレーム画像のみを表示するようにしてもよい。
なお、本変形例においても、上記実施形態と同様、選択するフレーム画像は必ずしも連続のフレーム画像である必要はない。また、サマリー表示領域１５Ａへの表示は、フレーム画像の選択後自動的に行ってもよいし、ユーザーの操作（例えばＯＫボタンをクリックする等）に基づいて行うようにしてもよい。

【0082】

また、異なる複数種類のフレーム画像の一覧からそれぞれフレーム画像を選択した場合は、図１４（ａ）に示したように重畳画像Ｉ_３の表示（ここでは、元フレーム画像と解析フレーム画像の重畳表示）の他、図１４（ｂ）に示した並列表示（ここでは、解析フレーム画像Ｉ_２、元フレーム画像Ｉ_０及び重畳画像Ｉ_３の並列表示）、図１４（ｃ）に示した切替表示（ここでは、解析フレーム画像Ｉ_２、元フレーム画像Ｉ_０及び重畳画像Ｉ_３の切替表示）を行う。
なお、異なる複数種類の動画の中からそれぞれフレーム画像を選択する際、それぞれのフレーム画像数は同一であっても異なっていても良い。

【0083】

本変形例のようにすれば、ユーザーにとって関心の高いフレーム画像のみがサマリー表示領域１５Ａに表示される為、ユーザーは読影に注力しやすくなり、読影精度が向上する。
特に、１呼吸分など呼吸周期の整数倍に相当するフレーム画像のみをサマリー表示領域１５Ａに表示することで、呼吸が正しく行われているかを評価しやすくなる。
また、１心拍分など、心拍周期の整数倍に相当するフレーム画像のみをサマリー表示領域１５Ａに表示することで、心臓が正常に動いているかを評価しやすくなる。
また、吸気始めあるいは、呼気始めに相当するフレーム画像のみをサマリー表示領域１５Ａに表示することで、呼吸器疾患症例に対しての、息の吸い難さや息の吐きにくさを評価しやすくなる。
また、心膨張あるいは心収縮に相当するフレーム画像のみをサマリー表示領域１５Ａに表示することで、心疾患症例に対しての、血液の拍入のし易さ、血液の拍出のし易さ、不整脈の有無や心肥大の有無を評価しやすくなる。
また、肺機能が正常あるいは異常と判断したフレーム画像群をサマリー表示領域１５Ａに表示させたうえで、表示部１５に元々表示されているフレーム画像の一覧と比較することで、表示部１５に元々表示されているフレーム画像群が総じて正常なのか異常なのかを判断しやすくなる。すなわち、画像表示部に元々表示されているフレーム画像群の読影効率が向上する。

【0084】

なお、図１５（ａ）に示したように、換気解析結果と、血流解析結果をサマリー表示領域１５Ａに重畳表示することで、容易に換気機能と血流機能を比較する画像（Ｖ／Ｑ画像）と同等の画像を得ることができる。
このようにすれば、ユーザーが選択したフレーム画像のみを用いて換気機能と血流機能の比較が行える為、例えばアーチファクトの影響をあまり受けていないＶ／Ｑ画像を得ることができ、診断精度が向上する。
また、図１５（ｂ）に示したように、換気解析結果と、血流解析結果、これらの重畳画像及び元画像の並列表示や、切替表示も可能である。

【0085】

（変形例２）
また、上記実施形態では、表示部１５に表示されたフレーム画像の一覧の中から、サマリー表示領域１５Ａに解析結果動画等を表示するためのフレーム画像を選択するようにしていたが、図１６に示したように、表示部１５において動画を再生し、動画から必要な画像を直接ドラックしサマリー表示領域１５Ａにドロップするようにしても良い。

【0086】

また、図１７に示したように、各フレーム画像にチェックボックスＢを設け、チェックＣを入れたフレーム画像がサマリー表示領域１５Ａに表示されるようにしても良い。このようにすれば、サマリー表示領域１５Ａに表示されたフレーム画像又は動画を見ながら、チェックＣを入れるフレーム画像を選べるため、ユーザーが選択したフレーム画像が妥当かどうか（診断にふさわしいフレーム画像かどうか）判断を素早く行うことが可能となる。

【0087】

（変形例３）
また、図１８に示したように、選択されサマリー表示領域１５Ａに表示されるフレーム画像については、選択されていないフレーム画像や動画とは異なるルックアップテーブル（Lookup table；ＬＵＴ）を適用した画像Ｉ_Ｌを表示してもよい。具体的には、濃度、コントラスト、カラーマップの変更し表示する。
このようにすれば、ユーザーが関心の高い画像が、ユーザーの関心の高い色づけでサマリー表示領域１５Ａに表示される為、ユーザーは読影に注力しやすくなり、読影精度が向上する。

【符号の説明】

【0088】

１，１Ａ，１Ｂ 放射線画像処理装置
１１ 制御部（画像取得手段、解析手段、選択手段、演算処理手段、画像表示手段）
１２ 通信部
１３，１３Ａ，１３Ｂ 記憶部
１４ 操作部
１５ 表示部
１５Ａ サマリー表示領域
１６ バス
Ｉ_０ 元動画、元フレーム画像
Ｉ_１ 解析結果動画、解析フレーム画像
Ｉ_２ 演算結果画像
Ｉ_３ 重畳画像
Ｒ 関心領域
ｒ_１〜ｒ_１０ （フレーム画像の）選択範囲
ｔ 時間
ｔ_１，ｔ_２ 所定時間

【図1】

【図2】

【図3】

【図4】

【図5】

【図6】

【図7】

【図8】

【図9】

【図10】

【図11】

【図12】

【図13】

【図14】

【図15】

【図16】

【図17】

【図18】