JP2021178111A - 動態画像解析装置、動態画像解析方法及びプログラム - Google Patents

Abstract

【課題】様々な観点から捉えた対象部位の状態の経時変化の情報を提供することのできる動態画像解析装置、動態画像解析方法及びプログラムを提供する。【解決手段】動態画像解析装置が、動態画像にタグ付けを行うとともに、動態画像を構成する複数のフレーム画像から特徴量を抽出し、この特徴量の時間方向における変化量を算出し、第１の動態画像と同じ対象部位を含む別撮影により得られた第２の動態画像が入力された場合に、タグから第２の動態画像と関連性があると判断される第１の動態画像である関連動態画像、及び第１の動態画像、第２の動態画像から抽出可能な特徴量を選定する制御部と、動態画像の一撮影を構成するフレーム画像間における経時変化を表示可能であるとともに、第２の動態画像及び第２の動態画像から抽出される特徴量と、選定された関連動態画像及び関連動態画像から抽出される特徴量との間の経時変化を表示可能である表示部と、を備える。【選択図】図５

Description

本発明は、動態画像解析装置、動態画像解析方法及びプログラムに関する。

近年、放射線撮影により、極短時間の連続撮影を行うことで生体における対象部位（検査および診断の対象となる部位）の動きを捉えた画像（動態画像という）を取得することが可能となっている。
例えば胸部を撮影した動態画像に対して画像解析を行うことにより、横隔膜や肺野辺縁の一撮影内における経時的変化を把握することのできる情報を取得することが可能である。

しかし、疾病の診断等を正確かつ効果的に行うためには、対象部位の状態について，術前術後の比較や治療開始からの経過観察等、ある程度の期間をおいての複数回の撮影における経時的な変化を知ることも重要である。
この点、例えば特許文献１には、骨密度分布を示す情報を測定日とともに記憶しておき、同じ患者の同一部位について測定日の異なる骨密度分布を比較し、骨密度の変化量を取得することが開示されている。

特許第５５２００７８号公報

しかしながら、特許文献１に記載の技術のように一般的な放射線撮影では、異なる時点に取得された瞬間的な撮影画像間での対象部位の状態の経時変化を把握することはできるが、一撮影内における経時変化を捉えることは困難である。
診断が進展するに連れて、診察内容も変わっていくものである。このため、診断等に有用な情報を医師等に提供するためには、こうした変化していく診察内容に応じて必要となる様々な観点から経時変化を捉えたいとの要望があるところ、既存の一般的な放射線撮影では、その画像の用途、検査目的等が確立され、診断に供される情報がある程度固定化してしまっているために、他の用途等に柔軟に適用し活用することが難しい。

本発明の課題は、様々な観点から捉えた対象部位の状態の経時変化の情報を提供することのできる動態画像解析装置、動態画像解析方法及びプログラムを提供することである。

上記課題を解決するため、請求項１に記載の発明の動態画像解析装置は、
生体における対象部位を含む被写体を放射線撮影することにより得られ一撮影のデータが複数のフレーム画像で構成される動態画像について、前記動態画像の類別に資する項目についてタグ付けを行うタグ設定部と、
前記動態画像を構成する１以上のフレーム画像から特徴量を抽出し、この特徴量の時間方向における変化量を算出する特徴変化量算出部と、
第１の前記動態画像と同じ前記対象部位を含む別撮影により得られた第２の前記動態画像が入力された場合に、前記タグ設定部により設定されたタグから前記第２の動態画像と関連性があると判断される前記第１の動態画像である関連動態画像、及び前記第１の動態画像、前記第２の動態画像から抽出可能な特徴量を選定するデータ選定部と、
前記動態画像の前記一撮影を構成する前記フレーム画像間における経時変化を表示可能であるとともに、前記第２の動態画像及び前記第２の動態画像から抽出される特徴量と、前記データ選定部によって選定された前記関連動態画像及び前記関連動態画像から抽出される特徴量との間の経時変化を表示可能である表示部と、
前記表示部の表示を制御する表示制御部と、
を備えることを特徴とする。

また、請求項９に記載の動態画像解析方法は、
生体における対象部位を含む被写体を放射線撮影することにより得られ一撮影のデータが複数のフレーム画像で構成される動態画像について、前記動態画像の類別に資する項目についてタグ付けを行うタグ設定工程と、
前記動態画像を構成する１以上のフレーム画像から特徴量を抽出し、この特徴量の時間方向における変化量を算出する特徴変化量算出工程と、
第１の前記動態画像と同じ前記対象部位を含む別撮影により得られた第２の前記動態画像が入力された場合に、前記タグ設定工程において設定されたタグから前記第２の動態画像と関連性があると判断される前記第１の動態画像である関連動態画像、及び前記第１の動態画像、前記第２の動態画像から抽出可能な特徴量を選定するデータ選定工程と、
表示部に、前記動態画像の前記一撮影を構成する前記フレーム画像間における経時変化を表示させるとともに、前記第２の動態画像及び前記第２の動態画像から抽出される特徴量と、前記データ選定工程において選定された前記関連動態画像及び前記関連動態画像から抽出される特徴量との間の経時変化を表示させる表示工程と、
を含むことを特徴とする。

