(57)【要約】【課題】パラメトリックイメージングにおいて立体視可能な画像を提供すること。【解決手段】実施形態のＸ線診断装置は、第１の方向から照射されるＸ線を用いて、所定の撮像条件で時系列に撮像した造影画像データと非造影画像データを差分して第１の差分画像データを生成し、第２の方向から照射されるＸ線を用いて、所定の撮像条件と略同一の撮像条件で時系列に撮像した造影画像データと非造影画像データを差分して第２の差分画像データを生成する。Ｘ線診断装置は、第１の差分画像データの各位置に、画素値の経時的変化に応じたカラーを割り当てた第１のカラー画像データを所定の処理条件で生成し、第２の差分画像データの各位置に、画素値の経時的変化に応じたカラーを割り当てた第２のカラー画像データを所定の処理条件と略同一の処理条件で生成する。Ｘ線診断装置は、第１のカラー画像データ及び第２のカラー画像データに基づく立体視画像を表示させる。【選択図】図１
2018000320.tif 000002
【特許請求の範囲】
【請求項１】
第１の方向から照射されるＸ線を用いて、所定の撮像条件で時系列に撮像した造影画像データと非造影画像データを差分して第１の差分画像データを生成し、第２の方向から照射されるＸ線を用いて、前記所定の撮像条件と略同一の撮像条件で時系列に撮像した造影画像データと非造影画像データを差分して第２の差分画像データを生成する差分画像生成部と、
前記第１の差分画像データの各位置に、画素値の経時的変化に応じたカラーを割り当てた第１のカラー画像データを所定の処理条件で生成するとともに、前記第２の差分画像データの各位置に、画素値の経時的変化に応じたカラーを割り当てた第２のカラー画像データを前記所定の処理条件と略同一の処理条件で生成するカラー画像生成部と、
前記第１のカラー画像データ及び前記第２のカラー画像データに基づく立体視画像を表示させる表示制御部と
を備える、Ｘ線診断装置。
【請求項２】
前記差分画像生成部は、血管の走行方向に応じて前記第１の方向と前記第２の方向との間の視差の角度方向を導出し、導出した前記第１の方向及び前記第２の方向のそれぞれで、前記第１の差分画像データ及び前記第２の差分画像データを生成する、
請求項１に記載のＸ線診断装置。
【請求項３】
前記カラー画像生成部は、第１の差分画像データの各位置における前記画素値の経時的変化に基づいて、被検体に注入される造影剤の流入時間を同定し、同定した流入時間に応じて前記カラーを各位置に割り当てて前記第１のカラー画像データを生成するとともに、第２の差分画像データの各位置における前記画素値の経時的変化に基づいて、被検体に注入される造影剤の流入時間を同定し、同定した流入時間に応じて前記カラーを各位置に割り当てて前記第２のカラー画像データを生成する、
請求項１又は２に記載のＸ線診断装置。
【請求項４】
前記差分画像生成部は、前記造影剤の注入条件、前記Ｘ線の照射条件、前記Ｘ線が照射される視野サイズ、前記Ｘ線を照射するＸ線管球の焦点と前記Ｘ線の検出器との間の距離、前記被検体が載置される寝台の位置、及び、前記Ｘ線管球及び前記検出器を支持する支持機の回転中心のうち、少なくとも一つが略同一である前記撮像条件を用いて、前記第１の差分画像データ及び前記第２の差分画像データを生成し、
前記カラー画像生成部は、前記流入時間の同定方法、及び、前記経時的変化に応じて割り当てられるカラーが規定されたカラーコードのうち、少なくとも一つが略同一である前記処理条件を用いて、前記第１のカラー画像データ及び前記第２のカラー画像データを生成する、
請求項３に記載のＸ線診断装置。
【請求項５】
前記差分画像生成部は、前記造影剤の注入開始から前記造影画像データの撮像開始までの時間、前記造影剤の注入速度、前記造影剤の注入器に加える圧力、及び、前記注入速度に到達する前の立ち上がり曲線のうち、少なくとも一つが略同一である前記注入条件を用いて、前記第１の差分画像データ及び前記第２の差分画像データを生成する、
請求項４に記載のＸ線診断装置。
【請求項６】
前記差分画像生成部は、前記Ｘ線による撮像のフレームレート、前記Ｘ線管球の管電圧、前記Ｘ線管球の管電流、前記Ｘ線のパルス幅、前記Ｘ線管球のコリメータの位置、前記Ｘ線管球のビームフィルタの種類及び厚み、及び、前記Ｘ線管球の補償フィルタの位置のうち、少なくとも一つが略同一である前記照射条件を用いて、前記第１の差分画像データ及び前記第２の差分画像データを生成する、
請求項４に記載のＸ線診断装置。
【請求項７】
前記カラー画像生成部は、前記第１の差分画像データ及び前記第２の差分画像データの同一位置から求まる前記流入時間を揃える補正処理を行って、前記第１のカラー画像データ及び前記第２のカラー画像データのうち少なくとも一方を再生成する、
請求項３〜６のいずれか一つに記載のＸ線診断装置。
【請求項８】
前記カラー画像生成部は、前記流入時間の変化に対する前記カラーの周期的な変化が規定された周期的カラーコードを用いる場合に、フレームの数、及び、前記周期的カラーコードにおけるフレーム間の時間変化量のうち、少なくとも一つが略同一である前記処理条件を用いて、時系列順に並んだ複数フレームの前記第１のカラー画像データ及び時系列順に並んだ複数フレームの前記第２のカラー画像データを生成し、
前記表示制御部は、同一時相のフレームの前記第１のカラー画像データ及び前記第２のカラー画像データに基づく立体視画像を時系列順に表示させる、
請求項３〜７のいずれか一つに記載のＸ線診断装置。
【請求項９】
前記表示制御部は、非同期再生の要求を受け付けると、前記複数フレームの第１のカラー画像データ及び前記複数フレームの第２のカラー画像データのうちいずれか一方について、時系列順に表示されるフレームの時相を変更する、
請求項８に記載のＸ線診断装置。
【請求項１０】
前記カラー画像生成部は、前記処理条件を変更する指示を受け付けた場合に、変更された処理条件を用いて、前記第１のカラー画像データ及び前記第２のカラー画像データを再生成する、
請求項１〜９のいずれか一つに記載のＸ線診断装置。
【請求項１１】
前記差分画像生成部は、予め設定された視差角度及び視差の角度方向のうち少なくとも一方を用いて前記第１の方向及び前記第２の方向を導出し、導出した前記第１の方向及び前記第２の方向のそれぞれで、前記第１の差分画像データ及び前記第２の差分画像データをそれぞれ生成する、
請求項１〜１０のいずれか一つに記載のＸ線診断装置。
【請求項１２】
前記差分画像生成部は、前記Ｘ線の照射を指示するためのＸ線トリガーボタンの半押しを受け付けた場合に、前記第１の方向から前記第２の方向へのＸ線照射方向の変更を実行する、
請求項１〜１１のいずれか一つに記載のＸ線診断装置。
【請求項１３】
前記差分画像生成部は、前記Ｘ線の照射方向の変更を指示するための装置駆動ボタンの押下を受け付けた場合に、前記第１の方向から前記第２の方向へのＸ線の照射方向の変更を実行する、
請求項１〜１１のいずれか一つに記載のＸ線診断装置。
【請求項１４】
注目領域を指定する操作を受け付けた場合に、撮像系機器に対する前記注目領域の位置を同定し、同定した前記注目領域の位置を撮像するように前記撮像系機器を制御する注目領域同定部を更に備える、
請求項１〜１３のいずれか一つに記載のＸ線診断装置。
【請求項１５】
前記差分画像生成部は、前記第１の差分画像データ及び前記第２の差分画像データのそれぞれを、２つの焦点を切り替えながらＸ線を照射するＸ線管球を用いて同一の期間に生成し、
前記カラー画像生成部は、前記２つの焦点を切り替える時間差を補正した上で、前記第１のカラー画像データ及び前記第２のカラー画像データを生成する、
請求項１〜９のいずれか一つに記載のＸ線診断装置。
【請求項１６】
前記撮像条件及び前記処理条件のうち少なくとも一つを、前記第１の差分画像データ及び前記第２の差分画像データの撮像、並びに、前記第１の差分画像データ及び前記第２の差分画像データを用いた処理において一致させることが予め設定された撮像設定情報を有する、
請求項１〜９のいずれか一つに記載のＸ線診断装置。
【請求項１７】
前記撮像設定情報を記憶する記憶回路と、
前記記憶回路に記憶された前記撮像設定情報を読み出して、読み出した撮像設定情報に基づいて、前記差分画像生成部及び前記カラー画像生成部の処理を実行する処理回路と
を更に備え、
前記処理回路は、
前記撮像設定情報において前記撮像条件及び前記処理条件を設定する操作を操作者から受け付けて、
受け付けた前記撮像条件を用いて、前記第１の方向からの前記第１の差分画像データの生成と、前記第２の方向からの前記第２の差分画像データの生成とを実行し、
受け付けた前記処理条件を用いて、前記第１のカラー画像データの生成と、前記第２のカラー画像データの生成とを実行する、
請求項１６に記載のＸ線診断装置。
【請求項１８】
第１の方向から所定の撮像条件で時系列に撮像した造影画像データと非造影画像データを差分して第１の差分画像データを生成し、第２の方向から前記所定の撮像条件と略同一の撮像条件で時系列に撮像した造影画像データと非造影画像データを差分して第２の差分画像データを生成する差分画像生成部と、
前記第１の差分画像データの各位置に、画素値の経時的変化に応じたカラーを割り当てた第１のカラー画像データを所定の処理条件で生成するとともに、前記第２の差分画像データの各位置に、画素値の経時的変化に応じたカラーを割り当てた第２のカラー画像データを前記所定の処理条件と略同一の処理条件で生成するカラー画像生成部と、
前記第１のカラー画像データ及び前記第２のカラー画像データに基づく立体視画像を表示させる表示制御部と
を備える、画像処理装置。
【請求項１９】
第１の方向から所定の撮像条件で時系列に撮像した造影画像データと非造影画像データを差分して第１の差分画像データを生成し、第２の方向から前記所定の撮像条件と略同一の撮像条件で時系列に撮像した造影画像データと非造影画像データを差分して第２の差分画像データを生成する差分画像生成ステップと、
前記第１の差分画像データの各位置に、画素値の経時的変化に応じたカラーを割り当てた第１のカラー画像データを所定の処理条件で生成するとともに、前記第２の差分画像データの各位置に、画素値の経時的変化に応じたカラーを割り当てた第２のカラー画像データを前記所定の処理条件と略同一の処理条件で生成するカラー画像生成ステップと、
前記第１のカラー画像データ及び前記第２のカラー画像データに基づく立体視画像を表示させる表示ステップと
を含む、画像診断支援方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明の実施形態は、Ｘ線診断装置、画像処理装置、及び画像診断支援方法に関する。
【背景技術】
【０００２】
従来、Ｘ線診断装置における血管撮像法の１つとして、ＤＳＡ（Digital  Subtraction  Angiography）が知られている。ＤＳＡは、被検体に対する造影剤の注入前後におけるＸ線画像データの差分を行うことで、造影剤によって染影された血管が選択的に描出された画像データを得る技術である。例えば、ＤＳＡでは、造影剤を注入する前に撮像することにより、造影剤非存在下におけるＸ線画像データをマスク画像データとして収集する。そして、造影剤を注入しながら撮像することにより、造影剤存在下におけるＸ線画像データをコントラスト画像データとして収集する。そして、マスク画像データとコントラスト画像データとの間における差分処理により、ＤＳＡ画像データが生成される。
【０００３】
また、上記のＤＳＡを利用して、造影剤の流入時間に関するパラメータを画像化するパラメトリックイメージングと呼ばれる技術も存在する。パラメトリックイメージングでは、例えば、ＤＳＡ画像データの各位置における画素値の変化を造影剤濃度の変化と見なし、画素値の時系列変化がピークあるいは特定の値となる時間を流入時間として算出する。そして、パラメトリックイメージングでは、算出した流入時間に応じたカラーを各位置にマッピングすることにより、パラメトリックイメージング画像データ（「パラメトリック画像データ」とも表記）を生成する。
【０００４】
また、Ｘ線診断装置では、立体視可能なＸ線画像データを提供する技術も各種提案されている。例えば、Ｃ型アームの角度を変更して右目用と左目用の視差画像を撮像する技術が知られている。
【先行技術文献】
【特許文献】
【０００５】
【特許文献１】特開２０１４−２００３３９号公報
【特許文献２】特開２０１４−２００５９３号公報
【発明の概要】
【発明が解決しようとする課題】
【０００６】
本発明が解決しようとする課題は、パラメトリックイメージングにおいて立体視可能な画像を提供することができるＸ線診断装置、画像処理装置、及び画像診断支援方法を提供することである。
【課題を解決するための手段】
【０００７】
実施形態のＸ線診断装置は、差分画像生成部と、カラー画像生成部と、表示制御部とを備える。差分画像生成部は、第１の方向から照射されるＸ線を用いて、所定の撮像条件で時系列に撮像した造影画像データと非造影画像データを差分して第１の差分画像データを生成し、第２の方向から照射されるＸ線を用いて、前記所定の撮像条件と略同一の撮像条件で時系列に撮像した造影画像データと非造影画像データを差分して第２の差分画像データを生成する。カラー画像生成部は、前記第１の差分画像データの各位置に、画素値の経時的変化に応じたカラーを割り当てた第１のカラー画像データを所定の処理条件で生成するとともに、前記第２の差分画像データの各位置に、画素値の経時的変化に応じたカラーを割り当てた第２のカラー画像データを前記所定の処理条件と略同一の処理条件で生成する。表示制御部は、前記第１のカラー画像データ及び前記第２のカラー画像データに基づく立体視画像を表示させる。
【図面の簡単な説明】
【０００８】
【図１】図１は、第１の実施形態に係るＸ線診断装置の構成の一例を示す図である。
【図２Ａ】図２Ａは、第１の実施形態に係るＣ型アームの回転方向について説明するための図である。
【図２Ｂ】図２Ｂは、第１の実施形態に係るＣ型アームの回転方向について説明するための図である。
【図３】図３は、第１の実施形態に係るＸ線診断装置の処理手順を示すフローチャートである。
【図４】図４は、第１の実施形態に係るＤＳＡ画像生成処理の処理手順を示すフローチャートである。
【図５】図５は、第１の実施形態に係るパラメトリック画像生成処理の処理手順を示すフローチャートである。
【図６】図６は、第１の実施形態に係るパラメトリック画像生成機能の処理を説明するための図である。
【図７Ａ】図７Ａは、第１の実施形態に係る静止画生成におけるカラーコード生成処理を説明するための図である。
【図７Ｂ】図７Ｂは、第１の実施形態に係る静止画生成におけるカラーコード生成処理を説明するための図である。
【図８】図８は、第１の実施形態に係る動画生成におけるカラーコード生成処理を説明するための図である。
【図９】図９は、第１の実施形態に係るパラメトリック画像生成機能の処理を説明するための図である。
【図１０】図１０は、第１の実施形態に係る表示制御機能による同期表示について説明するための図である。
【図１１】図１１は、第１の実施形態に係る非同期表示モードによるフレーム番号の補正処理について説明するための図である。
【図１２】図１２は、第１の実施形態に係るＸ線診断装置の効果を説明するための図である。
【図１３】図１３は、第２の実施形態に係るステレオＸ線管球を用いた撮像について説明するための図である。
【図１４】図１４は、第３の実施形態に係るＸ線診断装置の構成の一例を示す図である。
【図１５】図１５は、第３の実施形態に係る注目領域同定機能の処理を説明するための図である。
【図１６】図１６は、第３の実施形態に係る注目領域同定機能の処理を説明するための図である。
【図１７】図１７は、第３の実施形態に係る注目領域同定機能の処理を説明するための図である。
【図１８】図１８は、第３の実施形態に係る注目領域同定機能の処理を説明するための図である。
【図１９】図１９は、第３の実施形態に係る注目領域同定機能の処理を説明するための図である。
【発明を実施するための形態】
【０００９】
以下、図面を参照して、実施形態に係るＸ線診断装置、画像処理装置、及び画像診断支援方法を説明する。なお、実施形態は、以下の実施形態に限られるものではない。また、一つの実施形態に記載した内容は、原則として他の実施形態にも同様に適用される。
【００１０】
（第１の実施形態）
図１は、第１の実施形態に係るＸ線診断装置１の構成の一例を示す図である。図１に示すように、第１の実施形態に係るＸ線診断装置１は、Ｘ線撮像機構１０と、画像処理装置１００とを有する。
【００１１】
Ｘ線撮像機構１０は、Ｘ線管球１１と、検出器（Flat  Panel  Detector：ＦＰＤ）１２と、Ｃ型アーム１３と、寝台１４とを有し、インジェクター６０が接続される。
【００１２】
インジェクター６０は、被検体Ｐに挿入されたカテーテルから造影剤を注入するための装置である。ここで、インジェクター６０からの造影剤注入開始は、後述する画像処理装置１００を介して受信した注入開始指示に従って実行される場合であってもよいし、術者などの操作者が直接インジェクター６０に対して入力した注入開始指示に従って実行される場合であってもよい。
【００１３】
Ｃ型アーム１３は、Ｘ線管球１１と、Ｘ線管球１１から照射されたＸ線を検出する検出器１２とを支持する。Ｃ型アーム１３は、図示しないモータにより、寝台１４上に横臥する被検体Ｐの周りをプロペラのように高速回転する。ここで、Ｃ型アーム１３は、直交する３軸であるＸＹＺ軸に関してそれぞれ回転可能に支持され、図示しない駆動部によって各軸で個別に回転する。なお、Ｃ型アーム１３は、支持機の一例である。
【００１４】
図２Ａ及び図２Ｂは、第１の実施形態に係るＣ型アーム１３の回転方向について説明するための図である。図２Ａには、Ｙ軸の正方向から被検体Ｐを見た図を例示し、図２Ｂには、Ｘ軸の正方向から被検体Ｐを見た図を例示する。
【００１５】
図２Ａに示すように、被検体Ｐの正面を示すＺ軸を基準とし、被検体Ｐの左方向への回転は、ＬＡＯ（Left  Anterior  Oblique）と呼ばれ、ＬＡＯの反対の回転は、ＲＡＯ（Right  Anterior  Oblique）と呼ばれる。また、図２Ｂに示すように、被検体Ｐの正面を示すＺ軸を基準とし、被検体Ｐの頭頂方向への回転は、ＣＲＡ（CRANIAL）と呼ばれ、ＣＲＡの反対の回転は、ＣＡＵ（CAUDAL）と呼ばれる。例えば、Ｘ線管球１１がＺ軸の正方向（つまり、患者正面）に位置する場合、Ｃ型アーム１３の角度は、（ＲＡＯ／ＬＡＯ  ０，ＣＲＡ／ＣＡＵ  ０）と表される。これは、ＲＡＯ／ＬＡＯ角が０度であり、ＣＲＡ／ＣＡＵ角が０度であることを表す。
【００１６】
Ｘ線管球１１は、図示しない高電圧発生器から供給される高電圧を用いてＸ線を発生するＸ線源である。検出器１２は、被検体Ｐを透過したＸ線を検出するための複数のＸ線検出素子がマトリックス状に配列された装置である。この検出器１２が有する各Ｘ線検出素子は、被検体Ｐを透過したＸ線を後述するＡ／Ｄ変換器２１に出力する。
【００１７】
画像処理装置１００は、図１に示すように、Ａ／Ｄ（Analog／Digital）変換器２１と、画像メモリ２２と、サブトラクション回路２３と、フィルタリング回路２４と、アフィン変換回路２５と、ＬＵＴ（Look  Up  Table）２６と、撮像制御回路２７と、３次元再構成回路３１と、３次元画像処理回路３２と、処理回路３３と、ディスプレイ４０と、入力インターフェース５０を有する。
【００１８】
ディスプレイ４０は、画像処理装置１００によって処理された各種画像や、ＧＵＩ（Graphical  User  Interface）などの各種情報を表示する。例えば、ディスプレイ４０は、ＣＲＴ（Cathode  Ray  Tube）モニタや液晶モニタなどである。
【００１９】
ここで、ディスプレイ４０は、左目用の画像データと右目用の画像データとに基づいて立体視することが可能な立体視画像を表示させることができる立体視専用のディスプレイである。例えば、ディスプレイ４０は、表示面にかまぼこのような形状のレンズを幾つも並べたようなレンチキュラーシートや、ハエの目のような多数のレンズからなるハエの目レンズを張り付けた構造をしており、立体視専用メガネを使用しなくてもレンズで光の軌跡を変更することにより裸眼で立体画像を観察することができる。なお、ディスプレイ４０は、裸眼立体視専用のディスプレイでなくてもよい。かかる場合、ディスプレイ４０は、立体視専用メガネと同期するディスプレイであり、左目用の画像データを表示している時間はメガネの左側だけ透過し、右側は不透過となる。逆に右目用の画像データを表示している時間はメガネの右側だけ透過し、左側は不透過となる。あるいは、ディスプレイ４０は、表示面に偏光フィルタを張り付けた構造をしており、例えば偶数画素ラインには横偏光、奇数画素ラインには縦偏光を施す。立体視専用メガネの左目側では横偏光の光のみを、右目側では縦偏光の光のみを透過させるようになっており、偶数画素ラインに左目用の画像データ、奇数画素ラインに右目用の画像データを表示させる。このように立体視専用のメガネを用いることで、立体視可能なＸ線画像データを表示させる。
【００２０】
入力インターフェース５０は、例えば、マウス、キーボード、ボタン、パネルスイッチ、タッチコマンドスクリーン、フットスイッチ、トラックボール、ジョイスティック等の入力装置に対応する。入力インターフェース５０は、操作者からの各種指示を受け付け、受け付けた各種指示を画像処理装置１００の各回路に対して適宜転送する。
【００２１】
また、例えば、入力インターフェース５０には、Ｘ線の照射を指示するためのＸ線トリガーボタンが含まれる。Ｘ線トリガーボタンが操作者により押下されると、Ｘ線診断装置１は、Ｘ線画像データの撮像を開始する。また、例えば、入力インターフェース５０には、Ｘ線の照射方向の変更を指示するための装置駆動ボタンが含まれる。装置駆動ボタンが操作者により押下されると、Ｘ線診断装置１は、予め設定された方向にＣ型アーム１３を回転させることでＸ線の照射方向を変更する。
【００２２】
Ａ／Ｄ変換器２１は、検出器１２に接続され、検出器１２から入力されたアナログ信号をデジタル信号に変換し、変換したデジタル信号をＸ線画像データとして画像メモリ２２に格納する。
【００２３】
画像メモリ２２は、Ｘ線画像データを記憶する。また、画像メモリ２２は、後述する３次元再構成回路３１によって再構成された再構成データ（ボリュームデータ）や、３次元画像処理回路３２によって生成された３次元画像を記憶する。なお、画像メモリ２２は、コンピュータによって実行可能なプログラムを記憶可能である。
【００２４】
サブトラクション回路２３は、ＤＳＡ（Digital  Subtraction  Angiography）画像データなどの差分画像データを生成する。例えば、サブトラクション回路２３は、画像メモリ２２に記憶されたマスク画像データ及びコントラスト画像データを用いてＤＳＡ画像データを、又は２つのボリュームデータを用いて血管構造を有するボリュームデータを生成する。ここで、マスク画像データは、造影剤を注入される前に撮像されたＸ線画像データ（非造影画像データ）に対応する。また、コントラスト画像データは、造影剤を注入しながら撮像されたＸ線画像データ（造影画像データ）に対応する。
【００２５】
フィルタリング回路２４は、高周波強調フィルタリングなどを行う。アフィン変換回路２５は、画像の拡大や縮小、移動などを行う。ＬＵＴ２６は、諧調変換を行うためのテーブルを記憶する。
【００２６】
撮像制御回路２７は、後述する処理回路３３の制御のもと、Ｘ線撮像機構１０による撮像に係る各種処理を制御する。例えば、撮像制御回路２７は、Ｃ型アーム１３を回転させながら所定のフレームレートでＸ線画像を撮像する回転撮像を制御する。一例を挙げると、撮像制御回路２７は、インジェクター６０から造影剤注入開始時に出力される信号を契機として、単一の造影剤注入の後に複数回のＸ線画像データの回転撮像を制御する。ここで、撮像制御回路２７は、単一の造影剤の注入開始時刻を起点とした経過時間により複数回の回転撮像のスタートを制御することで、各回転撮像の対象に造影剤が到達するタイミングに合わせた回転撮像を行う。
【００２７】
また、撮像制御回路２７は、Ｃ型アーム１３を回転制御している間、図示しない高電圧発生器を制御してＸ線管球１１からＸ線を連続的又は断続的に発生させ、検出器１２によって被検体Ｐを透過したＸ線を検出させるように制御する。ここで、撮像制御回路２７は、後述する処理回路３３によって回転撮像ごとに設定されるＸ線の発生条件に基づいて、Ｘ線管球１１からＸ線を発生させる。
【００２８】
３次元再構成回路３１は、Ｘ線撮像機構１０による回転撮像によって収集されたＸ線画像から再構成データ（ボリュームデータ）を再構成する。例えば、３次元再構成回路３１は、サブトラクション回路２３によってマスク画像データとしての回転Ｘ線画像データとコントラスト画像としての回転Ｘ線画像の中でマスク画像データと角度が略一致するＸ線画像とが差分され、画像メモリ２２によって記憶されたサブトラクション後の投影データから血管構造を有するボリュームデータを再構成する。或いは、３次元再構成回路３１は、画像メモリ２２に記憶されたマスク画像データとしての回転Ｘ線画像データとコントラスト画像データとしての回転Ｘ線画像データとを用いて別々にボリュームデータを再構成し、２つのボリュームデータをサブトラクションすることで血管構造を有するボリュームデータを生成する。そして、３次元再構成回路３１は、再構成したボリュームデータを画像メモリ２２に格納する。
【００２９】
３次元画像処理回路３２は、画像メモリ２２によって記憶されたボリュームデータから３次元医用画像データを生成する。例えば、３次元画像処理回路３２は、ボリュームデータからボリュームレンダリング画像データや、ＭＰＲ（Multi  Planar  Reconstruction）画像データを生成する。そして、３次元画像処理回路３２は、生成した３次元医用画像データを画像メモリ２２に格納する。また、３次元画像処理回路３２は、ＬＵＴ２６を参照して、３次元医用画像データの階調変換を行う。
【００３０】
処理回路３３は、Ｘ線診断装置１全体を制御する。具体的には、処理回路３３は、Ｘ線撮像機構１０によるＸ線画像データの撮像、表示画像の生成、ディスプレイ４０における表示画像の表示などに係る各種処理を制御する。例えば、処理回路３３は、Ｘ線撮像機構１０による回転撮像や、回転撮像によって撮像されたＸ線画像データから３次元画像データを生成してディスプレイ４０に表示させる。
【００３１】
また、処理回路３３は、図１に示すように、パラメトリック画像生成機能３３２と、表示制御機能３３３とを実行する。ここで、例えば、図１に示す処理回路３３の構成要素であるパラメトリック画像生成機能３３２と、表示制御機能３３３とが実行する各処理機能は、コンピュータによって実行可能なプログラムの形態でＸ線診断装置１の記憶装置（例えば、画像メモリ２２）に記録されている。処理回路３３は、各プログラムを記憶装置から読み出し、実行することで各プログラムに対応する機能を実現するプロセッサである。換言すると、各プログラムを読み出した状態の処理回路３３は、図１の処理回路３３内に示された各機能を有することとなる。
【００３２】
なお、図１に図示した内容は、一例に過ぎない。例えば、図１には、サブトラクション回路２３、フィルタリング回路２４、アフィン変換回路２５、撮像制御回路２７、３次元再構成回路３１、３次元画像処理回路３２、及び処理回路３３の複数の回路（プロセッサ）を例示したが、これらの回路は必ずしも独立して構成されなくともよい。例えば、これらの回路のうち任意の回路を適宜組み合わせて構成されてもよい。
【００３３】
上記説明において用いた「プロセッサ」という文言は、例えば、ＣＰＵ（Central  Processing  Unit）、ＧＰＵ（Graphics  Processing  Unit）、或いは、特定用途向け集積回路（Application  Specific  Integrated  Circuit：ＡＳＩＣ）、プログラマブル論理デバイス（例えば、単純プログラマブル論理デバイス（Simple  Programmable  Logic  Device：ＳＰＬＤ）、複合プログラマブル論理デバイス（Complex  Programmable  Logic  Device：ＣＰＬＤ）、及びフィールドプログラマブルゲートアレイ（Field  Programmable  Gate  Array：ＦＰＧＡ））等の回路を意味する。プロセッサは記憶回路に保存されたプログラムを読み出し実行することで機能を実現する。なお、記憶回路にプログラムを保存する代わりに、プロセッサの回路内にプログラムを直接組み込むよう構成しても構わない。この場合、プロセッサは回路内に組み込まれたプログラムを読み出し実行することで機能を実現する。なお、本実施形態の各プロセッサは、プロセッサごとに単一の回路として構成される場合に限らず、複数の独立した回路を組み合わせて１つのプロセッサとして構成し、その機能を実現するようにしてもよい。更に、各図における複数の構成要素を１つのプロセッサへ統合してその機能を実現するようにしてもよい。
【００３４】
以上、第１の実施形態に係るＸ線診断装置１の構成の一例を説明した。このような構成を有するＸ線診断装置１は、パラメトリックイメージングにおいて立体視可能な画像を提供するために、処理回路３３は、ＤＳＡ画像生成機能３３１と、パラメトリック画像生成機能３３２と、表示制御機能３３３とを実行する。
【００３５】
すなわち、ＤＳＡ画像生成機能３３１は、視点が異なる２つのＤＳＡ画像データを、互いに略同一の撮像条件で時系列に沿って生成する。なお、ＤＳＡ画像生成機能３３１は、差分画像生成部の一例である。また、ＤＳＡ画像データは、差分画像データの一例である。
【００３６】
続いて、パラメトリック画像生成機能３３２は、生成された２つのＤＳＡ画像データを用いて、視点が異なる２つのパラメトリック画像データを、互いに略同一の処理条件で生成する。ここで、パラメトリック画像データは、造影剤の流入時間に関するパラメータを画像化したものである。例えば、パラメトリック画像データは、ＤＳＡ画像データの各位置に、画素値の経時的変化に応じたカラーを割り当てたカラー画像データである。パラメトリック画像生成機能３３２は、カラー画像生成部の一例である。
【００３７】
そして、表示制御機能３３３は、生成された２つのパラメトリック画像データに基づく立体視画像を表示させる。これにより、Ｘ線診断装置１は、パラメトリックイメージングにおいてステレオ観察可能な画像を提供することができる。
【００３８】
なお、処理回路３３が実行する各機能は、例えば、Ｘ線パラメトリックイメージング用撮像プログラムとして予め設定される。例えば、Ｘ線パラメトリックイメージング用撮像プログラムには、視差角度が予め登録されている。視差角度は、ステレオ観察を行うために設定される２つの視点の位置を決定するための角度である。厳密には、視差角度は、ステレオ観察を行う観察者の左右の眼球間の距離に依存して変化するが、一般的には３〜５度が適切と考えられている。本実施形態では、視差角度として「５度」が登録されている場合を説明するが、これに限らず、任意の角度が視差角度として設定可能である。また、視差角度は、予め設定される場合に限らず、例えば、撮像を行うごとに任意に設定されても良い。
【００３９】
また、Ｘ線パラメトリックイメージング用撮像プログラムには、例えば、造影画像データ及び非造影画像データを撮像するための撮像条件、及び、パラメトリック画像データを生成するための処理条件が予め登録されている。なお、撮像条件及び処理条件については、後述する。
【００４０】
図３は、第１の実施形態に係るＸ線診断装置１の処理手順を示すフローチャートである。図３に示す処理手順は、例えば、Ｘ線パラメトリックイメージング用撮像プログラムの処理を開始する旨の指示が操作者によって入力されることにより、開始される。
【００４１】
図３に示すように、処理回路３３は、撮像が開始されたか否かを判定する（ステップＳ１０１）。例えば、処理回路３３は、Ｘ線パラメトリックイメージング用撮像プログラムの処理を開始する旨の指示が操作者によって入力されると、撮像が開始されたと判定し、ステップＳ１０２以降の処理を実行する。また、このとき、処理回路３３は、インジェクター６０と通信を行って、造影剤注入の準備を行う。なお、ステップＳ１０１が否定される場合には、処理回路３３は、撮像を開始せず、以下の処理は待機状態である。
【００４２】
続いて、処理回路３３は、観察方向の変更を受け付ける（ステップＳ１０２）。ここで、観察方向とは、１回目の撮像と２回目の撮像との間の視差の角度方向に対応する。例えば、観察方向としては、「患者左右方向」や「患者上下方向」が設定される。「患者左右方向」は、例えば、１回目の撮像から２回目の撮像へ移行する際に、Ｘ線管球１１の位置が患者（被検体Ｐ）の左右方向（ＲＡＯ／ＬＡＯ）に移動するようにＣ型アーム１３を回転させることを表す。また、「患者上下方向」は、例えば、１回目の撮像から２回目の撮像へ移行する際に、Ｘ線管球１１の位置が患者の上下方向（ＣＲＡ／ＣＡＵ）に移動するようにＣ型アーム１３を回転させることを表す。
【００４３】
例えば、Ｘ線パラメトリックイメージング用撮像プログラムには、デフォルトの観察方向として「患者左右方向」が設定されている。そして、処理回路３３は、Ｘ線パラメトリックイメージング用撮像プログラムの処理を開始する旨の指示が操作者によって入力されてから数秒の間に観察方向の変更を受け付ける。この数秒（５秒〜１０秒程度）の間に、操作者がマウスやキーボードなどの入力装置を介して観察方向を「患者上下方向」に変更する旨の指示を入力すると、処理回路３３は、デフォルトの観察方向を操作者により指示された観察方向に変更する。また、数秒の間に操作者が観察方向を変更する旨の指示を入力しない場合には、処理回路３３は、デフォルトの観察方向を用いて以下の処理を行う。なお、ここでは観察方向として「患者左右方向」若しくは「患者上下方向」が設定される場合を説明するが、これに限らず、任意の方向が観察方向として設定可能である。
【００４４】
なお、観察方向（視差の角度方向）は、血管の走行方向に応じて決定されるのが好ましい。例えば、観察対象である血管が患者の上下方向に走行している場合には、観察方向は「患者左右方向」に設定されるのが好ましい。また、観察対象である血管が患者の左右方向に走行している場合には、観察方向は「患者上下方向」に設定されるのが好ましい。これは、観察方向と血管の走行方向が略平行である場合には、ステレオ観察用の立体視画像を表示しても立体に見え難い、或いはぼやけた見え方になってしまう場合があるからである。このため、観察方向は、血管の走行方向に対して平行ではなく、直交する方向に近いほど好ましい。なお、処理回路３３は、血管の走行方向に応じて観察方向を自動的に決定しても良いが、この処理回路３３の処理については後述する。
【００４５】
ステップＳ１０３において、ＤＳＡ画像生成機能３３１は、ＤＳＡ画像生成処理を実行する。ステップＳ１０３のＤＳＡ画像生成処理は、１回目の撮像におけるＤＳＡ画像データ（以下、「第１のＤＳＡ画像データ」とも表記する）を生成するための処理である。
【００４６】
例えば、ＤＳＡ画像生成機能３３１は、Ｘ線が照射される位置（照射位置）にＸ線管球１１を移動させ、ＤＳＡ画像生成処理を実行する。例えば、操作者は、Ｘ線の照射位置を患者側面（ＬＡＯ  ９０，ＣＲＡ／ＣＡＵ  ０）に指定する操作を行う。この操作により、ＤＳＡ画像生成機能３３１は、Ｃ型アーム１３を回転させてＸ線管球１１を患者側面（ＬＡＯ  ９０，ＣＲＡ／ＣＡＵ  ０）の位置に移動させ、ＤＳＡ画像生成処理を開始する。なお、ここでは、１回目の撮像におけるＸ線の照射位置が操作者により指定される場合を説明したが、これに限定されるものではない。例えば、Ｘ線の照射位置は、デフォルトの照射位置が予め設定される場合であってもよいし、撮像部位、視差角度、及び視差の角度方向などに基づいて決定される場合であってもよい。
【００４７】
図４は、第１の実施形態に係るＤＳＡ画像生成処理の処理手順を示すフローチャートである。なお、図４に示すＤＳＡ画像生成処理の処理手順は、図３に示したステップＳ１０３の処理に対応する。
【００４８】
図４に示すように、ＤＳＡ画像生成機能３３１は、Ｘ線トリガーボタンの押下を受け付ける（ステップＳ２０１）。例えば、ＤＳＡ画像生成機能３３１は、Ｘ線トリガーボタンが操作者により押下されると（ステップＳ２０１肯定）、Ｘ線管球１１にＸ線を照射させ、ステップＳ２０２以降の処理を開始する。なお、本実施形態では、ステップＳ２０１におけるＸ線トリガーボタンの押下が継続している間（つまり、長押し状態の間）に、Ｘ線の撮像が行われる場合を説明する。言い換えると、操作者は、Ｘ線トリガーボタンを長押しする時間を任意に変更することで、撮像期間及び撮像フレーム数を調整する。なお、Ｘ線トリガーボタンが押下されるまで（ステップＳ２０１否定）、ＤＳＡ画像生成機能３３１は、ステップＳ２０２以降の処理を開始しない。
【００４９】
なお、Ｘ線トリガーボタンの押下により開始される造影画像データ及び非造影画像データの撮像条件は、上述したＸ線パラメトリックイメージング用撮像プログラムに予め登録されている。ここで登録されている撮像条件としては、例えば、適切な画像レベルを決定するターゲット入射線量（所望の画像ノイズレベルを設定）、適切なＸ線収集条件を得るためのターゲット管電圧、管電圧の許容範囲、パルス幅の許容範囲、Ｘ線焦点サイズの選択肢など、従来の撮像プログラムに予め登録されている撮像条件が含まれる。
【００５０】
続いて、ＤＳＡ画像生成機能３３１は、非造影画像データを撮像する（ステップＳ２０２）。例えば、ＤＳＡ画像生成機能３３１は、予め決められた枚数（１枚、若しくは複数枚）の非造影画像データを撮像する。ＤＳＡ画像生成機能３３１は、撮像した非造影画像データを画像メモリ２２に格納する。
【００５１】
続いて、ＤＳＡ画像生成機能３３１は、造影剤の注入を開始させる（ステップＳ２０３）。例えば、ＤＳＡ画像生成機能３３１は、インジェクター６０に対して造影剤の注入を開始する旨の注入開始信号を送信する。インジェクター６０は、ＤＳＡ画像生成機能３３１から送信された注入開始信号を受信すると、造影剤の注入を開始する。
【００５２】
続いて、ＤＳＡ画像生成機能３３１は、造影画像データを撮像する（ステップＳ２０４）。例えば、ＤＳＡ画像生成機能３３１は、造影剤の注入が開始されてから予め設定された時間が経過した時点で、造影画像データの撮像を開始する。そして、ＤＳＡ画像生成機能３３１は、Ｘ線トリガーボタンがリリース（解放）されるまで（ステップＳ２０５否定）、造影画像データを時系列に沿って連続的に撮像する。ＤＳＡ画像生成機能３３１は、造影画像データを撮像するごとに、撮像した造影画像データを画像メモリ２２に格納する。Ｘ線トリガーボタンのリリース（解放）を受け付けると（ステップＳ２０５肯定）、ＤＳＡ画像生成機能３３１は、造影画像データの撮像を停止する（ステップＳ２０６）。
【００５３】
そして、ＤＳＡ画像生成機能３３１は、ＤＳＡ画像データを生成する（ステップＳ２０７）。例えば、ＤＳＡ画像生成機能３３１は、時系列に沿って連続的に撮像された複数フレームの造影画像データそれぞれから、非造影画像データをサブトラクションすることで、時系列に沿った複数フレームのＤＳＡ画像データを生成する。なお、サブトラクションに用いられる非造影画像データは、ノイズを低減させるため、複数枚の非造影画像データが加算平均（若しくはボトムトレース）された画像データであるのが好適であるが、任意の１枚の非造影画像データであってもよい。また、生成されたＤＳＡ画像データ（若しくは造影画像データ）は、ディスプレイ４０に適宜表示されてよい。
【００５４】
このように、ＤＳＡ画像生成機能３３１は、時系列に沿った複数フレームのＤＳＡ画像データ（第１のＤＳＡ画像データ）を生成する。そして、ＤＳＡ画像生成機能３３１は、生成した複数フレームのＤＳＡ画像データをパラメトリック画像生成機能３３２へ出力する。
【００５５】
なお、上述した図４の処理手順は、あくまで一例であり、図４の例に限定されるものではない。例えば、図４の処理手順は、必ずしも上述した順序で実行されなくてもよい。例えば、ＤＳＡ画像データを生成する処理（ステップＳ２０７）は、必ずしも撮像が停止された後（ステップＳ２０６）に実行されなくてもよい。例えば、ＤＳＡ画像生成機能３３１は、時系列に沿って造影画像データを撮像するごとに、撮像した造影画像データから非造影画像データをサブトラクションすることで、ＤＳＡ画像データを順次生成してもよい。
【００５６】
また、例えば、撮像の停止は、Ｘ線トリガーボタンのリリースに限らず、例えば、再びＸ線トリガーボタンを押下することで行われてもよい。この場合、ステップＳ２０１においてＸ線トリガーボタンが長押しされなくてもよい。また、撮像の停止は、撮像を停止するための専用のボタンを押下することで行われてもよい。
【００５７】
図３の説明に戻る。パラメトリック画像生成機能３３２は、パラメトリック画像生成処理を実行する（ステップＳ１０４）。ここで、ステップＳ１０４のパラメトリック画像生成処理は、１回目の撮像におけるパラメトリック画像データ（以下、「第１のパラメトリック画像データ」とも表記する）を生成するための処理である。
【００５８】
図５は、第１の実施形態に係るパラメトリック画像生成処理の処理手順を示すフローチャートである。なお、図５に示すパラメトリック画像生成処理の処理手順は、図３に示したステップＳ１０４の処理に対応する。
【００５９】
図５に示すように、パラメトリック画像生成機能３３２は、流入時間を同定する（ステップＳ３０１）。ここで、流入時間は、ＤＳＡ画像データの各位置における画素値の変化を造影剤濃度の変化と見なし、各画素値の時系列変化に基づいて定義されるパラメータである。流入時間の同定方法としては、ＴＴＰ（Time-to-Peak）及びＴＴＡ（Time-to-Arrival）のうち任意の同定方法が選択可能である。例えば、ＴＴＰは、画素値の時間変化が最大となる時間を流入時間として同定する方法である。また、ＴＴＡは、画素値の時間変化が所定値に到達する時間、若しくは、画素値の時間変化における最大値に対して所定割合に到達する時間を流入時間として同定する方法である。
【００６０】
図６は、第１の実施形態に係るパラメトリック画像生成機能３３２の処理を説明するための図である。図６には、ＤＳＡ画像データに含まれる任意の画素について解析された時間濃度プロファイルを例示する。図６において、横軸は時間（造影画像データの撮像時間）に対応し、縦軸は画素値に対応する。なお、図６では、流入時間の同定方法がＴＴＡである場合を説明する。
【００６１】
図６に示すように、例えば、パラメトリック画像生成機能３３２は、ＤＳＡ画像生成機能３３１によって生成されたＤＳＡ画像データに含まれる各画素について、時間濃度プロファイルを解析することで流入時間を同定する。図６の例では、パラメトリック画像生成機能３３２は、最大の画素値（Ｃｍａｘ）に対して２０％（０．２×Ｃｍａｘ）に到達する時間を流入時間として特定する。
【００６２】
なお、図６の内容はあくまで一例であり、図６の例に限定されるものではない。例えば、図６では、ＴＴＡにより流入時間を同定するための割合として「２０％」が設定される場合を説明したが、これに限らず、任意の割合（若しくは値）が設定されてもよい。また、流入時間の同定方法は、ＴＴＡに限らず、ＴＴＰであってもよい。
【００６３】
続いて、パラメトリック画像生成機能３３２は、カラーコードを生成する（ステップＳ３０２）。例えば、パラメトリック画像生成機能３３２は、静止画としてパラメトリック画像データを生成する場合には、１つのカラーコードを生成する。また、パラメトリック画像生成機能３３２は、動画としてパラメトリック画像データを生成する場合には、動画に含まれるフレーム数に対応する数のカラーコードを生成する。以下、静止画生成におけるカラーコード生成処理と、動画生成におけるカラーコード生成処理とを、順に説明する。なお、カラーコードを生成するための生成条件は、デフォルトの条件が予め設定されているが、操作者により適宜設定されてもよい。
【００６４】
図７Ａ及び図７Ｂは、第１の実施形態に係る静止画生成におけるカラーコード生成処理を説明するための図である。なお、図７Ａ及び図７Ｂでは、流入時間の同定方法がＴＴＡである場合を説明する。ＴＴＡでは、カラーコードの周期及び初期値は、デフォルトとして「Auto」が設定されている。なお、「Auto」は、ＤＳＡ画像データの収集時間（撮像期間）によって自動的に設定されることを表す。
【００６５】
図７Ａに示すように、例えば、パラメトリック画像生成機能３３２は、ＤＳＡ画像データの収集時間が０〜Ｔ秒である場合には、カラーコードの周期を「Ｔ」、カラーコードの初期値を「０」にそれぞれ設定する。このカラーコードでは、流入時間ｔが「０」のときに「赤」となり、「０」から「Ｔ／２」にかけて「赤」から「緑」に徐々に変化し、「Ｔ／２」のときに「緑」となることが規定される。また、このカラーコードでは、「Ｔ／２」から「Ｔ」にかけて「緑」から「青」に徐々に変化し、「Ｔ」のときに「青」となることが規定される。
【００６６】
また、他の例としては、図７Ｂに示すように、ＤＳＡ画像データに含まれる全画素の流入時間の最小値「Ｔｍｉｎ」をカラーコードの初期値としてもよい。このカラーコードでは、流入時間ｔが「Ｔｍｉｎ」のときに「赤」となり、「Ｔｍｉｎ」から「（Ｔ＋Ｔｍｉｎ）／２」にかけて「赤」から「緑」に徐々に変化し、「（Ｔ＋Ｔｍｉｎ）／２」のときに「緑」となることが規定される。また、このカラーコードでは、「（Ｔ＋Ｔｍｉｎ）／２」から「Ｔ」にかけて「緑」から「青」に徐々に変化し、「Ｔ」のときに「青」となることが規定される。なお、図７Ｂのカラーコードは、「０」から「Ｔｍｉｎ」にかけてはカラーを割り当てないので、流入時間の変化に対するカラーの変化が大きくなる。このため、図７Ｂのカラーコードは、流入時間の小さい変化を描出するのに適していると言える。
【００６７】
このように、パラメトリック画像生成機能３３２は、静止画用のカラーコードを生成する。なお、図７Ａ及び図７Ｂの内容はあくまで一例であり、図７Ａ及び図７Ｂの例に限定されるものではない。例えば、図７Ａ及び図７Ｂでは、流入時間の同定方法がＴＴＡである場合を説明したが、これに限らず、ＴＴＰであってもよい。また、カラーコードに規定されるカラーの順序は、図７Ａ及び図７Ｂの例に限定されるものではなく、任意に設定されてよい。
【００６８】
図８は、第１の実施形態に係る動画生成におけるカラーコード生成処理を説明するための図である。図８では、動画としてパラメトリック画像データを生成する生成方法がＣＣＣ（Circular  Color  Coding）である場合を説明する。なお、ＣＣＣにおける流入時間の同定方法は、ＴＴＡであってもＴＴＰであってもよい。
【００６９】
図８に示すように、例えば、パラメトリック画像生成機能３３２は、ＤＳＡ画像データの収集時間が０〜Ｔ秒である場合には、カラーコードの周期を「Ｌ（Ｌ＜Ｔ）」、カラーコードの初期値を「０」、カラーコードのステップを「ΔＴ」にそれぞれ設定する。そして、パラメトリック画像生成機能３３２は、動画に含まれるフレーム数に対応する数のカラーコードを生成する。例えば、動画のフレーム数が「Ｎ」であれば、パラメトリック画像生成機能３３２は、１フレーム目からＮフレーム目までのＮ個のカラーコードをそれぞれ生成する。
【００７０】
例えば、１フレーム目のカラーコードでは、流入時間ｔが「０」のときに「赤」となり、「０」から「Ｌ／３」にかけて「赤」から「緑」に徐々に変化し、「Ｌ／３」のときに「緑」となることが規定される。また、１フレーム目のカラーコードでは、「Ｌ／３」から「２Ｌ／３」にかけて「緑」から「青」に徐々に変化し、「２Ｌ／３」のときに「青」となることが規定される。また、１フレーム目のカラーコードでは、「２Ｌ／３」から「Ｌ」にかけて「青」から「赤」に徐々に変化し、「Ｌ」のときに「赤」に戻ることが規定される。このように、１フレーム目のカラーコードでは、流入時間ｔが「０」から「Ｌ」に変化する間に、「赤→緑→青→赤」と変化する。また、「Ｌ」以降については、「０」から「Ｌ」までのカラーの変化の繰り返しになる。つまり、１フレーム目のカラーコードでは、流入時間ｔが「Ｌ」から「２Ｌ」に変化する間に、「赤→緑→青→赤」と変化する。また、１フレーム目のカラーコードでは、流入時間ｔが「２Ｌ」から「３Ｌ」に変化する間に、「赤→緑→青→赤」と変化する。これ以降も同様に、１フレーム目のカラーコードでは、流入時間ｔが「Ｔ」に到達するまで、「赤→緑→青→赤」と変化する。
【００７１】
２フレーム目のカラーコードは、１フレーム目のカラーコードを「ΔＴ」ずらすことにより生成される。例えば、２フレーム目のカラーコードでは、流入時間ｔが「Δｔ」のときに「赤」となり、「Δｔ」から「Ｌ／３＋Δｔ」にかけて「赤」から「緑」に徐々に変化し、「Ｌ／３＋Δｔ」のときに「緑」となることが規定される。また、２フレーム目のカラーコードでは、「Ｌ／３＋Δｔ」から「２Ｌ／３＋Δｔ」にかけて「緑」から「青」に徐々に変化し、「２Ｌ／３＋Δｔ」のときに「青」となることが規定される。また、２フレーム目のカラーコードでは、「２Ｌ／３＋Δｔ」から「Ｌ＋Δｔ」にかけて「青」から「赤」に徐々に変化し、「Ｌ＋Δｔ」のときに「赤」に戻ることが規定される。なお、「Ｌ＋Δｔ」以降についても、１フレーム目のカラーコードと同様に、「ΔＴ」から「Ｌ＋ΔＴ」までのカラーの変化の繰り返しになる。
【００７２】
すなわち、Ｎフレーム目のカラーコードは、Ｎ−１フレーム目のカラーコードを「ΔＴ」ずらすことにより生成される。言い換えると、Ｎフレーム目のカラーコードは、１フレーム目のカラーコードを「ΔＴ×（Ｎ−１）」ずらすことにより生成される。
【００７３】
なお、図８の内容はあくまで一例であり、図８の例に限定されるものではない。例えば、カラーコードに規定されるカラーの順序は、図８の例に限定されるものではなく、任意に設定されてよい。
【００７４】
このように、パラメトリック画像生成機能３３２は、動画用のカラーコードとして、動画のフレーム数に対応する数のカラーコードを生成する。動画用のカラーコードは、流入時間の変化に対するカラーの周期的な変化が規定された周期的カラーコードであると言える。なお、パラメトリック画像生成処理における処理条件は、上述したＸ線パラメトリックイメージング用撮像プログラムに予め登録されている。ここで登録されている処理条件としては、例えば、流入時間の同定方法、カラーコードの周期、カラーコードの位相、カラーコードの初期値、及び、静止画であるか動画であるかを示す情報などが含まれる。
【００７５】
図５の説明に戻る。パラメトリック画像生成機能３３２は、パラメトリック画像データを生成する（ステップＳ３０３）。例えば、パラメトリック画像生成機能３３２は、ＤＳＡ画像データの各画素に、画素値の経時的変化に応じたカラーを割り当てることにより、パラメトリック画像データを生成する。
【００７６】
図９は、第１の実施形態に係るパラメトリック画像生成機能３３２の処理を説明するための図である。図９には、被検体Ｐの脳の血管が描出されたパラメトリック画像データを例示する。
【００７７】
図９に示すように、例えば、パラメトリック画像生成機能３３２は、ステップＳ３０２において生成したカラーコードを参照し、各画素の流入時間に対応するカラーを割り当てて、パラメトリック画像データを生成する。
【００７８】
具体的には、パラメトリック画像生成機能３３２は、静止画のパラメトリック画像データを生成する。例えば、パラメトリック画像生成機能３３２は、図７Ａ（若しくは図７Ｂ）に示した静止画用のカラーコードを参照し、各画素の流入時間に対応するカラーを割り当てて、パラメトリック画像データを生成する。
【００７９】
また、パラメトリック画像生成機能３３２は、動画のパラメトリック画像データを生成する。例えば、パラメトリック画像生成機能３３２は、図８に示した動画用のカラーコードを参照し、動画に含まれるフレーム数に対応する数のパラメトリック画像データをそれぞれ生成する。具体的には、パラメトリック画像生成機能３３２は、各画素の流入時間に対応するカラーを、１フレーム目のカラーコードに基づいて割り当てることで、１フレーム目のパラメトリック画像データを生成する。また、パラメトリック画像生成機能３３２は、各画素の流入時間に対応するカラーを、２フレーム目のカラーコードに基づいて割り当てることで、２フレーム目のパラメトリック画像データを生成する。このように、パラメトリック画像生成機能３３２は、各画素の流入時間に対応するカラーを、Ｎフレーム目のカラーコードに基づいて割り当てることで、Ｎフレーム目のパラメトリック画像データを生成する。
【００８０】
なお、図９では、流入時間に応じて割り当てられるカラーについて説明したが、更に、各画素（ピクセル）の最大の画素値に基づいて明度（明るさ）を調整するのが好適である。例えば、各画素の明度は、ＤＳＡ画像データに含まれる全画素のなかの最大の画素値「Ｄｍａｘ」に対する各画素の最大の画素値「Ｄ」の割合「Ｄ／Ｄｍａｘ」にするのが好適である。
【００８１】
このように、パラメトリック画像生成機能３３２は、静止画若しくは動画のパラメトリック画像データ（第１のパラメトリック画像データ）を生成する。そして、パラメトリック画像生成機能３３２は、生成したパラメトリック画像データを表示制御機能３３３へ出力する。
【００８２】
図３の説明に戻る。処理回路３３は、Ｘ線の照射方向を変更する（ステップＳ１０５）。例えば、ＤＳＡ画像生成機能３３１は、予め設定された視差角度及び視差の角度方向（観察方向）を用いて、２回目の撮像におけるＸ線の照射位置を導出する。ここで、例えば、視差角度は「５度」であり、視差の角度方向はステップＳ１０２の処理において設定された観察方向であり、例えば、「患者左右方向」である。
【００８３】
例えば、１回目の撮像におけるＸ線の照射位置が患者側面（ＬＡＯ  ９０，ＣＲＡ／ＣＡＵ  ０）であった場合には、２回目の撮像におけるＸ線の照射位置としては、（ＬＡＯ  ８５，ＣＲＡ／ＣＡＵ  ０）及び（ＬＡＯ  ９５，ＣＲＡ／ＣＡＵ  ０）の２通りの候補が挙げられる。この場合、ＤＳＡ画像生成機能３３１が特定のアルゴリズムに従って自動的に選択可能である。例えば、ＤＳＡ画像生成機能３３１は、（ＲＡＯ  ０）〜（ＲＡＯ  １８０）の間であれば右側（ＲＡＯの正方向）に、（ＬＡＯ  ０）〜（ＬＡＯ  １８０）の間であれば左側（ＬＡＯの正方向）に観察方向を変更する。ただし、ＤＳＡ画像生成機能３３１は、Ｃ型アーム１３が機械的に回転できない場合には、自動的に逆方向に観察方向を変更する。
【００８４】
また、Ｘ線の照射方向の変更は、例えば、Ｘ線トリガーボタンの半押しによって行われる。例えば、操作者がＸ線トリガーボタンの半押しを行っている間、ＤＳＡ画像生成機能３３１は、Ｃ型アーム１３を視差の角度方向「患者左右方向」に視差角度「５度」分回転させる。この結果、例えば、ＤＳＡ画像生成機能３３１は、Ｘ線管球１１を（ＬＡＯ  ９５，ＣＲＡ／ＣＡＵ  ０）の位置に移動させる。なお、このとき、ＤＳＡ画像生成機能３３１は、ＳＩＤ（Source  Image  Distance）、コリメータの位置、ビームフィルタの種類及び厚み、及び、補償フィルタの位置などを変化させずにＸ線管球１１を移動させる。なお、ＤＳＡ画像生成機能３３１は、Ｘ線トリガーボタンの半押しに限らず、例えば、装置駆動ボタンの押下によってＸ線管球１１を移動させてもよい。
【００８５】
なお、ここでは、撮像部位（注目領域）とＣ型アーム１３の回転中心とが略一致している場合を説明するが、実施形態はこれに限定されるものではない。例えば、撮像部位とＣ型アーム１３の回転中心とが異なる場合には、視差画像（パラメトリック画像データ）において撮像部位が描出される位置及び大きさが異なるが、例えば、視差画像上での画像処理（平行移動、拡大／縮小など）により一致させることが可能である。また、注目領域をＸ線管球１１に近づけてＸ線光学系での拡大率を上げることも可能であるが、この処理については、後述する。
【００８６】
続いて、ＤＳＡ画像生成機能３３１は、ＤＳＡ画像生成処理を実行する（ステップＳ１０６）。ステップＳ１０６のＤＳＡ画像生成処理は、２回目の撮像におけるＤＳＡ画像データ（以下、「第２のＤＳＡ画像データ」とも表記する）を生成するための処理である。
【００８７】
ここで、ＤＳＡ画像生成機能３３１は、ステップＳ１０５において変更されたＸ線の照射方向で撮像を行う点を除き、１回目のＤＳＡ画像生成処理における撮像条件と略同一の撮像条件で２回目のＤＳＡ画像生成処理を行う。すなわち、ＤＳＡ画像生成機能３３１は、造影剤の注入条件、Ｘ線の照射条件、Ｘ線が照射される視野サイズ（Field  Of  View：ＦＯＶ）、Ｘ線を照射するＸ線管球１１の焦点とＸ線の検出器と間の距離（ＳＩＤ）、寝台１４の位置、及び、Ｃ型アーム１３の回転中心（Iso-Center）が略同一である撮像条件を用いて、第１のＤＳＡ画像データ及び第２のＤＳＡ画像データを生成する。また、ＤＳＡ画像生成機能３３１は、造影剤の注入開始から造影画像データの撮像開始までの時間、造影剤の注入速度、注入剤の注入器に加える圧力、及び、注入速度に到達する前の立ち上がり曲線が略同一である注入条件を用いて、第１のＤＳＡ画像データ及び第２のＤＳＡ画像データを生成する。また、ＤＳＡ画像生成機能３３１は、Ｘ線による撮像のフレームレート、Ｘ線管球１１の管電圧、Ｘ線管球１１の管電流、Ｘ線のパルス幅、Ｘ線管球１１のコリメータの位置、Ｘ線管球１１のビームフィルタの種類及び厚み、及び、Ｘ線管球１１の補償フィルタの位置が略同一である照射条件を用いて、第１のＤＳＡ画像データ及び第２のＤＳＡ画像データを生成する。なお、上記の撮像条件は、必ずしも全ての条件が完全に一致していなくてもよい。例えば、注入剤の注入器に加える圧力が略同一であれば、注入速度に到達する前の立ち上がり曲線も略同一となる条件である。このため、注入剤の注入器に加える圧力、及び、注入速度に到達する前の立ち上がり曲線のうち、少なくとも一方が略同一であればよい。また、例えば、上記の条件は、パラメトリック画像データの生成に大きく影響を与えない範囲内で誤差を含んでいてもよい。すなわち、上記の撮像条件のうち、少なくとも一つが略同一であればよい。
【００８８】
このように、ＤＳＡ画像生成機能３３１は、ステップＳ１０５において変更されたＸ線の照射方向で撮像を行う点を除き、１回目の撮像における撮像条件と略同一の撮像条件で２回目の撮像を行う。なお、２回目の撮像においてＤＳＡ画像生成機能３３１が実行する処理は、図４に示したＤＳＡ画像生成処理と基本的に同様であるので、説明を省略する。
【００８９】
続いて、パラメトリック画像生成機能３３２は、パラメトリック画像生成処理を実行する（ステップＳ１０７）。ステップＳ１０７のパラメトリック画像生成処理は、２回目の撮像におけるパラメトリック画像データ（以下、「第２のパラメトリック画像データ」とも表記する）を生成するための処理である。
【００９０】
ここで、パラメトリック画像生成機能３３２は、１回目のパラメトリック画像生成処理における処理条件と略同一の処理条件で２回目のパラメトリック画像生成処理を行う。すなわち、パラメトリック画像生成機能３３２は、流入時間の同定方法、及び、経時的変化に応じて割り当てられるカラーが規定されたカラーコードが略同一である処理条件を用いて、第１のパラメトリック画像データ及び第２のパラメトリック画像データを生成する。また、パラメトリック画像生成機能３３２は、動画のパラメトリック画像データを生成する場合には、フレームの数、及び、動画用のカラーコード（周期的カラーコード）におけるフレーム間の時間変化量が略同一である処理条件を用いて、時系列順に並んだ複数フレームの第１のパラメトリック画像データ及び時系列順に並んだ複数フレームの第２のパラメトリック画像データを生成する。なお、上記の処理条件は、必ずしも全ての条件が完全に一致していなくてもよい。例えば、上記の処理条件は、パラメトリック画像データの生成に大きく影響を与えない範囲内で誤差を含んでいてもよい。すなわち、上記の処理条件のうち、少なくとも一つが略同一であればよい。
【００９１】
このように、パラメトリック画像生成機能３３２は、１回目のパラメトリック画像生成処理における処理条件と略同一の処理条件で２回目のパラメトリック画像生成処理を行う。なお、２回目のパラメトリック画像生成処理においてパラメトリック画像生成機能３３２が実行する処理は、図５に示したパラメトリック画像生成処理と基本的に同様であるが、同一のカラーコードを用いる場合には、カラーコードを生成する処理（ステップＳ３０２）は実行されなくてよい。また、流入時間を同定する処理（ステップＳ３０１）とパラメトリック画像データを生成する処理（ステップＳ３０３）については、同様であるので説明を省略する。
【００９２】
そして、表示制御機能３３３は、視点が異なる２つのパラメトリック画像データに基づいて、立体視画像を表示させる（ステップＳ１０８）。例えば、表示制御機能３３３は、ステップＳ１０４において生成された第１のパラメトリック画像データと、ステップＳ１０７において生成された第２のパラメトリック画像データとに基づいて、立体視画像を表示させる。
【００９３】
例えば、１回目の撮像におけるＸ線の照射位置が（ＬＡＯ  ９０，ＣＲＡ／ＣＡＵ  ０）であり、２回目の撮像におけるＸ線の照射位置が（ＬＡＯ  ９５，ＣＲＡ／ＣＡＵ  ０）である場合には、表示制御機能３３３は、第１のパラメトリック画像データを左目用とし、第２のパラメトリック画像データを右目用として、ステレオ表示を行う。具体的には、表示制御機能３３３は、第１のパラメトリック画像データ及び第２のパラメトリック画像データを、裸眼で立体視画像を観察可能なディスプレイ、若しくは立体視専用メガネと同期するディスプレイに表示させることで、ステレオ表示を行う。なお、ステレオ表示を行うディスプレイとしては、従来の如何なるステレオ表示用ディスプレイが適用されてもよい。
【００９４】
具体的には、表示制御機能３３３は、静止画の第１のパラメトリック画像データと、静止画の第２のパラメトリック画像データとに基づいて、静止画の立体視画像を表示させる。また、表示制御機能３３３は、動画の第１のパラメトリック画像データと、動画の第２のパラメトリック画像データとに基づいて、動画の立体視画像を同期表示させる。
【００９５】
図１０は、第１の実施形態に係る表示制御機能３３３による同期表示について説明するための図である。図１０には、再生時間の経過に伴って同期表示されるフレーム番号の順序を例示する。
【００９６】
図１０に示すように、動画の立体視画像を表示させる場合には、表示制御機能３３３は、同一時相のフレームの第１のパラメトリック画像データ及び第２のパラメトリック画像データに基づく立体視画像を時系列順に表示させる。具体的には、表示制御機能３３３は、１フレーム目の第１のパラメトリック画像データと、１フレーム目の第２のパラメトリック画像データとに基づいて、１フレーム目の立体視画像を表示させる。また、表示制御機能３３３は、２フレーム目の第１のパラメトリック画像データと、２フレーム目の第２のパラメトリック画像データとに基づいて、２フレーム目の立体視画像を表示させる。このように、表示制御機能３３３は、Ｎフレーム目の第１のパラメトリック画像データと、Ｎフレーム目の第２のパラメトリック画像データとに基づいて、Ｎフレーム目の立体視画像を表示させる。
【００９７】
このように、表示制御機能３３３は、視点が異なる２つのパラメトリック画像データに基づいて、立体視画像を表示させる。これによれば、Ｘ線診断装置１は、パラメトリックイメージングにおいて立体視可能な画像を提供することができる。
【００９８】
そして、処理回路３３は、補正処理を実行する（ステップＳ１０９）。例えば、処理回路３３は、パラメトリック画像データを再生成する補正処理、同一位置の流入時間を揃える補正処理、及び非同期表示モードによるフレーム番号の補正処理のうち、操作者の要求に応じて任意の補正処理を実行する。以下、各補正処理について順に説明する。
【００９９】
（パラメトリック画像データを再生成する補正処理）
パラメトリック画像データを再生成する補正処理は、パラメトリック画像生成処理の処理条件を変更する指示を受け付けた場合に実行される補正処理である。例えば、パラメトリック画像生成機能３３２は、パラメトリック画像生成処理の処理条件を変更する指示を受け付けると、変更された処理条件を用いて、第１のパラメトリック画像データ及び第２のパラメトリック画像データを再生成する。
【０１００】
例えば、処理条件の一つであるカラーコードの周期を「Ｔ」から「Ｔ／２」に変更する旨の指示を操作者から受け付けると、パラメトリック画像生成機能３３２は、変更されたカラーコードの周期「Ｔ／２」を用いて、カラーコードを再生成する。そして、パラメトリック画像生成機能３３２は、再生成されたカラーコードを用いて、第１のパラメトリック画像データ及び第２のパラメトリック画像データを再生成する。
【０１０１】
このように、パラメトリック画像生成機能３３２は、パラメトリック画像データを再生成する補正処理を実行する。
【０１０２】
（同一位置の流入時間を揃える補正処理）
同一位置の流入時間を揃える補正処理は、視点が異なる２つのパラメトリック画像データにおける同一位置を同定することで実行される補正処理である。例えば、パラメトリック画像生成機能３３２は、同一位置の流入時間を揃える補正処理を実行する旨の指示を操作者から受け付けると、第１のパラメトリック画像データ及び第２のパラメトリック画像データにおける同一位置を同定する。
【０１０３】
例えば、パラメトリック画像生成機能３３２は、２つのパラメトリック画像データの解剖学的な特徴を利用して、第１のパラメトリック画像データの各位置（画素）と、第２のパラメトリック画像データの各位置との対応づけを行う。そして、パラメトリック画像生成機能３３２は、同一位置における流入時間を比較して、流入時間を揃える補正を行う。
【０１０４】
具体的には、パラメトリック画像生成機能３３２は、第１のパラメトリック画像データ上で代表的な領域（例えば、血管の分岐を示す画素）を選択する。そして、パラメトリック画像生成機能３３２は、第１のパラメトリック画像データ上の代表的な領域と、解剖学的な構造が略一致する第２のパラメトリック画像データ上の領域を同定する。同定した領域において、１回目の撮像における流入時間と２回目の撮像における流入時間との間の誤差を算出する。そして、パラメトリック画像生成機能３３２は、算出した誤差が無くなるように、１回目の撮像又は２回目の撮像の全画素の流入時間を補正する。そして、パラメトリック画像生成機能３３２は、補正後の流入時間を用いて、流入時間が補正された側のパラメトリック画像データを再生成する。
【０１０５】
このように、パラメトリック画像生成機能３３２は、第１のＤＳＡ画像データ及び第２のＤＳＡ画像データの同一位置から求まる流入時間を揃える補正処理を行う。そして、パラメトリック画像生成機能３３２は、第１のパラメトリック画像データ及び第２のパラメトリック画像データのうち少なくとも一方を再生成する。
【０１０６】
なお、上記の説明はあくまで一例であり、上記の例に限定されるものではない。例えば、上記の説明では、代表的な領域を一つ選択する場合を説明したが、実施形態はこれに限定されるものではない。例えば、代表的な領域が複数選択される場合であってもよい。この場合、パラメトリック画像生成機能３３２は、複数の領域の流入時間の平均値を用いて誤差を算出するのが好適である。
【０１０７】
また、例えば、上記の説明では、同一領域における流入時間を用いて誤差を算出する場合を説明したが、実施形態はこれに限定されるものではない。例えば、パラメトリック画像生成機能３３２は、同一領域における時間濃度プロファイルを用いて、誤差を算出してもよい。
【０１０８】
また、例えば、上記の説明では、操作者から指示を受け付けることにより、同一位置の流入時間を揃える補正処理を実行する場合を説明したが、実施形態はこれに限定されるものではない。例えば、パラメトリック画像生成機能３３２は、操作者からの指示を受け付けなくとも、自動的に本補正処理を実行してもよい。なお、自動的に本補正処理を実行する場合には、本補正処理は、立体視画像を表示する処理（ステップＳ１０８）が実行される前に実行されてもよい。
【０１０９】
（非同期表示モードによるフレーム番号の補正処理）
非同期表示モードによるフレーム番号の補正処理は、動画としてパラメトリック画像データが生成される場合（ＣＣＣ）に、非同期表示モードにおいて実行される補正処理である。
【０１１０】
例えば、表示制御機能３３３は、非同期表示モードの要求を受け付けると、複数フレームの第１のパラメトリック画像データ及び複数フレームの第２のパラメトリック画像データのうちいずれか一方について、時系列順に表示されるフレームの時相を変更する。
【０１１１】
図１１は、第１の実施形態に係る非同期表示モードによるフレーム番号の補正処理について説明するための図である。図１１には、図１０の同期表示において非同期表示モードの要求が操作者によって行われた場合を例示する。
【０１１２】
図１１に示すように、非同期表示モードの要求を操作者から受け付けると、表示制御機能３３３は、第１のパラメトリック画像データ及び第２のパラメトリック画像データの同期表示を解除する。ここで、操作者がいずれか一方のフレーム番号を変更する指示を行うと、表示制御機能３３３は、指示された方のフレーム番号を変更する。例えば、表示制御機能３３３は、１フレーム目の第１のパラメトリック画像データと、１フレーム目の第２のパラメトリック画像データとが表示される場合に、第２のパラメトリック画像データのフレーム番号を繰り上げる旨の指示を操作者から受け付ける。この場合、表示制御機能３３３は、受け付けた指示に応じて、第２のパラメトリック画像データのフレーム番号を繰り上げる補正を行う。つまり、表示制御機能３３３は、１フレーム目の第１のパラメトリック画像データと、２フレーム目の第２のパラメトリック画像データとに基づいて、１フレーム目の立体視画像を表示させる。なお、操作者が非同期表示モードを終了する旨の指示を行うと、表示制御機能３３３は、非同期表示モードを終了し、同期表示モードに戻る。
【０１１３】
このように、表示制御機能３３３は、非同期表示モードの要求を受け付けると、複数フレームの第１のパラメトリック画像データ及び複数フレームの第２のパラメトリック画像データのうちいずれか一方について、時系列順に表示されるフレームの時相（フレーム番号）を変更する。
【０１１４】
このように、処理回路３３は、補正処理を実行する。そして、処理回路３３は、例えば、撮像を終了する旨の指示を操作者から受け付けると、図３の処理を終了する。
【０１１５】
なお、上述した図３の処理手順は、あくまで一例であり、図３の例に限定されるものではない。例えば、図３の処理手順は、必ずしも上述した順序で実行されなくてもよい。例えば、補正処理（ステップＳ１０９）は、必ずしも実行されなくてもよい。例えば、補正処理は、操作者が必要に応じて実行すればよい。また、例えば、同一位置の流入時間を揃える補正処理（ステップＳ１０９）が自動的に実行される場合には、補正処理は、立体視画像を表示する処理（ステップＳ１０８）が実行される前に実行されてもよい。さらに補正処理はフレーム単位での時相変更について説明したが、フレーム単位より細かく時間変更をしても良い。具体的にはΔｔ／５ずつ時間を変更しても良い。この場合Δｔ／５ずらした第２のパラメトリック画像データは存在しないので、パラメトリック画像データを作成し直す必要がある。
【０１１６】
また、図３の例では、専用の撮像プログラムとしてＸ線パラメトリックイメージング用撮像プログラムが用いられる場合を説明したが、実施形態はこれに限定されるものではない。例えば、従来のＤＳＡ画像撮像用プログラムを用いて２回の撮像を略同一の撮像条件で実行する場合であってもよい。この場合、撮像された２つのＤＳＡ画像データに対して、略同一の処理条件で従来のパラメトリックイメージング用プログラムを適用することにより、視差の異なる２つのパラメトリック画像データが生成される。
【０１１７】
上述してきたように、第１の実施形態に係るＸ線診断装置１は、視点が異なる２つのＤＳＡ画像データを、互いに略同一の撮像条件で時系列に沿って生成する。そして、Ｘ線診断装置１は、互いに略同一の処理条件で、生成した２つのＤＳＡ画像データから視点が異なる２つのパラメトリック画像データを生成する。そして、Ｘ線診断装置１は、生成した２つのパラメトリック画像データに基づいて、立体視画像を表示させる。これにより、Ｘ線診断装置１は、パラメトリックイメージングにおいてステレオ観察可能な画像を提供することができる。
【０１１８】
ところで、ステレオ表示を行うためには、左右の画像における同一位置が同一色で描出されていなければ立体視できないという問題がある。また、パラメトリックイメージングでは造影剤濃度の時系列変化に基づいて画像化を行うため、時系列変化を観察可能なレベルのフレームレートが要求される。すなわち、パラメトリックイメージングにおいてステレオ観察可能な画像を提供するためには、一定以上のフレームレートを維持しつつ、同一位置を同一色で描出した視差画像を生成することが要求される。そこで、第１の実施形態に係るＸ線診断装置１は、以下の処理により、一定以上のフレームレートを維持しつつ、同一位置を同一色で描出した視差画像（パラメトリック画像データ）を生成する。
【０１１９】
図１２は、第１の実施形態に係るＸ線診断装置１の効果を説明するための図である。図１２に示すように、例えば、第１の実施形態に係るＸ線診断装置１は、ステレオ表示用の左右のパラメトリック画像データを撮像するために、２回の撮像を行う。すなわち、Ｘ線診断装置１は、１回目の撮像を行ってからＸ線の照射方向を変更し、２回目の撮像を行うことで、視点が異なる２つのＤＳＡ画像データを撮像する。そして、Ｘ線診断装置１は、撮像した２つのＤＳＡ画像データにおいて造影剤の流入時間を算出し、算出した流入時間に応じたカラーを割り当てることで、視点が異なる２つのパラメトリック画像データを生成する。
【０１２０】
ここで、第１の実施形態に係るＸ線診断装置１は、２回の撮像において、それぞれの照射方向で固定した状態で時系列に沿った複数フレームのＤＳＡ画像データを撮像する。このため、Ｘ線診断装置１は、それぞれの照射方向におけるＤＳＡ画像データのフレームレートを一定以上に保つことができる。
【０１２１】
また、第１の実施形態に係るＸ線診断装置１は、２回の撮像において、Ｘ線の照射方向（Ｘ線管球１１の照射位置）が異なっている点を除き、略同一の撮像条件で２つのＤＳＡ画像データを撮像し、略同一の処理条件で２つのパラメトリック画像データを生成する。例えば、Ｘ線診断装置１は、Ｘ線の照射方向以外の照射条件や、造影剤を注入する注入条件、視野サイズ、ＳＩＤ、寝台１４の位置、及び、Iso-Centerなどの撮像条件を略同一として、視点が異なる２つのＤＳＡ画像データを撮像する。また、例えば、Ｘ線診断装置１は、流入時間の同定方法やカラーコード（若しくは動画用の周期的カラーコード）の処理条件を略同一として、視点が異なる２つのパラメトリック画像データを生成する。このため、Ｘ線診断装置１は、同一位置を同一色で描出した２つのパラメトリック画像データを生成することができる。
【０１２２】
このように、第１の実施形態に係るＸ線診断装置１は、一定以上のフレームレートを維持しつつ、同一位置を同一色で描出した２つのパラメトリック画像データを生成する。これによれば、Ｘ線診断装置１は、パラメトリックイメージングにおいてステレオ観察可能な画像を提供することができる。
【０１２３】
また、例えば、第１の実施形態に係るＸ線診断装置１は、パラメトリック画像生成処理の処理条件を変更する指示を受け付けて、視点が異なる２つのパラメトリック画像データを再生成する補正処理を実行する。これによれば、Ｘ線診断装置１は、例えば、操作者がパラメトリック画像生成処理の処理条件を変更しても、視点が異なる２つのパラメトリック画像データを常に同一の処理条件で生成／表示することができる。
【０１２４】
また、例えば、第１の実施形態に係るＸ線診断装置１は、第１のＤＳＡ画像データ及び第２のＤＳＡ画像データの同一位置から求まる流入時間を揃える補正処理を行う。そして、パラメトリック画像生成機能３３２は、第１のパラメトリック画像データ及び第２のパラメトリック画像データのうち少なくとも一方を再生成する。これによれば、Ｘ線診断装置１は、例えば、同一位置の流入時間に誤差が生じた場合に、この誤差を補正することができる。
【０１２５】
また、例えば、第１の実施形態に係るＸ線診断装置１は、動画表示において非同期表示モードの要求を受け付けると、複数フレームの第１のパラメトリック画像データ及び複数フレームの第２のパラメトリック画像データのうちいずれか一方について、時系列順に表示されるフレームの時相（フレーム番号あるいはフレーム時間）を変更する。これによれば、Ｘ線診断装置１は、動画であるパラメトリック画像データのカラーコードのステップの時間の単位あるいはそれより細かい時間単位で、パラメトリック画像データのフレームの時相を変更することができる。例えば、動画表示においては、カラーコードの変化が急峻であるので、流入時間に僅かな誤差があるだけで同一位置を同一色で描出できなくなってしまうことがある。具体例を挙げると、１回目の撮像と２回目の撮像との間において、医師が、被検体Ｐに挿入したカテーテル位置を動かしてしまった場合には、流入時間に誤差が生じてしまう。このような場合に、Ｘ線診断装置１は、カラーコードのステップの時間の単位（ΔＴ）あるいはそれより細かい時間単位でフレームの時相を変更するので、非常に細かいピッチでパラメトリック画像データの色を調整することができる。
【０１２６】
（第１の実施形態の変形例）
第１の実施形態において説明した観察方向は、血管の走行方向に応じて自動的に決定することができる。
【０１２７】
例えば、Ｘ線診断装置１は、部位ごとに、各部位における主要な血管の走行方向の情報を記憶する。例えば、Ｘ線診断装置１は、「脳」に対して、「患者上下方向」という情報を記憶装置に記憶している。
【０１２８】
そして、処理回路３３は、撮像部位を指定する旨の指示を操作者から受け付ける。そして、処理回路３３は、各部位における主要な血管の走行方向の情報を参照し、受け付けた撮像部位に対応する血管の走行方向の情報を取得する。例えば、撮像部位として「脳」が指定されると、処理回路３３は、「患者上下方向」の情報を記憶装置から取得する。そして、処理回路３３は、取得した血管の走行方向の情報に基づいて、観察方向を決定する。ここで、処理回路３３は、血管の走行方向に対して直交するように、観察方向を決定する。例えば、処理回路３３は、血管の走行方向が「患者上下方向」であれば、「患者左右方向」を観察方向として決定する。そして、処理回路３３は、決定した観察方向を用いて撮像を行う。なお、ここでは説明を簡単にするために血管の走行方向を「患者上下方向」、「患者左右方向」に限定して説明している。しかし血管の走行方向はその２方向に限定されず、例えばＸ軸方向からＹ軸方向に４５度傾いた走行をしている血管や、あるいはさらにＺ軸方向に４５度傾いた走行をしている血管などある。そのような場合、その走行方向に垂直な平面内で観察方向を変更することが、目的の血管を立体として観察するには最も望ましい。従って血管の走行方向を例えば３ＤＣＴ画像などから同定し、その同定した角度から観察方向を決定しても良い。あるいは立体的な位置関係がある程度分かれば良いような場合は、「患者左右方向」あるいは「患者上下方向」に観察角度を変更してもある程度立体構造は認識できる。この方法は複雑な角度指定が不要になり、例えば「患者左右方向」あるいは「患者上下方向」と言う単純な二者択一になると言うメリットがある。
【０１２９】
これによれば、Ｘ線診断装置１は、操作者が血管の走行方向を意識して観察方向を指定しなくとも、各部位の主要な血管の走行方向に応じた観察方向を決定することができる。なお、上記の内容はあくまで一例であり、上記の例に限定されるものではない。例えば、血管の走行方向の情報は、血管の種類ごとに記憶されてもよい。この場合、操作者が撮像対象となる血管を指定することにより、処理回路３３は、血管の種類に応じた観察方向を決定する。また、血管の走行方向ではなく、血管の走行方向に応じた観察方向が記憶されてもよい。この場合、操作者が撮像部位（若しくは血管）を指定することにより、処理回路３３は、血管の種類に応じた観察方向を記憶装置から読み出すことができる。
【０１３０】
（第２の実施形態）
第１の実施形態では、１つの焦点を有するＸ線管球１１が適用される場合を説明したが、実施形態はこれに限定されるものではない。例えば、Ｘ線管球１１として、２つの焦点を有するステレオＸ線管球が適用される場合であってもよい。ステレオＸ線管球は、２つの焦点を切り替えながらＸ線を照射するＸ線管球である。
【０１３１】
すなわち、第２の実施形態に係るＤＳＡ画像生成機能３３１は、第１のＤＳＡ画像データ及び第２のＤＳＡ画像データのそれぞれを、ステレオＸ線管球を用いて同一の期間に生成する。そして、パラメトリック画像生成機能３３２は、２つの焦点を切り替える時間差を補正した上で、第１のパラメトリック画像データ及び第２のパラメトリック画像データを生成する。
【０１３２】
図１３は、第２の実施形態に係るステレオＸ線管球を用いた撮像について説明するための図である。図１３には、焦点Ｆ１と焦点Ｆ２とを有するステレオＸ線管球を用いて撮像する場合を例示する。
【０１３３】
図１３に示すように、ＤＳＡ画像生成機能３３１は、Ｘ線トリガーボタンの押下を受け付けると、焦点Ｆ１と焦点Ｆ２とを切り替えながら時系列に沿った複数フレームのＤＳＡ画像データを生成する。つまり、ＤＳＡ画像生成機能３３１は、１回の造影剤注入で、複数フレームの第１のＤＳＡ画像データと、複数フレームの第２のＤＳＡ画像データとを撮像する。言い換えると、ＤＳＡ画像生成機能３３１は、造影剤の注入条件やＸ線の照射条件が同一である状況下で複数フレームの第１のＤＳＡ画像データと、複数フレームの第２のＤＳＡ画像データとを撮像する。なお、第２の実施形態に係るＤＳＡ画像生成機能３３１の処理は、１回の造影剤注入で視差のある２つのＤＳＡ画像データを生成する点を除き、第１の実施形態にて説明した内容を同様であるので、説明を省略する。
【０１３４】
ここで、焦点Ｆ１と焦点Ｆ２との切り替えにより、各焦点で撮像されるＤＳＡ画像データには時間差が生じる。例えば、「焦点Ｆ１→焦点Ｆ２→焦点Ｆ１→焦点Ｆ２→・・・」の順で切り替えながら撮像される場合には、焦点Ｆ２で撮像された各フレームの撮像時間は、焦点Ｆ１で撮像された各フレームの撮像時間より時間差（例えば、０．１秒）分遅くなる。
【０１３５】
そこで、パラメトリック画像生成機能３３２は、焦点の切り替えにかかる時間差を補正した上で、第１のパラメトリック画像データ及び第２のパラメトリック画像データを生成する。例えば、「焦点Ｆ１→焦点Ｆ２→焦点Ｆ１→焦点Ｆ２→・・・」の順で切り替えながら撮像される場合には、１フレーム目の画像は焦点Ｆ１と焦点Ｆ２で撮影した画像がペアで記録される。この時ペアで記録された画像には同じ撮像時間が記録されている。例えば１フレーム目の焦点Ｆ１の画像の撮像時間が０秒だったとし、焦点Ｆ１の画像を撮像してから焦点Ｆ２の画像を撮像するまでの時間差を０．１秒とすると、焦点Ｆ２の画像の撮像時間は０．１秒となる。しかし記録時には焦点Ｆ１の画像と焦点Ｆ２の画像は１ペアの画像として記録されるため、１つの撮像時間が記録される。その結果焦点Ｆ２の画像の撮像時間は実際の撮像時間と誤差を生じる。そこでパラメトリック画像生成機能３３２は、焦点Ｆ２で撮像された各フレームの撮像時間それぞれに焦点Ｆ１の画像を撮像してから焦点Ｆ２の画像を撮像するまでの時間差（今回の例では０．１秒）を加算する。これにより、焦点Ｆ１で撮像された各フレームの撮像時間と焦点Ｆ２で撮像された各フレームの撮像時間との時間差を補正する。そして、パラメトリック画像生成機能３３２は、各フレームの撮像時間の時間差が補正されたＤＳＡ画像データを用いて、パラメトリック画像データを生成する。なお、第２の実施形態に係るパラメトリック画像生成機能３３２の処理は、焦点Ｆ１の画像を撮像してから焦点Ｆ２の画像を撮像するまでの時間差を補正する点を除き、第１の実施形態にて説明した内容を同様であるので、説明を省略する。
【０１３６】
このように、２つの焦点を有するステレオＸ線管球が適用される場合であっても、Ｘ線診断装置１は、パラメトリックイメージングにおいて立体視可能な画像を提供することができる。
【０１３７】
なお、上記の説明はあくまで一例であり、上記の例に限定されるものではない。例えば、上記の説明では、ステレオＸ線管球が適用される場合を説明したが、これに限らず、例えば、焦点位置を電磁偏向により切り替えるフライングフォーカスが適用される場合であってもよい。つまり、第２の実施形態は、２つの焦点を切り替えながらＸ線を照射するＸ線管球を用いる場合に適用可能である。
【０１３８】
（第３の実施形態）
上記の実施形態では、撮像部位（注目領域）とＣ型アーム１３の回転中心（Iso-Center）とが略一致している場合を説明したが、実施形態はこれに限定されるものではない。例えば、注目領域をＸ線管球１１に近づけてＸ線光学系での拡大率を上げることも可能である。
【０１３９】
図１４は、第３の実施形態に係るＸ線診断装置１の構成の一例を示す図である。第３の実施形態に係るＸ線診断装置１は、図１に例示したＸ線診断装置１と同様の構成を備え、処理回路３３が注目領域同定機能３３４を更に備える点が相違する。そこで、第３の実施形態では、第３の実施形態と相違する点を中心に説明することとし、第１の実施形態において説明した構成と同様の機能を有する点については、図１と同一の符号を付し、説明を省略する。
【０１４０】
注目領域同定機能３３４は、注目領域を指定する操作を受け付けた場合に、撮像系機器に対する注目領域の位置を同定し、同定した注目領域の位置を撮像するように撮像系機器を制御する。なお、撮像系機器は、例えば、Ｘ線管球１１、検出器１２、及びＣ型アーム１３を含む。
【０１４１】
図１５〜図１９は、第３の実施形態に係る注目領域同定機能３３４の処理を説明するための図である。図１５〜図１９には、Ｘ線光学系での拡大率を上げる場合の注目領域及び撮像系機器の位置関係の遷移を例示する。なお、図１５では、ＳＩＤが１０００［ｍｍ］であり、ＦＯＶが６［ｉｎｃｈ］であり、Ｘ線管球１１の焦点からIso-Centerまでの距離が７００［ｍｍ］である場合を説明する。なお、Ｘ線光学系での拡大率を上げる場合には、Ｘ線画像におけるぼけを防止するために最小の焦点サイズを用いるのが好適である。
【０１４２】
図１５に示すように、注目領域同定機能３３４は、例えば、視点が異なる２つのＸ線透視画像上で注目領域を指定する操作を操作者から受け付けると、Ｘ線光学系での拡大率を決定する。具体的には、注目領域同定機能３３４は、Ｘ線透視画像上で注目領域の中心点及び注目領域の範囲を指定する操作を受け付けると、受け付けた注目領域の範囲を予め設定されたＦＯＶで撮影できるように、Ｘ線光学系での拡大率を決定する。例えば、操作者により指定された注目領域の大きさが５２．５［ｍｍ］であれば、ＳＩＤが１０００［ｍｍ］であり、ＦＯＶが６［ｉｎｃｈ］であるので、注目領域同定機能３３４は、Ｘ線光学系での拡大率を２．８６倍と算出する。この場合、Ｘ線管球１１の焦点から注目領域の中心までの距離は３５０［ｍｍ］となるので、Ｘ線管球１１の焦点からIso-Centerまでの距離（７００［ｍｍ］）と比較して非常に拡大率の大きい撮像となる。
【０１４３】
また、注目領域同定機能３３４は、Ｘ線透視画像上で指定された注目領域の中心点から、３次元空間（実空間）における注目領域の中心位置を同定する。つまり、注目領域同定機能３３４は、３次元空間における注目領域と撮像系機器との位置関係を同定する。
【０１４４】
続いて、図１６に示すように、視差角度を８度としてIso-Centerを中心としてＣ型アーム１３の回転が行われると、検出器（ＦＰＤ）１２上に投影される注目領域の位置がシフトする。そこで、注目領域同定機能３３４は、このシフト量を同定し、同定したシフト量に応じて撮像系機器を制御する。なお、図１６では、回転後におけるＸ線管球１１及び検出器１２の位置を実線で表し、回転前におけるＸ線管球１１及び検出器１２の位置を破線で表す。
【０１４５】
例えば、図１７及び図１８に示すように、注目領域同定機能３３４は、シフトした注目領域が検出器１２上に含まれる場合には、検出器１２上で画像を読み込む領域をシフトさせる。例えば、注目領域同定機能３３４は、検出器１２が１２［ｉｎｃｈ］であれば、６［ｉｎｃｈ］の画像読み込み領域を６９．９２６８［ｍｍ］シフトするように撮像系機器を制御する。なお、図１８では、回転（シフト）後における画像読み込み領域の位置を実線で表し、回転前における画像読み込み領域の位置を破線で表す。
【０１４６】
また、例えば、図１９に示すように、注目領域同定機能３３４は、シフトした注目領域が検出器１２上に含まれない場合には、シフトした注目領域が検出器１２上に含まれるようにＣ型アーム１３を水平移動させる。図１９の例では、注目領域同定機能３３４は、検出器１２の中心に注目領域の中心点が一致するように、Ｃ型アーム１３を水平移動させる。なお、図１９では、水平移動後におけるＸ線管球１１及び検出器１２の位置を実線で表し、水平移動前におけるＸ線管球１１及び検出器１２の位置を破線で表す。
【０１４７】
このように、第３の実施形態に係るＸ線診断装置１は、撮像系機器に対する注目領域の位置を同定し、同定した注目領域の位置を撮像するように撮像系機器を制御する。これによれば、Ｘ線診断装置１は、Ｘ線光学系での拡大率を上げた場合にも、注目領域の中心点が常に視差画像の中心となるように撮像することができる。
【０１４８】
なお、Ｘ線光学系での拡大率を上げた場合には、視差画像である２つのパラメトリック画像データに描出される注目領域の大きさが異なってしまう場合がある。この場合、パラメトリック画像データ上で画像処理として拡大／縮小を行うことにより、注目領域の大きさを揃えるのが好適である。なお、立体視に影響を与えない程度の違いであれば、画像処理を行わずにステレオ表示に用いてもよい。
【０１４９】
なお、上記の説明はあくまで一例であり、上記の例に限定されるものではない。例えば、図１７では、画像読み込み領域をシフトさせる場合を説明したが、これに限らず、例えば、検出器１２自体をシフトさせてもよい。具体的には、検出器１２をシフトさせるための移動機構をＣ型アーム１３が備える場合には、検出器１２自体をシフトさせることができる。
【０１５０】
また、例えば、図１９では、Ｃ型アーム１３の水平移動によって検出器１２の中心に注目領域の中心点を一致させる場合を説明したが、水平移動のみでは一致させることができない場合も考えられる。この場合、例えば、図１７で説明した画像読み込み領域のシフトを組み合わせてもよい。また、検出器１２の中心に注目領域の中心点が一致していなくても、検出器１２の範囲内にあれば、撮像後にパラメトリック画像データ上で画像処理として平行移動を行うことにより一致させてもよい。
【０１５１】
また、上記の説明では、Ｘ線透視画像上で注目領域の中心点及び注目領域の範囲が指定される場合を説明したが、３次元画像上で指定されてもよい。例えば、操作者は、３次元画像のＭＰＲ（Multi  Planar  Reconstructions）断面上に注目領域の立体的な形状パターンを配置することで、３次元画像上で注目領域の中心点及び注目領域の範囲を指定してもよい。
【０１５２】
（その他の実施形態）
上述した実施形態以外にも、種々の異なる形態にて実施されてもよい。
【０１５３】
（Ｘ線パラメトリックイメージング用撮像プログラム）
例えば、上記のＸ線パラメトリックイメージング用撮像プログラムでは、１度の条件設定（撮像条件、処理条件）により２方向の撮像を一連の処理として実行される。
【０１５４】
例えば、Ｘ線診断装置１は、撮像条件及び処理条件のうち少なくとも一つを、第１の差分画像データ及び第２の差分画像データの撮像、並びに、第１の差分画像データ及び第２の差分画像データを用いた処理において一致させることが予め設定された撮像設定情報を有する。この撮像設定情報は、例えば、Ｘ線診断装置１の記憶回路に記憶される。
【０１５５】
そして、処理回路３３は、記憶回路に記憶された撮像設定情報を読み出して、読み出した撮像設定情報に基づいて、ＤＳＡ画像生成機能３３１及びパラメトリック画像生成機能３３２を実行する。具体的には、処理回路３３は、撮像設定情報において撮像条件及び処理条件を設定する操作を操作者から受け付ける。そして、処理回路３３は、受け付けた撮像条件を用いて、第１の方向からの第１の差分画像データの生成と、第２の方向からの第２の差分画像データの生成とを実行する。また、処理回路３３は、受け付けた処理条件を用いて、第１のカラー画像データの生成と、第２のカラー画像データの生成とを実行する。これにより、例えば、操作者は、１度の条件設定で２方向の撮像を行うことが可能となる。
【０１５６】
（画像処理装置）
また、上記の実施形態にて説明した内容は、Ｘ線診断装置に限らず、任意の画像処理装置においても実現可能である。
【０１５７】
すなわち、ＤＳＡ画像生成機能３３１は、第１の方向から所定の撮像条件で時系列に撮影した造影画像データと非造影画像データを差分して第１の差分画像データを生成し、第２の方向から所定の撮像条件と略同一の撮像条件で時系列に撮影した造影画像データと非造影画像データを差分して第２の差分画像データを生成する。そして、パラメトリック画像生成機能３３２は、第１の差分画像データの各位置に、画素値の経時的変化に応じたカラーを割り当てた第１のカラー画像データを所定の処理条件で生成するとともに、第２の差分画像データの各位置に、画素値の経時的変化に応じたカラーを割り当てた第２のカラー画像データを所定の処理条件と略同一の処理条件で生成する。そして、表示制御機能３３３は、第１のカラー画像データ及び第２のカラー画像データに基づく立体視画像を表示させる。
【０１５８】
また、図示した各装置の各構成要素は機能概念的なものであり、必ずしも物理的に図示の如く構成されていることを要しない。すなわち、各装置の分散・統合の具体的形態は図示のものに限られず、その全部又は一部を、各種の負荷や使用状況等に応じて、任意の単位で機能的又は物理的に分散・統合して構成することができる。更に、各装置にて行なわれる各処理機能は、その全部又は任意の一部が、ＣＰＵ及び当該ＣＰＵにて解析実行されるプログラムにて実現され、或いは、ワイヤードロジックによるハードウェアとして実現され得る。
【０１５９】
また、上記の実施形態及び変形例において説明した各処理のうち、自動的に行なわれるものとして説明した処理の全部又は一部を手動的に行なうこともでき、或いは、手動的に行なわれるものとして説明した処理の全部又は一部を公知の方法で自動的に行なうこともできる。この他、上記文書中や図面中で示した処理手順、制御手順、具体的名称、各種のデータやパラメータを含む情報については、特記する場合を除いて任意に変更することができる。
【０１６０】
また、上記の実施形態及び変形例で説明した撮像方法（画像診断支援方法）は、予め用意された撮像プログラムをパーソナルコンピュータやワークステーション等のコンピュータで実行することによって実現することができる。この撮像方法は、インターネット等のネットワークを介して配布することができる。また、この撮像方法は、ハードディスク、フレキシブルディスク（ＦＤ）、ＣＤ−ＲＯＭ、ＭＯ、ＤＶＤ等のコンピュータで読み取り可能な記録媒体に記録され、コンピュータによって記録媒体から読み出されることによって実行することもできる。
【０１６１】
以上説明した少なくともひとつの実施形態によれば、パラメトリックイメージングにおいて立体視可能な画像を提供することができる。
【０１６２】
本発明のいくつかの実施形態を説明したが、これらの実施形態は、例として提示したものであり、発明の範囲を限定することは意図していない。これら実施形態は、その他の様々な形態で実施されることが可能であり、発明の要旨を逸脱しない範囲で、種々の省略、置き換え、変更を行うことができる。これら実施形態やその変形は、発明の範囲や要旨に含まれると同様に、特許請求の範囲に記載された発明とその均等の範囲に含まれるものである。
【符号の説明】
【０１６３】
１      Ｘ線診断装置
３３      処理回路
３３１      ＤＳＡ画像生成機能
３３２      パラメトリック画像生成機能
３３３      表示制御機能
【図１】
2018000320.tif 000003
【図２Ａ】
2018000320.tif 000004
【図２Ｂ】
2018000320.tif 000005
【図３】
2018000320.tif 000006
【図４】
2018000320.tif 000007
【図５】
2018000320.tif 000008
【図６】
2018000320.tif 000009
【図７Ａ】
2018000320.tif 000010
【図７Ｂ】
2018000320.tif 000011
【図８】
2018000320.tif 000012
【図９】
2018000320.tif 000013
【図１０】
2018000320.tif 000014
【図１１】
2018000320.tif 000015
【図１２】
2018000320.tif 000016
【図１３】
2018000320.tif 000017
【図１４】
2018000320.tif 000018
【図１５】
2018000320.tif 000019
【図１６】
2018000320.tif 000020
【図１７】
2018000320.tif 000021
【図１８】
2018000320.tif 000022
【図１９】
2018000320.tif 000023