ホーム > 特許ランキング > コニカミノルタ株式会社
知財求人 - 知財ポータルサイト「IP Force」
(19)【発行国】日本国特許庁(JP)
(12)【公報種別】特許公報(B1)
(11)【特許番号】
(24)【登録日】2024-01-22
(45)【発行日】2024-01-30
(54)【発明の名称】放射線画像解析装置及びプログラム
(51)【国際特許分類】
A61B 6/00 20240101AFI20240123BHJP
【ＦＩ】
A61B6/00 330A
A61B6/00 350B
A61B6/00 350C
【請求項の数】
8
(21)【出願番号】P 2023072945
(22)【出願日】2023-04-27
【審査請求日】2023-08-10
【新規性喪失の例外の表示】特許法第３０条第２項適用 令和４年９月２２日第２回日本びまん性肺疾患研究会プログラム抄録集において発表
(73)【特許権者】
【識別番号】000001270
【氏名又は名称】コニカミノルタ株式会社
(74)【代理人】
【識別番号】110001254
【氏名又は名称】弁理士法人光陽国際特許事務所
(72)【発明者】
【氏名】角森 昭教
(72)【発明者】
【氏名】野間 惠之
(72)【発明者】
【氏名】上山 維晋
【審査官】亀澤 智博
(56)【参考文献】
【文献】特開２０２０－１５１２３２（ＪＰ，Ａ）
【文献】特開２０２１－１４６０１２（ＪＰ，Ａ）
【文献】国際公開第２０２３／０６３３１８（ＷＯ，Ａ１）
【文献】特表２０２１－５０５２５８（ＪＰ，Ａ）
【文献】特開２０１９－１８７８６２（ＪＰ，Ａ）
【文献】国際公開第２０１７／１８３９５８（ＷＯ，Ａ１）
【文献】特開２０１６－２１４７２５（ＪＰ，Ａ）
【文献】国際公開第２００７／０４６２２０（ＷＯ，Ａ１）
(58)【調査した分野】(Int.Cl.，ＤＢ名)
Ａ６１Ｂ ６／００ － ６／１４
Ａ６１Ｂ ５／０８ － ５／０９７
Ｇ０６Ｔ １／００ ， ７／００
Ｇ１６Ｈ ３０／００ －３０／４０
(57)【特許請求の範囲】
【請求項1】
被検者の胸部の第１の放射線動態画像から抽出された第１の特徴量と、前記被検者に対する、前記第１の放射線動態画像が撮影された時期よりも過去の呼吸機能検査の検査結果情報とに基づいて、前記第１の放射線動態画像の撮影と同時期の前記被検者の呼吸機能を推定する推定手段、を備える放射線画像解析装置。
【請求項2】
前記第１の特徴量は、呼吸に伴って動く特定の構造物の領域の特徴量である、請求項１に記載の放射線画像解析装置。
【請求項3】
前記第１の特徴量は、肺野領域の面積、横隔膜の変位、気管径変化量又は肺野領域内の濃度変化率である、請求項１に記載の放射線画像解析装置。
【請求項4】
前記過去の呼吸機能検査は、前記被検者の初診時又は最新の呼吸機能検査である、請求項１に記載の放射線画像解析装置。
【請求項5】
前記呼吸機能検査は、ＦＶＣ又はＦＥＶ１．０を測定する検査を含み、前記推定手段は、前記呼吸機能としてＦＶＣ又はＦＥＶ１．０を推定する、請求項１に記載の放射線画像解析装置。
【請求項6】
前記推定手段は、さらに、前記呼吸機能検査と同時期に撮影された前記被検者の胸部の第２の放射線動態画像から抽出された第２の特徴量に基づいて、前記被検者の呼吸機能を推定する、請求項１に記載の放射線画像解析装置。
【請求項7】
前記推定手段は、さらに、前記被検者の属性情報に基づいて、前記被検者の呼吸機能を推定する、請求項１に記載の放射線画像解析装置。
【請求項8】
コンピューターを、被検者の胸部の第１の放射線動態画像から抽出された第１の特徴量と、前記被検者に対する、前記第１の放射線動態画像が撮影された時期よりも過去の呼吸機能検査の検査結果情報とに基づいて、前記第１の放射線動態画像の撮影と同時期の前記被検者の呼吸機能を推定する推定手段、
として機能させるためのプログラム。
【発明の詳細な説明】
【技術分野】
【0001】
本発明は、放射線画像解析装置及びプログラムに関する。
【背景技術】
【0002】
従来、呼吸機能を測定する方法としてスパイロメトリーが用いられている。しかし、スパイロメトリーは、患者への負担が大きい。特に、呼吸機能を経過観察するにあたり、スパイロメトリーを複数回行うとなると、患者への負担がより一層大きくなる。また、経過観察中に病状が悪化した場合には、スパイロメトリーでの呼吸機能検査が困難となる。
【0003】
スパイロメトリーに代わり、胸部の動態撮影で取得した動態画像から、呼吸機能を推定する方法が検討されている（例えば、特許文献１参照）。動態撮影は、スパイロメトリーに比べて侵襲性が低く、患者の負担が少ない。
【先行技術文献】
【特許文献】
【0004】
【文献】特開２０１９－１８７８６２号公報
【発明の概要】
【発明が解決しようとする課題】
【0005】
しかし、経過観察にあたり、例えば、術前術後などの呼吸機能の変化をとらえるには、呼吸機能をより正確に推定することが必要である。また、スパイロメトリーから動態撮影へ検査を置き換えるにあたり、検査の信頼性の観点から、推定精度を高める必要がある。
【0006】
本発明は、動態撮影で取得した動態画像から呼吸機能を推定するにあたり、推定精度を向上させることを課題とする。
【課題を解決するための手段】
【0007】
上記課題を解決するため、本発明の放射線画像解析装置は、
被検者の胸部の第１の放射線動態画像から抽出された第１の特徴量と、前記被検者に対する、前記第１の放射線動態画像が撮影された時期よりも過去の呼吸機能検査の検査結果情報とに基づいて、前記第１の放射線動態画像の撮影と同時期の前記被検者の呼吸機能を推定する推定手段、
を備える。
被検者の胸部の第１の放射線動態画像から抽出された第１の特徴量と、前記被検者に対する、前記第１の放射線動態画像が撮影された時期よりも過去の呼吸機能検査の検査結果情報とに基づいて、前記第１の放射線動態画像の撮影と同時期の前記被検者の呼吸機能を推定する推定手段、
を備える。
【0008】
また、本発明のプログラムは、
被検者の胸部の第１の放射線動態画像から抽出された第１の特徴量と、前記被検者に対する、前記第１の放射線動態画像が撮影された時期よりも過去の呼吸機能検査の検査結果情報とに基づいて、前記第１の放射線動態画像の撮影と同時期の前記被検者の呼吸機能を推定する推定手段、
として機能させる。
被検者の胸部の第１の放射線動態画像から抽出された第１の特徴量と、前記被検者に対する、前記第１の放射線動態画像が撮影された時期よりも過去の呼吸機能検査の検査結果情報とに基づいて、前記第１の放射線動態画像の撮影と同時期の前記被検者の呼吸機能を推定する推定手段、
として機能させる。
【発明の効果】
【0009】
本発明によれば、動態撮影で取得した動態画像から呼吸機能を推定するにあたり、推定精度を向上させることができる。
【図面の簡単な説明】
【0010】
【発明を実施するための形態】
【0011】
以下、図面を参照して本発明の実施形態を詳細に説明する。ただし、発明の範囲は、図示例に限定されない。
【0012】
〈放射線画像解析システム１００の構成〉
まず、本発明の実施形態の構成を説明する。
図１に、本実施形態における放射線画像解析システム１００の全体構成例を示す。
図１に示すように、放射線画像解析システム１００は、撮影装置１と、撮影用コンソール２と、解析装置３（放射線画像解析装置）と、を備えている。撮影装置１と撮影用コンソール２とは、通信ケーブル等により接続されている。撮影用コンソール２と解析装置３とは、ＬＡＮ（Local Area Network）等の通信ネットワークＮＴを介して接続されている。通信ネットワークＮＴには、図示しないＲＩＳ（Radiology Information Systems）、ＨＩＳ（Hospital Information Systems）、電子カルテシステム、ＰＡＣＳ（Picture Archiving and Communication System）等が接続されていることとしてもよい。放射線画像解析システム１００を構成する各装置は、ＤＩＣＯＭ（Digital Image and Communications in Medicine）規格に準じており、各装置間の通信は、ＤＩＣＯＭに則って行われる。
まず、本発明の実施形態の構成を説明する。
図１に、本実施形態における放射線画像解析システム１００の全体構成例を示す。
図１に示すように、放射線画像解析システム１００は、撮影装置１と、撮影用コンソール２と、解析装置３（放射線画像解析装置）と、を備えている。撮影装置１と撮影用コンソール２とは、通信ケーブル等により接続されている。撮影用コンソール２と解析装置３とは、ＬＡＮ（Local Area Network）等の通信ネットワークＮＴを介して接続されている。通信ネットワークＮＴには、図示しないＲＩＳ（Radiology Information Systems）、ＨＩＳ（Hospital Information Systems）、電子カルテシステム、ＰＡＣＳ（Picture Archiving and Communication System）等が接続されていることとしてもよい。放射線画像解析システム１００を構成する各装置は、ＤＩＣＯＭ（Digital Image and Communications in Medicine）規格に準じており、各装置間の通信は、ＤＩＣＯＭに則って行われる。
【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】
本実施形態において、制御部３１は、呼吸機能検査の検査結果情報として、FVC（努力肺活量）及びFEV1.0（一秒量）の検査結果情報を取得する。多くの場合、初診でスパイロメトリーなどの呼吸機能検査及び動態撮影を実施するため、制御部３１は、初診の呼吸機能検査の検査結果情報を過去の呼吸機能検査の検査結果情報として取得することができる。また、過去に複数回の呼吸機能検査が実施されている場合、制御部３１は、最新の呼吸機能検査の検査結果情報を取得することが好ましい。最新の呼吸機能検査における呼吸状態は、患者の現在の呼吸状態と最も近いと考えられるためである。
【0044】
次いで、制御部３１は、同一患者の過去の胸部の動態画像から抽出された注目領域の特徴量を取得する（ステップＳ１３）。
ステップＳ１３において、制御部３１は、ステップＳ１２で検査結果情報が取得された過去の呼吸機能検査と同時期に撮影された動態画像を記憶部３２から検索して取得する。一般的に、呼吸機能検査と動態撮影による検査は同日に行われる。そこで、制御部３１は、ステップＳ１２で検査結果情報が取得された過去の呼吸機能検査と同日に撮影された動態画像を記憶部３２から取得する。あるいは、制御部３１は、ステップＳ１２で検査結果情報が取得された過去の呼吸機能検査と同日に撮影された動態画像を図示しないＰＡＣＳから取得することとしてもよい。そして、制御部３１は、取得した動態画像から注目領域の特徴量を抽出する。
なお、制御部３１は、動態画像から抽出した特徴量の情報を動態画像に対応付けて記憶部３２に記憶しておくこととしてもよい。そして、制御部３１は、ステップＳ１２で検査結果情報が取得された過去の呼吸機能検査と同時期に撮影された動態画像に注目領域の特徴量の情報が対応付けられている場合には、その特徴量の情報を記憶部３２から読み出して取得することとしてもよい。
ステップＳ１３において、制御部３１は、ステップＳ１２で検査結果情報が取得された過去の呼吸機能検査と同時期に撮影された動態画像を記憶部３２から検索して取得する。一般的に、呼吸機能検査と動態撮影による検査は同日に行われる。そこで、制御部３１は、ステップＳ１２で検査結果情報が取得された過去の呼吸機能検査と同日に撮影された動態画像を記憶部３２から取得する。あるいは、制御部３１は、ステップＳ１２で検査結果情報が取得された過去の呼吸機能検査と同日に撮影された動態画像を図示しないＰＡＣＳから取得することとしてもよい。そして、制御部３１は、取得した動態画像から注目領域の特徴量を抽出する。
なお、制御部３１は、動態画像から抽出した特徴量の情報を動態画像に対応付けて記憶部３２に記憶しておくこととしてもよい。そして、制御部３１は、ステップＳ１２で検査結果情報が取得された過去の呼吸機能検査と同時期に撮影された動態画像に注目領域の特徴量の情報が対応付けられている場合には、その特徴量の情報を記憶部３２から読み出して取得することとしてもよい。
【0045】
次いで、制御部３１は、患者属性情報を取得する（ステップＳ１４）。
例えば、制御部３１は、通信部３５により図示しない電子カルテシステム等にアクセスして、患者属性情報を取得する。患者属性情報としては、例えば、性別、年齢、BMI（Body Mass Index）等が挙げられる。
例えば、制御部３１は、通信部３５により図示しない電子カルテシステム等にアクセスして、患者属性情報を取得する。患者属性情報としては、例えば、性別、年齢、BMI（Body Mass Index）等が挙げられる。
【0046】
次いで、制御部３１は、ステップＳ１１～Ｓ１４で取得した情報に基づいて、患者の呼吸機能を推定する（ステップＳ１５）。
例えば、制御部３１は、呼吸機能としてFVCを推定する。患者が男性である場合、制御部３１は、以下の式（１）によりFVCを推定する。患者が女性である場合、制御部３１は、以下の式（２）によりFVCを推定する。推定したFVCをFVCestとする。
ここで、Sinsは、最大肺野面積を示す。Sexpは、最小肺野面積を示す。Ageは、年齢を示す。「.clb」の識別子がついたパラメーター（Sins.clb、Sexp.clb、FVC.clb、FEV1.0.clb）は、過去の検査結果情報又は過去の撮影により取得された動態画像から抽出された特徴量であることを示す。「.clb」の識別子がついていないパラメーター（Sins、Sexp）は、今回の撮影により取得された動態画像から抽出された特徴量であることを示す。
すなわち、制御部３１は、患者が男性である場合は、式（１）のSinsに今回の動態画像から抽出した最大肺野面積を当てはめ、Sins.clbに過去の動態画像から抽出した最大肺野面積を当てはめ、Sexpに今回の動態画像から抽出した最小肺野面積を当てはめ、Sexp.clbに過去の動態画像から抽出した最小肺野面積を当てはめ、BMIに患者のBMIを当てはめ、Ageに患者の年齢を当てはめ、FVC.clbに過去の呼吸機能検査におけるFVCの値を当てはめ、FEV1.0.clbに過去の呼吸機能検査におけるFEV1.0を当てはめることにより、FVCest（FVCの推定値）を算出する。制御部３１は、患者が女性である場合は、式（２）のSinsに今回の動態画像から抽出した最大肺野面積を当てはめ、Sins.clbに過去の動態画像から抽出した最大肺野面積を当てはめ、Sexpに今回の動態画像から抽出した最小肺野面積を当てはめ、Sexp.clbに過去の動態画像から抽出した最小肺野面積を当てはめ、BMIに患者のBMIを当てはめ、Ageに患者の年齢を当てはめ、FVC.clbに過去の呼吸機能検査におけるFVCの値を当てはめ、FEV1.0.clbに過去の呼吸機能検査におけるFEV1.0を当てはめることにより、FVCest（FVCの推定値）を算出する。
例えば、制御部３１は、呼吸機能としてFVCを推定する。患者が男性である場合、制御部３１は、以下の式（１）によりFVCを推定する。患者が女性である場合、制御部３１は、以下の式（２）によりFVCを推定する。推定したFVCをFVCestとする。
【数1】
【数2】
すなわち、制御部３１は、患者が男性である場合は、式（１）のSinsに今回の動態画像から抽出した最大肺野面積を当てはめ、Sins.clbに過去の動態画像から抽出した最大肺野面積を当てはめ、Sexpに今回の動態画像から抽出した最小肺野面積を当てはめ、Sexp.clbに過去の動態画像から抽出した最小肺野面積を当てはめ、BMIに患者のBMIを当てはめ、Ageに患者の年齢を当てはめ、FVC.clbに過去の呼吸機能検査におけるFVCの値を当てはめ、FEV1.0.clbに過去の呼吸機能検査におけるFEV1.0を当てはめることにより、FVCest（FVCの推定値）を算出する。制御部３１は、患者が女性である場合は、式（２）のSinsに今回の動態画像から抽出した最大肺野面積を当てはめ、Sins.clbに過去の動態画像から抽出した最大肺野面積を当てはめ、Sexpに今回の動態画像から抽出した最小肺野面積を当てはめ、Sexp.clbに過去の動態画像から抽出した最小肺野面積を当てはめ、BMIに患者のBMIを当てはめ、Ageに患者の年齢を当てはめ、FVC.clbに過去の呼吸機能検査におけるFVCの値を当てはめ、FEV1.0.clbに過去の呼吸機能検査におけるFEV1.0を当てはめることにより、FVCest（FVCの推定値）を算出する。
【0047】
上記式（１）、式（２）は、FVCを目的変数、Sins 、Sins.clb、Sexp、Sexp.clb、BMI、Age、FVC.clb、FEV1.0.clbの各パラメーターを説明変数として、重回帰分析により生成された回帰式（重回帰モデル）である。
なお、例えば、FEV1.0を目的変数、上記と同様のパラメーター（Sins 、Sins.clb、Sexp、Sexp.clb、BMI、Age、FVC.clb、FEV1.0.clb）を説明変数として重回帰分析により生成された回帰式（重回帰モデル）を用いる等により、FEV1.0を推定することとしてもよい。
なお、例えば、FEV1.0を目的変数、上記と同様のパラメーター（Sins 、Sins.clb、Sexp、Sexp.clb、BMI、Age、FVC.clb、FEV1.0.clb）を説明変数として重回帰分析により生成された回帰式（重回帰モデル）を用いる等により、FEV1.0を推定することとしてもよい。
【0048】
次いで、制御部３１は、推定したFVCを表示部３４に表示するとともに、推定したFVCを今回の撮影により取得した動態画像に対応付けて記憶部３２に記憶させ（ステップＳ１６）、呼吸機能推定処理を終了する。
【0049】
図４は、上記式（１）又は式（２）により推定されたFVCest（推定FVC）と、FVCestの推定と同時期に同一患者について呼吸機能検査で測定されたFVC（実測FVC）と、の関係を示すグラフである。図４に示すグラフの縦軸は推定FVC（FVCest）、横軸は実測FCVである。図４に示すように、推定FVCは、実測FVCと非常に高い相関があることがわかる。すなわち、動態画像から抽出された注目領域の特徴量と、同一患者に対する過去の呼吸機能検査により取得した検査結果情報と、同一患者に対する過去の動態撮影により得られた動態画像から抽出した注目領域の特徴量の情報と、患者属性情報と、に基づいて、非常に精度よく呼吸機能を推定することができることがわかる。
【0050】
ここで、従来、動態画像から抽出された肺野面積などの呼吸に連動する注目領域の特徴量と呼吸機能とが関連することが示されている（例えば、参考文献１：N.Ohkura et al., “Chest Dynamic-Ventilatory Digital Radiography in Chronic Obstructive or Restrictive Lung Disease” International Journal of Chronic Obstructive Pulmonary Disease 2021:16 1393-1399, 参考文献２：M.Ueyama et al., ”Prediction of forced vital capacity with dynamic chest radiography in interstitial lung disease”, European Journal of Radiology, July21, 2021、特許文献１参照）。そこで、動態画像から得られた特徴量に基づいて呼吸機能を推定することが提案されている。しかし、呼吸機能は個人差が大きく、動態画像の注目領域の特徴量のみから呼吸機能を推定すると呼吸機能の個人差が反映されず、十分な推定精度が得られない場合があると考えられる。そこで、その患者の過去の呼吸機能検査結果を呼吸機能の推定に用いるパラメーターに加えることで、呼吸機能の個人差を呼吸機能の推定に反映させることが可能となり、従来よりも呼吸機能の推定精度を向上させることができると考えられる。
また、患者の過去の呼吸機能検査の検査結果情報だけでなく、さらに、この過去の呼吸機能検査と同時期に撮影された過去の動態画像における呼吸に連動する注目領域の特徴量の情報を加えることで、さらに呼吸機能の推定精度を向上させることができると考えられる。
また、年齢、性別、体格などは既に呼吸機能に影響することが広く知られている。よって、これらの患者属性情報を呼吸機能の推定に用いるパラメーターに加えることで、より高精度に呼吸機能を推定できると考えられる。
また、患者の過去の呼吸機能検査の検査結果情報だけでなく、さらに、この過去の呼吸機能検査と同時期に撮影された過去の動態画像における呼吸に連動する注目領域の特徴量の情報を加えることで、さらに呼吸機能の推定精度を向上させることができると考えられる。
また、年齢、性別、体格などは既に呼吸機能に影響することが広く知られている。よって、これらの患者属性情報を呼吸機能の推定に用いるパラメーターに加えることで、より高精度に呼吸機能を推定できると考えられる。
【0051】
（変形例１）
上記実施形態では、注目領域の特徴量を最大肺野面積及び最小肺野面積とした場合の、呼吸機能の推定に用いる回帰式の例を示した。しかし、注目領域の特徴量は、呼吸に伴って動く構造物の特徴量であれば、最大肺野面積及び最小肺野面積に限定されない。
上記実施形態では、注目領域の特徴量を最大肺野面積及び最小肺野面積とした場合の、呼吸機能の推定に用いる回帰式の例を示した。しかし、注目領域の特徴量は、呼吸に伴って動く構造物の特徴量であれば、最大肺野面積及び最小肺野面積に限定されない。
【0052】
例えば、注目領域を横隔膜、特徴量を呼吸による横隔膜の変位（Exc）として、下記の式（３）の回帰式を用いてFVCestを推定することとしてもよい。
ここで、ａ～ｊは、偏回帰係数（定数）を表す。Excは変位を示し、_Rは右肺、_Lは左肺を示す。「_clb」の識別子がついたパラメーター（Exc_R_clb、Exc_L_clb、FVC_clb、FEV1.0_clb）は、過去の呼吸機能検査の検査結果情報又は過去の動態画像から取得された特徴量であることを示す。「_clb」の識別子がついていないパラメーター（Exc_R、Exc_L）は、今回の撮影により取得された動態画像から得られた特徴量であることを示す。Ageは、年齢を示す。Genderは性別を示すパラメーターであり、例えば、男：1、女：0である。なお、式（３）、式（４）、式（５）におけるａ～ｈは異なる。
【数3】
【0053】
また、例えば、注目領域を気管、特徴量を呼吸による気管径変化量（Tra）として、下記の式（４）の回帰式を用いてFVCestを推定することとしてもよい。
ここで、ａ～ｈは、偏回帰係数（定数）を示す。Traは気管径変化量を示す。「_clb」の識別子がついたパラメーター（Tra_clb、FVC_clb、FEV1.0_clb）は、過去の呼吸機能検査の検査結果情報又は過去の動態画像から取得された特徴量であることを示す。「_clb」の識別子がついていないパラメーター（Tra）は、今回の撮影により取得された動態画像から得られた特徴量であることを示す。Ageは、年齢を示す。Genderは性別を示すパラメーターであり、例えば、男：1、女：0である。
【数4】
【0054】
また、例えば、注目領域を肺野、特徴量を呼吸による肺野領域内の濃度変化率（De）として、下記の式（５）の回帰式を用いてFVCestを推定することとしてもよい。なお、肺野領域の濃度変化率は、例えば、左右を含めた肺野領域の各画素値の濃度変化率の平均値である。
ここで、ａ～ｈは、偏回帰係数（定数）を示す。Deは肺野領域内の濃度変化率（平均濃度変化率）を示す。「_clb」の識別子がついたパラメーター（De_clb、FVC_clb、FEV1.0_clb）は、過去の呼吸機能検査の検査結果情報又は過去の動態画像から取得された特徴量であることを示す。「_clb」の識別子がついていないパラメーター（De）は、今回の撮影により取得された動態画像から得られた特徴量であることを示す。Ageは、年齢を示す。Genderは性別を示すパラメーターであり、例えば、男：1、女：0である。
【数5】
【0055】
なお、上記実施形態において、肺野面積は左右の肺野面積の合算としたが、左右の肺で別々に肺野面積を算出してそれぞれを重回帰式のパラメーターとしてもよい。同様に、上記変形例において、濃度変化率は左右を含めた肺野領域の各画素値の濃度変化率の平均値としたが、左右の肺で別々に濃度変化率を算出してそれぞれを重回帰式のパラメーターとしてもよい。また、例えば、FEV1.0を目的変数、式（３）～式（５）で用いたパラメーターを説明変数として重回帰分析により生成された回帰式を用いて、FEV1.0を推定することとしてもよい。
【0056】
（変形例２）
上記実施形態では、胸部正面の動態画像から抽出された注目領域の特徴量を、呼吸機能を推定するためのパラメーターとした。代わりに、胸部側面の動態画像から抽出された注目領域の特徴量を呼吸機能を推定するためのパラメーターとしてもよい。また、胸部正面と胸部側面の両方の動態画像から抽出された注目領域の特徴量を呼吸機能を推定するためのパラメーターとしてもよい。
上記実施形態では、胸部正面の動態画像から抽出された注目領域の特徴量を、呼吸機能を推定するためのパラメーターとした。代わりに、胸部側面の動態画像から抽出された注目領域の特徴量を呼吸機能を推定するためのパラメーターとしてもよい。また、胸部正面と胸部側面の両方の動態画像から抽出された注目領域の特徴量を呼吸機能を推定するためのパラメーターとしてもよい。
【0057】
（変形例３）
上記実施形態では、本発明を、FVC又はFEV1.0の推定に適用した場合を例にとり説明したが、例えば、VC（肺活量）、TLC（全肺気量）、RV（残気量）、FRC（機能的残気量）等の他の呼吸機能の予測に適用することとしてもよい。
上記実施形態では、本発明を、FVC又はFEV1.0の推定に適用した場合を例にとり説明したが、例えば、VC（肺活量）、TLC（全肺気量）、RV（残気量）、FRC（機能的残気量）等の他の呼吸機能の予測に適用することとしてもよい。
【0058】
以上説明したように、解析装置３の制御部３１は、被検者の胸部の動態画像から抽出された注目領域の特徴量と、その被検者に対する過去の呼吸機能検査の検査結果情報とに基づいて、その被検者の呼吸機能を推定する。
したがって、動態撮影で取得した動態画像から呼吸機能を推定するにあたり、推定精度を向上させることができる。
したがって、動態撮影で取得した動態画像から呼吸機能を推定するにあたり、推定精度を向上させることができる。
【0059】
また、制御部３１は、さらに、呼吸機能検査と同時期に撮影された被検者の胸部の動態画像から抽出された注目領域の特徴量に基づいて、その被検者の呼吸機能を推定することで、呼吸機能の推定精度を更に向上させることができる。
【0060】
また、制御部３１は、さらに、被検者の年齢、性別、体格などの属性情報に基づいて、被検者の呼吸機能を推定することで、呼吸機能の推定精度を更に向上させることができる。
【0061】
なお、本実施形態における記述は、本発明に係る好適な放射線画像解析装置及びプログラムの一例であり、これに限定されるものではない。
例えば、上記実施形態では、動態画像から注目領域の特徴量を算出する機能と、算出された特徴量を用いて呼吸機能を推定する機能が１つの装置に設けられている場合を例にとり説明した。しかし、動態画像から注目領域の特徴量を算出する機能は、呼吸機能を推定する機能を備える装置とは別の装置に設けられていることとしてもよい。
例えば、上記実施形態では、動態画像から注目領域の特徴量を算出する機能と、算出された特徴量を用いて呼吸機能を推定する機能が１つの装置に設けられている場合を例にとり説明した。しかし、動態画像から注目領域の特徴量を算出する機能は、呼吸機能を推定する機能を備える装置とは別の装置に設けられていることとしてもよい。
【0062】
また、例えば、上記の説明では、本発明に係るプログラムのコンピューター読み取り可能な媒体としてハードディスクや半導体の不揮発性メモリー等を使用した例を開示したが、この例に限定されない。その他のコンピューター読み取り可能な媒体として、ＣＤ-ＲＯＭ等の可搬型記録媒体を適用することが可能である。また、本発明に係るプログラムのデータを通信回線を介して提供する媒体として、キャリアウエーブ(搬送波)も適用される。
【0063】
その他、放射線画像解析システム１００を構成する各装置の細部構成及び細部動作に関しても、本発明の趣旨を逸脱することのない範囲で適宜変更可能である。
【符号の説明】
【0064】
１００ 放射線画像解析システム
１ 撮影装置
１１ 放射線源
１２ 放射線照射制御装置
１３ 放射線検出部
１４ 読取制御装置
２ 撮影用コンソール
２１ 制御部
２２ 記憶部
２３ 操作部
２４ 表示部
２５ 通信部
２６ バス
３ 解析装置
３１ 制御部
３２ 記憶部
３３ 操作部
３４ 表示部
３５ 通信部
３６ バス
１ 撮影装置
１１ 放射線源
１２ 放射線照射制御装置
１３ 放射線検出部
１４ 読取制御装置
２ 撮影用コンソール
２１ 制御部
２２ 記憶部
２３ 操作部
２４ 表示部
２５ 通信部
２６ バス
３ 解析装置
３１ 制御部
３２ 記憶部
３３ 操作部
３４ 表示部
３５ 通信部
３６ バス
【要約】
【課題】動態撮影で取得した動態画像から呼吸機能を推定するにあたり、推定精度を向上させる。
【解決手段】解析装置の制御部は、被検者の胸部の動態画像から抽出された注目領域の特徴量と、その被検者に対する過去の呼吸機能検査の検査結果情報とに基づいて、その被検者の呼吸機能を推定する。
【選択図】図３
【解決手段】解析装置の制御部は、被検者の胸部の動態画像から抽出された注目領域の特徴量と、その被検者に対する過去の呼吸機能検査の検査結果情報とに基づいて、その被検者の呼吸機能を推定する。
【選択図】図３
【図1】
【図2】
【図3】
【図4】