また、請求項１０に記載のプログラムは、
コンピューターに、
生体における対象部位を含む被写体を放射線撮影することにより得られ一撮影のデータが複数のフレーム画像で構成される動態画像について、前記動態画像の類別に資する項目についてタグ付けを行うタグ設定機能と、
前記動態画像を構成する１以上のフレーム画像から特徴量を抽出し、この特徴量の時間方向における変化量を算出する特徴変化量算出機能と、
第１の前記動態画像と同じ前記対象部位を含む別撮影により得られた第２の前記動態画像が入力された場合に、前記タグ設定機能により設定されたタグから前記第２の動態画像と関連性があると判断される前記第１の動態画像である関連動態画像、及び前記第１の動態画像、前記第２の動態画像から抽出可能な特徴量を選定するデータ選定機能と、
表示部に、前記動態画像の前記一撮影を構成する前記フレーム画像間における経時変化を表示させるとともに、前記第２の動態画像及び前記第２の動態画像から抽出される特徴量と、前記データ選定機能によって選定された前記関連動態画像及び前記関連動態画像から抽出される特徴量との間の経時変化を表示させる表示機能と、
を実現させることを特徴とする。

本発明によれば、様々な観点から捉えた対象部位の状態の経時変化の情報を提供することができる。

本実施形態における動態画像解析装置を含む医用画像システムのシステム構成を示す要部構成図である。 撮影プロトコルの一例を示す説明図である。 本実施形態における動態画像解析方法の全体的な手順を示すフローチャートである。 本実施形態における動態画像解析方法における経時変化の表示処理の手順を示すフローチャートである。 （ａ）は、一撮影内における経時変化を模式的に示す説明図であり、（ｂ）は、一撮影内における経時変化の特徴量の一例を示す説明図であり、（ｃ）は、異なる時点に撮影された画像の経時変化を模式的に示す説明図である。 複数時点に撮影された各画像から抽出される特徴量の一例を示す説明図である。 元画像と元画像から抽出される特徴量の例を示す説明図である。 横隔膜の最大変位量の一例を示す表である。 肺野面積の変位の一例を示す表である。 複数時点に撮影された画像から抽出された肺野面積の変化率を時系列的に示すグラフの一例である。 複数時点に撮影された画像から抽出された肺野面積の最大／最小の変化を時系列的に示すグラフの一例である。

以下、図１から図１１を参照して、本発明にかかる動態画像解析装置、動態画像解析方法及びプログラムの一実施形態について説明する。
ただし、以下に述べる実施形態には、本発明を実施するために技術的に好ましい種々の限定が付されているが、本発明の技術的範囲を以下の実施形態及び図示例に限定するものではない。

本実施形態における動態画像解析装置は、医用画像システム内に設けられて、各種のモダリティーによって取得された画像（医用画像）のうち特に動態画像Ｄ（図５参照）の解析を行うものであり、
図１は、本実施形態における医用画像システム１００の全体構成を示す図である。
図１に示すように、医用画像システム１００は、撮影装置１と、撮影用コンソール２とが通信ケーブル等により接続され、撮影用コンソール２と、診断用コンソール３とがＬＡＮ（Local Area Network）等の通信ネットワークＮＴを介して接続されて構成されている。医用画像システム１００を構成する各装置は、ＤＩＣＯＭ（Digital Image and Communications in Medicine）規格に準じており、各装置間の通信は、ＤＩＣＯＭに則って行われる。
後述するように、本実施形態では診断用コンソール３が動態画像解析装置として機能する。

〔撮影装置１の構成〕
撮影装置１は、例えば、呼吸運動に伴う肺の膨張及び収縮の形態変化、心臓の拍動等の、生体の動態を撮影する撮影手段である。
図１に示すように、撮影装置１は、放射線源１１、放射線照射制御装置１２、放射線検出部１３、読取制御装置１４を含んでいる。

ここで動態撮影とは、生体における対象部位（例えば胸部の肺野周辺や心臓等）を含む被写体を放射線撮影することにより画像を得るものであり、被写体に対し、Ｘ線等の放射線をパルス状にして所定時間間隔で繰り返し照射するか（パルス照射）、もしくは、低線量率にして途切れなく継続して照射する(連続照射)ことで、複数の画像を取得することをいう。動態撮影により得られた一連の画像を動態画像Ｄと呼ぶ。動態画像Ｄは、一撮影のデータが複数の画像で構成されており、動態画像Ｄを構成する画像をフレーム画像（図５においてフレーム画像ｄ１〜ｄｎ）と呼ぶ。
なお、以下の実施形態では、パルス照射により動態撮影を行う場合を例にとり説明する。また、以下の実施形態では、診断の対象となる対象部位を胸部の肺野周辺（肺野や横隔膜）とした場合を例にとり説明するが、これに限定されるものではない。例えば、胸部であればこのほか心臓、肺血管、肋間筋、胸郭、胃部であれば腸類や食道等、その他例えば膝、肘、首、背骨等、様々な整形部位や各種組織等を対象部位としてもよい。

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

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

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

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

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

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

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

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

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

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

〔診断用コンソール３の構成〕
診断用コンソール３は、撮影用コンソール２から動態画像Ｄを取得し、取得した動態画像Ｄを解析して解析結果としての画像や各種のデータを生成し、生成した解析結果を表示して医師の診断を支援する情報を提供する動態画像解析装置である。
診断用コンソール３は、図１に示すように、制御部３１、記憶部３２、操作部３３、表示部３４、通信部３５を備えて構成され、各部はバス３６により接続されている。

制御部３１は、ＣＰＵ、ＲＡＭ等により構成される。制御部３１のＣＰＵは、操作部３３の操作に応じて、記憶部３２に記憶されているシステムプログラムや、各種処理プログラムを読み出してＲＡＭ内に展開し、展開されたプログラムに従って、後述する画像解析処理を始めとする各種処理を実行し、診断用コンソール３各部の動作を集中制御する。

