(57)【特許請求の範囲】
【請求項１】
連続する２心拍分の時系列のＸ線画像を取得する取得部と、
前記連続する２心拍分の時系列のＸ線画像における造影剤の濃度の経時的な変化に基づいて、前記連続する２心拍のうち一方の心拍を特定する特定部と、
特定された前記一方の心拍に対応する時系列のＸ線画像における造影剤の濃度の経時的な変化に基づくパラメータを生成するパラメータ生成部と、
を備える医用画像処理装置。
【請求項２】
前記パラメータ生成部は、前記時系列のＸ線画像における各位置について、特定された前記一方の心拍に対応する時系列のＸ線画像における造影剤の濃度の経時的な変化に基づくパラメータを生成し、生成した当該パラメータをそれぞれ対応する位置に割り当てたパラメータマップを生成する、請求項１に記載の医用画像処理装置。
【請求項３】
前記パラメータ生成部は、特定された前記一方の心拍に対応する時系列のＸ線画像に対して設定された領域内の画素の統計値の経時的な変化に基づくパラメータを生成する、請求項１に記載の医用画像処理装置。
【請求項４】
前記取得部は、造影剤注入後に時系列に収集されたＸ線画像から、前記造影剤が撮影対象部位に到達後の連続する２心拍分の時系列のＸ線画像を取得する、請求項１〜３のいずれか一つに記載の医用画像処理装置。
【請求項５】
前記特定部は、前記Ｘ線画像に設定された関心領域における時間濃度曲線を導出し、当該時間濃度曲線において連続する複数フレーム間での変化量が最大となる期間を含んだ心拍を前記一方の心拍として特定する、請求項１〜４のいずれか一つに記載の医用画像処理装置。
【請求項６】
前記特定部は、心拍の影響を受ける臓器に設定された前記関心領域中の画素値の合計値或いは平均値の時間濃度曲線を導出して前記一方の心拍を特定する、請求項５に記載の医用画像処理装置。
【請求項７】
前記特定部は、心拍の影響を受ける臓器に設定された前記関心領域中の各画素の変化を示す時間濃度曲線をそれぞれ導出し、導出した各時間濃度曲線を加算した時間濃度曲線を導出して前記一方の心拍を特定する、請求項５に記載の医用画像処理装置。
【請求項８】
前記特定部は、前記関心領域として左肺に設定された第１の関心領域及び右肺に設定された第２の関心領域それぞれに対応する領域中の画素値の合計値或いは平均値の時間濃度曲線を導出して前記一方の心拍を特定する、請求項５に記載の医用画像処理装置。
【請求項９】
前記特定部は、前記関心領域として左肺に設定された第１の関心領域及び右肺に設定された第２の関心領域それぞれに対応する領域中の各画素の変化を示す時間濃度曲線をそれぞれ導出し、導出した各時間濃度曲線を加算した時間濃度曲線を導出して前記一方の心拍を特定する、請求項５に記載の医用画像処理装置。
【請求項１０】
前記特定部は、前記関心領域として左肺に設定された第１の関心領域及び右肺に設定された第２の関心領域それぞれに対応する領域中の各画素の変化を示す時間濃度曲線をそれぞれ導出して前記一方の心拍を特定する、請求項５に記載の医用画像処理装置。
【請求項１１】
前記特定部は、左右の肺を含んだ前記関心領域中の画素値の合計値或いは平均値の時間濃度曲線を導出して前記一方の心拍を特定する、請求項５に記載の医用画像処理装置。
【請求項１２】
前記特定部は、左右の肺を含んだ前記関心領域中の各画素の変化を示す時間濃度曲線をそれぞれ導出し、導出した各時間濃度曲線を加算した時間濃度曲線を導出して前記一方の心拍を特定する、請求項５に記載の医用画像処理装置。
【請求項１３】
前記パラメータ生成部は、前記パラメータに基づいて、左右の肺のパラメータの比を更に算出する、請求項１〜１２のいずれか一つに記載の医用画像処理装置。
【請求項１４】
Ｘ線を照射するＸ線管と、
前記Ｘ線管から照射され被検体を透過したＸ線を検出して検出信号を出力するＸ線検出器と、
前記検出信号に基づいてＸ線画像を時系列に沿って生成する画像生成部と、
連続する２心拍分の時系列のＸ線画像を取得する取得部と、
前記連続する２心拍分の時系列のＸ線画像における造影剤の濃度の経時的な変化に基づいて、前記連続する２心拍のうち一方の心拍を特定する特定部と、
特定された前記一方の心拍に対応する時系列のＸ線画像における造影剤の濃度の経時的な変化に基づくパラメータを生成するパラメータ生成部と、
を備えるＸ線診断装置。
【請求項１５】
連続する２心拍分の時系列のＸ線画像を取得する取得部と、
前記連続する２心拍のうち一方の心拍を特定する特定部と、
特定された前記一方の心拍に対応する時系列のＸ線画像における造影剤の濃度の経時的な変化に基づくパラメータを生成するパラメータ生成部と、
を備え、
前記特定部は、前記Ｘ線画像に設定された関心領域における時間濃度曲線を導出し、当該時間濃度曲線において連続する複数フレーム間での変化量が最大となる期間を含んだ心拍を前記一方の心拍として特定する医用画像処理装置。
【請求項１６】
連続する２心拍分の時系列のＸ線画像を取得する取得部と、
前記連続する２心拍分の時系列のＸ線画像におけるカテーテルの向きに応じて、前記連続する２心拍のうち一方の心拍を特定する特定部と、
特定された前記一方の心拍に対応する時系列のＸ線画像における造影剤の濃度の経時的な変化に基づくパラメータを生成するパラメータ生成部と、
を備える医用画像処理装置。
【請求項１７】
連続する２心拍分の時系列のＸ線画像を取得する取得部と、
前記連続する２心拍のうち一方の心拍を特定する特定部と、
特定された前記一方の心拍に対応する時系列のＸ線画像における造影剤の濃度の経時的な変化に基づくパラメータを生成し、前記パラメータに基づいて、左右の肺のパラメータの比を算出するパラメータ生成部と、
を備える医用画像処理装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明の実施形態は、医用画像処理装置及びＸ線診断装置に関する。
【背景技術】
【０００２】
血管に造影剤を投与してＸ線二次元投影画像を取得し、投影画像における時間濃度曲線から得られるパラメータを新たな画像として表示するパラメトリックイメージング（Parametric  imaging）技術が知られている。
【先行技術文献】
【特許文献】
【０００３】
【特許文献１】特開２０１４−１２１３３号公報
【発明の概要】
【発明が解決しようとする課題】
【０００４】
本発明が解決しようとする課題は、肺血管造影において血流を正確に測定することである。
【課題を解決するための手段】
【０００５】
実施形態の医用画像処理装置は、取得部と、特定部と、パラメータ生成部とを備える。取得部は、連続する２心拍分の時系列のＸ線画像を取得する。特定部は、前記連続する２心拍のうち一方の心拍を特定する。パラメータ生成部は、特定された前記一方の心拍に対応する期間における造影剤の濃度の経時的な変化に基づくパラメータを生成する。
【図面の簡単な説明】
【０００６】
【図１】図１は、第１の実施形態に係る医用画像処理システムの構成の一例を示す図である。
【図２】図２は、肺血管造影を説明するための図である。
【図３】図３は、第１の実施形態に係る医用画像処理装置による処理手順を示すフローチャートである。
【図４】図４は、第１の実施形態を説明するための図である。
【図５】図５は、第１の実施形態を説明するための図である。
【図６】図６は、第１の実施形態を説明するための図である。
【図７】図７は、第１の実施形態の変形例を説明するための図である。
【図８】図８は、第２の実施形態を説明するための図である。
【図９】図９は、第２の実施形態を説明するための図である。
【図１０】図１０は、第３の実施形態を説明するための図である。
【発明を実施するための形態】
【０００７】
以下、図面を参照して、実施形態に係る医用画像処理装置及びＸ線診断装置を説明する。また、以下では、医用画像処理装置及びＸ線診断装置を有する医用画像処理システムについて説明する。なお、実施形態は、以下の実施形態に限られるものではない。また、一つの実施形態に記載した内容は、原則として他の実施形態にも同様に適用される。例えば、医用画像処理システムにＰＡＣＳ（Picture  Archiving  and  Communication  System）が導入されている場合、各装置は、医用画像データに付帯情報を付与したＤＩＣＯＭ（Digital  Imaging  and  Communications  in  Medicine）形式の医用画像データ等を相互に送受信する。ここで、付帯情報には、例えば、患者を識別する患者ＩＤ（Identifier）、検査を識別する検査ＩＤ、各装置を識別する装置ＩＤ、各装置による１回の撮影を識別するシリーズＩＤなどが含まれる。
【０００８】
（第１の実施形態）
図１は、第１の実施形態に係る医用画像処理システム１の構成の一例を示す図である。図１に示すように、第１の実施形態に係る医用画像処理システム１は、Ｘ線診断装置１００と医用画像処理装置２００とを有する。
【０００９】
Ｘ線診断装置１００は、高電圧発生器１１と、Ｘ線管１２と、コリメータ１３と、天板１４と、Ｃアーム１５と、Ｘ線検出器１６とを備える。また、第１の実施形態に係るＸ線診断装置１００は、Ｃアーム回転・移動機構１７と、天板移動機構１８と、Ｃアーム・天板機構制御回路１９と、絞り制御回路２０と、処理回路２１と、入力インターフェース２２と、ディスプレイ２３とを備える。また、第１の実施形態に係るＸ線診断装置１００は、画像データ生成回路２４と、記憶回路２５と、画像処理回路２６と、通信インターフェース２７とを備える。また、Ｘ線診断装置１００は、インジェクター３０及び心電計４０と接続される。そして、Ｘ線診断装置１００は、図１に示すように、各回路が相互に接続され、各回路間で種々の電気信号を送受信したり、インジェクター３０と電気信号を送受信したりする。
【００１０】
インジェクター３０は、被検体Ｐに挿入されたカテーテルから造影剤を注入するための装置である。ここで、インジェクター３０からの造影剤注入は、後述する処理回路２１を介して受信した注入指示に従って実行される。具体的には、インジェクター３０は、後述する処理回路２１から受信する造影剤の注入開始指示や、注入停止指示、さらに、注入速度などを含む造影剤注入条件に応じた造影剤注入を実行する。なお、インジェクター３０は、操作者が直接インジェクター３０に対して入力した注入指示に従って注入開始や、注入停止を実行することも可能である。
【００１１】
心電計４０は、図示しない端子が取り付けられた被検体Ｐの心電波形（ＥＣＧ：Electrocardiogram）を取得し、取得した心電波形を、時間情報とともに、処理回路２１、画像データ生成回路２４および画像処理回路２６に送信する。
【００１２】
図１に示すＸ線診断装置１００においては、各処理機能がコンピュータによって実行可能なプログラムの形態で記憶回路２５へ記憶されている。Ｃアーム・天板機構制御回路１９、絞り制御回路２０、処理回路２１、画像データ生成回路２４、及び、画像処理回路２６は、記憶回路２５からプログラムを読み出して実行することで各プログラムに対応する機能を実現するプロセッサである。換言すると、各プログラムを読み出した状態の各回路は、読み出したプログラムに対応する機能を有することとなる。
【００１３】
高電圧発生器１１は、処理回路２１による制御の下、高電圧を発生し、発生した高電圧をＸ線管１２に供給する。Ｘ線管１２は、高電圧発生器１１から供給される高電圧を用いて、Ｘ線を発生する。
【００１４】
コリメータ１３は、絞り制御回路２０による制御の下、Ｘ線管１２が発生したＸ線を、被検体Ｐの関心領域に対して選択的に照射されるように絞り込む。例えば、コリメータ１３は、スライド可能な４枚の絞り羽根を有する。コリメータ１３は、絞り制御回路２０による制御の下、これらの絞り羽根をスライドさせることで、Ｘ線管１２が発生したＸ線を絞り込んで被検体Ｐに照射させる。天板１４は、被検体Ｐを載せるベッドであり、図示しない寝台の上に配置される。なお、被検体Ｐは、Ｘ線診断装置１００に含まれない。
【００１５】
Ｘ線検出器１６は、被検体Ｐを透過したＸ線を検出する。例えば、Ｘ線検出器１６は、マトリックス状に配列された検出素子を有する。各検出素子は、被検体Ｐを透過したＸ線を電気信号に変換して蓄積し、蓄積した電気信号を画像データ生成回路２４に送信する。
【００１６】
Ｃアーム１５は、Ｘ線管１２、コリメータ１３及びＸ線検出器１６を保持する。Ｘ線管１２及びコリメータ１３とＸ線検出器１６とは、Ｃアーム１５により被検体Ｐを挟んで対向するように配置される。なお、図１では、Ｘ線診断装置１００がシングルプレーンの場合を例に挙げて説明しているが、実施形態はこれに限定されるものではなく、バイプレーンの場合であってもよい。
【００１７】
Ｃアーム回転・移動機構１７は、Ｃアーム１５を回転及び移動させるための機構であり、天板移動機構１８は、天板１４を移動させるための機構である。Ｃアーム・天板機構制御回路１９は、処理回路２１による制御の下、Ｃアーム回転・移動機構１７及び天板移動機構１８を制御することで、Ｃアーム１５の回転や移動、天板１４の移動を調整する。絞り制御回路２０は、処理回路２１による制御の下、コリメータ１３が有する絞り羽根の開度を調整することで、被検体Ｐに対して照射されるＸ線の照射範囲を制御する。
【００１８】
画像データ生成回路２４は、Ｘ線検出器１６によってＸ線から変換された電気信号を用いて画像データを生成し、生成した画像データを記憶回路２５に格納する。例えば、画像データ生成回路２４は、Ｘ線検出器１６から受信した電気信号に対して、電流・電圧変換やＡ（Analog）／Ｄ（Digital）変換、パラレル・シリアル変換を行い、画像データを生成する。一例を挙げると、画像データ生成回路２４は、造影剤が注入されていない状態で撮像された画像データ（マスク画像）及び造影剤が注入された状態で撮像された画像データ（コントラスト画像）を生成する。そして、画像データ生成回路２４は、生成したマスク画像及びコントラスト画像を記憶回路２５に格納する。
【００１９】
記憶回路２５は、画像データ生成回路２４によって生成された画像データを受け付けて記憶する。例えば、記憶回路２５は、造影剤が投与される前後の被検体Ｐの画像データを記憶する。また、記憶回路２５は、後述する差分画像を記憶する。
【００２０】
画像処理回路２６は、記憶回路２５が記憶する画像データに対して各種画像処理を行う。例えば、画像処理回路２６は、記憶回路２５が記憶するマスク画像とコントラスト画像とを読み出し、サブトラクション（Ｌｏｇサブ）することで差分画像を生成する。
【００２１】
なお、画像処理回路２６は、造影剤投与直前の１フレームをマスク画像として用いることで、体動による位置合わせ（レジストレーション）のミスを最小限に抑えることが可能である。また、画像処理回路２６は、移動平均（平滑化）フィルタ、ガウシアンフィルタ、メディアンフィルタなどの画像処理フィルタによるノイズ低減処理を実行することも可能である。すなわち、画像処理回路２６は、造影剤を用いて経時的に撮影された複数のＸ線画像群それぞれに対して、位置ずれ補正及びノイズ除去を含む前処理を実行することができる。
【００２２】
入力インターフェース２２は、所定の領域（例えば、差分画像内で補正処理の対象領域）などの設定などを行うためのトラックボール、スイッチボタン、マウス、キーボード等によって実現される。入力インターフェース２２は、処理回路２１に接続されており、操作者から受け取った入力操作を電気信号へ変換し処理回路２１へと出力する。
【００２３】
ディスプレイ２３は、操作者の指示を受け付けるためのＧＵＩ（Graphical  User  Interface）や、画像処理回路２６によって生成された差分画像、パラメトリックイメージングによるカラー画像などを表示する。
【００２４】
通信インターフェース２７は、処理回路２１に接続され、ネットワークを介して接続された医用画像処理装置２００等の他装置との間で行われる各種データの伝送及び通信を制御する。例えば、通信インターフェース２７は、ネットワークカードやネットワークアダプタ、ＮＩＣ（Network  Interface  Controller）等によって実現される。図１では、通信インターフェース２７は、記憶回路２５に記憶されているＸ線画像を医用画像処理装置２００に送信する。
【００２５】
処理回路２１は、Ｘ線診断装置１００全体の動作を制御する。例えば、処理回路２１は、装置全体を制御するための制御機能に対応するプログラムを記憶回路２５から読み出して実行することにより、種々の処理を実行する。例えば、処理回路２１は、入力インターフェース２２から転送された操作者の指示に従って高電圧発生器１１を制御し、Ｘ線管１２に供給する電圧を調整することで、被検体Ｐに対して照射されるＸ線量やＯＮ／ＯＦＦを制御する。また、例えば、処理回路２１は、操作者の指示に従ってＣアーム・天板機構制御回路１９を制御し、Ｃアーム１５の回転や移動、天板１４の移動を調整する。また、例えば、処理回路２１は、操作者の指示に従って絞り制御回路２０を制御し、コリメータ１３が有する絞り羽根の開度を調整することで、被検体Ｐに対して照射されるＸ線の照射範囲を制御する。
【００２６】
また、処理回路２１は、操作者の指示に従って、画像データ生成回路２４による画像データ生成処理や、画像処理回路２６による画像処理、あるいは解析処理などを制御する。また、処理回路２１は、操作者の指示を受け付けるためのＧＵＩや記憶回路２５が記憶する画像などを、ディスプレイ２３に表示するように制御する。また、処理回路２１は、インジェクター３０に対して、造影剤注入開始及び終了の信号を送信することで、造影剤の注入タイミングを制御する。また、処理回路２１は、心電計４０から受信するＥＣＧの時間情報と画像データ生成回路２４によって生成された画像データの時間情報とを対応付ける。
【００２７】
医用画像処理装置２００は、入力インターフェース２０１と、ディスプレイ２０２と、記憶回路２０３と、処理回路２０４と、通信インターフェース２０５とを有する。
【００２８】
通信インターフェース２０５は、処理回路２０４に接続され、ネットワークを介して接続された各種の医用画像診断装置等の他装置との間で行われる各種データの伝送及び通信を制御する。例えば、通信インターフェース２０５は、ネットワークカードやネットワークアダプタ、ＮＩＣ（Network  Interface  Controller）等によって実現される。図１では、通信インターフェース２０５は、Ｘ線診断装置１００からＸ線画像を受信し、受信したＸ線画像を処理回路２０４に出力する。
【００２９】
記憶回路２０３は、処理回路２０４に接続され、各種データを記憶する。例えば、記憶回路２０３は、ＲＡＭ（Random  Access  Memory）、フラッシュメモリ等の半導体メモリ素子や、ハードディスク、光ディスク等によって実現される。記憶回路２０３は、Ｘ線診断装置１００から受信したＸ線画像を記憶する。例えば、記憶回路２０３は、Ｘ線画像として、マスク画像及びコントラスト画像を記憶する。或いは、記憶回路２０３は、Ｘ線画像として、差分画像を記憶する。また、或いは、記憶回路２０３は、Ｘ線画像として、マスク画像、コントラスト画像及び差分画像を記憶する。また、記憶回路２０３は、処理回路２０４による処理結果を記憶する。
【００３０】
入力インターフェース２０１は、処理回路２０４に接続され、操作者から受け付けた入力操作を電気信号に変換して処理回路２０４に出力する。例えば、入力インターフェース２０１は、トラックボール、スイッチボタン、マウス、キーボード、タッチパネル等によって実現される。
【００３１】
ディスプレイ２０２は、処理回路２０４に接続され、処理回路２０４から出力される各種情報及び各種画像データを表示する。例えば、ディスプレイ２０２は、液晶モニタやＣＲＴ（Cathode  Ray  Tube）モニタ、タッチパネル等によって実現される。
【００３２】
処理回路２０４は、入力インターフェース２０１を介して操作者から受け付けた入力操作に応じて、医用画像処理装置２００が有する各構成要素を制御する。例えば、処理回路２０４は、プロセッサによって実現される。処理回路２０４は、通信インターフェース２０５から出力されるＸ線画像を記憶回路２０３に記憶させる。また、処理回路２０４は、記憶回路２０３からＸ線画像データを読み出し、ディスプレイ２０２に表示する。
【００３３】
また、処理回路２０４は、図１に示すように、取得機能２０４ａと、特定機能２０４ｂと、パラメータ生成機能２０４ｃとを実行する。ここで、例えば、図１に示す処理回路２０４の構成要素である取得機能２０４ａと、特定機能２０４ｂと、パラメータ生成機能２０４ｃとが実行する各処理機能は、コンピュータによって実行可能なプログラムの形態で記憶回路２０３に記録されている。処理回路２０４は、各プログラムを記憶回路２０３から読み出し、実行することで各プログラムに対応する機能を実現するプロセッサである。換言すると、各プログラムを読み出した状態の処理回路２０４は、図１の処理回路２０４内に示された各機能を有することとなる。
【００３４】
以上、医用画像処理システム１の構成の一例について説明した。かかる構成のもと、本実施形態に係る医用画像処理システム１では、医用画像処理装置２００は、肺血管造影にパラメトリックイメージングを適用する。例えば、医用画像処理装置２００は、心臓疾患の手術直後の判断や心臓疾患の手術前における手術計画に、パラメトリックイメージングを利用する。
【００３５】
例えば、肺血管造影は、肺動脈にカテーテルを留置し、造影剤を投与する。ここで、造影剤は、血液が心臓を出て肺を通り心臓に戻る肺循環に沿って流れる。すなわち、心拍に同期して右心室から肺動脈へと血液が流れる際に、造影剤も左右の肺動脈に分かれて流れる。そして、肺動脈に流れた造影剤は、次第に肺胞に達し、その後肺静脈を経由して左心房に戻ってくる。なお、肺循環に要する時間は成人で４−５秒と言われており、小児の場合の肺循環は、４−５秒より短い。また、正常な心臓では、左右の肺に入る血液の量は、略１対１になる。
【００３６】
ここで、心臓疾患には、完全大血管転位症（ＴＧＡ：Transposition  of  the  great  arteries）という指定難病がある。正常な心臓では大動脈は左心室から、肺動脈は右心室から起始するが、ＴＧＡでは、大動脈と肺動脈が転位し、大動脈が右心室から、肺動脈が左心室から起始する。上述したように、正常な心臓では、左右の肺に入る血液の量は、略１対１になるが、ＴＧＡでは、左右の肺に入る血液の量が１対１にはならない場合がある。このようなことから、例えば、ＴＧＡの手術直後やＴＧＡの手術前における手術計画の判断に、肺血管造影にパラメトリックイメージング技術を適用して肺の組織に入っている造影剤量の左右の比を観察することが考えられる。
【００３７】
しかしながら、肺血管造影にパラメトリックイメージング技術を適用する場合、造影剤が肺の微小循環系内に停留する時間が１秒以内と短く、また心拍運動の影響を受けるため、単純に時間濃度曲線をパラメータ化しただけでは正しく血流を測定できない。
【００３８】
図２は、肺血管造影を説明するための図である。図２では、造影剤存在下の肺血管を時系列に沿って撮影した複数のＸ線画像を示す。また、図２では、上段から下段にかけて時間経過を示し、また、各段に置いて左側から右側にかけて時間経過を示す。また、図２では、造影剤投与直後から１秒までのＸ線画像を示す。図２に示すように、肺動脈から肺胞への造影剤流入量を、画像の時間濃度曲線から推定する場合、造影剤投与直後から１秒程度は、肺胞での造影剤は増えていく一方である。すなわち、図２の上段に示すＸ線画像では、肺動脈に造影剤が蓄積していく様子を示している。そして、図２の中段左端に示すＸ線画像は、肺動脈から肺胞へ造影剤が流れ込むタイミングのＸ線画像である。図２の中段及び下段に示すＸ線画像では、時間の経過とともに肺胞が染まっていく様子を示している。
【００３９】
その後、１秒以降では、肺胞に到達した造影剤は流れ去っていく。よって、造影剤投与直後から１秒程度は増えている量のみが観測されるが、１秒以降は増えていく量と減っていく量の差分が観測されることになる。また、肺胞から流れ出た造影剤は肺静脈を通るため二次元画像上で重なってしまうため、観察される画素値は肺胞の画素値に加えて肺静脈の画素値を含むことになる。したがって測定をシンプルにするためには、造影剤投与後１秒程度のデータのみを用いることが望ましい。
【００４０】
更に、心臓は拍動する。肺動脈に注入した造影剤は、徐々に肺胞に向かって流れるのではなく、心拍運動に応じて流れる。すなわち、心臓の拡張期には、造影剤は肺動脈に蓄積され、心臓の収縮期に造影剤が肺動脈から肺胞へと一斉に流れ出す。したがって、任意のタイミングで血流を測定しても正確ではなく、心拍運動を考慮した所定のタイミングで血流を測定する必要がある。これらの理由により、心拍運動にともなって観測されるＸ線画像における濃度の時間濃度曲線を観察し、所定のタイミングのデータのみから計算をおこなうことが望ましいと考えられる。例えば、比較例として最初の心拍（第１の心拍とも言う）のデータのみから計算をおこなう場合について説明する。
【００４１】
比較例において、造影剤を投与するタイミングが心拍運動のタイミングとずれてしまった場合、造影剤投与途中で肺に造影剤が流れ込む。このような流れ込みがあると、造影剤の投与量が少なくなるので、第１の心拍のデータでは正確に血流を測定することができなくなる。この場合、第１の心拍による流れを測定するのではなく、第２の心拍による流れを測定した方が正確な測定となる。
【００４２】
なお、第３の心拍を使うことは好ましくない。第３の心拍では、造影剤が肺胞から肺静脈に流れ出すフェーズになってしまうことがある。また肺静脈が二次元画像上では重なってしまうため、測定される信号値に余分な要素が入り込んでしまう。
【００４３】
このように、比較例では、第１の心拍及び第２の心拍のどちらのタイミングで血流を測定すれば良いのかを特定できず、正確に血流を測定することが困難であった。そこで、第１の実施形態では、医用画像処理装置２００は、肺血管造影において血流を正確に測定するために、以下に示すパラメータ生成処理を実行する。すなわち、医用画像処理装置２００は、時間濃度曲線を参照し、第１の心拍のデータ及び第２の心拍のデータを取得し、取得したデータそれぞれの特性を調べ、第１の心拍のデータもしくは第２の心拍のデータのいずれか一方を特定する。そして、医用画像処理装置２００は、特定した心拍のデータを用いてパラメータを生成する。パラメータ生成処理は、処理回路２０４が、取得機能２０４ａ、特定機能２０４ｂ及びパラメータ生成機能２０４ｃを実行することにより実現される。以下では、図３から図７を用いて、パラメータ生成処理について説明する。図３は、第１の実施形態に係る医用画像処理装置２００による処理手順を示すフローチャートであり、図４から図６は、第１の実施形態を説明するための図である。
【００４４】
図３では、医用画像処理装置２００の処理回路２０４の動作を説明するフローチャートを示し、各構成要素がフローチャートのどのステップに対応するかを説明する。ステップＳ１は、取得機能２０４ａに対応するステップである。処理回路２０４が記憶回路２０３から取得機能２０４ａに対応する所定のプログラムを呼び出し実行することにより、取得機能２０４ａが実現されるステップである。ステップＳ１では、取得機能２０４ａは、連続する２心拍分の時系列のＸ線画像を取得する。なお、以下では、取得機能２０４ａは、Ｘ線画像として時系列の差分画像を取得する場合について説明するが、実施形態はこれに限定されるものではない。例えば、取得機能２０４ａは、Ｘ線画像としてマスク画像と、時系列のコントラスト画像とを取得し、取得したマスク画像と、時系列のコントラスト画像とをサブトラクションすることで、時系列の差分画像を取得してもよい。
【００４５】
ここで、取得機能２０４ａは、造影剤注入後に時系列に収集されたＸ線画像から、造影剤が撮影対象部位に到達後の連続する２心拍分の時系列のＸ線画像を取得する。取得機能２０４ａは、以下に示す第１の方法から第３の方法のうちいずれかの方法を用いて、連続する２心拍分の時系列のＸ線画像を取得する。
【００４６】
第１の方法では、取得機能２０４ａは、差分画像の画素値に基づいて連続する２心拍分の時系列のＸ線画像を取得する。例えば、取得機能２０４ａは、時系列に沿って生成された差分画像を記憶回路２０３から取得する。そして、取得機能２０４ａは、各差分画像の画素値を算出し、算出した画素値が第１の閾値以上となった差分画像を特定する。また、取得機能２０４ａは、画素値が第１の閾値以上となった差分画像の中から、算出した画素値が時系列順で最初に極大値となった差分画像を特定する。そして、取得機能２０４ａは、画素値が時系列順で最初に第１の閾値以上となった差分画像から画素値が時系列順で最初に極大値となった差分画像までの間のＸ線画像を第１の心拍のＸ線画像とする。同様にして、取得機能２０４ａは、算出した画素値が第２の閾値以上となった差分画像を特定する。また、取得機能２０４ａは、画素値が第２の閾値以上となった差分画像の中から、算出した画素値が時系列順で最初に極大値となった差分画像を特定する。そして、取得機能２０４ａは、画素値が時系列順で最初に第２の閾値以上となった差分画像から画素値が時系列順で最初に極大値となった差分画像までの間のＸ線画像を第２の心拍のＸ線画像とする。
【００４７】
第２の方法では、取得機能２０４ａは、心電波形に基づいて、連続する２心拍分の時系列のＸ線画像を特定する。かかる場合、取得機能２０４ａは、造影剤を注入した時間を基準時間とし、この基準時間に対応するＸ線画像から心電波形の２心拍に対応するＸ線画像を特定し、連続する２心拍分の時系列のＸ線画像を取得する。
【００４８】
第３の方法では、取得機能２０４ａは、操作者からＸ線画像の指定を受け付け、指定されたＸ線画像を、連続する２心拍分の時系列のＸ線画像として取得する。なお、取得機能２０４ａは、第１の方法から第３の方法を任意に組み合わせて、連続する２心拍分の時系列のＸ線画像を取得してもよい。
【００４９】
ステップＳ２からステップＳ４は、特定機能２０４ｂに対応するステップである。処理回路２０４が記憶回路２０３から特定機能２０４ｂに対応する所定のプログラムを呼び出し実行することにより、特定機能２０４ｂが実現されるステップである。ステップＳ２では、特定機能２０４ｂは、関心領域（ＲＯＩ：Region  of  Interest）の設定を受け付ける。
【００５０】
ここで、操作者は、ステップＳ１で取得したＸ線画像のうち任意のＸ線画像を選択し、選択したＸ線画像において関心領域を設定する。ここで、操作者は、例えば、入力インターフェース２０１を介して、右の肺と左の肺とにそれぞれ関心領域を設定する。これにより、特定機能２０４ｂは、例えば図４に示すように、右の肺と左の肺とにそれぞれ関心領域の設定を受け付ける。図４では、特定機能２０４ｂが、右の肺に関心領域５１の設定を受け付け、左の肺に関心領域５２の設定を受け付けた場合を示す。
【００５１】
図３に戻る。ステップＳ３では、特定機能２０４ｂは、時間濃度曲線（Time  Density  Curves：ＴＤＣ）を導出する。例えば、特定機能２０４ｂは、関心領域として左肺に設定された第１の関心領域及び右肺に設定された第２の関心領域それぞれに対応する領域中の画素値の合計値或いは平均値の時間濃度曲線を導出する。図５では、横軸がフレーム数を示し、縦軸が画素値を示す。また、図５では、各関心領域中における画素値の合計値の時間濃度曲線を導出する場合について説明する。特定機能２０４ｂは、図５に示すように、右の肺に設定された関心領域５１内の画素値の合計値の時間濃度曲線６１と、左の肺に設定された関心領域５２内の画素値の合計値の時間濃度曲線６２とを導出する。
【００５２】
ステップＳ４では、特定機能２０４ｂは、連続する２心拍のうち一方の心拍を特定する。例えば、特定機能２０４ｂは、Ｘ線画像に設定された関心領域における時間濃度曲線を導出し、当該時間濃度曲線において連続する複数フレーム間での変化量が最大となる期間を含んだ心拍を一方の心拍として特定する。より具体的には、特定機能２０４ｂは、ステップＳ３で導出した２つの時間濃度曲線を加算する。そして、特定機能２０４ｂは、加算後の時間濃度曲線において、第１フレーム目から順に最終フレーム（第ｎフレームとする）まで、画素値の変化量を算出する。
【００５３】
そして、特定機能２０４ｂは、算出した変化量の中から最大値を特定する。図６では、横軸がフレーム数を示し、縦軸が画素値を示す。また、図６では、右の肺と左の肺とにそれぞれの関心領域から導出した時間濃度曲線を加算した時間濃度曲線６３を導出した場合について説明する。特定機能２０４ｂは、例えば、図６において、第ｉフレームから第ｉ＋６フレームまでの期間における画素値の変化量が最大値であると特定する。そして、特定機能２０４ｂは、第１の心拍及び第２の心拍のうち、特定した期間に対応する心拍を特定する。
【００５４】
ステップＳ５からステップＳ９は、パラメータ生成機能２０４ｃに対応するステップである。処理回路２０４が記憶回路２０３からパラメータ生成機能２０４ｃに対応する所定のプログラムを呼び出し実行することにより、パラメータ生成機能２０４ｃが実現されるステップである。ステップＳ５では、パラメータ生成機能２０４ｃは、パラメータを生成する。
【００５５】
例えば、パラメータ生成機能２０４ｃは、時系列のＸ線画像における各位置について、特定された一方の心拍に対応する期間における造影剤の濃度の経時的な変化に基づくパラメータを生成し、生成した当該パラメータをそれぞれ対応する位置に割り当てたパラメータマップを生成する。より具体的には、パラメータ生成機能２０４ｃは、特定された一方の心拍に対応する時系列の差分画像の画素ごとの造影剤の時間濃度曲線を用いて各種パラメータ値を算出する。
【００５６】
ステップＳ６では、パラメータ生成機能２０４ｃは、パラメータマップを生成する。例えば、パラメータ生成機能２０４ｃは、ステップＳ５で算出したパラメータ値に基づいてパラメトリックイメージングを実行する。
【００５７】
ここで、パラメータ生成機能２０４ｃは、観察者である操作者が所望する血流情報に応じたパラメータ値を画素ごとに算出し、算出したパラメータ値に対応する色で画像上の各画素をカラー化した画像をパラメータマップとして生成する。例えば、パラメータ生成機能２０４ｃは、時間濃度曲線を用いて画素値の時系列変化がピークあるいは特定の値となる時間を流入時間として算出し、算出した時間にカラー情報を割り当てる。なお、パラメータ生成機能２０４ｃは、既存の技術を適用してパラメータマップを生成する。また、パラメータ生成機能２０４ｃは、静止画としてパラメータマップを生成してもよいし、動画としてパラメータマップを生成してもよい。
【００５８】
そして、ステップＳ７では、パラメータ生成機能２０４ｃは、ステップＳ６で生成したパラメータマップをディスプレイ２０２に表示させる。ここで、パラメータ生成機能２０４ｃは、静止画としてパラメータマップを生成した場合、時系列のパラメータマップを所定枚数分並べてディスプレイ２０２に表示させる。また、パラメータ生成機能２０４ｃは、動画としてパラメータマップを生成した場合、パラメータマップをディスプレイ２０２に動画表示させる。これにより、操作者は、左右の肺への血流量を視認することが可能になる。
【００５９】
ステップＳ８では、パラメータ生成機能２０４ｃは、パラメータに基づいて、左右の肺のパラメータの比を算出する。例えば、パラメータ生成機能２０４ｃは、パラメータマップにおいて左右の肺それぞれの関心領域の設定を操作者から受け付け、受け付けた関心領域におけるパラメータの比を左右の肺のパラメータの比として算出する。
【００６０】
ステップＳ９では、パラメータ生成機能２０４ｃは、算出した左右の肺のパラメータの比をディスプレイ２０２に表示させる。これにより、操作者は、左右の肺への血流量をより正確に把握することが可能になる。
【００６１】
なお、パラメータ生成機能２０４ｃは、ステップＳ５において、特定された一方の心拍に対応する期間における、時系列のＸ線画像に対して設定された領域内の画素の統計値の経時的な変化に基づくパラメータを生成するようにしてもよい。かかる場合、パラメータ生成機能２０４ｃは、図３に示すステップＳ６及びステップＳ７の処理を省略し、ステップＳ８を実行する。すなわち、パラメータ生成機能２０４ｃは、ステップＳ８において、ステップＳ５で生成した継時的な変化に基づくパラメータの比を算出する。そして、パラメータ生成機能２０４ｃは、ステップＳ９において、算出したパラメータの比をディスプレイ２０２に表示させる。
【００６２】
上述したように、第１の実施形態では、医用画像処理装置２００は、連続する２心拍のうち一方の心拍を特定し、特定した一方の心拍に対応する期間における造影剤の濃度の経時的な変化に基づくパラメータを生成する。これにより、第１の実施形態によれば、肺血管造影において血流を正確に測定することができる。
【００６３】
（第１の実施形態の変形例）
上述した第１の実施形態では、連続する複数フレーム間での変化量が最大となる期間を含んだ心拍を一方の心拍として特定する場合について説明したが、実施形態はこれに限定されるものではない。例えば、特定機能２０４ｂは、第１の心拍及び第２の心拍について、造影剤が肺動脈から肺胞に流れ出す瞬間の時刻Ａと、時間濃度曲線の傾きがゆるやかになる時刻Ｂとを検出し、時刻Ａから時刻Ｂまでの間の時間濃度曲線の増加量を算出する。そして、第１の心拍の増加量と、第２の心拍の増加量を比較して、増加量が大きい心拍を一方の心拍として特定する。
【００６４】
また、上述した実施形態では、特定機能２０４ｂは、関心領域中の画素値の合計値或いは平均値の時間濃度曲線を導出する場合について説明したが、実施形態はこれに限定されるものではない。例えば、特定機能２０４ｂは、関心領域として左肺に設定された第１の関心領域及び右肺に設定された第２の関心領域それぞれに対応する領域中の各画素の変化を示す時間濃度曲線をそれぞれ導出し、導出した各時間濃度曲線を加算した時間濃度曲線を導出して一方の心拍を特定してもよい。
【００６５】
また、上述した実施形態では、左右の肺の時間濃度曲線を加算して、一方の心拍を特定するものとして説明したが、実施形態はこれに限定されるものではない。例えば、特定機能２０４ｂは、関心領域として左肺に設定された第１の関心領域及び右肺に設定された第２の関心領域それぞれに対応する領域中の各画素の変化を示す時間濃度曲線をそれぞれ導出して一方の心拍を特定してもよい。
【００６６】
また、上述した実施形態では、特定機能２０４ｂは、左右の肺それぞれに関心領域の設定を受け付ける場合について説明したが、実施形態はこれに限定されるものではない。例えば、特定機能２０４ｂは、左右の肺を含んだ関心領域の設定を受け付けてもよい。
【００６７】
図７は、第１の実施形態の変形例を説明するための図である。図７では、左右の肺が描出されたＸ線画像を示す。そして、操作者は、Ｘ線画像において、入力インターフェース２０１を介して、左右の肺を含んだ関心領域を設定する。これにより、特定機能２０４ｂは、図７に示すように、Ｘ線画像において左右の肺を含んだ関心領域５３の設定を操作者から受け付ける。
【００６８】
そして、かかる場合、特定機能２０４ｂは、左右の肺を含んだ関心領域中の画素値の合計値或いは平均値の時間濃度曲線を導出して一方の心拍を特定する。或いは、特定機能２０４ｂは、左右の肺を含んだ関心領域中の各画素の変化を示す時間濃度曲線をそれぞれ導出し、導出した各時間濃度曲線を加算した時間濃度曲線を導出して一方の心拍を特定する。
【００６９】
（第２の実施形態）
上述した第１の実施形態では、特定機能２０４ｂは、時間濃度曲線の変化量に基づいて、第１の心拍のデータ及び第２の心拍のデータそれぞれの特性を調べ、第１の心拍のデータもしくは第２の心拍のデータのいずれか一方を特定するものとして説明した。
【００７０】
ところで、医用画像処理システム１では、心電計４０により心電波形を取得可能である。ここで、造影剤が肺動脈から肺胞に流れ込むタイミングは、心拍運動に同期する。このため、例えば、第１の心拍の際に、造影剤が肺動脈に充分量蓄積されてから肺胞に流れ込む場合は、第１の心拍のデータをパラメトリックイメージング技術に利用すればよい。また、一方で、第１の心拍の際に、造影剤が肺動脈に充分量蓄積されずに肺胞に流れ込む場合は、第２の心拍のデータをパラメトリックイメージング技術に利用すればよい。
【００７１】
このようなことから、特定機能２０４ｂは、心電波形を利用して、造影剤が肺動脈から肺胞に流れるタイミングを判定し、パラメトリックイメージング技術に利用する心拍を特定してもよいものである。そこで、第２の実施形態では、心電波形に基づいて一方の心拍を特定する場合について説明する。なお、第２の実施形態に係る医用画像処理システム１の構成は、医用画像処理装置２００の処理回路が実行する特定機能２０４ｂの一部が異なる点を除いて、図１に示す第１の実施形態に係る医用画像処理システム１の構成と同様である。このため、第２の実施形態では、特定機能２０４ｂについてのみ説明する。
【００７２】
第２の実施形態に係る特定機能２０４ｂは、心電波形に基づいて一方の心拍を特定する。ここで、心房へ血液が流入し、心室から血液が駆出されるタイミングで、造影剤が肺動脈から肺胞へと流れ出す。このタイミングは、心電波形においてＲ波の頂点付近からＴ波の終わり付近までとなる。そこで、特定機能２０４ｂは、例えば、第１の心拍におけるＲ波の頂点付近のＸ線画像を選択して、肺動脈に蓄積された造影剤量が、パラメトリックイメージングに充分な量であるか否かを判定する。
【００７３】
より具体的には、特定機能２０４ｂは、肺動脈に蓄積している造影剤に対応する画素数を算出し、算出した画素数が所定の閾値以上である場合に、パラメトリックイメージングに充分な量の造影剤が肺動脈に蓄積されていると判定する。一方、特定機能２０４ｂは、算出した画素数が所定の閾値未満である場合に、パラメトリックイメージングに充分な量の造影剤が肺動脈に蓄積していないと判定する。
【００７４】
図８及び図９は、第２の実施形態を説明するための図である。図８及び図９では、心拍のタイミングであるＲ波の頂点付近のＸ線画像を示す。また、図８及び図９では、左右の肺に加えて、肺動脈に留置されたカテーテルＣを図示している。図８では、第１の心拍の際に充分量の造影剤が肺動脈に蓄積されてから肺胞に流れ込む場合を示す。例えば、図８に示すように、肺動脈には黒色で示す造影剤が蓄積している。ここで、例えば、肺動脈に蓄積している造影剤に対応する画素数が所定の閾値以上であるとする。かかる場合、特定機能２０４ｂは、パラメトリックイメージングに充分な量の造影剤が肺動脈に蓄積されていると判定する。
【００７５】
図９では、第１の心拍の際に充分量の造影剤が肺動脈に蓄積されずに肺胞に流れ込む場合を示す。例えば、図９に示すように、肺動脈にはドット状で示す造影剤が蓄積している。ここで、例えば、肺動脈に蓄積している造影剤に対応する画素数が所定の閾値未満であるとする。かかる場合、特定機能２０４ｂは、パラメトリックイメージングに充分な量の造影剤が肺動脈に蓄積されていないと判定する。
【００７６】
上述したように、第２の実施形態では、医用画像処理装置２００は、心電波形を利用して、造影剤が肺動脈から肺胞に流れ込むタイミングで充分量の造影剤が蓄積されていたか否かを推測して、一方の心拍を特定する。これにより、第２の実施形態によれば、肺血管造影において血流を正確に測定することができる。
【００７７】
（第３の実施形態）
ところで、カテーテルの向きが偏っていた場合、カテーテルの向きに応じた造影剤の流れが観測されてしまう。この場合も、第１の心拍による流れを測定するのではなく、第２の心拍による流れを測定した方が正確な測定となる。
【００７８】
このようなことから、特定機能２０４ｂは、カテーテルの向きを考慮して、パラメトリックイメージング技術に利用する心拍を特定してもよいものである。そこで、第３の実施形態では、カテーテルの向きに応じて、一方の心拍を特定する場合について説明する。なお、第３の実施形態に係る医用画像処理システム１の構成は、医用画像処理装置２００の処理回路が実行する特定機能２０４ｂの一部が異なる点を除いて、図１に示す第１の実施形態に係る医用画像処理システム１の構成と同様である。このため、第３の実施形態では、特定機能２０４ｂについてのみ説明する。
【００７９】
第３の実施形態に係る特定機能２０４ｂは、カテーテルの向きに応じて、一方の心拍を特定する。図１０は、第３の実施形態を説明するための図である。図１０では、造影剤を注入時のＸ線画像を示す。また、図１０には、カテーテルＣが描出されている。例えば、特定機能２０４ｂは、事前にカテーテル先端部の形状を記憶したテーブルを参照して、Ｘ線画像からカテーテル先端部を抽出する。そして、特定機能２０４ｂは、抽出したカテーテル先端部の位置を特定し、カテーテルの向きを判定する。例えば、特定機能２０４ｂは、図１０に示すように、カテーテルＣの先端部が右肺に向いていると判定する。かかる場合、特定機能２０４ｂは、一方の心拍として第２の心拍を特定する。
【００８０】
上述したように、第３の実施形態によれば、医用画像処理装置２００は、カテーテルの向きに応じて、一方の心拍を特定する。これにより、第３の実施形態によれば、肺血管造影において血流を正確に測定することができる。
【００８１】
（第３の実施形態の変形例）
なお、造影剤投与した瞬間に、肺動脈から肺胞への流れ込みがない場合は、投与した造影剤は肺動脈内にいったん充満し、次の心拍のタイミングで一斉に肺胞に流れ込む。このため、カテーテルの向きの影響を受けなくなる。このようなことから、特定機能２０４ｂは、カテーテルの向きが偏っていた場合に、第２の実施形態と同様にして心電波形を考慮して造影剤が肺動脈から肺胞に流れ込むタイミングを判定し、一方の心拍を特定してもよい。
【００８２】
（第４の実施形態）
上述した実施形態では、医用画像処理装置２００によって第１の心拍と第２の心拍の特性を調べて、一方の心拍を特定する場合について説明した。第４の実施形態では、操作者が目視によって一方の心拍を選択し、操作者によって選択された心拍を受け付けること、一方の心拍を特定する場合について説明する。なお、第４の実施形態に係る医用画像処理システム１の構成は、医用画像処理装置２００の処理回路が実行する特定機能２０４ｂの一部が異なる点を除いて、図１に示す第１の実施形態に係る医用画像処理システム１の構成と同様である。このため、第４の実施形態では、特定機能２０４ｂについてのみ説明する。
【００８３】
第４の実施形態に係る特定機能２０４ｂは、操作者から指定されたＸ線画像を含んだ心拍を、一方の心拍として特定する。かかる場合、例えば、特定機能２０４ｂは、取得機能２０４ａによって取得された連続する２心拍分の時系列のＸ線画像をディスプレイ２０２に表示させる。そして、操作者は、入力インターフェース２０１を介して、第１の心拍のＸ線画像又は第２の心拍のＸ線画像を選択する。ここで、操作者は、複数枚の第１の心拍のＸ線画像又は複数枚の第２の心拍のＸ線画像の一方を選択しても良い。
【００８４】
より具体的には、操作者は、第１の心拍のＸ線画像又は第２の心拍のＸ線画像において、肺胞の第１分岐に造影剤が入ったＸ線画像から、肺胞の外側が染まったＸ線画像までを選択する。そして、特定機能２０４ｂは、操作者から指定されたＸ線画像を含んだ心拍を、一方の心拍として特定する。なお、操作者は、１枚の第１の心拍のＸ線画像又は１枚の第２の心拍のＸ線画像の一方を選択しても良い。かかる場合も、特定機能２０４ｂは、選択された１枚のＸ線画像を含んだ心拍を特定する。
【００８５】
上述したように、第４の実施形態では、操作者が連続する２心拍分の時系列のＸ線画像を目視して一方の心拍を特定する。これにより、第４の実施形態によれば、肺血管造影において血流を正確に測定することができる。
【００８６】
（その他の実施形態）
実施形態は、上述した実施形態に限られるものではない。
【００８７】
上述した実施形態では、肺動脈に造影剤を注入する場合について説明したが、実施形態はこれに限定されるものではない。例えば、肺血管造影において下肢部から造影剤を注入してもよい。かかる場合、取得機能２０４ａは、造影剤注入後に時系列に収集されたＸ線画像から、造影剤が撮影対象部位に到達後の連続する２心拍分の時系列のＸ線画像を取得する。かかる場合、取得機能２０４ａは、第１の実施形態にて説明した、第１の方法から第３の方法を用いて、連続する２心拍分の時系列のＸ線画像を取得する。なお、第２の方法を用いる場合、取得機能２０４ａは、造影剤注入後に造影剤注入部位から撮影対象部位まで造影剤が到達する時間を事前に記憶しておき、この時間を更に用いて、連続する２心拍分の時系列のＸ線画像を取得する。
【００８８】
また、上述した実施形態では、医用画像処理装置２００は、特定した心拍に対応する期間における造影剤の濃度の経時的な変化に基づくパラメータを、各位置に割り当てたパラメータを生成するものとして説明したが、実施形態はこれに限定されるものではない。例えば、医用画像処理装置２００は、特定した期間のＸ線画像だけを用いて、パラメータを生成してもよい。
【００８９】
また、上述した実施形態では、肺血管造影において血流を測定する場合について説明したが、実施形態はこれに限定されるものではない。例えば、心拍の影響を受ける臓器であれば、肺以外の臓器において血流を測定する場合にも、上述した実施形態が適用可能である。かかる場合、特定機能２０４ｂは、心拍の影響を受ける臓器に設定された関心領域中の画素値の合計値或いは平均値の時間濃度曲線を導出して一方の心拍を特定する。或いは、特定機能２０４ｂは、心拍の影響を受ける臓器に設定された関心領域中の各画素の変化を示す時間濃度曲線をそれぞれ導出し、導出した各時間濃度曲線を加算した時間濃度曲線を導出して一方の心拍を特定する。
【００９０】
また、上述した実施形態では、医用画像処理装置２００は、Ｘ線診断装置１００により生成されたＸ線画像を取得して、パラメータ生成処理を実行する場合につて説明したが、実施形態はこれに限定されるものではない。例えば、医用画像処理装置２００は、Ｘ線ＣＴ装置によって生成されたＣＴ画像を取得して、第１の実施形態から第４の実施形態にて説明したパラメータ生成処理を実行してもよい。
【００９１】
また、上述した実施形態では、医用画像処理装置２００において、パラメータ生成処理を実行する場合につて説明したが、実施形態はこれに限定されるものではない。例えば、Ｘ線診断装置１００において、パラメータ生成処理を実行してもよい。かかる場合、Ｘ線診断装置１００の処理回路２１は、上述した第１の実施形態から第４の実施形態に係る医用画像処理装置２００の処理回路２０４と同様のパラメータ生成処理を実行する。すなわち、Ｘ線診断装置１００において、パラメータ生成処理を実行する場合、処理回路２１は、取得機能２１ａと、特定機能２１ｂと、パラメータ生成機能２１ｃとを実行する。ここで、例えば、取得機能２１ａは、上述した取得機能２０４ａと同様の処理を実行する。特定機能２１ｂは、上述した特定機能２０４ｂと同様の処理を実行する。パラメータ生成機能２１ｃは、上述したパラメータ生成機能２０４ｃと同様の処理を実行する。
【００９２】
また、上述した実施形態では、医療画像処理システム１には、インジェクター３０及び心電計４０が含まれる場合について説明したが、実施形態はこれに限定されるものではない。例えば、医療画像処理システム１には、インジェクター３０及び心電計４０が含まれなくてもよい。かかる場合、造影剤は、術者によって被検体Ｐに注入される。
【００９３】
上記説明において用いた「プロセッサ」という文言は、例えば、ＣＰＵ（Central  Processing  Unit）、ＧＰＵ(Graphics  Processing  Unit)、或いは、特定用途向け集積回路（Application  Specific  Integrated  Circuit：ＡＳＩＣ）、プログラマブル論理デバイス（例えば、単純プログラマブル論理デバイス（Simple  Programmable  Logic  Device：ＳＰＬＤ）、複合プログラマブル論理デバイス（Complex  Programmable  Logic  Device：ＣＰＬＤ）、及びフィールドプログラマブルゲートアレイ（Field  Programmable  Gate  Array：ＦＰＧＡ））等の回路を意味する。プロセッサは記憶回路に保存されたプログラムを読み出し実行することで機能を実現する。なお、記憶回路にプログラムを保存する代わりに、プロセッサの回路内にプログラムを直接組み込むよう構成しても構わない。この場合、プロセッサは回路内に組み込まれたプログラムを読み出し実行することで機能を実現する。なお、本実施形態の各プロセッサは、プロセッサごとに単一の回路として構成される場合に限らず、複数の独立した回路を組み合わせて１つのプロセッサとして構成し、その機能を実現するようにしてもよい。さらに、図１における複数の構成要素を１つのプロセッサへ統合してその機能を実現するようにしてもよい。
【００９４】
上記の実施形態の説明において、図示した各装置の各構成要素は機能概念的なものであり、必ずしも物理的に図示の如く構成されていることを要しない。すなわち、各装置の分散・統合の具体的形態は図示のものに限られず、その全部又は一部を、各種の負荷や使用状況等に応じて、任意の単位で機能的又は物理的に分散・統合して構成することができる。さらに、各装置にて行なわれる各処理機能は、その全部または任意の一部が、ＣＰＵおよび当該ＣＰＵにて解析実行されるプログラムにて実現され、或いは、ワイヤードロジックによるハードウェアとして実現され得る。
【００９５】
また、上記の実施形態で説明した制御方法は、予め用意された制御プログラムをパーソナルコンピュータやワークステーション等のコンピュータで実行することによって実現することができる。この制御プログラムは、インターネット等のネットワークを介して配布することができる。また、この制御プログラムは、ハードディスク、フレキシブルディスク（ＦＤ）、ＣＤ−ＲＯＭ、ＭＯ、ＤＶＤ等のコンピュータで読み取り可能な記録媒体に記録され、コンピュータによって記録媒体から読み出されることによって実行することもできる。
【００９６】
以上説明した少なくとも一つの実施形態によれば、肺血管造影において血流を正確に測定することができる。
【００９７】
本発明のいくつかの実施形態を説明したが、これらの実施形態は、例として提示したものであり、発明の範囲を限定することは意図していない。これら実施形態は、その他の様々な形態で実施されることが可能であり、発明の要旨を逸脱しない範囲で、種々の省略、置き換え、変更を行うことができる。これら実施形態やその変形は、発明の範囲や要旨に含まれると同様に、特許請求の範囲に記載された発明とその均等の範囲に含まれるものである。
【符号の説明】
【００９８】
１  医用画像処理システム
２００  医用画像処理装置
２０４  処理回路
２０４ａ  取得機能
２０４ｂ  特定機能
２０４ｃ  パラメータ生成機能
【図１】
0006987550.tif 000002
【図２】
0006987550.tif 000003
【図３】
0006987550.tif 000004
【図４】
0006987550.tif 000005
【図５】
0006987550.tif 000006
【図６】
0006987550.tif 000007
【図７】
0006987550.tif 000008
【図８】
0006987550.tif 000009
【図９】
0006987550.tif 000010
【図１０】
0006987550.tif 000011