ホーム > 特許ランキング > ザイオソフト株式会社
特開2024-151639医用画像処理装置、医用画像処理方法、及び医用画像処理プログラム
(19)【発行国】日本国特許庁(JP)
(12)【公報種別】公開特許公報(A)
(11)【公開番号】P2024151639
(43)【公開日】2024-10-25
(54)【発明の名称】医用画像処理装置、医用画像処理方法、及び医用画像処理プログラム
(51)【国際特許分類】
A61B 6/03 20060101AFI20241018BHJP
【ＦＩ】
A61B6/03 360E
【審査請求】未請求
【請求項の数】13
【出願形態】ＯＬ
(21)【出願番号】P 2023065151
(22)【出願日】2023-04-12
(71)【出願人】
【識別番号】500109320
【氏名又は名称】ザイオソフト株式会社
(74)【代理人】
【識別番号】110002000
【氏名又は名称】弁理士法人栄光事務所
(72)【発明者】
【氏名】茅野 秀介
(72)【発明者】
【氏名】那須 悠
【テーマコード（参考）】
4C093
【Ｆターム（参考）】
4C093AA22
4C093CA35
4C093DA02
4C093FB12
4C093FF42
4C093FF46
4C093FG04
4C093FG13
(57)【要約】
【課題】ＣＰＲ画像上での二重描画の発生を抑制でき、ＣＰＲ面上の様子を把握し易くできる医用画像処理装置を提供する。
【解決手段】医用画像処理装置は、ユーザインタフェースが、第１の方向ベクトルを回転させる回転操作を受け付け、プロセッサが、被検体のボリュームデータを取得し、ボリュームデータにおける３次元空間上のパスを取得し、前記第１の方向ベクトルを回転させる場合に、第１の方向ベクトルに垂直な第１の投影面に対して第１の方向ベクトルに沿ってパスを投影することで得られる第１の投影パスが第１の投影面上で重複する第１の重複度合いが、所定の閾値以下である第１の推奨方向ベクトルを決定し、パスと第１の推奨方向ベクトルとに基づいて、パス上の各点を通過し第１の推奨方向ベクトルと平行な直線の集合によって定義される第１の曲面を決定し、第１の曲面上のボリュームデータに基づいて、第１の曲面を可視化した第１のＣＰＲ画像を生成する。
【選択図】図１０
【解決手段】医用画像処理装置は、ユーザインタフェースが、第１の方向ベクトルを回転させる回転操作を受け付け、プロセッサが、被検体のボリュームデータを取得し、ボリュームデータにおける３次元空間上のパスを取得し、前記第１の方向ベクトルを回転させる場合に、第１の方向ベクトルに垂直な第１の投影面に対して第１の方向ベクトルに沿ってパスを投影することで得られる第１の投影パスが第１の投影面上で重複する第１の重複度合いが、所定の閾値以下である第１の推奨方向ベクトルを決定し、パスと第１の推奨方向ベクトルとに基づいて、パス上の各点を通過し第１の推奨方向ベクトルと平行な直線の集合によって定義される第１の曲面を決定し、第１の曲面上のボリュームデータに基づいて、第１の曲面を可視化した第１のＣＰＲ画像を生成する。
【選択図】図１０
【特許請求の範囲】
【請求項1】
プロセッサとユーザインタフェースとを備える医用画像処理装置であって、
前記ユーザインタフェースは、第１の方向ベクトルを回転させる回転操作を受け付け、
前記プロセッサは、
被検体のボリュームデータを取得し、
前記ボリュームデータにおける３次元空間上のパスを取得し、
前記第１の方向ベクトルを回転させうる範囲内で、前記第１の方向ベクトルに垂直な第１の投影面に対して前記第１の方向ベクトルに沿って前記パスを投影することで得られる第１の投影パスが前記第１の投影面上で重複する第１の重複度合いが、所定の第１の閾値以下である第１の推奨方向ベクトルを決定し、
前記パスと前記第１の推奨方向ベクトルとに基づいて、前記パス上の各点を通過し前記第１の推奨方向ベクトルと平行な直線の集合によって定義される第１の曲面を決定し、
前記第１の曲面上の前記ボリュームデータに基づいて、前記第１の曲面を可視化した第１のＣＰＲ画像を生成する、
医用画像処理装置。
前記ユーザインタフェースは、第１の方向ベクトルを回転させる回転操作を受け付け、
前記プロセッサは、
被検体のボリュームデータを取得し、
前記ボリュームデータにおける３次元空間上のパスを取得し、
前記第１の方向ベクトルを回転させうる範囲内で、前記第１の方向ベクトルに垂直な第１の投影面に対して前記第１の方向ベクトルに沿って前記パスを投影することで得られる第１の投影パスが前記第１の投影面上で重複する第１の重複度合いが、所定の第１の閾値以下である第１の推奨方向ベクトルを決定し、
前記パスと前記第１の推奨方向ベクトルとに基づいて、前記パス上の各点を通過し前記第１の推奨方向ベクトルと平行な直線の集合によって定義される第１の曲面を決定し、
前記第１の曲面上の前記ボリュームデータに基づいて、前記第１の曲面を可視化した第１のＣＰＲ画像を生成する、
医用画像処理装置。
【請求項2】
プロセッサとユーザインタフェースとを備える医用画像処理装置であって、
前記ユーザインタフェースは、第１の方向ベクトルを回転させる回転操作を受け付け、
前記プロセッサは、
被検体のボリュームデータを取得し、
前記ボリュームデータにおける３次元空間上のパスを取得し、
前記第１の方向ベクトルを回転させうる範囲内で、前記第１の方向ベクトルに垂直な第１の投影面に対して前記第１の方向ベクトルに沿って前記パスを投影することで得られる第１の投影パスが前記第１の投影面上で重複する第１の重複度合いが、所定の第１の基準方向ベクトルに垂直な第１の基準投影面に対して前記第１の基準方向ベクトルに沿って前記パスを投影することで得られる第１の基準投影パスが前記第１の基準投影面上で重複する度合いである第１の基準重複度合いよりも小さい、第１の推奨方向ベクトルを決定し、
前記パスと前記第１の推奨方向ベクトルとに基づいて、前記パス上の各点を通過し前記第１の推奨方向ベクトルと平行な直線の集合によって定義される第１の曲面を決定し、
前記第１の曲面上の前記ボリュームデータに基づいて、前記第１の曲面を可視化した第１のＣＰＲ画像を生成する、
医用画像処理装置。
前記ユーザインタフェースは、第１の方向ベクトルを回転させる回転操作を受け付け、
前記プロセッサは、
被検体のボリュームデータを取得し、
前記ボリュームデータにおける３次元空間上のパスを取得し、
前記第１の方向ベクトルを回転させうる範囲内で、前記第１の方向ベクトルに垂直な第１の投影面に対して前記第１の方向ベクトルに沿って前記パスを投影することで得られる第１の投影パスが前記第１の投影面上で重複する第１の重複度合いが、所定の第１の基準方向ベクトルに垂直な第１の基準投影面に対して前記第１の基準方向ベクトルに沿って前記パスを投影することで得られる第１の基準投影パスが前記第１の基準投影面上で重複する度合いである第１の基準重複度合いよりも小さい、第１の推奨方向ベクトルを決定し、
前記パスと前記第１の推奨方向ベクトルとに基づいて、前記パス上の各点を通過し前記第１の推奨方向ベクトルと平行な直線の集合によって定義される第１の曲面を決定し、
前記第１の曲面上の前記ボリュームデータに基づいて、前記第１の曲面を可視化した第１のＣＰＲ画像を生成する、
医用画像処理装置。
【請求項3】
前記ユーザインタフェースは、第２の方向ベクトルを回転させる回転操作を受け付け、
前記プロセッサは、
前記第１の方向ベクトルに対して所定の角度を維持して第２の方向ベクトルを回転させ、
前記第１の方向ベクトル及び前記第２の方向ベクトルを回転させる場合に、前記第１の重複度合いと、前記第２の方向ベクトルに垂直な第２の投影面に対して前記第２の方向ベクトルに沿って前記パスを投影することで得られる第２の投影パスが前記第２の投影面上で重複する第２の重複度合いと、を統合した統合重複度合いが、所定の第２の閾値以下である場合、前記第２の閾値以下の前記統合重複度合いが得られた前記第１の重複度合いが得られた前記第１の方向ベクトルを前記第１の推奨方向ベクトルに決定し、前記第２の閾値以下の前記統合重複度合いが得られた前記第２の重複度合いが得られた前記第２の方向ベクトルを第２の推奨方向ベクトルに決定し、
前記パスと前記第２の推奨方向ベクトルとに基づいて、前記パス上の各点を通過し前記第２の推奨方向ベクトルと平行な直線の集合によって定義される第２の曲面を決定し、
前記第２の曲面上のボリュームデータに基づいて、前記第２の曲面を可視化した第２のＣＰＲ画像を生成し、
前記第１のＣＰＲ画像と前記第２のＣＰＲ画像とを同時に表示させる、
請求項１に記載の医用画像処理装置。
前記プロセッサは、
前記第１の方向ベクトルに対して所定の角度を維持して第２の方向ベクトルを回転させ、
前記第１の方向ベクトル及び前記第２の方向ベクトルを回転させる場合に、前記第１の重複度合いと、前記第２の方向ベクトルに垂直な第２の投影面に対して前記第２の方向ベクトルに沿って前記パスを投影することで得られる第２の投影パスが前記第２の投影面上で重複する第２の重複度合いと、を統合した統合重複度合いが、所定の第２の閾値以下である場合、前記第２の閾値以下の前記統合重複度合いが得られた前記第１の重複度合いが得られた前記第１の方向ベクトルを前記第１の推奨方向ベクトルに決定し、前記第２の閾値以下の前記統合重複度合いが得られた前記第２の重複度合いが得られた前記第２の方向ベクトルを第２の推奨方向ベクトルに決定し、
前記パスと前記第２の推奨方向ベクトルとに基づいて、前記パス上の各点を通過し前記第２の推奨方向ベクトルと平行な直線の集合によって定義される第２の曲面を決定し、
前記第２の曲面上のボリュームデータに基づいて、前記第２の曲面を可視化した第２のＣＰＲ画像を生成し、
前記第１のＣＰＲ画像と前記第２のＣＰＲ画像とを同時に表示させる、
請求項１に記載の医用画像処理装置。
【請求項4】
前記ユーザインタフェースは、第２の方向ベクトルを回転させる回転操作を受け付け、
前記プロセッサは、
前記第１の方向ベクトルに対して所定の角度を維持して第２の方向ベクトルを回転させ、
前記第１の方向ベクトル及び前記第２の方向ベクトルを回転させる場合に、前記第１の重複度合いと、前記第２の方向ベクトルに垂直な第２の投影面に対して前記第２の方向ベクトルに沿って前記パスを投影することで得られる第２の投影パスが前記第２の投影面上で重複する第２の重複度合いと、を統合した統合重複度合いが、前記第１の基準重複度合いと、所定の第２の基準方向ベクトルに垂直な第２の基準投影面に対して前記第２の基準方向ベクトルに沿って前記パスを投影することで得られる第２の基準投影パスが前記第２の基準投影面上で重複する度合いである第２の基準重複度合いと、を統合した基準統合重複度合いよりも小さい場合、前記基準統合重複度合いよりも小さい前記統合重複度合いが得られた前記第１の重複度合いが得られた前記第１の方向ベクトルを前記第１の推奨方向ベクトルに決定し、前記基準統合重複度合いよりも小さい前記統合重複度合いが得られた前記第２の重複度合いが得られた前記第２の方向ベクトルを第２の推奨方向ベクトルに決定し、
前記パスと前記第２の推奨方向ベクトルとに基づいて、前記パス上の各点を通過し前記第２の推奨方向ベクトルと平行な直線の集合によって定義される第２の曲面を決定し、
前記第２の曲面上のボリュームデータに基づいて、前記第２の曲面を可視化した第２のＣＰＲ画像を生成し、
前記第１のＣＰＲ画像と前記第２のＣＰＲ画像とを同時に表示させる、
請求項２に記載の医用画像処理装置。
前記プロセッサは、
前記第１の方向ベクトルに対して所定の角度を維持して第２の方向ベクトルを回転させ、
前記第１の方向ベクトル及び前記第２の方向ベクトルを回転させる場合に、前記第１の重複度合いと、前記第２の方向ベクトルに垂直な第２の投影面に対して前記第２の方向ベクトルに沿って前記パスを投影することで得られる第２の投影パスが前記第２の投影面上で重複する第２の重複度合いと、を統合した統合重複度合いが、前記第１の基準重複度合いと、所定の第２の基準方向ベクトルに垂直な第２の基準投影面に対して前記第２の基準方向ベクトルに沿って前記パスを投影することで得られる第２の基準投影パスが前記第２の基準投影面上で重複する度合いである第２の基準重複度合いと、を統合した基準統合重複度合いよりも小さい場合、前記基準統合重複度合いよりも小さい前記統合重複度合いが得られた前記第１の重複度合いが得られた前記第１の方向ベクトルを前記第１の推奨方向ベクトルに決定し、前記基準統合重複度合いよりも小さい前記統合重複度合いが得られた前記第２の重複度合いが得られた前記第２の方向ベクトルを第２の推奨方向ベクトルに決定し、
前記パスと前記第２の推奨方向ベクトルとに基づいて、前記パス上の各点を通過し前記第２の推奨方向ベクトルと平行な直線の集合によって定義される第２の曲面を決定し、
前記第２の曲面上のボリュームデータに基づいて、前記第２の曲面を可視化した第２のＣＰＲ画像を生成し、
前記第１のＣＰＲ画像と前記第２のＣＰＲ画像とを同時に表示させる、
請求項２に記載の医用画像処理装置。
【請求項5】
前記第１の重複度合いは、前記第１の投影面上で前記第１の投影パスが交差又は所定距離以内に接近する度合いである、
請求項１又は２に記載の医用画像処理装置。
請求項１又は２に記載の医用画像処理装置。
【請求項6】
前記ユーザインタフェースは、前記第１の推奨方向ベクトルの決定のタイミングを指示する操作を受け付け、
前記プロセッサは、
前記ユーザインタフェースにより前記第１の推奨方向ベクトルの決定のタイミンングを指示する操作を受け付けた際、前記第１の推奨方向ベクトルを決定する、
請求項１又は２に記載の医用画像処理装置。
前記プロセッサは、
前記ユーザインタフェースにより前記第１の推奨方向ベクトルの決定のタイミンングを指示する操作を受け付けた際、前記第１の推奨方向ベクトルを決定する、
請求項１又は２に記載の医用画像処理装置。
【請求項7】
前記プロセッサは、前記第１の推奨方向ベクトルを前記第１の方向ベクトルの初期値として決定する、
請求項１又は２に記載の医用画像処理装置。
請求項１又は２に記載の医用画像処理装置。
【請求項8】
前記プロセッサは、
前記第１のＣＰＲ画像において、前記投影面上で前記投影パスが交差、重複又は接近していると判定された領域であり、前記ボリュームデータ内の同一の点が複数回描画される領域である二重描画領域を表示デバイスに表示させる、
請求項１又は２に記載の医用画像処理装置。
前記第１のＣＰＲ画像において、前記投影面上で前記投影パスが交差、重複又は接近していると判定された領域であり、前記ボリュームデータ内の同一の点が複数回描画される領域である二重描画領域を表示デバイスに表示させる、
請求項１又は２に記載の医用画像処理装置。
【請求項9】
前記プロセッサは、
前記パスを複数のパス部分に分割し、
前記複数のパス部分のうちの第１のパス部分を指定し、
前記第１の重複度合いは、前記第１のパス部分を投影することで得られる第１の投影パス部分が前記第１の投影面上で重複する度合いである、
請求項１又は２に記載の医用画像処理装置。
前記パスを複数のパス部分に分割し、
前記複数のパス部分のうちの第１のパス部分を指定し、
前記第１の重複度合いは、前記第１のパス部分を投影することで得られる第１の投影パス部分が前記第１の投影面上で重複する度合いである、
請求項１又は２に記載の医用画像処理装置。
【請求項10】
被検体のボリュームデータを取得するステップと、
前記ボリュームデータにおける３次元空間上のパスを取得するステップと、
第１の方向ベクトルを回転させる回転操作を受け付けるユーザインタフェースを介して前記第１の方向ベクトルを回転させる場合に、前記第１の方向ベクトルに垂直な第１の投影面に対して前記第１の方向ベクトルに沿って前記パスを投影することで得られる第１の投影パスが前記第１の投影面上で重複する第１の重複度合いが、所定の閾値以下である第１の推奨方向ベクトルを決定するステップと、
前記パスと前記第１の推奨方向ベクトルとに基づいて、前記パス上の各点を通過し前記第１の推奨方向ベクトルと平行な直線の集合によって定義される第１の曲面を決定し、
前記第１の曲面上の前記ボリュームデータに基づいて、前記第１の曲面を可視化した第１のＣＰＲ画像を生成するステップと、
を有する医用画像処理方法。
前記ボリュームデータにおける３次元空間上のパスを取得するステップと、
第１の方向ベクトルを回転させる回転操作を受け付けるユーザインタフェースを介して前記第１の方向ベクトルを回転させる場合に、前記第１の方向ベクトルに垂直な第１の投影面に対して前記第１の方向ベクトルに沿って前記パスを投影することで得られる第１の投影パスが前記第１の投影面上で重複する第１の重複度合いが、所定の閾値以下である第１の推奨方向ベクトルを決定するステップと、
前記パスと前記第１の推奨方向ベクトルとに基づいて、前記パス上の各点を通過し前記第１の推奨方向ベクトルと平行な直線の集合によって定義される第１の曲面を決定し、
前記第１の曲面上の前記ボリュームデータに基づいて、前記第１の曲面を可視化した第１のＣＰＲ画像を生成するステップと、
を有する医用画像処理方法。
【請求項11】
被検体のボリュームデータを取得するステップと、
前記ボリュームデータにおける３次元空間上のパスを取得するステップと、
第１の方向ベクトルを回転させる回転操作を受け付けるユーザインタフェースを介して前記第１の方向ベクトルを回転させる場合に、前記第１の方向ベクトルに垂直な第１の投影面に対して前記第１の方向ベクトルに沿って前記パスを投影することで得られる第１の投影パスが前記第１の投影面上で重複する第１の重複度合いが、所定の第１の基準方向ベクトルに垂直な第１の基準投影面に対して前記第１の基準方向ベクトルに沿って前記パスを投影することで得られる第１の基準投影パスが前記第１の基準投影面上で重複する度合いよりも小さい、第１の推奨方向ベクトルを決定するステップと、
前記パスと前記第１の推奨方向ベクトルとに基づいて、前記パス上の各点を通過し前記第１の推奨方向ベクトルと平行な直線の集合によって定義される第１の曲面を決定するステップと、
前記第１の曲面上の前記ボリュームデータに基づいて、前記第１の曲面を可視化した第１のＣＰＲ画像を生成するステップと、
を有する医用画像処理方法。
前記ボリュームデータにおける３次元空間上のパスを取得するステップと、
第１の方向ベクトルを回転させる回転操作を受け付けるユーザインタフェースを介して前記第１の方向ベクトルを回転させる場合に、前記第１の方向ベクトルに垂直な第１の投影面に対して前記第１の方向ベクトルに沿って前記パスを投影することで得られる第１の投影パスが前記第１の投影面上で重複する第１の重複度合いが、所定の第１の基準方向ベクトルに垂直な第１の基準投影面に対して前記第１の基準方向ベクトルに沿って前記パスを投影することで得られる第１の基準投影パスが前記第１の基準投影面上で重複する度合いよりも小さい、第１の推奨方向ベクトルを決定するステップと、
前記パスと前記第１の推奨方向ベクトルとに基づいて、前記パス上の各点を通過し前記第１の推奨方向ベクトルと平行な直線の集合によって定義される第１の曲面を決定するステップと、
前記第１の曲面上の前記ボリュームデータに基づいて、前記第１の曲面を可視化した第１のＣＰＲ画像を生成するステップと、
を有する医用画像処理方法。
【請求項12】
請求項１０に記載の医用画像処理方法をコンピュータに実行させるための医用画像処理プログラム。
【請求項13】
請求項１１に記載の医用画像処理方法をコンピュータに実行させるための医用画像処理プログラム。
【発明の詳細な説明】
【技術分野】
【0001】
本開示は、医用画像処理装置、医用画像処理方法、及び医用画像処理プログラムに関する。
【背景技術】
【0002】
従来、医用画像処理装置がＣＰＲ（ＣｕｒｖｅｄＰｌａｎａｒＲｅｃｏｎｓｔｒｕｃｔｉｏｎ）画像を表示することが知られている。ＣＰＲ画像は、ＣＴ装置やＭＲ（ＭａｇｎｅｔｉｃＲｅｓｏｎａｎｃｅ）装置での撮像により得られたボリュームデータにおける任意の曲線をＣＰＲ方向ベクトルの方向に平行移動させて得られる曲面（ＣＰＲ面）上のボクセルで構成される画像である。従来、ＣＰＲ面を規定するためのＣＰＲ方向ベクトルの方向をユーザ操作により指定する画像処理装置が知られている（特許文献１参照）。
【先行技術文献】
【特許文献】
【0003】
【特許文献1】特許第４３４１４３４号公報
【発明の概要】
【発明が解決しようとする課題】
【0004】
特許文献１の画像処理装置では、ユーザ所望の方向にＣＰＲ方向ベクトルを指定可能である。しかし、パスとＣＰＲ方向ベクトルとの組み合わせによっては、ＣＰＲ画像上でボリュームデータ内の同一の点が複数回描画されること（二重描画）が多くなり、ユーザがＣＰＲ面上の様子を把握し難くなることがあった。
【0005】
本開示は、上記事情に鑑みてされたものであって、ＣＰＲ画像上での二重描画の発生を抑制でき、ＣＰＲ面上の様子を把握し易くできる医用画像処理装置、医用画像処理方法、及び医用画像処理プログラムを提供する。
【課題を解決するための手段】
【0006】
本開示の一態様は、プロセッサを備える医用画像処理装置であって、前記プロセッサは、被検体のボリュームデータを取得し、前記ボリュームデータにおける３次元空間上のパスを取得し、前記パス上の各点を通過する第１の方向ベクトルを複数仮定し、仮定された複数の第１の方向ベクトルのうち、前記第１の方向ベクトルに垂直な第１の投影面に対して前記第１の方向ベクトルに沿って前記パスを投影することで得られる第１の投影パスが前記第１の投影面上で重複する第１の重複度合いが、所定の閾値以下である第１の推奨方向ベクトルを決定し、前記パスと前記第１の推奨方向ベクトルとに基づいて、前記パス上の各点を通過し前記第１の推奨方向ベクトルと平行な直線の集合によって定義される第１の曲面を決定し、前記第１の曲面上の前記ボリュームデータに基づいて、前記第１の曲面を可視化した第１のＣＰＲ画像を生成する、医用画像処理装置である。
【0007】
本開示の一態様は、被検体のボリュームデータを取得するステップと、前記ボリュームデータにおける３次元空間上のパスを取得するステップと、前記パス上の各点を通過する第１の方向ベクトルを複数仮定するステップと、仮定された複数の第１の方向ベクトルのうち、前記第１の方向ベクトルに垂直な第１の投影面に対して前記第１の方向ベクトルに沿って前記パスを投影することで得られる第１の投影パスが前記第１の投影面上で重複する第１の重複度合いが、所定の閾値以下である第１の推奨方向ベクトルを決定するステップと、前記パスと前記第１の推奨方向ベクトルとに基づいて、前記パス上の各点を通過し前記第１の推奨方向ベクトルと平行な直線の集合によって定義される第１の曲面を決定し、前記第１の曲面上の前記ボリュームデータに基づいて、前記第１の曲面を可視化した第１のＣＰＲ画像を生成するステップと、を有する医用画像処理方法である。
【0008】
本開示の一態様は、ＣＰＲ画像上での二重描画の発生を抑制でき、ＣＰＲ面上の様子を把握し易くなる。
【発明の効果】
【0009】
本開示によれば、ＣＰＲ画像上での二重描画の発生を抑制でき、ＣＰＲ面上の様子を把握し易くできる。
【図面の簡単な説明】
【0010】
【図1】第１の実施形態における医用画像処理装置のハードウェア構成例を示すブロック図
【図2】医用画像処理装置の機能構成例を示すブロック図
【図3】パスとＣＰＲ方向ベクトルとＣＰＲ面の関係を説明するための図
【図4】ＣＰＲ軸とＣＰＲ方向ベクトルの回転とを説明するための図
【図5】基準方向ベクトルの一例を示す図
【図6】任意の臓器に含まれる血管をパスとしたＣＰＲ面の一例を示す模式図
【図7】血管をパスとしたＣＰＲ面の一例を示す画像例を示す図
【図8A】拡張投影パスが重複領域を有しない場合の一例を示す図
【図8B】投影パスが交差を有し、拡張投影パスが重複領域を有する場合の一例を示す図
【図8C】投影パスＰＨＰが交差を有さずに接近しており、拡張投影パスが重複領域を有する場合の一例を示す図
【図9A】投影パスの一例を示す図（その１）
【図9B】投影パスの一例を示す図（その２）
【図10】医用画像処理装置の動作例を示すフローチャート
【図11】ＣＰＲ画像内での重複領域ＴＲの表示例を示す図
【図12】ＣＰＲ画像内での二重描画領域の表示例を示す図
【図13】パスと複数のＣＰＲ方向ベクトルと複数のＣＰＲ面の関係を説明するための図
【図14A】２つの基準方向ベクトルに基づく２つのＣＰＲ画像の表示例を示す図
【図14B】２つの推奨方向ベクトルに基づく２つのＣＰＲ画像の表示例を示す図
【図15】パス部分や投影パス部分を説明するための図
【図16】ＣＰＲ方向ベクトルの投影面における投影パスの重複とＣＰＲ面との関係を説明するための図（その１）
【図17】ＣＰＲ方向ベクトルの投影面における投影パスの重複とＣＰＲ面との関係を説明するための図（その２）
【発明を実施するための形態】
【0011】
以下、本開示の実施形態について、図面を用いて説明する。
【0012】
（本開示の一形態を得るに至った経緯）
図１６及び図１７は、ＣＰＲ画像において二重描画が発生する仕組みを、ＣＰＲ方向ベクトルの投影面における投影パスの重複とＣＰＲ面との関係を用いて説明するための図である。
図１６及び図１７は、ＣＰＲ画像において二重描画が発生する仕組みを、ＣＰＲ方向ベクトルの投影面における投影パスの重複とＣＰＲ面との関係を用いて説明するための図である。
【0013】
図１６では、ボリュームデータに含まれる所定の平面ＨＭＸ上にＵ字型のパスＰＨＸがあるとする。ＣＰＲ方向ベクトルＶＴＸの方向がこの平面ＨＭＸに対して垂直な方向である場合、パスＰＨＸが投影される投影面ＳＰＸにおいて、投影されたパス（投影パスＰＨＰＸ）は所定の平面ＨＭＸ上と同様にＵ字型のパスとなる。このような場合、投影パスＰＨＰＸには重複箇所が無く、よって投影パスＰＨＰＸ上の点に対応するＣＰＲ画像上の部位はそれぞれ一度ずつしか描画されないので、二重描画の無い良好なＣＰＲ画像が得られる。
【0014】
図１７では、ボリュームデータに含まれる所定の平面ＨＭＸ上にＵ字型のパスＰＨＸがあるとする。ＣＰＲ方向ベクトルＶＴＸの方向がこの平面ＨＭＸに対して平行な方向である場合、パスＰＨＸが投影される投影面ＳＰＸにおいて、投影パスＰＨＰＸは、一方向に進んだ後、その逆方向に進む二重線となる。このような場合、投影パスＰＨＰＸは重複し、よって投影パスＰＨＰＸ上の点に対応するＣＰＲ画像上の部位はそれぞれ２回描画される。このように、パスＰＨＸとＣＰＲ方向ベクトルＶＴＸとの組み合わせによっては、ＣＰＲ画像上でパスＰＨＸ上の同一の点に対応するＣＰＲ面上の線が複数回描画され、ＣＰＲ画像に二重描画が発生し、ユーザがＣＰＲ面ＳＣ上の様子を把握し難くなることがあった。
【0015】
以下の実施形態では、ＣＰＲ画像上での二重描画の発生を抑制でき、ＣＰＲ面上の様子を把握し易くできる医用画像処理装置、医用画像処理方法、及び医用画像処理プログラムについて説明する。
【0016】
（第１の実施形態）
図１は、第１の実施形態における医用画像処理装置１００の構成例を示すブロック図である。医用画像処理装置１００は、ポート１１０、ＵＩ１２０、ディスプレイ１３０、プロセッサ１４０、及びメモリ１５０を備える。
図１は、第１の実施形態における医用画像処理装置１００の構成例を示すブロック図である。医用画像処理装置１００は、ポート１１０、ＵＩ１２０、ディスプレイ１３０、プロセッサ１４０、及びメモリ１５０を備える。
【0017】
医用画像処理装置１００には、ＣＴ装置２００が接続される。医用画像処理装置１００は、ＣＴ装置２００からボリュームデータを取得し、取得されたボリュームデータに対して処理を行う。医用画像処理装置１００は、ＰＣとＰＣに搭載されたソフトウェアにより構成されてもよい。
【0018】
ＣＴ装置２００は、被検体へＸ線を照射し、体内の組織によるＸ線の吸収の違いを利用して、画像（ＣＴ画像）を撮像する。被検体は、生体、人体、又は動物等を含んでよい。ＣＴ装置２００は、Ｘ線検出器からシノグラムを取得し、シノグラムに基づいて、画像再構成によって被検体の断層画像（スライス画像、スライスデータともいう）を生成する。ＣＴ装置２００は、スライスデータに基づいて、例えばスライスデータを積層して、ボリュームデータを生成する。スライスデータ及びボリュームデータは、被検体内部の任意の箇所の情報を含む。ＣＴ装置２００は、ＣＴ画像としてのボリュームデータを医用画像処理装置１００へ、有線回線又は無線回線を介して送信する。ＣＴ画像の撮像には、ＣＴ撮像に関する撮像条件や造影剤の投与に関する造影条件が考慮されてよい。なお、造影は、血管の他、消化器官、胆管等に対して行われてよい。造影は、臓器の特性に応じて異なるタイミングで複数回実施されてよい。
【0019】
医用画像処理装置１００内のポート１１０は、通信ポート、外部装置接続ポート、又は組み込みデバイスへの接続ポート等を含み、ＣＴ画像から得られたボリュームデータを取得する。取得されたボリュームデータは、直ぐにプロセッサ１４０に送られて各種処理されてもよいし、メモリ１５０において保管された後、必要時にプロセッサ１４０へ送られて各種処理されてもよい。また、ボリュームデータは、記録媒体や記録メディアを介して取得されてもよい。また、ボリュームデータは中間データ、圧縮データ、シノグラム、又はスライスデータ等の形で取得されてもよい。また、ボリュームデータは医用画像処理装置１００に取り付けられたセンサデバイスからの情報から取得されてもよい。
【0020】
ＵＩ１２０は、タッチパネル、ポインティングデバイス、キーボード、又はマイクロホンを含んでよい。ＵＩ１２０は、医用画像処理装置１００のユーザから、任意の入力操作を受け付ける。ユーザは、医師、放射線技師、学生、又はその他医療従事者（Paramedic Staff）を含んでよい。
【0021】
ＵＩ１２０は、各種操作を受け付ける。例えば、ボリュームデータやボリュームデータに基づく画像（例えば後述する３次元画像、２次元画像）における、関心領域（ＲＯＩ）の指定や輝度条件（例えばウィンドウ情報）の設定等の操作を受け付ける。関心領域は、各種組織（例えば、血管、気管支、臓器、器官、骨、脳）の領域を含んでよい。組織は、病変組織、正常組織、又は腫瘍組織等を含んでよい。ウィンドウ情報は、ウィンドウ幅（ＷＷ：Window Width）及びウィンドウレベル（ＷＬ：Window Level）の少なくとも一方を含み、「ＷＷ／ＷＬ」とも記載する。ウィンドウ情報は、表示される画像の輝度を調整する情報である。また、ＵＩ１２０は、後述するＣＰＲ方向ベクトルを回転させる操作、つまりＣＰＲ面及びＣＰＲ画像を回転させる操作を受け付け可能である。
【0022】
ディスプレイ１３０は、例えばＬＣＤを含んでよく、各種情報を表示する。各種情報は、ボリュームデータから得られる３次元画像や２次元画像を含んでよい。３次元画像は、ボリュームレンダリング画像、サーフェスレンダリング画像、仮想内視鏡画像、仮想超音波画像、又はＣＰＲ画像等を含んでもよい。ボリュームレンダリング画像は、レイサム（RaySum）画像、ＭＩＰ画像、ＭｉｎＩＰ画像、平均値画像、又はレイキャスト画像等を含んでもよい。２次元画像は、アキシャル画像、サジタル画像、コロナル画像、ＭＰ画像、又はＣＰＲ画像等を含んでよい。
【0023】
メモリ１５０は、各種ＲＯＭやＲＡＭの一次記憶装置を含む。メモリ１５０は、ＨＤＤやＳＳＤの二次記憶装置を含んでもよい。メモリ１５０は、ＵＳＢメモリやＳＤカードの三次記憶装置を含んでもよい。メモリ１５０は、各種情報やプログラムを記憶する。各種情報は、ポート１１０により取得されたボリュームデータ、プロセッサ１４０により生成された画像、プロセッサ１４０により設定された設定情報、各種プログラムを含んでもよい。メモリ１５０は、プログラムが記録される非一過性の記録媒体の一例である。
【0024】
プロセッサ１４０は、ＣＰＵ、ＤＳＰ、又はＧＰＵ等を含んでもよい。プロセッサ１４０は、各種集積回路（例えばＬＳＩ又はＦＰＧＡ）により構成されてもよい。プロセッサ１４０は、メモリ１５０に記憶された医用画像処理プログラムを実行することにより、各種処理や制御を行う処理部１６０として機能する。
【0025】
図２は、処理部１６０の機能構成例を示すブロック図である。
【0026】
処理部１６０は、領域処理部１６１、パス取得部１６２、ＣＰＲ処理部１６３、画像生成部１６４、及び表示制御部１６５を備える。処理部１６０は、医用画像処理装置１００の各部を統括する。処理部１６０は、例えば、ＣＰＲを用いたボリュームデータの少なくとも一部の可視化に関する処理を行う。なお、処理部１６０に含まれる各部は、１つのハードウェアにより異なる機能として実現されてもよいし、複数のハードウェアにより異なる機能として実現されてもよい。また、処理部１６０に含まれる各部は、専用のハードウェア部品により実現されてもよい。
【0027】
領域処理部１６１は、例えばポート１１０を介して、被検体のボリュームデータを取得する。領域処理部１６１は、ボリュームデータに含まれる任意の領域を抽出する。領域処理部１６１は、例えばボリュームデータのボクセル値に基づいて、自動で関心領域を指定し、関心領域を抽出してよい。領域処理部１６１は、例えばＵＩ１２０を介して、手動で関心領域を指定し、関心領域を抽出してよい。本実施形態では、関心領域は、例えば管状組織（例えば血管、腸、又は気管支）を含む。
【0028】
パス取得部１６２は、ボリュームデータに含まれる３次元空間上のパスＰＨ、つまり３次元のパスＰＨの走行に関する情報を取得する。パスＰＨは、被検体における管状組織（例えば血管、腸、又は気管支）の中心線であってよい。パスＰＨは、被検体の脊椎や大腿骨に沿った線であってよい。パスＰＨは、例えばＵＩ１２０を介してユーザが任意に設定した線であってよい。パスＰＨは、ボリュームデータ内で様々な方向に延びていることがあり、様々な方向に曲がっていることがある。パスＰＨに沿うように、ＣＰＲ面が決定される。したがって、ＣＰＲ面には、パスＰＨにおける各点が含まれる。ＣＰＲ面は、曲面であり、パスＰＨと、後述するＣＰＲ方向ベクトルＶＴによって一意に定まる。ＣＰＲ面は、パスＰＨ上の各点を通過し任意のＣＰＲ方向ベクトルＶＴと平行な直線の集合によって定義される曲面である。パスＰＨ上の各点を通過する各ＣＰＲ方向ベクトルＶＴの方向は同一である。各ＣＰＲ方向ベクトルＶＴの方向が変化すると、異なるＣＰＲ面が得られる。よって、各ＣＰＲ方向ベクトルＶＴの方向が回転すると、得られるＣＰＲ面が回転する。
【0029】
ＣＰＲ処理部１６３は、パスＰＨとともにＣＰＲ面を規定するＣＰＲ方向ベクトルＶＴを複数仮定する。この場合、ＣＰＲ処理部１６３は、ＵＩ１２０を介して手動で複数の任意の方向のＣＰＲ方向ベクトルＶＴを仮定してもよいし、演算により自動的に複数の任意の方向のＣＰＲ方向ベクトルＶＴを仮定してもよい。また、ＣＰＲ処理部１６３は、ＵＩ１２０を介してＣＰＲ方向ベクトルＶＴを回転されることによって複数の任意の方向のＣＰＲ方向ベクトルＶＴを仮定してもよい。また、複数の任意のＣＰＲ方向ベクトルＶＴは、仮定できれば十分であり、実際に作成される必要は無い。ＣＰＲ処理部１６３は、仮定された複数のＣＰＲ方向ベクトルＶＴのうち、このＣＰＲ方向ベクトルＶＴにより規定されるＣＰＲ面を可視化したＣＰＲ画像において二重描画がされることを極力抑制した１つのＣＰＲ方向ベクトルＶＴを、推奨方向ベクトルＲＶとして決定する。なお、推奨方向ベクトルＲＶの候補として該当する候補が複数存在する場合、複数のうちいずれか１つの推奨方向ベクトルが決定されてよい。
【0030】
ここで、ＣＰＲ画像において二重描画が生じることは、ＣＰＲ面上に３次元上の同一の点が生じることを示す。また、３次元上の同一の点は、ＣＰＲ面上に直線として生じるので、ＣＰＲ画像においての二重描画は、折り返しが発生したかのようにみえる。つまり、ＣＰＲ処理部１６３は、仮定された複数のＣＰＲ方向ベクトルＶＴのうち、ＣＰＲ面上に３次元上の同一の点が生じることを極力抑制した１つのＣＰＲ方向ベクトルＶＴを、推奨方向ベクトルＲＶとして決定する。なお、推奨方向ベクトルＲＶの決定例の詳細については後述する。ＣＰＲ処理部１６３は、パスＰＨと推奨方向ベクトルＲＶとに基づいて、ＣＰＲ面を決定し、ＣＰＲ面を生成する。また、ＣＰＲ面を生成する前にＣＰＲ面上に３次元上の同一の点が生じることの判定方法についても後述する。
【0031】
画像生成部１６４は、各種画像を生成する。画像生成部１６４は、取得されたボリュームデータの少なくとも一部（例えば抽出された領域のボリュームデータ）に基づいて、３次元画像や２次元画像や断層画像を生成する。画像生成部１６４は、各種レンダリング（例えばボリュームレンダリング又はサーフィスレンダリング）を行って、画像を生成してよい。画像生成部１６４は、例えば、ＣＰＲ面上に位置するボリュームデータ（ボクセル）に基づいて、決定されたＣＰＲ面を可視化したＣＰＲ画像を生成する。
【0032】
表示制御部１６５は、各種データ、情報、又は画像をディスプレイ１３０に表示させる。画像は、被検体内の組織の一部を表現した画像であり、例えば、画像生成部１６４で生成された画像、所定の断面（例えばＣＰＲ面又はＭＰＲ面）の断面画像（例えばＣＰＲ画像又はＭＰＲ画像）、所定の断層の断層画像を含んでよい。
【0033】
次に、パスＰＨとＣＰＲ方向ベクトルＶＴとＣＰＲ面ＳＣの関係について説明する。
図３は、パスＰＨとＣＰＲ方向ベクトルＶＴとＣＰＲ面ＳＣの関係を説明するための図である。
図３は、パスＰＨとＣＰＲ方向ベクトルＶＴとＣＰＲ面ＳＣの関係を説明するための図である。
【0034】
ＣＰＲ処理部１６３は、ボリュームデータＶＤにパスＰＨが含まれており、取得（抽出）されたパスに対して、複数の任意の方向のＣＰＲ方向ベクトルＶＴを仮定する。ＣＰＲ処理部１６３は、仮定されたＣＰＲ方向ベクトルＶＴに垂直な投影面ＳＰを仮定する。仮定されたＣＰＲ方向ベクトルＶＴは複数存在するので、ＣＰＲ方向ベクトルＶＴ毎に投影面ＳＰが存在し、つまり投影面ＳＰが複数存在する。ＣＰＲ処理部１６３は、仮定された投影面ＳＰ毎に、３次元のパスＰＨを投影することで投影パスＰＨＰを生成する。よって、投影パスＰＨＰは二次元平面上のパスである。この投影パスＰＨＰはＣＰＲ方向ベクトルＶＴによって形状が異なる。ＣＰＲ処理部１６３は、投影パスＰＨＰが交差ないし二重線となるような重複の度合い（重複度合い）に基づいて、投影パスＰＨＰの重複度合いが比較的小さくなるＣＰＲ方向ベクトルＶＴを、推奨方向ベクトルＲＶとして決定する。画像生成部１６４は、パスＰＨと推奨方向ベクトルＲＶとにより定まるＣＰＲ面ＳＣ上にパスＰＨを描画し、ＣＰＲ画像ＧＣを生成する。図３では、図３のように仮定されたＣＰＲ方向ベクトルＶＴが推奨方向ベクトルＲＶであったことを例示している。投影パスＰＨＰに重複があるとＣＰＲ画像において二重描画が生じるからである。
【0035】
次に、ＣＰＲ軸ＡＣとＣＰＲ方向ベクトルＶＴの回転とについて説明する。
図４は、ＣＰＲ軸ＡＣとＣＰＲ方向ベクトルＶＴの回転とを説明するための図である。ＣＰＲ軸ＡＣは、ＣＰＲ方向ベクトルＶＴの回転中心となる軸である。
図４は、ＣＰＲ軸ＡＣとＣＰＲ方向ベクトルＶＴの回転とを説明するための図である。ＣＰＲ軸ＡＣは、ＣＰＲ方向ベクトルＶＴの回転中心となる軸である。
【0036】
ＣＰＲ方向ベクトルＶＴの回転の自由度は、２自由度である。具体的には、図４では、ＣＰＲ方向ベクトルＶＴの方向を規定するためのθ及びφの角度が、任意に調整可能である。ＣＰＲ方向ベクトルＶＴの回転の２自由度は、１自由度に制限され、例えばθの角度のみ調整可能であってもよい。これにより、ユーザがＵＩ１２０を介して手動でＣＰＲ軸ＡＣに対してＣＰＲ方向ベクトルＶＴを回転させる場合に、操作が容易になる。ＣＰＲ処理部１６３は、ＵＩ１２０を介して手動で２自由度又は１自由度を有するＣＰＲ方向ベクトルＶＴの方向（回転角度、向き）を指示してもよいし、演算によりＣＰＲ方向ベクトルＶＴの方向を決定してもよい。ＣＰＲ処理部１６３は、ＣＰＲ軸ＡＣを回転中心の軸として、ＣＰＲ方向ベクトルＶＴを回転させてＣＰＲ方向ベクトルＶＴの方向を変更可能である。
【0037】
なお、φは、例えば９０°であるが、他の角度であってもよい。φ＝９０°である場合には、θが変更されると、ＣＰＲ方向ベクトルは平面ＨＭ０上で回転する。φ≠９０°である場合には、ＣＰＲ軸ＡＣに対してＣＰＲ方向ベクトルＶＴが傾斜しており、φの角度が固定されたままθが変更されると、ＣＰＲ方向ベクトルＶＴは円錐面上を回転する。
【0038】
ＣＰＲ軸ＡＣは、パスＰＨの両端点を通る直線であってよい。また、ＣＰＲ軸ＡＣは、所定の基準方向ベクトルＳＶに沿った直線であってよい。基準方向ベクトルは、例えば、（０，０，１）、（０，１，０）、又は（０，０，１）等で規定される方向（向き）のベクトルであってよい。（０，０，１）と（０，１，０）と（０，０，１）とは、それぞれ単位ベクトルである。また、ＣＰＲ軸ＡＣは、パスＰＨ上の所定の注目点における接線としての直線であってもよい。パス上の注目点は、例えばＵＩ１２０を介して指定されてよいし、演算により自動で指定されてもよい。
【0039】
図５は、基準方向ベクトルＳＶの一例を示す図である。
【0040】
基準方向ベクトルＳＶは、例えば、ボリュームデータの座標系（Ｘ，Ｙ，Ｚ）に対して、３つの座標軸のそれぞれに沿う（１，０，０）、（０，１，０）、又は（０，０，１）の単位ベクトルで示されてよい。また、基準方向ベクトルＳＶは、ＣＰＲ軸ＡＣに沿うベクトルと、３つの座標軸に沿う（０，０，１）、（０，１，０）、及び（０，０，１）のうちのいずれかの単位ベクトルと、の外積で得られるベクトルで示されてもよい。言い換えると、基準方向ベクトルＳＶは、ＣＰＲ軸ＡＣに対して垂直な任意の方向のベクトルのいずれかであってもよい。また、基準方向ベクトルＳＶは、パスＰＨからみた縦隔方向であってもよい。
【0041】
図６は、任意の臓器ＺＫに含まれる血管をパスＰＨとしたＣＰＲ面の一例を示す模式図である。
【0042】
図６では、パス取得部１６２は、ＵＩ１２０を介してユーザが観察を希望する血管を指定することで、この血管のパスＰＨを取得している。ＣＰＲ処理部１６３は、血管の中心線をパスＰＨとし、パスＰＨと推奨方向ベクトルＲＶとに基づいてＣＰＲ面ＳＣを決定している。図６では、パスＰＨの曲がり具合に沿ってＣＰＲ面ＳＣが曲がっていることが分かる。画像生成部１６４は、決定されたＣＰＲ面ＳＣ上のボクセルの情報（例えばボクセル値）を基に描画したＣＰＲ画像を生成する。
【0043】
図７は、血管をパスＰＨとしたＣＰＲ面の一例を示す画像例を示す図である。
図７では、臓器（例えば心臓）にパスＰＨ（例えば血管）が３次元空間上の様々な方向に曲がっている。領域αでは、ＣＰＲ面ＳＣにおいてパスＰＨ上の同じ点が複数回出現されているので、ＣＰＲ画像ＧＣではパスＰＨ上の同じ点が複数回描画されることになる。図７のＣＰＲ面ＳＣは、例えば、基準方向ベクトルＳＶに基づいて決定されたＣＰＲ面（基準ＣＰＲ面）である。
図７では、臓器（例えば心臓）にパスＰＨ（例えば血管）が３次元空間上の様々な方向に曲がっている。領域αでは、ＣＰＲ面ＳＣにおいてパスＰＨ上の同じ点が複数回出現されているので、ＣＰＲ画像ＧＣではパスＰＨ上の同じ点が複数回描画されることになる。図７のＣＰＲ面ＳＣは、例えば、基準方向ベクトルＳＶに基づいて決定されたＣＰＲ面（基準ＣＰＲ面）である。
【0044】
次に、投影パスＰＨＰの重複度合いの評価例について説明する。前述のように、投影パスＰＨＰに重複があるとＣＰＲ画像に二重描画が生じるからである。
【0045】
ＣＰＲ処理部１６３は、各ＣＰＲ方向ベクトルＶＴの方向に垂直な平面に所定のパスＰＨを投影して、投影パスＰＨＰを生成する。投影パスＰＨＰ上の点からＣＰＲ方向ベクトルＶＴの方向の線上にあるボクセルがＣＰＲ画像に描画されるので、投影パスＰＨＰで重複が生じると、重複した点からＣＰＲ方向ベクトルＶＴの方向の線上のボクセルがＣＰＲ画像上で二重描画されることとなる。よって、投影パスＰＨＰの重複度合いは、投影パスＰＨＰ内の各点が複数回描画される頻度の高さを示す。ＣＰＲ処理部１６３は、投影パスＰＨＰの重複度合いを評価し、評価結果に基づいて推奨方向ベクトルＲＶを決定する。この場合、ＣＰＲ処理部１６３は、投影パスＰＨＰの重複度合いを評価する評価関数に基づいて、推奨方向ベクトルＲＶを決定してよい。評価関数は、例えば、投影パスＰＨＰの投影された面上で面積を用いた評価関数や、投影パスＰＨＰ上の各点間の距離を用いた評価関数が考えられる。
【0046】
ＣＰＲ処理部１６３は、推奨方向ベクトルＲＶとして、仮定された（ＵＩ１２０によって指定可能な）複数のＣＰＲ方向ベクトルＶＴのうち、投影パスＰＨＰの重複度合いが所定の閾値ＴＨ１よりも小さいＣＰＲ方向ベクトルＶＴ、つまりＣＰＲ画像において二重描画の少ないＣＰＲ方向ベクトルＶＴを選択してよい。閾値ＴＨ１は、例えばＣＰＲ画像において二重描画がされる範囲が観察対象の管状組織の径より小さくなるような値である。また、閾値ＴＨ１は、例えばＣＰＲ画像において二重描画がされる範囲がパスＰＨの最小半径より小さくなるような値である。一例として、ＣＰＲ処理部１６３は、仮定された複数のＣＰＲ方向ベクトルＶＴのうち、重複度合いが最小である（評価値が最大である）と評価されたＣＰＲ方向ベクトルＶＴを、推奨方向ベクトルＲＶとして決定してよい。また、ＣＰＲ処理部１６３は、仮定された複数のＣＰＲ方向ベクトルＶＴの評価を全て行い、得られた全ての重複度合いの平均値よりも重複度合いが小さい（全ての評価値の平均値よりも評価値が大きい）ＣＰＲ方向ベクトルＶＴを、推奨方向ベクトルＲＶとして決定してよい。これにより、医用画像処理装置１００は、二重描画が極力抑制されたＣＰＲ面ＳＣを決定可能であり、ＣＰＲ画像の描画精度を向上できる。
【0047】
また、ＣＰＲ処理部１６３は、基準方向ベクトルＳＶ（例えば（１，０，０）の単位ベクトル）に垂直な投影面ＳＰ（基準投影面ＳＰＳともいう）に対して、基準方向ベクトルＳＶに沿ってパスＰＨを投影することで得られる投影パスＰＨＰ（基準投影パスＰＨＳともいう）が、基準投影面ＳＰＳ上で重複する重複度合い（基準重複度合いともいう）を評価する。ＣＰＲ処理部１６３は、推奨方向ベクトルＲＶとして、仮定された複数のＣＰＲ方向ベクトルＶＴのうち、投影パスＰＨＰが投影面ＳＰ上で重複する重複度合いが、評価された基準重複度合いよりも小さい推奨方向ベクトルＲＶを決定してもよい。これにより、医用画像処理装置１００は、座標軸に沿う画一的な基準方向ベクトルＳＶの場合よりも二重描画が少ないＣＰＲ面ＳＣを決定可能であり、ＣＰＲ画像の描画精度を向上できる。なお、基準重複度合いは、予め導出されてメモリ１５０に保持されており、ＣＰＲ処理部１６３が必要時にメモリ１５０から取得するようにしてもよい。
【0048】
ＣＰＲ方向ベクトルＶＴの評価例つまり推奨方向ベクトルＲＶの決定例について、以下の第１評価例～第３評価例を例示する。
【0049】
（第１評価例）
第１評価例では、ＣＰＲ処理部１６３は、投影面ＳＰに投影された投影パスＰＨＰの太さを太くして、拡張投影パスＰＨＥを生成する。投影パスＰＨＰは、ＣＰＲ方向ベクトルＶＴに沿う方向を視線方向として、仮想二次元平面としての投影面ＳＰにパスＰＨを投影して得られるパスである。第１評価例では、ＣＰＲ処理部１６３は、投影面ＳＰにおける拡張投影パスＰＨＥが存在する領域であるパス領域ＲＰの面積に基づいて、重複度合いを算出する。例えば、拡張投影パスＰＨＥ内で同一の点が複数回描画された重複領域ＴＲの面積が大きい程、拡張投影パスＰＨＥの面積が小さくなる。この場合、重複度合いが大きくなる。例えば、拡張投影パスＰＨＥ内で重複領域ＴＲの面積が小さい程、拡張投影パスＰＨＥの面積が大きくなる。この場合、重複度合いが小さくなる。また、投影パスＰＨＰが投影パスＰＨＰ内で交差する点を有していない場合でも、投影パスＰＨＰ内の複数の点の距離が所定距離以下である場合、つまり投影パスＰＨＰが交差せずに接近している場合でも、拡張投影パスＰＨＥ内に重複領域ＴＲが発生することがある。この場合には、重複度合いが大きくなる。
第１評価例では、ＣＰＲ処理部１６３は、投影面ＳＰに投影された投影パスＰＨＰの太さを太くして、拡張投影パスＰＨＥを生成する。投影パスＰＨＰは、ＣＰＲ方向ベクトルＶＴに沿う方向を視線方向として、仮想二次元平面としての投影面ＳＰにパスＰＨを投影して得られるパスである。第１評価例では、ＣＰＲ処理部１６３は、投影面ＳＰにおける拡張投影パスＰＨＥが存在する領域であるパス領域ＲＰの面積に基づいて、重複度合いを算出する。例えば、拡張投影パスＰＨＥ内で同一の点が複数回描画された重複領域ＴＲの面積が大きい程、拡張投影パスＰＨＥの面積が小さくなる。この場合、重複度合いが大きくなる。例えば、拡張投影パスＰＨＥ内で重複領域ＴＲの面積が小さい程、拡張投影パスＰＨＥの面積が大きくなる。この場合、重複度合いが小さくなる。また、投影パスＰＨＰが投影パスＰＨＰ内で交差する点を有していない場合でも、投影パスＰＨＰ内の複数の点の距離が所定距離以下である場合、つまり投影パスＰＨＰが交差せずに接近している場合でも、拡張投影パスＰＨＥ内に重複領域ＴＲが発生することがある。この場合には、重複度合いが大きくなる。
【0050】
図８Ａは、拡張投影パスＰＨＥが重複領域ＴＲを有しない場合の一例を示す図である。図８Ｂは、投影パスＰＨＰが交差を有し、拡張投影パスＰＨＥが重複領域ＴＲを有する場合の一例を示す図である。図８Ｃは、投影パスＰＨＰが交差を有さずに接近しており、拡張投影パスＰＨＥが重複領域ＴＲを有する場合の一例を示す図である。
【0051】
ＣＰＲ処理部１６３は、投影面ＳＰにおいて、投影パスＰＨＰに沿って半径ｒの円を動かした場合に、この円が通る領域の面積をＳとしてよい。つまり、この面積Ｓは、パス領域ＲＰの面積に相当する。また、ＣＰＲ処理部１６３は、投影パスＰＨＰと同じ長さの直線の線分に沿って半径ｒの円を動かした場合に、つまり重複領域ＴＲが無い状態で、この円が通る領域の面積をＳ＿０としてよい。つまり、この面積Ｓ＿０は、重複領域ＴＲが無いパス領域ＲＰの面積に相当する。また、投影パスＰＨＰに沿う半径ｒの円の集合は、拡張投影パスＰＨＥに相当する。拡張投影パスＰＨＥの長さが同じなら、拡張投影パスＰＨＥ内での各点の重なりが大きい程、面積Ｓは小さくなる。したがって、ＣＰＲ処理部１６３は、例えば、「Ｓ＿０－Ｓ」の値や（Ｓ＿０－Ｓ）／Ｓ＿０の値を重複度合いとして用いてよい。つまり、ＣＰＲ処理部１６３は、上記の２つの面積の差分や上記の２つの面積の比率に基づいて、重複度合いを評価し、この評価結果に基づいて推奨方向ベクトルＲＶを決定してよい。
【0052】
（第２評価例）
第２評価例では、ＣＰＲ処理部１６３は、投影面ＳＰに投影された投影パスＰＨＰ上の各点について、投影パスＰＨＰ上での各点間の距離と投影面ＳＰ上での各点間の直線距離とに基づいて、ＣＰＲ方向ベクトルＶＴを評価する。例えば、ＣＰＲ処理部１６３は、投影パスＰＨＰを辿ると相互にある程度離れてる（例えば閾値ＴＨ２以上の距離である）２点が、投影面ＳＰ上でのこの２点の直線距離が近い（例えば閾値ＴＨ３未満の距離である）場合、投影パスＰＨＰの重複度合いが大きいと判定する。一方、ＣＰＲ処理部１６３は、投影パスＰＨＰ上での距離を辿ると相互にある程度離れてる２点が、投影面ＳＰ上でのこの２点の直線距離でも離れている（例えば閾値ＴＨ３以上の距離である）場合、投影パスＰＨＰの重複度合いが小さいと判定する。
第２評価例では、ＣＰＲ処理部１６３は、投影面ＳＰに投影された投影パスＰＨＰ上の各点について、投影パスＰＨＰ上での各点間の距離と投影面ＳＰ上での各点間の直線距離とに基づいて、ＣＰＲ方向ベクトルＶＴを評価する。例えば、ＣＰＲ処理部１６３は、投影パスＰＨＰを辿ると相互にある程度離れてる（例えば閾値ＴＨ２以上の距離である）２点が、投影面ＳＰ上でのこの２点の直線距離が近い（例えば閾値ＴＨ３未満の距離である）場合、投影パスＰＨＰの重複度合いが大きいと判定する。一方、ＣＰＲ処理部１６３は、投影パスＰＨＰ上での距離を辿ると相互にある程度離れてる２点が、投影面ＳＰ上でのこの２点の直線距離でも離れている（例えば閾値ＴＨ３以上の距離である）場合、投影パスＰＨＰの重複度合いが小さいと判定する。
【0053】
ＣＰＲ処理部１６３は、例えば以下の（式１）で示される評価関数Ｆ１を用いて、投影パスＰＨＰ（式中では単に「Ｐ」ともいう）上の点ｐについて、評価関数Ｆ１を計算して、評価値Ｆ１（ｐ）（評価値Ｅ１ともいう）を得る。つまり、投影パスＰＨＰ上で点ｐと点ｑとを投影パスＰＨＰの全体又は一部において順次変更しながら、点ｐ毎に評価値Ｆ（ｐ）を算出する。
【0054】
【数1】
【0055】
なお、点ｑは、点ｐと同様に、投影パスＰＨＰに含まれる点である。｜｜ｐ－ｑ｜｜２は、点ｐと点ｑとのユークリッド距離つまり平面上（投影面ＳＰ）での直線距離である。ｄｉｓｔａｎｃｅｐ（ｐ，ｑ）は、投影パスＰＨＰに沿った、投影パスＰＨＰ上の点ｐと点ｑとの距離である。ｔ１は、投影パスＰＨＰ上である程度離れていることの基準を示す閾値ＴＨ２であり、例えば１ｍｍ又は２ｍｍである。つまり、ｔ１は、ｄｉｓｔａｎｃｅｐ（ｐ，ｑ）の下限値を示す。なお、例えば、（式１）において、更に、ｄｉｓｔａｎｃｅｐ（ｐ，ｑ）の上限値を示す閾値ＴＨ２５を設けて、ｄｉｓｔａｎｃｅｐ（ｐ，ｑ）の範囲を更に制限してもよい。
【0056】
つまり、評価関数Ｆ１（Ｆ（ｐ））は、投影パスＰＨＰ上の点ｐと点ｐ以外の点ｑとの直線距離について、最短の直線距離を評価値Ｅ１として返す。ＣＰＲ処理部１６３は、評価関数Ｆ１は、投影パスＰＨＰ上の点ｐ毎に、最短の直線距離（ｐｑ間の最短直線距離ともいう）を評価値Ｅ１として導出する。よって、評価値Ｅ１が大きい程、投影パスＰＨＰの重複度合いが小さく、ＣＰＲ方向ベクトルＶＴの評価が高く、推奨方向ベクトルＲＶとして決定され易くなる。評価値Ｅ１が小さい程、投影パスＰＨＰの重複度合いが大きく、ＣＰＲ方向ベクトルＶＴの評価が低く、推奨方向ベクトルＲＶとして決定され難くなる。
【0057】
また、ＣＰＲ処理部１６３は、例えば以下の（式２）を用いて、評価関数Ｆ２を計算する。
【0058】
【数2】
【0059】
なお、ｔ２は、距離に関する閾値であり、評価関数Ｆ１で得られた評価値Ｅ１としてのｐｑ間の最短直線距離と比較され、ｐｑ間の最短直線距離に応じた値を評価関数Ｆ２により得られる値（評価値Ｅ２）に影響を与えるかどうかを定める指標である。例えば、０＜ｔ２＜ｔ１を満たす。（式２）では、ｔ２を用いることで、各点ｐについて、ｔ２よりもｐｑ間の最短直線距離が小さい場合には、ｔ２とｐｑ間の最短直線距離との差分を蓄積しし、ｔ２よりもｐｑ間の最短直線距離が大きい場合には、ｔ２とｐｑ間の最短直線距離との差分を蓄積しない（０を蓄積する）。そして、ＣＰＲ処理部１６３は、蓄積された差分を加算することで、評価関数Ｆ２の評価値Ｅ１を得る。また、ｔ２が可変とされる場合には、例えば、投影面ＳＰ上で点ｐからどの程度近くに最短直線距離となる点ｑが存在する場合に、上記の差分を蓄積するかを調整できる。
【0060】
評価値Ｅ２が大きい程、重複度合いが小さく、ＣＰＲ方向ベクトルＶＴの評価が高く、推奨方向ベクトルＲＶとして決定され易くなる。評価値Ｅ２が小さい程、重複度合いが大きく、ＣＰＲ方向ベクトルＶＴの評価が低く、推奨方向ベクトルＲＶとして決定され難くなる。
【0061】
なお、評価関数Ｆ２は、（式２）に限られない。例えば、ＣＰＲ処理部１６３は、（式２）のような積分の式の代わりに、投影パスＰＨＰ上の点を所定の間隔でサンプリングし、離散的に計算することで重複度合いを算出してもよい。
【0062】
【0063】
例えば、評価関数Ｆ１（Ｆ（ｐ））の算出時つまり各点ｐでのｐｑ間の最短直線距離の算出時に、少なくとも１回、評価関数Ｆ１（Ｆ（ｐ））＝０となる点ｐと点ｑとの集合が、重複領域ＴＲとなる。したがって、ＣＰＲ処理部１６３は、評価関数Ｆ１又は評価関数Ｆ１及びＦ２を用いて、投影パスＰＨＰにおいて重複が発生している位置を識別できる。なお、ＣＰＲ処理部１６３は、点ｐについての評価値Ｅ１（ｐｑ間の最短直線距離）が０でなくても０に近い（所定閾値以下の）値であれば、この位置で投影パスＰＨＰにおいて重複が発生していると判定してもよい。
【0064】
図９Ａでは、評価関数Ｆ１に基づいて重複していると判定された点ｐ及び点ｑの組み合わせが存在しておらず、重複領域ＴＲが存在しないことを例示している。図９Ｂでは、評価関数Ｆ１に基づいて重複していると判定された点ｐ及び点ｑの組み合わせが存在しており、重複領域ＴＲが含まれていることを示している。
【0065】
このような第２評価例によれば、点ｐ及び点ｑを用いて評価関数Ｆ１，Ｆ２を計算することで、投影パスＰＨＰの重複度合いを容易に算出可能である。なお、第１評価例の面積を用いた重複度合いの導出と比較すると、第２評価例の評価関数を用いた重複度合いの導出の方が、プロセッサ１４０が高速に演算できる。
【0066】
なお、第２評価例では、評価関数Ｆ１及び評価関数Ｆ２は、投影パスＰＨＰ（例えば血管）の太さに応じて変化するパラメータを追加してもよい。これにより、例えば、投影パスＰＨＰが太い血管である程、重複合いが大きく、投影パスＰＨＰが細い血管である程、重複度合いが小さくなるようにされてもよい。
【0067】
次に、医用画像処理装置１００の動作例について説明する。
図１０は、医用画像処理装置１００の動作例を示すフローチャートである。
図１０は、医用画像処理装置１００の動作例を示すフローチャートである。
【0068】
まず、ポート１１０は、被検体（例えば患者）のボリュームデータを取得する（Ｓ１１）。領域処理部１６１は、例えばＵＩ１２０を介して指定された関心領域（例えば血管や腸などの管状組織）を抽出する。パス取得部１６２は、抽出された管状組織を観察対象のパスＰＨとして取得する（Ｓ１２）。ＣＰＲ処理部１６３は、ＣＰＲ軸ＡＣを設定し、ＵＩ１２０によって回転できる（回転させ得る）ＣＰＲ方向ベクトルの範囲（例えば、ＣＰＲ軸ＡＣに対してφ＝９０°、θ＝０°～３６０°として）を設定する（Ｓ１３）。ＣＰＲ処理部１６３は、例えばＣＰＲ軸ＡＣに垂直な（つまりφ＝９０°として）ＣＰＲ方向ベクトルＶＴｉを３６０°分を作成（仮定）し、作成された各ＣＰＲ方向ベクトルＶＴｉを配列ＶＡに格納する（Ｓ１４）。配列ＶＡは、メモリ１５０に保持される。ＣＰＲ処理部１６３は、配列ＶＡに格納された各ＣＰＲ方向ベクトルＶＴｉについて、例えば評価関数Ｆ１，Ｆ２を用いて、各ＣＰＲ方向ベクトルＶＴｉに対応する各投影パスＰＨＰでの各重複量Ｄｉ（重複度合いの一例）を算出し、配列ＤＡに格納する（Ｓ１５）。配列ＤＡは、メモリ１５０に保持される。
【0069】
ＣＰＲ処理部１６３は、配列ＤＡに含まれる各重複量Ｄｉのうち、最小の重複量Ｄｍｉｎに対応するＣＰＲ方向ベクトルＶＴｍｉｎとしての推奨方向ベクトルＲＶを算出する（Ｓ１６）。ＣＰＲ処理部１６３は、パスＰＨと推奨方向ベクトルＲＶとを用いてＣＰＲ面ＳＣを決定する（Ｓ１７）。画像生成部１６４は、ＣＰＲ面上のボリュームデータ（ボクセル）によって描画してＣＰＲ画像ＧＣを生成する。表示制御部１６５は、ＣＰＲ画像ＧＣをディスプレイ１３０に表示させる（Ｓ１７）。また、表示制御部１６５は、推奨方向ベクトルＲＶを表現する情報（例えば推奨方向ベクトルＲＶの方向（角度）を示す情報）を表示させる（Ｓ１８）。ＣＰＲ処理部１６３は、ＵＩ１２０に推奨方向ベクトルＲＶを設定し、以降、ＵＩ１２０を介してＣＰＲ方向ベクトルを推奨方向ベクトルＲＶから回転させる操作を受け付ける（Ｓ１９）。
【0070】
このように、医用画像処理装置１００は、ユーザが任意にＣＰＲ方向ベクトルを選択するのではなく、投影パスＰＨＰ内での重複を加味して推奨方向ベクトルを自動的に決定できる。よって、客観的にユーザにとって観察し易いＣＰＲ画像を得ることができる。また、医用画像処理装置１００は、推奨方向ベクトルＲＶに基づいてがＣＰＲ画像を生成することで、ユーザがＵＩ１２０を介してＣＰＲ方向ベクトルの回転させる前に、一見して見易いＣＰＲ画像を生成できる。そして、ユーザは、一見して見易いＣＰＲ画像を確認した後、ＣＰＲ方向ベクトルを回転させてパスＰＨの表現が異なるＣＰＲ画像を順次確認でき、パスＰＨ等の観察対象の全体を把握できる。
【0071】
図１１は、ＣＰＲ画像の表示例を示す図である。
図１１では、ＣＰＲ画像ＧＣ１１，ＧＣ１２は、冠動脈のＲＣＡ（Ｒｉｇｈｔ Ｃｏｒｏｎａｒｙ Ａｒｔｅｒｙ）を示す画像である。ＣＰＲ画像ＧＣ１１内の二重描画領域ＤＲを示す情報は、パスＰＨ内での二重描画の発生、つまり基準投影パスＰＨＳでの重複の発生をユーザへ警告するための警告情報の一例である。
図１１では、ＣＰＲ画像ＧＣ１１，ＧＣ１２は、冠動脈のＲＣＡ（Ｒｉｇｈｔ Ｃｏｒｏｎａｒｙ Ａｒｔｅｒｙ）を示す画像である。ＣＰＲ画像ＧＣ１１内の二重描画領域ＤＲを示す情報は、パスＰＨ内での二重描画の発生、つまり基準投影パスＰＨＳでの重複の発生をユーザへ警告するための警告情報の一例である。
【0072】
ＣＰＲ画像ＧＣ１１は、所定の基準方向ベクトルＳＶ（θ＝０，φ＝９０°）を用いてＣＰＲ面ＳＣ（基準ＣＰＲ面ＳＳＣ）を決定し、この基準ＣＰＲ面ＳＳＣを可視化した画像である。この場合、ＣＰＲ処理部１６３は、基準方向ベクトルＳＶに垂直な基準投影面ＳＰＳに対して基準方向ベクトルＳＶに沿ってパスＰＨを投影することで、基準投影パスＰＨＳを得てよい。そして、ＣＰＲ処理部１６３は、基準投影パスＰＨＳの重複度合い（基準重複度合い）を算出し、メモリ１５０に保持しておく。
【0073】
図１１では、ＣＰＲ画像ＧＣ１１が表示された表示面がＸ座標とＹ座標との二次元平面であるとすると、Ｘ軸の中心の座標ｘ０とこの中心ｘ０からｘ軸負側に距離Ｌ１の位置（－ｘＬ１）とで囲まれた範囲α１、が折り返されて、Ｘ方向の中心座標ｘ０とこの中心座標ｘ０からｘ軸正側に距離Ｌ１の位置（＋ｘＬ１）とで囲まれた範囲α２、が出現している。つまり、二重描画領域ＤＲは、範囲α１と範囲α２とを含む。
【0074】
表示制御部１６５は、二重描画領域ＤＲを示す情報を表示させる。図１１では、二重描画領域ＤＲを示す情報として、範囲α１と範囲α２との境界線ＴＬ１と、二重描画領域ＤＲの境界線ＴＬ１と反対側の外縁を示す外縁線ＴＬ２と、Ｘ軸上での境界線ＴＬ１の位置を示すマークＭＫ１と、Ｘ軸上での外縁線ＴＬ２の位置を示すマークＭＫ２と、が表示されている。なお、二重描画領域ＤＲを示す情報はこれに限られない。なお、境界線ＴＬ１と外縁線ＴＬ２との表示が省略され、マークＭＫ１，ＭＫ２のみが表示されるようにしてもよい。二重描画領域ＤＲを示す情報が表示されることで、ユーザが二重描画領域ＤＲの位置を把握し易くなる。
【0075】
また、ＣＰＲ画像ＧＣ１２は、推奨方向ベクトルＲＶ（θ＝１１９°、φ＝９０°）を用いてＣＰＲ面ＳＣを決定し、このＣＰＲ面ＳＣを可視化した画像である。図１１では、ＣＰＲ画像ＧＣ１２において二重描画が発生していない。よって、ＣＰＲ画像ＧＣ１２は、ＣＰＲ画像ＧＣ１１と比較すると、ユーザがＣＰＲ面ＳＣの様子を把握し易い画像であると言える。
【0076】
なお、図１１では、処理部１６０は、基準方向ベクトルに対応するＣＰＲ画像ＧＣ１１において、二重描画領域ＤＲを示す情報を表示させているが、推奨方向ベクトルに対応するＣＰＲ画像ＧＣ１２において二重描画領域ＤＲが存在する場合には、同様に、二重描画領域ＤＲを示す情報を表示させてよい。二重描画領域ＤＲを示す情報の内容は、上述の内容と同様でよい。
【0077】
なお、推奨方向ベクトルＲＶが決定された後であっても、処理部１６０は、ＵＩ１２０を介してＣＰＲ方向ベクトルＶＴを回転させてＣＰＲ画像ＧＣを回転させてもよい。
【0078】
【0079】
また、ＣＰＲ画像ＧＣ１３は、推奨方向ベクトルＲＶ（θ＝－１４６°、φ＝９０°）を用いてＣＰＲ面ＳＣを決定し、このＣＰＲ面ＳＣを可視化した画像である。ＣＰＲ画像ＧＣ１３では、基準方向ベクトルに対応して得られたＣＰＲ画像と比較して、二重描画領域ＤＲの大きさが低減されているが、残存している状態である。ＣＰＲ画像ＧＣ１３では、図１１のＣＰＲ画像ＧＣ１１と同様に、二重描画領域ＤＲを示す情報が警告情報の一例として表示されている。警告情報は、パスＰＨ内での二重描画の発生、つまり投影パスＰＨＰでの重複の発生をユーザへ警告するための情報である。
【0080】
また、図１２では、ＣＰＲ処理部１６３は、パスＰＨの周辺では境界線ＴＬ１と外縁線ＴＬ２とを非表示とすることで、ユーザにとって、パスＰＨの近傍の組織を見易くしつつ、二重描画領域ＤＲの位置を容易に把握できる。
【0081】
（複数のＣＰＲ画像の同時表示）
次に、相互に向きが異なる推奨方向ベクトルＲＶを決定し、複数の推奨方向ベクトルＲＶに対応する複数のＣＰＲ画像を生成及び表示することについて説明する。
次に、相互に向きが異なる推奨方向ベクトルＲＶを決定し、複数の推奨方向ベクトルＲＶに対応する複数のＣＰＲ画像を生成及び表示することについて説明する。
【0082】
ＣＰＲ処理部１６３は、複数のＣＰＲ面ＳＣを決定するべく、同一のパスＰＨに対して相互に異なる方向の複数（例えば２つ）のＣＰＲ方向ベクトルＶＴを仮定する。ＣＰＲ処理部１６３は、例えば、θ＝θ１，φ＝９０°のＣＰＲ方向ベクトルＶＴ１と、θ＝θ１＋９０°，φ＝９０°のＣＰＲ方向ベクトルＶＴ２と、を仮定する。そして、θ１の角度を順次変更することで、複数のＣＰＲ方向ベクトルＶＴ１を仮定し、複数のＣＰＲ方向ベクトルＶＴ２を同時に仮定する。この場合、同時に仮定されるＣＰＲ方向ベクトルＶＴ１とＣＰＲ方向ベクトルＶＴ２とは、垂直である。ＣＰＲ処理部１６３は、ＣＰＲ方向ベクトルＶＴ１とＣＰＲ方向ベクトルＶＴ２とを９０°異なる状態でθ１を同じ角度ずつ回転させることで、ＣＰＲ方向ベクトルＶＴ１とＣＰＲ方向ベクトルＶＴ２との角度の関係を維持した状態で、角度θを変更していくことで２つのＣＰＲ方向ベクトルＶＴ１，ＶＴ２を回転させてよい。
【0083】
図１３は、パスＰＨと複数のＣＰＲ方向ベクトルＶＴ１，ＶＴ２と複数のＣＰＲ面ＳＣ１，ＳＣ２の関係を説明するための図である。
【0084】
ＣＰＲ処理部１６３は、仮定されたＣＰＲ方向ベクトルＶＴ１に垂直な投影面ＳＰ１を仮定する。仮定されたＣＰＲ方向ベクトルＶＴ１は複数存在するので、ＣＰＲ方向ベクトルＶＴ１毎に投影面ＳＰ１が存在し、つまり投影面ＳＰ１が複数存在する。ＣＰＲ処理部１６３は、仮定された投影面ＳＰ１毎に、３次元のパスＰＨを投影することで投影パスＰＨＰ１を生成する。よって、投影パスＰＨＰ１は二次元平面上のパスである。この投影パスＰＨＰ１はＣＰＲ方向ベクトルＶＴ１によって形状が異なる。ＣＰＲ処理部１６３は、投影パスＰＨＰ１が交差ないし二重線となるような重複の度合い（重複度合い）に基づいて、投影パスＰＨＰ１の重複度合いが比較的小さくなるＣＰＲ方向ベクトルＶＴ１を、推奨方向ベクトルＲＶ１として決定する。図１３では、図１３のように仮定されたＣＰＲ方向ベクトルＶＴ１が推奨方向ベクトルＲＶ１であったことを例示している。
【0085】
ＣＰＲ処理部１６３は、仮定されたＣＰＲ方向ベクトルＶＴ２に垂直な投影面ＳＰ２を仮定する。仮定されたＣＰＲ方向ベクトルＶＴ２は複数存在するので、ＣＰＲ方向ベクトルＶＴ２毎に投影面ＳＰ２が存在し、つまり投影面ＳＰ２が複数存在する。また、仮定された複数のＣＰＲ方向ベクトルＶＴ２は、仮定された複数のＣＰＲ方向ベクトルＶＴ１にそれぞれ対応する。ＣＰＲ処理部１６３は、仮定された投影面ＳＰ２毎に、３次元のパスＰＨを投影することで投影パスＰＨＰ２を生成する。よって、投影パスＰＨＰ２は二次元平面上のパスである。この投影パスＰＨＰ２はＣＰＲ方向ベクトルＶＴ２によって形状が異なる。ＣＰＲ処理部１６３は、投影パスＰＨＰ２が交差ないし二重線となるような重複の度合い（重複度合い）に基づいて、投影パスＰＨＰ２の重複度合いが比較的小さくなるＣＰＲ方向ベクトルＶＴ２を、推奨方向ベクトルＲＶ２として決定する。図１３では、図１３のように仮定されたＣＰＲ方向ベクトルＶＴ２が推奨方向ベクトルＲＶ２であったことを例示している。
【0086】
処理部１６０は、前述のように評価関数等を用いて、ＣＰＲ方向ベクトルＶＴ１に基づいて推奨方向ベクトルＲＶ１及びＣＰＲ面ＳＣ１を決定し、ＣＰＲ画像ＧＣ１を生成して表示させる。処理部１６０は、前述のように評価関数等を用いて、ＣＰＲ方向ベクトルＶＴ２に基づいて推奨方向ベクトルＲＶ２及びＣＰＲ面ＳＣ２を決定し、ＣＰＲ画像ＧＣ２を生成して表示させる。この際、投影パスＰＨＰ１の重複度合いと対応する投影パスＰＨＰ２の重複度合いとを統合した統合重複度合いＩＯに基づいて、推奨方向ベクトルＲＶ１と推奨方向ベクトルＲＶ２とを決定してよい。
【0087】
なお、複数のＣＰＲ方向ベクトルＶＴ（例えばＶＴ１，ＶＴ２）の成す角度（例えばθの角度）は、９０°でなくてもよい。この場合でも、ＣＰＲ処理部１６３は、ＣＰＲ方向ベクトルＶＴ１とＣＰＲ方向ベクトルＶＴ２とを９０°以外の角度（例えばθの角度）で異なる状態で同じ角度ずつ回転させることで、ＣＰＲ方向ベクトルＶＴ１とＣＰＲ方向ベクトルＶＴ２との角度の関係を維持した状態で、角度θを変更していくことで、２つのＣＰＲ方向ベクトルＶＴ１，ＶＴ２を回転させてよい。
【0088】
ユーザは、互いに例えば垂直なＣＰＲ面ＳＣ１のＣＰＲ画像ＧＣ１とＣＰＲ面ＳＣ２のＣＰＲ画像ＧＣ２とを同時に観察することがある。この２つの平面的な表示により三次元空間におけるパスＰＨの様子を想起し易いためである。この場合に、ＣＰＲ画像ＧＣ１とＣＰＲ画像ＧＣ２においていずれにも二重描画領域ＤＲが小さい又はないことが望ましい。
【0089】
ＣＰＲ処理部１６３は、仮定された２つのＣＰＲ方向ベクトルＶＴ１，ＶＴ２のそれぞれについて、例えば第１評価例又は第２評価例のように投影パスＰＨＰの占める面積又は評価関数Ｆ１，Ｆ２に基づいて、２つの推奨方向ベクトルＲＶ１，ＲＶ２のそれぞれを算出してよい。
【0090】
例えば、ＣＰＲ処理部１６３は、仮定されたＣＰＲ方向ベクトルＶＴ１に対応する投影面ＳＰ１に投影された投影パスＰＨＰ１の面積ＡＲ１と、仮定されたＣＰＲ方向ベクトルＶＴ２に対応する投影面ＳＰ２に投影された投影パスＰＨＰ２の面積ＡＲ２と、に基づいて、推奨方向ベクトルＲＶ１及び推奨方向ベクトルＲＶ２を生成してよい。この場合、ＣＰＲ処理部１６３は、例えば、面積ＡＲ１と面積ＡＲ２とを単純に合計し（統合の一例）、合計面積を得てよい。合計面積が小さい程、重複度合い（統合重複度合いＩＯの一例）が小さく、よって仮定されたＣＰＲ方向ベクトルＣＴ１，ＣＴ２が推奨方向ベクトルＲＶ１，ＲＶ２として決定され易くてよい。一方、合計面積が大きい程、重複度合い（統合重複度合いＩＯの一例）が大きく、よって仮定されたＣＰＲ方向ベクトルＣＴ１，ＣＴ２が推奨方向ベクトルＲＶ１，ＲＶ２として決定され難くてよい。なお、ここでは複数の面積を単純に合計することで面積を統合することを例示したが、これに限られない。例えば、処理部１６０は、複数の面積を乗算したり、平均したり、又は相乗平均したりすることで、複数の面積を統合してもよい。
【0091】
例えば、ＣＰＲ処理部１６３は、投影面ＳＰ１に投影された投影パスＰＨＰ１上の各点Ｐ及び各点Ｑを用いて、上記の評価関数Ｆ１，Ｆ２を用いて、評価関数Ｆ１の評価値Ｅ１１（評価値Ｅ１の一例）や評価関数Ｆ２の評価値Ｅ２１（評価値Ｅ２の一例）を算出してよい。ＣＰＲ処理部１６３は、投影面ＳＰ２に投影された投影パスＰＨＰ２上の各点Ｐ及び各点Ｑを用いて、上記の評価関数Ｆ１，Ｆ２を用いて、評価関数Ｆ１の評価値Ｅ１２（評価値Ｅ１の一例）や評価関数Ｆ２の評価値Ｅ２２（評価値Ｅ２の一例）を算出してよい。この場合、ＣＰＲ処理部１６３は、例えば、評価値Ｅ２１と評価値Ｅ２２とを単純に合計し（統合の一例）、合計評価値（統合重複度合いＩＯの一例）を得てよい。合計評価値が大きい程、２つの投影パスＰＨＰ１，ＰＨＰ２での重複を加味した統合重複度合いＩＯが小さく、よって仮定されたＣＰＲ方向ベクトルＣＴ１，ＣＴ２が推奨方向ベクトルＲＶ１，ＲＶ２として決定され易くてよい。一方、合計評価値が小さい程、統合重複度合いＩＯが大きく、よって仮定されたＣＰＲ方向ベクトルＣＴ１，ＣＴ２が推奨方向ベクトルＲＶ１，ＲＶ２として決定され難くてよい。なお、ここでは複数の評価値を単純に合計することで評価値を統合することを例示したが、これに限られない。例えば、処理部１６０は、複数の評価値を乗算したり、平均したり、又は相乗平均したりすることで、複数の評価値を統合してもよい。
【0092】
例えば、ＣＰＲ処理部１６３は、統合重複度合いＩＯ（例えば合計面積や合計評価値）が所定の閾値ＴＨ１１以下である場合、閾値ＴＨ１１以下の統合重複度合いＩＯが得られた際の投影パスＰＨＰ１の重複度合いが得られた際のＣＰＲ方向ベクトルＶＴ１を、推奨方向ベクトルＲＶ１に決定してよい。例えば、ＣＰＲ処理部１６３は、閾値ＴＨ１１以下である場合、閾値ＴＨ１１以下の統合重複度合いＩＯが得られた際の投影パスＰＨＰ２の重複度合いが得られた際のＣＰＲ方向ベクトルＶＴ２を、推奨方向ベクトルＲＶ２に決定してよい。
【0093】
【0094】
ＣＰＲ画像ＧＣ２１は、一の基準方向ベクトルＳＶ１（例えば、θ＝０，φ＝９０°つまり（１，０，０）で示される単位ベクトル）を用いてＣＰＲ面ＳＣ１（基準ＣＰＲ面ＳＳＣ１）を決定し、この基準ＣＰＲ面ＳＳＣ１を可視化した画像である。ＣＰＲ画像ＧＣ２２は、他の基準方向ベクトル（例えば、θ＝９０，φ＝９０°つまり（０，１，０）で示される単位ベクトル）を用いてＣＰＲ面ＳＣ２（基準ＣＰＲ面ＳＳＣ２）を決定し、この基準ＣＰＲ面ＳＳＣ２を可視化した画像である。よって、ＣＰＲ画像ＧＣ２１の基準ＣＰＲ面ＳＳＣ１とＣＰＲ画像ＧＣ２２の基準ＣＰＲ面ＳＳＣ２とは垂直である。また、ＣＰＲ画像ＧＣ２１には重複領域ＴＲ２１が存在し、ＣＰＲ画像ＧＣ２２には重複領域が存在しない。
【0095】
【0096】
ＣＰＲ画像ＧＣ３１は、推奨方向ベクトルＲＶ１（例えば、θ＝１２４°，φ＝９０°で示されるベクトル）を用いてＣＰＲ面ＳＣ１を決定し、このＣＰＲ面ＳＣ１を可視化した画像である。ＣＰＲ画像ＧＣ３２は、推奨方向ベクトルＲＶ２（例えば、θ＝－１４６°，φ＝９０°で示されるベクトル）を用いてＣＰＲ面ＳＣ２を決定し、このＣＰＲ面ＳＣ２を可視化した画像である。よって、ＣＰＲ画像ＧＣ３１のＣＰＲ面ＳＣ１とＣＰＲ画像ＧＣ３２のＣＰＲ面ＳＣ２とは垂直である。また、ＣＰＲ画像ＧＣ２１には重複領域が存在せず、ＣＰＲ画像ＧＣ２２には重複領域ＴＲ２２が存在する。しかし、重複領域ＴＲ２２は、重複領域ＴＲ２１よりも小さい。つまり、２つの推奨方向ベクトルＲＶ１，ＲＶ２を用いた２つのＣＰＲ面ＳＣ１，ＳＣ２が描画された２つのＣＰＲ画像ＧＣ３１，ＧＣ３２では、２つの基準方向ベクトルを用いた２つの基準ＣＰＲ面ＳＳＣ１，ＳＳＣ２が描画された２つのＣＰＲ画像ＧＣ２１，ＧＣ２２よりも、２つのＣＰＲ画像の重複領域ＴＲの合計の大きさが小さくなっている。
【0097】
なお、上述の説明では、統合重複度合いＩＯと閾値ＴＨ１１とを比較して、この比較結果に基づいて推奨方向ベクトルＲＶ１，ＲＶ２を決定することを例示したが、これに限られない。ＣＰＲ処理部１６３は、複数（例えば２つ）の異なる基準方向ベクトルＳＶを用いて、基準統合重複度合いＳＩＯと統合重複度合いＩＯとを比較し、この比較結果に基づいて推奨方向ベクトルＲＶ１，ＲＶ２を決定してもよい。ここでの統合重複度合いＩＯは、上述と同様に、仮定されたＣＰＲ方向ベクトルＶＴ１，ＶＴ２に基づく統合重複度合いＩＯでよい。基準統合重複度合いＳＩＯは、複数の基準方向ベクトルＳＶのそれぞれが複数の基準投影面ＳＰＳのそれぞれに投影された基準投影パスＰＨＳの重複度合いを統合した重複度合いでよい。
【0098】
例えば、ＣＰＲ処理部１６３は、基準方向ベクトルＳＶ１に垂直な投影面ＳＰ１（基準投影面ＳＰＳ１）に対して基準方向ベクトルＳＶ１に沿ってパスＰＨを投影することで、投影パスＰＨＰ１（基準投影パスＰＨＳ１）を得てよい。ＣＰＲ処理部１６３は、基準方向ベクトルＳＶ２に垂直な投影面ＳＰ２（基準投影面ＳＰＳ２）に対して基準方向ベクトルＳＶ２に沿ってパスＰＨを投影することで、投影パスＰＨＰ２（基準投影パスＰＨＳ２）を得てよい。そして、ＣＰＲ処理部１６３は、基準投影パスＰＨＳ１の重複度合い（第１基準重複度合い）と基準投影パスＰＨＳ２の重複度合い（第２基準重複度合い）とを算出し、第１基準重複度合いと第２基準重複度合いを統合した基準統合重複度合いＳＩＯを算出し、メモリ１５０に保持しておいてよい。基準統合重複度合いＳＩＯは、統合重複度合いＩＯの場合と同様に、基準となる合計面積や合計評価値であってよい。
【0099】
例えば、ＣＰＲ処理部１６３は、統合重複度合いＩＯが基準統合重複度合いＳＩＯよりも小さい場合、基準統合重複度合いＳＩＯよりも小さい統合重複度合いＩＯが得られた際の投影パスＰＨＰ１の重複度合いが得られた際のＣＰＲ方向ベクトルＶＴ１を、推奨方向ベクトルＲＶ１に決定してよい。例えば、ＣＰＲ処理部１６３は、統合重複度合いＩＯが基準統合重複度合いＳＩＯよりも小さい場合、基準統合重複度合いＳＩＯよりも小さい統合重複度合いＩＯが得られた際の投影パスＰＨＰ２の重複度合いが得られた際のＣＰＲ方向ベクトルＶＴ２を、推奨方向ベクトルＲＶ２に決定してよい。
【0100】
なお、推奨方向ベクトルＲＶ１，ＲＶ２が決定された後であっても、処理部１６０は、ＵＩ１２０を介してＣＰＲ方向ベクトルを回転させてＣＰＲ画像を回転させてもよい。この場合、処理部１６０は、ＵＩ１２０を介して、推奨方向ベクトルＲＶ１，ＲＶ２との成す角度を同じ角度で維持して、推奨方向ベクトルＲＶ１，ＲＶ２を回転させてよい。
【0101】
次に、本実施形態のバリエーションについて説明する。
【0102】
本実施形態では、ＣＰＲ画像ＧＣは、Ｓ－ＣＰＲ（Straight-CPR又はStretch-CPR）画像であってもよい。この場合、画像生成部１６４は、ＣＰＲ面ＳＣ上のパスＰＨを平面に直線的に引き延ばした状態の画像を生成する。Ｓ－ＣＰＲ画像を用いることで、ユーザが、例えば管状組織内のつまりや管状組織の太さが変化している箇所を把握し易くなる。
【0103】
本実施形態では、ＣＰＲ面ＳＣの厚さを考慮しないＣＰＲ面ＳＣを生成することを主に例示したが、これに限られない。ＣＰＲ処理部１６３は、ＣＰＲ面ＳＣの厚さを考慮した厚み付きＣＰＲ面を生成してもよい。よって、画像生成部１６４は、ＣＰＲ面に垂直な奥行き方向に厚みを持たせた範囲のボリュームデータを用いてＣＰＲ画像を生成してよい。この場合、ユーザは、厚み付きＣＰＲ面に対して垂直な方向に沿って移動する対象を確認しやすくなる。厚み方向において移動する観察対象の部位が、厚み付きＣＰＲ画像に移り込むためである。厚み付きＣＰＲ画像は、例えば、所定のボリュームレンダリング手法（例えばＭＩＰ、例えば平均値法）によって描画されてよい。また、厚み付きＣＰＲ画像では、厚みが厚い程、ＣＰＲ画像に二重描画が生じやすい。また、投影パスＰＨＰにおいては、厚みが厚い程、投影パスＰＨＰの太さが太い、閾値ＴＨ３が大きいとして表現でき、これによって同一の投影パスＰＨＰに対して重複度合いが大きく評価されやすい。
【0104】
これに対し、本実施形態の医用画像処理装置１００は、厚み付きＣＰＲ画像での二重描画の発生を極力低減しつつ、パスＰＨの様子を観察し易くできる。また、例えば第２評価例の評価関数Ｆ１又はＦ２において、重複が発生したとする投影パスＰＨＰの各点間の距離と比較される（式１）の閾値ｔ１や（式２）の閾値ｔ２が調整可能であってもよい。この調整は、処理部１６０により実施されてもよい。この場合、医用画像処理装置１００は、投影パスＰＨＰ内での重複の判定し易さを調整できる。
【0105】
本実施形態では、ＣＰＲ処理部１６３は、パスＰＨから作成されるＣＰＲ画像の全体を二重描画の度合いの評価対象とすること、つまり投影パスＰＨＰの全体を重複度合いの評価対象とすることを主に例示したが、これに限られない。ＣＰＲ処理部１６３は、パスＰＨの一部を二重描画の度合いの評価対象としてもよく、つまりパスＰＨから作成されるＣＰＲ画像の一部を二重描画の度合いの評価対象としてもよく、つまり投影パスＰＨＰの部分重複度合いを算出して評価してもよい。
【0106】
図１５は、パス部分ＰＨｐや投影パス部分ＰＨＰｐを説明するための図である。パスＰＨの全体とは、例えば、図３のパスＰＨにおける位置ｐ１から位置ｐ２までの範囲である。パスＰＨの一部とは、例えば、図３及び図１５のパスＰＨにおける位置ｐ１から位置ｐ３までの範囲や、パスＰＨにおける位置ｐ２から位置ｐ３までの範囲である。
【0107】
パスＰＨ上の場所によってはＣＰＲ画像での二重描画が診断に影響を与える度合いが異なる場合や、パスＰＨ上の場所によっては便宜的に表示していても診断に用いない場合が考えられる。この場合、ＣＰＲ処理部１６３は、パスＰＨの一部に対応する投影パスＰＨＰを複数の投影パス部分ＰＨＰｐに分割してよい。ＣＰＲ処理部１６３は、例えばＵＩ１２０を介して、複数の投影パス部分ＰＨＰｐのうちの重複度合い評価対象の投影パス部分ＰＨＰｐ１を指定してよい。図１５では、一例として、投影パス部分ＰＨＰｐ１の範囲は、投影パスＰＨＰにおける、パスＰＨ上の位置ｐ１に対応する位置ｐ１１と、パスＰＨ上の位置ｐ３に対応する位置ｐ１３と、の間の範囲である。部分重複度合いは、指定された投影パス部分ＰＨＰｐ１が投影面ＳＰ上で重複する度合であってよい。なお、投影パス部分ＰＨＰｐ１は、投影パスＰＨＰにおいて複数指定されてもよい。なお、投影パス部分ＰＨＰｐの重複度合いの評価方法は、上述したような投影パスＰＨＰ全体の評価（例えば面積を用いた評価、評価関数を用いた評価）と同様でよい。
【0108】
また、ＣＰＲ処理部１６３は、パスＰＨ自体を複数のパス部分ＰＨｐに分割して、重複度合い評価対象のパス部分ＰＨｐ１を指定して、パス部分ＰＨｐ１に対応する投影パス部分ＰＨＰｐ１を指定してもよい。この場合、ＣＰＲ処理部１６３は、例えばＵＩ１２０を介して、複数のパス部分ＰＨｐのうちの重複度合い評価対象のパス部分ＰＨｐ１を指定してよい。図１５では、一例として、パス部分ＰＨｐ１の範囲は、パスＰＨにおける位置ｐ１１と位置ｐ３との間の範囲である。ＣＰＲ処理部１６３は、指定されたパス部分ＰＨｐ１を投影面ＳＰに投影することで、パス部分ＰＨｐ１に対応する投影パス部分ＰＨＰｐ１を指定してもよい。
【0109】
例えば、ＣＰＲ処理部１６３は、投影パスＰＨＰに沿って動く半径ｒの円の移動可能範囲を、投影パス部分ＰＨＰｐ１に限定してよい。また、ＣＰＲ処理部１６３は、評価関数Ｆ１，Ｆ２における点ｐ及び点ｑの取りうる値の範囲を、投影パス部分ＰＨＰｐ１の範囲に限定してよい。
【0110】
例えば、分割して得られた投影パス部分ＰＨＰｐに応じて重複度合いの評価が異なることが予想される場合、ＣＰＲ処理部１６３は、所定の投影パス部分ＰＨＰｐを重複度合いの評価対象から除外し、この投影パス部分ＰＨＰｐに重複が発生することを許容してもよい。例えば、冠動脈のＬＣＸ(Left circumflex artery)の根本のように解剖学的に旋回することが知られている投影パス部分ＰＨＰｐは、予め重複度合いの評価対象から除外してもよい。つまり、ＣＰＲ処理部１６３は、投影パスＰＨＰにおける冠動脈のＬＣＸ(Left circumflex artery)の根本以外の部分を、重複度合いの評価対象の投影パス部分ＰＨＰｐ１として指定してよい。この場合、医用画像処理装置１００は、重複を抑制することが困難な投影パス部分ＰＨＰｐに対して過大なリソースを使用することを抑制できる。なお、冠動脈のＬＣＸの根本部分が重複度合いの評価対象から除外されてよいのは、冠動脈のＬＣＸ(Left circumflex artery)がＬＡＤ(left anterior descending artery)と根本部分が共通しているのでＬＡＤを観察すれば十分だからである。また、例えば、冠動脈を評価する場合は、末梢部分にはステント留置などの手術の対象とならない部位が存在するので、末梢部分を重複度合いの評価対象から除外しても良い。
【0111】
本実施形態では、評価関数Ｆ１及び評価関数Ｆ２の少なくとも一方は、同一の点が複数回描画される回数を重複数としてパラメータとして含んでよい。例えば、重複数が大きい程、重複度合いが大きいとしてもよい。一方、重複数が小さい程、重複度合いが小さいとしてもよい。
【0112】
本実施形態では、処理部１６０は、任意のタイミングで推奨方向ベクトルＲＶを決定してよい。例えば、医用画像処理装置１００の起動時やＣＰＲ画像を表示するアプリケーションの起動時に、初期値（初期設定）として、推奨方向ベクトルＲＶを決定してよい。この場合、初期の状態から、好適なＣＰＲ面ＳＣが決定され、ＣＰＲ画像が表示され得る。また、例えば、ＵＩ１２０が推奨方向ベクトルＲＶを決定するタイミングを指定するＵＩであるタイミング指定ＵＩ１２０Ａ（例えば最適化ボタン）を有していてよい。ＣＰＲ処理部１６３は、タイミング指定ＵＩ１２０Ａを介して推奨方向ベクトルＲＶを決定するタイミングの指定を受けた際に、上述した推奨方向ベクトルＲＶの決定を行ってよい。
【0113】
これにより、医用画像処理装置１００は、初期設定とは異なるユーザ所望のタイミングで、推奨方向ベクトルＲＶを決定し、推奨方向ベクトルＲＶに基づいてＣＰＲ画像を生成できる。よって、ユーザがＵＩ１２０を介してＣＰＲ方向ベクトルの回転させている最中でも、一見して見易いＣＰＲ画像にリセットされるようにＣＰＲ画像を生成できる。そして、ユーザは、一見して見易いＣＰＲ画像を確認した後、ＣＰＲ方向ベクトルを回転させてパスＰＨの表現が異なるＣＰＲ画像を順次確認でき、パスＰＨ等の観察対象の全体を把握できる。
【0114】
また、相互に異なる複数のＣＰＲ方向ベクトルＶＴを用いて複数のＣＰＲ画像を生成する場合、ＣＰＲ処理部１６３は、３つ以上のＣＰＲ方向ベクトルＶＴを決定し、３つ以上のＣＰＲ画像を生成し、同時に表示させてもよい。この場合、例えば、上記のθが６０°ずつ異なるようにして３つのＣＰＲ方向ベクトルＶＴ１１，ＶＴ１２，ＶＴ１３を決定してもよい。
【0115】
このように、本実施形態の医用画像処理装置１００は、容易に好適なＣＰＲ方向ベクトルＶＴつまり推奨方向ベクトルＲＶを決定でき、ユーザが見易いＣＰＲ画像を表示できる。よって、医用画像処理装置１００は、３次元空間でのＣＰＲ方向ベクトルＶＴの指定に不慣れなユーザであっても、ＣＰＲ方向ベクトルＶＴを適切に決定できる。また、医用画像処理装置１００は、演算により推奨方向ベクトルＲＶを決定するので、ユーザがＵＩ１２０を介して手動でユーザが見易いＣＰＲ画像を作成したときと比べて、再現性が高い推奨方向ベクトルＲＶを安定して得ることができる。よって、医用画像処理装置１００は、ユーザによってＣＰＲ方向ベクトルＶＴの自動決定時の方向が大きく異なることを抑制できる。
【0116】
また、医用画像処理装置１００は、任意のタイミングで、例えば初期設定の際に、又はＵＩ１２０（例えばタイミング指定ＵＩ１２０Ａ）を介して推奨方向ベクトルＲＶの設定が指示された際に、推奨方向ベクトルＲＶを決定できる。また、医用画像処理装置１００は、推奨方向ベクトルＲＶの決定後であっても、ＵＩ１２０を介して推奨方向ベクトルＲＶを回転可能とすることで、ＣＰＲ画像を回転させながら観察対象の管状組織の様子を詳細に観察可能である。ＵＩ１２０による推奨方向ベクトルＲＶの回転では、例えば２自由度を有してよく、上述のθ及びφの角度を調整可能でよい。これにより、推奨方向ベクトルＲＶが自動決定されても、ユーザが注目したいパスＰＨ上の箇所を様々な方向から確認可能である。
【0117】
以上、図面を参照しながら各種の実施形態について説明したが、本開示はかかる例に限定されないことは言うまでもない。当業者であれば、特許請求の範囲に記載された範疇内において、各種の変更例または修正例に想到し得ることは明らかであり、それらについても当然に本開示の技術的範囲に属するものと了解される。
【0118】
また、医用画像処理装置１００は、少なくともプロセッサ１４０及びメモリ１５０を備えてよい。ポート１１０、ＵＩ１２０、及びディスプレイ１３０は、医用画像処理装置１００に対して外付けであってもよい。
【0119】
また、撮像されたＣＴ画像としてのボリュームデータは、ＣＴ装置２００から医用画像処理装置１００へ送信されることを例示した。この代わりに、ボリュームデータが一旦蓄積されるように、ネットワーク上のサーバ（例えば画像データサーバ（ＰＡＣＳ）（不図示））等へ送信され、保管されてもよい。この場合、必要時に医用画像処理装置１００のポート１１０が、ボリュームデータを、有線回線又は無線回線を介してサーバ等から取得してもよいし、任意の記憶媒体（不図示）を介して取得してもよい。
【0120】
また、撮像されたＣＴ画像としてのボリュームデータは、ＣＴ装置２００から医用画像処理装置１００へポート１１０を経由して送信されることを例示した。これは、実質的にＣＴ装置２００と医用画像処理装置１００とを併せて一製品として成立している場合も含まれるものとする。また、医用画像処理装置１００がＣＴ装置２００のコンソールとして扱われている場合も含む。
【0121】
また、ＣＴ装置２００により画像を撮像し、被検体内部の情報を含むボリュームデータを生成することを例示したが、他の装置により画像を撮像し、ボリュームデータを生成してもよい。他の装置は、ＭＲＩ（Magnetic Resonance Imaging）装置、ＰＥＴ（Positron Emission Tomography）装置、血管造影装置（Angiography装置）、又はその他のモダリティ装置を含む。また、ＰＥＴ装置は、他のモダリティ装置と組み合わせて用いられてもよい。
【0122】
また、医用画像処理装置１００における動作が規定された医用画像処理方法として表現可能である。また、コンピュータに医用画像処理方法の各ステップを実行させるためのプログラムとして表現可能である。
【0123】
（上記実施形態の概要）
以上により、本開示には少なくとも以下の事項が記載されている。なお、括弧内には、上記した実施形態において対応する構成要素等を例示しているが、これに限定されるものではない。
以上により、本開示には少なくとも以下の事項が記載されている。なお、括弧内には、上記した実施形態において対応する構成要素等を例示しているが、これに限定されるものではない。
【0124】
［項目１］
プロセッサ（プロセッサ１４０）とユーザインタフェース（ＵＩ１２０）とを備える医用画像処理装置（医用画像処理装置１００）であって、
前記ユーザインタフェースは、第１の方向ベクトル（ＣＰＲ方向ベクトルＶＴ，ＶＴ１）を回転させる回転操作を受け付け、
前記プロセッサは、
被検体のボリュームデータ（ボリュームデータＶＤ）を取得し、
前記ボリュームデータにおける３次元空間上のパス（パスＰＨ）を取得し、
前記第１の方向ベクトルを回転させる場合に、前記第１の方向ベクトルに垂直な第１の投影面（投影面ＳＰ）に対して前記第１の方向ベクトルに沿って前記パスを投影することで得られる第１の投影パス（投影パスＰＨＰ，ＰＨＰ１）が前記第１の投影面上で重複する第１の重複度合いが、所定の第１の閾値（閾値ＴＨ１）以下である第１の推奨方向ベクトル（推奨方向ベクトルＲＶ，ＲＶ１）を決定し、
前記パスと前記第１の推奨方向ベクトルとに基づいて、前記パス上の各点を通過し前記第１の推奨方向ベクトルと平行な直線の集合によって定義される第１の曲面（ＣＰＲ面ＳＣ，ＳＣ１）を決定し、
前記第１の曲面上の前記ボリュームデータに基づいて、前記第１の曲面を可視化した第１のＣＰＲ画像（ＣＰＲ画像ＧＣ，ＧＣ１２，ＧＣ３１）を生成する、
医用画像処理装置。
プロセッサ（プロセッサ１４０）とユーザインタフェース（ＵＩ１２０）とを備える医用画像処理装置（医用画像処理装置１００）であって、
前記ユーザインタフェースは、第１の方向ベクトル（ＣＰＲ方向ベクトルＶＴ，ＶＴ１）を回転させる回転操作を受け付け、
前記プロセッサは、
被検体のボリュームデータ（ボリュームデータＶＤ）を取得し、
前記ボリュームデータにおける３次元空間上のパス（パスＰＨ）を取得し、
前記第１の方向ベクトルを回転させる場合に、前記第１の方向ベクトルに垂直な第１の投影面（投影面ＳＰ）に対して前記第１の方向ベクトルに沿って前記パスを投影することで得られる第１の投影パス（投影パスＰＨＰ，ＰＨＰ１）が前記第１の投影面上で重複する第１の重複度合いが、所定の第１の閾値（閾値ＴＨ１）以下である第１の推奨方向ベクトル（推奨方向ベクトルＲＶ，ＲＶ１）を決定し、
前記パスと前記第１の推奨方向ベクトルとに基づいて、前記パス上の各点を通過し前記第１の推奨方向ベクトルと平行な直線の集合によって定義される第１の曲面（ＣＰＲ面ＳＣ，ＳＣ１）を決定し、
前記第１の曲面上の前記ボリュームデータに基づいて、前記第１の曲面を可視化した第１のＣＰＲ画像（ＣＰＲ画像ＧＣ，ＧＣ１２，ＧＣ３１）を生成する、
医用画像処理装置。
【0125】
これにより、医用画像処理装置は、第１の投影パスの重複度合いが少ない推奨方向ベクトルを決定できる。よって、医用画像処理装置は、様々な方向にＣＰＲ方向ベクトルを指定可能である場合でも、与えられたパスに対して好適なＣＰＲ方向ベクトルを推奨することによって、ＣＰＲ画像上でボリュームデータ内の同一の点が複数回描画されること（二重描画）を抑制できる。よって、ユーザがＣＰＲ面上の様子を把握し易くなる。また、これによって、一つのＣＰＲ方向ベクトルが推奨されるのでユーザに関わらず同一のＣＰＲ画像が作成可能となり、診療の客観性再現性が増す。
【0126】
［項目２］
プロセッサとユーザインタフェースとを備える医用画像処理装置であって、
前記ユーザインタフェースは、第１の方向ベクトルを回転させる回転操作を受け付け、
前記プロセッサは、
被検体のボリュームデータを取得し、
前記ボリュームデータにおける３次元空間上のパスを取得し、
前記第１の方向ベクトルを回転させる場合に、仮定された複数の第１の方向ベクトルのうち、前記第１の方向ベクトルに垂直な第１の投影面に対して前記第１の方向ベクトルに沿って前記パスを投影することで得られる第１の投影パスが前記第１の投影面上で重複する第１の重複度合いが、所定の第１の基準方向ベクトル（基準方向ベクトルＳＶ，ＳＶ１）に垂直な第１の基準投影面（基準投影面ＳＰＳ，ＳＰＳ１）に対して前記第１の基準方向ベクトルに沿って前記パスを投影することで得られる第１の基準投影パス（基準投影パスＰＨＳ，ＰＨＳ１）が前記第１の基準投影面上で重複する度合いである第１基準重複度合いよりも小さい、第１の推奨方向ベクトルを決定し、
前記パスと前記第１の推奨方向ベクトルとに基づいて、前記パス上の各点を通過し前記第１の推奨方向ベクトルと平行な直線の集合によって定義される第１の曲面を決定し、
前記第１の曲面上のボリュームデータに基づいて、前記第１の曲面を可視化した第１のＣＰＲ画像を生成する、
医用画像処理装置。
プロセッサとユーザインタフェースとを備える医用画像処理装置であって、
前記ユーザインタフェースは、第１の方向ベクトルを回転させる回転操作を受け付け、
前記プロセッサは、
被検体のボリュームデータを取得し、
前記ボリュームデータにおける３次元空間上のパスを取得し、
前記第１の方向ベクトルを回転させる場合に、仮定された複数の第１の方向ベクトルのうち、前記第１の方向ベクトルに垂直な第１の投影面に対して前記第１の方向ベクトルに沿って前記パスを投影することで得られる第１の投影パスが前記第１の投影面上で重複する第１の重複度合いが、所定の第１の基準方向ベクトル（基準方向ベクトルＳＶ，ＳＶ１）に垂直な第１の基準投影面（基準投影面ＳＰＳ，ＳＰＳ１）に対して前記第１の基準方向ベクトルに沿って前記パスを投影することで得られる第１の基準投影パス（基準投影パスＰＨＳ，ＰＨＳ１）が前記第１の基準投影面上で重複する度合いである第１基準重複度合いよりも小さい、第１の推奨方向ベクトルを決定し、
前記パスと前記第１の推奨方向ベクトルとに基づいて、前記パス上の各点を通過し前記第１の推奨方向ベクトルと平行な直線の集合によって定義される第１の曲面を決定し、
前記第１の曲面上のボリュームデータに基づいて、前記第１の曲面を可視化した第１のＣＰＲ画像を生成する、
医用画像処理装置。
【0127】
これにより、医用画像処理装置は、重複度合いが少ない推奨方向ベクトルを決定できる。よって、医用画像処理装置は、様々な方向にＣＰＲ方向ベクトルを指定可能である場合でも、与えられたパスに対して好適なＣＰＲ方向ベクトルを推奨することによって、ＣＰＲ画像上でボリュームデータ内の同一の点が複数回描画されること（二重描画）を抑制できる。よって、ユーザがＣＰＲ面上の様子を把握し易くなる。また、これによって、一つのＣＰＲ方向ベクトルが推奨されるのでユーザに関わらず同一のＣＰＲ画像が作成可能となり、診療の客観性再現性が増す。
【0128】
［項目３］
前記ユーザインタフェースは、第２の方向ベクトル（ＣＰＲ方向ベクトルＶＴ２）を回転させる回転操作を受け付け、
前記プロセッサは、
前記第１の方向ベクトルに対して所定の角度を維持して第２の方向ベクトルを回転させ、
前記第１の方向ベクトル及び前記第２の方向ベクトルを回転させる場合に、前記第１の重複度合いと、仮定された複数の第２の方向ベクトルのうち、前記第２の方向ベクトルに垂直な第２の投影面（投影面ＳＰ２）に対して前記第２の方向ベクトルに沿って前記パスを投影することで得られる第２の投影パス（投影パスＰＨＰ２）が前記第２の投影面上で重複する第２の重複度合いと、を統合した統合重複度合い（統合重複度合いＩＯ）が、所定の第２の閾値（閾値ＴＨ１１）以下である場合、前記第２の閾値以下の前記統合重複度合いが得られた前記第１の重複度合いが得られた前記第１の方向ベクトルを前記第１の推奨方向ベクトルに決定し、前記第２の閾値以下の前記統合重複度合いが得られた前記第２の重複度合いが得られた前記第２の方向ベクトルを第２の推奨方向ベクトル（推奨方向ベクトルＲＶ２）に決定し、
前記パスと前記第２の推奨方向ベクトルとに基づいて、前記パス上の各点を通過し前記第２の推奨方向ベクトルと平行な直線の集合によって定義される第２の曲面（ＣＰＲ面ＳＣ２）を決定し、
前記第２の曲面上のボリュームデータに基づいて、前記第２の曲面を可視化した第２のＣＰＲ画像（ＣＰＲ画像ＧＣ３２）を生成し、
前記第１のＣＰＲ画像と前記第２のＣＰＲ画像とを同時に表示させる、
項目１に記載の医用画像処理装置。
前記ユーザインタフェースは、第２の方向ベクトル（ＣＰＲ方向ベクトルＶＴ２）を回転させる回転操作を受け付け、
前記プロセッサは、
前記第１の方向ベクトルに対して所定の角度を維持して第２の方向ベクトルを回転させ、
前記第１の方向ベクトル及び前記第２の方向ベクトルを回転させる場合に、前記第１の重複度合いと、仮定された複数の第２の方向ベクトルのうち、前記第２の方向ベクトルに垂直な第２の投影面（投影面ＳＰ２）に対して前記第２の方向ベクトルに沿って前記パスを投影することで得られる第２の投影パス（投影パスＰＨＰ２）が前記第２の投影面上で重複する第２の重複度合いと、を統合した統合重複度合い（統合重複度合いＩＯ）が、所定の第２の閾値（閾値ＴＨ１１）以下である場合、前記第２の閾値以下の前記統合重複度合いが得られた前記第１の重複度合いが得られた前記第１の方向ベクトルを前記第１の推奨方向ベクトルに決定し、前記第２の閾値以下の前記統合重複度合いが得られた前記第２の重複度合いが得られた前記第２の方向ベクトルを第２の推奨方向ベクトル（推奨方向ベクトルＲＶ２）に決定し、
前記パスと前記第２の推奨方向ベクトルとに基づいて、前記パス上の各点を通過し前記第２の推奨方向ベクトルと平行な直線の集合によって定義される第２の曲面（ＣＰＲ面ＳＣ２）を決定し、
前記第２の曲面上のボリュームデータに基づいて、前記第２の曲面を可視化した第２のＣＰＲ画像（ＣＰＲ画像ＧＣ３２）を生成し、
前記第１のＣＰＲ画像と前記第２のＣＰＲ画像とを同時に表示させる、
項目１に記載の医用画像処理装置。
【0129】
これにより、医用画像処理装置は、２つの異なるＣＰＲ方向ベクトルに対応する２つの異なるＣＰＲ画像を表示する場合でも、第１の投影パスの第１の重複度合い及び第２の投影パスの第２の重複度合いを加味した統合重複度合いが極力少ない第１の推奨方向ベクトル及び第２の推奨方向ベクトルを決定できる。よって、医用画像処理装置は、様々な方向に第１のＣＰＲ方向ベクトル及び第１のＣＰＲ方向ベクトルを指定可能である場合でも、与えられたパスに対して好適なＣＰＲ方向ベクトルを推奨することによって、ＣＰＲ画像上でボリュームデータ内の同一の点が複数回描画されること（二重描画）を抑制できる。よって、ユーザがＣＰＲ面上の様子を把握し易くなる。また、これによって、２つのＣＰＲ方向ベクトルが推奨されるのでユーザに関わらず同一のＣＰＲ画像の組が作成可能となり、診療の客観性再現性が増す。
【0130】
［項目４］
前記ユーザインタフェースは、第２の方向ベクトルを回転させる回転操作を受け付け、
前記プロセッサは、
前記第１の方向ベクトルに対して所定の角度を維持して第２の方向ベクトルを回転させ、
前記第１の方向ベクトル及び前記第２の方向ベクトルを回転させる場合に、前記第１の重複度合いと、仮定された複数の第２の方向ベクトルのうち、前記第２の方向ベクトルに垂直な第２の投影面に対して前記第２の方向ベクトルに沿って前記パスを投影することで得られる第２の投影パスが前記第２の投影面上で重複する第２の重複度合いと、を統合した統合重複度合いが、前記第１の基準重複度合いと、所定の第２の基準方向ベクトル（基準方向ベクトルＳＶ２）に垂直な第２の基準投影面（基準投影面ＳＰＳ２）に対して前記第２の基準方向ベクトルに沿って前記パスを投影することで得られる第２の基準投影パス（基準投影パスＰＨＳ２）が前記第２の基準投影面上で重複する度合いである第２の基準重複度合いと、を統合した基準統合重複度合い（基準統合重複度合いＳＩＯ）よりも小さい場合、前記基準統合重複度合いよりも小さい前記統合重複度合いが得られた前記第１の重複度合いが得られた前記第１の方向ベクトルを前記第１の推奨方向ベクトルに決定し、前記基準統合重複度合いよりも小さい前記統合重複度合いが得られた前記第２の重複度合いが得られた前記第２の方向ベクトルを第２の推奨方向ベクトルに決定し、
前記パスと前記第２の推奨方向ベクトルとに基づいて、前記パス上の各点を通過し前記第２の推奨方向ベクトルと平行な直線の集合によって定義される第２の曲面を決定し、
前記第２の曲面上のボリュームデータに基づいて、前記第２の曲面を可視化した第２のＣＰＲ画像を生成し、
前記第１のＣＰＲ画像と前記第２のＣＰＲ画像とを同時に表示させる、
項目２に記載の医用画像処理装置。
前記ユーザインタフェースは、第２の方向ベクトルを回転させる回転操作を受け付け、
前記プロセッサは、
前記第１の方向ベクトルに対して所定の角度を維持して第２の方向ベクトルを回転させ、
前記第１の方向ベクトル及び前記第２の方向ベクトルを回転させる場合に、前記第１の重複度合いと、仮定された複数の第２の方向ベクトルのうち、前記第２の方向ベクトルに垂直な第２の投影面に対して前記第２の方向ベクトルに沿って前記パスを投影することで得られる第２の投影パスが前記第２の投影面上で重複する第２の重複度合いと、を統合した統合重複度合いが、前記第１の基準重複度合いと、所定の第２の基準方向ベクトル（基準方向ベクトルＳＶ２）に垂直な第２の基準投影面（基準投影面ＳＰＳ２）に対して前記第２の基準方向ベクトルに沿って前記パスを投影することで得られる第２の基準投影パス（基準投影パスＰＨＳ２）が前記第２の基準投影面上で重複する度合いである第２の基準重複度合いと、を統合した基準統合重複度合い（基準統合重複度合いＳＩＯ）よりも小さい場合、前記基準統合重複度合いよりも小さい前記統合重複度合いが得られた前記第１の重複度合いが得られた前記第１の方向ベクトルを前記第１の推奨方向ベクトルに決定し、前記基準統合重複度合いよりも小さい前記統合重複度合いが得られた前記第２の重複度合いが得られた前記第２の方向ベクトルを第２の推奨方向ベクトルに決定し、
前記パスと前記第２の推奨方向ベクトルとに基づいて、前記パス上の各点を通過し前記第２の推奨方向ベクトルと平行な直線の集合によって定義される第２の曲面を決定し、
前記第２の曲面上のボリュームデータに基づいて、前記第２の曲面を可視化した第２のＣＰＲ画像を生成し、
前記第１のＣＰＲ画像と前記第２のＣＰＲ画像とを同時に表示させる、
項目２に記載の医用画像処理装置。
【0131】
これにより、医用画像処理装置は、２つの異なるＣＰＲ方向ベクトルに対応する２つの異なるＣＰＲ画像を表示する場合でも、画一的な基準方向ベクトルが選択された場合の２つの基準重複度合いを加味した基準統合重複度合いよりも、第１の投影パスの第１の重複度合い及び第２の投影パスの第２の重複度合いを加味した統合重複度合いが極力少ない、第１の推奨方向ベクトル及び第２の推奨方向ベクトルを決定できる。よって、医用画像処理装置は、様々な方向に第１のＣＰＲ方向ベクトル及び第２のＣＰＲ方向ベクトルを指定可能である場合でも、与えられたパスに対して好適なＣＰＲ方向ベクトルを推奨することによって、ＣＰＲ画像上でボリュームデータ内の同一の点が複数回描画されること（二重描画）を抑制できる。よって、ユーザがＣＰＲ面上の様子を把握し易くなる。また、これによって、２つのＣＰＲ方向ベクトルが推奨されるのでユーザに関わらず同一のＣＰＲ画像の組が作成可能となり、診療の客観性再現性が増す。
【0132】
［項目５］
前記第１の重複度合いは、前記第１の投影面上で前記第１の投影パスが交差又は所定距離以内に接近する度合いである、
項目１から４のいずれか１項に記載の医用画像処理装置。
前記第１の重複度合いは、前記第１の投影面上で前記第１の投影パスが交差又は所定距離以内に接近する度合いである、
項目１から４のいずれか１項に記載の医用画像処理装置。
【0133】
これにより、医用画像処理装置は、ＣＰＲ画像上で二重描画がされるとみなせる状態にある場合に、第１の重複度合いの少ない第１の推奨方向ベクトルを決定できる。
【0134】
［項目６］
前記ユーザインタフェースは、前記第１の推奨方向ベクトルの決定のタイミングを指示する操作を受け付け、
前記プロセッサは、
前記ユーザインタフェースにより前記方第１の推奨向ベクトルの決定のタイミンングを指示する操作を受け付けた際、前記第１の推奨方向ベクトルを決定する、
項目１から５のいずれか１項に記載の医用画像処理装置。
前記ユーザインタフェースは、前記第１の推奨方向ベクトルの決定のタイミングを指示する操作を受け付け、
前記プロセッサは、
前記ユーザインタフェースにより前記方第１の推奨向ベクトルの決定のタイミンングを指示する操作を受け付けた際、前記第１の推奨方向ベクトルを決定する、
項目１から５のいずれか１項に記載の医用画像処理装置。
【0135】
これにより、医用画像処理装置は、第１の方向ベクトルを回転させる場合に、ユーザ所望のタイミングで第１の推奨方向ベクトルを決定できる。また、ユーザが適切なＣＰＲ方向ベクトルを発見できなかった場合に第１の推奨方向ベクトルを決定できる。
【0136】
［項目７］
前記プロセッサは、前記推奨方向ベクトルを前記第１の方向ベクトルの初期値として決定する、
項目１から６のいずれか１項に記載の医用画像処理装置。
前記プロセッサは、前記推奨方向ベクトルを前記第１の方向ベクトルの初期値として決定する、
項目１から６のいずれか１項に記載の医用画像処理装置。
【0137】
これにより、医用画像処理装置は、第１の方向ベクトルを回転させる場合に、ユーザ所望のタイミングで第１の推奨方向ベクトルを決定できる。よって、第１の方向ベクトルの回転操作に不慣れなユーザであっても、好適な第１の方向ベクトルである第１の推奨方向ベクトル（例えば関心領域の代表的な状態を観察し易い第１のＣＰＲ画像が得られる第１の推奨方向ベクトル）を起点として、様々な方向にＣＰＲ画像を回転できる。第１の推奨方向ベクトルから第１の方向ベクトルを回転させることで、ユーザは、関心領域の全体を観察し易くなる。また、これによって、ＣＰＲ方向ベクトルが推奨されるのでユーザに関わらず同一のＣＰＲ画像が作成可能となり、診療の客観性再現性が増す。
【0138】
［項目８］
前記プロセッサは、
前記第１のＣＰＲ画像において、前記投影面上で前記投影パスが交差、重複又は接近していると判定された領域であり、前記ボリュームデータ内の同一の点が複数回描画される領域である二重描画領域を表示デバイス（ディスプレイ１３０）に表示させる、
項目１から７のいずれか１項に記載の医用画像処理装置。
前記プロセッサは、
前記第１のＣＰＲ画像において、前記投影面上で前記投影パスが交差、重複又は接近していると判定された領域であり、前記ボリュームデータ内の同一の点が複数回描画される領域である二重描画領域を表示デバイス（ディスプレイ１３０）に表示させる、
項目１から７のいずれか１項に記載の医用画像処理装置。
【0139】
これにより、医用画像処理装置は、第１のＣＰＲ画像における二重描画領域の位置、大きさ又は範囲等を確認し易くできる。
【0140】
［項目９］
前記プロセッサは、
前記パスを複数のパス部分（パス部分ＰＨｐ）に分割し、
前記複数のパス部分のうちの第１のパス部分（パス部分ＰＨｐ１）を指定し、
前記第１の重複度合いは、前記第１のパス部分を投影することで得られる第１の投影パス部分（投影パス部分ＰＨＰｐ１）が前記第１の投影面上で重複する度合いである、
項目１から８のいずれか１項に記載の医用画像処理装置。
前記プロセッサは、
前記パスを複数のパス部分（パス部分ＰＨｐ）に分割し、
前記複数のパス部分のうちの第１のパス部分（パス部分ＰＨｐ１）を指定し、
前記第１の重複度合いは、前記第１のパス部分を投影することで得られる第１の投影パス部分（投影パス部分ＰＨＰｐ１）が前記第１の投影面上で重複する度合いである、
項目１から８のいずれか１項に記載の医用画像処理装置。
【0141】
これにより、医用画像処理装置は、第１の重複度合いを導出する投影パスの部分を制限することで、所定の投影パス部分（ここでは第１の投影パス部分）の重複の状態を重視した第１の重複度合いを得ることができる。
【0142】
［項目１０］
被検体のボリュームデータを取得するステップと、
前記ボリュームデータにおける３次元空間上のパスを取得するステップと、
第１の方向ベクトルを回転させる回転操作を受け付けるユーザインタフェースを介して前記第１の方向ベクトルを回転させる場合に、前記第１の方向ベクトルに垂直な第１の投影面に対して前記第１の方向ベクトルに沿って前記パスを投影することで得られる第１の投影パスが前記第１の投影面上で重複する第１の重複度合いが、所定の閾値以下である第１の推奨方向ベクトルを決定するステップと、
前記パスと前記第１の推奨方向ベクトルとに基づいて、前記パス上の各点を通過し前記第１の推奨方向ベクトルと平行な直線の集合によって定義される第１の曲面を決定し、
前記第１の曲面上の前記ボリュームデータに基づいて、前記第１の曲面を可視化した第１のＣＰＲ画像を生成するステップと、
を有する医用画像処理方法。
被検体のボリュームデータを取得するステップと、
前記ボリュームデータにおける３次元空間上のパスを取得するステップと、
第１の方向ベクトルを回転させる回転操作を受け付けるユーザインタフェースを介して前記第１の方向ベクトルを回転させる場合に、前記第１の方向ベクトルに垂直な第１の投影面に対して前記第１の方向ベクトルに沿って前記パスを投影することで得られる第１の投影パスが前記第１の投影面上で重複する第１の重複度合いが、所定の閾値以下である第１の推奨方向ベクトルを決定するステップと、
前記パスと前記第１の推奨方向ベクトルとに基づいて、前記パス上の各点を通過し前記第１の推奨方向ベクトルと平行な直線の集合によって定義される第１の曲面を決定し、
前記第１の曲面上の前記ボリュームデータに基づいて、前記第１の曲面を可視化した第１のＣＰＲ画像を生成するステップと、
を有する医用画像処理方法。
【0143】
これにより、医用画像処理方法は、項目１と同様の効果が得られる。
【0144】
［項目１１］
被検体のボリュームデータを取得するステップと、
前記ボリュームデータにおける３次元空間上のパスを取得するステップと、
第１の方向ベクトルを回転させる回転操作を受け付けるユーザインタフェースを介して前記第１の方向ベクトルを回転させる場合に、前記第１の方向ベクトルに垂直な第１の投影面に対して前記第１の方向ベクトルに沿って前記パスを投影することで得られる第１の投影パスが前記第１の投影面上で重複する第１の重複度合いが、所定の第１の基準方向ベクトルに垂直な第１の基準投影面に対して前記第１の基準方向ベクトルに沿って前記パスを投影することで得られる第１の基準投影パスが前記第１の基準投影面上で重複する度合いよりも小さい、第１の推奨方向ベクトルを決定するステップと、
前記パスと前記第１の推奨方向ベクトルとに基づいて、前記パス上の各点を通過し前記第１の推奨方向ベクトルと平行な直線の集合によって定義される第１の曲面を決定するステップと、
前記第１の曲面上のボリュームデータに基づいて、前記第１の曲面を可視化した第１のＣＰＲ画像を生成するステップと、
を有する医用画像処理方法。
被検体のボリュームデータを取得するステップと、
前記ボリュームデータにおける３次元空間上のパスを取得するステップと、
第１の方向ベクトルを回転させる回転操作を受け付けるユーザインタフェースを介して前記第１の方向ベクトルを回転させる場合に、前記第１の方向ベクトルに垂直な第１の投影面に対して前記第１の方向ベクトルに沿って前記パスを投影することで得られる第１の投影パスが前記第１の投影面上で重複する第１の重複度合いが、所定の第１の基準方向ベクトルに垂直な第１の基準投影面に対して前記第１の基準方向ベクトルに沿って前記パスを投影することで得られる第１の基準投影パスが前記第１の基準投影面上で重複する度合いよりも小さい、第１の推奨方向ベクトルを決定するステップと、
前記パスと前記第１の推奨方向ベクトルとに基づいて、前記パス上の各点を通過し前記第１の推奨方向ベクトルと平行な直線の集合によって定義される第１の曲面を決定するステップと、
前記第１の曲面上のボリュームデータに基づいて、前記第１の曲面を可視化した第１のＣＰＲ画像を生成するステップと、
を有する医用画像処理方法。
【0145】
これにより、医用画像処理方法は、項目２と同様の効果が得られる。
【0146】
［項目１２］
項目１０に記載の医用画像処理方法をコンピュータに実行させるための医用画像処理プログラム。
項目１０に記載の医用画像処理方法をコンピュータに実行させるための医用画像処理プログラム。
【0147】
これにより、医用画像処理プログラムは、項目１と同様の効果が得られる。
【0148】
［項目１３］
項目１１に記載の医用画像処理方法をコンピュータに実行させるための医用画像処理プログラム。
項目１１に記載の医用画像処理方法をコンピュータに実行させるための医用画像処理プログラム。
【0149】
これにより、医用画像処理プログラムは、項目２と同様の効果が得られる。
【産業上の利用可能性】
【0150】
本開示は、ＣＰＲ画像上での二重描画の発生を抑制でき、ＣＰＲ面上の様子を把握し易くできる医用画像処理装置、医用画像処理方法、及び医用画像処理プログラム等に有用である。
【符号の説明】
【0151】
１００ 医用画像処理装置
１１０ ポート
１２０ ユーザインタフェース（ＵＩ）
１３０ ディスプレイ
１４０ プロセッサ
１５０ メモリ
１６０ 処理部
１６１ 領域処理部
１６２ パス取得部
１６３ ＣＰＲ処理部
１６４ 画像生成部
１６５ 表示制御部
２００ ＣＴ装置
ＧＣ ＣＰＲ画像
ＰＨ パス
ＰＨｐ パス部分
ＰＨＥ 拡張投影パス
ＰＨＰ，ＰＨＰ１，ＰＨＰ２ 投影パス
ＰＨＰｐ，ＰＨＰｐ１ 投影パス部分
ＲＶ，ＲＶ１，ＲＶ２ 推奨方向ベクトル
ＳＣ，ＳＣ１，ＳＣ２ ＣＰＲ面
ＳＰ，ＳＰ１，ＳＰ２ 投影面
ＶＴ，ＶＴ１，ＶＴ２ ＣＰＲ方向ベクトル
１１０ ポート
１２０ ユーザインタフェース（ＵＩ）
１３０ ディスプレイ
１４０ プロセッサ
１５０ メモリ
１６０ 処理部
１６１ 領域処理部
１６２ パス取得部
１６３ ＣＰＲ処理部
１６４ 画像生成部
１６５ 表示制御部
２００ ＣＴ装置
ＧＣ ＣＰＲ画像
ＰＨ パス
ＰＨｐ パス部分
ＰＨＥ 拡張投影パス
ＰＨＰ，ＰＨＰ１，ＰＨＰ２ 投影パス
ＰＨＰｐ，ＰＨＰｐ１ 投影パス部分
ＲＶ，ＲＶ１，ＲＶ２ 推奨方向ベクトル
ＳＣ，ＳＣ１，ＳＣ２ ＣＰＲ面
ＳＰ，ＳＰ１，ＳＰ２ 投影面
ＶＴ，ＶＴ１，ＶＴ２ ＣＰＲ方向ベクトル
【図1】
【図2】
【図3】
【図4】
【図5】
【図6】
【図7】
【図8A】
【図8B】
【図8C】
【図9A】
【図9B】
【図10】
【図11】
【図12】
【図13】
【図14A】
【図14B】
【図15】
【図16】
【図17】