本実施形態の制御部３１は、後述するように、動態画像Ｄの類別に資する項目についてタグ付けを行うタグ設定部として機能する。
また制御部３１は、動態画像Ｄを構成する１以上のフレーム画像から特徴量を抽出し、この特徴量の時間方向における変化量を算出する特徴変化量算出部として機能する。
また制御部３１は、過去に撮影されたある動態画像Ｄ（これを「第１の動態画像」という）と同じ対象部位を含む別撮影により得られた新たな動態画像Ｄ（これを「第２の動態画像」という）が入力された場合に、タグ設定部として機能する制御部３１により設定されたタグから第２の動態画像Ｄと関連性があると判断される第１の動態画像Ｄである関連動態画像（すなわち関連する過去の動態画像Ｄ）、及び第１の動態画像Ｄ、第２の動態画像Ｄから抽出可能な特徴量を選定するデータ選定部として機能する。
さらに制御部３１は、表示部３４の表示を制御する表示制御部としても機能する。
なお、制御部３１の上記各機能については、後に詳述する。

記憶部３２は、不揮発性の半導体メモリーやハードディスク等により構成される。記憶部３２は、制御部３１で画像解析処理等を実行するためのプログラムを始めとする各種プログラムやプログラムにより処理の実行に必要なパラメーター、或いは処理結果等のデータを記憶する。これらの各種プログラムは、読取可能なプログラムコードの形態で格納され、制御部３１は、当該プログラムコードに従った動作を逐次実行する。
また記憶部３２には、過去に取得された動態画像Ｄや動態画像Ｄを解析することで得られた特徴量、特徴量の変化量等の各種データが格納される。なお、動態画像Ｄや動態画像Ｄを解析することで得られる各種データは、診断用コンソール３内の記憶部３２に記憶されている場合に限定されない。例えば診断用コンソール３の外部に設けられた記憶部やサーバー等に記憶されていてもよい。

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

表示部３４は、ＬＣＤやＣＲＴ等のモニターにより構成され、制御部３１から入力される表示信号の指示に従って、各種表示を行う。
本実施形態の表示部３４は、動態画像Ｄの一撮影を構成するフレーム画像ｄ１〜ｄｎ間における経時変化を表示可能であるとともに、第２の動態画像Ｄ及び第２の動態画像Ｄから抽出される特徴量と、データ選定部としての制御部３１によって選定された関連動態画像及び関連動態画像から抽出される特徴量との間の経時変化等を表示可能となっている。
なお、表示部３４における具体的な表示については、後に詳述する。

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

〔動態画像解析装置としての診断用コンソールの動作〕
次に、上記動態画像解析装置としての診断用コンソール３の動作について説明する。

まず前提として、動態画像を取得するための撮影装置１、撮影用コンソール２による撮影動作について説明する。

撮影を行う際には、まず撮影実施者により撮影用コンソール２の操作部２３が操作され、被検者（被写体Ｍ）の患者情報（患者の氏名、身長、体重、年齢、性別等）や検査情報等の入力が行われる。
検査情報としては、例えば、撮影対象となる対象部位（例えば肺野や横隔膜等）や、対象部位の動作・組織変化の種類・方法（拡縮動作なのか、上下動なのか等）、撮影時のポジショニング（正面、側面等）、撮影時の撮影条件（管電圧、照射角度、撮影時間等）である。
次に、放射線照射条件が記憶部２２から読み出されて放射線照射制御装置１２に設定されるとともに、画像読取条件が記憶部２２から読み出されて読取制御装置１４に設定され、操作部２３の操作による放射線照射の指示が待機される状態となる。
ここで、撮影実施者は、被写体Ｍを放射線源１１と放射線検出部１３の間に配置してポジショニングを行う。

また、例えば肺野の拡縮動作や横隔膜の上下動を動態画像Ｄとしてとらえる場合には、呼吸状態下で撮影を行うため、被検者（被写体Ｍ）に楽にするように指示し、安静呼吸を促す。撮影準備が整った時点で、操作部２３を操作して放射線照射指示を入力する。
なお本実施形態では、後述するように、新たな撮影が行われた際に過去画像との比較評価が行われる。この比較評価においては経時変化が起こる条件を同一にしないと誤診に繋がるおそれがある。このため、検査時の対象部位や、体位・方向といったポジショニングの条件（例えば胸部側面を撮影する場合には被検者（被写体Ｍ）を所定の位置に立たせる等）、撮影時間等を一定に揃えることが好ましい他、呼吸等の撮影時の動作方法についても所定の撮影プロトコルに則って行うことで、撮影時の条件同一性を確保した状態で撮影を行うことが好ましい。
このように撮影時の条件を揃えておくことにより、画像解析について将来新たなアルゴリズムが登場した場合に、過去に撮影された画像であっても当該新たなアルゴリズムを用いて再処理を行うことができ、解析データを新たなアルゴリズムに基づいたものに最新化することが可能となる。

例えば、図３は、肺野の拡縮動作や横隔膜の上下動を動態画像としてとらえる際の撮影プロトコルの一例を模式的に示した説明図である。
図３に示す例では、自動の音声案内（図中「オートボイス」）によって息を吸うタイミング、息を止めるタイミング、息を吐くタイミングの誘導を行う。
そして、音声による呼吸タイミングの誘導に合わせてＸ線等の放射線をパルス照射することにより動態撮影を行う。
このようにすることで、被検者（被写体Ｍ）が大きく息を吸った後の息止めから息を吐いていく「呼気相」、息を吐き切った後の息止めから息を吸っていく「吸気相」という、呼吸における一サイクルを一撮影内に収めることができる。

操作部２３により放射線照射指示が入力されると、放射線照射制御装置１２及び読取制御装置１４に撮影開始指示が出力され、動態撮影が開始される。即ち、放射線照射制御装置１２に設定されたパルス間隔で放射線源１１により放射線が照射され、放射線検出部１３によりフレーム画像ｄ１〜ｄｎが取得される。

