(57)【要約】【課題】マーカを用いずに被写体を放射線撮影した動態画像において、対象構造物の境界が不鮮明な場合であっても、各フレーム画像における対象構造物の位置を精度良く計測できるようにする。【解決手段】診断用コンソール３によれば、制御部３１は、動態画像の複数のフレーム画像のそれぞれにおける所定範囲内の各位置について、横隔膜との類似度を示す評価値を算出し、算出された評価値に基づいて、複数のフレーム画像のそれぞれから横隔膜位置候補を一つ以上抽出し、複数のフレーム画像のそれぞれから抽出された横隔膜位置候補を一つずつ時系列に繋いで横隔膜が移動したルート候補としてルート保持部３２１に記憶させる。そして、ルート保持部３２１に記憶されているルート候補の中から一つのルート候補を横隔膜の移動ルートとして決定し、決定したルートに含まれる横隔膜位置候補を複数のフレーム画像のそれぞれにおける横隔膜の位置として決定する。【選択図】図３
2018068400.tif 000002
【特許請求の範囲】
【請求項１】
被写体に放射線を照射して動態撮影することにより得られた複数のフレーム画像における所定の構造物の位置を計測する動態画像処理装置であって、
前記複数のフレーム画像のそれぞれにおける所定範囲内の各位置について、前記構造物との類似度を示す評価値を算出する評価値算出部と、
前記評価値算出部により算出された評価値に基づいて、前記複数のフレーム画像のそれぞれから前記構造物の位置候補を一つ以上抽出し、少なくとも一つの前記フレーム画像からは前記構造物の位置候補を複数抽出する位置候補抽出部と、
前記複数のフレーム画像のそれぞれから抽出された前記構造物の位置候補を時系列に繋いで前記構造物が移動したルート候補として複数のルート候補を記憶する記憶部と、
前記記憶部に記憶されている複数のルート候補の中から一つのルート候補を前記構造物の移動ルートとして決定し、決定したルートに含まれる前記構造物の位置候補を前記複数のフレーム画像のそれぞれにおける前記構造物の位置として決定するルート決定部と、
を備える動態画像処理装置。
【請求項２】
前記位置候補抽出部は、前記複数のフレーム画像のそれぞれにおいて、前記評価値算出部により算出された評価値が極大となる位置を前記構造物の位置候補として抽出する請求項１に記載の動態画像処理装置。
【請求項３】
前記位置候補抽出部は、前記構造物の位置候補の抽出範囲を前記構造物の移動速度及び／又は加速度に基づいて制限する請求項２に記載の動態画像処理装置。
【請求項４】
前記位置候補抽出部は、前記構造物の移動方向が前記複数のフレーム画像のそれぞれを撮影したときの呼吸状態に合致するか否かに基づいて、抽出する前記構造物の位置候補を制限する請求項２又は３に記載の動態画像処理装置。
【請求項５】
前記ルート決定部は、前記記憶部に記憶されている各ルート候補に含まれる前記構造物の位置候補のそれぞれについて算出された評価値に基づいて、前記構造物の移動ルートを決定する請求項１〜４の何れか一項に記載の動態画像処理装置。
【請求項６】
前記ルート決定部は、前記ルート候補のそれぞれにおける、前記構造物の位置候補の評価値の合計値又は平均値、前記構造物の位置候補の評価値が所定の閾値以上又は所定の閾値未満のフレーム画像数又はフレーム画像数の割合、前記複数のルート候補間で同一のフレーム画像から算出された前記構造物の位置候補の評価値を比較したときの高さ順位の合計値又は平均値、又は、前記複数のルート候補間で同一のフレーム画像から算出された前記構造物の位置候補の評価値を比較したときの高さ順位が所定の閾値以上又は所定の閾値未満のフレーム画像数又はフレーム画像の割合に基づいて、前記構造物の移動ルートを決定する請求項５に記載の動態画像処理装置。
【請求項７】
前記ルート決定部は、前記記憶部に記憶されている各ルート候補における前記構造物の移動量及び／又は移動方向に基づいて、前記構造物の移動ルートを決定する請求項１〜６の何れか一項に記載の動態画像処理装置。
【請求項８】
前記構造物の位置を計測するための初期設定をユーザーが入力するための初期設定入力部を備える請求項１〜７の何れか一項に記載の動態画像処理装置。
【請求項９】
前記構造物の位置の計測結果を出力する出力部を備える請求項１〜８の何れか一項に記載の動態画像処理装置。
【請求項１０】
前記記憶部に記憶されている複数のルート候補をより少ない複数のルート候補に限定するルート限定部を備える請求項１〜９の何れか一項に記載の動態画像処理装置。
【請求項１１】
前記ルート限定部は、前記記憶部に記憶されている各ルート候補に含まれる前記構造物の位置候補のそれぞれについて算出された評価値に基づいて、前記ルート候補を限定する請求項１０に記載の動態画像処理装置。
【請求項１２】
前記ルート限定部は、前記評価値が所定の閾値以上の前記構造物の位置候補が存在した場合に、その位置候補を含むルート候補に限定する請求項１１に記載の動態画像処理装置。
【請求項１３】
前記ルート限定部は、前記記憶部に記憶されている各ルート候補における前記構造物の移動量及び／又は移動方向に基づいて、前記ルート候補を限定する請求項１０〜１２の何れか一項に記載の動態画像処理装置。
【請求項１４】
前記ルート限定部は、前記ルート候補のそれぞれにおける前記構造物の動きが前記動態撮影時の呼吸運動に対応しているか否かに基づいて、前記ルート候補を限定する請求項１０〜１３の何れか一項に記載の動態画像処理装置。
【請求項１５】
前記所定の構造物は、横隔膜、心壁、肋骨、又は血管である請求項１〜１４の何れか一項に記載の動態画像処理装置。
【請求項１６】
前記評価値算出部は、テンプレートマッチングを用いて前記評価値を算出する請求項１〜１５の何れか一項に記載の動態画像処理装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、動態画像処理装置に関する。
【背景技術】
【０００２】
近年、Ｘ線動画ＦＰＤ（Flat Panel Detector）の開発により、Ｘ線動態画像を撮影することが可能となった。ＣＴ（Computed Tomography）やＭＲＩ（Magnetic Resonance Imaging）は高価で設置できる医療機関も限られるが、Ｘ線動態撮影システムは比較的安価で多くの医療機関に導入できるものであり、医師は容易に体内の構造物の動きを観察し、診断することができるようになってきている。
【０００３】
ところで、医用画像から体内の所定の構造物の位置を抽出する手法として、テンプレートマッチングを用いた抽出処理が知られている。テンプレートマッチングは、抽出したい対象構造物の基準画像をテンプレート画像として用意しておき、テンプレート画像を抽出対象の画像の中で移動させながら、テンプレート画像と重なる部分との相関値を算出していく処理である。一般的に、算出された相関値が最大の箇所が対象構造物の位置として抽出される。
【０００４】
例えば、特許文献１には、３次元超音波動画像において、被写体の心壁の動きをテンプレートマッチングを用いて計測することが記載されている。
また、特許文献２には、２次元のＸ線透視画像において、体内の腫瘍付近に埋め込んだマーカの位置をテンプレートマッチングを用いて探索し、腫瘍とマーカの位置関係から腫瘍の位置を特定することが記載されている。
【先行技術文献】
【特許文献】
【０００５】
【特許文献１】特開２０１４−７６３３１号公報
【特許文献２】特許第５６６７４８９号公報
【発明の概要】
【発明が解決しようとする課題】
【０００６】
特許文献１では、上述のように、３次元画像にテンプレートマッチングを適用している。しかしながら、２次元のＸ線動態画像では、構造物同士の重なりが生じるため、対象の構造物の境界が曖昧になり、テンプレートマッチングを適用した場合、対象構造物の位置の相関値が低くなり精度よく位置を検出できない場合がある。図１１に、Ｘ線動態画像の例を示す。図１１の左図は横隔膜の境界が明瞭だが、右図では心臓と肋骨が横隔膜に重なり、境界が曖昧になっている。
【０００７】
また、特許文献２では、Ｘ線透視画像内の体内に埋め込んだ球形のマーカの位置・動きをテンプレートマッチングを用いて計測しているが、マーカを利用するとコスト・医療技術ともに大がかりなものとなり、利用現場が限定される。マーカを用いずに横隔膜等の臓器を対象構造物としてテンプレートマッチングを行うと、横隔膜等の臓器はマーカと比べて画像上に映りにくく、また時間とともに変形するため、相関値が最大の位置を抽出してしまうと高い精度で位置を検出することは困難である。
【０００８】
本発明の課題は、マーカを用いずに被写体を放射線撮影した動態画像において、対象構造物の境界が不鮮明な場合であっても、各フレーム画像における対象構造物の位置を精度良く計測できるようにすることである。
【課題を解決するための手段】
【０００９】
上記課題を解決するため、請求項１に記載の発明は、
被写体に放射線を照射して動態撮影することにより得られた複数のフレーム画像における所定の構造物の位置を計測する動態画像処理装置であって、
前記複数のフレーム画像のそれぞれにおける所定範囲内の各位置について、前記構造物との類似度を示す評価値を算出する評価値算出部と、
前記評価値算出部により算出された評価値に基づいて、前記複数のフレーム画像のそれぞれから前記構造物の位置候補を一つ以上抽出し、少なくとも一つの前記フレーム画像からは前記構造物の位置候補を複数抽出する位置候補抽出部と、
前記複数のフレーム画像のそれぞれから抽出された前記構造物の位置候補を時系列に繋いで前記構造物が移動したルート候補として複数のルート候補を記憶する記憶部と、
前記記憶部に記憶されている複数のルート候補の中から一つのルート候補を前記構造物の移動ルートとして決定し、決定したルートに含まれる前記構造物の位置候補を前記複数のフレーム画像のそれぞれにおける前記構造物の位置として決定するルート決定部と、
を備える。
【００１０】
請求項２に記載の発明は、請求項１に記載の発明において、
前記位置候補抽出部は、前記複数のフレーム画像のそれぞれにおいて、前記評価値算出部により算出された評価値が極大となる位置を前記構造物の位置候補として抽出する。
【００１１】
請求項３に記載の発明は、請求項２に記載の発明において、
前記位置候補抽出部は、前記構造物の位置候補の抽出範囲を前記構造物の移動速度及び／又は加速度に基づいて制限する。
【００１２】
請求項４に記載の発明は、請求項２又は３に記載の発明において、
前記位置候補抽出部は、前記構造物の移動方向が前記複数のフレーム画像のそれぞれを撮影したときの呼吸状態に合致するか否かに基づいて、抽出する前記構造物の位置候補を制限する。
【００１３】
請求項５に記載の発明は、請求項１〜４の何れか一項に記載の発明において、
前記ルート決定部は、前記記憶部に記憶されている各ルート候補に含まれる前記構造物の位置候補のそれぞれについて算出された評価値に基づいて、前記構造物の移動ルートを決定する。
【００１４】
請求項６に記載の発明は、請求項５に記載の発明において、
前記ルート決定部は、前記ルート候補のそれぞれにおける、前記構造物の位置候補の評価値の合計値又は平均値、前記構造物の位置候補の評価値が所定の閾値以上又は所定の閾値未満のフレーム画像数又はフレーム画像数の割合、前記複数のルート候補間で同一のフレーム画像から算出された前記構造物の位置候補の評価値を比較したときの高さ順位の合計値又は平均値、又は、前記複数のルート候補間で同一のフレーム画像から算出された前記構造物の位置候補の評価値を比較したときの高さ順位が所定の閾値以上又は所定の閾値未満のフレーム画像数又はフレーム画像の割合に基づいて、前記構造物の移動ルートを決定する。
【００１５】
請求項７に記載の発明は、請求項１〜６の何れか一項に記載の発明において、
前記ルート決定部は、前記記憶部に記憶されている各ルート候補における前記構造物の移動量及び／又は移動方向に基づいて、前記構造物の移動ルートを決定する。
【００１６】
請求項８に記載の発明は、請求項１〜７の何れか一項に記載の発明において、
前記構造物の位置を計測するための初期設定をユーザーが入力するための初期設定入力部を備える。
【００１７】
請求項９に記載の発明は、請求項１〜８の何れか一項に記載の発明において、
前記構造物の位置の計測結果を出力する出力部を備える。
【００１８】
請求項１０に記載の発明は、請求項１〜９の何れか一項に記載の発明において、
前記記憶部に記憶されている複数のルート候補をより少ない複数のルート候補に限定するルート限定部を備える。
【００１９】
請求項１１に記載の発明は、請求項１０に記載の発明において、
前記ルート限定部は、前記記憶部に記憶されている各ルート候補に含まれる前記構造物の位置候補のそれぞれについて算出された評価値に基づいて、前記ルート候補を限定する。
【００２０】
請求項１２に記載の発明は、請求項１１に記載の発明において、
前記ルート限定部は、前記評価値が所定の閾値以上の前記構造物の位置候補が存在した場合に、その位置候補を含むルート候補に限定する。
【００２１】
請求項１３に記載の発明は、請求項１０〜１２の何れか一項に記載の発明において、
前記ルート限定部は、前記記憶部に記憶されている各ルート候補における前記構造物の移動量及び／又は移動方向に基づいて、前記ルート候補を限定する。
【００２２】
請求項１４に記載の発明は、請求項１０〜１３の何れか一項に記載の発明において、
前記ルート限定部は、前記ルート候補のそれぞれにおける前記構造物の動きが前記動態撮影時の呼吸運動に対応しているか否かに基づいて、前記ルート候補を限定する。
【００２３】
請求項１５に記載の発明は、請求項１〜１４の何れか一項に記載の発明において、
前記所定の構造物は、横隔膜、心壁、肋骨、又は血管である。
【００２４】
請求項１６に記載の発明は、請求項１〜１５の何れか一項に記載の発明において、
前記評価値算出部は、テンプレートマッチングを用いて前記評価値を算出する。
【発明の効果】
【００２５】
本発明によれば、マーカを用いずに被写体を放射線撮影した動態画像において、対象構造物の境界が不鮮明な場合であっても、各フレーム画像における対象構造物の位置を精度良く計測することが可能となる。
【図面の簡単な説明】
【００２６】
【図１】本発明の実施形態における位置計測システムの全体構成を示す図である。
【図２】図１の撮影用コンソールの制御部により実行される撮影制御処理を示すフローチャートである。
【図３】図１の診断用コンソールの制御部により実行される位置計測処理を示すフローチャートである。
【図４】初期設定入力画面の一例を示す図である。
【図５】横隔膜位置候補の抽出手法を説明するための図である。
【図６】ルート保持部のデータ格納例を示す図である。
【図７】ルート候補の限定手法を説明するための図である。
【図８】計測結果画面の一例を示す図である。
【図９】結果修正画面の一例を示す図である。
【図１０】結果修正画面の他の例を示す図である。
【図１１】胸部の動態画像における横隔膜の境界を示す図である。
【発明を実施するための形態】
【００２７】
以下、図面を参照して本発明の実施形態を詳細に説明する。ただし、発明の範囲は、図示例に限定されない。
【００２８】
〔位置計測システム１００の構成〕
まず、構成を説明する。
図１に、本実施形態における位置計測システム１００の全体構成を示す。
図１に示すように、位置計測システム１００は、撮影装置１と、撮影用コンソール２とが通信ケーブル等により接続され、撮影用コンソール２と、診断用コンソール３とがＬＡＮ（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）等を備え、通信ネットワークＮＴに接続された各装置との間のデータ送受信を制御する。
【００４０】
〔診断用コンソール３の構成〕
診断用コンソール３は、撮影用コンソール２から送信された動態画像を解析することにより医師の診断を支援するための動態画像処理装置である。
診断用コンソール３は、図１に示すように、制御部３１、記憶部３２、操作部３３、表示部３４、通信部３５を備えて構成され、各部はバス３６により接続されている。
【００４１】
制御部３１は、ＣＰＵ、ＲＡＭ等により構成される。制御部３１のＣＰＵは、操作部３３の操作に応じて、記憶部３２に記憶されているシステムプログラムや、各種処理プログラムを読み出してＲＡＭ内に展開し、展開されたプログラムに従って、後述する位置計測処理を始めとする各種処理を実行し、診断用コンソール３各部の動作を集中制御する。  
【００４２】
記憶部３２は、不揮発性の半導体メモリーやハードディスク等により構成される。記憶部３２は、制御部３１で位置計測処理を実行するためのプログラムを始めとする各種プログラムやプログラムにより処理の実行に必要なパラメーター、或いは処理結果等のデータを記憶する。これらの各種プログラムは、読取可能なプログラムコードの形態で格納され、制御部３１は、当該プログラムコードに従った動作を逐次実行する。
【００４３】
また、記憶部３２は、撮影用コンソール２から受信された動態画像及び当該動態画像における位置計測結果を対応付けて記憶する。
また、記憶部３２には、位置計測処理におけるルート候補を記憶するためのルート保持部３２１が設けられている（図６参照）。ルート保持部３２１については詳細を後述する。
【００４４】
操作部３３は、カーソルキー、数字入力キー、及び各種機能キー等を備えたキーボードと、マウス等のポインティングデバイスを備えて構成され、ユーザーによるキーボードに対するキー操作やマウス操作により入力された指示信号を制御部３１に出力する。また、操作部３３は、表示部３４の表示画面にタッチパネルを備えても良く、この場合、手指やタッチペン等によりタッチパネルを介して入力された指示信号を制御部３１に出力する。
【００４５】
表示部３４は、ＬＣＤやＣＲＴ等のモニターにより構成され、制御部３１から入力される表示信号の指示に従って、各種表示を行う。
【００４６】
通信部３５は、ＬＡＮアダプターやモデムやＴＡ等を備え、通信ネットワークＮＴに接続された各装置との間のデータ送受信を制御する。
【００４７】
〔位置計測システム１００の動作〕
次に、上記位置計測システム１００における動作について説明する。
【００４８】
（撮影装置１、撮影用コンソール２の動作）
まず、撮影装置１、撮影用コンソール２による撮影動作について説明する。
図２に、撮影用コンソール２の制御部２１において実行される撮影制御処理を示す。撮影制御処理は、制御部２１と記憶部２２に記憶されているプログラムとの協働により実行される。
【００４９】
まず、撮影実施者により撮影用コンソール２の操作部２３が操作され、被検者（被写体Ｍ）の患者情報（患者の氏名、身長、体重、年齢、性別等）や検査情報（被写体部位（ここでは、胸部））の入力が行われる（ステップＳ１）。
【００５０】
次いで、放射線照射条件が記憶部２２から読み出されて放射線照射制御装置１２に設定されるとともに、画像読取条件が記憶部２２から読み出されて読取制御装置１４に設定される（ステップＳ２）。
【００５１】
次いで、操作部２３の操作による放射線照射の指示が待機される（ステップＳ３）。ここで、撮影実施者は、被写体Ｍを放射線源１１と放射線検出部１３の間に配置してポジショニングを行う。また、本実施形態においては、呼吸状態下で撮影を行うため、被検者（被写体Ｍ）に楽にするように指示し、安静呼吸を促す。或いは、「吸って、吐いて」等の深呼吸の誘導を行うこととしてもよい。撮影準備が整った時点で、操作部２３を操作して放射線照射指示を入力する。
【００５２】
操作部２３により放射線照射指示が入力されると（ステップＳ３；ＹＥＳ）、放射線照射制御装置１２及び読取制御装置１４に撮影開始指示が出力され、動態撮影が開始される（ステップＳ４）。即ち、放射線照射制御装置１２に設定されたパルス間隔で放射線源１１により放射線が照射され、放射線検出部１３によりフレーム画像が取得される。
【００５３】
予め定められたフレーム数の撮影が終了すると、制御部２１により放射線照射制御装置１２及び読取制御装置１４に撮影終了の指示が出力され、撮影動作が停止される。撮影されるフレーム数は、少なくとも１呼吸サイクルが撮影できる枚数である。
【００５４】
撮影により取得されたフレーム画像は順次撮影用コンソール２に入力され、撮影順を示す番号（フレーム番号）と対応付けて記憶部２２に記憶されるとともに（ステップＳ５）、表示部２４に表示される（ステップＳ６）。撮影実施者は、表示された動態画像によりポジショニング等を確認し、撮影により診断に適した画像が取得された（撮影ＯＫ）か、再撮影が必要（撮影ＮＧ）か、を判断する。そして、操作部２３を操作して、判断結果を入力する。
【００５５】
操作部２３の所定の操作により撮影ＯＫを示す判断結果が入力されると（ステップＳ７；ＹＥＳ）、動態撮影で取得された一連のフレーム画像のそれぞれに、動態画像を識別するための識別ＩＤや、患者情報、検査情報、放射線照射条件、画像読取条件、撮影順を示す番号（フレーム番号）等の情報が付帯され（例えば、ＤＩＣＯＭ形式で画像データのヘッダ領域に書き込まれ）、通信部２５を介して診断用コンソール３に送信される（ステップＳ８）。そして、本処理は終了する。一方、操作部２３の所定の操作により撮影ＮＧを示す判断結果が入力されると（ステップＳ７；ＮＯ）、記憶部２２に記憶された一連のフレーム画像が削除され（ステップＳ９）、本処理は終了する。この場合、再撮影が必要となる。
【００５６】
なお、被検者に呼吸センサー等を装着し、撮影用コンソール２において、動態画像を撮影する際に動態撮影と同期した状態で呼吸センサーから被検者の呼吸運動の情報（例えば、呼吸運動を示す波形等）を取得し、制御部２１により、取得した情報に基づいて、動態撮影時の呼吸運動を示す情報（例えば、呼吸回数、各呼吸における呼気区間と吸気区間の時間、各フレーム画像が撮影されたタイミングの呼吸状態（例えば、呼気区間、安静（最大）呼気位、吸気区間、安静（最大）吸気位等）等））を判別して判別結果を各フレーム画像に付帯させることが好ましい。
【００５７】
（診断用コンソール３の動作）
次に、診断用コンソール３における動作について説明する。
診断用コンソール３においては、通信部３５を介して撮影用コンソール２から動態画像の一連のフレーム画像が受信されると、受信された動態画像の一連のフレーム画像が記憶部３２に記憶される。記憶部３２に記憶されている動態画像の中から操作部３３により動態画像が選択され、構造物の位置計測の実行が指示されると、制御部３１と記憶部３２に記憶されているプログラムとの協働により、選択された動態画像に対して図３に示す位置計測処理が実行される。位置計測処理の各ステップにおけるバリエーションの処理や表示される各画面が備える機能についても、制御部３１と記憶部３２に記憶されているプログラムとの協働により実現される。
【００５８】
ここで、本実施形態における位置計測処理は、あるフレーム画像（ここでは、１番目のフレーム画像）における位置計測対象の構造物の位置（計測対象点という）をユーザーが指定し、指定された位置の画像をテンプレート画像としてテンプレートマッチングを行って構造物の位置を計測する処理である。本実施形態では、横隔膜の位置計測を行う場合を例にとり説明する。
【００５９】
以下、図３を参照して位置計測処理の流れについて説明する。
まず、選択された動態画像が記憶部３２から取得される（ステップＳ１１）。
【００６０】
次いで、表示部３４に初期設定入力画面３４１が表示され、ユーザーによる操作部３３の操作に応じて計測対象点の位置を含む初期設定が入力される（ステップＳ１２：初期設定入力部）。
【００６１】
図４は、初期設定入力画面３４１の一例を示す図である。初期設定入力画面３４１は、ユーザーが位置計測のためのテンプレート画像の作成と横隔膜位置候補の抽出範囲の設定に必要な情報を入力するための画面である。図４に示すように、初期設定入力画面３４１には、画像表示領域３４１ａと、横隔膜位置候補の抽出範囲の設定に必要なパラメーターを入力するための初期設定入力領域３４１ｂと、階調補正ボタン３４１ｃと、決定ボタン３４１ｄと、が設けられている。
【００６２】
画像表示領域３４１ａには、動態画像の１番目フレーム画像が表示されている。ユーザーは、画像表示領域３４１ａに表示されたフレーム画像上から操作部３３のタッチペンやマウスなどのポインティングデバイスにより計測対象点を指定することができる。例えば、画像表示領域３４１ａに表示されているフレーム画像における横隔膜の境界が操作部３３のポインティングデバイスにより指定（例えば、クリック）されると、指定された点の位置が計測対象点として、制御部３１のＲＡＭに保持される。
【００６３】
画像表示領域３４１ａに表示されている画像の階調は、ユーザー操作に応じて調整することができる。例えば、操作部３３により階調補正ボタン３４１ｃが押下されると、押下された回数に応じて画像表示領域３４１ａに表示されている画像の階調が調整される。これにより、ユーザーが最も横隔膜が見えやすい状態で計測対象点を指定することができる。あるいは、制御部３１により画像表示領域３４１ａに表示するフレーム画像に対して自動で階調補正処理を行い、階調補正処理されたフレーム画像を画像表示領域３４１ａに表示して計測対象点の指定を受け付けてもよい。 
【００６４】
なお、ステップＳ１２において、横隔膜の計測対象点の指定入力はユーザーが手動で行うこととして説明したが、表示されているフレーム画像から制御部３１が画像処理により自動で横隔膜の位置を認識して指定してもよい。
【００６５】
また、計測対象点を指定するためのフレーム画像は１番目のフレーム画像に限らず、ユーザーが所望するフレーム画像を用いることとしてもよい。例えば、ユーザーが所望するフレーム画像を選択する方法として、１枚ずつフレーム画像を見てユーザーが手動で選択するものが考えられる。あるいは、ユーザーが所望のフレーム番号を入力するためのユーザインターフェイスを備えてもよい。これにより、ユーザーの裁量で自由に計測対象点を決定することができる。また、例えば、初期設定入力画面３４１に、画像表示領域３４１ａに表示するフレーム画像を選択するためのプルダウンメニュー（例えば、安静呼気位、安静吸気位、安静呼気位→安静吸気位の中間、安静吸気位→安静呼気位の中間等）を設けてユーザーが所望のフレーム画像を選択できるようにしてもよい。また、横隔膜の動きや形状は安静呼気位がもっとも自然状態に近いため、安静呼気位のフレーム画像を画像表示領域３４１ａに自動で表示してユーザーが計測対象点を指定できるようにしてもよい。
【００６６】
また、初期設定入力画面３４１には、ユーザーによる計測対象点の指定を補助するための補助機能を設けることが好ましい。
例えば、画像表示領域３４１ａに表示されているフレーム画像上に、横隔膜の形態の基準になるような画像を半透明などにして重ねて表示する。例えば、安静呼気位のフレーム画像の横隔膜の領域を認識し、認識した横隔膜の形態の画像を画像表示領域３４１ａに表示されているフレーム画像に重ねることで、横隔膜の形態の特徴的な部分を表す適切な計測対象点をユーザーが把握した状態で指定することが可能となる。
【００６７】
別の方法としては、予め手動または自動で計測対象点のＸ座標を決定しておき、そのＸ座標をライン等で画像表示領域３４１ａに表示されている１番目のフレーム画像上に表示しておく。これにより、ユーザーが指定すべき計測対象点を把握しやすくなる。また、計測対象点のＸ座標が予め決定されている場合は、ユーザーの入力はＹ座標のみを受け付けるようにし、どのＸ座標が入力されても無視するようにしてもよい。これにより、ユーザーはＹ座標の調整のみに集中することができる。 ここで、フレーム画像の水平方向（左右方向）をＸ座標、垂直方向（上下方向）をＹ座標とする。
【００６８】
また、画像表示領域３４１ａに表示されているフレーム画像の横隔膜と肺野の境界が不明瞭な場合、ユーザーは小さな輝度勾配を目視で探して計測対象点を指定する必要がある。その補助機能として、ユーザー操作に応じて適宜画像を所望の大きさに拡大する機能を備えることとしてもよい。例えば、画像表示領域３４１ａに表示されているフレーム画像において操作部３３のマウス等の操作により指定された位置の周辺を拡大表示することとしてもよい。或いは、同じ目的で、ユーザー操作に応じて、フレーム画像にエッジ強調処理を施したエッジ強調画像を画像表示領域３４１ａに表示することとしてもよい。エッジ強調処理は、例えば、一次微分フィルター、ソーベルフィルター、プレヴィットフィルター等による処理が挙げられる。これにより、ユーザーが指定すべき計測対象点を把握しやすくなる。 
【００６９】
また、操作部３３により計測対象点の位置を指定する際、カーソル自体が目視の妨げとなって、指定しようとする計測対象点と実際の画像上での横隔膜の位置が僅かにずれることがある。それを防ぐために、カーソルを半透明で表示する。或いは十字型などのカーソルの位置を指す点を見えやすい形で表示する。画像上の輝度情報をもとにカーソルの色を変えてもよい。例えば、輝度が低い領域では明るい色のカーソルを表示する。これにより、横隔膜位置とカーソルの目視を補助することができる。また、ここまではユーザーが操作部３３のマウス等によりクリックした瞬間に計測対象点の位置を決定するＧＵＩを前提に説明してきたが、クリックした際には、その位置を計測対象点の位置を仮置きで表示するようにし、ドラッグ中は、計測対象点の仮置き位置をマウス位置に応じて移動して表示させるようにし、マウスを離した際に、そのマウスを離した位置を計測対象点として決定するようにしてもよい。 
【００７０】
また、ユーザーが手動で計測対象点を指定すると、指定ミスを起こす場合がある。その場合の修正手段として、自動で計測対象点の位置を調整する機能又は手動で計測対象点を修正・削除する機能を備えることが好ましい。
自動で計測対象点の位置を調整する機能としては、例えば、計測対象点付近で輝度勾配の大きい箇所に計測対象点を移動させる。輝度勾配の大きい箇所とは、エッジ強調画像上での画素値が高い箇所と読み替えてもよい。輝度勾配の大きい箇所が指定した計測対象点付近に見つからない場合は、入力ミスとして捉え、エラーを返して、再度計測対象点を入力するように促してもう１度計測対象点の指定を受け付ける。エラーの返し方として、例えば、エラーとなった計測対象点の色を変えること、エラー内容を文字で表現すること、エラーとなった計測対象点を別の形・記号に変えること、点滅させること、あるいは音を鳴らすことなどが挙げられる。横隔膜が動き得るエリアを予め決めて、そのエリア外に計測対象点を指定した場合に上記のエラーを返すことも考えられる。
【００７１】
手動で計測対象点を調整する機能として、大まかな位置を操作部３３のマウスで指定した後の次のステップとして、指定された計測対象点の位置を、例えばマウスホイール等の操作に応じて少しずつ移動させて微調整可能とすることが望ましい。もちろん、大きく調整するために、マウスのドラッグ操作等に応じて計測対象点の位置を大きく移動させてもかまわない。調整中の計測対象点は、点のサイズや色を変えたり、点滅させたりして表示することが好ましい。また、一度指定された計測対象点を、例えばDeleteキー等により削除して指定し直せる機能を備えることが好ましい。 
【００７２】
また、ユーザーの所望の計測対象点で計測を行うために、上記の修正なども含めて最終的に指定された点をユーザーが承諾する機能を備えることが好ましい。例えば、画像表示領域３４１ａの近傍に承諾ボタンを設け、操作部３３により承諾ボタンが押下されたり、指定された点をダブルクリックしたりする等、規定された承諾のためのアクションが行われた場合に、制御部３１は、指定された点を計測対象点として決定する。
【００７３】
計測対象点は、複数の構造物について指定してもかまわない。例えば、横隔膜と肺尖の点を指定することで、肺のサイズの時間変化を計測することができる。 
【００７４】
一方、１つの構造物（ここでは、横隔膜）について複数の計測対象点を指定してもかまわない。もしくは計測対象点を線で指定してもよい。これにより、後述するルート候補の限定において、複数点の空間的な関係に基づいてルート候補を限定することができ、位置測定の安定性が増す。また、計測対象点を複数又は線で指定してその複数点や線の位置を追跡的に測定することで、横隔膜についての形状変化を把握できるため、ＣＯＰＤなどの病態に特有の横隔膜の特徴的な変形の動きなどを捉えることができ、ＣＯＰＤの兆候のある患者を早期に検出できる。
計測対象点を複数点や線で指定する場合においても、上記の１点で指定する場合の機能は含んでいるものとする。 
【００７５】
複数の計測対象点を指定する場合、指定する点の個数は特に限定されない。また、複数のフレーム画像に計測対象点を打ってもよい。これにより、横隔膜の形状の移り変わりを捉えることが可能となる。例えば、肋横角を二次元的に追跡し、横隔膜の指定された各点をＹ方向に追跡し、それを複合して、横隔膜の境界全体の位置の追跡結果を作成してもよい。これにより、横隔膜の境界について、二次元的な位置移動を抽出することができる。
【００７６】
計測対象点を線で指定する（横隔膜の位置を線で指定する）ことができるようにするために、初期設定入力画面３４１には以下の機能を備えることが好ましい。例えば、基本機能としては、ユーザーが操作部３３のポインティングデバイスにより、画像表示領域３４１ａの所望の位置に線を引く機能を備える。ただし、この機能のみでは、線を正確に引くのは難しく、ユーザーの負担が大きい。そこで、ユーザーが容易に横隔膜の線を指定するために、操作部３３により横隔膜上に２点以上が指定された場合、指定された点の間の横隔膜位置を自動で埋める（補間する）機能を備える。指定された点の間を自動で埋める方法として、例えば、動的計画法やテンプレートマッチング、あるいは輝度勾配が強い箇所をなぞる等の方法が考えられる。また、指定された点の間に自動で線分を引いて、その線分付近の輝度勾配の大きい箇所に各点を調整して正しい横隔膜の線を形成してもよい。
他にも、横隔膜が含まれる領域をユーザーが操作部３３により図形で囲んで、その中で輝度勾配の大きい箇所を横隔膜として検知する機能を備えてもよい。ここでいう図形とはたとえば矩形、円形、もしくはペンツールなどで囲んだ任意の図形を指す。これらの機能により、ユーザーが容易に横隔膜の線を指定することが可能となる。
【００７７】
横隔膜の線を間違えて指定した場合、その修正をする必要がある。例えば、複数の点をユーザーが指定した場合については、ユーザーが操作部３３により、既に指定された点と点の間に新たに点を指定すると、制御部３１は、点間の距離を縮めて線分を再検知する。これにより横隔膜の線の指定の正確さが向上する。 或いは、誤検知した点を操作部３３により矩形などで囲んで削除できるようにし、操作部３３によりその区間に新たな点が指定されると、新たに指定された点の間を上述の方法で自動で埋めることとしてもよい。これにより、点群として誤検知した場合にまとめて削除して修正を行うことが可能となる。 
【００７８】
また、計測対象点を同じＸ座標に複数のフレーム画像に対して指定してもよい。これにより、指定点の位置情報を、ルート候補を限定する際や計測結果を修正する際の制限条件として用いることができる。例えば、ルート候補が指定点付近を通過するか否かをルート候補の限定条件として用いることができる。
【００７９】
初期設定入力領域３４１ｂには、図４に示すように、横隔膜位置候補の抽出範囲の設定に必要な情報、例えば、横隔膜の速度（移動速度）や加速度の閾値、テンプレートマッチングにおける探索方向等のパラメーターを入力するための入力欄が設けられている。ここで、横隔膜は基本的には上下方向に移動するが、想定される移動速度や加速度には限りがある。そこで、本実施形態では、横隔膜位置候補を抽出する際の抽出範囲を横隔膜の想定される移動速度及び加速度により制限することができ、初期設定入力領域３４１ｂでは、この横隔膜の速度及び加速度の閾値を入力することができる。また、横隔膜の移動は基本的には上下方向であるため、デフォルトの探索方向は上下方向であるが、初期設定入力領域３４１ｂでは、ユーザーが所望する探索方向を指定することができる。なお、本実施形態では、探索方向が上下方向である場合を例にとり説明する。
【００８０】
初期設定入力画面３４１において、操作部３３により決定ボタン３４１ｄが押下されると、初期設定入力画面３４１から入力された内容がＲＡＭに設定され、処理はステップＳ１３に移行する。
なお、初期設定入力画面３４１の表示をせずに、位置計測処理に必要なパラメーターのデフォルト値を予め設定しても構わない。これにより、必要な操作が減り、ユーザーの負担を軽減することができる。
【００８１】
ステップＳ１３においては、変数ｎに１が設定され（ステップＳ１３）、ｎ番目のフレーム画像に前処理が実行される（ステップＳ１４）。
ステップＳ１４においては、まず、ｎ番目のフレーム画像に対してノイズ除去処理及び階調補正処理が行われる。ノイズ除去処理としては、例えば、メディアンフィルター、移動平均フィルター、ガウシアン平滑フィルター等によるフィルター処理が挙げられる。次いで、エッジ強調処理が行われる。エッジ強調処理としては、例えば、一次微分フィルター、ソーベルフィルター、プレヴィットフィルター等によるフィルター処理が挙げられる。 これらの前処理をフレーム画像に施すことにより、そのフレーム画像に含まれるノイズ情報を除去して横隔膜を際立たせることができ、横隔膜位置の計測の精度を向上させることができる。
【００８２】
次いで、ｎ＝１であるか否かが判断される（ステップＳ１５）。
ｎ＝１であると判断された場合（ステップＳ１５；ＹＥＳ）、１番目のフレーム画像の計測対象点の周囲の領域からテンプレート画像が作成される（ステップＳ１６）。例えば、１番目のフレーム画像における計測対象点を中心とするｗ画素×ｈ画素（ｗ、ｈは正の整数）の領域がテンプレート画像として取得される。次いで、計測対象点のＹ座標がフレーム１の横隔膜位置としてルート保持部３２１に記憶され（ステップＳ１７）、ステップＳ２２の処理に移行する。
【００８３】
一方、ｎ＝１ではないと判断された場合（ステップＳ１５；ＮＯ）、ｎ番目のフレーム画像に対し、ステップＳ１６で作成されたテンプレート画像を用いたテンプレートマッチングが実行され、計測対象点と同じＸ座標上の予め定められたＹ座標の範囲内の各位置を中心とするｗ画素×ｈ画素の領域とテンプレート画像との相関値（例えば、相互相関係数）がそれぞれ算出される（ステップＳ１８：評価値算出部）。ステップＳ１８で算出される相関値は、横隔膜と各位置との類似度を表す評価値である。
【００８４】
なお、評価値の算出方法は相互相関（NCC; Normalized Cross-Correlation）が望ましいが、相互相関以外にも、SSD(Sum of Squared Difference)、SAD(Sum of Absolute Difference) 、ZNCC(Zero-means Normalized Cross-Correlation) を用いてもよい。
【００８５】
次いで、初期設定入力画面３４１で設定された速度及び加速度で制限した抽出範囲内で相関値が極大となる位置が横隔膜位置候補として抽出される（ステップＳ１９：位置候補抽出部）。
ここで、相関値が最大となる位置を横隔膜位置として決定するのではなく、相関値が極大となる位置を横隔膜位置候補として抽出しておくことで、絶対的に相関値が高くなくても空間的に見て相対的に相関値が高い位置を候補として残すことができ、大きな別構造物などと重なったケースなど絶対的な相関値が低くなってしまう条件下においても、横隔膜の位置を漏らさず位置候補として抽出することができる。なお、ここでは抽出範囲は一次元であるが、二次元での探索処理の場合、極大点は二次元にて抽出すればよい。また、ステップＳ１９においては、少なくとも一つのフレーム画像から複数の横隔膜位置候補が抽出されることを想定している。 
【００８６】
また、抽出する横隔膜位置候補を予め設定した横隔膜の速度及び加速度の範囲内の相関値が極大の位置に限定して抽出することで、横隔膜としてありえない位置候補を大幅に削減することができ、計算量を削減できる。ここで、横隔膜の速度は、フレーム画像間の横隔膜の移動量であり、例えば、１フレーム前の画像と当該フレーム画像との間の横隔膜位置のＹ座標の差分により算出できる。また、横隔膜の加速度は、フレーム画像間の横隔膜の移動量の変化であり、時間領域が重複しない隣接する２つの速度（前述）の差分により算出できる。例えば、２フレーム前の画像と１フレーム前の画像との間の横隔膜位置のＹ座標の差分と、１フレーム前の画像と当該フレーム画像との間の横隔膜位置のＹ座標の差分との差分により算出することができる。
【００８７】
図５に、ステップＳ１９における横隔膜位置候補の抽出を模式的に示す図である。図５に示すグラフの縦軸は、右のエッジ強調画像のＹ座標と対応している。時間Ｔ−２、時間Ｔ−１のフレーム画像において○で示した極大点が横隔膜位置候補として抽出されていた場合、時間Ｔのフレーム画像については、相関値が極大となる点のうち、設定した速度及び加速度の範囲内の点（○で囲んだ点）が横隔膜位置候補として抽出される。
【００８８】
なお、ステップＳ１８におけるテンプレートマッチングの探索範囲を初期設定入力画面３４１で設定された速度及び加速度に基づいて予め絞り込んでおくと、無駄な範囲の相関値算出を行わずに済み、計算量を削減することができるので好ましい。
【００８９】
ここで、１番目のフレーム画像の横隔膜位置は、指定された計測対象点の位置であるが、２番目以降のフレーム画像は、複数の横隔膜位置候補が抽出される可能性がある。例えば、２番目のフレーム画像からＰ個の横隔膜位置候補が抽出された場合、横隔膜が１番目のフレーム画像から２番目のフレーム画像までに移動したルートとしては、１フレーム目の計測対象点から２フレーム目のＰ個の横隔膜位置候補のそれぞれまでのＰ通りの候補が考えられる。３番目のフレーム画像から横隔膜位置候補を抽出する際には、テンプレートマッチングで得られた極大点のそれぞれをＰ通りのルート候補のそれぞれに追加した場合に速度及び加速度が予め定められた範囲内となるか否かを判断し、速度及び加速度が予め定められた範囲内となる極大点をそのルート候補の横隔膜位置候補とする。それ以降のフレーム画像についても同様に、ルート候補毎に横隔膜位置候補の抽出を行う。なお、ルート候補間で速度・加速度により制限する範囲が重複する場合は、まとめて１回だけ位置候補抽出を行ってもよい。これにより、計算の重複を避けることができる。
【００９０】
なお、ここでは、横隔膜位置候補の抽出範囲を横隔膜の速度及び加速度により制限することとして説明しているが、速度又は加速度の何れか一方により制限することとしてもよい。また、ｎ番目のフレーム画像に付帯されている呼吸状態を示す情報に基づいて、そのフレーム画像撮影時の呼吸状態（例えば、呼気区間、安静（最大）呼気位、吸気区間、安静（最大）吸気位）を判定し、横隔膜の移動方向がフレーム画像の撮影時の呼吸状態に合致しているか否かに基づいて、抽出する横隔膜位置候補を制限することとしてもよい。これらの制限は組み合わせて用いてもよい。
【００９１】
次いで、ｎ番目のフレーム画像の横隔膜位置候補を追加したルート候補が作成され、ルート保持部３２１に保存される（ステップＳ２０）。
図６に、ルート保持部３２１の一例を示す。ルート保持部３２１は、１〜Ｎ番目のフレーム画像のそれぞれから抽出された横隔膜位置候補を一つずつ時系列に繋いで横隔膜が移動したルート候補の情報として記憶するものである。図６に示すように、ルート保持部３２１は、ルート候補を識別する情報（例えば、ルート１、ルート２、ルート３、・・・）と、そのルート候補における各フレーム画像の横隔膜位置候補のＹ座標とを対応付けて記憶する。ここで、図６におけるＹ２１とは、２番目のフレーム画像の１番目の横隔膜位置候補であることを示している。（未計算）とは、まだ横隔膜位置候補が抽出されていないことを示す。なお、未計算であるか否かを示す情報は、フレーム画像毎に別途保持してもよい。また、ルート保持部３２１は、図６に示したものに限定されず、例えば、ルート３とルート４はフレーム１〜３の範囲で同一であるため、これを圧縮して（例えば、ルート３の欄にＹ座標の値を格納し、ルート４にはルート３と同じであることを示すフラグ等をたてて）保持させてもよい。また、横隔膜位置候補の”Ｙ座標”だけでなく”相関値”を併せて記憶することとしてもよい。
【００９２】
ステップＳ２０においては、抽出された横隔膜位置候補が１つだったルート候補については、ルート保持部３２１のｎ番目のフレーム画像の欄にその横隔膜位置候補のＹ座標を追加して保存する。抽出された横隔膜位置候補が２以上だったルート候補については、ルート保持部３２１に新たなルート候補を格納する行を追加して記憶する。ここで、追加するルート候補の内容は、ｎ番目のフレーム画像の横隔膜位置候補のみが異なり他のフレーム画像は同じ横隔膜位置候補の情報をコピーして記憶する。
【００９３】
次いで、ルート保持部３２１に記憶されたルート候補をより少ない複数のルート候補に限定する（ステップＳ２１：ルート限定部）。
横隔膜位置のテンプレート画像との相関値は、他の臓器との重なり等により局所的な時間で見ると最大にならない場合がある。しかし、時間方向に大域的に見ると横隔膜位置のテンプレート画像との相関値は最大になることが多い。そこで、ルート保持部３２１に記憶されているルート候補を大域的な時間範囲の中で相関値が高いルート候補に限定する。ルート候補を限定するタイミングとしては、全フレーム画像までのルート候補が出揃った時点が考えられるが、本実施形態のように逐次処理中にルート候補の限定を行うことも可能である。逐次処理中にルート候補の数を限定していくことで、計算量が減り、処理速度が向上する。ルート候補を保持するために使用される記憶容量も大幅に低減することができる。
【００９４】
ステップＳ２１における具体的な方法としては、例えば、或るフレーム画像までの横隔膜位置候補を抽出した時点でのルート候補数が、予め設定された値よりも大きくなった時に、その時点でルート候補を限定する。ルート候補を限定する方法としては、例えば、図７に示すように、それまでに算出された横隔膜位置候補のテンプレート画像との相関値（以下、横隔膜位置候補の相関値と呼ぶ）の合計値または平均値が所定の閾値以上のルート候補に限定する。或いは、横隔膜位置候補の相関値の合計値または平均値が高いものを上位から任意の個数だけに限定するという方法でもよい。また、或るフレーム画像時点での横隔膜位置候補の相関値に予め設定された閾値以上のものが存在する時に、その横隔膜位置候補を含むルート候補に限定してもよい。また、横隔膜位置候補の相関値が最大のフレーム画像数が多いルート候補を上位から任意の個数に限定することとしてもよい。
その他、それまでに算出された横隔膜位置候補のうち、上述の相関値が所定の閾値以上又は所定の閾値未満の横隔膜位置候補をもつフレーム画像の数又はその割合、複数のルート候補間で同一のフレーム画像から算出された横隔膜位置候補の相関値を比較したときの高さ順位の合計値や平均値、複数のルート候補間で同一のフレーム画像から算出された横隔膜位置候補の相関値を比較したときの高さ順位が所定の閾値以上又は所定の閾値未満のフレーム画像の数又はその割合等に基づいて、ルート候補の限定を行うこととしてもよい。
また、複数のフレーム画像の同じＸ座標に対して計測対象点が指定されている場合、指定された各計測対象点付近を通過していることをルート候補の限定条件として用いることとしてもよい。
【００９５】
また、ルート候補を限定するための評価指標として、相関値ではなく、横隔膜位置候補の移動量及び／又は移動方向を用いることもできる。横隔膜は、呼吸運動に伴い、ある決まった範囲内で上下の周期運動を行うので、移動量は限定されており、また移動方向も断続的に変化しないからである。例えば、横隔膜位置候補の移動量が小さいルート候補を上位から任意の個数に限定することとしてもよい。或いは、横隔膜の動きが呼吸運動に合致しているか否かに基づいてルート候補の限定を行ってもよい。例えば、フレーム画像間の横隔膜の速度の正負が予め設定した個数以上連続して同じルート候補に限定してもよい。
【００９６】
また、上記のルート候補の限定の方法はすべて、横隔膜位置候補の抽出が終わった全フレーム画像の情報を用いるものだが、他にも分岐したフレーム画像以降のフレーム画像の情報を用いる方法も考えられる。例えば、図７では、分岐した６フレーム目以降のフレーム画像の情報を用いる。これによって、ルート候補限定の判断材料となる時間領域のみを見ることができるので、情報の分散が大きくなり、ルート候補の限定精度が向上する。また、計算量を低減させることができ、処理速度も向上する。具体的な方法としては、分岐した時点のフレーム画像から最終フレーム画像までの情報を用いる。或いは、分岐した時点のフレーム画像から、横隔膜位置候補が１つに収束されるフレーム画像までの情報のみを用いてもよい。更に、分岐したフレーム画像の前後のフレーム画像の情報を加えてもよい。これにより、分岐前後での時間的な一貫性をもとに、ルート候補の妥当性を判断できる。例えば、１秒程度の範囲のフレーム画像の情報を加えることで、各ルート候補が横隔膜らしい緩慢な動きを示しているか否か、移動方向の変化が呼吸運動に即したものであるか等を判断することが可能となる。 
【００９７】
また、ステップＳ１２において、複数の計測対象点（複数点）が指定された場合は、更に、各フレーム画像における複数点の空間的な関係に基づいて、ルート候補を限定することができる。例えば、複数点の位置候補から想定される横隔膜の境界線が断線したり、滑らかでなくなったり、想定される形状ではなくなる場合に、その原因となる横隔膜位置候補を含むルート候補を削除する。複数点の位置候補から想定される横隔膜の境界線が想定される形状であるか否かは、例えば、２次関数や３次関数で曲線近似したときのズレの合計値で判断することができる。或いは、複数のフレーム画像において複数点の空間的な関係に基づいて、ルート候補を限定してもよい。横隔膜においては、どの箇所でも呼気と吸気で動く方向が一方向となるため、複数のフレーム画像における複数点の動く方向が一致していない場合、他と違う動きをしている横隔膜位置候補を含むルート候補を削除してもよい。複数点のＹ座標の動き幅が肋横角（体側）に近いほど大きく脊椎側に近いほど小さくなっているか否かに基づいてルート候補を限定することもできる。 
【００９８】
なお、ルート候補の限定に用いる評価指標は単独だけでなく、複数の評価指標を組み合わせて用いてもよい。
例えば、相関値の合計と横隔膜位置候補の移動量を評価指標として用いる場合、それぞれを独立に用いてルート候補を限定し、両評価指標において所定の条件を満たしたものをルート候補として残すことができる。これにより、より一層ルート候補の数を少なくすることができる。或いは、２つの評価指標に係数により重みをつけて合算することで１つの評価指標として用いることもできる。これにより、複数の評価指標であっても１次元的にルート候補を限定することができる。
【００９９】
ステップＳ２１におけるルート候補の限定は、複数のルート候補がある場合にのみ実行される。また、ルート保持部３２１に記憶されているルート候補が所定の条件に合致した場合（例えば、ルート候補数が予め設定された値よりも大きくなった場合）に実行される。
【０１００】
次いで、全フレーム画像についてステップＳ１４〜Ｓ２１の処理が終了したか否かが判断される（ステップＳ２２）。
全フレーム画像についてステップＳ１４〜Ｓ２１の処理が終了していないと判断した場合（ステップＳ２２；ＮＯ）、変数ｎが１インクリメントされ（ステップＳ２３）、ステップＳ１４に戻る。
【０１０１】
全フレーム画像についてステップＳ１４〜Ｓ２１の処理が終了したと判断された場合（ステップＳ２２；ＹＥＳ）、一つのルート候補が横隔膜の移動を示す横隔膜の移動ルートとして決定され、決定されたルートに含まれる各フレーム画像の横隔膜位置候補が各フレーム画像の横隔膜位置として決定される（ステップＳ２４：ルート決定部）。
ステップＳ２４におけるルートの決定方法としては、ステップＳ２１でルート候補を限定したときに説明した方法を適用することができ、ルート候補の限定数を任意の個数ではなく１つに決定する。また、動態画像の付帯情報として、動態撮影時の呼吸運動を示す情報が付帯されている場合、ルート候補のそれぞれにおける横隔膜の移動が動態撮影時の呼吸運動に対応しているか否かに基づいて、ルート候補を限定することが好ましい。例えば、動態画像の付帯情報として呼吸回数が設定されていた場合に、撮影全体の横隔膜候補位置により想定される呼吸回数が付帯情報の呼吸回数と一致するものにルート候補を限定してからルート決定を行うことが望ましい。
【０１０２】
次いで、横隔膜位置の計測結果が表示部３４に表示される（ステップＳ２５：結果出力部）。
図８に、ステップＳ２５において表示部３４に表示される計測結果画面３４２の一例を示す。図８に示すように、計測結果画面３４２には、各フレーム画像の横隔膜位置に点が付与された動態画像を表示するための画像表示領域３４２ａと、横隔膜位置のＹ座標を時間軸にプロットした横隔膜波形を示すグラフ３４２ｂと、計測結果の妥当性の判断結果３４２ｃと、承諾ボタン３４２ｄと、修正ボタン３４２ｅと、やり直しボタン３４２ｆと、が設けられている。なお、グラフ３４２ｂの縦軸は、放射線検出部１３上の実際の画素の位置で表現しても構わない。
【０１０３】
画像表示領域３４２ａに表示される動態画像上の横隔膜位置に表示する点のサイズや線の太さは問わない。点の形状は、例えば、円や正方形などが想定される。また、点のサイズや点の形状は、操作部３３によりユーザーが指定してもデフォルトで設定されていても構わない。また、ユーザーが目視レベルで計測対象点の位置変化の特性を分析することを想定して、例えば、横隔膜位置の速度や加速度によって点の色を変化させる機能を備えることが好ましい。
【０１０４】
グラフ３４２ｂにおいては、フレーム画像毎の横隔膜位置を点で表現してもよいし、点を結んだ線で表現してもよい。フレーム画像毎の横隔膜位置を点で表現する場合は、点のサイズは問わない。点の形状についても円や正方形等が想定されるが、特に限定されない。点のサイズや点の形状は、操作部３３によりユーザーが指定してもデフォルトで設定されていても構わない。
一方、フレーム画像毎の横隔膜位置を線として表現する場合は、線の種類として望ましくは折れ線や近似曲線が考えられる。また点の場合と同様に、線の太さについては特に限定されない。
複数の計測対象点を同時に指定した場合、計測結果は１枚のグラフにまとめて表示される。この場合、計測対象点ごとにグラフの色分けをする機能を備えることが好ましい。横隔膜位置を点で表現する場合は、計測対象点によって色ではなく形状を変えることで区別してもよい。 
【０１０５】
また、ユーザーが参照したいフレーム画像を効率よく表示できるようにするために、操作部３３によりグラフ３４２ｂ上の１点が指定されると、画像表示領域３４２ａにおいて、対応するフレーム画像を瞬時に表示する機能を備えることが好ましい。
【０１０６】
ユーザーが横隔膜の位置変化を詳細に分析する場合、グラフではなく、横隔膜の位置座標を数値として表示することも必要となる。そこで、計測結果画面３４２には、図８の表示に加え、フレーム画像毎の横隔膜位置の座標を表として表示することが好ましい。この時、望ましくは、フレーム画像毎の横隔膜位置座標を表示するだけでなく、表のセルを速度や加速度によって異なる色で塗りつぶし、カラーマップとして表示する機能も備えることが好ましい。これにより、複数の計測対象点を指定して計測した際、計測対象点毎の横隔膜の移動特性をユーザーが判断しやすくなる。
また、操作部３３により表の値が指定されると、画像表示領域３４２ａに、対応するフレーム画像を瞬時に表示する機能を備えることが好ましい。これにより、ユーザーが参照したいフレーム画像を効率よく表示することが可能となる。
また、フレーム画像毎の横隔膜位置座標の数値データをＣＳＶファイルなどの編集可能な形式で通信部３５により外部のコンピュータ装置等に出力することが望ましい。これにより、ユーザーが横隔膜位置座標の数値データを用いて、横隔膜の位置変化を定量的に分析することができる。
【０１０７】
また、ユーザーが計測結果を容易に診断に役立てることができるようにするために、計測された横隔膜位置に基づいて、制御部３１により各種特徴量を算出して表示部３４の画面上に画面出力、ならびにファイル出力をする機能を備えることが好ましい。算出する特徴量としては、例えば以下の１〜６が挙げられる。 
【０１０８】
１．横隔膜移動量の最大値
横隔膜移動量の最大値は、計測された横隔膜位置のＹ座標の最大値と最小値の差の絶対値である。 
２．横隔膜速度の最大値
横隔膜速度の最大値は、単位時間あたり、または単位フレームあたりの横隔膜位置の移動量の最大値である。基本的には動態画像の全フレーム画像の計測結果から算出するが、呼気時の計測結果から呼気時の横隔膜速度の最大値を、吸気時の計測結果から吸気時の横隔膜速度の最大値をそれぞれ算出することとしてもよい。 
３．呼吸時間
呼吸時間とは、呼気にかかる時間、吸気にかかる時間、又は１回の呼吸（呼気＋吸気）にかかる時間のことである。呼気にかかる時間は、横隔膜のＹ座標が上方向に移動している時間を算出することにより求めることができる。吸気にかかる時間は、横隔膜のＹ座標が下方向に移動している時間を算出することにより求めることができる。１回の呼吸に係る時間は、呼気時間と吸気時間を加算することで求めることができる。呼気にかかる時間と吸気にかかる時間とを比較する場合は、例えば、呼気にかかる時間と吸気にかかる時間との比や差を求める。 
４．呼吸のばらつき
複数回の呼吸中の動態画像を撮影する場合、横隔膜の移動量や速度、呼吸時間は、呼吸間でばらつきがある。そのばらつき度合いを定量化する方法として、たとえば、複数回の呼吸において上述の横隔膜移動量の最大値、横隔膜速度の最大値、呼吸時間等の特徴量をそれぞれ算出し、それぞれの分散や標準偏差、特徴量の（最大値−最小値）等を求めることが考えられる。呼吸毎の横隔膜移動量でなく、呼吸毎の最大の吸気位、最大の呼気位の横隔膜位置座標でそれぞれ同様の計算をしてもよい。呼吸時間に関しても、呼気にかかる時間、吸気にかかる時間でそれぞれ同様の計算をしてもよい。 
５．横隔膜と任意の基準点との距離の最大値
横隔膜と任意の基準点との距離の最大値とは、横隔膜位置が或る基準点から最も離れたときの横隔膜位置と基準点の距離である。ある基準点とは、望ましくは肺尖部を想定する。他にも基準点として、鎖骨、胸郭と鎖骨の交差する箇所、肺門、気管分岐部などが考えられるが、これに限ったものではない。 
６．横隔膜移動量の最大値、横隔膜速度の最大値を正規化したもの
横隔膜移動量の最大値、横隔膜速度の最大値の特徴量を他人の胸部の動態画像と比較して用いる場合、被写体間で特徴量を正規化する必要がある。そのために、横隔膜移動量の最大値、横隔膜速度の最大値の特徴量を被写体固有の大きさを表す値で除算する。被写体固有の大きさを表す値とは、例えば、横隔膜と任意の基準点との距離の最大値の特徴量や身長、肺野の面積、胸郭の幅などが考えられるが、これに限ったものではない。また横隔膜速度の最大値の特徴量に関しては、その被写体の横隔膜の可動域で正規化する方法も考えられる。即ち、横隔膜速度の最大値を横隔膜移動量の最大値で除算する。
【０１０９】
上記の横隔膜位置計測結果から得られる特徴量を数値データとして出力してもかまわないが、グラフや動態画像とあわせて表示をさせてもよい。数値データを表示させる際の表示方法は以下のとおりである。
例えば、１つの特徴量のみを表示する場合、計測結果画面３４２に、動態画像やグラフ３４２ｂに加えて、算出した特徴量を記載した特徴量表示用のテーブルを表示することが好ましい。このとき、特徴量から病気を示唆できる場合は、その特徴量から分かる病気の重度によって区別できるような表示にする。特に、病気の重度が高いものほど目立つように表示する。例えば、正常、注意、異常の３段階で、正常な場合は青や緑のような寒色に近い色で、注意や異常な場合は黄色や赤のような暖色に近い色で文字や表のセルを表示することが考えられる。同様に、異常に近づくほど文字サイズを大きくしたり、あるいは文字サイズを太くしたりしてもよい。この表示形式は後述のグラフや動態画像に表示する場合も同様である。 
【０１１０】
また、特徴量は、グラフ３４２ｂに書き込んで表示してもよい。１つの特徴量のみをグラフに書き込む際の表示方法について記す。
グラフ３４２ｂに特徴量を書き込む場合は、グラフ３４２ｂにおいてその特徴量の算出に用いた値の位置付近に表示することが望ましい。例えば、横隔膜速度の最大値の特徴量をグラフ３４２ｂに書き込む場合は、横隔膜の速度が最大となった時刻、位置の付近に書き込む。ただし、表示位置はこの限りではない。 
画像表示領域３４２ａに表示された動態画像に書き込む場合は、例えば、動態画像の右上等の横隔膜の目視の妨げにならない場所に表示することが望ましいが、この限りではない。また、その特徴量が算出されたフレーム画像における特徴量表示は他と区別できるようにすることが望ましい。例えば、横隔膜速度の最大値の特徴量の場合、横隔膜の速度が最大となったフレーム画像で特徴量表示のサイズを大きくすることが考えられる。特徴量表示の色を変えることとしてもよい。
【０１１１】
また、複数の特徴量を表示する場合、個々の特徴量の表示方法は１つの特徴量表示の際と同様だが、複数個表示するとなると画面が見にくくなるという問題が発生する。そこで、表示する特徴量を限定することを考える。例えば、病気を示唆する値を示した特徴量のみを表示することが考えられる。また、ユーザーの見たい特徴量のみを表示するという方法も考えられる。その場合、ユーザーの所望の特徴量をユーザー自身が操作部３３により選択できるユーザインターフェイスを備える。あるいは、診断に用いられる頻度の高い特徴量をデフォルトとして表示し、ユーザーがそれに特徴量を増減できるようなユーザインターフェイスを備えてもよい。これらの方法は、特徴量を数値データとして表示する際も、グラフ３４２ｂや画像表示領域３４２ａに表示された動態画像に表示する際にも共通する。
【０１１２】
計測結果の妥当性の判断結果３４２ｃには、妥当性が高いことを示すＧｏｏｄ又は妥当性が低いことを示すＰｏｏｒが表示される。計測結果の妥当性が低いとは、例えば、膨大（予め定められた所定数以上）の数のルート候補が作成され、その中から一つのルートに決定された場合や、決定されたルート上に相関値の高い（所定の閾値以上の）横隔膜位置候補が少なかった（予め定められた所定数を下回った）場合、あるいは相関値の低い（所定の閾値を下回った）横隔膜位置候補が多い（予め定められた所定数以上である）ことなどが挙げられる。膨大な数のルート候補が作成されたということは、１フレーム内で横隔膜位置を判断できないフレーム画像が多いということなので、ステップＳ２４で誤ったルートを採用する可能性がある。また、相関値は計測対象点を含む横隔膜領域との類似度を示す値なので、それが低いフレーム画像数の多いルートは横隔膜を計測できたという確からしさが低いと言える。そのため、ユーザーに計測結果を注意深く見るように促すために警告として、計測結果の妥当性の低さを通知する必要がある。計測結果の妥当性が低い場合、その旨を示す「Ｐｏｏｒ」が表示されるが、例えば、ルート候補の数が膨大である等、妥当性の低い具体的な理由を表示してもよい。
【０１１３】
また、計測結果画面３４２には、最終的にユーザーが計測結果を承諾するのか、修正するのか、あるいは計測をやり直すのかを判断した判断結果を入力するためのユーザインターフェイスとして、承諾ボタン３４２ｄ、修正ボタン３４２ｅ、やり直しボタン３４２ｆが設けられている。修正ボタン３４２ｅを押下すると、計測結果を修正するための結果修正画面３４３に移行する。やり直しボタン３４２ｆを押下すると、初期設定入力画面３４１に移行して、計測をやり直すことができる。
【０１１４】
次いで、操作部３３により計測結果画面３４２の承諾ボタン３４２ｄ（又は、後述する承諾ボタン３４３ｄ、３４４ｄ）が押下されたか否かが判断される（ステップＳ２６）。
【０１１５】
操作部３３により承諾ボタン３４２ｄ（又は、後述する承諾ボタン３４３ｄ、３４４ｄ）が押下されていないと判断された場合（ステップＳ２６；ＮＯ）、操作部３３により修正ボタン３４２ｅ（又は、後述する修正ボタン３４３ｅ、３４４ｅ）が押下されたか否かが判断される（ステップＳ２７）。
修正ボタン３４２ｅ（又は、後述する修正ボタン３４３ｅ、３４４ｅ）が押下されたと判断された場合（ステップＳ２７；ＹＥＳ）、表示部３４に結果修正画面が表示され（ステップＳ２８）、操作部３３の操作に応じて計測結果が修正され（ステップＳ２９）、処理はステップＳ２６に戻る。
【０１１６】
計測結果の修正方法として、動態画像の全フレーム画像において、計測ミスをした点を正しい位置にユーザーが操作部３３により手動で修正（例えば、ドラッグする）というものが考えられる。ただし、計測ミスしたフレーム画像数が多い場合、計測ミスをした全てのフレーム画像をユーザーが手動で修正するのは困難である。そのため、結果修正画面においては、望ましくは、ユーザー操作に応じて動態画像の各フレーム画像（決定された横隔膜位置に点等のマークを付したもの）を切り替え可能に表示し、操作部３３によりユーザーが横隔膜位置を修正した（ドラッグした）フレーム画像の横隔膜位置情報をもとにその周囲のフレーム画像についても自動で修正する機能を備えることが望ましい。このときの自動修正方法としては、例えば、手動で修正されたフレーム画像の横隔膜位置と周囲のフレーム画像において計測された横隔膜位置とのフレーム間平均による補間処理を用いることができるが、この限りではない。 
【０１１７】
また、上述のように、複数のルート候補から一つのルートが決定されるが、その結果が残念ながらユーザーの所望の結果と異なる場合も存在する。その場合、ユーザーが適切なルートを容易に選択できるようにするために、一つのルートが決定される前の複数のルート候補を表示し、表示した複数のルート候補の中から最も適するものをユーザーが選択するためのユーザインターフェイス（例えば、結果修正画面３４３）を表示部３４に表示することとしてもよい。
【０１１８】
図９に、上述の結果修正画面３４３の一例を示す。図９に示すように、結果修正画面３４３には、画像表示領域３４３ａと、グラフ３４３ｂと、計測結果の妥当性の判断結果３４３ｃと、承諾ボタン３４３ｄと、修正ボタン３４３ｅと、やり直しボタン３４３ｆと、が表示されている。画像表示領域３４３ａには、図９に示すように、一つのルートが決定される前の複数のルート候補の横隔膜位置候補Ｐ（図９ではＰ１〜Ｐ３の３つ）が表示される。また、グラフ３４３ｂには、複数のルート候補の横隔膜位置候補のＹ座標の時間変化Ｒ（図９ではＲ１〜Ｒ３）が表示される。操作部３３により複数の横隔膜位置候補Ｐの中から何れかが選択された場合、又は、操作部３３により複数のルート候補Ｒの中から何れかが選択された場合、選択された点Ｐに対応するルート候補、又は選択されたルート候補Ｒが横隔膜位置のルートとして決定される。
【０１１９】
また、新たに計測対象点を入力可能なユーザインターフェイス（例えば、結果修正画面３４４）を表示することとしてもよい。
図１０に、結果修正画面３４４の一例を示す。図１０に示すように、結果修正画面３４４には、画像表示領域３４４ａと、グラフ３４４ｂと、計測結果の妥当性の判断結果３４４ｃと、承諾ボタン３４４ｄと、修正ボタン３４４ｅと、やり直しボタン３４４ｆと、が表示されている。
例えば、図１０に示すように、画像表示領域３４４ａには、始めに計測対象点（計測対象点Ｐ０とする）を指定したフレーム画像が表示される。操作部３３により、計測対象点Ｐ０付近に新たに計測対象点Ｐ１が指定されると、制御部３１により、新たに指定された計測対象点Ｐ１に基づいて横隔膜の移動ルートが計測される。そして、各フレーム画像について、始めに指定した計測対象点Ｐ０に基づいて計測した複数のルート候補と新たに指定した計測対象点Ｐ１に基づく横隔膜の移動ルートとの位置関係から始めに計測した複数のルート候補から妥当性の高いルートを選択する。或いは、始めに計測対象点Ｐ０を指定したフレーム画像と異なるフレーム画像を画像表示領域３４４ａに表示することにより、始めに指定された計測対象点Ｐ０と同じＸ座標の計測対象点の指定を異なるフレーム画像上で受け付け、すでに作成されている複数のルート候補のうち、新たに指定された計測対象点を通るルート候補に絞ることとしてもよい。これにより、精度よくルートを決定することが可能となる。
【０１２０】
ステップＳ２７において、修正ボタン３４２ｅ（又は、後述する修正ボタン３４３ｅ、３４４ｅ）が押下されず（ステップＳ２７；ＮＯ）、やり直しボタン３４２ｆ（又は、後述するやり直しボタン３４３ｆ、３４４ｆ）が押下されたと判断した場合（ステップＳ２８；ＹＥＳ）、処理はステップＳ１２に戻り、表示部３４に初期設定入力画面３４１が表示され、位置計測のやり直しが実行される。
【０１２１】
一方、操作部３３により承諾ボタン３４２ｄ（又は、３４３ｄ、３４４ｄ）が押下されたと判断された場合（ステップＳ２６；ＹＥＳ）、記憶部３２に記憶されている動態画像に位置計測結果が対応付けて記憶され（ステップＳ３１）、位置計測処理は終了する。
【０１２２】
以上説明したように、診断用コンソール３によれば、制御部３１は、胸部の動態画像の複数のフレーム画像のそれぞれにおける所定範囲内の各位置について、横隔膜との類似度を示す評価値を算出し、算出された評価値に基づいて、複数のフレーム画像のそれぞれから横隔膜位置候補を一つ以上抽出し、複数のフレーム画像のそれぞれから抽出された横隔膜位置候補を一つずつ時系列に繋いで横隔膜が移動したルート候補として複数のルート候補をルート保持部３２１に記憶させる。そして、ルート保持部３２１に記憶されている複数のルート候補の中から一つのルート候補を横隔膜の移動ルートとして決定し、決定したルートに含まれる横隔膜位置候補を複数のフレーム画像のそれぞれにおける横隔膜の位置として決定する。
従って、例えば、他の構造物との重なりが生じ、局所的な時間で見ると実際の横隔膜位置の評価値が最大とならないケースであっても、横隔膜位置候補として抽出し、横隔膜位置候補を時系列に繋いで大域的な時間範囲の中で横隔膜位置を決定することができるので、マーカを用いずに撮影した胸部の動態画像において、横隔膜の境界（輪郭）が不鮮明であっても、各フレーム画像における横隔膜位置を精度良く計測することが可能となる。
【０１２３】
例えば、制御部３１は、複数のフレーム画像のそれぞれにおいて、算出された評価値が極大となる位置を横隔膜位置候補として抽出する。従って、絶対的に評価値が高くなくても、空間的に見て相対的に評価値が高い位置を横隔膜位置候補として残すことができ、他の構造物と重なったケース等、ある空間範囲において絶対的な評価値が低くなってしまう条件下においても、漏らさず横隔膜位置候補として抽出することができる。
【０１２４】
また、例えば、制御部３１は、横隔膜位置候補の抽出範囲を横隔膜の移動速度及び／又は加速度に基づいて制限するので、横隔膜としてあり得ない横隔膜位置候補が抽出されることを低減することができる。また、計算量を削減することもできる。
【０１２５】
また、例えば、制御部３１は、横隔膜の移動方向が複数のフレーム画像のそれぞれを撮影したときの呼吸状態に合致するか否かに基づいて、各フレーム画像から抽出する横隔膜位置候補を制限するので、横隔膜としてあり得ない横隔膜位置候補が抽出されることを低減することができる。計算量を削減することもできる。
【０１２６】
また、例えば、制御部３１は、ルート保持部３２１に記憶されている各ルート候補に含まれる横隔膜位置候補のそれぞれについて算出された評価値に基づいて、横隔膜の移動ルートを決定することができる。例えば、ルート候補のそれぞれにおける、横隔膜位置候補の評価値の合計値又は平均値、横隔膜位置候補の評価値が所定の閾値以上又は所定の閾値未満のフレーム画像数又はフレーム画像数の割合、複数のルート候補間で同一のフレーム画像から算出された横隔膜位置候補の評価値を比較したときの高さ順位の合計値又は平均値、又は、複数のルート候補間で同一のフレーム画像から算出された横隔膜位置候補の評価値を比較したときの高さ順位が所定の閾値以上又は所定の閾値未満のフレーム画像数又はフレーム画像の割合に基づいて、横隔膜の移動ルートを決定することで、精度良く横隔膜の移動ルートを決定できる。
【０１２７】
また、例えば、制御部３１は、ルート保持部３２１に記憶されている各ルート候補における横隔膜の移動量及び／又は移動方向に基づいて、横隔膜の移動ルートを決定することで、横隔膜の想定される周期運動に基づいて、精度良く横隔膜の移動ルートを決定することができる。
【０１２８】
また、制御部３１は、横隔膜位置を計測するための初期設定をユーザーが入力するための初期設定入力画面を表示部３４に表示することで、横隔膜位置の計測に係る設定情報をユーザーが入力することが可能となる。
【０１２９】
また、制御部３１は、横隔膜位置の計測結果を出力することで、ユーザーに横隔膜位置の計測結果を提供することができる。
【０１３０】
また、制御部３１は、ルート保持部３２１に記憶されている複数のルート候補をより少ない複数のルート候補に限定することで、ルート保持部３２１で必要とする記憶容量を低減することができる。
【０１３１】
例えば、制御部３１は、ルート保持部３２１に記憶されている各ルート候補に含まれる横隔膜位置候補のそれぞれについて算出された評価値に基づいて、ルート候補を限定する。例えば、他の構造物との重なりのないフレーム画像においては、横隔膜位置と合致した横隔膜位置候補の評価値が明らかに高くなるため、制御部３１は、評価値が所定の閾値以上の横隔膜位置候補が存在した場合に、その位置候補を含むルート候補に限定する。これにより、あり得ないルート候補（明らかに横隔膜位置に合致している横隔膜位置候補を通らないルート候補）を除外することができる。
【０１３２】
また、例えば、制御部３１は、ルート保持部３２１に記憶されている各ルート候補における横隔膜の移動量及び／又は移動方向に基づいて、ルート候補を限定することで、横隔膜としてあり得ない動きとなるルート候補を除外することができる。
【０１３３】
また、例えば、制御部３１は、ルート候補のそれぞれにおける横隔膜の動きが動態撮影時の呼吸運動に対応しているか否かに基づいて、ルート候補を限定することで、動態撮影時の呼吸運動に対応していない横隔膜の動きを示すルート候補を除外することができる。
【０１３４】
なお、本実施形態における記述は、本発明に係る好適な位置計測システムの一例であり、これに限定されるものではない。
【０１３５】
例えば、上記実施形態においては、被写体部位が胸部である場合を例にとり説明したが、これに限定されず、他の部位を撮影した動態画像としてもよい。
【０１３６】
また、本実施形態においては、横隔膜位置を計測する場合を例にとり説明したが、本発明は、例えば、心壁、肋骨、血管等の他の構造物に対しても適用可能である。ただし、構造物によっては、位置の探索方向を調整する必要がある。例えば、計測対象の構造物が心壁である場合、探索方向を左右方向（Ｘ方向）とすることが好ましい。肋骨、血管の場合は、横隔膜と同様に、Ｙ方向とすることが好ましい。また、構造物によって、Ｙ方向やＸ方向ではなく、斜め一方向を探索方向としてもよい。また、２次元の位置変化を探索して計測することとしてもよい。
【０１３７】
また、上記実施形態においては、１番目のフレーム画像から順次時間後方に向かって横隔膜位置候補を抽出していく場合を例にとり説明したが、これに限定されず、例えば、或る基準となるフレーム画像から時間前方と時間後方のそれぞれに向かって順次横隔膜位置候補を抽出してもよい。例えば、１〜１５番目のフレーム画像が存在する場合、７番目のフレーム画像が基準フレーム画像として指定された場合、７番目のフレーム画像から１番目のフレーム画像まで、順次さかのぼって横隔膜位置候補を抽出し、７番目のフレーム画像から１５番目のフレーム画像まで、順次横隔膜位置候補を抽出することとしてもよい。
【０１３８】
また、上記実施形態においては、動態画像の全フレーム画像の横隔膜位置を計測する場合を例にとり説明したが、動態画像のフレーム画像のうち、位置計測を行うフレーム画像の範囲が指定された場合は、その範囲で位置計測を行うこととしてもよい。
【０１３９】
また、上記実施形態においては、ユーザーが横隔膜の速度、加速度の閾値を設定することとしたが、制御部３１により自動的に設定してもよい。また、本実施形態では、速度、加速度の上限を設定する場合を示しているが、上限だけでなく下限も設定することとしてもよい。
自動的に横隔膜の速度、加速度の閾値を設定する場合、例えば、以下の方法により設定を行うことができる。
例えば、さまざまな胸部の動態画像の横隔膜の動きから算出した速度及び加速度の範囲の最大値を上限の閾値、最小値を下限の閾値に設定する。或いは、さまざまな胸部の動態画像の横隔膜の動きから算出した速度及び加速度が取った値の範囲にマージンを与えた範囲を横隔膜が取り得る速度及び加速度の範囲とし、その最大値を上限の閾値、最小値を下限の閾値に設定することとしてもよい。この閾値は、動態画像を収集する過程で速度及び加速度が取った範囲をアップデートし、その結果を閾値に反映することで、より精度のよい横隔膜の計測結果を得ることが可能となる。
【０１４０】
また、上記実施形態においては、動態画像においてユーザーが指定した計測対象点を中心とするｗ画素×ｈ画素の領域をテンプレート画像としたが、横隔膜を模した幾何学図形や、横隔膜の画像を別のサンプルデータから抽出した画像や、多数のＸ線画像から作成した一般的な横隔膜の画像等をテンプレート画像としてもよい。
【０１４１】
また、上記実施形態においては、動態画像の複数のフレーム画像のそれぞれにおける所定範囲内の各位置と横隔膜（所定の構造物）との類似度を示す評価値をテンプレートマッチングにより求める場合を例にとり説明したが、評価値の算出手法はこれに限定されない。例えば、計測対象点を中心とするｗ画素×ｈ画素の領域を二値化したときの白の画素と黒の画素の比の値を横隔膜の特徴を表す基準値とし、横隔膜を探索する所定範囲内の各位置（各画素を中心とするｗ画素×ｈ画素の領域）について算出した比の値と基準値との差分値等を各位置と横隔膜との類似度を表す評価値として用いることとしてもよい。
【０１４２】
また、上記実施形態においては、ルート候補を逐次限定していくこととしたが、例えば、ルート保持部３２１として大容量の記憶容量が確保できるような場合等においては、ルート候補の限定を省略しても構わない。
【０１４３】
また、上記の説明では、本発明に係るプログラムのコンピュータ読み取り可能な媒体としてハードディスクや半導体の不揮発性メモリー等を使用した例を開示したが、この例に限定されない。その他のコンピュータ読み取り可能な媒体として、ＣＤ-ＲＯＭ等の可
搬型記録媒体を適用することが可能である。また、本発明に係るプログラムのデータを通信回線を介して提供する媒体として、キャリアウエーブ(搬送波)も適用される。
【０１４４】
その他、位置計測システム１００を構成する各装置の細部構成及び細部動作に関しても、本発明の趣旨を逸脱することのない範囲で適宜変更可能である。
【符号の説明】
【０１４５】
１００  位置計測システム
１  撮影装置
１１  放射線源
１２  放射線照射制御装置
１３  放射線検出部
１４  読取制御装置
２  撮影用コンソール
２１  制御部
２２  記憶部
２３  操作部
２４  表示部
２５  通信部
２６  バス
３  診断用コンソール
３１  制御部
３２  記憶部
３２１  ルート保持部
３３  操作部
３４  表示部
３５  通信部
３６  バス
【図１】
2018068400.tif 000003
【図２】
2018068400.tif 000004
【図３】
2018068400.tif 000005
【図４】
2018068400.tif 000006
【図５】
2018068400.tif 000007
【図６】
2018068400.tif 000008
【図７】
2018068400.tif 000009
【図８】
2018068400.tif 000010
【図９】
2018068400.tif 000011
【図１０】
2018068400.tif 000012
【図１１】
2018068400.tif 000013