特許 7073661 動態解析装置及び動態解析システム

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

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

(11)【特許番号】

(24)【登録日】2022-05-16

(45)【発行日】2022-05-24

(54)【発明の名称】動態解析装置及び動態解析システム

(51)【国際特許分類】

   A61B 6/00 20060101AFI20220517BHJP

【ＦＩ】

A61B6/00 350C

A61B6/00 330A

【請求項の数】 8

(21)【出願番号】P 2017185574

(22)【出願日】2017-09-27

(65)【公開番号】P2019058368

(43)【公開日】2019-04-18

【審査請求日】2020-03-18

(73)【特許権者】

【識別番号】000001270

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

(74)【代理人】

【識別番号】110001254

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

(72)【発明者】

【氏名】松谷 哲嗣

【審査官】佐藤 秀樹

(56)【参考文献】

【文献】米国特許出願公開第２００８／０２６９５９２（ＵＳ，Ａ１）

【文献】特開２０１２－１９２２５５（ＪＰ，Ａ）

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

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

Ａ６１Ｂ ６／００－６／１４

Ｇ０６Ｔ ７／００－７／９０

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

【請求項1】

被写体における血管を含む部位を放射線撮影することにより得られた動態画像内の複数の領域において前記血管の位置や形状に関する指標値の経時変化を測定する指標値測定手段と、
前記指標値測定手段が測定した各領域の指標値の経時変化に基づいて、各領域における血管の伸長率を算出する伸長率算出手段と、
前記伸長率算出手段が算出した各領域における伸長率を一覧表示する伸長率表示手段と、を備え、
前記血管を含む部位は肺野であり、
前記指標値測定手段は、肺野全体に分布するように配置され、前記動態画像に含まれる各フレーム画像において同一部位に付された複数のマーカーのうち、肺野内に設定された大きさ及び座標が一定の領域内にあるマーカーの数を指標値として測定することを特徴とする動態解析装置。

【請求項2】

被写体における血管を含む部位を放射線撮影することにより得られた動態画像内の複数の領域において前記血管の位置や形状に関する指標値の経時変化を測定する指標値測定手段と、
前記指標値測定手段が測定した各領域の指標値の経時変化に基づいて、各領域における血管の伸長率を算出する伸長率算出手段と、
前記伸長率算出手段が算出した各領域における伸長率を一覧表示する伸長率表示手段と、
前記動態画像に所定の画像処理を施すことにより肺野内の肺血管を検出する血管検出手段と、
を備え、
前記血管を含む部位は肺野であり、
前記指標値測定手段は、肺野全体に分布するように配置され、前記動態画像に含まれる各フレーム画像において同一部位に付された複数のマーカーのうち、肺野内に設定された大きさ及び座標が一定の領域内にあるマーカーの数を指標値として測定し、
前記指標値として、前記領域に占める肺血管領域の密度を測定することを特徴とする動態解析装置。

【請求項3】

前記動態画像に所定の画像処理を施すことにより肺野内の肺血管を検出する血管検出手段を備え、
前記指標値測定手段は、前記指標値として、一本の肺血管上に設定された二点の距離を測定することを特徴とする請求項１から請求項２のいずれか一項に記載の動態解析装置。

【請求項4】

前記動態画像に所定の画像処理を施すことにより肺野内の肺血管を検出する血管検出手段を備え、
前記指標値測定手段は、前記指標値として、二本の異なる肺血管上にそれぞれ設定された二点の距離を測定することを特徴とする請求項１から請求項３のいずれか一項に記載の動態解析装置。

【請求項5】

前記動態画像に所定の画像処理を施すことにより肺野内の肺血管を検出する血管検出手段を備え、
前記指標値測定手段は、前記指標値として、肺血管の周辺に設定された三つ以上の点を直線で結ぶことにより構成される多角形の面積を測定することを特徴とする請求項１から請求項４のいずれか一項に記載の動態解析装置。

【請求項6】

前記伸長率算出手段は、二時点における指標値の差、指標値の増加速度、加速度及び増加開始タイミングのいずれかに基づいて前記伸長率を算出することを特徴とする請求項１から５のいずれか一項に記載の動態解析装置。

【請求項7】

前記伸長率表示手段は、伸長率に対応する色で領域を塗り分けた撮影部位の画像を表示可能であることを特徴とする請求項１から６のいずれか一項に記載の動態解析装置。

【請求項8】