予め定められたフレーム数の撮影が終了すると、制御部２１により放射線照射制御装置１２及び読取制御装置１４に撮影終了の指示が出力され、撮影動作が停止される。撮影されるフレーム数は特に限定されないが、例えば図３に示すような呼吸による部位の変化を撮影する場合であれば、少なくとも１呼吸サイクルが撮影できる枚数からなる動態画像Ｄを取得する。

次に、本実施形態における動態画像Ｄの解析方法について図３及び図４等を参照しつつ説明する。
撮影装置１により動態画像Ｄが撮影されると、動態画像Ｄを構成するフレーム画像ｄ１〜ｄｎが順次撮影用コンソール２に入力され、撮影順を示す番号（フレーム番号）と対応付けられる。
そして、通信部３５を介して撮影用コンソール２から動態画像Ｄの一連のフレーム画像ｄ１〜ｄｎが診断用コンソール３において受信される。これによって動態画像Ｄ（第１の動態画像）が診断用コンソール３に取得され（ステップＳ１）、制御部３１と記憶部３２に記憶されているプログラムとの協働により図３及び図４に示す動態画像解析処理が実行される。

具体的にはまず、取得した動態画像Ｄについて、制御部３１がタグ設定部として、動態画像Ｄの類別に資する項目についてタグ付けを行う（ステップＳ２）。
ここで「類別に資する項目」は、記憶部３２等に蓄積されていく動態画像Ｄを後日検索することを可能とするためのインデックスとしての機能を果たすものである。
タグ設定部としての制御部によって設定されるタグは、例えば撮影時のポジショニング（正面か側面か、右方向か左方向か等）、対象部位（例えば肘や膝等の整形部位や胸部、腹部等）の動作・組織変化の種類や方法（例えば呼吸による肺野の拡縮、横隔膜の上下動、腸の蠕動、心臓や血管の拍動、脈動、関節等の曲げ伸ばし、部位の左右動等）や実際の組織の内容、各種の検査条件（例えば撮影時間、撮影間隔（フレームレート）、照射角度、線量、管球位置等の各種撮影条件等、撮影機器の型番等の情報）のうち少なくともいずれかの項目について設定されることが好ましい。
この他、例えば患者情報（例えば患者氏名や患者ＩＤ、患者の性別、年齢等）、撮影日、撮影時刻等、ＤＩＣＯＭ形式で画像データのヘッダ領域に書き込まれる付帯情報についてもタグ付けされる。

なお、タグ設定部としての制御部３１は、これら動態画像Ｄの「類別に資する項目」について、単一の項目についてタグ付けを行ってもよいし、複数の項目の組み合わせについてタグ付けを行ってもよい。
例えば、ある動態画像Ｄについて、胸部の画像であるとの単一のタグを付してもよいし、胸部の画像かつ肺野の拡縮を捉えた画像であるとの複数項目の組み合わせでタグを付してもよい。
タグが設定されると、制御部３１は動態画像Ｄをタグと対応付けて記憶部３２に記憶させる（ステップＳ３）。

そして、新たに動態画像Ｄ（第２の動態画像）が取得されると（ステップＳ４）、当該動態画像Ｄ（第２の動態画像）についても同様に、タグ設定部としての制御部３１が各種の「類別に資する項目」についてタグ付けを行い（ステップＳ５）、動態画像Ｄをタグと対応付けて記憶部３２に記憶させる（ステップＳ６）。

そして制御部３１は、一撮影で得られた動態画像Ｄ内の経時変化の表示を行うか否かを随時判断する（ステップＳ７）。なおこの判断は、新たな動態画像Ｄ（第２の動態画像）が取得されるごとに行ってもよいし、それ以外の所定のタイミングで行ってもよい。
一撮影内での経時変化の表示を行う場合（ステップＳ７；ＹＥＳ）には、動態画像Ｄを構成するフレーム画像ｄ１〜ｄｎ間の経時変化を表示部３４に表示させる（ステップＳ８）。

他方、一撮影内での経時変化の表示を行わない場合（ステップＳ７；ＮＯ）及び一撮影内での経時変化の表示を行う場合にはさらに、異なる時点で撮影された複数の動態画像Ｄ間における経時変化の表示を行うか否かを判断する（ステップＳ９）。
異なる時点で撮影された動態画像間での経時変化の表示を行う場合（ステップＳ９；ＹＥＳ）には、動態画像Ｄ間における経時変化の比較評価を行い、その結果等を表示部３４に表示させる（ステップＳ１０）。

動態画像Ｄの一撮影内での経時変化の表示と、異なる時点で撮影された動態画像Ｄ間での経時変化の比較評価及びその結果の表示とを行う場合の手順について、図４等を参照しつつ説明する。
この場合、制御部３１はデータ選定部として機能し、タグの情報に基づいて、記憶部３２内の動態画像Ｄ（過去に撮影された第１の動態画像）のうちから、第２の動態画像Ｄと関連性があると判断されるものを関連動態画像として選定・抽出する（ステップＳ２１）。
例えば、制御部３１は、第２の動態画像Ｄと同じ被検者（被写体Ｍ）に対する胸部の動態画像で、正面を向いたポジショニングで撮影されたもの、といったタグで検索をかけることにより、既に取得されている第１の動態画像Ｄの中から関連動態画像を選定・抽出する。
なお、検索に用いるタグは医師等が自分の求めるデータに応じて入力してもよいし、第２の動態画像に付されているタグに基づいて、自動的に選定するようにしてもよい。

