(57)【要約】【課題】適切なデータを簡便に選択することができるＸ線ＣＴ装置、医用情報処理装置及び医用情報処理プログラムを提供すること。【解決手段】実施形態に係るＸ線ＣＴ装置は、収集部と、選択部と、解析部とを備える。収集部は、心臓の冠動脈の少なくとも一部を含む、複数時相に係る画像データを収集する。選択部は、複数の前記画像データに関する画像指標を取得し、各画像データにおける画像指標に基づいて、複数の画像データのうち、所定の時間間隔よりも大きい時間間隔を有する画像データの組み合わせの中から画像データの組を抽出する。解析部は、抽出された画像データの組に基づいて、冠動脈に関する流体解析を行い、冠動脈に関する流体パラメータを求める。【選択図】図２
2017217459.tif 000002
【特許請求の範囲】
【請求項１】
心臓の冠動脈の少なくとも一部を含む、複数時相に係る画像データを収集する収集部と、
複数の前記画像データに関する画像指標を取得し、各画像データにおける画像指標に基づいて、複数の前記画像データのうち、所定の時間間隔よりも大きい時間間隔を有する画像データの組み合わせの中から画像データの組を抽出する選択部と、
前記抽出された画像データの組に基づいて、前記冠動脈に関する流体解析を行い、前記冠動脈に関する流体パラメータを求める解析部と、
を備える、Ｘ線ＣＴ装置。
【請求項２】
前記選択部は、前記画像データの組を、１心拍の心拍時相中におけるWave  Free  Periodとなる時相中から抽出する、請求項１記載のＸ線ＣＴ装置。
【請求項３】
前記選択部は、前記画像データの組を、１心拍の心拍時相中の、７０%以降の心時相から抽出する、請求項１記載のＸ線ＣＴ装置。
【請求項４】
前記選択部は、前記画像データの組み合わせの中から、前記各画像データにおける画質指標が所定の閾値以上である画像データの組を抽出する、請求項１〜３のいずれか１つに記載のＸ線ＣＴ装置。
【請求項５】
前記選択部は、前記各画像データに含まれる所定の部位の抽出結果に基づいて、前記画像データの組み合わせの中から前記画像データの組を抽出する、請求項１〜３のいずれか１つに記載のＸ線ＣＴ装置。
【請求項６】
前記所定の部位は、前記冠動脈の血管壁又は前記冠動脈の芯線である、請求項５記載のＸ線ＣＴ装置。
【請求項７】
前記選択部は、前記画像指標に基づいて、Wave  Free  Periodに含まれる時相の各画像データの前記心臓の拍動に伴う動きを評価し、前記拍動に伴う動きが相対的に小さい時相の画像データの組を抽出する、請求項１記載のＸ線ＣＴ装置。
【請求項８】
前記選択部は、前記画像データに含まれる冠動脈領域における鮮明度に基づいて、前記拍動に伴う動きを評価する、請求項７記載のＸ線ＣＴ装置。
【請求項９】
前記選択部は、前記Wave  Free  Periodを複数の部分範囲に分割し、各部分範囲について、前記部分範囲に含まれる時相の中から前記拍動に伴う動きが相対的に小さい時相に対応する画像データの組を抽出する、請求項７又は８記載のＸ線ＣＴ装置。
【請求項１０】
前記選択部は、前記複数の部分範囲の幅がそれぞれ近似するように、前記Wave  Free  Periodを前記部分範囲に分割する、請求項９記載のＸ線ＣＴ装置。
【請求項１１】
前記選択部は、前記画像データの組として、少なくとも前記Wave  Free  Periodにおける上端付近の時相及び下端付近の時相に対応する画像データを抽出する、請求項７又は８記載のＸ線ＣＴ装置。
【請求項１２】
前記選択部は、前記Wave  Free  Periodに含まれる時相において、前記時相間の時間間隔が所定の閾値以上となるように、前記複数の時相に対応する画像データを抽出する、請求項７又は８記載のＸ線ＣＴ装置。
【請求項１３】
前記選択部は、前記画像データの組において隣接する時相間の各時間間隔が大きく、かつ、相互に近似するように、前記画像データの組を抽出する、請求項１〜１２のいずれか１つに記載のＸ線ＣＴ装置。
【請求項１４】
前記選択部は、前記画像データの転送時に前記画像データの組を抽出する、請求項１〜１３のいずれか１つに記載のＸ線ＣＴ装置。
【請求項１５】
前記選択部は、前記流体解析のアプリケーションの起動時に前記画像データの組を抽出する、請求項１〜１３のいずれか１つに記載のＸ線ＣＴ装置。
【請求項１６】
前記選択部が抽出した前記画像データの組が所定の条件に適合しない場合に、操作者に対して警告を提示する提示部をさらに備える、請求項１〜１５のいずれか１つに記載のＸ線ＣＴ装置。
【請求項１７】
心臓の冠動脈の少なくとも一部を含む、複数時相に係る画像データに関する画像指標を取得し、各画像データにおける画像指標に基づいて、複数の前記画像データのうち、所定の時間間隔よりも大きい時間間隔を有する画像データの組み合わせの中から画像データの組を抽出する選択部と、
前記抽出された画像データの組に基づいて、前記冠動脈に関する流体解析を行い、前記冠動脈に関する流体パラメータを求める解析部と、
を備える、
医用情報処理装置。
【請求項１８】
心臓の冠動脈の少なくとも一部を含む、複数時相に係る画像データに関する画像指標を取得し、各画像データにおける画像指標に基づいて、複数の前記画像データのうち、所定の時間間隔よりも大きい時間間隔を有する画像データの組み合わせの中から画像データの組を抽出する選択ステップと、
前記抽出された画像データの組に基づいて、前記冠動脈に関する流体解析を行い、前記冠動脈に関する流体パラメータを求める解析ステップと、
をコンピュータに実行させる、医用情報処理プログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明の実施形態は、Ｘ線ＣＴ装置、医用情報処理装置及び医用情報処理プログラムに関する。
【背景技術】
【０００２】
従来、臓器の虚血性疾患の原因には、大別して血行障害と臓器そのものの機能障害とがあることが知られている。例えば、冠動脈の血行障害の一例である狭窄は、虚血性心疾患に至る重大な病変であるが、そのような虚血性心疾患では、薬物治療を行うべきか、ステント治療を行うべきか等を判断する必要がある。近年では、冠動脈の血行性虚血評価を行う診断として、カテーテルによる冠動脈造影検査（ＣＡＧ：Coronary  Angiography）においてプレッシャーワイヤを用いて心筋血流予備量比（ＦＦＲ：Fractional  Flow  Reserve）を計測する手法が推奨されつつある。
【０００３】
これに対し、例えば、Ｘ線ＣＴ（Computed  Tomography）装置やＭＲＩ（Magnetic  Resonance  Imaging）装置、超音波診断装置等の医用画像診断装置によって収集された心臓の医用画像を用いて冠動脈の血行性虚血評価を非侵襲的に行う手法も知られている。近年、このような種々の手法により血行性虚血評価が行われ、評価に応じた治療が行われる。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特開２０１４−１０８２０８号公報
【特許文献２】特開２０１４−０８７６３５号公報
【特許文献３】特開２０１５−１１０１５５号公報
【発明の概要】
【発明が解決しようとする課題】
【０００５】
本発明が解決しようとする課題は、適切なデータを簡便に選択することができるＸ線ＣＴ装置、医用情報処理装置及び医用情報処理プログラムを提供することである。
【課題を解決するための手段】
【０００６】
実施形態に係るＸ線ＣＴ装置は、収集部と、選択部と、解析部とを備える。収集部は、心臓の冠動脈の少なくとも一部を含む、複数時相に係る画像データを収集する。選択部は、複数の前記画像データに関する画像指標を取得し、各画像データにおける画像指標に基づいて、複数の前記画像データのうち、所定の時間間隔よりも大きい時間間隔を有する画像データの組み合わせの中から画像データの組を抽出する。解析部は、前記抽出された画像データの組に基づいて、前記冠動脈に関する流体解析を行い、前記冠動脈に関する流体パラメータを求める。
【図面の簡単な説明】
【０００７】
【図１】図１は、第１の実施形態に係る医用情報処理システムの構成の一例を示す図である。
【図２】図２は、第１の実施形態に係るＸ線ＣＴ装置の構成の一例を示す図である。
【図３】図３は、第１の実施形態に係る解析機能による処理の一例を説明するための図である。
【図４】図４は、第１の実施形態に係る選択機能による時相の選択を説明するための図である。
【図５】図５は、第１の実施形態に係る各心位相の画質の変化の一例を示す図である。
【図６】図６は、第１の実施形態に係る選択機能による選択処理の一例を説明するための図である。
【図７】図７は、第１の実施形態に係る選択機能による最適心位相の検索処理の一例を説明するための図である。
【図８】図８は、第１の実施形態に係る選択機能による閾値処理の一例を説明するための図である。
【図９】図９は、第１の実施形態に係る選択機能による選択処理の一例を説明するための図である。
【図１０】図１０は、第１の実施形態に係るＸ線ＣＴ装置による処理手順を示すフローチャートである。
【図１１】図１１は、第２の実施形態に係る選択機能による選択処理の一例を説明するための図である。
【図１２】図１２は、第２の実施形態に係るＸ線ＣＴ装置による処理手順を示すフローチャートである。
【図１３】図１３は、第３の実施形態に係る選択機能による時相の選択を説明するための図である。
【図１４】図１４は、第３の実施形態に係るＸ線ＣＴ装置による処理手順を示すフローチャートである。
【図１５】図１５は、第３の実施形態に係る医用情報処理装置の構成の一例を示す図である。
【発明を実施するための形態】
【０００８】
以下に添付図面を参照して、本願に係るＸ線ＣＴ装置、医用情報処理装置及び医用情報処理プログラムの実施形態を詳細に説明する。なお、本願に係るＸ線ＣＴ装置、医用情報処理装置及び医用情報処理プログラムは、以下に示す実施形態によって限定されるものではない。
【０００９】
（第１の実施形態）
まず、第１の実施形態について説明する。第１の実施形態では、本願が開示する技術をＸ線ＣＴ装置に適用した場合の例を説明する。なお、以下、Ｘ線ＣＴ装置を含む医用情報処理システムを例に挙げて説明する。また、以下では、一例として、心臓の血管を解析対象とした場合の例を説明する。
【００１０】
図１は、第１の実施形態に係る医用情報処理システムの構成の一例を示す図である。図１に示すように、第１の実施形態に係る医用情報処理システムは、Ｘ線ＣＴ（Computed  Tomography）装置１００と、画像保管装置２００と、医用情報処理装置３００とを備える。
【００１１】
例えば、第１の実施形態に係るＸ線ＣＴ装置１００は、図１に示すように、ネットワーク４００を介して、画像保管装置２００と、医用情報処理装置３００に接続される。なお、医用情報処理システムは、ネットワーク４００を介して、ＭＲＩ装置や超音波診断装置、ＰＥＴ（Positron  Emission  Tomography）装置等の他の医用画像診断装置にさらに接続されてもよい。
【００１２】
画像保管装置２００は、各種の医用画像診断装置によって収集された画像データを保管する。例えば、画像保管装置２００は、サーバ装置等のコンピュータ機器によって実現される。本実施形態では、画像保管装置２００は、ネットワーク４００を介してＸ線ＣＴ装置１００からＣＴ画像データ（ボリュームデータ）を取得し、取得したＣＴ画像データを装置内又は装置外に設けられた記憶回路に記憶させる。
【００１３】
医用情報処理装置３００は、ネットワーク４００を介して各種の医用画像診断装置から画像データを取得し、取得した画像データを処理する。例えば、医用情報処理装置３００は、ワークステーション等のコンピュータ機器によって実現される。本実施形態では、医用情報処理装置３００は、ネットワーク４００を介してＸ線ＣＴ装置１００又は画像保管装置２００からＣＴ画像データを取得し、取得したＣＴ画像データに対して各種画像処理を行う。そして、医用情報処理装置３００は、画像処理を行う前又は行った後のＣＴ画像データをディスプレイ等に表示する。
【００１４】
Ｘ線ＣＴ装置１００は、被検体のＣＴ画像データを収集する。具体的には、Ｘ線ＣＴ装置１００は、被検体を略中心にＸ線管及びＸ線検出器を旋回移動させ、被検体を透過したＸ線を検出して投影データを収集する。そして、Ｘ線ＣＴ装置１００は、収集された投影データに基づいて、時系列の３次元ＣＴ画像データを生成する。図２は、第１の実施形態に係るＸ線ＣＴ装置１００の構成の一例を示す図である。図２に示すように、第１の実施形態に係るＸ線ＣＴ装置１００は、架台１０と、寝台装置２０と、コンソール３０とを有する。
【００１５】
架台１０は、被検体Ｐ（患者）にＸ線を照射し、被検体Ｐを透過したＸ線を検出して、コンソール３０に出力する装置であり、Ｘ線照射制御回路１１と、Ｘ線発生装置１２と、検出器１３と、データ収集回路（ＤＡＳ：Data  Acquisition  System）１４と、回転フレーム１５と、架台駆動回路１６とを有する。
【００１６】
回転フレーム１５は、Ｘ線発生装置１２と検出器１３とを被検体Ｐを挟んで対向するように支持し、後述する架台駆動回路１６によって被検体Ｐを中心とした円軌道にて高速に回転する円環状のフレームである。
【００１７】
Ｘ線照射制御回路１１は、高電圧発生部として、Ｘ線管１２ａに高電圧を供給する装置であり、Ｘ線管１２ａは、Ｘ線照射制御回路１１から供給される高電圧を用いてＸ線を発生する。Ｘ線照射制御回路１１は、後述するスキャン制御回路３３の制御により、Ｘ線管１２ａに供給する管電圧や管電流を調整することで、被検体Ｐに対して照射されるＸ線量を調整する。
【００１８】
また、Ｘ線照射制御回路１１は、ウェッジ１２ｂの切り替えを行う。また、Ｘ線照射制御回路１１は、コリメータ１２ｃの開口度を調整することにより、Ｘ線の照射範囲（ファン角やコーン角）を調整する。なお、本実施形態は、複数種類のウェッジを、操作者が手動で切り替える場合であっても良い。
【００１９】
Ｘ線発生装置１２は、Ｘ線を発生し、発生したＸ線を被検体Ｐへ照射する装置であり、Ｘ線管１２ａと、ウェッジ１２ｂと、コリメータ１２ｃとを有する。
【００２０】
Ｘ線管１２ａは、図示しない高電圧発生部により供給される高電圧により被検体ＰにＸ線ビームを照射する真空管であり、回転フレーム１５の回転にともなって、Ｘ線ビームを被検体Ｐに対して照射する。Ｘ線管１２ａは、ファン角及びコーン角を持って広がるＸ線ビームを発生する。例えば、Ｘ線照射制御回路１１の制御により、Ｘ線管１２ａは、フル再構成用に被検体Ｐの全周囲でＸ線を連続曝射したり、ハーフ再構成用にハーフ再構成可能な曝射範囲（１８０度＋ファン角）でＸ線を連続曝射したりすることが可能である。また、Ｘ線照射制御回路１１の制御により、Ｘ線管１２ａは、予め設定された位置（管球位置）でＸ線（パルスＸ線）を間欠曝射したりすることが可能である。また、Ｘ線照射制御回路１１は、Ｘ線管１２ａから曝射されるＸ線の強度を変調させることも可能である。例えば、Ｘ線照射制御回路１１は、特定の管球位置では、Ｘ線管１２ａから曝射されるＸ線の強度を強くし、特定の管球位置以外の範囲では、Ｘ線管１２ａから曝射されるＸ線の強度を弱くする。
【００２１】
ウェッジ１２ｂは、Ｘ線管１２ａから曝射されたＸ線のＸ線量を調節するためのＸ線フィルタである。具体的には、ウェッジ１２ｂは、Ｘ線管１２ａから被検体Ｐへ照射されるＸ線が、予め定められた分布になるように、Ｘ線管１２ａから曝射されたＸ線を透過して減衰するフィルタである。例えば、ウェッジ１２ｂは、所定のターゲット角度や所定の厚みとなるようにアルミニウムを加工したフィルタである。なお、ウェッジは、ウェッジフィルタ（wedge  filter）や、ボウタイフィルタ（bow-tie  filter）とも呼ばれる。
【００２２】
コリメータ１２ｃは、後述するＸ線照射制御回路１１の制御により、ウェッジ１２ｂによってＸ線量が調節されたＸ線の照射範囲を絞り込むためのスリットである。
【００２３】
架台駆動回路１６は、回転フレーム１５を回転駆動させることによって、被検体Ｐを中心とした円軌道上でＸ線発生装置１２と検出器１３とを旋回させる。
【００２４】
検出器１３は、被検体Ｐを透過したＸ線を検出する２次元アレイ型検出器（面検出器）であり、複数チャンネル分のＸ線検出素子を配してなる検出素子列が被検体Ｐの体軸方向（図２に示すＺ軸方向）に沿って複数列配列されている。具体的には、第１の実施形態における検出器１３は、被検体Ｐの体軸方向に沿って３２０列など多列に配列されたＸ線検出素子を有し、例えば、被検体Ｐの肺や心臓を含む範囲など、広範囲に被検体Ｐを透過したＸ線を検出することが可能である。
【００２５】
データ収集回路１４は、ＤＡＳであり、検出器１３が検出したＸ線の検出データから、投影データを収集する。例えば、データ収集回路１４は、検出器１３により検出されたＸ線強度分布データに対して、増幅処理やＡ／Ｄ変換処理、チャンネル間の感度補正処理等を行なって投影データを生成し、生成した投影データを後述するコンソール３０に送信する。例えば、回転フレーム１５の回転中に、Ｘ線管１２ａからＸ線が連続曝射されている場合、データ収集回路１４は、全周囲分（３６０度分）の投影データ群を収集する。また、データ収集回路１４は、収集した各投影データに管球位置を対応付けて、後述するコンソール３０に送信する。管球位置は、投影データの投影方向を示す情報となる。なお、チャンネル間の感度補正処理は、後述する前処理回路３４が行なっても良い。
【００２６】
寝台装置２０は、被検体Ｐを載せる装置であり、図２に示すように、寝台駆動装置２１と、天板２２とを有する。寝台駆動装置２１は、天板２２をＺ軸方向へ移動して、被検体Ｐを回転フレーム１５内に移動させる。天板２２は、被検体Ｐが載置される板である。
【００２７】
なお、架台１０は、例えば、天板２２を移動させながら回転フレーム１５を回転させて被検体Ｐをらせん状にスキャンするヘリカルスキャンを実行する。または、架台１０は、天板２２を移動させた後に被検体Ｐの位置を固定したままで回転フレーム１５を回転させて被検体Ｐを円軌道にてスキャンするコンベンショナルスキャンを実行する。または、架台１０は、天板２２の位置を一定間隔で移動させてコンベンショナルスキャンを複数のスキャンエリアで行うステップアンドシュート方式を実行する。
【００２８】
コンソール３０は、操作者によるＸ線ＣＴ装置１００の操作を受け付けるとともに、架台１０によって収集された投影データを用いてＣＴ画像データ（ボリュームデータ）を再構成する装置である。コンソール３０は、図２に示すように、入力回路３１と、ディスプレイ３２と、スキャン制御回路３３と、前処理回路３４と、記憶回路３５と、画像再構成回路３６と、処理回路３７とを有する。
【００２９】
入力回路３１は、Ｘ線ＣＴ装置１００の操作者が各種指示や各種設定の入力に用いるマウスやキーボード、トラックボール、スイッチ、ボタン、ジョイスティック等を有し、操作者から受け付けた指示や設定の情報を、処理回路３７に転送する。例えば、入力回路３１は、操作者から、ＣＴ画像データの撮影条件や、ＣＴ画像データを再構成する際の再構成条件、ＣＴ画像データに対する画像処理条件等を受け付ける。また、入力回路３１は、被検体Ｐに対する検査を選択するための操作を受け付ける。また、入力回路３１は、画像上の部位を指定するための指定操作を受け付ける。
【００３０】
ディスプレイ３２は、操作者によって参照されるモニタであり、処理回路３７による制御のもと、ＣＴ画像データから生成された画像データを操作者に表示したり、入力回路３１を介して操作者から各種指示や各種設定等を受け付けるためのＧＵＩ（Graphical  User  Interface）を表示したりする。また、ディスプレイ３２は、スキャン計画の計画画面や、スキャン中の画面などを表示する。
【００３１】
スキャン制御回路３３は、処理回路３７による制御のもと、Ｘ線照射制御回路１１、架台駆動回路１６、データ収集回路１４及び寝台駆動装置２１の動作を制御することで、架台１０における投影データの収集処理を制御する。具体的には、スキャン制御回路３３は、位置決め画像（スキャノ画像）を収集する撮影及び診断に用いる画像を収集する本撮影（スキャン）における投影データの収集処理をそれぞれ制御する。
【００３２】
前処理回路３４は、データ収集回路１４によって生成された投影データに対して、対数変換処理と、オフセット補正、感度補正及びビームハードニング補正等の補正処理とを行なって、補正済みの投影データを生成する。具体的には、前処理回路３４は、データ収集回路１４によって生成された位置決め画像の投影データ及び本撮影によって収集された投影データのそれぞれについて、補正済みの投影データを生成して、記憶回路３５に格納する。
【００３３】
記憶回路３５は、前処理回路３４により生成された投影データを記憶する。具体的には、記憶回路３５は、前処理回路３４によって生成された、位置決め画像の投影データ及び本撮影によって収集される診断用の投影データを記憶する。また、記憶回路３５は、後述する画像再構成回路３６によって再構成されたＣＴ画像データなどを記憶する。また、記憶回路３５は、後述する処理回路３７による処理結果を適宜記憶する。
【００３４】
画像再構成回路３６は、記憶回路３５が記憶する投影データを用いてＣＴ画像データを再構成する。具体的には、画像再構成回路３６は、位置決め画像の投影データ及び診断に用いられる画像の投影データから、ＣＴ画像データをそれぞれ再構成する。ここで、再構成方法としては、種々の方法があり、例えば、逆投影処理が挙げられる。また、逆投影処理としては、例えば、ＦＢＰ（Filtered  Back  Projection）法による逆投影処理が挙げられる。或いは、画像再構成回路３６は、逐次近似法を用いてＣＴ画像データを再構成することもできる。
【００３５】
また、画像再構成回路３６は、ＣＴ画像データに対して各種画像処理を行うことで、画像データを生成する。そして、画像再構成回路３６は、再構成したＣＴ画像データや、各種画像処理により生成した画像データを記憶回路３５に格納する。
【００３６】
処理回路３７は、架台１０、寝台装置２０及びコンソール３０の動作を制御することによって、Ｘ線ＣＴ装置１００の全体制御を行う。具体的には、処理回路３７は、スキャン制御回路３３を制御することで、架台１０で行なわれるＣＴスキャンを制御する。また、処理回路３７は、画像再構成回路３６を制御することで、コンソール３０における画像再構成処理や画像生成処理を制御する。また、処理回路３７は、記憶回路３５が記憶する各種画像データを、ディスプレイ３２に表示するように制御する。
【００３７】
このような構成のもと、本実施形態に係るＸ線ＣＴ装置１００は、適切なデータを簡便に選択することを可能にする。具体的には、Ｘ線ＣＴ装置１００は、血管を含む医用画像（例えば、３次元ＣＴ画像データ等）を用いた流体解析によって血流に関する指標値を算出する際のデータとして適切なデータを簡便に選択することを可能にする。例えば、Ｘ線ＣＴ装置１００は、経時的に収集した複数時相のＣＴ画像データから、流体解析に用いるデータとして適切な時相のＣＴ画像データを選択する。
【００３８】
ここで、血流に関する指標値としては、例えば、心筋血流予備量比（ＦＦＲ：Fractional  Flow  Reserve）や、血管内の力学的な指標、血液の流量に関する指標などが挙げられる。ＦＦＲは、血管において心臓から近い近位部における圧力と、心臓から遠い遠位部における圧力との比であり、例えば、「ＦＦＲ＝Ｐｄ（遠位部の圧力）／Ｐａ（近位部の圧力）」で表される。例えば、血管に狭窄（治療対象部位）が生じている場合、狭窄により遠位部の圧力が低下するため、ＦＦＲの値は低下する。流体解析は、このようなＦＦＲの値などを算出することで、治療の要否などの判断に利用される。なお、血管内の力学的な指標としては、例えば、圧力、ベクトル、せん断応力などが挙げられる。また、血液の流量に関する指標としては、流量及び流速などが挙げられる。
【００３９】
以下、第１の実施形態では、冠動脈を対象とした流体解析に際して、適切なＣＴ画像データを選択する場合を例に挙げて説明する。本実施形態に係る処理回路３７は、図２に示すように、制御機能３７ａ、解析機能３７ｂ、選択機能３７ｃ、取得機能３７ｄ及び提示機能３７ｅを実行する。ここで、例えば、図２に示す処理回路３７の構成要素である制御機能３７ａ、解析機能３７ｂ、選択機能３７ｃ、取得機能３７ｄ及び提示機能３７ｅが実行する各処理機能は、コンピュータによって実行可能なプログラムの形態で記憶回路３５に記録されている。処理回路３７は、各プログラムを記憶回路３５から読み出し、実行することで各プログラムに対応する機能を実現するプロセッサである。換言すると、各プログラムを読み出した状態の処理回路３７は、図２の処理回路３７内に示された各機能を有することとなる。なお、本実施形態で説明する解析機能３７ｂは、特許請求の範囲に記載した解析部の一例である。また、選択機能３７ｃは、特許請求の範囲に記載した選択部の一例である。また、取得機能３７ｄは、取得部の一例である。また、提示機能３７ｅは、特許請求の範囲に記載した提示部の一例である。
【００４０】
制御機能３７ａは、Ｘ線ＣＴ装置１００の全体制御を実行する。例えば、制御機能３７ａは、上述したスキャンに係る制御、ＣＴ画像データの再構成及びＣＴ画像データから各種画像の生成に係る制御、及び、生成した各種画像の表示制御などを実行する。
【００４１】
解析機能３７ｂは、ＣＴ画像データに基づいて流体解析を実行する。具体的には、解析機能３７ｂは、３次元のＣＴ画像データから血管の形状を表す時系列の血管形状データを抽出する。例えば、解析機能３７ｂは、記憶回路３５から経時的に収集された複数時相のＣＴ画像データを読み出し、読み出した複数時相のＣＴ画像データに対して画像処理を行うことで、時系列の血管形状データを抽出する。
【００４２】
ここで、解析機能３７ｂは、ＣＴ画像データに含まれる血管領域に指標値を算出する対象領域を設定する。具体的には、解析機能３７ｂは、操作者による入力回路３１を介した指示又は画像処理によって、血管領域に対象領域を設定する。そして、解析機能３７ｂは、設定した対象領域の血管形状データとして、例えば、血管の芯線（芯線の座標情報）、芯線に垂直な断面での血管及び内腔の断面積、芯線に垂直な断面での円柱方向の、芯線から内壁までの距離及び芯線から外壁までの距離などをＣＴ画像データから抽出する。なお、解析機能３７ｂは、解析手法に応じて、その他種々の血管形状データを抽出することができる。
【００４３】
さらに、解析機能３７ｂは、流体解析の解析条件を設定する。具体的には、解析機能３７ｂは、解析条件として、血液の物性値、反復計算の条件、解析の初期値などを設定する。例えば、解析機能３７ｂは、血液の物性値として、血液の粘性、密度などを設定する。また、解析機能３７ｂは、反復計算の条件として、反復計算における最大反復回数、緩和係数、残差の許容値などを設定する。また、解析機能３７ｂは、解析の初期値として、流量、圧力、流体抵抗、圧力境界の初期値などを設定する。なお、解析機能３７ｂによって用いられる各種値は、システムに予め組み込んでおいてもよいし、操作者が対話的に定義してもよい。
【００４４】
また、解析機能３７ｂは、画像データ中の血管における治療対象部位を設定する。具体的には、解析機能３７ｂは、血管における治療対象部位を、手動又は自動で設定する。例えば、解析機能３７ｂは、入力回路３１を介して受け付けられた範囲を治療対象部位として設定する。かかる場合には、入力回路３１が解析条件を変更する範囲（治療対象部位）を受け付け、解析機能３７ｂが受け付けられた範囲を治療対象部位として設定する。また、解析機能３７ｂは、対象領域内の形状に基づいて、治療対象部位を自動で設定する。例えば、解析機能３７ｂは、対象領域内の形状に基づいて狭窄部分を抽出し、抽出した狭窄部分のうち一定の狭窄度以上となる狭窄部分を治療対象部位として設定する。なお、狭窄部分の抽出は、任意の手法を用いることができる。
【００４５】
そして、解析機能３７ｂは、血管を含む画像データを用いた流体解析により血管の血流に関する指標値を算出する。具体的には、解析機能３７ｂは、血管形状データと解析条件とを用いた流体解析を実行し、血管の対象領域における血流に関する指標値（流体パラメータ）を算出する。例えば、解析機能３７ｂは、血管の内腔や外壁の輪郭、血管の断面積及び芯線などの血管形状データと、血液の物性値、反復計算の条件及び解析の初期値などの設定条件に基づいて、血管の所定の位置ごとに、圧力、血液の流量、血液の流速、ベクトル及びせん断応力などの指標値を算出する。そして、解析機能３７ｂは、算出した指標値からさらにＦＦＲなどの指標値を算出する。
【００４６】
図３は、第１の実施形態に係る解析機能３７ｂによる処理の一例を説明するための図である。図３に示すように、例えば、解析機能３７ｂは、大動脈及び冠動脈を含む３次元のＣＴ画像データから、対象領域であるＬＡＤについて、芯線の座標や断面情報を含む血管形状データを抽出する。さらに、解析機能３７ｂは、抽出されたＬＡＤを対象とする解析の解析条件を設定する。そして、解析機能３７ｂは、抽出されたＬＡＤの血管形状データ及び設定された条件を用いて流体解析を行うことで、例えば、対象領域ＬＡＤの入口の境界から出口の境界まで、芯線に沿った所定の位置ごとに圧力、血液の流量、血液の流速、ベクトル及びせん断応力などの指標値を算出する。すなわち、解析機能３７ｂは、対象領域について、圧力、血液の流量、血液の流速、ベクトル及びせん断応力などの分布を算出する。そして、解析機能３７ｂは、例えば、算出した圧力分布に基づいて、対象領域における各位置のＦＦＲを算出する。
【００４７】
上述したように、解析機能３７ｂは、経時的に収集された複数時相のＣＴ画像データからそれぞれ血管形状データを抽出し、抽出した複数時相の血管形状データと解析条件とを用いた流体解析を行うことで、血流に関する指標値を算出する。ここで、このような冠動脈の流体解析を実行する際に精度の高い解析結果を得るためには、冠動脈の形状の変化（例えば、断面積の変化）をできるだけ含むような複数時相であり、かつ、各時相のＣＴ画像データができるだけシャープなものであることが望ましい。すなわち、時系列の血管形状データが、血液が流入して冠動脈の面積が最大となる時相から血液が流出して面積が最小となる時相までの一連の変化をできるだけ含みつつ、拍動による画像の動き（ぼけ）が少ない時相のＣＴ画像データを用いることが望ましい。そこで、第１の実施形態では、選択機能３７ｃが、経時的に収集されたＣＴ画像データの各時相の中から、上記したＣＴ画像データに対応する時相を選択する。
【００４８】
選択機能３７ｃは、経時的に収集された複数時相のＣＴ画像データの中から、冠動脈の形状の変化を含みつつ、拍動による画像のぼけが少ない時相を選択する。具体的には、選択機能３７ｃは、複数のＣＴ画像データに対応する各時相から、拍動に伴う動きが相対的に小さい時相であり、かつ、時相間の時間間隔が大きくなるように複数の時相を選択する。より具体的には、選択機能３７ｃは、複数のＣＴ画像データに対応する各時相において、拍動に伴う動きが相対的に小さい時相の範囲を抽出し、抽出した範囲に含まれる時相に対応する画像データにおける拍動に伴う動きの評価に基づいて、範囲に含まれる時相の中から拍動に伴う動きが相対的に小さい時相であり、かつ、時相間の時間間隔が大きくなるように複数の時相を選択する。
【００４９】
ここで、選択機能３７ｃは、各ＣＴ画像データにおける画像指標に基づいて、心臓の心拍に伴う動きを評価し、拍動に伴う動きが相対的に小さい時相であり、かつ、時相間の時間間隔が大きくなるように複数の時相を選択する。すなわち、選択機能３７ｃは、各ＣＴ画像データについて、画像指標をそれぞれ取得する。そして、選択機能３７ｃは、各画像データにおける画像指標に基づいて、複数のＣＴ画像データのうち、所定の時間間隔よりも大きい時間間隔を有するＣＴ画像データの組み合わせの中からＣＴ画像データの組を抽出する。例えば、選択機能３７ｃは、ＣＴ画像データの組み合わせの中から、各ＣＴ画像データにおける画質指標が所定の閾値以上であるＣＴ画像データの組を抽出する。ここで、画質指標は、画像指標の一例である。以下、画質指標に基づく、複数時相の選択（ＣＴ画像データの組の選択）の例を説明する。
【００５０】
図４は、第１の実施形態に係る選択機能３７ｃによる時相の選択を説明するための図である。図４においては、上段に心拍を示し、中段に心臓の動きを示し、下段に冠動脈の面積を示す。また、図４は、横方向が時間を示し、心拍、心臓の動き及び冠動脈の面積の時間変化を対応付けて示す。例えば、選択機能３７ｃは、心位相７０％〜９９％の範囲に含まれる心位相を選択する。ここで、心位相７０％〜９９％は、図４に示すように、心臓の動きがあまりなく、冠動脈の面積の変化が大きい時相である。心臓は収縮と拡張により動き、図４の中段に示すように、拡張期後半（心位相７０％〜９９％）で動きが安定する。選択機能３７ｃは、この動きが安定した心位相７０％〜９９％に含まれる心位相を選択することにより、拍動に伴う動きが相対的に小さい時相を選択する。なお、心位相７０％〜９９％は、Wave  Free  Periodに相当する。Wave  Free  Periodは、脈から冠動脈を下降する駆出圧波と、左室圧により形成され冠動脈を逆行性に伝わる反射圧の影響が無い、拡張中期の時相である。
【００５１】
また、図４の下段に示すように、冠動脈の面積は、心位相７０％付近で最大となり、９９％付近で最小となる。これは、心位相７０％付近で冠動脈に血液が流入し始め、その後、９９％に進むにつれて血液が流出していくためである。選択機能３７ｃは、この冠動脈の面積の変化をできるだけ含むように、心位相７０％〜９９％の範囲において、時相間の時間間隔ができるだけ大きくなるように複数の時相を選択する。
【００５２】
ここで、心位相７０〜９９％では、心臓の動きがあまりなく安定しているが、実際に収集したＣＴ画像データでは画質が変動する。また、心位相９９％付近では、心臓の次の収縮の影響を受けてしまい、ぼけが生じる場合もある。図５は、第１の実施形態に係る各心位相の画質の変化の一例を示す図である。図５は、横軸に時間を示し、縦軸に画質を示す。例えば、ＣＴ画像データの画質は、図５に示すように、心位相７０％付近に向けて徐々に向上し、ある画質付近で変動しながら推移する。そして、ＣＴ画像データの画質は、心位相９９％付近で急激に低下する。このように、ＣＴ画像データの画質は、心位相７０％〜９９％の範囲内でも一定ではなく、大きく変化する。
【００５３】
そこで、選択機能３７ｃは、拍動に伴う動きが相対的に小さい心位相７０％〜９９％に対応する各ＣＴ画像データの拍動に伴う動きをそれぞれ評価することにより、心位相７０％〜９９％の中から拍動に伴う動きがより小さい心位相を選択する。ここで、選択機能３７ｃは、心位相７０％〜９９％の中から拍動に伴う動きがより小さい心位相で、かつ、時間間隔がより大きくなるように心位相を選択する。以下、心位相７０％〜９９％の中から４つの心位相を選択する場合を例に挙げて、選択機能３７ｃの処理を説明する。
【００５４】
かかる場合には、選択機能３７ｃは、複数の時相において隣接する心位相間の各時間間隔が大きく、かつ、相互に近似するように、複数の心位相を選択する。例えば、選択機能３７ｃは、心位相７０％〜９９％の中から４つの心位相を選択する場合に、各心位相におけるＣＴ画像データが拍動に伴う動きがより小さく、かつ、隣り合う心位相間での時間間隔が等間隔に近くなるように、４つの心位相を選択する。これは、冠動脈の流体解析が、上述した冠動脈の面積がどのように変動するか（変動の傾向）を含めて解析するため、流体解析に適用するＣＴ画像データの心位相はできるだけ等間隔のものが望ましいためである。
【００５５】
以下、選択機能３７ｃによる４つの心位相の選択処理の一例を説明する。例えば、選択機能３７ｃは、範囲（Wave  Free  Period）を複数の部分範囲に分割し、各部分範囲について、部分範囲に含まれる時相の中から拍動に伴う動きが相対的に小さい時相をそれぞれ選択する。すなわち、選択機能３７ｃは、心位相７０％〜９９％を４つの領域に分割し、分割した各領域から１つずつ心位相を選択する。図６は、第１の実施形態に係る選択機能３７ｃによる選択処理の一例を説明するための図である。例えば、選択機能３７ｃは、心位相７０％〜９９％を、「７０〜７６％」、「７７％〜８３％」、「８４％〜９１％」、「９２〜９９％」の４つの領域に分割し、各領域で最適位相を検索して選択する。そして、選択機能３７ｃによって各領域から選択された４つの心位相に対応する各心位相データから最適な４つの心位相のＣＴ画像データが再構成される。
【００５６】
なお、図６に示す各領域における最適心位相の検索は、例えば、拍動に伴う動きがより小さい心位相を検索することである。例えば、選択機能３７ｃは、心電図同期で収集された心位相７０％〜９０％の投影データにおける１％ごとの心位相データを用いて画像内のぼけ（鮮明度）の程度を評価し、画像内のぼけが最も少ない心位相データを再構成する心位相データとして選択する。
【００５７】
図７は、第１の実施形態に係る選択機能３７ｃによる最適心位相の検索処理の一例を説明するための図である。図７においては、「７０〜７６％」の最適心位相を検索する場合について示す。例えば、選択機能３７ｃは、複数心拍分の投影データから冠動脈を含む１スライスのサイノグラムをそれぞれ抽出する。そして、選択機能３７ｃは、図７に示すように、抽出した各サイノグラムから「７０％」の心位相データ、「７１％」の心位相データ、「７２％」の心位相データ、「７３％」の心位相データ、「７４％」の心位相データ、「７５％」の心位相データ及び「７６％」の心位相データをそれぞれ抽出して再構成を行うことで、各心位相における冠動脈を含む１スライスの画像を生成する。
【００５８】
そして、選択機能３７ｃは、生成した各心位相の画像についてＳＤ値をそれぞれ算出し、算出したＳＤ値に基づいて、ぼけの少ない（鮮明度の高い）画像を抽出し、抽出した画像の心位相を最適心位相として選択する。例えば、選択機能３７ｃは、図７に示すように、各画像の画素ごとの信号強度（画素値）の分布において裾の幅が狭い画像の心位相「７４％」を、「７０〜７６％」の領域における最適心位相として選択する。なお、上述したＳＤ値とは、画像に含まれる画素ごとの信号強度の標準偏差（Standard  Deviation）であり、ノイズの特性を示す指標値である。
【００５９】
選択機能３７ｃは、「７７％〜８３％」、「８４％〜９１％」、「９２〜９９％」の３つの領域に対しても、上述した選択処理を実行して、各領域から最適心位相を選択する。ここで、上述したように各領域から最適心位相を選択すると、各領域で選択された最適心位相の時間間隔が短くなる場合がある。例えば、「７０〜７６％」の領域の最適心位相を「７４％」と選択し、「７７％〜８３％」の領域の最適心位相を「７７％」と選択した場合、時間間隔が「３％」分となる。上述したように、流体解析に用いるＣＴ画像データでは、心位相間の時間間隔が大きく、均等であるほうが望ましい。
【００６０】
そこで、選択機能３７ｃは、範囲に含まれる心位相において、心位相間の時間間隔が所定の閾値以上となるように、複数の心位相を選択する。図８は、第１の実施形態に係る選択機能３７ｃによる閾値処理の一例を説明するための図である。例えば、選択機能３７ｃは、図８に示すように、時間間隔の閾値として「５％」を用いて選択処理を実行する。一例を挙げると、選択機能３７ｃは、「７０〜７６％」の領域の最適心位相を「７４％」と選択した場合、「７４％」から「５％」以上間隔をあけた「７９％」を下端として「７７〜８３％」の領域の最適心位相を検索する。すなわち、選択機能３７ｃは、「７９〜８３％」の領域で最適心位相を検索する。なお、時間間隔の閾値は、任意に設定される。
【００６１】
上述したように、選択機能３７ｃが、「７０〜７６％」、「７７％〜８３％」、「８４％〜９１％」、「９２〜９９％」の４つの領域から最適心位相を選択すると、取得機能３７ｄは、選択機能３７ｃによって選択された複数の時相に対応する心位相データを取得する。具体的には、取得機能３７ｄは、選択された各最適心位相のＣＴ画像データを再構成するための心位相データを取得する。すなわち、取得機能３７ｄは、選択された各最適心位相について、体軸方向の各列における心位相データを取得する。画像再構成回路３６は、取得された４つの最適心位相に対応する心位相データを用いて、４つの最適心位相の３次元のＣＴ画像データを再構成する。解析機能３７ｂは、再構成された４つの心位相のＣＴ画像データを用いて上述した流体解析を実行する。
【００６２】
ここで、第１の実施形態に係るＸ線ＣＴ装置１００では、選択機能３７ｃによって選択された最適心位相が所定の条件に合わない場合に、警告を提示することも可能である。具体的には、提示機能３７ｅは、選択機能３７ｃによって選択された複数の心位相が所定の条件に適合しない場合に、操作者に対して警告を提示する。例えば、提示機能３７ｅは、選択された最適心位相において、上端及び下端の心位相間の時間間隔が所定の閾値よりも短い場合に、ディスプレイ３２に警告を表示させる。この場合、提示機能３７ｅは、警告を表示させるとともに、代替案を提示することもできる。例えば、提示機能３７ｅは、上端及び下端の心位相間の時間間隔が所定の閾値よりも長くなる心位相で、拍動に伴う動きが最も小さい心位相を代替案として提示する。
【００６３】
上述した実施形態では、心位相「７０％〜９９％」の範囲を４つの領域に分割して処理を実行する場合について説明した。しかしながら、実際のスキャンにおいて心位相「７０％〜９９％」の範囲の心位相データが収集されない場合もある。そこで、選択機能３７ｃは、実際に収集された心位相の範囲を分割することもできる。図９は、第１の実施形態に係る選択機能による選択処理の一例を説明するための図である。例えば、選択機能３７ｃは、心電同期スキャンにおいて心位相７０％〜９９％で収集された投影データのうち、実際に収集された心位相の範囲を４つの領域に分割し、各領域で最適位相を検索して選択する。そして、選択機能３７ｃによって各領域から選択された４つの心位相に対応する各心位相データから最適な４つの心位相のＣＴ画像データが再構成される。
【００６４】
一例を挙げると、実際に収集された心位相の範囲が「６８〜９５％」の場合、選択機能３７ｃは、「６８〜９５％」を４つの領域に分割する。ここで、選択機能３７ｃは、各領域に含まれる心位相が「７０％〜９９％」の範囲を超えないように領域を分割する。例えば、選択機能３７ｃは、「６８〜９５％」のうち、「６８％」及び「６９％」を対象とせずに、「６８〜９５％」を「７０〜７６％」、「７７〜８２％」、「８３〜８８％」及び「８９〜９５％」の４つの領域に分割する。
【００６５】
また、上述した実施形態では、心位相データから１スライス分の画像を生成して最適心位相検索を行う場合を例に挙げて説明した。しかしながら、実施形態はこれに限定されるものではなく、例えば、各心位相に対応するＣＴ画像データを再構成して最適心位相検索を行う場合であってもよい。この場合、例えば、選択機能３７ｃは、再構成したＣＴ画像データから冠動脈を含むスライスの画像を生成して、生成したスライスのＳＤ値を比較することにより最適心位相を検索することができる。また、例えば、選択機能３７ｃは、再構成したＣＴ画像データから冠動脈の３次元領域を抽出し、抽出した冠動脈の３次元領域（各ボクセル）についてＳＤ値を比較することにより最適心位相を検索することもできる。
【００６６】
また、上述した実施形態では、各領域からそれぞれ最適心位相を選択する場合を例に挙げて説明した。しかしながら、実施形態はこれに限定されるものではなく、例えば、範囲の上端及び下端が固定で選択されてもよい。一例を挙げると、選択機能３７ｃは、心位相「７０％〜９９％」の範囲において、心位相「７０％」及び「９９％」を最適心位相として固定で選択する場合であってもよい。或いは、選択機能３７ｃは、複数の心位相として、少なくとも範囲における上端付近の心位相及び下端付近の心位相を選択する場合であってもよい。一例を挙げると、選択機能３７ｃは、心位相「７０％〜９９％」の範囲において、心位相「７０％」付近及び「９９％」付近を最適心位相として優先的に選択する場合であってもよい。
【００６７】
次に、第１の実施形態に係るＸ線ＣＴ装置１００による処理の手順について説明する。図１０は、第１の実施形態に係るＸ線ＣＴ装置１００による処理手順を示すフローチャートである。なお、図１０では、心位相７０〜９９％の範囲を分割する場合の処理の手順について示す。
【００６８】
ここで、図１０におけるステップＳ１０１は、例えば、処理回路３７が制御機能３７ａに対応するプログラムを記憶回路３５から呼び出して実行することにより実現される。また、ステップＳ１０２、ステップＳ１０３及びステップＳ１０７は、例えば、処理回路３７が選択機能３７ｃに対応するプログラムを記憶回路３５から呼び出して実行することにより実現される。また、ステップＳ１０４及びステップＳ１０６は、例えば、処理回路３７が提示機能３７ｅに対応するプログラムを記憶回路３５から呼び出して実行することにより実現される。また、ステップＳ１０５は、例えば、画像再構成回路３６が、対応するプログラムを記憶回路３５から読み出して実行することにより実現される。また、ステップＳ１０８は、例えば、処理回路３７が解析機能３７ｂに対応するプログラムを記憶回路３５から呼び出して実行することにより実現される。
【００６９】
本実施形態に係るＸ線ＣＴ装置１００では、まず、処理回路３７が、投影データを収集する（ステップＳ１０１）。そして、処理回路３７が、心位相７０〜９９％を複数の領域に分割し（ステップＳ１０２）、各領域において、冠動脈の動きが少ない心位相を選択する（ステップＳ１０３）。続いて、処理回路３７は、選択された心位相が条件を満たすか否かを判定する（ステップＳ１０４）。ここで、選択された心位相が条件を満たす場合（ステップＳ１０４肯定）、画像再構成回路３６が、選択された心位相に対応するボリュームデータを再構成する（ステップＳ１０５）。
【００７０】
一方、選択された心位相が条件を満たさない場合（ステップＳ１０４否定）、処理回路３７は、警告を表示する（ステップＳ１０６）。そして、処理回路３７は、再選択の操作を受け付けたか否かを判定する（ステップＳ１０７）。ここで、再選択の操作を受け付けた場合（ステップＳ１０７肯定）、処理回路３７は、ステップＳ１０３に戻って、冠動脈の動きが少ない心位相を選択する。一方、再選択の操作を受け付けなかった場合（ステップＳ１０７否定）、画像再構成回路３６が、選択された心位相に対応するボリュームデータを再構成する（ステップＳ１０５）。ステップＳ１０５において、各心位相に対応するボリュームデータが再構成されると、処理回路３７は、再構成された各心位相のボリュームデータを用いて流体解析を実行する（ステップＳ１０８）。
【００７１】
上述したように、第１の実施形態によれば、データ収集回路１４は、経時的な複数の投影データを収集する。選択機能３７ｃは、複数の投影データに対応する各時相から、拍動に伴う動きが相対的に小さい心位相であり、かつ、時相間の時間間隔が大きくなるように複数の心位相を選択する。取得機能３７ｄは、選択機能３７ｃによって選択された複数の心位相に対応する投影データを取得する。従って、第１の実施形態に係るＸ線ＣＴ装置１００は、流体解析を実行する際の適切なデータを簡便に取得することを可能にする。
【００７２】
また、第１の実施形態によれば、選択機能３７ｃは、複数の投影データに対応する各時相において拍動に伴う動きが相対的に小さい心位相の範囲を抽出し、抽出した範囲に含まれる心位相に対応する画像データにおける拍動に伴う動きの評価に基づいて、範囲に含まれる心位相の中から拍動に伴う動きが相対的に小さい心位相であり、かつ、心位相間の時間間隔が大きくなるように複数の心位相を選択する。従って、第１の実施形態に係るＸ線ＣＴ装置１００は、流体解析を実行する際のより適切なデータを取得することを可能にする。
【００７３】
また、第１の実施形態によれば、選択機能３７ｃは、投影データに含まれる冠動脈領域における鮮明度に基づいて、拍動に伴う動きを評価する。従って、第１の実施形態に係るＸ線ＣＴ装置１００は、拍動に伴う動きを容易に評価することを可能にする。
【００７４】
また、第１の実施形態によれば、選択機能３７ｃは、範囲を複数の部分範囲に分割し、各部分範囲について、部分範囲に含まれる心位相の中から拍動に伴う動きが相対的に小さい心位相をそれぞれ選択する。また、選択機能３７ｃは、複数の部分範囲の幅がそれぞれ近似するように、範囲を部分範囲に分割する。従って、第１の実施形態に係るＸ線ＣＴ装置１００は、ある程度の間隔を持った心位相を容易に選択することを可能にする。
【００７５】
また、第１の実施形態によれば、選択機能３７ｃは、複数の心位相として、少なくとも範囲における上端付近の心位相及び下端付近の心位相を選択する。従って、第１の実施形態に係るＸ線ＣＴ装置１００は、冠動脈の面積が最大となる心時相及び最小となる心時相に近い心位相を容易に選択することを可能にする。
【００７６】
また、第１の実施形態によれば、選択機能３７ｃは、範囲に含まれる心位相において、心位相間の時間間隔が所定の閾値以上となるように、複数の心位相を選択する。従って、第１の実施形態に係るＸ線ＣＴ装置１００は、心位相間の時間間隔を等間隔に近づけることを可能にする。
【００７７】
また、第１の実施形態によれば、選択機能３７ｃは、複数の心位相において隣接する心位相間の各時間間隔が大きく、かつ、相互に近似するように、複数の心位相を選択する。従って、第１の実施形態に係るＸ線ＣＴ装置１００は、心位相間の時間間隔をより大きく等間隔に近づけることを可能にする。
【００７８】
また、第１の実施形態によれば、提示機能３７ｅは、選択機能３７ｃによって選択された複数の心位相が所定の条件に適合しない場合に、操作者に対して警告を提示する。従って、第１の実施形態に係るＸ線ＣＴ装置１００は、あまり適切ではないＣＴ画像データを用いた流体解析の実行を事前に抑止することを可能にする。
【００７９】
（第２の実施形態）
次に、第２の実施形態について説明する。なお、本実施形態に係るＸ線ＣＴ装置１００の構成は、基本的には、図２に示したＸ線ＣＴ装置１００の構成と同じである。そのため、以下では、第１の実施形態に係るＸ線ＣＴ装置１００と異なる点を中心に説明することとし、図２に示した構成要素と同様の役割を果たす構成要素については同じ符号を付すこととして詳細な説明を省略する。
【００８０】
上述した第１の実施形態では、心位相の範囲を複数の領域に分割して最適心位相を選択する場合を例に挙げて説明した。第２の実施形態では、心位相の範囲を分割せずに最適心位相を選択する場合について説明する。例えば、第２の実施形態に係るＸ線ＣＴ装置１００は、全心位相を対象に拍動に伴う動きをそれぞれ評価し、評価結果に基づいて最適心位相を選択する。
【００８１】
図１１は、第２の実施形態に係る選択機能３７ｃによる選択処理の一例を説明するための図である。図１１においては、心位相７０〜９９％の範囲から最適心位相を選択する場合を示す。例えば、選択機能３７ｃは、心位相７０〜９９％の範囲における１％ごとの心位相について、上述したＳＤ値をそれぞれ算出する。そして、選択機能３７ｃは、図１１に示すように、各心位相に対してＳＤ値が低い順に順位を付けて最適心位相を選択する。ここで、選択機能３７ｃは、単純に順位の高い心位相から順に選択するのではなく、心位相間の時間間隔も考慮して最適心位相を選択する。一例を挙げると、心位相間の時間間隔が「５％」以上離れ、かつ、ＳＤ値による評価の順位が高くなるように最適心位相を選択する。
【００８２】
なお、第２の実施形態においても、範囲(例えば、７０〜９９％)の上端付近及び下端付近を優先的に最適心位相として選択する場合であってもよい。この場合には、選択機能３７ｃは、まず、範囲の上端付近及び下端付近でそれぞれ順位が高い心位相を最適心位相として選択する。そして、選択機能３７ｃは、選択した最適心位相から所定の時間間隔（例えば、「５％」）以上離れた心位相からさらに最適心位相を選択する。
【００８３】
次に、第２の実施形態に係るＸ線ＣＴ装置１００による処理の手順について説明する。図１２は、第２の実施形態に係るＸ線ＣＴ装置１００による処理手順を示すフローチャートである。なお、図１２では、心位相７０〜９９％の範囲から最適心位相を選択する場合の処理の手順について示す。
【００８４】
ここで、図１２におけるステップＳ２０１は、例えば、処理回路３７が制御機能３７ａに対応するプログラムを記憶回路３５から呼び出して実行することにより実現される。また、ステップＳ２０２及びステップＳ２０３は、例えば、処理回路３７が選択機能３７ｃに対応するプログラムを記憶回路３５から呼び出して実行することにより実現される。また、ステップＳ２０４は、例えば、画像再構成回路３６が、対応するプログラムを記憶回路３５から読み出して実行することにより実現される。また、ステップＳ２０５は、例えば、処理回路３７が解析機能３７ｂに対応するプログラムを記憶回路３５から呼び出して実行することにより実現される。
【００８５】
本実施形態に係るＸ線ＣＴ装置１００では、まず、処理回路３７が、投影データを収集する（ステップＳ２０１）。そして、処理回路３７が、心位相７０〜９９％の各心位相における動きを評価する（ステップＳ２０２）。そして、処理回路３７が、動きの少なさと時間間隔に基づいて、心位相を選択する（ステップＳ２０３）。その後、画像再構成回路３６が、選択された心位相に対応するボリュームデータを再構成する（ステップＳ２０４）。ステップＳ２０４において、各心位相に対応するボリュームデータが再構成されると、処理回路３７は、再構成された各心位相のボリュームデータを用いて流体解析を実行する（ステップＳ２０５）。
【００８６】
上述したように、第２の実施形態によれば、選択機能３７ｃは、拍動に伴う動きが相対的に小さい心位相の範囲に含まれる各心位相について拍動に伴う動きを評価して、最適心位相を選択する。従って、第２の実施形態に係るＸ線ＣＴ装置１００は、流体解析を実行する際の適切なデータを取得することを可能にする。
【００８７】
（第３の実施形態）
さて、これまで第１及び第２の実施形態について説明したが、上述した第１及び第２の実施形態以外にも、種々の異なる形態にて実施されてよいものである。
【００８８】
上述した第１及び第２の実施形態で説明した最適心位相の選択処理は、任意のタイミングで実行することができる。例えば、選択機能３７ｃは、データ転送時に複数の最適心位相を選択する。一例を挙げると、選択機能３７ｃは、再構成済みのボリュームデータを医用情報処理装置３００などに転送する際に、ボリュームデータを用いて最適心位相を選択する。或いは、選択機能３７ｃは、流体解析のアプリケーションの起動時に複数の最適心位相を選択する。一例を挙げると、選択機能３７ｃは、入力回路３１を介して流体解析のアプリケーションを起動させるための操作を受け付けた場合に、最適心位相を選択する。この場合、選択機能３７ｃは、心位相データを用いた選択処理を実行する場合であってもよく、ボリュームデータ（ＣＴ画像データ）を用いた選択処理を実行する場合であってもよい。
【００８９】
また、上述した実施形態では、画像指標として、画質指標を用いる場合を例に挙げて説明した。しかしながら、実施形態はこれに限定されるものではなく、Ｘ線ＣＴ装置１００は、種々の画像指標を用いることができる。例えば、Ｘ線ＣＴ装置１００は、画像指標として部位の抽出能を用いることができる。Ｘ線ＣＴ装置１００は、上述した血管形状データなどの部位の形状の抽出処理を行うことができる。ここで、ＣＴ画像データからの部位の抽出の精度は、収集されたＣＴ画像データの状態に応じて変化する。例えば、心臓の拍動に伴う動きが大きい場合、部位の抽出の精度が低下する場合がある。すなわち、第３の実施形態に係るＸ線ＣＴ装置１００は、このような画像データに起因する部位の抽出能を画像指標として用いる。
【００９０】
具体的には、第３の実施形態に係る選択機能３７ｃは、各ＣＴ画像データに含まれる所定の部位の抽出結果に基づいて、ＣＴ画像データの組み合わせの中から流体解析に用いるＣＴ画像データの組を抽出する。ここで、選択機能３７ｃは、所定の部位として、例えば、冠動脈の血管壁又は冠動脈の芯線を用いる。以下、これらの部位の抽出能に基づく、複数時相の選択（ＣＴ画像データの組の選択）の例を説明する。
【００９１】
ここで、まず、第３の実施形態に係るＸ線ＣＴ装置１００では、部位の抽出に際して、各時相のＣＴ画像データを再構成する。すなわち、画像再構成回路３６は、収集された複数心拍分の投影データから各心位相に対応するＣＴ画像データを再構成する。例えば、画像再構成回路３６は、心位相７０％〜９９％の各心位相データから各心位相のＣＴ画像データをそれぞれ再構成する。
【００９２】
各心位相のＣＴ画像データが再構成されると、解析機能３７ｂが、各心位相のＣＴ画像データから部位を抽出する。ここで、解析機能３７ｂは、種々の手法により、各心位相のＣＴ画像データから部位を抽出することができる。例えば、解析機能３７ｂは、各ＣＴ画像データに対してそれぞれ抽出処理を施すことで、血管壁や芯線を抽出する。すなわち、解析機能３７ｂは、画像再構成回路３６によって再構成された全てのＣＴ画像データに対して部位抽出を行うことで、全てのＣＴ画像データからそれぞれ部位を抽出する。なお、解析機能３７ｂは、部位抽出に用いる抽出アルゴリズムとして、既存のアルゴリズムを適宜利用することができる。
【００９３】
また、解析機能３７ｂは、１つの心位相のＣＴ画像データから部位を抽出し、部位を抽出したＣＴ画像データと他の心位相のＣＴ画像データとを位置合わせすることで、他の心位相のＣＴ画像データにおける部位を抽出すうこともできる。かかる場合には、解析機能３７ｂは、まず、１つの心位相のＣＴ画像データに対して抽出処理を施すことで、血管壁や芯線を抽出する。次に、解析機能３７ｂは、血管壁や芯線を抽出したＣＴ画像データと他の心位相のＣＴ画像データとを位置合わせする。例えば、解析機能３７ｂは、位置合わせ処理を行うＣＴ画像データからそれぞれ特徴点を抽出し、抽出した特徴点を一致させるために変形させる非剛体の位置合わせ処理を実行する。
【００９４】
一例を挙げると、解析機能３７ｂは、血管壁や芯線を抽出したＣＴ画像データと他の心位相のＣＴ画像データから、対応する特徴点（解剖学的な特徴点）をそれぞれ抽出する。そして、解析機能３７ｂは、抽出した各特徴点を一致させるための座標変換情報を算出する。例えば、解析機能３７ｂは、他の心位相のＣＴ画像データにおける特徴点を、血管壁や芯線を抽出したＣＴ画像データにおける特徴点に一致させるための３次元座標系における座標変換情報を算出する。すなわち、解析機能３７ｂは、他の心位相のＣＴ画像データにおいて、血管壁や芯線を抽出したＣＴ画像データの各位置に対応する位置を算出するための変換情報を算出する。解析機能３７ｂは、抽出した血管壁や芯線の座標に対して変換情報を用いることで、他の心位相のＣＴ画像データにおける血管壁や芯線の座標を特定することできる。
【００９５】
上述したように、部位が抽出されると、選択機能３７ｃは、抽出結果に基づいて、複数時相（ＣＴ画像データの組）を選択する。以下、冠動脈の血管壁の抽出結果を用いる場合の例と、冠動脈の芯線の抽出結果を用いる場合の例とを順に説明する。例えば、選択機能３７ｃは、解析機能３７ｂによって抽出された血管壁の断面積の変動に基づいて複数時相を選択する。図１３は、第３の実施形態に係る選択機能３７ｃによる時相の選択を説明するための図である。ここで、図１３においては、横軸に時相を示し、横軸に断面積を示す。例えば、選択機能３７ｃは、図１３のグラフ中の点に示すように、解析機能３７ｂによって抽出された血管壁における所定の断面の断面積を、各ＣＴ画像データについてそれぞれ算出する。そして、選択機能３７ｃは、冠動脈の断面積の一般的な変化を示す直線Ｌ１の傾きと、算出した各心位相の断面積とを比較して、断面積の経時変化として不適切な心位相のＣＴ画像データを選択の対象から除外する。例えば、選択機能３７ｃは、図１３に示す心位相Ｐ１のＣＴ画像データを選択の対象から除外する。すなわち、選択機能３７ｃは、隣接する断面積の値を相互に比較して、隣接する断面積の値からの乖離度が大きい心位相のＣＴ画像データを選択の対象から除外する。
【００９６】
一例を挙げると、選択機能３７ｃは、まず、「７０〜７６％」、「７７％〜８３％」、「８４％〜９１％」、「９２〜９９％」の４つの領域から、心位相間の時間間隔が所定の閾値以上となるように、それぞれ心位相を選択する。そして、選択機能３７ｃは、選択した各心位相における断面積を算出して相互に比較する。さらに、選択機能３７ｃは、直線Ｌ１の傾きに基づく許容範囲を超えて断面積が変化しているか否かを判定する。ここで、選択した各心位相における断面積の変化が直線Ｌ１の傾きに基づく許容範囲を超えていない場合、選択機能３７ｃは、「７０〜７６％」、「７７％〜８３％」、「８４％〜９１％」、「９２〜９９％」の４つの領域から選択した心位相を最適心位相として決定する。
一方、選択した各心位相における断面積の変化が直線Ｌ１の傾きに基づく許容範囲を超えている場合、選択機能３７ｃは、所定の条件に基づいて、心位相を再選択する。例えば、「７７％〜８３％」から選択した心位相における断面積の変化が許容範囲を超えている場合、選択機能３７ｃは、「７７％〜８３％」から再度心位相を選択する。ここで、選択機能３７ｃは、隣接する領域「７０〜７６％」及び「８４％〜９１％」において選択済みの心位相との時間間隔をできるだけ維持するように、「７７％〜８３％」から心位相を再選択する。
【００９７】
例えば、選択機能３７ｃは、心位相を再選択する場合に、不適切と判定した心位相の隣接する心位相を選択する。一例を挙げると、最初の選択において「７９％」を選択していた場合、選択機能３７ｃは、「７８％」又は「８０％」の心位相を選択する。
【００９８】
或いは、選択機能３７ｃは、心位相を再選択する場合に、所定の数の位相（例えば、３％分）を隔てた心位相を選択する。一例を挙げると、最初の選択において「７９％」を選択していた場合、選択機能３７ｃは、「８２％」の心位相を選択する。上述したように、ＣＴ画像データからの部位の抽出の精度は、動きが大きい場合に低下する場合がある。従って、不適切と判定した心位相の隣接する心位相を再選択した場合、再選択した心位相も不適切となる可能性がある。そこで、選択機能３７ｃは、上述したように所定の数の位相を隔てた心位相を選択することで、精度よく再選択を行うことができる。
【００９９】
なお、選択機能３７ｃは、心位相を再選択した場合に、４つの心位相間の時間間隔を再調整する場合であってもよい。例えば、「７９％」の心位相を不適切として、「８２％」の心位相を再選択した場合、選択機能３７ｃは、「８４％〜９１％」及び「９２〜９９％」から選択済みの心位相を再選択する。すなわち、選択機能３７ｃは、心位相間の時間間隔が広くなるように、「８４％〜９１％」における心位相と「９２〜９９％」における心位相を再選択する。
【０１００】
また、上述した例では、４つの領域から選択した各心位相における冠動脈の断面積を算出して、相互に比較する場合について説明した。しかしながら、実施形態はこれに限定されるものではなく、全ての心位相における冠動脈の断面線を算出して、相互に比較する場合であってもよい。かかる場合には、例えば、選択機能３７ｃは、全ての心位相における冠動脈の断面線を算出して、断面積の変化が直線Ｌ１の傾きに基づく許容範囲を超えていない心位相を、心位相間の時間間隔がより広くなるように選択する。
【０１０１】
また、上述した例では、冠動脈の血管壁の抽出について、断面積を用いて判定する場合について説明した。しかしながら、実施形態はこれに限定されるものではなく、例えば、冠動脈の血管壁が抽出されたか否かを用いる場合であってもよい。かかる場合には、例えば、選択機能３７ｃは、血管壁が抽出されなかったＣＴ画像データを選択の対象から除外する。また、例えば、選択機能３７ｃは、血管壁の輪郭の類似度を用いて血管壁の抽出の精度を判定することもできる。かかる場合には、例えば、選択機能３７ｃは、判定対象の心位相における血管壁の輪郭を両隣の心位相の血管壁の輪郭と比較して、類似度が共に所定の閾値を下回った場合に、判定対象の心位相を選択の対象から除外する。
【０１０２】
次に、冠動脈の芯線の抽出結果を用いる場合について説明する。かかる場合には、例えば、選択機能３７ｃは、解析機能３７ｂによって抽出された冠動脈の芯線の長さを比較することで、芯線の抽出の精度を判定する。例えば、全てのＣＴ画像データから芯線をそれぞれ抽出する場合、芯線の抽出が失敗して短くなる場合がある。選択機能３７ｃは、各心位相のＣＴ画像データから抽出された各芯線の長さをそれぞれ算出する。そして、選択機能３７ｃは、算出した芯線の長さを各心位相間で比較して、芯線の長さが他よりも短い心位相を選択の対象から除外する。
【０１０３】
上述した例では、抽出した血管壁や芯線に基づいて、選択機能３７ｃが自動で再選択を行う場合について説明した。しかしながら、実施形態はこれに限定されるものではなく、例えば、抽出結果を操作者に提示して、再選択を行う指示を受け付ける場合であってもよい。かかる場合には、例えば、提示機能３７ｅは、解析機能３７ｂによる部位の抽出処理の結果をディスプレイ３２に表示させる。操作者は、ディスプレイ３２に表示された抽出処理の結果を参照して、再選択を行うか否かを判断する。そして、操作者は、入力回路３１を介して再選択を行うか否かの指示を入力する。ここで、入力回路３１が再選択を行う旨の指示を受け付けた場合、選択機能３７ｃは、心位相を再選択する。一方、入力回路３１が再選択を行わない旨の指示を受け付けた場合、解析機能３７ｂは、選択済みの心位相のＣＴ画像データを用いた流体解析を実施する。
【０１０４】
以上、画像指標として部位の抽出能を用いる場合の例について説明した。次に、第３の実施形態に係るＸ線ＣＴ装置１００による処理の手順について説明する。図１４は、第３の実施形態に係るＸ線ＣＴ装置１００による処理手順を示すフローチャートである。なお、図１４では、画像指標として冠動脈の血管壁の抽出能を用いる場合の処理の手順について示す。
【０１０５】
ここで、図１４におけるステップＳ３０１は、例えば、処理回路３７が制御機能３７ａに対応するプログラムを記憶回路３５から呼び出して実行することにより実現される。また、ステップＳ３０２は、例えば、画像再構成回路３６が、対応するプログラムを記憶回路３５から読み出して実行することにより実現される。また、ステップＳ３０３及びステップ３０７は、例えば、処理回路３７が解析機能３７ｂに対応するプログラムを記憶回路３５から呼び出して実行することにより実現される。また、ステップＳ３０４、ステップＳ３０５及びステップＳ３１０は、例えば、処理回路３７が選択機能３７ｃに対応するプログラムを記憶回路３５から呼び出して実行することにより実現される。また、ステップＳ３０６及びステップＳ３０８は、例えば、処理回路３７が提示機能３７ｅに対応するプログラムを記憶回路３５から呼び出して実行することにより実現される。
【０１０６】
本実施形態に係るＸ線ＣＴ装置１００では、まず、処理回路３７が、投影データを収集する（ステップＳ３０１）。そして、画像再構成回路３６が、各心位相のボリュームデータを再構成する（ステップＳ３０２）。続いて、処理回路３７が、再構成された各心位相のボリュームデータから冠動脈の血管壁をそれぞれ抽出する（ステップＳ３０３）。そして、処理回路３７が、心位相７０〜９９％を複数の領域に分割し（ステップＳ３０４）、各領域において、心位相を選択する（ステップＳ３０５）。
【０１０７】
続いて、処理回路３７は、選択された心位相が条件を満たすか否かを判定する（ステップＳ３０６）。ここで、選択された心位相が条件を満たす場合（ステップＳ３０６肯定）、処理回路３７は、選択した各心位相のボリュームデータを用いて流体解析を実行する（ステップＳ３０７）。
【０１０８】
一方、選択された心位相が条件を満たさない場合（ステップＳ３０６否定）、処理回路３７は、警告を表示して（ステップＳ３０８）、再選択の操作を受け付けたか否かを判定する（ステップＳ３０９）。ここで、再選択の操作を受け付けた場合（ステップＳ３０９肯定）、処理回路３７は、条件に基づいて位相を再選択する（ステップＳ３１０）。一方、再選択の操作を受け付けなかった場合（ステップＳ３０９否定）、処理回路３７は、選択した各心位相のボリュームデータを用いて流体解析を実行する（ステップＳ３０７）。
【０１０９】
また、上述した実施形態では、Ｘ線ＣＴ装置１００が各種処理を実行する場合について説明した。しかしながら、実施形態はこれに限定されるものではなく、例えば、医用情報処理装置３００において各種処理が実行される場合であってもよい。図１３は、第３の実施形態に係る医用情報処理装置３００の構成の一例を示す図である。
【０１１０】
例えば、図１３に示すように、医用情報処理装置３００は、Ｉ／Ｆ（インターフェース）回路３１０と、記憶回路３２０と、入力回路３３０と、ディスプレイ３４０と、処理回路３５０とを有する。
【０１１１】
Ｉ／Ｆ回路３１０は、処理回路３５０に接続され、ネットワーク４００を介して接続された各種のＸ線ＣＴ装置１００又は画像保管装置２００との間で行われる各種データの伝送及び通信を制御する。例えば、Ｉ／Ｆ回路３１０は、ネットワークカードやネットワークアダプタ、ＮＩＣ（Network  Interface  Controller）等によって実現される。本実施形態では、Ｉ／Ｆ回路３１０は、Ｘ線ＣＴ装置１００又は画像保管装置２００からＣＴ画像データを受信し、受信したＣＴ画像データを処理回路３５０に出力する。
【０１１２】
記憶回路３２０は、処理回路３５０に接続され、各種データを記憶する。例えば、記憶回路３２０は、ＲＡＭ（Random  Access  Memory）、フラッシュメモリ等の半導体メモリ素子や、ハードディスク、光ディスク等によって実現される。本実施形態では、記憶回路３２０は、Ｘ線ＣＴ装置１００又は画像保管装置２００から受信したＣＴ画像データを記憶する。
【０１１３】
入力回路３３０は、処理回路３５０に接続され、操作者から受け付けた入力操作を電気信号に変換して処理回路３５０に出力する。例えば、入力回路３３０は、トラックボール、スイッチボタン、マウス、キーボード、タッチパネル等によって実現される。
【０１１４】
ディスプレイ３４０は、処理回路３５０に接続され、処理回路３５０から出力される各種情報及び各種画像データを表示する。例えば、ディスプレイ３４０は、液晶モニタやＣＲＴ（Cathode  Ray  Tube）モニタ、タッチパネル等によって実現される。
【０１１５】
処理回路３５０は、入力回路３３０を介して操作者から受け付けた入力操作に応じて、画像処理装置３００が有する各構成要素を制御する。例えば、処理回路３５０は、プロセッサによって実現される。本実施形態では、処理回路３５０は、Ｉ／Ｆ回路３１０から出力されるＣＴ画像データを記憶回路３２０に記憶させる。また、処理回路３５０は、記憶回路３２０からＣＴ画像データを読み出し、ディスプレイ３４０に表示する。
【０１１６】
そして、処理回路３５０は、図１３に示すように、制御機能３５１と、選択機能３５２と、取得機能３５３と、提示機能３５４と、解析機能３５５とを実行する。制御機能３５１は、医用情報処理装置３００の全体を制御する。選択機能３５２は、上述した選択機能３７ｃと同様の処理を実行する。取得機能３５３は、上述した取得機能３７ｄと同様の処理を実行する。提示機能３５４は、上述した提示機能３７ｅと同様の処理を実行する。解析機能３５５は、上述した解析機能３７ｂと同様の処理を実行する。
【０１１７】
上述した実施形態では、単一の処理回路（処理回路３７及び処理回路３５０）によって各処理機能が実現される場合の例を説明したが、実施形態はこれに限られない。例えば、処理回路３７及び処理回路３５０は、複数の独立したプロセッサを組み合わせて構成され、各プロセッサが各プログラムを実行することにより各処理機能を実現するものとしても構わない。また、処理回路３７及び処理回路３５０が有する各処理機能は、単一又は複数の処理回路に適宜に分散又は統合されて実現されてもよい。
【０１１８】
上述した各実施形態の説明で用いた「プロセッサ」という文言は、例えば、ＣＰＵ（Central  Processing  Unit）、ＧＰＵ(Graphics  Processing  Unit)、或いは、特定用途向け集積回路（Application  Specific  Integrated  Circuit：ＡＳＩＣ）、プログラマブル論理デバイス（例えば、単純プログラマブル論理デバイス（Simple  Programmable  Logic  Device：ＳＰＬＤ）、複合プログラマブル論理デバイス（Complex  Programmable  Logic  Device：ＣＰＬＤ）、及びフィールドプログラマブルゲートアレイ（Field  Programmable  Gate  Array：ＦＰＧＡ））等の回路を意味する。ここで、記憶回路にプログラムを保存する代わりに、プロセッサの回路内にプログラムを直接組み込むように構成しても構わない。この場合には、プロセッサは回路内に組み込まれたプログラムを読み出し実行することで機能を実現する。また、本実施形態の各プロセッサは、プロセッサごとに単一の回路として構成される場合に限らず、複数の独立した回路を組み合わせて一つのプロセッサとして構成され、その機能を実現するようにしてもよい。
【０１１９】
ここで、プロセッサによって実行されるプログラムは、ＲＯＭ（Read  Only  Memory）や記憶部等に予め組み込まれて提供される。なお、このプログラムは、これらの装置にインストール可能な形式又は実行可能な形式のファイルでＣＤ（Compact  Disk）−ＲＯＭ、ＦＤ（Flexible  Disk）、ＣＤ−Ｒ（Recordable）、ＤＶＤ（Digital  Versatile  Disk）等のコンピュータで読み取り可能な記憶媒体に記録されて提供されてもよい。また、このプログラムは、インターネット等のネットワークに接続されたコンピュータ上に格納され、ネットワーク経由でダウンロードされることにより提供又は配布されてもよい。例えば、このプログラムは、各機能部を含むモジュールで構成される。実際のハードウェアとしては、ＣＰＵが、ＲＯＭ等の記憶媒体からプログラムを読み出して実行することにより、各モジュールが主記憶装置上にロードされて、主記憶装置上に生成される。
【０１２０】
以上説明した少なくとも一つの実施形態によれば、適切なデータを簡便に選択することができる。
【０１２１】
本発明のいくつかの実施形態を説明したが、これらの実施形態は、例として提示したものであり、発明の範囲を限定することは意図していない。これら実施形態は、その他の様々な形態で実施されることが可能であり、発明の要旨を逸脱しない範囲で、種々の省略、置き換え、変更を行うことができる。これら実施形態やその変形は、発明の範囲や要旨に含まれると同様に、特許請求の範囲に記載された発明とその均等の範囲に含まれるものである。
【符号の説明】
【０１２２】
１００  Ｘ線ＣＴ装置
３７ｂ、３５５  解析機能
３７ｃ、３５２  選択機能
３７ｄ、３５３  取得機能
３７ｅ、３５４  提示機能
３００  医用情報処理装置
【図１】
2017217459.tif 000003
【図２】
2017217459.tif 000004
【図３】
2017217459.tif 000005
【図４】
2017217459.tif 000006
【図５】
2017217459.tif 000007
【図６】
2017217459.tif 000008
【図７】
2017217459.tif 000009
【図８】
2017217459.tif 000010
【図９】
2017217459.tif 000011
【図１０】
2017217459.tif 000012
【図１１】
2017217459.tif 000013
【図１２】
2017217459.tif 000014
【図１３】
2017217459.tif 000015
【図１４】
2017217459.tif 000016
【図１５】
2017217459.tif 000017