請求項１から７のいずれか一項に記載の動態解析装置と、
被写体における血管を含む部位を放射線撮影して動態画像を生成し、前記動態解析装置へ出力する撮影装置と、を備えることを特徴とする動態解析システム。

【発明の詳細な説明】

【技術分野】

【0001】

本発明は、動態解析装置、及びこの装置を備える動態解析システムに関する。

【背景技術】

【0002】

従来のフィルム／スクリーンや輝尽性蛍光体プレートを用いた放射線（Ｘ線）の静止画撮影及び診断に対し、ＦＰＤ（flat panel detector）等の半導体イメージセンサーを利用して診断対象の部位（対象部位という）の動態画像を撮影し、診断に応用する試みがなされるようになってきている。具体的には、半導体イメージセンサーの画像データの読取・消去の応答性の早さを利用し、半導体イメージセンサーの読取・消去のタイミングと合わせて放射源からパルス状の放射線を連続的に照射し、１秒間に複数回の撮影を行って、対象部位の動態を撮影する。撮影により取得された一連の複数枚の画像を順次表示することにより、医師は対象部位の一連の動きを観察することが可能となる。

【0003】

肺の診断においては、肺の機能（換気機能や肺血流機能）が低下している箇所がないかを観察することは重要である。しかし、医師が動態画像を観察して目視で機能の異常個所を認識することは難しい。特に、肺の呼吸運動や心臓の拍動には個人差があり、個人差を考慮しながら換気機能や肺血流機能の異常箇所を視認することは困難である。
そこで、動態撮影で得られた動態画像を解析して診断支援情報を生成し、早期診断に向けて医師に提供することが提案されている。例えば、特許文献１には、呼吸動態画像のフレーム間差分画像を作成し、フレーム間差分画像の対応画素毎に絶対値が最大の差分画素を求めて差分画素からなる最大値画像を作成し、その画像を特定の呼吸位相の画像に重ねて表示することが記載されている。このような技術を用いれば、呼吸状態に関連した画像の濃度変化を基に、肺疾患の局所的な位置を検出することができる。

【先行技術文献】

【特許文献】

【0004】

【文献】特許第４４０４２９１号公報

【発明の概要】

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

【0005】

特許文献１に記載されたような画像の濃度変化に基づいて肺疾患を検出するシステムでは、呼吸に伴う肋骨の移動や体動、体厚の変化により、画像の濃度変化にアーティファクトが生じてしまうことがあった。
具体的には、呼吸に伴って肋骨が所定の関心領域（ＲＯＩ）内へ出入りすることで、濃度変化量が大きくなり、疾患部位を健常な状態であると誤判断してしまう要因となっていた。

【0006】

本発明の課題は、上記課題に鑑みてなされたもので、動態画像を解析して肺疾患を検出する動態解析システムにおいて、肋骨の移動や体動、体厚の変化の影響を受けることなく、肺疾患を正確に検出できるようにすることを目的とする。

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

【0007】

上記課題を解決するため、本発明に係る動態解析装置は、
被写体における血管を含む部位を放射線撮影することにより得られた動態画像内の複数の領域において前記血管の位置や形状に関する指標値の経時変化を測定する指標値測定手段と、
前記指標値測定手段が測定した各領域の指標値の経時変化に基づいて、各領域における血管の伸長率を算出する伸長率算出手段と、
前記伸長率算出手段が算出した各領域における伸長率を一覧表示する伸長率表示手段と、を備え、
前記血管を含む部位は肺野であり、
前記指標値測定手段は、肺野全体に分布するように配置され、前記動態画像に含まれる各フレーム画像において同一部位に付された複数のマーカーのうち、肺野内に設定された大きさ及び座標が一定の領域内にあるマーカーの数を指標値として測定することを特徴とする。

【発明の効果】

【0008】

本発明によれば、肋骨の移動や体動、体厚の変化の影響を受けることなく、肺疾患を正確に検出することができる。

【図面の簡単な説明】

【0009】

【図1】本発明の実施形態における動態解析システムの全体構成を示す図である。

【図2】図１の動態解析システムを構成する撮影用コンソールが実行する撮影制御処理を表すフローチャートである。

【図3】図１の動態解析システムを構成する診断用コンソールが実行する画像解析処理を表すフローチャートである。

【図4】指標値を測定するためのマーカーの配置例である。

【図5】指標値を測定するためのマーカーの配置例である。

【図6】指標値を測定するための画像の処理例である。