なお、選定・抽出する範囲は関連性があると判断された動態画像Ｄ全てでもよいし、例えば直近２年以内に取得された画像等、範囲を限定してもよい。また本実施形態では、新たに取得された（すなわち最新の）動態画像Ｄを第２の動態画像とし、これに関連する動態画像Ｄを過去画像から選定する場合を例示するが、第２の動態画像Ｄは新たな動態画像に限定されず、例えばある手術の直後に撮影された動態画像Ｄがある場合に、これを第２の動態画像として、当該動態画像Ｄに関連する過去から直近までの（すなわち第２の動態画像の前後の）動態画像Ｄを関連動態画像として選定してもよい。

そしてデータ選定部としての制御部３１は、動態画像Ｄ（第２の動態画像及び第１の動態画像（関連動態画像））を構成する複数のフレーム画像ｄ１〜ｄｎ間の経時変化（一撮影内の経時変化）を解析し、比較評価可能なデータを得る（ステップＳ２２）。
例えば図５（ａ）は、呼吸動作を行いながら動態撮影を行うことで一撮影内に呼吸によって生じる横隔膜の上下動（経時的な変化）の様子を取得した例を説明する図である。
図５（ａ）に示すように、動態画像Ｄを構成するフレーム画像ｄ１〜ｄｎ間における横隔膜の上下動を解析してグラフに表すと、図５（ａ）の右端の図に示すような、一撮影内の経時変化を比較評価可能なグラフとなる。

動態画像Ｄからどのような比較評価可能なデータを解析・抽出できるかは、動態画像Ｄに適用されるアルゴリズムによって異なる。また本実施形態では、前述のように、撮影条件等を揃えて撮影を行うことにより、過去の画像に新たなアルゴリズムを適用することで、抽出可能な特徴量のデータをアップデートすることができるように構成されている。
図６は、動態画像に適用されるアルゴリズムによって解析・抽出することのできる比較評価可能なデータの違いを示す説明図である。
図６に示す例は、２０１８年１２月時点に適用されていたアルゴリズムでは、動態画像Ｄから強調画像Ｄｈ、画像信号値のデータ（信号値グラフＧａ）、左右の横隔膜の位置（高さ位置）を追跡したデータ（横隔膜追跡グラフＧｂ）を解析・抽出可能であったことを示している。これに対して、２０１９年１２月時点のアルゴリズムを動態画像Ｄに適用すると、強調画像Ｄｈ、左右の横隔膜の位置（高さ位置）を追跡したデータ（横隔膜追跡グラフＧｂ）、肺野辺縁検出データＤｏｌを解析・抽出することができるようになっている。

このため、２０１９年１２月時点のアルゴリズムが適用されるようになる前に取得され、記憶されていた動態画像Ｄについても、画像取得時とは異なる２０１９年１２月時点のアルゴリズムを適用することにより、例えば画像取得時には得られなかった肺野辺縁検出データＤｏｌを過去に取得された動態画像Ｄから抽出することが可能となる。
これにより、関連動態画像として選定した画像が古いアルゴリズムの下で解析されたものであって、新たな動態画像（第２の動態画像）と共通する比較評価可能なデータがない場合でも、新たな動態画像と同じアルゴリズムを関連動態画像に適用することで、共通する比較評価可能なデータを得ることが可能となる。
このように、本実施形態では、新たなアルゴリズムによってデータをアップデートすることにより、過去の動態画像Ｄのデータも有効に活用することができる。

各動態画像Ｄの一撮影内の経時変化を解析すると、特徴変化量算出部としての制御部３１は、動態画像Ｄ（第２の動態画像及び第１の動態画像（関連動態画像））を構成する複数のフレーム画像ｄ１〜ｄｎから動態画像Ｄ同士を比較評価することが可能な特徴量を抽出し、各動態画像Ｄをサマライズする（ステップＳ２３）。
本実施形態では、動態画像Ｄ内の経時変化に関する各種のデータについて、データの収集、保存、表示等、あらゆる場面においてサマライズ（集約）することが可能となっている。

例えば図７に示すように、ある動態画像Ｄから横隔膜追跡グラフＧｂ及び肺野辺縁検出データＤｏｌを解析・抽出することができた場合、左右の横隔膜の最大移動量（最大の移動地点を図中破線丸印で示す）を抜き出すことで、複数のフレーム画像ｄ１〜ｄｎからなる動態画像Ｄに関するデータを１つの値に集約することができる。
すなわち、複数のフレーム画像ｄ１〜ｄｎ間で経時的に変化しばらつきのある数値を、図８に示すように、右の横隔膜の最大移動量３２ｍｍ、左の横隔膜の最大移動量４９ｍｍという単一の数値にまとめることができる。
同様に、肺野辺縁検出データＤｏｌについて、肺野辺縁の検出座標から最大／最小面積を抜き出すことで、複数のフレーム画像ｄ１〜ｄｎからなる動態画像Ｄに関するデータを１つの値に集約することができる。また肺野の面積は経時変化するため、面積の最大／最小という変化をさらに定量化して観察するということもできる。

例えば図５（ａ）に示すような、一撮影内における横隔膜の上下動の経時変化からは、何に着目するかによって様々な特徴量を抽出することができる。
例えば、グラフの始点（左端）から終点（右端）までの振幅の平均値をとっていく手法もある。図５（ｂ）に示す例では、横隔膜が最も高い位置に来た時の値（最高位ｈｐ）と横隔膜が最も低い位置に来た時の値（最低位ｌｐ）とを特徴量として抽出する。このように一撮影内における最高位ｈｐと最低位ｌｐとを抽出することで、複数のフレーム画像ｄ１〜ｄｎからなる動態画像Ｄを他の動態画像Ｄと比較評価することが可能な１つの値に集約（サマライズ）する。

一撮影内の動態画像Ｄの経時変化をすべてデータとして抽出して、これをそのまま動態画像Ｄ同士の比較に用いると、比較対象が明確とならず、動態画像Ｄ同士を正しく比較評価することができない。
この点、各動態画像Ｄを着目すべきポイントを絞ってサマライズ（集約）することで、意味のある比較評価を行うことができるようになる。
すなわち、各動態画像Ｄ（第２の動態画像及び第１の動態画像（関連動態画像））が１つの値に集約（サマライズ）されることによって、特徴変化量算出部としての制御部３１は、容易に各値の変化量を算出する（ステップＳ２４）ことが可能となる。
特徴変化量算出部としての制御部３１は、１以上のフレーム画像ｄ１〜ｄｎ間の時間方向における変化に基づいて各値の変化量を算出する。
なお、制御部３１は、動態画像Ｄを構成する１以上のフレーム画像ｄ１〜ｄｎから変化量を算出すればよく、必ずしも動態画像Ｄに含まれるデータがサマライズ（集約）されている場合に限定されない。

例えば図５（ｃ）に示す例では、異なる時点（図５（ｃ）では、右から左に向かって順に最新のデータ、１か月前、２か月前、３か月前、４か月前）に撮影された複数の動態画像Ｄ同士を１つのグラフ等の上に比較評価可能に表現している。
同様に、肺野辺縁検出データＤｏｌから、肺野辺縁の検出座標に基づき、図９に示すように、右の肺野の最大／最小面積、左の肺野の最大／最小面積、左右の肺野を合計した場合の最大／最小面積、右の肺野の面積変化率、左の肺野の面積変化率、左右の肺野の面積変化率を特徴量として得ることができる。

肺野辺縁検出データＤｏｌから抽出されるこれらの値の経時的な変化を動態画像Ｄ同士で比較評価することで、図１０、図１１のようなグラフを得ることができる。
このうち図１０は、肺野面積の変化率の、動態画像Ｄ間での経時的な変化を示すグラフであり、グラフの左から右に行くにつれて新しいデータとなっている。図１０では、２０１８年３月１日から２０１９年の７月２６日までの約１６か月の間の肺野面積の変化率の推移を見ることができる。
また、図１１は、肺野面積の最大／最小値の、動態画像Ｄ間での経時的な変化を示すグラフであり、グラフの左から右に行くにつれて新しいデータとなっている。図１１では、２０１８年３月１日から２０１９年の７月２６日までの約１６か月の間の肺野面積の最大／最小値の推移を見ることができる。

制御部３１は表示制御部として、表示部３４の表示を制御し、一撮影内の経時変化を表示部３４に表示させる（ステップＳ２５）。
具体的には、動態画像Ｄの一撮影を構成する各フレーム画像ｄ１〜ｄｎや、図５（ａ）に示すグラフ、図６に示す強調画像Ｄｈ、信号値グラフＧａ、横隔膜追跡グラフＧｂ、肺野辺縁検出データＤｏｌ、数値を一覧とした図８、図９に示すような各種の表等、動態画像Ｄを解析することで得られる一撮影内の経時変化を示す各種データを表示部３４に表示させる。

さらに制御部３１は表示制御部として、表示部３４の表示を制御し、撮影時点の異なる複数の動態画像Ｄ間（すなわち、第２の動態画像とデータ選定部としての制御部３１によって選定された関連動態画像との間）における経時変化を表示部３４に表示させる（ステップＳ２６）。
具体的には、図５（ｃ）に示すグラフ、図１０、図１１に示すグラフのように、各動態画像Ｄ（第２の動態画像及び第１の動態画像（関連動態画像））についてサマライズ（集約）された数値同士を比較評価可能な状態で表示部３４に表示させる。また、第２の動態画像や関連動態画像の画像そのものを表示させてもよい。

また、表示部３４に表示される解析結果は医学的な評価に用いられるものであるため、医師等が見たときに医学的な観点から違和感を覚えることのないような表示を行う必要がある。
この点、例えば横隔膜の周辺には一般的に左側に心臓という大きな部位が配置される。このために、常に右側が低く左側が高い、というように画像上の絶対値としては必ず左右差が生じてしまう。しかし、左右の横隔膜がどの程度上下動するか、という機能的な観点からは、図７に示す横隔膜追跡グラフＧｂのように、左右の最低値を揃えて、横隔膜が機能的にどこまで上昇するかを表現することが好ましい。
そこで本実施形態では、表示制御部としての制御部３１は、対象部位の機能的な変化における基点を揃えて（すなわち、横隔膜の上下動を見る場合であれば左右の横隔膜の最低値をゼロ等に揃えて）特徴量の変化を表示部３４に表示させるようになっている。
また例えば、図１０、図１１のようなグラフにおいて、ある時点をクリックすると、当該時点のデータと、過去のデータや最新のデータとを比較可能に表示させてもよい。

このように、動態画像Ｄの一撮影内の経時変化と、撮影時点の異なる複数の動態画像Ｄ間における経時変化とを表示部３４に表示させることで、診断に役立つ資料として医師等に提供することができる。