【図7】測定する指標値の経時変化の例である。

【図8】伸長量の比較方法を表す模式図である。

【図9】伸長量の一覧表示の例である。

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

【0010】

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

【0011】

〔動態解析システム１００の構成〕
まず、構成を説明する。
図１に、本実施形態における動態解析システム１００の全体構成を示す。
図１に示すように、動態解析システム１００は、撮影装置１と、撮影用コンソール２とが通信ケーブル等により接続され、撮影用コンソール２と、診断用コンソール３とがＬＡＮ（Local Area Network）等の通信ネットワークＮＴを介して接続されて構成されている。動態解析システム１００を構成する各装置は、ＤＩＣＯＭ（Digital Image and Communications in Medicine）規格に準じており、各装置間の通信は、ＤＩＣＯＭに則って行われる。

【0012】

なお、撮影装置１と撮影用コンソール２とが通信ネットワークを介して接続されていてもよいし、撮影用コンソール２と診断用コンソール３が通信ケーブルで接続されていてもよい。
また、図１には、撮影用コンソール２と診断用コンソール３を別々に設けたものを例示したが、これらは一体（一の装置が撮影用コンソールと診断用コンソールの機能を兼ね備えたもの）になっていてもよい。

【0013】

〔撮影装置１の構成〕
撮影装置１は、例えば、呼吸運動に伴う肺の膨張及び収縮の形態変化、心臓の拍動等の、周期性（サイクル）を持つ動態を撮影する撮影手段である。動態撮影とは、被写体に対し、Ｘ線等の放射線をパルス状にして所定時間間隔で繰り返し照射するか（パルス照射）、もしくは、低線量率にして途切れなく継続して照射する（連続照射）ことで、動態を示す複数の画像を取得することをいう。動態撮影により得られた一連の画像を動態画像と呼ぶ。また、動態画像を構成する複数の画像のそれぞれをフレーム画像と呼ぶ。なお、以下の実施形態では、パルス照射により動態撮影を行う場合を例にとり説明する。

【0014】

放射線源１１は、被写体Ｍを挟んで放射線検出部１３と対向する位置に配置され、放射線照射制御装置１２の制御に従って、被写体Ｍに対し放射線（Ｘ線）を照射する。
放射線照射制御装置１２は、撮影用コンソール２に接続されており、撮影用コンソール２から入力された放射線照射条件に基づいて放射線源１１を制御して放射線撮影を行う。撮影用コンソール２から入力される放射線照射条件は、例えば、パルスレート、パルス幅、パルス間隔、１撮影あたりの撮影フレーム数、Ｘ線管電流の値、Ｘ線管電圧の値、付加フィルター種等である。パルスレートは、１秒あたりの放射線照射回数であり、後述するフレームレートと一致している。パルス幅は、放射線照射１回当たりの放射線照射時間である。パルス間隔は、１回の放射線照射開始から次の放射線照射開始までの時間であり、後述するフレーム間隔と一致している。

【0015】

放射線検出部１３は、ＦＰＤ等の半導体イメージセンサーにより構成される。ＦＰＤは、例えば、ガラス基板等を有しており、基板上の所定位置に、放射線源１１から照射されて少なくとも被写体Ｍを透過した放射線をその強度に応じて検出し、検出した放射線を電気信号に変換して蓄積する複数の検出素子（画素）がマトリックス状に配列されている。各画素は、例えばＴＦＴ（Thin Film Transistor）等のスイッチング部を備えて構成されている。ＦＰＤにはＸ線をシンチレーターを介して光電変換素子により電気信号に変換する間接変換型、Ｘ線を直接的に電気信号に変換する直接変換型があるが、何れを用いてもよい。
放射線検出部１３は、被写体Ｍを挟んで放射線源１１と対向するように設けられている。

【0016】

読取制御装置１４は、撮影用コンソール２に接続されている。読取制御装置１４は、撮影用コンソール２から入力された画像読取条件に基づいて放射線検出部１３の各画素のスイッチング部を制御して、当該各画素に蓄積された電気信号の読み取りをスイッチングしていき、放射線検出部１３に蓄積された電気信号を読み取ることにより、画像データを取得する。この画像データがフレーム画像である。そして、読取制御装置１４は、取得したフレーム画像を撮影用コンソール２に出力する。画像読取条件は、例えば、フレームレート、フレーム間隔、画素サイズ、画像サイズ（マトリックスサイズ）等である。フレームレートは、１秒あたりに取得するフレーム画像数であり、パルスレートと一致している。フレーム間隔は、１回のフレーム画像の取得動作開始から次のフレーム画像の取得動作開始までの時間であり、パルス間隔と一致している。