以上説明したように、本実施形態における動態画像解析装置である診断用コンソール３によれば、動態画像Ｄにタグ付けを行うとともに、動態画像Ｄを構成する複数のフレーム画像ｄ１〜ｄｎから特徴量を抽出し、この特徴量の時間方向における変化量を算出し、第１の動態画像Ｄと同じ対象部位を含む別撮影により得られた第２の動態画像Ｄが入力された場合に、タグから第２の動態画像Ｄと関連性があると判断される第１の動態画像Ｄである関連動態画像、及び第１の動態画像Ｄ、第２の動態画像Ｄから抽出可能な特徴量を選定する制御部３１と、動態画像Ｄの一撮影を構成するフレーム画像ｄ１〜ｄｎ間における経時変化を表示可能であるとともに、第２の動態画像Ｄ及び第２の動態画像Ｄから抽出される特徴量と、選定された関連動態画像及び関連動態画像から抽出される特徴量との間の経時変化を表示可能である表示部３４と、を備える。
これにより、一撮影内において経時変化を捉えることのできる動態画像Ｄを用いて、動態画像Ｄ間における経時変化を捉えることが可能となる。
このため、例えば呼吸の一サイクル等、一撮影で得られるごく短い時間内での対象部位の経時的変化（例えば横隔膜の移動変化、肺野面積の拡縮変化等）と、こうした一撮影における経時変化を定量化して、複数時点での撮影で得られた動態画像Ｄ同士を比較したときの経時変化との両方を診断のための情報として提供することができ、両者を組み合わせた診断を可能とする。
例えば、肺等の活動状況の比較を行うためには、術前術後や治療開始から現在までの経過等、比較的長いスパンでの比較検討、状況把握が必要であるとともに、いつの時点ではどのような動き方をしていたのか、といったある時点での動作状況を適切に把握することも重要である。本実施形態における動態画像解析装置によれば、こうした両者を組み合わせた視点で患者の経過観察を行うことが可能となり、優れた臨床価値を実現する。

また本実施形態では、タグ設定部としての制御部３１が、撮影時のポジショニング、対象部位の動作・組織変化の種類、方法、撮影時の撮影条件のうち少なくともいずれかの項目について動態画像Ｄにタグ付けを行う。
このため、多角的な視点で画像を分類することができ、適切な関連動態画像を選定することが可能となる。

また本実施形態では、タグ設定部としての制御部３１が、動態画像Ｄの類別に資する項目について、単一で又は複数の項目の組み合わせでタグ付けを行う。
このように、単一のタグのみならず複数のタグを組み合わせるため、より適切な関連動態画像を選定することが可能となる。

また本実施形態では、特徴変化量算出部としての制御部３１が、１以上のフレーム画像ｄ１〜ｄｎから変化量を算出する。
このため、適切に定量化された動態画像Ｄから変化量を算出することができる。

また本実施形態では、特徴変化量算出部としての制御部３１が、１以上のフレーム画像間の時間方向における変化に基づいて変化量を算出する。
これにより、動態画像Ｄ内における経時的変化を適切に捉えることができる。

また本実施形態では、特徴変化量算出部としての制御部３１が、過去に取得された動態画像Ｄについて画像取得時とは異なるアルゴリズムを適用した上で特徴量を抽出し、この特徴量の時間方向における変化量を算出することが可能となっている。
これにより、新たに発見され適用されるアルゴリズムがあり、当該アルゴリズムを用いれば過去の画像取得時には抽出できなかった新たなデータ等がある場合に、当該アルゴリズムを過去に取得された動態画像Ｄに適用することでデータをアップデートすることができる。
これにより、過去に取得された動態画像Ｄについても有効に利用することができる。

また本実施形態では、表示制御部としての制御部３１が、動態画像Ｄの一撮影を構成するフレーム画像ｄ１〜ｄｎ間における経時変化を表示部３４に表示させるとともに、第２の動態画像及び第２の動態画像から抽出される特徴量と、データ選定部によって選定された関連動態画像及び関連動態画像から抽出される特徴量との間の経時変化を比較評価可能な状態で表示部３４表示させる。
これにより、一撮影で得られるごく短い時間内での対象部位の経時的変化と、複数時点での撮影で得られた動態画像Ｄ同士を比較したときの経時変化との両方を、診断を支援する情報として提供することができる。

また本実施形態では、表示制御部としての制御部３１が、対象部位の機能的な変化における基点を揃えて特徴量の変化を表示部３４表示させる。
このため、絶対値で見たときには医学的な見地から必ず差異が生じるような場合に、違和感なく分かりやすい表示を行うことができる。

なお、以上本発明の実施形態について説明したが、本発明は、かかる実施形態に限定されず、その要旨を逸脱しない範囲で、種々変形が可能であることは言うまでもない。

例えば、本実施形態では、医用画像システム１００内の診断用のコンソール３が動態画像解析装置として機能する場合を例示したが、動態画像解析装置はこれに限定されない。例えば医用画像システム１００外に、独立した動態画像解析ワークステーションが設けられ、各種ネットワーク等を介して医用画像システム１００から動態画像Ｄを取得する構成としてもよい。

また例えば、動態画像解析装置において解析処理が行われる動態画像Ｄは、医用画像システム１００内で取得されたものではなく、各種ネットワークを介して外部の装置から提供されたものであってもよい。
このように構成することで、１人の患者が複数の医療機関を受診しているような場合に、過去に受診した別の医療機関で撮影された動態画像Ｄも有効に活用することが可能となる。
なお、この場合、経時的変化を適切に比較することができるように、同一の撮影プロトコル等、共通のルールに則って撮影時の条件の同一性が確保された撮影された動態画像Ｄであることが好ましい。

また例えば、表示部３４において経時変化をグラフ等で表示する場合に、正常値を示しているデータと異常値となっているデータとで色を変える等により、変化の状況、傾向等を区別して表示させてもよい。
これにより、症状が改善傾向であるか否か等を視覚的に分かりやすく伝えることができる。

その他、動態画像解析装置（本実施形態では診断用コンソール３）を構成する各部の細部構成及び細部動作に関しても、本発明の趣旨を逸脱することのない範囲で適宜変更可能である。