【0017】

ここで、放射線照射制御装置１２と読取制御装置１４は互いに接続され、互いに同期信号をやりとりして放射線照射動作と画像の読み取りの動作を同調させるようになっている。

【0018】

〔撮影用コンソール２の構成〕
撮影用コンソール２は、放射線照射条件や画像読取条件を撮影装置１に出力して撮影装置１による放射線撮影及び放射線画像の読み取り動作を制御するとともに、撮影装置１により取得された動態画像を撮影技師等の撮影実施者によるポジショニングの確認や診断に適した画像であるか否かの確認用に表示する。
撮影用コンソール２は、図１に示すように、制御部２１、記憶部２２、操作部２３、表示部２４、通信部２５を備えて構成され、各部はバス２６により接続されている。

【0019】

制御部２１は、ＣＰＵ（Central Processing Unit）、ＲＡＭ（Random Access Memory）等により構成される。制御部２１のＣＰＵは、操作部２３の操作に応じて、記憶部２２に記憶されているシステムプログラムや各種処理プログラムを読み出してＲＡＭ内に展開し、展開されたプログラムに従って後述する撮影制御処理を始めとする各種処理を実行し、撮影用コンソール２各部の動作や、撮影装置１の放射線照射動作及び読み取り動作を集中制御する。

【0020】

記憶部２２は、不揮発性の半導体メモリーやハードディスク等により構成される。記憶部２２は、制御部２１で実行される各種プログラムやプログラムにより処理の実行に必要なパラメーター、或いは処理結果等のデータを記憶する。例えば、記憶部２２は、図２に示す撮影制御処理を実行するためのプログラムを記憶している。また、記憶部２２は、撮影部位に対応付けて放射線照射条件及び画像読取条件を記憶している。各種プログラムは、読取可能なプログラムコードの形態で格納され、制御部２１は、当該プログラムコードに従った動作を逐次実行する。

【0021】

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

【0022】

表示部２４は、ＬＣＤ（Liquid Crystal Display）やＣＲＴ（Cathode Ray Tube）等のモニターにより構成され、制御部２１から入力される表示信号の指示に従って、操作部２３からの入力指示やデータ等を表示する。

【0023】

通信部２５は、ＬＡＮアダプターやモデムやＴＡ（Terminal Adapter）等を備え、通信ネットワークＮＴに接続された各装置との間のデータ送受信を制御する。

【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】

指標値の差は、図７（ａ）に示したように、複数得られる測定値を、時間（フレーム番号）を横軸、測定値を縦軸とするグラフにしたときの、所定時点における測定値（極大値、極小値又は中間値）と所定時点よりも後の時点における測定値極大値、極小値又は中間値）との差である。なお、それらの比を用いてもよい。
増加速度は、図７（ｂ）に示したように、測定値が増加する期間における、単位時間当たりの測定値の増加量（グラフの傾き）である。
加速度は、図７（ｃ）に示したように、測定値が増加しているある時点における測定値の増加量（グラフの接線の傾き）である。
伸長タイミングは、図７（ｄ）に示したように、増加を始める（グラフの接線の傾きが所定以上となる）タイミングである。

【0054】

得られた測定値は、指標値の測定にマーカーＭを用いた場合には、マーカーＭ上に記憶し、文字・色情報により、動態画像上にマッピングする。
一方、指標値の測定にマーカーＭを用いなかった場合には、測定した位置（座標）毎に測定値を記憶し、文字・色情報により、動態画像上にマッピングするようにする。
以上のような伸長率を算出する本実施形態の診断用コンソール３は、本発明における伸長率算出手段ということになる。

【0055】

伸長率を算出した後は、算出した各領域における伸長率を一覧表示する（ステップＳ１３）。一覧表示の仕方としては、表の形であってもよいし、撮影部位の画像（静止画でも動画でもよい）の対応する箇所に表示する形であってもよい。
ユーザーは、一覧表示を見て、ある測定対象領域Ｒの伸長率を、他の領域Ｒ２における伸長率と比較する。具体的には、図８（ａ）に示したように、当該測定対象領域Ｒの周囲の領域Ｒ２における伸長率と比較する。又は、図８（ｂ）に示したように、左右一方側の肺Ｌにおけるある測定対象領域Ｒの伸長率を、他方側の肺Ｌにおける対応する（左右対称の位置にある）領域Ｒ２の伸長率と比較する。