１００ 医用画像システム
１ 撮影装置
１１ 放射線源
１２ 放射線照射制御装置
１３ 放射線検出部
１４ 読取制御装置
２ 撮影用コンソール
２１ 制御部
２２ 記憶部
２３ 操作部
２４ 表示部
２５ 通信部
２６ バス
３ 診断用コンソール
３１ 制御部
３２ 記憶部
３３ 操作部
３４ 表示部
３５ 通信部
３６ バス

Claims (10)

1. 生体における対象部位を含む被写体を放射線撮影することにより得られ一撮影のデータが複数のフレーム画像で構成される動態画像について、前記動態画像の類別に資する項目についてタグ付けを行うタグ設定部と、
   前記動態画像を構成する１以上のフレーム画像から特徴量を抽出し、この特徴量の時間方向における変化量を算出する特徴変化量算出部と、
   第１の前記動態画像と同じ前記対象部位を含む別撮影により得られた第２の前記動態画像が入力された場合に、前記タグ設定部により設定されたタグから前記第２の動態画像と関連性があると判断される前記第１の動態画像である関連動態画像、及び前記第１の動態画像、前記第２の動態画像から抽出可能な特徴量を選定するデータ選定部と、
   前記動態画像の前記一撮影を構成する前記フレーム画像間における経時変化を表示可能であるとともに、前記第２の動態画像及び前記第２の動態画像から抽出される特徴量と、前記データ選定部によって選定された前記関連動態画像及び前記関連動態画像から抽出される特徴量との間の経時変化を表示可能である表示部と、
   前記表示部の表示を制御する表示制御部と、
   を備えることを特徴とする動態画像解析装置。
2. 前記タグ設定部は、撮影時のポジショニング、前記対象部位の動作・組織変化の種類、方法、撮影時の撮影条件のうち少なくともいずれかの項目について前記動態画像にタグ付けを行うことを特徴とする請求項１に記載の動態画像解析装置。
3. 前記タグ設定部は、前記動態画像の類別に資する項目について、単一で又は複数の項目の組み合わせでタグ付けを行うことを特徴とする請求項１又は請求項２に記載の動態画像解析装置。
4. 前記特徴変化量算出部は、前記１以上のフレーム画像から変化量を算出することを特徴とする請求項１から請求項３のいずれか一項に記載の動態画像解析装置。
5. 前記特徴変化量算出部は、前記１以上のフレーム画像間の時間方向における変化に基づいて変化量を算出することを特徴とする請求項１から請求項３のいずれか一項に記載の動態画像解析装置。
6. 前記特徴変化量算出部は、過去に取得された前記動態画像について画像取得時とは異なるアルゴリズムを適用した上で前記特徴量を抽出し、この特徴量の時間方向における変化量を算出することが可能であることを特徴とする請求項１から請求項５のいずれか一項に記載の動態画像解析装置。
7. 前記表示制御部は、前記動態画像の前記一撮影を構成する前記フレーム画像間における経時変化を表示させるように前記表示部の表示を制御するとともに、前記第２の動態画像及び前記第２の動態画像から抽出される特徴量と、前記データ選定部によって選定された前記関連動態画像及び前記関連動態画像から抽出される特徴量との間の経時変化を比較評価可能な状態で表示させるように前記表示部の表示を制御することを特徴とする請求項１から請求項６のいずれか一項に記載の動態画像解析装置。
8. 前記表示制御部は、前記対象部位の機能的な変化における基点を揃えて前記特徴量の変化を表示させるように前記表示部の表示を制御することを特徴とする請求項１から請求項７のいずれか一項に記載の動態画像解析装置。
9. 生体における対象部位を含む被写体を放射線撮影することにより得られ一撮影のデータが複数のフレーム画像で構成される動態画像について、前記動態画像の類別に資する項目についてタグ付けを行うタグ設定工程と、
   前記動態画像を構成する１以上のフレーム画像から特徴量を抽出し、この特徴量の時間方向における変化量を算出する特徴変化量算出工程と、
   第１の前記動態画像と同じ前記対象部位を含む別撮影により得られた第２の前記動態画像が入力された場合に、前記タグ設定工程において設定されたタグから前記第２の動態画像と関連性があると判断される前記第１の動態画像である関連動態画像、及び前記第１の動態画像、前記第２の動態画像から抽出可能な特徴量を選定するデータ選定工程と、
   表示部に、前記動態画像の前記一撮影を構成する前記フレーム画像間における経時変化を表示させるとともに、前記第２の動態画像及び前記第２の動態画像から抽出される特徴量と、前記データ選定工程において選定された前記関連動態画像及び前記関連動態画像から抽出される特徴量との間の経時変化を表示させる表示工程と、
   を含むことを特徴とする動態画像解析方法。
10. コンピューターに、
    生体における対象部位を含む被写体を放射線撮影することにより得られ一撮影のデータが複数のフレーム画像で構成される動態画像について、前記動態画像の類別に資する項目についてタグ付けを行うタグ設定機能と、
    前記動態画像を構成する１以上のフレーム画像から特徴量を抽出し、この特徴量の時間方向における変化量を算出する特徴変化量算出機能と、
    第１の前記動態画像と同じ前記対象部位を含む別撮影により得られた第２の前記動態画像が入力された場合に、前記タグ設定機能により設定されたタグから前記第２の動態画像と関連性があると判断される前記第１の動態画像である関連動態画像、及び前記第１の動態画像、前記第２の動態画像から抽出可能な特徴量を選定するデータ選定機能と、
    表示部に、前記動態画像の前記一撮影を構成する前記フレーム画像間における経時変化を表示させるとともに、前記第２の動態画像及び前記第２の動態画像から抽出される特徴量と、前記データ選定機能によって選定された前記関連動態画像及び前記関連動態画像から抽出される特徴量との間の経時変化を表示させる表示機能と、
    を実現させることを特徴とするプログラム。

Images