【0056】

健全な（疾患のない）肺において、ある特定の部位の伸長率が周囲に比べて大きく異なるということは考えにくいため、ある測定対象領域Ｒの伸長率が、周囲の領域の伸長率から大きくかい離している（小さい、遅い）場合、肺野におけるその測定対象領域Ｒに対応する部位に疾患等の異常が存在する可能性が高いということになる。
以上のような伸長率を一覧表示する診断用コンソール３は、本発明における伸長率表示手段ということになる。
なお、他の領域の伸長率と比較するのではなく、所定の基準値を予め保持しておき、その基準値と比較するようにしてもよい。

【0057】

また、本実施形態の診断用コンソール３は、図９に示したように、伸長率に対応する色で領域を塗り分けた撮影部位の静止画像（カラーマッピング）を表示部３４に表示することが可能となっている。相対的に伸長率が低い測定対象領域Ｒを周囲と異なる色（ここでは白）にすることで、数値をそのまま一覧表示するよりも伸長率が他からかい離している領域を見つけやすくなる。
その際、伸長率の数値を色による表示と併せて行うようにしてもよい。伸長率は、常に全小領域Ｒ１の値を一覧表示することとしてもよいし、カラーマッピング表示の画像上で操作部３３により指定された（クリックされた又はカーソルを合わせた）領域の伸長率だけを表示することとしてもよい。

【0058】

以上説明したように、診断用コンソール３は、被写体における肺野（血管Ｖを含む部位）を放射線撮影することにより得られた動態画像の画像データを入力する入力手段と、入力手段により入力された動態画像内の複数の測定対象領域Ｒにおいて肺血管Ｖの位置や形状に関する指標値の経時変化を測定する指標値測定手段と、指標値測定手段が測定した各測定対象領域Ｒの指標値の経時変化に基づいて、各測定対象領域Ｒにおける肺血管Ｖの伸長率を算出する伸長率算出手段と、血管伸長率計算部が算出した各測定対象領域Ｒにおける伸長率を一覧表示する伸長率表示手段と、を備える。

【0059】

これにより、診断用コンソール３は、画像の濃度変化を用いることなく、肺野における硬い領域を検出することができるため、肋骨の移動や体動、体厚の変化の影響を受けることなく、肺疾患を正確に検出することができる。

【0060】

なお、本実施形態における記述は、本発明に係る動態解析システムの一例であり、これに限定されるものではない。
例えば、上記実施形態においては、人体の胸部の動態を撮影することにより得られた胸部動態画像から対象部位として肺や肺血管Ｖの柔らかさを評価する場合を例にとり説明したが、これに限定されず、他の部位を撮影した動態画像から他の部位の柔らかさを評価する場合に適用することとしてもよい。

【0061】

また、例えば、上記の説明では、本発明に係るプログラムのコンピュータ読み取り可能な媒体としてハードディスクや半導体の不揮発性メモリー等を使用した例を開示したが、この例に限定されない。その他のコンピュータ読み取り可能な媒体として、ＣＤ-ＲＯＭ
等の可搬型記録媒体を適用することが可能である。また、本発明に係るプログラムのデータを通信回線を介して提供する媒体として、キャリアウエーブ（搬送波）も適用される。

【0062】

その他、動態解析システム１００を構成する各装置の細部構成及び細部動作に関しても、本発明の趣旨を逸脱することのない範囲で適宜変更可能である。

【符号の説明】

【0063】

１００ 動態解析システム
１ 撮影装置
１１ 放射線源
１２ 放射線照射制御装置
１３ 放射線検出部
１４ 読取制御装置
２ 撮影用コンソール
２１ 制御部
２２ 記憶部
２３ 操作部
２４ 表示部
２５ 通信部
２６ バス
３ 診断用コンソール（動態解析装置）
３１ 制御部
３２ 記憶部
３３ 操作部
３４ 表示部
３５ 通信部
３６ バス
Ｍ マーカー
Ｒ 測定対象領域
Ｒ１ 小領域
Ｒ２ 他の領域
Ｒｖ 血管領域
Ｒｎ 非血管領域
Ｖ 血管

【図1】

【図2】

【図3】

【図4】

【図5】

【図6】

【図7】

【図8】

【図9】