知財求人 - 知財ポータルサイト「IP Force」
特許7180123医用画像処理装置、医用画像処理方法、プログラム、及びデータ作成方法
(19)【発行国】日本国特許庁(JP)
(12)【公報種別】特許公報(B2)
(11)【特許番号】
(24)【登録日】2022-11-21
(45)【発行日】2022-11-30
(54)【発明の名称】医用画像処理装置、医用画像処理方法、プログラム、及びデータ作成方法
(51)【国際特許分類】
A61B 6/03 20060101AFI20221122BHJP
G06T 15/08 20110101ALI20221122BHJP
G06T 1/00 20060101ALI20221122BHJP
【ＦＩ】
A61B6/03 360G
G06T15/08
G06T1/00 290B
【請求項の数】
28
(21)【出願番号】P 2018103767
(22)【出願日】2018-05-30
(65)【公開番号】P2019205791
(43)【公開日】2019-12-05
【審査請求日】2021-03-23
(73)【特許権者】
【識別番号】000002897
【氏名又は名称】大日本印刷株式会社
(74)【代理人】
【識別番号】100096091
【弁理士】
【氏名又は名称】井上 誠一
(72)【発明者】
【氏名】茂出木 敏雄
【審査官】松岡 智也
(56)【参考文献】
【文献】特開２０１１－０１９６３３（ＪＰ，Ａ）
【文献】国際公開第２０１２／０７３７６９（ＷＯ，Ａ１）
【文献】特開２００６－１３６６１９（ＪＰ，Ａ）
(58)【調査した分野】(Int.Cl.，ＤＢ名)
Ａ６１Ｂ ６／００－６／１４
(57)【特許請求の範囲】
【請求項1】
複数の断層画像を取得する画像取得手段と、信号値と色値及び不透明度との対応関係を定義するカラーマップを取得するカラーマップ取得手段と、
前記カラーマップを参照することで、前記断層画像の各画素の信号値を、信号値に応じた色値及び不透明度に変換し、色値及び不透明度を保持するボクセルの集合であるボクセルデータを作成するボクセル作成手段と、
前記ボクセルデータに基づいてボリュームレンダリング像を生成するレンダリング手段と、
前記断層画像毎に、断層画像の被写体領域の外郭を所定の幾何形状で近似し、当該幾何形状の情報を取得する幾何情報取得手段と、を備え、
前記ボクセル作成手段は、前記断層画像の各画素の信号値に応じた前記不透明度を、当該断層画像に対応する幾何形状の情報及び当該断層画像のスライス位置の情報に基づいて補正する
医用画像処理装置。
【請求項2】
前記断層画像毎に、断層画像の各画素に対応する不透明度の補正倍率を、当該断層画像に対応する幾何形状の情報、及び当該断層画像のスライス位置の情報に基づいて算出する補正倍率算出手段と、を更に備え、前記ボクセル作成手段は、前記断層画像の各画素の信号値に応じた前記不透明度に対して、前記補正倍率算出手段により算出された補正倍率を乗算する
請求項１に記載の医用画像処理装置。
【請求項3】
前記補正倍率算出手段は、算出した補正倍率を格納する補正テーブルを作成するテーブル作成手段と、を更に備え、前記ボクセル作成手段は、前記断層画像の各画素の信号値に応じた前記不透明度に対して、前記補正テーブルを参照して得られる補正倍率を乗算する
請求項２に記載の医用画像処理装置。
【請求項4】
前記補正倍率算出手段は、前記補正倍率を、前記幾何形状の幾何中心から遠ざかるにつれ小さくし、かつ、スライス位置が中央から末端に位置するほど小さくする請求項２または請求項３に記載の医用画像処理装置。
【請求項5】
前記幾何形状は楕円であり、前記幾何情報取得手段は、前記幾何形状の情報として、楕円の中心座標、楕円の横方向のサイズ、及び楕円の縦方向のサイズを取得し、
前記スライス位置の情報は、スライス総数、及びスライス順位であり、
前記補正倍率算出手段は、前記中心座標、前記横方向のサイズ、前記縦方向のサイズ、前記スライス総数、及び前記スライス順位に基づいて、補正倍率を算出する
請求項２から請求項４のいずれかに記載の医用画像処理装置。
【請求項6】
前記補正倍率算出手段は、前記中心座標における前記補正倍率を最大とし、前記中心座標と前記横方向のサイズ及び前記縦方向のサイズの比が同一となる前記楕円内部の各同心楕円上の前記補正倍率を、各同心楕円の横方向のサイズ及び縦方向のサイズが大きくなるにつれ小さくし、かつ、前記スライス総数／２とスライス順位との差が大きくなるにつれ前記補正倍率を小さくする請求項５に記載の医用画像処理装置。
【請求項7】
前記補正倍率算出手段は、前記中心座標に所定のオフセットを加算して前記中心座標を補正し、前記横方向のサイズ及び前記縦方向のサイズに所定の横方向のサイズに対する倍率及び縦方向のサイズに対する倍率を乗算して前記横方向のサイズ及び前記縦方向のサイズを補正し、前記スライス総数に所定のスライス総数に対する倍率を乗算して補正し、前記スライス順位を所定のスライスオフセットを加算して補正し、補正された前記中心座標、前記横方向のサイズ、前記縦方向のサイズ、前記スライス総数、及び前記スライス順位に基づいて、前記補正倍率を算出する請求項５または請求項６に記載の医用画像処理装置。
【請求項8】
前記オフセット、前記横方向のサイズに対する倍率、前記縦方向のサイズに対する倍率、前記スライス総数に対する倍率、及び前記スライスオフセットをユーザに調整させる調整手段と、を更に備える請求項７に記載の医用画像処理装置。
【請求項9】
複数の断層画像を取得する画像取得手段と、信号値と色値及び不透明度との対応関係を定義するカラーマップを取得するカラーマップ取得手段と、
前記カラーマップを参照することで、前記断層画像の各画素の信号値を、信号値に応じた色値及び不透明度に変換し、色値及び不透明度を保持するボクセルの集合であるボクセルデータを作成するボクセル作成手段と、
前記ボクセルデータに基づいてボリュームレンダリング像を生成するレンダリング手段と、
前記断層画像毎に、断層画像の被写体領域の外郭を所定の幾何形状で近似し、当該幾何形状の情報を取得する幾何情報取得手段と、
前記断層画像毎に、断層画像の各画素に対応する不透明度の補正倍率を、当該断層画像に対応する幾何形状の情報に基づいて算出する補正倍率算出手段と、を備え、
前記補正倍率算出手段は、前記被写体領域の中心部から離れた骨領域や体表を透過させ、中心部に近い臓器を鮮明に可視化するように前記補正倍率を、前記幾何形状の幾何中心から遠ざかるにつれ小さくし、
前記ボクセル作成手段は、
前記断層画像の各画素の信号値に応じた前記不透明度を、当該断層画像に対応する幾何形状の情報に基づいて補正し、かつ前記断層画像の各画素の信号値に応じた前記不透明度に対して、前記補正倍率算出手段により算出された補正倍率を乗算する
医用画像処理装置。
【請求項10】
前記補正倍率算出手段は、算出した補正倍率を格納する補正テーブルを作成するテーブル作成手段と、を更に備え、前記ボクセル作成手段は、前記断層画像の各画素の信号値に応じた前記不透明度に対して、前記補正テーブルを参照して得られる補正倍率を乗算する
請求項９に記載の医用画像処理装置。
【請求項11】
前記幾何形状は楕円であり、前記幾何情報取得手段は、前記幾何形状の情報として、楕円の中心座標、楕円の横方向のサイズ、及び楕円の縦方向のサイズを取得し、
前記補正倍率算出手段は、前記中心座標、前記横方向のサイズ、及び前記縦方向のサイズに基づいて、前記補正倍率を算出する
請求項９または請求項１０に記載の医用画像処理装置。
【請求項12】
前記補正倍率算出手段は、前記中心座標における前記補正倍率を最大とし、前記中心座標と前記横方向のサイズ及び前記縦方向のサイズの比が同一となる前記楕円内部の各同心楕円上の前記補正倍率を、各同心楕円の横方向のサイズ及び縦方向のサイズが大きくなるにつれ小さくする請求項１１に記載の医用画像処理装置。
【請求項13】
前記補正倍率算出手段は、前記中心座標に所定のオフセットを加算して前記中心座標を補正し、前記横方向のサイズ及び前記縦方向のサイズに所定の横方向のサイズに対する倍率及び縦方向のサイズに対する倍率を乗算して前記横方向のサイズ及び前記縦方向のサイズを補正し、補正された前記中心座標、前記横方向のサイズ及び前記縦方向のサイズに基づいて、前記補正倍率を算出する請求項１１または請求項１２に記載の医用画像処理装置。
【請求項14】
前記オフセット、前記横方向のサイズに対する倍率、及び前記縦方向のサイズに対する倍率をユーザに調整させる調整手段と、を更に備える請求項１３に記載の医用画像処理装置。
【請求項15】
前記断層画像の幾何形状の近似精度を所定の基準で判断し、所定の基準に満たない前記断層画像の幾何形状の情報を、所定の基準を満たす他の前記断層画像の幾何形状の情報に基づいて補正する幾何情報補正手段と、を更に備える請求項１から請求項１４のいずれかに記載の医用画像処理装置。
【請求項16】
前記幾何形状は楕円であり、前記幾何情報取得手段は、前記幾何形状の情報として、楕円の中心座標、楕円の横方向のサイズ及び楕円の縦方向のサイズを取得し、
前記幾何情報補正手段は、楕円の横方向のサイズが画像の横方向のサイズと一致する前記断層画像における楕円の横方向のサイズを、楕円の横方向のサイズが画像の横方向のサイズより小さい他の前記断層画像における楕円の横方向のサイズ及び縦方向のサイズの比と同一となるように、補正する
請求項１５記載の医用画像処理装置。
【請求項17】
信号値と所定の閾値との差分値に応じて当該信号値に対応する不透明度を減衰させるように前記カラーマップを調整するカラーマップ調整手段と、を更に備え、前記ボクセル作成手段は、調整された前記カラーマップを参照して、ボクセルデータを作成する
請求項１から請求項１６のいずれかに記載の医用画像処理装置。
【請求項18】
代表的な信号値と色値及び不透明度との対応関係に基づいて、所定範囲の信号値と色値及び不透明度との対応関係を定義するカラーマップを作成するカラーマップ作成手段と、を更に備え、前記カラーマップ取得手段は、作成した前記カラーマップを取得する
請求項１から請求項１７のいずれかに記載の医用画像処理装置。
【請求項19】
前記レンダリング手段は、前記ボクセルデータを生成するボリュームレンダリング像に投影変換した座標系を視点座標系とすると、視点座標系において、前記ボリュームレンダリング像の各画素よりＺ軸方向に沿って、Ｚ軸の上限値より下限値に向けて、視点座標系のボクセル座標毎に座標変換を行って前記ボクセルデータより不透明度を取得しながら、不透明ボクセルを探索し、最初に見つかった不透明ボクセルの視点座標系におけるＺ座標を、前記ボリュームレンダリング像の画素毎に記録した探索制御マスクを作成する探索制御マスク作成手段と、
前記ボリュームレンダリング像の画素毎に、前記探索制御マスクからＺ座標を取得し、取得したＺ座標よりＺ軸の下限値に向けてＺ軸方向に沿って、所定の光強度をもつ仮想光線を照射する際、視点座標系のボクセル座標毎に座標変換を行って前記ボクセルデータより不透明度を取得し、不透明ボクセルが見つかった場合、当該ボクセル座標に対して座標変換を行って前記ボクセルデータより色値を取得し、当該ボクセルの不透明度に基づいて前記光強度を減衰させるとともに、当該ボクセルの不透明度及び色値並びに前記減衰させた光強度に基づいて累積輝度値を算出する処理を繰り返し、算出された累積輝度値に基づいて、前記ボリュームレンダリング像の当該画素に対応する画素値として与えるレイキャスティング手段と、
を備える請求項１から請求項１８のいずれかに記載の医用画像処理装置。
【請求項20】
前記探索制御マスク作成手段及び前記レイキャスティング手段は、前記座標変換を行って前記ボクセルデータより不透明度または色値を取得する際、所定の回転を定義した回転行列、ＸＹＺ軸各方向のオフセット値、ＸＹＺ軸方向の拡大又は縮小倍率、Ｚ軸方向の変倍率、注視点から視点までの距離を含む前記所定の座標変換のパラメータを取得し、
前記視点座標系のボクセルの整数値の座標を、前記パラメータに基づいて前記ボクセルデータの座標系に変換を行って、前記ボクセルデータの実数値の座標を算出し、
算出した実数値の座標の近傍の複数の整数値の座標に対応する前記ボクセルデータの複数のボクセルを特定し、
特定した複数のボクセルの不透明度または色値に基づいて前記ボクセルデータより取得される不透明度または色値として算出する
請求項１９に記載の医用画像処理装置。
【請求項21】
前記レンダリング手段は、前記ボクセルデータに基づいて３Ｄテクスチャを生成する３Ｄテクスチャ生成手段と、
前記３Ｄテクスチャに対して所定の座標変換を行って変換後３Ｄテクスチャを生成する座標変換手段と、
３次元空間のＸＹ座標面上の四角形をＺ軸方向に並べた積層四角形を設定する積層四角形設定手段と、
所定の視点からＺ軸方向に平行な視線上の前記四角形のＸＹ座標に対応する前記変換後３Ｄテクスチャのボクセルの色値を前記ボクセルの不透明度に基づいて前記視点から遠い四角形の順にアルファブレンディングして取得し、前記ボリュームレンダリング像の画素値として与える画素値算出手段と、
を備える請求項１から請求項１８のいずれかに記載の医用画像処理装置。
【請求項22】
前記座標変換手段は、所定の回転を定義した回転行列、視野角度、視点位置、クリッピング位置、ＸＹＺ軸各方向のオフセット値、ＸＹＺ軸方向の拡大又は縮小倍率、Ｚ軸方向の変倍率を含む所定の座標変換のパラメータを取得し、
前記３Ｄテクスチャに対して、前記取得したパタメータを用いた前記所定の座標変換を行って前記変換後３Ｄテクスチャを生成する
請求項２１に記載の医用画像処理装置。
【請求項23】
前記座標変換手段及び前記画素値算出手段は、ビデオカードに搭載されたＧＰＵ及びフレームメモリを用いて実行する請求項２２に記載の医用画像処理装置。
【請求項24】
コンピュータの制御部が、複数の断層画像を取得する画像取得ステップと、
信号値と色値及び不透明度との対応関係を定義するカラーマップを取得するカラーマップ取得ステップと、
前記カラーマップを参照することで、前記断層画像の各画素の信号値を、信号値に応じた色値及び不透明度に変換し、色値及び不透明度を保持するボクセルの集合であるボクセルデータを作成するボクセル作成ステップと、
前記ボクセルデータに基づいてボリュームレンダリング像を生成するレンダリングステップと、
前記断層画像毎に、断層画像の被写体領域の外郭を所定の幾何形状で近似し、当該幾何形状の情報を取得する幾何情報取得ステップと、を備え、
前記ボクセル作成ステップは、前記断層画像の各画素の信号値に応じた前記不透明度を、当該断層画像に対応する幾何形状の情報及び当該断層画像のスライス位置の情報に基づいて補正する
医用画像処理方法。
【請求項25】
コンピュータの制御部が、
複数の断層画像を取得する画像取得ステップと、
信号値と色値及び不透明度との対応関係を定義するカラーマップを取得するカラーマップ取得ステップと、
前記カラーマップを参照することで、前記断層画像の各画素の信号値を、信号値に応じた色値及び不透明度に変換し、色値及び不透明度を保持するボクセルの集合であるボクセルデータを作成するボクセル作成ステップと、
前記ボクセルデータに基づいてボリュームレンダリング像を生成するレンダリングステップと、
前記断層画像毎に、断層画像の被写体領域の外郭を所定の幾何形状で近似し、当該幾何形状の情報を取得する幾何情報取得ステップと、
前記断層画像毎に、断層画像の各画素に対応する不透明度の補正倍率を、当該断層画像に対応する幾何形状の情報に基づいて算出する補正倍率算出ステップと、を備え、
前記補正倍率算出ステップは、前記被写体領域の中心部から離れた骨領域や体表を透過させ、中心部に近い臓器を鮮明に可視化するように前記補正倍率を、前記幾何形状の幾何中心から遠ざかるにつれ小さくし、
前記ボクセル作成ステップは、
前記断層画像の各画素の信号値に応じた前記不透明度を、当該断層画像に対応する幾何形状の情報に基づいて補正し、かつ前記断層画像の各画素の信号値に応じた前記不透明度に対して、前記補正倍率算出ステップにより算出された補正倍率を乗算する
医用画像処理方法。
【請求項26】
コンピュータを、複数の断層画像を取得する画像取得手段、
信号値と色値及び不透明度との対応関係を定義するカラーマップを取得するカラーマップ取得手段、
前記カラーマップを参照することで、前記断層画像の各画素の信号値を、信号値に応じた色値及び不透明度に変換し、色値及び不透明度を保持するボクセルの集合であるボクセルデータを作成するボクセル作成手段、
前記ボクセルデータに基づいてボリュームレンダリング像を生成するレンダリング手段、
前記断層画像毎に、断層画像の被写体領域の外郭を所定の幾何形状で近似し、当該幾何形状の情報を取得する幾何情報取得手段、として機能させ、
前記ボクセル作成手段は、前記断層画像の各画素の信号値に応じた前記不透明度を、当該断層画像に対応する幾何形状の情報及び当該断層画像のスライス位置の情報に基づいて補正するプログラム。
【請求項27】
コンピュータを、
複数の断層画像を取得する画像取得手段と
信号値と色値及び不透明度との対応関係を定義するカラーマップを取得するカラーマップ取得手段、
前記カラーマップを参照することで、前記断層画像の各画素の信号値を、信号値に応じた色値及び不透明度に変換し、色値及び不透明度を保持するボクセルの集合であるボクセルデータを作成するボクセル作成手段、
前記ボクセルデータに基づいてボリュームレンダリング像を生成するレンダリング手段、
前記断層画像毎に、断層画像の被写体領域の外郭を所定の幾何形状で近似し、当該幾何形状の情報を取得する幾何情報取得手段、
前記断層画像毎に、断層画像の各画素に対応する不透明度の補正倍率を、当該断層画像に対応する幾何形状の情報に基づいて算出する補正倍率算出手段、として機能させ、
前記補正倍率算出手段は、前記被写体領域の中心部から離れた骨領域や体表を透過させ、中心部に近い臓器を鮮明に可視化するように前記補正倍率を、前記幾何形状の幾何中心から遠ざかるにつれ小さくし、
前記ボクセル作成手段は、
前記断層画像の各画素の信号値に応じた前記不透明度を、当該断層画像に対応する幾何形状の情報に基づいて補正し、かつ前記断層画像の各画素の信号値に応じた前記不透明度に対して、前記補正倍率算出手段により算出された補正倍率を乗算するプログラム。
【請求項28】
コンピュータの制御部が、複数の断層画像を取得する画像取得ステップと、
前記断層画像毎に、断層画像の被写体領域の外郭を所定の幾何形状で近似し、当該幾何形状の情報を取得する幾何情報取得ステップと、
前記断層画像の各画素の信号値に応じたオパシティカーブにより規定される不透明度を、当該断層画像に対応する幾何形状の情報及び当該断層画像のスライス位置の情報に基づいて補正する不透明度補正ステップと、
前記断層画像の各画素に対応する補正後の不透明度を格納したデータを作成するデータ作成ステップと、
を含むデータ作成方法。
【発明の詳細な説明】
【技術分野】
【0001】
本開示は、医用画像処理装置、医用画像処理方法、プログラム、及びデータ作成方法に関する。
【背景技術】
【0002】
医用画像診断において、特定の人体組織の観察に適したボリュームレンダリング像を得たい場面がある。例えば、胸部や頭部にある内臓や脳等の所定の人体組織の観察を行いたい場合に、内臓や脳等は骨に囲まれており、骨領域はむしろ診断の妨げになる。ＣＴ画像から生成されるボリュームレンダリング像では一般に骨が鮮明に描写され、内臓や血管は隠れてしまうことがあるため、可視化にあたっては切断を行うなど工夫が必要となることがある。骨領域は比較的高信号であるが、信号値分布が広範で造影血管と被るため、カラーマップを工夫して骨領域だけを透明化したオパシティカーブを設計することは容易ではなかった。
【0003】
たとえば、下記の特許文献１では、ＣＴ画像の信号値に基づいて領域の抽出とラベリングを行い、ラベリングされた各肋骨領域が適正か否かを、別途構築した統計上の肋骨の体軸方向の高さ幅情報を有する統計データベースを用いた手法が開示されている。これにより、骨領域（肋骨領域）を検出して除去（透明化）させることができるが、特許文献１の方法では、肋骨に先天性の奇形があったり、骨折を伴う場合には判定誤りが発生し、適切に除去できなくなる。
【0004】
また、領域拡張法（リージョングローイング法）を用いて非観察対象としたい骨領域を抽出した３次元マスクを作成し、３次元マスクを参照しながらマスク処理により非観察対象としたい骨領域を除去したボリュームレンダリング像を生成する方法が開示されている（特許文献２、３参照）。領域拡張法とは、非観察対象領域の画素を指定し、その画素を開始点（拡張開始点）として、近傍画素を次々と抽出する方法である。
【先行技術文献】
【特許文献】
【0005】
【文献】特開２００９－２４０５６９号公報
【文献】特許４０８７５１７号公報
【文献】特許５２５７９５８号公報
【発明の概要】
【発明が解決しようとする課題】
【0006】
しかしながら、領域拡張法を用いた３次元マスクの作成は、ユーザによる複数の拡張開始点の指示が必須で、拡張の打ち切り段階もユーザが指示する必要があるため、自動化が難しく、ユーザごとに結果にバラつきが発生するという問題がある。また領域拡張法に限らず、３次元マスクの作成は、ユーザのスキルや経験が必要であり、その作成に時間や手間がかかる。
【0007】
本開示は、前述した問題点に鑑みてなされたものであり、その目的とすることは、特定の領域を鮮明に可視化することが可能な、医用画像処理装置等を提供することである。
【課題を解決するための手段】
【0008】
本開示の一実施例によると、複数の断層画像を取得する画像取得手段と、信号値と色値及び不透明度との対応関係を定義するカラーマップを取得するカラーマップ取得手段と、前記カラーマップを参照することで、前記断層画像の各画素の信号値を、信号値に応じた色値及び不透明度に変換し、色値及び不透明度を保持するボクセルの集合であるボクセルデータを作成するボクセル作成手段と、前記ボクセルデータに基づいてボリュームレンダリング像を生成するレンダリング手段と、前記断層画像毎に、断層画像の被写体領域の外郭を所定の幾何形状で近似し、当該幾何形状の情報を取得する幾何情報取得手段と、を備え、前記ボクセル作成手段は、前記断層画像の各画素の信号値に応じた前記不透明度を、当該断層画像に対応する幾何形状の情報及び当該断層画像のスライス位置の情報に基づいて補正する医用画像処理装置が提供される。
【0009】
本開示の一実施形態によれば、カラーマップを参照して断層画像の各画素の信号値を色及び不透明度に変換してボクセルデータを作成する際、カラーマップ（オパシティカーブ）に定義された不透明度を、断層画像の被写体領域の外郭を近似した幾何形状の情報及び当該断層画像のスライス位置の情報に基づいて補正する。このように、被写体形状（幾何形状）及びスライス位置の情報に基づいて不透明度をコントロールすることで、容易に特定の領域を鮮明に可視化させたレンダリング像を得ることができる。
【0010】
また、前記断層画像毎に、断層画像の各画素に対応する不透明度の補正倍率を、当該断層画像に対応する幾何形状の情報、及び当該断層画像のスライス位置の情報に基づいて算出する補正倍率算出手段と、を更に備え、前記ボクセル作成手段は、前記断層画像の各画素の信号値に応じた前記不透明度に対して、前記補正倍率算出手段により算出された補正倍率を乗算するようにしても良い。これにより、断層画像の各画素の信号値に応じた不透明度に対して、幾何形状の情報及びスライス位置の情報に基づいて算出された補正倍率を乗算することで、不透明度を補正することができる。
【0011】
また、前記補正倍率算出手段は、算出した補正倍率を格納する補正テーブルを作成するテーブル作成手段と、を更に備え、前記ボクセル作成手段は、前記断層画像の各画素の信号値に応じた前記不透明度に対して、前記補正テーブルを参照して得られる補正倍率を乗算するようにしても良い。これにより、不透明度の補正倍率を格納した補正テーブルを参照して、不透明度を補正することができる。
【0012】
また、前記補正倍率算出手段は、前記補正倍率を、前記幾何形状の幾何中心（重心）から遠ざかるにつれ小さくし、かつ、スライス位置が中央から末端に位置するほど小さくするようにしても良い。これにより、被写体領域の中心から遠ざかるつれ不透明度が減衰するので、被写体領域の中心部から離れた骨領域や体表が透過され、被写体領域の中心部に近い臓器等が鮮明に可視化される。また、スライス位置が中央から遠ざかるつれ不透明度が減衰するので、特にスライス末端近くで被写体領域に占める骨領域（頭頂骨や顎骨）が増大する頭部ＣＴにおいて、骨領域を効果的に透過させることができる。
【0013】
また、前記幾何形状は楕円であり、前記幾何情報取得手段は、前記幾何形状の情報として、楕円の中心座標、楕円の横方向のサイズ、及び楕円の縦方向のサイズを取得し、前記スライス位置の情報は、スライス総数、及びスライス順位であり、前記補正倍率算出手段は、前記中心座標、前記横方向のサイズ、前記縦方向のサイズ、前記スライス総数、及び前記スライス順位に基づいて、前記補正倍率を算出するようにしても良い。これにより、被写体領域の外郭を楕円で近似し、楕円パラメータ（中心座標、横方向のサイズ、縦方向のサイズ）、及びスライスパラメータ（スライス総数、スライス順位）に基づいて、不透明度の補正倍率を算出することができる。なお本開示において、横方向とは、断層画像の座標系を規定する２つの座標軸方向（Ｘ軸方向、Ｙ軸方向）のうち、該断層画像に映る被写体の左右軸の方向と近いほうの座標軸方向のことを指す。また、縦方向とは、断層画像の座標系を規定する２つの座標軸方向（Ｘ軸方向、Ｙ軸方向）のうち、該断層画像に映る被写体の背腹軸（前後軸）の方向と近いほうの座標軸方向のことを指す。ここで、断層画像のＸ軸およびＹ軸は、例えば、モダリティでの撮影時に適宜設定されるが、本開示では、被写体の左右軸と方向が近いほうの座標軸をＸ軸、被写体の背腹軸（前後軸）と方向が近いほうの座標軸をＹ軸とする。すなわち、本開示では、横方向は断層画像のＸ軸方向に相当し、縦方向は断層画像のＹ軸方向に相当するものとして説明する。
【0014】
また、前記補正倍率算出手段は、前記中心座標における前記補正倍率を最大とし、前記中心座標と前記横方向のサイズ及び前記縦方向のサイズの比が同一となる前記楕円内部の各同心楕円上の前記補正倍率を、各同心楕円の横方向のサイズ及び縦方向のサイズが大きくなるにつれ小さくし、かつ、前記スライス総数／２とスライス順位との差が大きくなるにつれ前記補正倍率を小さくするようにしても良い。これにより、楕円内部の不透明度の補正倍率がスライス方向を考慮して好適に算出される。
【0015】
また、前記補正倍率算出手段は、前記中心座標に所定のオフセットを加算して前記中心座標を補正し、前記横方向のサイズ及び前記縦方向のサイズに所定の横方向のサイズに対する倍率及び縦方向のサイズに対する倍率を乗算して前記横方向のサイズ及び前記縦方向のサイズを補正し、前記スライス総数に所定のスライス総数に対する倍率を乗算して補正し、前記スライス順位を所定のスライスオフセットを加算して補正し、補正された前記中心座標、前記横方向のサイズ、前記縦方向のサイズ、前記スライス総数、及び前記スライス順位に基づいて、前記補正倍率を算出するようにしても良い。これにより、楕円パラメータ（中心座標、横方向のサイズ、縦方向のサイズ）、及びスライスパラメータ（スライス総数、スライス順位）を補正したうえで、不透明度の補正倍率を算出することができる。
【0016】
また、前記オフセット、前記横方向のサイズに対する倍率、及び前記縦方向のサイズに対する倍率、前記スライス総数に対する倍率、及び前記スライスオフセットをユーザに調整させる調整手段と、を更に備えるようにしても良い。これにより、楕円パラメータ（中心座標、横方向のサイズ、縦方向のサイズ）、及びスライスパラメータ（スライス総数、スライス順位）の補正量をユーザが調整することができる。
【0017】
また、本開示の一実施態様によると、複数の断層画像を取得する画像取得手段と、信号値と色値及び不透明度との対応関係を定義するカラーマップを取得するカラーマップ取得手段と、前記カラーマップを参照することで、前記断層画像の各画素の信号値を、信号値に応じた色値及び不透明度に変換し、色値及び不透明度を保持するボクセルの集合であるボクセルデータを作成するボクセル作成手段と、前記ボクセルデータに基づいてボリュームレンダリング像を生成するレンダリング手段と、前記断層画像毎に、断層画像の被写体領域の外郭を所定の幾何形状で近似し、当該幾何形状の情報を取得する幾何情報取得手段と、前記断層画像毎に、断層画像の各画素に対応する不透明度の補正倍率を、当該断層画像に対応する幾何形状の情報に基づいて算出する補正倍率算出手段と、を備え、前記補正倍率算出手段は、前記被写体領域の中心部から離れた骨領域や体表を透過させ、中心部に近い臓器を鮮明に可視化するように前記補正倍率を、前記幾何形状の幾何中心から遠ざかるにつれ小さくし、前記ボクセル作成手段は、前記断層画像の各画素の信号値に応じた前記不透明度を、当該断層画像に対応する幾何形状の情報に基づいて補正し、かつ前記断層画像の各画素の信号値に応じた前記不透明度に対して、前記補正倍率算出手段により算出された補正倍率を乗算する医用画像処理装置が提供される。
【0018】
本開示の一実施形態によれば、カラーマップを参照して断層画像の各画素の信号値を色及び不透明度に変換してボクセルデータを作成する際、カラーマップ（オパシティカーブ）に定義された不透明度を、断層画像の被写体領域の外郭を近似した幾何形状の情報に基づいて補正する。このように、被写体形状（幾何形状）に基づいて不透明度をコントロールすることで、容易に特定の領域を鮮明に可視化させたレンダリング像を得ることができる。また、補正倍率算出手段によって、前記断層画像毎に、断層画像の各画素に対応する不透明度の補正倍率を、当該断層画像に対応する幾何形状の情報に基づいて算出する際に、被写体領域の中心部から離れた骨領域や体表を透過させ、中心部に近い臓器を鮮明に可視化するように前記補正倍率を、前記幾何形状の幾何中心から遠ざかるにつれ小さくする。断層画像の各画素の信号値に応じた不透明度に対して、幾何形状の情報に基づいて算出された補正倍率を乗算することで、不透明度を補正することができる。また、被写体領域の中心から遠ざかるつれ不透明度が減衰するので、被写体領域の中心部から離れた骨領域や体表が透過され、被写体領域の中心部に近い臓器等が鮮明に可視化される。
【0019】
また、前記補正倍率算出手段は、算出した補正倍率を格納する補正テーブルを作成するテーブル作成手段と、を更に備え、前記ボクセル作成手段は、前記断層画像の各画素の信号値に応じた前記不透明度に対して、前記補正テーブルを参照して得られる補正倍率を乗算するようにしても良い。これにより、不透明度の補正倍率を格納した補正テーブルを参照して、不透明度を補正することができる。
【0021】
また、前記幾何形状は楕円であり、前記幾何情報取得手段は、前記幾何形状の情報として、楕円の中心座標、楕円の横方向のサイズ、及び楕円の縦方向のサイズを取得し、前記補正倍率算出手段は、前記中心座標、前記横方向のサイズ、及び前記縦方向のサイズに基づいて、補正倍率を算出するようにしても良い。これにより、楕円パラメータ（中心座標、横方向のサイズ、縦方向のサイズ）に基づいて、不透明度の補正倍率を算出することができる。
【0022】
また、前記補正倍率算出手段は、前記中心座標における前記補正倍率を最大とし、前記中心座標と前記横方向のサイズ及び前記縦方向のサイズの比が同一となる前記楕円内部の各同心楕円上の前記補正倍率を、各同心楕円の横方向のサイズ及び縦方向のサイズが大きくなるにつれ小さくしても良い。これにより、楕円内部の不透明度の補正倍率が好適に算出される。
【0023】
また、前記補正倍率算出手段は、前記中心座標に所定のオフセットを加算して前記中心座標を補正し、前記横方向のサイズ及び前記縦方向のサイズに所定の横方向のサイズに対する倍率及び縦方向のサイズに対する倍率を乗算して前記横方向のサイズ及び前記縦方向のサイズを補正し、補正された前記中心座標、前記横方向のサイズ、及び前記縦方向のサイズに基づいて、前記補正倍率を算出するようにしても良い。これにより、楕円パラメータ（中心座標、横方向のサイズ、縦方向のサイズ）を補正したうえで、不透明度の補正倍率を算出することができる。
【0024】
また、前記オフセット、前記横方向のサイズに対する倍率、前記縦方向のサイズに対する倍率をユーザに調整させる調整手段と、を更に備えるようにしても良い。これにより、楕円パラメータの補正量をユーザが調整することができる。
【0025】
また、前記断層画像の幾何形状の近似精度を所定の基準で判断し、所定の基準に満たない前記断層画像の幾何形状の情報を、所定の基準を満たす他の前記断層画像の幾何形状の情報に基づいて補正する幾何情報補正手段と、を更に備えるようにしても良い。これにより、近似精度が良好でない断層画像の幾何形状の情報（幾何パラメータ）を補正することができる。
【0026】
また、前記幾何形状は楕円であり、前記幾何情報取得手段は、前記幾何形状の情報として、楕円の中心座標、楕円の横方向のサイズ及び楕円の縦方向のサイズを取得し、前記幾何情報補正手段は、楕円の横方向のサイズが画像の横方向のサイズと一致する前記断層画像における楕円の横方向のサイズを、楕円の横方向のサイズが画像の横方向のサイズより小さい他の前記断層画像における楕円の横方向のサイズ及び縦方向のサイズの比と同一となるように、補正するようにしても良い。これにより、被写体領域が画像範囲に収まらない断層画像であっても、幾何パラメータ（楕円パラメータ）を良好に算出することができる。
【0027】
また、信号値と所定の閾値との差分値に応じて当該信号値に対応する不透明度を減衰させるように前記カラーマップを調整するカラーマップ調整手段と、を更に備え、前記ボクセル作成手段は、調整された前記カラーマップを参照して、ボクセルデータを作成するようにしても良い。これにより、カラーマップ（オパシティカーブ）の不透明度を信号値に応じて調整することができる。
【0028】
また、代表的な信号値と色値及び不透明度との対応関係に基づいて、所定範囲の信号値と色値及び不透明度との対応関係を定義するカラーマップを作成するカラーマップ作成手段と、を更に備え、前記カラーマップ取得手段は、作成した前記カラーマップを取得するようにしても良い。これにより、例えば、代表的な信号値（例えば、１０箇所程度の特徴的な信号値）と当該信号値に対する色及び不透明度をユーザが設定するだけで、自動でカラーマップを作成することができる。
【0029】
また、前記レンダリング手段は、前記ボクセルデータを生成するボリュームレンダリング像に投影変換した座標系を視点座標系とすると、視点座標系において、前記ボリュームレンダリング像の各画素よりＺ軸方向に沿って、Ｚ軸の上限値より下限値に向けて、視点座標系のボクセル座標毎に座標変換を行って前記ボクセルデータより不透明度を取得しながら、不透明ボクセルを探索し、最初に見つかった不透明ボクセルの視点座標系におけるＺ座標を、前記ボリュームレンダリング像の画素毎に記録した探索制御マスクを作成する探索制御マスク作成手段と、前記ボリュームレンダリング像の画素毎に、前記探索制御マスクからＺ座標を取得し、取得したＺ座標よりＺ軸の下限値に向けてＺ軸方向に沿って、所定の光強度をもつ仮想光線を照射する際、視点座標系のボクセル座標毎に座標変換を行って前記ボクセルデータより不透明度を取得し、不透明ボクセルが見つかった場合、当該ボクセル座標に対して座標変換を行って前記ボクセルデータより色値を取得し、当該ボクセルの不透明度に基づいて前記光強度を減衰させるとともに、当該ボクセルの不透明度及び色値並びに前記減衰させた光強度に基づいて累積輝度値を算出する処理を繰り返し、算出された累積輝度値に基づいて、前記ボリュームレンダリング像の当該画素に対応する画素値として与えるレイキャスティング手段と、を備えるようにしても良い。これにより、レイキャスティング法によるレンダリング処理が実行される。
【0030】
また、前記探索制御マスク作成手段及び前記レイキャスティング手段は、前記座標変換を行って前記ボクセルデータより不透明度または色値を取得する際、
所定の回転を定義した回転行列、ＸＹＺ軸各方向のオフセット値、ＸＹＺ軸方向の拡大又は縮小倍率、Ｚ軸方向の変倍率、注視点から視点までの距離を含む前記所定の座標変換のパラメータを取得し、前記視点座標系のボクセルの整数値の座標を、前記パラメータに基づいて前記ボクセルデータの座標系に変換を行って、前記ボクセルデータの実数値の座標を算出し、算出した実数値の座標の近傍の複数の整数値の座標に対応する前記ボクセルデータの複数のボクセルを特定し、特定した複数のボクセルの不透明度または色値に基づいて前記ボクセルデータより取得される不透明度または色値として算出するようにしても良い。これにより、レンダリング処理中の座標変換処理が実行される。
所定の回転を定義した回転行列、ＸＹＺ軸各方向のオフセット値、ＸＹＺ軸方向の拡大又は縮小倍率、Ｚ軸方向の変倍率、注視点から視点までの距離を含む前記所定の座標変換のパラメータを取得し、前記視点座標系のボクセルの整数値の座標を、前記パラメータに基づいて前記ボクセルデータの座標系に変換を行って、前記ボクセルデータの実数値の座標を算出し、算出した実数値の座標の近傍の複数の整数値の座標に対応する前記ボクセルデータの複数のボクセルを特定し、特定した複数のボクセルの不透明度または色値に基づいて前記ボクセルデータより取得される不透明度または色値として算出するようにしても良い。これにより、レンダリング処理中の座標変換処理が実行される。
【0031】
また、前記レンダリング手段は、前記ボクセルデータに基づいて３Ｄテクスチャを生成する３Ｄテクスチャ生成手段と、前記３Ｄテクスチャに対して所定の座標変換を行って変換後３Ｄテクスチャを生成する座標変換手段と、３次元空間のＸＹ座標面上の四角形をＺ軸方向に並べた積層四角形を設定する積層四角形設定手段と、所定の視点からＺ軸方向に平行な視線上の前記四角形のＸＹ座標に対応する前記変換後３Ｄテクスチャのボクセルの色値を前記ボクセルの不透明度に基づいて前記視点から遠い四角形の順にアルファブレンディングして取得し、前記ボリュームレンダリング像の画素値として与える画素値算出手段と、を備えるようにしても良い。これにより、３Ｄテクスチャマッピング法によるレンダリング処理が実行される。
【0032】
また、前記座標変換手段は、所定の回転を定義した回転行列、視野角度、視点位置、クリッピング位置、ＸＹＺ軸各方向のオフセット値、ＸＹＺ軸方向の拡大又は縮小倍率、Ｚ軸方向の変倍率を含む所定の座標変換のパラメータを取得し、前記３Ｄテクスチャに対して、前記取得したパタメータを用いた前記所定の座標変換を行って前記変換後３Ｄテクスチャを生成するようにしても良い。これにより、レンダリング処理中の座標変換処理が実行される。
【0033】
また、前記座標変換手段及び前記画素値算出手段は、ビデオカードに搭載されたＧＰＵ及びフレームメモリを用いて実行するようにしても良い。これにより、レンダリング処理中の座標変換処理及び画像値算出処理が、高速に実行される。
【0034】
また、本開示の一実施形態によると、コンピュータの制御部が、複数の断層画像を取得する画像取得ステップと、信号値と色値及び不透明度との対応関係を定義するカラーマップを取得するカラーマップ取得ステップと、前記カラーマップを参照することで、前記断層画像の各画素の信号値を、信号値に応じた色値及び不透明度に変換し、色値及び不透明度を保持するボクセルの集合であるボクセルデータを作成するボクセル作成ステップと、前記ボクセルデータに基づいてボリュームレンダリング像を生成するレンダリングステップと、前記断層画像毎に、断層画像の被写体領域の外郭を所定の幾何形状で近似し、当該幾何形状の情報を取得する幾何情報取得ステップと、を備え、前記ボクセル作成ステップは、前記断層画像の各画素の信号値に応じた前記不透明度を、当該断層画像に対応する幾何形状の情報及び当該断層画像のスライス位置の情報に基づいて補正する医用画像処理方法が提供される。
本開示の一実施形態によれば、カラーマップを参照して断層画像の各画素の信号値を色及び不透明度に変換してボクセルデータを作成する際、カラーマップ（オパシティカーブ）に定義された不透明度を、断層画像の被写体領域の外郭を近似した幾何形状の情報及び当該断層画像のスライス位置の情報に基づいて補正する。このように、被写体形状（幾何形状）及び当該断層画像のスライス位置の情報に基づいて不透明度をコントロールすることで、従来のように３次元マスクを作成することなく、特定の領域を鮮明に可視化させたレンダリング像を得ることができる。
本開示の一実施形態によれば、カラーマップを参照して断層画像の各画素の信号値を色及び不透明度に変換してボクセルデータを作成する際、カラーマップ（オパシティカーブ）に定義された不透明度を、断層画像の被写体領域の外郭を近似した幾何形状の情報及び当該断層画像のスライス位置の情報に基づいて補正する。このように、被写体形状（幾何形状）及び当該断層画像のスライス位置の情報に基づいて不透明度をコントロールすることで、従来のように３次元マスクを作成することなく、特定の領域を鮮明に可視化させたレンダリング像を得ることができる。
【0035】
また、本開示の一実施形態によると、コンピュータの制御部が、複数の断層画像を取得する画像取得ステップと、信号値と色値及び不透明度との対応関係を定義するカラーマップを取得するカラーマップ取得ステップと、前記カラーマップを参照することで、前記断層画像の各画素の信号値を、信号値に応じた色値及び不透明度に変換し、色値及び不透明度を保持するボクセルの集合であるボクセルデータを作成するボクセル作成ステップと、前記ボクセルデータに基づいてボリュームレンダリング像を生成するレンダリングステップと、前記断層画像毎に、断層画像の被写体領域の外郭を所定の幾何形状で近似し、当該幾何形状の情報を取得する幾何情報取得ステップと、前記断層画像毎に、断層画像の各画素に対応する不透明度の補正倍率を、当該断層画像に対応する幾何形状の情報に基づいて算出する補正倍率算出ステップと、を備え、前記補正倍率算出ステップは、前記被写体領域の中心部から離れた骨領域や体表を透過させ、中心部に近い臓器を鮮明に可視化するように前記補正倍率を、前記幾何形状の幾何中心から遠ざかるにつれ小さくし、前記ボクセル作成ステップは、前記断層画像の各画素の信号値に応じた前記不透明度を、当該断層画像に対応する幾何形状の情報に基づいて補正し、かつ前記断層画像の各画素の信号値に応じた前記不透明度に対して、前記補正倍率算出ステップにより算出された補正倍率を乗算する医用画像処理方法が提供される。
本開示の一実施形態によれば、カラーマップを参照して断層画像の各画素の信号値を色及び不透明度に変換してボクセルデータを作成する際、カラーマップ（オパシティカーブ）に定義された不透明度を、断層画像の被写体領域の外郭を近似した幾何形状の情報に基づいて補正する。このように、被写体形状（幾何形状）
に基づいて不透明度をコントロールすることで、従来のように３次元マスクを作成することなく、特定の領域を鮮明に可視化させたレンダリング像を得ることができる。また、前記断層画像毎に、断層画像の各画素に対応する不透明度の補正倍率を、当該断層画像に対応する幾何形状の情報に基づいて算出する際に、被写体領域の中心部から離れた骨領域や体表を透過させ、中心部に近い臓器を鮮明に可視化するように前記補正倍率を、前記幾何形状の幾何中心から遠ざかるにつれ小さくする。断層画像の各画素の信号値に応じた不透明度に対して、幾何形状の情報に基づいて算出された補正倍率を乗算することで、不透明度を補正することができる。また、被写体領域の中心から遠ざかるつれ不透明度が減衰するので、被写体領域の中心部から離れた骨領域や体表が透過され、被写体領域の中心部に近い臓器等が鮮明に可視化される。
本開示の一実施形態によれば、カラーマップを参照して断層画像の各画素の信号値を色及び不透明度に変換してボクセルデータを作成する際、カラーマップ（オパシティカーブ）に定義された不透明度を、断層画像の被写体領域の外郭を近似した幾何形状の情報に基づいて補正する。このように、被写体形状（幾何形状）
に基づいて不透明度をコントロールすることで、従来のように３次元マスクを作成することなく、特定の領域を鮮明に可視化させたレンダリング像を得ることができる。また、前記断層画像毎に、断層画像の各画素に対応する不透明度の補正倍率を、当該断層画像に対応する幾何形状の情報に基づいて算出する際に、被写体領域の中心部から離れた骨領域や体表を透過させ、中心部に近い臓器を鮮明に可視化するように前記補正倍率を、前記幾何形状の幾何中心から遠ざかるにつれ小さくする。断層画像の各画素の信号値に応じた不透明度に対して、幾何形状の情報に基づいて算出された補正倍率を乗算することで、不透明度を補正することができる。また、被写体領域の中心から遠ざかるつれ不透明度が減衰するので、被写体領域の中心部から離れた骨領域や体表が透過され、被写体領域の中心部に近い臓器等が鮮明に可視化される。
【0036】
また、本開示の一実施形態によると、コンピュータを、複数の断層画像を取得する画像取得手段、信号値と色値及び不透明度との対応関係を定義するカラーマップを取得するカラーマップ取得手段、前記カラーマップを参照することで、前記断層画像の各画素の信号値を、信号値に応じた色値及び不透明度に変換し、色値及び不透明度を保持するボクセルの集合であるボクセルデータを作成するボクセル作成手段、前記ボクセルデータに基づいてボリュームレンダリング像を生成するレンダリング手段、前記断層画像毎に、断層画像の被写体領域の外郭を所定の幾何形状で近似し、当該幾何形状の情報を取得する幾何情報取得手段、として機能させ、前記ボクセル作成手段は、前記断層画像の各画素の信号値に応じた前記不透明度を、当該断層画像に対応する幾何形状の情報及び当該断層画像のスライス位置の情報に基づいて補正するプログラムが提供される。
また、本開示の一実施形態によると、コンピュータを、複数の断層画像を取得する画像取得手段、信号値と色値及び不透明度との対応関係を定義するカラーマップを取得するカラーマップ取得手段、前記カラーマップを参照することで、前記断層画像の各画素の信号値を、信号値に応じた色値及び不透明度に変換し、色値及び不透明度を保持するボクセルの集合であるボクセルデータを作成するボクセル作成手段、前記ボクセルデータに基づいてボリュームレンダリング像を生成するレンダリング手段、前記断層画像毎に、断層画像の被写体領域の外郭を所定の幾何形状で近似し、当該幾何形状の情報を取得する幾何情報取得手段、前記断層画像毎に、断層画像の各画素に対応する不透明度の補正倍率を、当該断層画像に対応する幾何形状の情報に基づいて算出する補正倍率算出手段、として機能させ、前記補正倍率算出手段は、前記被写体領域の中心部から離れた骨領域や体表を透過させ、中心部に近い臓器を鮮明に可視化するように前記補正倍率を、前記幾何形状の幾何中心から遠ざかるにつれ小さくし、前記ボクセル作成手段は、前記断層画像の各画素の信号値に応じた前記不透明度を、当該断層画像に対応する幾何形状の情報に基づいて補正し、かつ前記断層画像の各画素の信号値に応じた前記不透明度に対して、前記補正倍率算出手段により算出された補正倍率を乗算するプログラムが提供される。
本開示の一実施形態に係るプログラムを汎用のコンピュータにインストールすることによって、本開示の一実施形態に係る医用画像処理装置を得ることができる。
また、本開示の一実施形態によると、コンピュータを、複数の断層画像を取得する画像取得手段、信号値と色値及び不透明度との対応関係を定義するカラーマップを取得するカラーマップ取得手段、前記カラーマップを参照することで、前記断層画像の各画素の信号値を、信号値に応じた色値及び不透明度に変換し、色値及び不透明度を保持するボクセルの集合であるボクセルデータを作成するボクセル作成手段、前記ボクセルデータに基づいてボリュームレンダリング像を生成するレンダリング手段、前記断層画像毎に、断層画像の被写体領域の外郭を所定の幾何形状で近似し、当該幾何形状の情報を取得する幾何情報取得手段、前記断層画像毎に、断層画像の各画素に対応する不透明度の補正倍率を、当該断層画像に対応する幾何形状の情報に基づいて算出する補正倍率算出手段、として機能させ、前記補正倍率算出手段は、前記被写体領域の中心部から離れた骨領域や体表を透過させ、中心部に近い臓器を鮮明に可視化するように前記補正倍率を、前記幾何形状の幾何中心から遠ざかるにつれ小さくし、前記ボクセル作成手段は、前記断層画像の各画素の信号値に応じた前記不透明度を、当該断層画像に対応する幾何形状の情報に基づいて補正し、かつ前記断層画像の各画素の信号値に応じた前記不透明度に対して、前記補正倍率算出手段により算出された補正倍率を乗算するプログラムが提供される。
本開示の一実施形態に係るプログラムを汎用のコンピュータにインストールすることによって、本開示の一実施形態に係る医用画像処理装置を得ることができる。
【0037】
また、本開示の一実施形態によると、コンピュータの制御部が、複数の断層画像を取得する画像取得ステップと、前記断層画像毎に、断層画像の被写体領域の外郭を所定の幾何形状で近似し、当該幾何形状の情報を取得する幾何情報取得ステップと、前記断層画像の各画素の信号値に応じたオパシティカーブにより規定される不透明度を、当該断層画像に対応する幾何形状の情報及び当該断層画像のスライス位置の情報に基づいて補正する不透明度補正ステップと、前記断層画像の各画素に対応する補正後の不透明度を格納したデータを作成するデータ作成ステップと、を含むデータ作成方法が提供される。本開示の一実施形態によれば、断層画像の各画素に適用される不透明度を参照用のデータとして作成しておくことができるので、断層画像からボクセルデータを作成する際、不透明度を都度算出する必要がなくなり、レンダリング結果を得るまでの処理時間を短縮できる。
【発明の効果】
【0038】
本開示により、特定の領域を鮮明に可視化させることできる。
【図面の簡単な説明】
【0039】
【図1】医用画像処理装置１のハードウェア構成を示す図
【図2】ボリュームレンダリング処理の概要を示す図
【図3】医用画像処理装置１の機能構成を示す図
【図4】補正テーブル作成部２６の内部構成を示す図
【図5】レンダリング部２８の内部構成を示す図
【図6】医用画像処理装置１の全体動作を示すフローチャート
【図7】調整されたカラーマップＣｍａｐの性質を示す図
【図8】補正テーブル作成処理を示すフローチャート
【図9】楕円パラメータＰ（ｚ）を取得する処理を説明する図
【図10】楕円パラメータＰ（ｚ）の補正について説明する図
【図11】断層画像Ｄｏ（ｚ）毎に、被写体形状の外郭が楕円形状で近似されていく様子を示す図
【図12】補正倍率Ｓαの分布パターンを示す図
【図13】ＸＹ方向の振幅倍率と補正倍率との関係を示す図
【図14】ボクセル作成処理を示すフローチャート
【図15】レイキャスティング法の概要を示す図
【図16】レイキャスティング法によるレンダリング処理を示すフローチャート
【図17】探索制御マスク算出処理を示すフローチャート
【図18】探索制御マスク算出処理を説明する図
【図19】不透明ボクセル探索処理を示すフローチャート
【図20】レイキャスティング処理を示すフローチャート
【図21】３Ｄテクスチャマッピング法ついて説明する図
【図22】色補正処理を示すフローチャート
【図23】３Ｄテクスチャマッピング法によるレンダリング処理を示すフローチャート
【図24】投影画面設定処理を示すフローチャート
【図25】レンダリング処理を示すフローチャート
【図26】従来手法により生成したレンダリング像（正面）
【図27】提案手法により生成したレンダリング像（正面）
【図28】提案手法により生成したレンダリング像（背面）
【図29】提案手法により生成したレンダリング像（側面）
【図30】第２の実施の形態における補正倍率Ｓα（ｘ、ｙ、ｚ）の分布パターンを示す図
【図31】ＸＹＺ方向の振幅倍率と補正倍率との関係を示す図
【図32】第３の実施形態における医用画像処理装置１ａの機能構成を示す図
【図33】第３の実施形態における医用画像処理装置１ａの全体動作を示すフローチャート
【図34】パラメータ調整画面４０の表示例を示す図
【図35】第４の実施形態における医用画像処理装置１ｂの機能構成を示す図
【図36】体軸断面、冠状断面、矢状断面について説明する図
【図37】第４の実施形態における医用画像処理装置１ｂの全体動作を示すフローチャート
【図38】ＭＰＲ像の表示例を示す図
【図39】ＭＰＲ像の表示例を示す図
【発明を実施するための形態】
【0040】
以下図面に基づいて、本開示の実施の形態を詳細に説明する。本実施の形態では、Ｘ線ＣＴ装置により撮影されたＣＴ画像に基づいてボリュームレンダリング像を生成する場合について説明するが、本開示はＣＴ画像以外の医用画像（ＭＲＩ画像やＰＥＴ画像等）に基づいてボリュームレンダリング像を生成する場合にも適用可能である。
【0041】
［第１の実施の形態］
【0042】
図１は、本実施の形態における医用画像処理装置１のハードウェア構成を示すブロック図である。図１に示すように、医用画像処理装置１は、制御部１１、記憶部１２、メディア入出力部１３、通信制御部１４、入力部１５、表示部１６、周辺機器Ｉ／Ｆ部１７等が、バス１８を介して接続される汎用のコンピュータで実現される。但し、これに限ることなく、用途、目的に応じて様々な構成を採ることが可能である。
【0043】
制御部１１は、ＣＰＵ（Central Processing Unit）、ＧＰＵ（Graphics Processing Unit）、ＲＯＭ（Read Only
Memory）、ＲＡＭ（Random Access Memory）、フレームメモリ(Frame Memory)等によって構成される。ＣＰＵ、ＧＰＵは、記憶部１２、ＲＯＭ、記録媒体等に格納されるプログラムをＲＡＭ、フレームメモリ上のワークメモリ領域に呼び出して実行し、バス１８を介して接続された各装置を駆動制御し、医用画像処理装置１が行う後述する処理を実現する。
Memory）、ＲＡＭ（Random Access Memory）、フレームメモリ(Frame Memory)等によって構成される。ＣＰＵ、ＧＰＵは、記憶部１２、ＲＯＭ、記録媒体等に格納されるプログラムをＲＡＭ、フレームメモリ上のワークメモリ領域に呼び出して実行し、バス１８を介して接続された各装置を駆動制御し、医用画像処理装置１が行う後述する処理を実現する。
【0044】
ＧＰＵ及びフレームメモリは、医用画像処理装置１のビデオカードに搭載されている。ＲＯＭは、不揮発性メモリであり、コンピュータのブートプログラムやＢＩＯＳ等のプログラム、データ等を恒久的に保持している。ＲＡＭ、フレームメモリは、揮発性メモリであり、記憶部１２、ＲＯＭ、記録媒体等からロードしたプログラム、データ等を一時的に保持するとともに、制御部１１が各種処理を行う為に使用するワークエリアを備える。
【0045】
記憶部１２は、ＨＤＤ（Hard Disk Drive）等であり、制御部１１が実行するプログラム、プログラム実行に必要なデータ、ＯＳ（Operating System）等が格納される。プログラムに関しては、ＯＳに相当する制御プログラムや、後述する処理をコンピュータに実行させるためのアプリケーションプログラムが格納されている。これらの各プログラムコードは、制御部１１により必要に応じて読み出されてＲＡＭ、フレームメモリに移され、ＣＰＵ、ＧＰＵに読み出されて各種の手段として実行される。
【0046】
メディア入出力部１３（ドライブ装置）は、データの入出力を行い、例えば、ＣＤドライブ（－ＲＯＭ、－Ｒ、－ＲＷ等）、ＤＶＤドライブ（－ＲＯＭ、－Ｒ、－ＲＷ等）等のメディア入出力装置を有する。通信制御部１４は、通信制御装置、通信ポート等を有し、コンピュータとネットワーク間の通信を媒介する通信インタフェースであり、ネットワークを介して、他のコンピュータ間との通信制御を行う。ネットワークは、有線、無線を問わない。
【0047】
入力部１５は、データの入力を行い、例えば、キーボード、マウス等のポインティングデバイス、テンキー等の入力装置を有する。入力部１５を介して、コンピュータに対して、操作指示、動作指示、データ入力等を行うことができる。表示部１６は、液晶パネル等のディスプレイ装置、ディスプレイ装置と連携してコンピュータのビデオ機能を実現するための論理回路等（ビデオアダプタ等）を有する。なお、入力部１５及び表示部１６は、タッチパネルディスプレイのように、一体となっていてもよい。
【0048】
周辺機器Ｉ／Ｆ（Interface）部１７は、コンピュータに周辺機器を接続させるためのポートであり、周辺機器Ｉ／Ｆ部１７を介してコンピュータは周辺機器とのデータの送受信を行う。周辺機器Ｉ／Ｆ部１７は、ＵＳＢ（Universal Serial Bus）やＩＥＥＥ１３９４やＲＳ－２３２Ｃ等によって構成されており、通常複数の周辺機器Ｉ／Ｆを有する。周辺機器との接続形態は有線、無線を問わない。バス１８は、各装置間の制御信号、データ信号等の授受を媒介する経路である。
【0049】
医用画像処理装置１は、１台のコンピュータで構成されてもよいし、複数のコンピュータがネットワークを介して構成されてもよい。例えば、医用画像処理装置１が、サーバとクライアント端末で構成される場合、クライアント端末においてデータの入力を受け付けて、サーバが各種の処理を行い、クライアント端末が処理結果を表示するようにしてもよい。以下の説明では、簡素な構成例として、医用画像処理装置１が１台のコンピュータで構成された例を説明する。
【0050】
図２は、医用画像処理装置１が実行するボリュームレンダリング処理の概要を示す図である。図２に示すように、複数の断層画像（図の例では、５１２×５１２ピクセルの３７０枚の胸部ＣＴ画像）に基づいてボリュームレンダリング処理を実行し、ボリュームレンダリング像として、遠位の視点から被写体を観察する全体レンダリング像（図左下）や視点を被写体内に自由に移動させて気管支や大腸等を観察する仮想内視鏡像（図右下）を生成し表示する。本実施の形態では、主に、全体レンダリング像（図左下）を生成する場合について説明する。
【0051】
図３は、医用画像処理装置１の機能構成を示す図である。図３に示すように、医用画像処理装置１は、画像取得部２１、階調圧縮部２２、領域指定部２３、カラーマップ取得部２４、カラーマップ調整部２５、補正テーブル作成部２６、ボクセル作成部２７、レンダリング部２８、及びデータ出力部２９を備える。
【0052】
画像取得部２１は、ＣＴ装置により撮影された断層画像群Ｄｏを取得する。断層画像群Ｄｏは、被写体（人体）を頭尾軸に沿って所定のスライス間隔で連続的に撮影した複数の断層画像（ＣＴ画像）からなる。各断層画像はＤＩＣＯＭ形式の２次元の画像データである。ＤＩＣＯＭ形式は、１ファイルにヘッダ部と画像データ部を含む医療画像で一般的に用いられる画像フォーマットであり、画像撮影時のパラメータや診断情報を保存しておくことができる。
【0053】
１つの断層画像は、例えば、５１２×５１２ピクセルの画像（ＣＴ画像）である。断層画像の各画素には、信号値ｖが付与されており、ＣＴ画像の場合、信号値ｖはＣＴ値である。本実施の形態では、信号値ｖ（ＣＴ値）は１６ビット（－３２７６８≦ｖ≦３２７６７）のデータとする（但し、信号値ｖのビット数は特に限定されない）。
【0054】
ＣＴ値は、水を基準として表現した組織のＸ線減弱係数であり、ＣＴ値により組織や病変の種類等が判断できるようになっている(単位はＨＵ（Hounsfield Unit））。ＣＴ値は、水と空気のＸ線減弱係数で標準化されており、水のＣＴ値は０、空気のＣＴ値を－１０００である。また、脂肪のＣＴ値は－１２０～－１００程度であり、通常組織のＣＴ値は０～１２０程度であり、骨のＣＴ値は１０００前後を示す。
【0055】
断層画像群Ｄｏは、ＸＹの２次元データである断層画像をＺ軸方向に積層したものであり、ＸＹＺの３次元データとして表現可能である。例えば断層画像群Ｄｏは、以下のように、ＸＹＺの３次元データ（Ｄｏ（ｘ、ｙ、ｚ））として定義される。尚、本開示において、Ｘ軸は人体の左右軸、Ｙ軸は人体の背腹軸（前後軸）、Ｚ軸は人体の頭尾軸（上下軸）に相当するものとする。ＸＹ軸は、例えば、ＣＴ装置での撮影時に設定される。また、本開示において、Ｘ軸方向を横方向、Ｙ軸方向を縦方向と呼ぶ場合がある。
【0056】
（式１）
－３２７６８≦Ｄｏ（ｘ、ｙ、ｚ）≦３２７６７
０≦ｘ≦Ｓｘ－１、０≦ｙ≦Ｓｙ－１、０≦ｚ≦Ｓｚ－１
解像度：Ｒｘｙ、解像度Ｒｚ
－３２７６８≦Ｄｏ（ｘ、ｙ、ｚ）≦３２７６７
０≦ｘ≦Ｓｘ－１、０≦ｙ≦Ｓｙ－１、０≦ｚ≦Ｓｚ－１
解像度：Ｒｘｙ、解像度Ｒｚ
【0057】
（式１）において、Ｓｘは断層画像の横方向（Ｘ軸方向）の画素数、Ｓｙは断層画像の縦方向（Ｙ軸方向）の画素数、Ｓｚはスライス枚数（断層画像の枚数）を表す。ＲｘｙはＸ軸方向及びＹ軸方向の断層画像の解像度であり、画素の間隔の逆数、すなわち単位距離あたりの画素数を示す。Ｒｚはスライスの解像度であり、断層画像のスライス間隔（例えば、０．５ｍｍや１ｍｍ）の逆数、すなわち単位距離あたりのスライス枚数を表す。
【0058】
尚、断層画像群Ｄｏのうち、所定のｚ番目（０≦ｚ≦Ｓｚ－１）の断層画像を表す場合は、「断層画像Ｄｏ（ｚ）」と表記することがある。
【0059】
階調圧縮部２２は、画像取得部２１により取得した断層画像群Ｄｏの信号値を、必要に応じて８ビットに階調圧縮する（後述の（式６）参照）。８ビットに階調圧縮された断層画像Ｄｏ（ｘ、ｙ、ｚ）は以下のように定義される。
【0060】
（式２）
０≦Ｄｏ（ｘ、ｙ、ｚ）≦２５５
０≦ｘ≦Ｓｘ－１、０≦ｙ≦Ｓｙ－１、０≦ｚ≦Ｓｚ－１
解像度：Ｒｘｙ、解像度Ｒｚ
０≦Ｄｏ（ｘ、ｙ、ｚ）≦２５５
０≦ｘ≦Ｓｘ－１、０≦ｙ≦Ｓｙ－１、０≦ｚ≦Ｓｚ－１
解像度：Ｒｘｙ、解像度Ｒｚ
【0061】
このように階調圧縮をすることで、断層画像群を保持するためのメモリ容量を半分に抑えることができる。たとえ信号値の階調が１６ビットあっても、カラーマップＣｍａｐにより、変換される色値（具体的にはＲＧＢ値）及び不透明度の階調はディスプレイの階調に合わせて８ビットに制限されるため、階調圧縮に伴う画質劣化は殆ど生じない。
【0062】
領域指定部２３は、ユーザから被写体の関心領域（ＲＯＩ：Region of
Interest）の指定を受け付ける。例えば、関心領域ＲＯＩを直方体で指定する場合は、以下のように、Ｘ方向ＲＯＩ（Ｘ方向の開区間）、Ｙ方向ＲＯＩ（Ｙ方向の開区間）、Ｚ方向ＲＯＩ（Ｚ方向の開区間）を指定する。
Interest）の指定を受け付ける。例えば、関心領域ＲＯＩを直方体で指定する場合は、以下のように、Ｘ方向ＲＯＩ（Ｘ方向の開区間）、Ｙ方向ＲＯＩ（Ｙ方向の開区間）、Ｚ方向ＲＯＩ（Ｚ方向の開区間）を指定する。
【0063】
（式３）
Ｘ方向ＲＯＩ＝（Ｘｓ、Ｘｅ）
Ｙ方向ＲＯＩ＝（Ｙｓ、Ｙｅ）
Ｚ方向ＲＯＩ＝（Ｚｓ、Ｚｅ）
尚、関心領域を指定しない場合は、Ｘｓ＝０、Ｘｅ＝Ｓｘ－１、Ｙｓ＝０、Ｙｅ＝Ｓｙ－１、Ｚｓ＝０、Ｚｅ＝Ｓｚ－１となる。
Ｘ方向ＲＯＩ＝（Ｘｓ、Ｘｅ）
Ｙ方向ＲＯＩ＝（Ｙｓ、Ｙｅ）
Ｚ方向ＲＯＩ＝（Ｚｓ、Ｚｅ）
尚、関心領域を指定しない場合は、Ｘｓ＝０、Ｘｅ＝Ｓｘ－１、Ｙｓ＝０、Ｙｅ＝Ｓｙ－１、Ｚｓ＝０、Ｚｅ＝Ｓｚ－１となる。
【0064】
カラーマップ取得部２４は、断層画像群Ｄｏに適用するカラーマップＣｍａｐを取得する。カラーマップ取得部２４は、取得したカラーマップＣｍａｐを更に調整する場合には、カラーマップ調整部２５に出力し、カラーマップＣｍａｐの調整を行わない場合には、ボクセル作成部２７に出力する。
【0065】
カラーマップＣｍａｐは、信号値ｖと色値（具体的にはＲＧＢ値）及び不透明度(α値)との対応関係を定義するものであり、信号値ｖを２４ビットの色（ＲＧＢ値）及び８ビットの不透明度（α値）に変換する関数（実体的には２次元のデータテーブル）として表現可能である。例えば、１６ビットの断層画像群Ｄｏ（（式１）参照）に適用されるカラーマップＣｍａｐは次のように定義される。
【0066】
（式４）
０≦Ｃｍａｐ（ｖ、ｎ）≦２５５
－３２７６８≦ｖ≦３２７６７
ｎ＝０（Ｒ）、１（Ｇ）、２（Ｂ）、３（α）
０≦Ｃｍａｐ（ｖ、ｎ）≦２５５
－３２７６８≦ｖ≦３２７６７
ｎ＝０（Ｒ）、１（Ｇ）、２（Ｂ）、３（α）
【0067】
また、８ビットの断層画像群Ｄｏ（（式２）参照）に適用されるカラーマップＣｍａｐ（ｖ、ｎ）は次のように定義される。
【0068】
（式５）
０≦Ｃｍａｐ（ｖ、ｎ）≦２５５
０≦ｖ≦２５５
ｎ＝０（Ｒ）、１（Ｇ）、２（Ｂ）、３（α）
０≦Ｃｍａｐ（ｖ、ｎ）≦２５５
０≦ｖ≦２５５
ｎ＝０（Ｒ）、１（Ｇ）、２（Ｂ）、３（α）
【0069】
（式４）、（式５）に示すカラーマップＣｍａｐ（ｖ、ｎ）（０≦ｎ≦３）のうち、特にカラーマップＣｍａｐ（ｖ、ｎ）（０≦ｎ≦２）は、信号値ｖを色値（ＲＧＢ値）に変換する関数に相当する。本開示では、信号値ｖを色値に変換するカラーマップＣｍａｐ（ｖ、ｎ）（０≦ｎ≦２）を「カラーパレット」と称する。
【0070】
また、（式４）、（式５）に示すカラーマップＣｍａｐ（ｖ、ｎ）（０≦ｎ≦３）のうち、特にカラーマップＣｍａｐ（ｖ、３）は、信号値ｖを不透明度に変換する関数に相当し、一般的に「オパシティカーブ」と呼ばれる。
【0071】
カラーマップＣｍａｐ（ｖ、３）（オパシティカーブ）は、例えば、通常組織（１６ビットの信号値ｖ＝０～１２０程度）の不透明度が２５５に設定（不透明度＝２５５は不透明（光がすべて反射）であることを示す）され、骨領域（１６ビットの信号値ｖ＝１０００前後）の不透明度が１～２５４の中間値に設定（不透明度＝１～２５４は半透明（入射光の一部が反射され、その他は透過）であることを示す）され、その他の組織等の部位を不透明度が０に設定（不透明度＝０は透明（入射光の全てが透過）であることを示す）される。
【0072】
カラーマップ調整部２５は、必要に応じて、カラーマップＣｍａｐ（ｖ、ｎ）のうち、信号値ｖと不透明度との対応関係を定義するＣｍａｐ（ｖ、３）（オパシティカーブともいう）を信号値に応じて調整し（後述の（式７）（式８）参照）、調整したカラーマップＣｍａｐをボクセル作成部２７に出力する。
【0073】
補正テーブル作成部２６は、断層画像群Ｄｏの各画素に対応する不透明度の補正倍率を格納した補正テーブルＳαを作成し、ボクセル作成部２７に出力する。
【0074】
補正テーブル作成部２６は、図３に示すように、幾何情報取得部２６ａ、幾何情報補正部２６ｂ、及び補正倍率算出部２６ｃから構成される。
幾何情報取得部２６ａは、断層画像Ｄｏ（ｚ）毎に、断層画像Ｄｏ（ｚ）の被写体領域の外郭を所定の幾何形状で近似し、幾何形状の情報（幾何パラメータＰ（ｚ）、後述の（式９）参照）を取得する。
幾何情報取得部２６ａは、断層画像Ｄｏ（ｚ）毎に、断層画像Ｄｏ（ｚ）の被写体領域の外郭を所定の幾何形状で近似し、幾何形状の情報（幾何パラメータＰ（ｚ）、後述の（式９）参照）を取得する。
【0075】
幾何情報補正部２６ｂは、幾何形状の近似精度を所定の基準で判断し、近似精度が良好でない断層画像Ｄｏ（ｚ）の幾何パラメータＰ（ｚ）を、近似精度が良好な他の断層画像Ｄｏ（ｚ）の幾何パラメータＰ（ｚ）に基づいて補正する（後述の（式１０）参照）。
【0076】
補正倍率算出部２６ｃは、断層画像Ｄｏ（ｚ）毎に、断層画像Ｄｏ（ｚ）に対応する幾何パラメータＰ（ｚ）に基づいて、断層画像Ｄｏ（ｚ）の各画素（ｘ、ｙ、ｚ）に対応する不透明度の補正倍率Ｓα（ｘ、ｙ、ｚ）を算出し、補正倍率Ｓα（ｘ、ｙ、ｚ）を格納した３次元の補正テーブルＳαを作成する（後述の（式１１）参照）。
【0077】
この際、補正倍率Ｓα（ｘ、ｙ、ｚ）は、断層画像Ｄｏ（ｚ）の被写体領域の外郭を近似した幾何形状の幾何中心（後述する実施形態では、楕円の中心）から遠ざかるにつれその倍率が小さくなるように算出される。
【0078】
尚、補正倍率Ｓα（ｘ、ｙ、ｚ）は３次元のデータ配列であり、そのままデータテーブル（補正テーブル）として扱えるため、補正倍率と補正テーブルは実質的に同一である。このため、補正倍率と補正テーブルは同一の符号「Ｓα」で表す。
【0079】
ボクセル作成部２７は、断層画像群Ｄｏ、カラーマップ取得部２４により取得されたカラーマップＣｍａｐ（又はカラーマップ調整部２５により調整されたカラーマップＣｍａｐ）、補正テーブル作成部２６により作成された補正テーブルＳα、及び領域指定部２３により指定された関心領域ＲＯＩに基づいて、ボクセル構造体Ｖ（ボクセルデータ）を作成し、レンダリング部２８に出力する。
【0080】
具体的には、ボクセル作成部２７は、カラーマップＣｍａｐを参照することで、断層画像群Ｄｏの各画素（ｘ、ｙ、ｚ）の信号値ｖを、信号値ｖに応じた色値及び不透明度に変換し、色値及び不透明度を保持するボクセルの集合であるボクセル構造体Ｖ（ボクセルデータ）を作成する。
【0081】
特に本開示では、ボクセル作成部２７は、カラーマップＣｍａｐを参照して断層画像群Ｄｏの各画素（ｘ、ｙ、ｚ）の信号値ｖを色値及び不透明度に変換してボクセル構造体Ｖ（ボクセルデータ）を作成する際、カラーマップＣｍａｐ（ｖ、３）（オパシティカーブ）に定義された不透明度をそのまま用いるのではなく、各画素（ｘ、ｙ、ｚ）に対応する補正倍率Ｓα（ｘ、ｙ、ｚ）を乗算することで不透明度を補正する（後述の（式１３）参照）。
【0082】
前述したように、補正倍率Ｓα（ｘ、ｙ、z）は、被写体領域の外郭を近似した幾何形状の幾何中心（重心）から遠ざかるにつれ倍率が小さくなるので、被写体領域の中心から離れたボクセルほど不透明度が減衰する。これにより、骨領域（胸骨や頭蓋骨）や皮膚、外部の寝台、固定用治具等が透明化され、被写体領域の中心付近にある臓器等が鮮明に可視化されたレンダリング像を得ることができる。
【0083】
レンダリング部２８は、ボクセル作成部２７により作成されたボクセル構造体Ｖ（ボクセルデータ）に基づいて、ボリュームレンダリング処理を実行し、ボリュームレンダリング像（３次元画像）を生成し表示する。レンダリング部２８は、図５に示すように、第１レンダリング部２８ａ及び第２レンダリング部２８ｂを含む。
【0084】
第１レンダリング部２８ａは、ＣＰＵによりレンダリング処理を実行するレンダリング部であり、第２レンダリング部２８ｂは、コンピュータグラフィックスＡＰＩのオープン標準規格であるＯｐｅｎＧＬ（Open Graphics Library）を用いてＧＰＵによりレンダリング処理を実行するレンダリング部である。レンダリング処理は、第１レンダリング部２８ａ又は第２レンダリング部２８ｂにより実行される。
【0085】
データ出力部２９は、ユーザからデータの保存操作を受け付けることによって、各種データを出力して記憶部１２に保存する。例えば、データ出力部２９は、レンダリング部２８により生成し表示されたボリュームレンダリング像（３次元画像）を出力して保存する。
【0086】
また、データ出力部２９は、ボクセル構造体Ｖ（ボクセルデータ）を作成する際に参照可能なデータ（例えば、補正後の不透明度を格納した不透明度テーブルＴα等）を出力しても良い。これにより、同じ断層画像群Ｄｏに基づいて再度レンダリングを行う際、いくつかの処理（例えば、補正テーブルＳαの作成処理等）を省略することができ、レンダリング像を得るまでの処理時間を短縮できる。
【0087】
【0088】
まず、医用画像処理装置１の制御部１１（画像取得部２１）は、レンダリング対象となるＤＩＣＯＭ形式の複数の断層画像（断層画像群）Ｄｏを取得する（図６のステップＳ１）。断層画像群Ｄｏの取得方法は任意であるが、例えば、断層画像群Ｄｏが予め記憶部１２に記憶されている場合は、制御部１１は、記憶部１２から断層画像群Ｄｏを読込み取得することができる。
【0089】
また、断層画像群ＤｏがＵＳＢ等の記憶媒体に記憶されている場合は、制御部１１は、周辺機器Ｉ／Ｆ部１７に接続した当該記憶媒体から断層画像群Ｄｏを読込み取得することができる。或いは、サーバに断層画像群Ｄｏが記憶されている場合は、制御部１１は、通信制御部１４を介して当該サーバから断層画像群Ｄｏをダウンロードして取得することができる。
【0090】
続いて、制御部１１（階調圧縮部２２）は、必要に応じて、断層画像群Ｄｏに対して階調圧縮処理を施す（図６のステップＳ２）。具体的には、制御部１１は、以下の手順で１６ビットの断層画像群Ｄｏ（（式１）参照）を８ビットの断層画像群Ｄｏ（（式２）参照）に変換する。
【0091】
（式６）
（１）断層画像群ＤｏのＳｚ／２番目の中間スライスＤｏ（ｚ／２）における全ての画素の最小値Ｄｍｉｎ、最大値Ｄｍａｘを算出する。
（２）下限値Ｌｍｉｎ＝（Ｄｍａｘ－Ｄｍｉｎ）・γ＋Ｄｍｉｎ、上限値Ｌｍａｘ＝（Ｄｍａｘ－Ｄｍｉｎ）・（１－γ）＋Ｄｍｉｎを設定する。ここで、γは階調圧縮画像のコントラスト調整幅で、０に近いほどコントラストが増大する（但し、輝度が小さくなる）。通常はγ＝０．１に設定する。
（３）断層画像群Ｄｏを以下の計算式で８ビットの断層画像群Ｄｏに変換する（下限値Ｌｍｉｎ～上限値Ｌｍａｘの範囲で２５６段階に圧縮する）。
Ｄｏ（ｘ、ｙ、ｚ）
＝（Ｄｏ（ｘ、ｙ、ｚ）－Ｌｍｉｎ）・２５５／（Ｌｍａｘ－Ｌｍｉｎ）
但し、Ｄｏ（ｘ、ｙ、ｚ）＞２５５の場合はＤｏ（ｘ、ｙ、ｚ）＝２５５、Ｄｏ（ｘ、ｙ、ｚ）＜０の場合はＤｏ（ｘ、ｙ、ｚ）＝０に飽和させる。
（１）断層画像群ＤｏのＳｚ／２番目の中間スライスＤｏ（ｚ／２）における全ての画素の最小値Ｄｍｉｎ、最大値Ｄｍａｘを算出する。
（２）下限値Ｌｍｉｎ＝（Ｄｍａｘ－Ｄｍｉｎ）・γ＋Ｄｍｉｎ、上限値Ｌｍａｘ＝（Ｄｍａｘ－Ｄｍｉｎ）・（１－γ）＋Ｄｍｉｎを設定する。ここで、γは階調圧縮画像のコントラスト調整幅で、０に近いほどコントラストが増大する（但し、輝度が小さくなる）。通常はγ＝０．１に設定する。
（３）断層画像群Ｄｏを以下の計算式で８ビットの断層画像群Ｄｏに変換する（下限値Ｌｍｉｎ～上限値Ｌｍａｘの範囲で２５６段階に圧縮する）。
Ｄｏ（ｘ、ｙ、ｚ）
＝（Ｄｏ（ｘ、ｙ、ｚ）－Ｌｍｉｎ）・２５５／（Ｌｍａｘ－Ｌｍｉｎ）
但し、Ｄｏ（ｘ、ｙ、ｚ）＞２５５の場合はＤｏ（ｘ、ｙ、ｚ）＝２５５、Ｄｏ（ｘ、ｙ、ｚ）＜０の場合はＤｏ（ｘ、ｙ、ｚ）＝０に飽和させる。
【0092】
続いて、制御部１１（領域指定部２３）は、必要に応じて、ユーザから被写体の関心領域ＲＯＩの指定を受け付ける（図６のステップＳ３）。例えば、前記した（式３）のように、関心領域ＲＯＩを直方体等で指定する。
【0093】
続いて、制御部１１（カラーマップ取得部２４）は、断層画像群Ｄｏに適用するカラーマップＣｍａｐを取得する（ステップＳ４）。ここで、断層画像群Ｄｏの信号値が１６ビットの場合（（式１）参照）は、１６ビット対応のカラーマップＣｍａｐ（（式４）参照）を取得し、断層画像群Ｄｏの信号値が８ビットに圧縮されている場合（（式２）参照）は、８ビット対応のカラーマップＣｍａｐ（（式５）参照）を取得する。
【0094】
カラーマップＣｍａｐの取得方法は特に限定されず、例えば、制御部１１は、予め用意されている複数のカラーマップＣｍａｐの中から、断層画像群Ｄｏに適用するカラーマップＣｍａｐをユーザに選択させることで取得することができる。また、制御部１１は、前回レンダリング時に使用したカラーマップＣｍａｐを記憶部１２から読込み、断層画像群Ｄｏに適用するカラーマップＣｍａｐとして取得しても良い。
【0095】
また、制御部１１は、カラーマップ作成画面（不図示）においてユーザが作成したカラーマップＣｍａｐを取得しても良い。この場合、制御部１１は、カラーマップ作成画面において、代表的な信号値ｖ（例えば、１０箇所程度の特徴的な信号値）と当該信号値ｖに対する色値（ＲＧＢ値）及び不透明度（α値）をユーザに設定させ、ユーザが設定した代表的な信号値ｖと色値及び不透明度との対応関係に基づいて所定範囲の信号値ｖを色値及び不透明度に変換するカラーマップＣｍａｐを自動生成する。
【0096】
続いて、制御部１１は、ステップＳ４において取得したカラーマップＣｍａｐを更に調整する場合（ステップＳ５；「Ｙｅｓ」）、ステップＳ６へ移行し、カラーマップＣｍａｐの調整を行う。具体的には、制御部１１（カラーマップ調整部２５）は、以下のように、カラーマップＣｍａｐ（ｖ、ｎ）のうち、信号値ｖと不透明度との対応関係を定義するカラーマップＣｍａｐ（ｖ、３）（オパシティカーブ）を調整する。
【0097】
断層画像群Ｄｏの信号値が１６ビットの場合（（式１）参照）は、制御部１１は、次の手順で、１６ビット対応のカラーマップＣｍａｐ（ｖ、３）（（式４）参照）を調整する。
【0098】
（式７）
（１）断層画像群Ｄｏの中で、Ｓｚ／２番目の中間スライスＤｏ（ｚ／２）における全ての画素の最小値Ｄｍｉｎ、最大値Ｄｍａｘを算出する。
（２）下限値Ｌｍｉｎ＝（Ｄｍａｘ－Ｄｍｉｎ）・γ＋Ｄｍｉｎ、上限値Ｌｍａｘ＝（Ｄｍａｘ－Ｄｍｉｎ）・（１－γ）＋Ｄｍｉｎを設定する。ここで、γは階調圧縮画像のコントラスト調整幅で、０に近いほどコントラストが増大する（但し、輝度が小さくなる）。通常はγ＝０．１に設定する。
（３）１６ビットの信号値ｖを以下の計算式で８ビットの信号値ｖ８に変換する（下限値Ｌｍｉｎ～上限値Ｌｍａｘの範囲で２５６段階に圧縮する）。
ｖ８＝（ｖ－Ｌｍｉｎ）・２５５／（Ｌｍａｘ－Ｌｍｉｎ）
但し、ｖ８＞２５５の場合はｖ８＝２５５、ｖ８＜０の場合はｖ８＝０に飽和させる。
（４）次のように、１６ビット対応のカラーマップＣｍａｐ（ｖ、３）（（式４）参照）を調整する。
ｖ８＞Ｓｖの場合：Ｃｍａｐ（ｖ、３）←Ｃｍａｐ（ｖ、３）（ｖ８－Ｓｖ）δｖ８≦Ｓｖの場合：調整しない
（１）断層画像群Ｄｏの中で、Ｓｚ／２番目の中間スライスＤｏ（ｚ／２）における全ての画素の最小値Ｄｍｉｎ、最大値Ｄｍａｘを算出する。
（２）下限値Ｌｍｉｎ＝（Ｄｍａｘ－Ｄｍｉｎ）・γ＋Ｄｍｉｎ、上限値Ｌｍａｘ＝（Ｄｍａｘ－Ｄｍｉｎ）・（１－γ）＋Ｄｍｉｎを設定する。ここで、γは階調圧縮画像のコントラスト調整幅で、０に近いほどコントラストが増大する（但し、輝度が小さくなる）。通常はγ＝０．１に設定する。
（３）１６ビットの信号値ｖを以下の計算式で８ビットの信号値ｖ８に変換する（下限値Ｌｍｉｎ～上限値Ｌｍａｘの範囲で２５６段階に圧縮する）。
ｖ８＝（ｖ－Ｌｍｉｎ）・２５５／（Ｌｍａｘ－Ｌｍｉｎ）
但し、ｖ８＞２５５の場合はｖ８＝２５５、ｖ８＜０の場合はｖ８＝０に飽和させる。
（４）次のように、１６ビット対応のカラーマップＣｍａｐ（ｖ、３）（（式４）参照）を調整する。
ｖ８＞Ｓｖの場合：Ｃｍａｐ（ｖ、３）←Ｃｍａｐ（ｖ、３）（ｖ８－Ｓｖ）δｖ８≦Ｓｖの場合：調整しない
【0099】
尚、Ｃｍａｐ（ｖ、ｎ）（０≦Ｃｍａｐ（ｖ、ｎ）≦２５５、－３２７６８≦ｖ≦３２７６７、ｎ＝０（Ｒ）、１（Ｇ）、２（Ｂ）、３（α））であり、８ビット換算の閾値をＳｖとし、減衰係数をδとする（δは負の実数で、通常はδ＝－０．５）。閾値Ｓｖは、不透明度を減衰させる信号値の閾値であり、骨領域を減衰させる場合には、Ｓｖ＝１２８程度に設定する。
【0100】
一方、断層画像群Ｄｏの信号値が８ビットの場合（（式２）参照）は、制御部１１は、次のように、８ビット対応のカラーマップＣｍａｐ（（式５）参照）を調整する。
【0101】
（式８）
ｖ＞Ｓｖの場合：Ｃｍａｐ（ｖ、３）←Ｃｍａｐ（ｖ、３）（ｖ－Ｓｖ）δ
ｖ≦Ｓｖの場合：調整しない
ｖ＞Ｓｖの場合：Ｃｍａｐ（ｖ、３）←Ｃｍａｐ（ｖ、３）（ｖ－Ｓｖ）δ
ｖ≦Ｓｖの場合：調整しない
【0102】
尚、Ｃｍａｐ（ｖ、ｎ）（０≦Ｃｍａｐ（ｖ、ｎ）≦２５５、０≦ｖ≦２５５、ｎ＝０（Ｒ）、１（Ｇ）、２（Ｂ）、３（α））であり、閾値をＳｖとし、減衰係数をδとする（δは負の実数で、通常はδ＝－０．５）。閾値Ｓｖは、不透明度を減衰させる信号値の閾値であり、骨領域を減衰させる場合には、Ｓｖ＝１２８程度に設定する。
【0103】
図７（ａ）は、（式７）においてＣｍａｐ（ｖ、３）（オパシティカーブ）に乗算される（ｖ８－Ｓｖ）δの減衰特性を示す図である。図に示すように、信号値ｖ８＞閾値Ｓｖの場合、信号値ｖ８と閾値Ｓｖとの差分値（ｖ８－Ｓｖ）に応じた減衰比率（ｖ８－Ｓｖ）δで不透明度が減衰する。例えば、図７（ｂ）に示すように、内臓領域と骨領域の中心からの距離に差が無い場合、両者の不透明度に対して後述のステップＳ７の補正テーブルに基づく補正を加えても、内臓領域は骨領域に被ってしまう。そこで同図に示すように、内臓領域と骨領域の信号値の差を利用し、閾値Ｓｖを内臓領域と骨領域の信号値の境界に設定することで、骨領域の不透明度を効果的に減衰させることができ、内臓領域が鮮明に可視化される。
【0104】
一方、カラーマップＣｍａｐを調整しない場合には（ステップＳ５；「Ｎｏ」）、制御部１１は、上記したカラーマップＣｍａｐの調整を省略し、ステップＳ７へ移行する。
【0105】
【0106】
図８は、補正テーブルＳαを作成する処理を示すフローチャートである。
まず、制御部１１（幾何情報取得部２６ａ）は、断層画像Ｄｏ（ｚ）（０≦ｚ≦Ｓｚ－１）毎に、断層画像Ｄｏ（ｚ）の被写体形状の外郭を所定の幾何形状で近似し、幾何パラメータＰ（ｚ）を取得する（図８のステップＳ２１）。本実施の形態では、被写体形状の外郭を楕円で近似し、幾何パラメータＰ（ｚ）として楕円のパラメータ（楕円パラメータＰ（ｚ））を取得するものとする。
まず、制御部１１（幾何情報取得部２６ａ）は、断層画像Ｄｏ（ｚ）（０≦ｚ≦Ｓｚ－１）毎に、断層画像Ｄｏ（ｚ）の被写体形状の外郭を所定の幾何形状で近似し、幾何パラメータＰ（ｚ）を取得する（図８のステップＳ２１）。本実施の形態では、被写体形状の外郭を楕円で近似し、幾何パラメータＰ（ｚ）として楕円のパラメータ（楕円パラメータＰ（ｚ））を取得するものとする。
【0107】
図９は、断層画像Ｄｏ（ｚ）の楕円パラメータＰ（ｚ）を取得する処理を説明する図である。
制御部１１は、図９（ａ）に示すように、断層画像Ｄｏ（ｚ）に対して、信号値が所定の閾値以上の画素を「白」、閾値未満の画素を「黒」で塗りつぶして２値化し、被写体領域を抽出する。図９（ａ）の例では、２値化された断層画像Ｄｏ（ｚ）の「白」の領域が被写体領域である。閾値としては、例えば、空気以外を被写体領域として抽出する場合には、信号値ｖが１６ビットの場合は「－７００」、信号値ｖが８ビットの場合には「１０」などに設定する。
制御部１１は、図９（ａ）に示すように、断層画像Ｄｏ（ｚ）に対して、信号値が所定の閾値以上の画素を「白」、閾値未満の画素を「黒」で塗りつぶして２値化し、被写体領域を抽出する。図９（ａ）の例では、２値化された断層画像Ｄｏ（ｚ）の「白」の領域が被写体領域である。閾値としては、例えば、空気以外を被写体領域として抽出する場合には、信号値ｖが１６ビットの場合は「－７００」、信号値ｖが８ビットの場合には「１０」などに設定する。
【0108】
続いて、制御部１１は、図９（ｂ）に示すように、抽出した被写体領域（白領域）を全て含む最小の矩形領域Ｒ（ｚ）を特定し、図９（ｃ）のように、矩形領域Ｒ（ｚ）に内接する楕円Ｅ（ｚ）を、被写体形状の外郭を近似する楕円として求める。そして、制御部１１は、図９（ｃ）に示すように、楕円Ｅ（ｚ）の楕円パラメータＰ（ｚ）として、楕円Ｅ（ｚ）の中心座標（Ｃｘ（ｚ）、Ｃｙ（ｚ））、横方向のサイズＷ（ｚ）、及び縦方向のサイズＨ（ｚ）を以下のように算出する。尚、図９（ｃ）に示すように、矩形領域Ｒ（ｚ）の最小座標を（Ｘｍｉｎ、Ｙｍｉｎ）、最大座標を（Ｘｍａｘ、Ｙｍａｘ）とする。
【0109】
（式９）
中心座標Ｃｘ（ｚ）＝（Ｘｍａｘ＋Ｘｍｉｎ）／２
Ｃｙ（ｚ）＝（Ｙｍａｘ＋Ｙｍｉｎ）／２
横方向のサイズＷ（ｚ）＝Ｘｍａｘ－Ｘｍｉｎ＋１
縦方向のサイズＨ（ｚ）＝Ｙｍａｘ－Ｙｍｉｎ＋１
中心座標Ｃｘ（ｚ）＝（Ｘｍａｘ＋Ｘｍｉｎ）／２
Ｃｙ（ｚ）＝（Ｙｍａｘ＋Ｙｍｉｎ）／２
横方向のサイズＷ（ｚ）＝Ｘｍａｘ－Ｘｍｉｎ＋１
縦方向のサイズＨ（ｚ）＝Ｙｍａｘ－Ｙｍｉｎ＋１
【0110】
制御部１１は、全ての断層画像Ｄｏ（ｚ）（０≦ｚ≦Ｓｚ－１）について、被写体領域の外郭を近似する楕円Ｅ（ｚ）を求め、楕円パラメータＰ（ｚ）を算出する。
【0111】
続いて、制御部１１（幾何情報補正部２６ｂ）は、幾何形状（楕円）の近似精度が良好でない断層画像Ｄｏ（ｚ）について、幾何パラメータＰ（ｚ）（楕円パラメータＰ（ｚ））の補正を行う（図８のステップＳ２２）。
【0112】
図１０は、断層画像Ｄｏ（ｚ）の楕円パラメータＰ（ｚ）の補正について説明する図である。図１０（ａ）の断層画像Ｄｏ（ｚ）は、楕円近似が良好に行えない断層画像の例である。図１０（ａ）の断層画像は、被写体領域（幅方向）の両端が切れており、被写体領域が画像範囲に収まっていない（撮影時の条件によってこのような画像が得られる場合がある）。このような場合、被写体領域を正確に楕円近似できない。
【0113】
すなわち、図１０（ｂ）に示すように、断層画像から得られる楕円Ｅ（ｚ）の横方向のサイズＷ（ｚ）が、断層画像の横方向のサイズＳｘと同一となり（Ｗ（ｚ）＝Ｓｘ）、実際の被写体領域の幅（Ｓｘより大きい幅）を正確に反映したものでなくなる。このため、本実施の形態では、被写体領域（幅方向）の両端が切れた断層画像Ｄｏ（ｚ）について横方向のサイズＷ（ｚ）の補正を行う。
【0114】
具体的には、制御部１１は、まず、ステップＳ２１により取得された楕円の横方向のサイズＷ（ｚ）がＷ（ｚ）＝Ｓｘとなる１以上の断層画像Ｄｏ（ｚ）を、被写体領域の両端が切れた断層画像（幾何形状の近似精度が良好でない断層画像）として特定する。
【0115】
次に、制御部１１は、特定された断層画像Ｄｏ（ｚ）毎に、断層画像Ｄｏ（ｚ）とスライス位置が最も近く、かつ、ステップＳ２１により取得された楕円の横方向のサイズＷ（ｚ）がＷ（ｚ）＜Ｓｘを満たす断層画像（「断層画像Ｄｏ（ｚ’）」と表記）を特定する。
【0116】
断層画像Ｄｏ（ｚ’）は、楕円Ｅ（ｚ’）の横方向のサイズＷ（ｚ’）が画像サイズＳｘより小さいため、被写体領域（幅方向）が画像範囲に収まっている断層画像（幾何形状の近似精度が良好な断層画像）である。そして、制御部１１は、被写体領域の両端が切れた断層画像Ｄｏ（ｚ）から得られる楕円Ｅ（ｚ）の横方向のサイズＷ（ｚ）を、断層画像Ｄｏ（ｚ’） から得られる楕円Ｅ（ｚ’）の横方向のサイズＷ（ｚ’）及び縦方向のサイズＨ（ｚ’）の比と同一となるように、以下の式で補正する。
【0117】
（式１０）
Ｗ（ｚ）＝Ｈ（ｚ）・Ｗ（ｚ’）／Ｈ（ｚ’）
Ｗ（ｚ）＝Ｈ（ｚ）・Ｗ（ｚ’）／Ｈ（ｚ’）
【0118】
これにより、図１０（ｃ）に示すように、横方向のサイズＷ（ｚ）が、真の被写体領域の外郭を近似するような楕円Ｅ（ｚ）の横方向のサイズＷ（ｚ）に補正され、被写体領域（幅方向）が画像範囲に収まらない場合でも、楕円パラメータＰ（ｚ）を良好に算出することができる。
【0119】
以上、図８のステップＳ２１及びＳ２２の処理により、図１１に示すように、断層画像Ｄｏ（ｚ）毎に、被写体形状の外郭が楕円形状で近似され、全ての断層画像Ｄｏ（ｚ）（０≦ｚ≦Ｓｚ－１）について楕円パラメータＰ（ｚ）（０≦ｚ≦Ｓｚ－１）が得られる。
【0120】
そして、制御部１１（補正倍率算出部２６ｃ）は、断層画像Ｄｏ（ｚ）毎に、断層画像Ｄｏ（ｚ）に対応する幾何パラメータＰ（ｚ）（楕円パラメータＰ（ｚ））に基づいて、断層画像Ｄｏ（ｚ）の各画素（ｘ、ｙ、ｚ）に対応する不透明度の補正倍率Ｓα（ｘ、ｙ、ｚ）を算出し、補正テーブルＳαを作成する（図８のステップＳ２３）。
例えば、ｚ番目の断層画像Ｄｏ（ｚ）の各画素（ｘ、ｙ、ｚ）の補正倍率Ｓα（ｘ、ｙ、ｚ）は、以下のように算出される。
例えば、ｚ番目の断層画像Ｄｏ（ｚ）の各画素（ｘ、ｙ、ｚ）の補正倍率Ｓα（ｘ、ｙ、ｚ）は、以下のように算出される。
【0121】
（式１１）
Ｓα（ｘ、ｙ、ｚ）＝１－ｋ・ｒ（ｘ、ｙ、ｚ）・ａｍｐ
但し、０≦Ｓα（ｘ、ｙ、ｚ）≦１
ｒ（ｘ、ｙ、ｚ）＝｛（ｘ－Ｃｘ（ｚ）－ｘｏｆｆｓｅｔ）／（ａ・Ｗ（ｚ）／２）｝２＋｛（ｙ－Ｃｙ（ｚ）－ｙｏｆｆｓｅｔ）／（ｂ・Ｈ（ｚ）／２）｝２
Ｓα（ｘ、ｙ、ｚ）＝１－ｋ・ｒ（ｘ、ｙ、ｚ）・ａｍｐ
但し、０≦Ｓα（ｘ、ｙ、ｚ）≦１
ｒ（ｘ、ｙ、ｚ）＝｛（ｘ－Ｃｘ（ｚ）－ｘｏｆｆｓｅｔ）／（ａ・Ｗ（ｚ）／２）｝２＋｛（ｙ－Ｃｙ（ｚ）－ｙｏｆｆｓｅｔ）／（ｂ・Ｈ（ｚ）／２）｝２
【0122】
（式１１）において、ｋは減衰係数であり初期値は１．０である。ｋ＝１．０の場合、骨及び寝台が不可視となる。ａは横方向のサイズＷ（ｚ）の補正係数、ｂは縦方向のサイズＨ（ｚ）の補正係数である。補正しない場合はａ＝ｂ＝１．０に設定する。
ｘｏｆｆｓｅｔ、ｙｏｆｆｓｅｔはＸＹ面におけるオフセット値であり、初期値はｘｏｆｆｓｅｔ＝ｙｏｆｆｓｅｔ＝０（単位：ボクセル）である。
ｘｏｆｆｓｅｔ、ｙｏｆｆｓｅｔはＸＹ面におけるオフセット値であり、初期値はｘｏｆｆｓｅｔ＝ｙｏｆｆｓｅｔ＝０（単位：ボクセル）である。
【0123】
図１２は、（式１１）により算出される補正倍率Ｓα（ｘ、ｙ、ｚ）の分布パターン（同心楕円分布）を示す図である。図に示すように、楕円の中心（Ｃｘ（ｚ）、Ｃｙ（ｚ））の補正倍率Ｓαは１（最大）となり、楕円内部は中心から遠ざかるにつれ補正倍率Ｓαが小さくなる。具体的には、楕円の中心（Ｃｘ（ｚ）、Ｃｙ（ｚ））と横方向のサイズＷ（ｚ）及び縦方向のサイズＨ（ｚ）の比が同一となる各同心楕円上の補正倍率Ｓαが、各同心楕円の径（横方向のサイズ及び縦方向のサイズ）が大きい程小さくなる（径の２乗に反比例して小さくなる）。楕円外部の補正倍率Ｓαは一律に０とする。これにより、被写体領域の中心から外側にいくにつれ不透明度が減衰するので、骨領域（胸骨や頭蓋骨）や皮膚、外部の寝台、固定用治具等が透明化され、被写体領域の中心付近にある臓器等が鮮明に可視化される。
【0124】
尚、（式１１）のａｍｐはＸＹ平面の振幅倍率であり、初期値はａｍｐ＝１．０である。図１３に示すように、被写体中心部の補正倍率を平坦にし、飽和領域を広くしたい場合は、ａｍｐ＞１．０に設定する（但し、陰影感が弱くなる）。
【0125】
制御部１１は、全ての断層画像Ｄｏ（ｚ）について、補正倍率Ｓα（ｘ、ｙ、ｚ）（０≦ｘ≦Ｓｘ－１、０≦ｙ≦Ｓｙ－１、０≦ｚ≦Ｓｚ－１）を算出し、算出された補正倍率Ｓα（ｘ、ｙ、ｚ）を格納する３次元の補正テーブルＳαを作成する。尚、前述したように、補正倍率Ｓα（ｘ、ｙ、ｚ）は３次元のデータ配列であり、そのままデータテーブル（補正テーブル）として扱えるため、補正倍率と補正テーブルは実質的に同一である。
【0126】
【0127】
制御部１１（ボクセル作成部２７）は、ステップＳ１において取得された断層画像群Ｄｏ（又はステップＳ２において階調圧縮された断層画像群Ｄｏ）、ステップＳ４において取得されたカラーマップＣｍａｐ（又はステップＳ６において調整されたカラーマップＣｍａｐ）、ステップＳ７において作成された補正テーブルＳα、及びステップＳ３において指定された関心領域ＲＯＩに基づいて、ボクセル構造体Ｖ（ボクセルデータ）を作成する（図６のステップＳ８）。
【0128】
本実施の形態では、ボクセル構造体Ｖ（ボクセルデータ）は、２４ビットの色値（ＲＧＢ値）を保持するボクセル構造体Ｖｃと、８ビットの不透明度（α値）を保持するボクセル構造体Ｖαとの２つのボクセル構造体から構成されるものとする。このようにボクセル構造体を２つに分けて作成することで、色値及び不透明度に関するボクセル処理を個別に実行できるため、メモリアクセス等の向上が図られる。
【0129】
図１４は、ボクセル構造体Ｖ（Ｖｃ、Ｖα）を作成する処理を示すフローチャートである。
まず、制御部１１は、カラーマップＣｍａｐ（ｖ、ｎ）（０≦ｎ≦２）（＝カラーパレット）、及び関心領域ＲＯＩに基づいて、２４ビットのＲＧＢ値を保持するボクセル構造体Ｖｃ（ｘ、ｙ、ｚ、ｎ）（０≦Ｖ（ｘ、ｙ、ｚ、ｎ）≦２５５、０≦ｎ≦２、０≦ｘ≦Ｓｘ－１、０≦ｙ≦Ｓｙ－１、０≦ｚ≦Ｓｚ－１；解像度：Ｒｘｙ、Ｒｚ）を以下のように作成する（図１４のステップＳ８０１）。
まず、制御部１１は、カラーマップＣｍａｐ（ｖ、ｎ）（０≦ｎ≦２）（＝カラーパレット）、及び関心領域ＲＯＩに基づいて、２４ビットのＲＧＢ値を保持するボクセル構造体Ｖｃ（ｘ、ｙ、ｚ、ｎ）（０≦Ｖ（ｘ、ｙ、ｚ、ｎ）≦２５５、０≦ｎ≦２、０≦ｘ≦Ｓｘ－１、０≦ｙ≦Ｓｙ－１、０≦ｚ≦Ｓｚ－１；解像度：Ｒｘｙ、Ｒｚ）を以下のように作成する（図１４のステップＳ８０１）。
【0130】
（式１２）
（１）Ｘｓ＜ｘ＜Ｘｅ、Ｙｓ＜ｙ＜Ｙｅ、及びＺｓ＜ｚ＜Ｚｅを全て満たす場合（関心領域ＲＯＩ内の場合）
Ｖｃ（ｘ、ｙ、ｚ、ｎ）＝Ｃｍａｐ（Ｄｏ（ｘ、ｙ、ｚ）、ｎ） （０≦ｎ≦２）
（２）ｘ≦Ｘｓ、ｘ≧Ｘｅ、ｙ≦Ｙｓ、ｙ≧Ｙｅ、ｚ≦Ｚｓ、又はｚ≧Ｚｅのいずれかを満たす場合（関心領域ＲＯＩ外の場合）
Ｖｃ（ｘ、ｙ、ｚ、ｎ）＝０ （０≦ｎ≦２）
（１）Ｘｓ＜ｘ＜Ｘｅ、Ｙｓ＜ｙ＜Ｙｅ、及びＺｓ＜ｚ＜Ｚｅを全て満たす場合（関心領域ＲＯＩ内の場合）
Ｖｃ（ｘ、ｙ、ｚ、ｎ）＝Ｃｍａｐ（Ｄｏ（ｘ、ｙ、ｚ）、ｎ） （０≦ｎ≦２）
（２）ｘ≦Ｘｓ、ｘ≧Ｘｅ、ｙ≦Ｙｓ、ｙ≧Ｙｅ、ｚ≦Ｚｓ、又はｚ≧Ｚｅのいずれかを満たす場合（関心領域ＲＯＩ外の場合）
Ｖｃ（ｘ、ｙ、ｚ、ｎ）＝０ （０≦ｎ≦２）
【0131】
次に、制御部１１は、カラーマップＣｍａｐ（ｖ、３）（＝オパシティカーブ）、補正テーブルＳα、及び関心領域ＲＯＩに基づいて、８ビットの不透明度を保持するボクセル構造体Ｖα（ｘ、ｙ、ｚ）（０≦Ｖα（ｘ、ｙ、ｚ）≦２５５、０≦ｘ≦Ｓｘ－１、０≦ｙ≦Ｓｙ－１、０≦ｚ≦Ｓｚ－１；解像度：Ｒｘｙ、Ｒｚ）を以下のように作成する（図１４のステップＳ８０２）。
【0132】
（式１３）
（１）Ｘｓ＜ｘ＜Ｘｅ、Ｙｓ＜ｙ＜Ｙｅ、及びＺｓ＜ｚ＜Ｚｅを全て満たす場合（関心領域ＲＯＩ内の場合）
Ｖα（ｘ、ｙ、ｚ）＝Ｃｍａｐ（Ｄｏ（ｘ、ｙ、ｚ）、３）
Ｖα（ｘ、ｙ、ｚ）＝Ｖα（ｘ、ｙ、ｚ）・Ｓα（ｘ、ｙ、ｚ）
（２）ｘ≦Ｘｓ、ｘ≧Ｘｅ、ｙ≦Ｙｓ、ｙ≧Ｙｅ、ｚ≦Ｚｓ、又はｚ≧Ｚｅのいずれかを満たす場合（関心領域ＲＯＩ外の場合）
Ｖα（ｘ、ｙ、ｚ）＝０
（１）Ｘｓ＜ｘ＜Ｘｅ、Ｙｓ＜ｙ＜Ｙｅ、及びＺｓ＜ｚ＜Ｚｅを全て満たす場合（関心領域ＲＯＩ内の場合）
Ｖα（ｘ、ｙ、ｚ）＝Ｃｍａｐ（Ｄｏ（ｘ、ｙ、ｚ）、３）
Ｖα（ｘ、ｙ、ｚ）＝Ｖα（ｘ、ｙ、ｚ）・Ｓα（ｘ、ｙ、ｚ）
（２）ｘ≦Ｘｓ、ｘ≧Ｘｅ、ｙ≦Ｙｓ、ｙ≧Ｙｅ、ｚ≦Ｚｓ、又はｚ≧Ｚｅのいずれかを満たす場合（関心領域ＲＯＩ外の場合）
Ｖα（ｘ、ｙ、ｚ）＝０
【0133】
（式１３）のＶα（ｘ、ｙ、ｚ）＝Ｖα（ｘ、ｙ、ｚ）・Ｓα（ｘ、ｙ、ｚ）は、本開示の要旨に係る処理であり、これにより、カラーマップＣｍａｐに定義された不透明度に対して補正倍率Ｓαが乗算され、不透明度が補正される。特に、補正倍率Ｓα（ｘ、ｙ、z）は、被写体領域の外郭を近似した楕円の中心から遠ざかるにつれその倍率が小さくなるので、被写体領域の中心から離れたボクセルほど不透明が減衰する。これにより、骨領域（胸骨や頭蓋骨）や皮膚、外部の寝台、固定用治具等が透明化され、被写体領域の中心付近にある臓器等が鮮明に可視化されたレンダリング像を得ることができるようになる。
【0134】
尚、メモリアクセスの向上等を目的に、３２ビットのボクセル構造体Ｖを色値を保持するボクセル構造体Ｖｃと不透明度を保持するボクセル構造体Ｖαとに分けているが、ボクセル構造体を分けなくてもよい。この場合、制御部１１は、以下のように、３２ビットの色値及び不透明度を保持する１つのボクセル構造体Ｖ（ｘ、ｙ、ｚ、ｎ）（０≦Ｖ（ｘ、ｙ、ｚ、ｎ）≦２５５、０≦ｎ≦３、０≦ｘ≦Ｓｘ－１、０≦ｙ≦Ｓｙ－１、０≦ｚ≦Ｓｚ－１；解像度：Ｒｘｙ、Ｒｚ）を作成する。
【0135】
（式１４）
（１）Ｘｓ＜ｘ＜Ｘｅ、Ｙｓ＜ｙ＜Ｙｅ、及びＺｓ＜ｚ＜Ｚｅを全て満たす場合（関心領域ＲＯＩ内の場合）
Ｖ（ｘ、ｙ、ｚ、ｎ）＝Ｃｍａｐ（Ｄｏ（ｘ、ｙ、ｚ）、ｎ） （０≦ｎ≦３）
Ｖ（ｘ、ｙ、ｚ、３）＝Ｖ（ｘ、ｙ、ｚ、３）・Ｓα（ｘ、ｙ、ｚ）
（２）ｘ≦Ｘｓ、ｘ≧Ｘｅ、ｙ≦Ｙｓ、ｙ≧Ｙｅ、ｚ≦Ｚｓ、又はｚ≧Ｚｅのいずれかを満たす場合（関心領域ＲＯＩ外の場合）
Ｖ（ｘ、ｙ、ｚ、ｎ）＝０ （０≦ｎ≦３）
（１）Ｘｓ＜ｘ＜Ｘｅ、Ｙｓ＜ｙ＜Ｙｅ、及びＺｓ＜ｚ＜Ｚｅを全て満たす場合（関心領域ＲＯＩ内の場合）
Ｖ（ｘ、ｙ、ｚ、ｎ）＝Ｃｍａｐ（Ｄｏ（ｘ、ｙ、ｚ）、ｎ） （０≦ｎ≦３）
Ｖ（ｘ、ｙ、ｚ、３）＝Ｖ（ｘ、ｙ、ｚ、３）・Ｓα（ｘ、ｙ、ｚ）
（２）ｘ≦Ｘｓ、ｘ≧Ｘｅ、ｙ≦Ｙｓ、ｙ≧Ｙｅ、ｚ≦Ｚｓ、又はｚ≧Ｚｅのいずれかを満たす場合（関心領域ＲＯＩ外の場合）
Ｖ（ｘ、ｙ、ｚ、ｎ）＝０ （０≦ｎ≦３）
【0136】
上記した（式１２）～（式１４）に示すボクセル構造体を作成する処理は、断層画像群Ｄｏの信号値が１６ビット（（式１）参照）の場合は、１６ビット対応のカラーマップＣｍａp（（式４）参照）を用いて実行され、断層画像群Ｄｏの信号値が８ビット（（式２）参照）に階調圧縮されている場合は、８ビット対応のカラーマップＣｍａp（（式５）参照）を用いて実行される。
【0137】
制御部１１は、ボクセル構造体Ｖ（Ｖｃ、Ｖα）を作成した後、必要に応じてボクセル構造体Ｖ（Ｖｃ、Ｖα）に対して、公知のスムーシング処理を施したり、公知の陰影処理を施しても良い。
【0138】
続いて、制御部１１（第１レンダリング部２８ａ又は第２レンダリング部２８ｂ）は、ボクセル構造体Ｖに基づいて、ボリュームレンダリング像を生成する（図６のステップＳ９）。以降、ＣＰＵ（第１レンダリング部２８ａ）により実行されるレイキャスティング法によるレンダリング処理と、ＧＰＵ（第２レンダリング部２８ｂ）により実行される３Ｄテクスチャマッピング法によるレンダリング処理について順に説明する。
【0139】
最初に、図１５～図２０を参照しながら、レイキャスティング法によるレンダリング処理について説明する。
図１５は、レイキャスティング法の概要を示す図である。図１５（ａ）に示すように、投影面（レンダリング像）の各点からボクセルデータに対して任意の方向にレイを発射し、最深点（レイ方向で光が届く最も深い点）よりその軌跡を逆にたどり、投影面上の点の輝度値を算出する。
また、図１５（ｂ）に示すように、計算を簡単にするため、ボクセルデータに座標変換を施し、投影面の各点よりＺ軸負方向の最深点までレイを進め、同時に累積輝度計算を行う。
図１５は、レイキャスティング法の概要を示す図である。図１５（ａ）に示すように、投影面（レンダリング像）の各点からボクセルデータに対して任意の方向にレイを発射し、最深点（レイ方向で光が届く最も深い点）よりその軌跡を逆にたどり、投影面上の点の輝度値を算出する。
また、図１５（ｂ）に示すように、計算を簡単にするため、ボクセルデータに座標変換を施し、投影面の各点よりＺ軸負方向の最深点までレイを進め、同時に累積輝度計算を行う。
【0140】
図１６は、レンダリング処理を示すフローチャートである。まず、制御部１１は、座標変換パラメータを設定する（ステップＳ３１）。後述する探索制御マスクの算出（ステップＳ３２）、レイキャスティング処理（ステップＳ３３）の各処理で、視点座標系からボクセル座標系への座標変換処理を行う。そのため、制御部１１は、以下のように、座標変換処理に共通する座標変換パラメータを事前に設定・算出しておく。
【0141】
＜座標変換パラメータ＞
回転パラメータ行列Ｒ：
Ｒ＝［Ｒ１１ Ｒ１２ Ｒ１３；
Ｒ２１ Ｒ２２ Ｒ２３；
Ｒ３１ Ｒ３２ Ｒ３３］
（ボクセル座標系から視点座標系への座標変換を行うための３×３の行列の逆行列、ＧＵＩ側はボクセル座標系から視点座標系への座標変換を指示するが、レンダリング側は視点座標系からボクセル座標系に座標変換を行う。）
Ｘ軸方向のオフセットＸｏｆｆ
Ｙ軸方向のオフセットＹｏｆｆ
Ｚ軸方向のオフセットＺｏｆｆ
Ｘ軸方向のＲＯＩ：Ｘｓ－Ｘｅ（０≦Ｘｓ、Ｘｅ≦Ｓｘ－１）
Ｙ軸方向のＲＯＩ：Ｙｓ－Ｙｅ（０≦Ｙｓ、Ｙｅ≦Ｓｙ－１）
Ｚ軸方向のＲＯＩ：Ｚｓ－Ｚｅ（０≦Ｚｓ、Ｚｅ≦Ｓｚ－１）
透視変換パラメータ：原点からＺ軸方向の距離Ｄｉｓｔ（正の実数値、平行投影の場合はＤｉｓｔ＝０）
拡大縮小倍率Ｓｃａｌｅ（ＸＹＺ軸方向で同一）
Ｚ方向変倍率Ｓｃｚ＝Ｒｘｙ／Ｒｚ
Ｚ方向変倍後のＺ方向サイズＳｚ’＝Ｓｚ・Ｓｃｚ
座標変換サブサンプル・オフセット：Ｘ軸方向ｄｘ、Ｙ軸方向ｄｙ、Ｚ軸方向ｄｚ
回転パラメータ行列Ｒ：
Ｒ＝［Ｒ１１ Ｒ１２ Ｒ１３；
Ｒ２１ Ｒ２２ Ｒ２３；
Ｒ３１ Ｒ３２ Ｒ３３］
（ボクセル座標系から視点座標系への座標変換を行うための３×３の行列の逆行列、ＧＵＩ側はボクセル座標系から視点座標系への座標変換を指示するが、レンダリング側は視点座標系からボクセル座標系に座標変換を行う。）
Ｘ軸方向のオフセットＸｏｆｆ
Ｙ軸方向のオフセットＹｏｆｆ
Ｚ軸方向のオフセットＺｏｆｆ
Ｘ軸方向のＲＯＩ：Ｘｓ－Ｘｅ（０≦Ｘｓ、Ｘｅ≦Ｓｘ－１）
Ｙ軸方向のＲＯＩ：Ｙｓ－Ｙｅ（０≦Ｙｓ、Ｙｅ≦Ｓｙ－１）
Ｚ軸方向のＲＯＩ：Ｚｓ－Ｚｅ（０≦Ｚｓ、Ｚｅ≦Ｓｚ－１）
透視変換パラメータ：原点からＺ軸方向の距離Ｄｉｓｔ（正の実数値、平行投影の場合はＤｉｓｔ＝０）
拡大縮小倍率Ｓｃａｌｅ（ＸＹＺ軸方向で同一）
Ｚ方向変倍率Ｓｃｚ＝Ｒｘｙ／Ｒｚ
Ｚ方向変倍後のＺ方向サイズＳｚ’＝Ｓｚ・Ｓｃｚ
座標変換サブサンプル・オフセット：Ｘ軸方向ｄｘ、Ｙ軸方向ｄｙ、Ｚ軸方向ｄｚ
【0142】
制御部１１は、上記した回転パラメータ行列Ｒに対して、初期値として単位行列を設定する。すなわち、Ｒ１１＝１、Ｒ１２＝０、Ｒ１３＝０、Ｒ２１＝０、Ｒ２２＝１、Ｒ２３＝０、Ｒ３１＝０、Ｒ３２＝０、Ｒ３３＝１と設定する。
そして、ＧＵＩの指示に従い、Ｘ軸中心回転Ｒｘ、Ｙ軸中心回転Ｒｙ、Ｚ軸中心回転Ｒｚ（角度単位：ラジアン）のいずれかを逐次指定し、以下のように、各々回転行列Ａを生成して回転パラメータ行列Ｒに右から乗算して、回転パラメータ行列Ｒを更新する。これにより、ＧＵＩの指示により生成されるボクセル座標系から視点座標系への回転行列の逆行列が算出される。
そして、ＧＵＩの指示に従い、Ｘ軸中心回転Ｒｘ、Ｙ軸中心回転Ｒｙ、Ｚ軸中心回転Ｒｚ（角度単位：ラジアン）のいずれかを逐次指定し、以下のように、各々回転行列Ａを生成して回転パラメータ行列Ｒに右から乗算して、回転パラメータ行列Ｒを更新する。これにより、ＧＵＩの指示により生成されるボクセル座標系から視点座標系への回転行列の逆行列が算出される。
【0143】
回転行列Ａを
Ａ＝［Ａ１１ Ａ１２ Ａ１３；
Ａ２１ Ａ２２ Ａ２３；
Ａ３１ Ａ３２ Ａ３３］
とすると、
Ｘ軸中心回転Ｒｘの場合の回転行列Ａの各要素は、
Ａ１１＝１、Ａ１２＝０、Ａ１３＝０
Ａ２１＝０、Ａ２２＝ｃｏｓＲｘ、Ａ２３＝ｓｉｎＲｘ
Ａ３１＝０、Ａ３２＝ｓｉｎＲｘ、Ａ３３＝ｃｏｓＲｘ
Ｙ軸中心回転Ｒｙの場合の回転行列Ａの各要素は、
Ａ１１＝ｃｏｓＲｙ、Ａ１２＝０、Ａ１３＝ｓｉｎＲｙ
Ａ２１＝０、Ａ２２＝１、Ａ２３＝０
Ａ３１＝－ｓｉｎＲｙ、Ａ３２＝０、Ａ３３＝ｃｏｓＲｙ
Ｚ軸中心回転Ｒｚの場合の回転行列Ａの各要素は、
Ａ１１＝ｃｏｓＲｚ、Ａ１２＝ｓｉｎＲｚ、Ａ１３＝０
Ａ２１＝－ｓｉｎＲｚ、Ａ２２＝ｃｏｓＲｚ、Ａ２３＝０
Ａ３１＝０、Ａ３２＝０、Ａ３３＝１
となる。
回転パラメータ行列Ｒは、Ｒ←Ｒ×Ａと更新される。
Ａ＝［Ａ１１ Ａ１２ Ａ１３；
Ａ２１ Ａ２２ Ａ２３；
Ａ３１ Ａ３２ Ａ３３］
とすると、
Ｘ軸中心回転Ｒｘの場合の回転行列Ａの各要素は、
Ａ１１＝１、Ａ１２＝０、Ａ１３＝０
Ａ２１＝０、Ａ２２＝ｃｏｓＲｘ、Ａ２３＝ｓｉｎＲｘ
Ａ３１＝０、Ａ３２＝ｓｉｎＲｘ、Ａ３３＝ｃｏｓＲｘ
Ｙ軸中心回転Ｒｙの場合の回転行列Ａの各要素は、
Ａ１１＝ｃｏｓＲｙ、Ａ１２＝０、Ａ１３＝ｓｉｎＲｙ
Ａ２１＝０、Ａ２２＝１、Ａ２３＝０
Ａ３１＝－ｓｉｎＲｙ、Ａ３２＝０、Ａ３３＝ｃｏｓＲｙ
Ｚ軸中心回転Ｒｚの場合の回転行列Ａの各要素は、
Ａ１１＝ｃｏｓＲｚ、Ａ１２＝ｓｉｎＲｚ、Ａ１３＝０
Ａ２１＝－ｓｉｎＲｚ、Ａ２２＝ｃｏｓＲｚ、Ａ２３＝０
Ａ３１＝０、Ａ３２＝０、Ａ３３＝１
となる。
回転パラメータ行列Ｒは、Ｒ←Ｒ×Ａと更新される。
【0144】
以上の座標変換処理は、探索制御マスクの算出（ステップＳ３２）、レイキャスティング処理（ステップＳ３３）の各処理の中で逐次実行される。
【0145】
次に、図１７のフローチャートを参照して、図１６のステップＳ３２における探索制御マスクの算出処理について説明する。
探索制御マスクＭ（ｘ、ｙ、ｌ）（０≦ｘ≦Ｓｘ－１、０≦ｙ≦Ｓｙ－１、０≦ｌ≦Ｌ－１）は、算出するレンダリング像の画素（ｘ、ｙ）毎に先頭の不透明ボクセルのＺ座標を記憶したものである。座標変換時にサブサンプリング（ボリュームレンダリング像のジャギーを除去するため、アンチエリアシングともいわれる）を行う場合は、制御部１１は、以下の処理をＬ回（０≦ｌ≦Ｌ－１）実行する。
探索制御マスクＭ（ｘ、ｙ、ｌ）（０≦ｘ≦Ｓｘ－１、０≦ｙ≦Ｓｙ－１、０≦ｌ≦Ｌ－１）は、算出するレンダリング像の画素（ｘ、ｙ）毎に先頭の不透明ボクセルのＺ座標を記憶したものである。座標変換時にサブサンプリング（ボリュームレンダリング像のジャギーを除去するため、アンチエリアシングともいわれる）を行う場合は、制御部１１は、以下の処理をＬ回（０≦ｌ≦Ｌ－１）実行する。
【0146】
まず、制御部１１は、処理回数ｌ＝０とし、サブサンプル・オフセット値をｄｘ＝ｄｙ＝ｄｚ＝０と初期化する（ステップＳ４１）。
【0147】
続いて、制御部１１は、ｘ＝２ｍ及びｙ＝２ｎ（ｍ、ｎは整数）の値をとる画素（ｘ、ｙ）に対して、Ｚｓ＝Ｓｚ’－１から座標変換を行いながら先頭の不透明ボクセルを探索し、そのＺ座標をＭ（ｘ、ｙ、ｌ）に記録する。ｘ＝Ｓｘ－１又はｙ＝Ｓｙ－１のいずれかの値をもつ画素（ｘ、ｙ）に対しても同様の探索を行う（ステップＳ４２）。これにより、図１８（ａ）のように、縦横偶数番目の画素について不透明ボクセルのＺ座標が算出される。不透明度ボクセルを探索する処理は、図１９にて後述する。
【0148】
続いて、制御部１１は、ｘ＝２ｍ＋１（ｘ＜Ｓｘ－１）及びｙ＝２ｎ（ｍ、ｎは整数）の値をとる画素（ｘ、ｙ）に対して、Ｍ（ｘ－１、ｙ、ｌ）＝Ｍ（ｘ＋１、ｙ、ｌ）の場合、Ｍ（ｘ、ｙ、ｌ）＝Ｍ（ｘ－１、ｙ、ｌ）を与え、Ｍ（ｘ－１、ｙ、ｌ）≠Ｍ（ｘ＋１、ｙ、ｌ）の場合、上記同様Ｚｓ＝Ｓｚ’－１から座標変換を行いながら先頭の不透明ボクセルを探索（図１９参照）し、そのＺ座標をＭ（ｘ、ｙ、ｌ）に記録する（ステップＳ４３）。これにより、図１８（ｂ）のように、Ｘ方向に奇数番目の画素について不透明ボクセルのＺ座標が算出される。隣接する左右画素の不透明ボクセルのＺ座標が同一の場合、そのＺ座標で補間する（不透明ボクセルの探索は行わない）。
【0149】
続いて、制御部１１は、全てのｘ（０≦ｘ≦Ｓｘ－１）及びｙ＝２ｎ＋１（ｎは整数）の値をとる画素（ｘ、ｙ）に対して、Ｍ（ｘ、ｙ－１、ｌ）＝Ｍ（ｘ、ｙ＋１、ｌ）の場合、Ｍ（ｘ、ｙ、ｌ）＝Ｍ（ｘ、ｙ－１、ｌ）を与え、Ｍ（ｘ、ｙ－１、ｌ）≠Ｍ（ｘ、ｙ＋１、ｌ）の場合、上記同様Ｚｓ＝Ｓｚ’－１から座標変換を行いながら先頭の不透明ボクセルを探索（図１９参照）し、そのＺ座標をＭ（ｘ、ｙ、ｌ）に記録する（ステップＳ４４）。これにより、図１８（ｃ）のように、Ｙ方向に奇数番目の画素について不透明ボクセルのＺ座標が算出される。隣接する上下画素の不透明ボクセルのＺ座標が同一の場合、そのＺ座標で補間する（不透明ボクセルの探索は行わない）。
【0150】
続いて、制御部１１は、処理回数をｌ←ｌ＋１に更新し、サブサンプル・オフセット値を、ｄｘ←ｄｘ＋１／Ｌ、ｄｙ←ｄｙ＋１／Ｌ、ｄｚ←ｄｚ＋１／Ｌに更新する（ステップＳ４５）。
制御部１１は、処理回数ｌ＝Ｌ－１となるまで（ステップＳ４６；Ｙｅｓ）、上記したステップＳ４２～Ｓ４５の処理を繰り返す。
制御部１１は、処理回数ｌ＝Ｌ－１となるまで（ステップＳ４６；Ｙｅｓ）、上記したステップＳ４２～Ｓ４５の処理を繰り返す。
【0151】
次に、図１９のフローチャートを参照して、探索制御マスク算出処理（図１７）において実行される不透明度ボクセルの探索処理について説明する。
まず、制御部１１は、ｚをｚ＝Ｚｓと初期化し、探索対象画素（ｘ、ｙ）を入力する（ステップＳ５１）。続いて、制御部１１は、３次元座標（ｘ、ｙ、ｚ）に対して座標変換を行い、ボクセルα値（ｎ＝３）を算出する（ステップＳ５２）。座標変換を行い、ボクセルα値を算出する処理（座標変換処理）は後述する。
まず、制御部１１は、ｚをｚ＝Ｚｓと初期化し、探索対象画素（ｘ、ｙ）を入力する（ステップＳ５１）。続いて、制御部１１は、３次元座標（ｘ、ｙ、ｚ）に対して座標変換を行い、ボクセルα値（ｎ＝３）を算出する（ステップＳ５２）。座標変換を行い、ボクセルα値を算出する処理（座標変換処理）は後述する。
【0152】
ステップＳ５２において算出されたボクセルα値がα＝０の場合、制御部１１は、ｚ←ｚ－ｍ（例えばｍ＝８）に更新する（ステップＳ５３）。そして、ｚ＜０の場合、制御部１１は、Ｍ（ｘ、ｙ、ｌ）＝－１に設定し（ステップＳ５４）、処理を終了する。ｚ≧０の場合、制御部１１は、ステップＳ５２に戻る。
【0153】
一方、ステップＳ５２において算出されたボクセルα値がα＞０の場合、制御部１１は、まず、ｉをｉ＝０と初期化する（ステップＳ５５）。続いて、制御部１１は、ｉ←ｉ＋１に更新し（ステップＳ５６）、ｉ＜ｍかつｚ＋ｉ＜Ｚｓの場合、３次元座標（ｘ、ｙ、ｚ＋ｉ）に対して座標変換を行い、ボクセルα値を算出する（ステップＳ５７）。座標変換を行い、ボクセルα値を算出する処理（座標変換処理）は後述する。
【0154】
一方、ｉ≧ｍまたはｚ＋ｉ≧Ｚｓの場合、制御部１１は、Ｍ（ｘ、ｙ、ｌ）＝ｚ＋ｉ－１に設定し（ステップＳ５８）、処理を終了する。
【0155】
また、ステップＳ５７において算出されたボクセルα値がα＝０の場合、制御部１１は、Ｍ（ｘ、ｙ、ｌ）＝ｚ＋ｉ－１に設定し（ステップＳ５８）、処理を終了する。
【0156】
一方、ステップＳ５７において算出されたボクセルα値がα＞０の場合、ステップＳ５６に戻り、ステップＳ５８において不透明ボクセルのＺ座標が設定されるまで、ステップＳ５６～Ｓ５７の処理を繰り返す。
【0157】
不透明度ボクセル探索処理（図１９）のステップＳ５２及びステップＳ５７において実行される座標変換処理について説明する。
座標変換処理は、視点座標系をボクセル座標系に変換する処理であり、ＧＵＩ側の変換処理とは逆になる。ＧＵＩ側では関心領域ＲＯＩによるクリッピング、スケーリング、Ｚ方向変倍処理、オフセット（ＸＹＺ方向同時）、回転、透視変換の順に行うものと仮定し、制御部１１は、与えられた視点座標系の３次元座標値（ｘ、ｙ、ｚ）（整数値）に対応する８ビットのボクセル構造体Ｖα（ｘ、ｙ、ｚ）又は２４ビットのボクセル構造体Ｖｃ（ｘ、ｙ、ｚ、ｎ）の実数の座標値（ｘｒ、ｙｒ、ｚｒ）を以下のように算出する。
座標変換処理は、視点座標系をボクセル座標系に変換する処理であり、ＧＵＩ側の変換処理とは逆になる。ＧＵＩ側では関心領域ＲＯＩによるクリッピング、スケーリング、Ｚ方向変倍処理、オフセット（ＸＹＺ方向同時）、回転、透視変換の順に行うものと仮定し、制御部１１は、与えられた視点座標系の３次元座標値（ｘ、ｙ、ｚ）（整数値）に対応する８ビットのボクセル構造体Ｖα（ｘ、ｙ、ｚ）又は２４ビットのボクセル構造体Ｖｃ（ｘ、ｙ、ｚ、ｎ）の実数の座標値（ｘｒ、ｙｒ、ｚｒ）を以下のように算出する。
【0158】
制御部１１は、視点座標系の座標値（ｘ、ｙ、ｚ）（整数値）を次のように実数値（ｘｘ、ｙｙ、ｚｚ）に変換する。
ｘｘ＝ｘ－Ｓｘ／２＋ｄｘ
ｙｙ＝ｙ－Ｓｙ／２＋ｄｙ
ｚｚ＝ｚ－Ｓｚ’／２＋ｄｚ
ｘｘ＝ｘ－Ｓｘ／２＋ｄｘ
ｙｙ＝ｙ－Ｓｙ／２＋ｄｙ
ｚｚ＝ｚ－Ｓｚ’／２＋ｄｚ
【0159】
続いて、原点からＺ軸方向の距離Ｄｉｓｔ＞０の場合、制御部１１は、以下のように透視変換を行う。
ｚｚ’＝Ｄｉｓｔ・ｚｚ／（Ｄｉｓｔ＋ｚｚ）
ｘｘ’＝ｘｘ・（Ｄｉｓｔ－ｚｚ’）／Ｄｉｓｔ
ｙｙ’＝ｙｙ・（Ｄｉｓｔ－ｚｚ’）／Ｄｉｓｔ
透視変換後の（ｘｘ’、ｙｙ’、ｚｚ’）を（ｘｘ、ｙｙ、ｚｚ）とする。
ｚｚ’＝Ｄｉｓｔ・ｚｚ／（Ｄｉｓｔ＋ｚｚ）
ｘｘ’＝ｘｘ・（Ｄｉｓｔ－ｚｚ’）／Ｄｉｓｔ
ｙｙ’＝ｙｙ・（Ｄｉｓｔ－ｚｚ’）／Ｄｉｓｔ
透視変換後の（ｘｘ’、ｙｙ’、ｚｚ’）を（ｘｘ、ｙｙ、ｚｚ）とする。
【0160】
続いて、制御部１１は、以下のように、回転パラメータ行列Ｒを用いて回転処理を行う。
ｘｘ’＝Ｒ１１・ｘｘ＋Ｒ１２・ｙｙ＋Ｒ１３・ｚｚ
ｙｙ’＝Ｒ２１・ｘｘ＋Ｒ２２・ｙｙ＋Ｒ２３・ｚｚ
ｚｚ’＝Ｒ３１・ｘｘ＋Ｒ３２・ｙｙ＋Ｒ３３・ｚｚ
回転処理後の（ｘｘ’、ｙｙ’、ｚｚ’）を（ｘｘ、ｙｙ、ｚｚ）とする。
ｘｘ’＝Ｒ１１・ｘｘ＋Ｒ１２・ｙｙ＋Ｒ１３・ｚｚ
ｙｙ’＝Ｒ２１・ｘｘ＋Ｒ２２・ｙｙ＋Ｒ２３・ｚｚ
ｚｚ’＝Ｒ３１・ｘｘ＋Ｒ３２・ｙｙ＋Ｒ３３・ｚｚ
回転処理後の（ｘｘ’、ｙｙ’、ｚｚ’）を（ｘｘ、ｙｙ、ｚｚ）とする。
【0161】
制御部１１は、スケーリング、Ｚ方向変倍処理、及びオフセット（ＸＹＺ方向同時）を同時に行い、次のようにボクセル構造体の対応するボクセルの座標値（ｘｒ、ｙｒ、ｚｒ）（実数値）を取得する。
ｘｒ＝ｘｘ／Ｓｃａｌｅ＋Ｓｘ／２－Ｘｏｆｆ
ｙｒ＝ｙｙ／Ｓｃａｌｅ＋Ｓｙ／２－Ｙｏｆｆ
ｚｒ＝ｚｚ／Ｓｃａｌｅ／Ｓｃｚ＋Ｓｚ／２－Ｚｏｆｆ
ｘｒ＝ｘｘ／Ｓｃａｌｅ＋Ｓｘ／２－Ｘｏｆｆ
ｙｒ＝ｙｙ／Ｓｃａｌｅ＋Ｓｙ／２－Ｙｏｆｆ
ｚｒ＝ｚｚ／Ｓｃａｌｅ／Ｓｃｚ＋Ｓｚ／２－Ｚｏｆｆ
【0162】
続いて、制御部１１は、算出したボクセルの座標値（ｘｒ、ｙｒ、ｚｒ）（実数値）に対して、小数点以下を切り捨て整数化した座標値を（ｘｉ、ｙｉ、ｚｉ）（整数値）とし、切り捨てた小数点以下の端数を（ｗｘ、ｗｙ、ｗｚ）（０≦ｗｘ、ｗｙ、ｗｚ＜１）とする（すなわち、ｘｒ＝ｘｉ＋ｗｘ、ｙｒ＝ｙｉ＋ｗｙ、ｚｒ＝ｚｉ＋ｗｚ）。
【0163】
そして、制御部１１は、Ｘ軸・Ｙ軸・Ｚ軸方向のＲＯＩを考慮して、視点座標系の３次元座標値（ｘ、ｙ、ｚ）（整数値）に対応するボクセルのα値（Ｖα）を次のように３通りの場合に分けて決定する。
【0164】
１）ｘｉ＜Ｘｓ、ｘｉ＞Ｘｅ、ｙｉ＜Ｙｓ、ｙｉ＞Ｙｅ、ｚｉ＜Ｚｓ、又はｚｉ＞Ｚｅのいずれかを満たす場合（クリッピング範囲）
Ｖα＝０
２）上記１）の条件を満たさない場合において、ｘｉ＋１＞Ｘｅ、ｙｉ＋１＞Ｙｅ、又はｚｉ＋１＞Ｚｅのいずれかを満たす場合（補間しない）
Ｖα＝Ｖα（ｘｉ、ｙｉ、ｚｉ）
Ｖα＝０
２）上記１）の条件を満たさない場合において、ｘｉ＋１＞Ｘｅ、ｙｉ＋１＞Ｙｅ、又はｚｉ＋１＞Ｚｅのいずれかを満たす場合（補間しない）
Ｖα＝Ｖα（ｘｉ、ｙｉ、ｚｉ）
【0165】
３）上記１）２）の条件を満たさない場合（補間する）
Ｖα＝（１－ｗｚ）（１－ｗｙ）（１－ｗｘ）・Ｖα（ｘｉ、ｙｉ、ｚｉ）＋（１－ｗｚ）（１－ｗｙ）・ｗｘ・Ｖα（ｘｉ＋１、ｙｉ、ｚｉ）＋（１－ｗｚ）・ｗｙ・（１－ｗｘ）・Ｖα（ｘｉ、ｙｉ＋１、ｚｉ）＋（１－ｗｚ）・ｗｙ・ｗｘ・Ｖα（ｘｉ＋１、ｙｉ＋１、ｚｉ）＋ｗｚ・（１－ｗｙ）（１－ｗｘ）・Ｖα（ｘｉ、ｙｉ、ｚｉ＋１）＋ｗｚ・（１－ｗｙ）・ｗｘ・Ｖα（ｘｉ＋１、ｙｉ、ｚｉ＋１）＋ｗｚ・ｗｙ・（１－ｗｘ）・Ｖα（ｘｉ、ｙｉ＋１、ｚｉ＋１）＋ｗｚ・ｗｙ・ｗｘ・Ｖα（ｘｉ＋１、ｙｉ＋１、ｚｉ＋１）
Ｖα＝（１－ｗｚ）（１－ｗｙ）（１－ｗｘ）・Ｖα（ｘｉ、ｙｉ、ｚｉ）＋（１－ｗｚ）（１－ｗｙ）・ｗｘ・Ｖα（ｘｉ＋１、ｙｉ、ｚｉ）＋（１－ｗｚ）・ｗｙ・（１－ｗｘ）・Ｖα（ｘｉ、ｙｉ＋１、ｚｉ）＋（１－ｗｚ）・ｗｙ・ｗｘ・Ｖα（ｘｉ＋１、ｙｉ＋１、ｚｉ）＋ｗｚ・（１－ｗｙ）（１－ｗｘ）・Ｖα（ｘｉ、ｙｉ、ｚｉ＋１）＋ｗｚ・（１－ｗｙ）・ｗｘ・Ｖα（ｘｉ＋１、ｙｉ、ｚｉ＋１）＋ｗｚ・ｗｙ・（１－ｗｘ）・Ｖα（ｘｉ、ｙｉ＋１、ｚｉ＋１）＋ｗｚ・ｗｙ・ｗｘ・Ｖα（ｘｉ＋１、ｙｉ＋１、ｚｉ＋１）
【0166】
尚、後述するレイキャスティング処理（図２０参照）のようにボクセルのＲＧＢ値を取得する場合は、制御部１１は、視点座標系の３次元座標値（ｘ、ｙ、ｚ）（整数値）に対応するボクセルのＲＧＢ値（Ｖｃ（ｎ））（０≦ｎ≦２）を次のように３通りの場合に分けて決定する。
【0167】
１）ｘｉ＜Ｘｓ、ｘｉ＞Ｘｅ、ｙｉ＜Ｙｓ、ｙｉ＞Ｙｅ、ｚｉ＜Ｚｓ、又はｚｉ＞Ｚｅのいずれかを満たす場合（クリッピング範囲）
Ｖｃ（ｎ）＝０ （０≦ｎ≦２）
２）上記１）の条件を満たさない場合において、ｘｉ＋１＞Ｘｅ、ｙｉ＋１＞Ｙｅ、又はｚｉ＋１＞Ｚｅのいずれかを満たす場合（補間しない）
Ｖｃ（ｎ）＝Ｖｃ（ｘｉ、ｙｉ、ｚｉ、ｎ） （０≦ｎ≦２）
Ｖｃ（ｎ）＝０ （０≦ｎ≦２）
２）上記１）の条件を満たさない場合において、ｘｉ＋１＞Ｘｅ、ｙｉ＋１＞Ｙｅ、又はｚｉ＋１＞Ｚｅのいずれかを満たす場合（補間しない）
Ｖｃ（ｎ）＝Ｖｃ（ｘｉ、ｙｉ、ｚｉ、ｎ） （０≦ｎ≦２）
【0168】
３）上記１）２）の条件を満たさない場合（補間する）
Ｖｃ（ｎ）＝（１－ｗｚ）（１－ｗｙ）（１－ｗｘ）・Ｖｃ（ｘｉ、ｙｉ、ｚｉ、ｎ）＋（１－ｗｚ）（１－ｗｙ）・ｗｘ・Ｖｃ（ｘｉ＋１、ｙｉ、ｚｉ、ｎ）＋（１－ｗｚ）・ｗｙ・（１－ｗｘ）・Ｖｃ（ｘｉ、ｙｉ＋１、ｚｉ、ｎ）＋（１－ｗｚ）・ｗｙ・ｗｘ・Ｖｃ（ｘｉ＋１、ｙｉ＋１、ｚｉ、ｎ）＋ｗｚ・（１－ｗｙ）（１－ｗｘ）・Ｖｃ（ｘｉ、ｙｉ、ｚｉ＋１、ｎ）＋ｗｚ・（１－ｗｙ）・ｗｘ・Ｖｃ（ｘｉ＋１、ｙｉ、ｚｉ＋１、ｎ）＋ｗｚ・ｗｙ・（１－ｗｘ）・Ｖｃ（ｘｉ、ｙｉ＋１、ｚｉ＋１、ｎ）＋ｗｚ・ｗｙ・ｗｘ・Ｖｃ（ｘｉ＋１、ｙｉ＋１、ｚｉ＋１、ｎ） （０≦ｎ≦２）
Ｖｃ（ｎ）＝（１－ｗｚ）（１－ｗｙ）（１－ｗｘ）・Ｖｃ（ｘｉ、ｙｉ、ｚｉ、ｎ）＋（１－ｗｚ）（１－ｗｙ）・ｗｘ・Ｖｃ（ｘｉ＋１、ｙｉ、ｚｉ、ｎ）＋（１－ｗｚ）・ｗｙ・（１－ｗｘ）・Ｖｃ（ｘｉ、ｙｉ＋１、ｚｉ、ｎ）＋（１－ｗｚ）・ｗｙ・ｗｘ・Ｖｃ（ｘｉ＋１、ｙｉ＋１、ｚｉ、ｎ）＋ｗｚ・（１－ｗｙ）（１－ｗｘ）・Ｖｃ（ｘｉ、ｙｉ、ｚｉ＋１、ｎ）＋ｗｚ・（１－ｗｙ）・ｗｘ・Ｖｃ（ｘｉ＋１、ｙｉ、ｚｉ＋１、ｎ）＋ｗｚ・ｗｙ・（１－ｗｘ）・Ｖｃ（ｘｉ、ｙｉ＋１、ｚｉ＋１、ｎ）＋ｗｚ・ｗｙ・ｗｘ・Ｖｃ（ｘｉ＋１、ｙｉ＋１、ｚｉ＋１、ｎ） （０≦ｎ≦２）
【0169】
次に、図２０のフローチャートを参照して、図１６のステップＳ３３におけるレイキャスティング処理について説明する。
制御部１１は、まず、生成する２４ビット（ＲＧＢ）のレンダリング画像Ｉｍａｇｅ（ｘ、ｙ、ｎ）の初期値を全て０に設定する（Ｉｍａｇｅ（ｘ、ｙ、ｎ）＝０、ｎ＝０（Ｒ）、１（Ｇ）、２（Ｂ））。そして、サブサンプル回数Ｌとして、各２次元座標（ｘ、ｙ）（０≦ｘ≦Ｓｘ－１、０≦ｙ≦Ｓｙ－１）に対して、以下の処理を実行する。
制御部１１は、まず、生成する２４ビット（ＲＧＢ）のレンダリング画像Ｉｍａｇｅ（ｘ、ｙ、ｎ）の初期値を全て０に設定する（Ｉｍａｇｅ（ｘ、ｙ、ｎ）＝０、ｎ＝０（Ｒ）、１（Ｇ）、２（Ｂ））。そして、サブサンプル回数Ｌとして、各２次元座標（ｘ、ｙ）（０≦ｘ≦Ｓｘ－１、０≦ｙ≦Ｓｙ－１）に対して、以下の処理を実行する。
【0170】
まず、制御部１１は、光源Ｌｉｇｈｔ（ｎ）（ｎ＝０（Ｒ）、１（Ｇ）、２（Ｂ）、０≦Ｌｉｇｈｔ（ｎ）≦２５５）を設定し、サブサンプル・オフセット値をｄｘ＝ｄｙ＝ｄｚ＝０、処理回数ｌ＝０に初期化する（ステップＳ６１）。
【0171】
続いて、制御部１１は、ｚ＝Ｍ（ｘ、ｙ、ｌ）とし（ステップＳ６２）、ｚ＜０の場合、ステップＳ７０へ移行し、画素（ｘ、ｙ）におけるＲＧＢ画素値を決定する。
【0172】
一方、ｚ≧０の場合、制御部１１は、仮想光強度Ｔｒａｎｓ＝１．０、累積輝度値Ｅｎｅｒｇｙ（ｎ）＝０．０（０≦ｎ≦２）に初期化し（ステップＳ６３）、３次元座標（ｘ、ｙ、ｚ）を座標変換してボクセルα値（Ｖα）を取得し、Ａｌｐｈａ＝Ｖα／２５５を算出する（ステップＳ６４）。
【0173】
ステップＳ６４において算出したＡｌｐｈａがＡｌｐｈａ＝０かつｚ＝０の場合、ステップＳ７０へ移行し、画素（ｘ、ｙ）におけるＲＧＢ画素値を決定する。
【0174】
一方、ステップＳ６４において算出したＡｌｐｈａがＡｌｐｈａ＞０の場合、制御部１１は、３次元座標（ｘ、ｙ、ｚ）を座標変換してボクセルＲＧＢ値Ｖｃ（ｎ）を取得する（ステップＳ６７）。ボクセルのＲＧＢ値Ｖｃ（ｎ）を取得は、前述した座標変換処理により行われる。
【0175】
続いて、制御部１１は、累積輝度をＥｎｅｒｇｙ（ｎ）＝Ｅｎｅｒｇｙ（ｎ）＋Ｔｒａｎｓ／Ａｌｐｈａ・Ｖｃ（ｎ）／２５５、透過光強度をＴｒａｎｓ＝Ｔｒａｎｓ・（１．０－Ａｌｐｈａ）と更新する（ステップＳ６８）。Ａｌｐｈａ＝１．０又はＴｒａｎｓ＜０．００１の場合、ステップＳ７０へ移行し、画素（ｘ、ｙ）におけるＲＧＢ画素値を決定する。それ以外の場合、制御部１１は、ｚ＝ｚ－１に更新し（ステップＳ６９）、ｚ≧０の場合、ステップＳ６４に戻り、ｚ≦０の場合、ステップＳ７０へ移行し、画素（ｘ、ｙ）におけるＲＧＢ画素値を決定する。
【0176】
一方、ステップＳ６４において算出したＡｌｐｈａがＡｌｐｈａ＝０かつｚ＞０の場合、制御部１１は、Ｚｓ＝ｚ－１から座標変換して先頭の不透明ボクセルのＺ座標ｚ’を探索する（ステップＳ６５）。不透明ボクセルの探索は、前述した不透明ボクセル探索処理（図１９参照）により実行される。
【0177】
探索されたｚ’がｚ’＜０の場合、ステップＳ７０へ移行し、画素（ｘ、ｙ）におけるＲＧＢ画素値を決定する。ｚ’≧０の場合、制御部１１は、ｚ＝ｚ'に更新し（ステップＳ６６）、ステップＳ６４の処理に戻る。
【0178】
ステップＳ７０において、制御部１１は、画素（ｘ、ｙ）におけるＲＧＢ画素値を次のように決定する。
Ｉｍａｇｅ（ｘ、ｙ、ｎ）＝Ｉｍａｇｅ（ｘ、ｙ、ｎ）＋ｋ・Ｅｎｅｒｇｙ（ｎ）・Ｌｉｇｈｔ（ｎ）／Ｌ
ここで、ｋは強度倍率であり、初期値はｋ＝１．０に設定される。
Ｉｍａｇｅ（ｘ、ｙ、ｎ）＝Ｉｍａｇｅ（ｘ、ｙ、ｎ）＋ｋ・Ｅｎｅｒｇｙ（ｎ）・Ｌｉｇｈｔ（ｎ）／Ｌ
ここで、ｋは強度倍率であり、初期値はｋ＝１．０に設定される。
【0179】
続いて、制御部１１は、処理回数をｌ←ｌ＋１に更新し、サブサンプル・オフセット値を、ｄｘ←ｄｘ＋１／Ｌ、ｄｙ←ｄｙ＋１／Ｌ、ｄｚ←ｄｚ＋１／Ｌに更新する（ステップＳ７１）。制御部１１は、処理回数ｌ＝Ｌ－１となるまで（ステップＳ７２；Ｙｅｓ）、上記したステップＳ６２～Ｓ７１の処理を繰り返し、最終的なレンダリング画像Ｉｍａｇｅ（ｘ、ｙ、ｎ）（ｎ≦０≦２）を得る。
【0180】
【0181】
【0182】
ＧＰＵによるレンダリングはＯｐｅｎＧＬシステムを用いて実行するため、制御部１１は、作成したボクセル構造体Ｖを３ＤテクスチャとしてＯｐｅｎＧＬシステムへ登録する。具体的には、制御部１１は、２４ビットのボクセル構造体Ｖｃ（（式１２）参照）と８ビットのボクセル構造体Ｖα（式（１３）参照）を３２ビットのボクセル構造体Ｖ（（式１４）参照）に統合し、ｇｌＴｅｘＩｍａｇｅ３Ｄ関数を用いて３Ｄテクスチャを生成し、３ＤテクスチャをＯｐｅｎＧＬシステムに登録する。
【0183】
また、制御部１１は、図２１に示すように、３Ｄテクスチャを構成する各スライスを四角形で定義し、３次元空間のＸＹ座標面上のＺ軸方向に並べた積層四角形を設定し、テクスチャ座標と四角形座標との対応付けを行い、３Ｄテクスチャを各四角形に貼り付ける。３Ｄテクスチャマッピング法では、四角形に貼り付けた３Ｄテクスチャマップを視線方向に合わせて座標変換することで、任意方向のレンダリング像を生成する。
【0184】
尚、ＧＰＵ（第２レンダリング部２８ｂ）によるレンダリング処理では、レンダリング時にモアレが発生する場合がある。制御部１１は、レンダリング時のモアレの発生を抑えるため、次のように、関心領域ＲＯＩ境界面の透明ボクセルの色補正処理を事前に行う。
【0185】
図２２は色補正処理のフローチャートである。制御部１１は、境界面ボクセル（ｘ、ｙ、ｚ）を抽出する（ステップＳ８１）。境界面ボクセル（ｘ、ｙ、ｚ）とは、ｘ＝Ｘｓ、ｘ＝Ｘｅ、ｙ＝Ｙｓ、ｙ＝Ｙｅ、ｚ＝Ｚｓ、又はｚ＝Ｚｅのいずれかを満たすボクセルである。
【0186】
続いて、制御部１１は、抽出した境界面ボクセルの２７近傍ボクセルより、不透明ボクセルを抽出し、抽出した不透明ボクセルの平均色（Ｒａ、Ｇａ、Ｂａ）を算出する（ステップＳ８２）。不透明ボクセルとは、Ｖα（ｘ＋ｉ、ｙ＋ｊ、ｚ＋ｋ）＞０（－１≦ｉ≦１、－１≦ｊ≦１、－１≦ｋ≦１）の条件を満たすボクセルである。尚、不透明ボクセルが抽出されない場合は、制御部１１は、Ｒａ＝Ｇａ＝Ｂａ＝０とする。
【0187】
そして、制御部１１は、以下のように、境界面ボクセルの不透明度を０に設定し、ＲＧＢ値を算出した不透明ボクセルの平均色に設定する（ステップＳ８３）。
Ｖα（ｘ、ｙ、ｚ）＝０
Ｖｃ（ｘ、ｙ、ｚ、１）＝Ｒａ、Ｖｃ（ｘ、ｙ、ｚ、２）＝Ｇａ、Ｖｃ（ｘ、ｙ、ｚ、３）＝Ｂａ
Ｖα（ｘ、ｙ、ｚ）＝０
Ｖｃ（ｘ、ｙ、ｚ、１）＝Ｒａ、Ｖｃ（ｘ、ｙ、ｚ、２）＝Ｇａ、Ｖｃ（ｘ、ｙ、ｚ、３）＝Ｂａ
【0188】
制御部１１は、全ての境界面ボクセル（ｘ、ｙ、ｚ）について、上記した色補正処理（ステップＳ８１～８３）を実行する。
【0189】
次に、図２３のフローチャートを参照して、３Ｄテクスチャマッピング法によるレンダリング処理の流れを説明する。まず、制御部１１は、ＯｐｅｎＧＬの透視変換パラメータを設定する（ステップＳ９１）。仮想内視鏡モードでレンダリングを行う場合は、制御部１１は、ＯｐｅｎＧＬシステムに対して、次の透視変換パラメータを設定する。
【0190】
＜透視変換パラメータ＞
・カメラの視野角度（単位：度、平行投影の場合は０、焦点距離に相当）
・視点座標（カメラの位置、通常はＺ軸上に設定）
・注視点座標（被写体の凝視点、通常はＺ軸上で原点ｚ＝０に固定）
・近方クリッピング位置（視点からの距離、視点より若干距離を置く）
・遠方クリッピング位置（視点からの距離、通常は１０００など無限大に固定しクリッピングしない）
・カメラの視野角度（単位：度、平行投影の場合は０、焦点距離に相当）
・視点座標（カメラの位置、通常はＺ軸上に設定）
・注視点座標（被写体の凝視点、通常はＺ軸上で原点ｚ＝０に固定）
・近方クリッピング位置（視点からの距離、視点より若干距離を置く）
・遠方クリッピング位置（視点からの距離、通常は１０００など無限大に固定しクリッピングしない）
【0191】
一方、通常のレンダリングを行う際は、カメラの視野角度を０に設定し、平行投影モードに設定する。
【0192】
続いて、制御部１１は、ＯｐｅｎＧＬの投影画面設定を行う（ステップＳ９２）。
図２４は、投影画面設定の処理を示すフローチャートである。図２４に示すように、制御部１１は、レンダリング画像のスクリーンサイズ（縦横画素数、縦横アスペクト比率）を設定し（ステップＳ１０１）、平行投影（通常の外観レンダリング）又は透視投影（内視鏡モード）の設定を行う（ステップＳ１０２）。そして、透視投影が設定された場合は、制御部１１は、透視投影パラメータ（カメラの視野角度（焦点距離）、視点座標、注視点座標、近方クリッピング位置、遠方クリッピング位置）の設定を行う（ステップＳ１０３）。
図２４は、投影画面設定の処理を示すフローチャートである。図２４に示すように、制御部１１は、レンダリング画像のスクリーンサイズ（縦横画素数、縦横アスペクト比率）を設定し（ステップＳ１０１）、平行投影（通常の外観レンダリング）又は透視投影（内視鏡モード）の設定を行う（ステップＳ１０２）。そして、透視投影が設定された場合は、制御部１１は、透視投影パラメータ（カメラの視野角度（焦点距離）、視点座標、注視点座標、近方クリッピング位置、遠方クリッピング位置）の設定を行う（ステップＳ１０３）。
【0193】
図２３のフローチャートの説明に戻る。
続いて、制御部１１は、以下のように、ＯｐｅｎＧＬの座標変換パラメータ（回転、スケール、移動、Ｚ方向変倍パラメータＲｘｙ／Ｒｚ）を設定する（ステップＳ９３）。
続いて、制御部１１は、以下のように、ＯｐｅｎＧＬの座標変換パラメータ（回転、スケール、移動、Ｚ方向変倍パラメータＲｘｙ／Ｒｚ）を設定する（ステップＳ９３）。
【0194】
＜ＯｐｅｎＧＬの座標変換パラメータ＞
回転パラメータ行列Ｒ：
Ｒ＝［Ｒ１１ Ｒ１２ Ｒ１３ Ｒ１４；
Ｒ２１ Ｒ２２ Ｒ２３ Ｒ２４；
Ｒ３１ Ｒ３２ Ｒ３３ Ｒ３４；
Ｒ４１ Ｒ４２ Ｒ４３ Ｒ４４］
（ボクセル座標系から視点座標系への座標変換を行うための４×４の行列で、レンダリング側もＧＵＩ側と同一方向の座標変換をＯｐｅｎＧＬに指示する。実際は３×３の座標変換行列で定義できるが、ＯｐｅｎＧＬの仕様上、４×４の座標変換行列に拡張）
Ｘ軸方向のオフセットＸｏｆｆ
Ｙ軸方向のオフセットＹｏｆｆ
Ｚ軸方向のオフセットＺｏｆｆ
拡大縮小倍率Ｓｃａｌｅ（ＸＹＺ軸方向で同一）
Ｚ方向変倍率Ｓｃｚ＝Ｒｘｙ／Ｒｚ
Ｚ方向変倍後のＺ方向サイズＳｚ’＝Ｓｚ・Ｓｃｚ
回転パラメータ行列Ｒ：
Ｒ＝［Ｒ１１ Ｒ１２ Ｒ１３ Ｒ１４；
Ｒ２１ Ｒ２２ Ｒ２３ Ｒ２４；
Ｒ３１ Ｒ３２ Ｒ３３ Ｒ３４；
Ｒ４１ Ｒ４２ Ｒ４３ Ｒ４４］
（ボクセル座標系から視点座標系への座標変換を行うための４×４の行列で、レンダリング側もＧＵＩ側と同一方向の座標変換をＯｐｅｎＧＬに指示する。実際は３×３の座標変換行列で定義できるが、ＯｐｅｎＧＬの仕様上、４×４の座標変換行列に拡張）
Ｘ軸方向のオフセットＸｏｆｆ
Ｙ軸方向のオフセットＹｏｆｆ
Ｚ軸方向のオフセットＺｏｆｆ
拡大縮小倍率Ｓｃａｌｅ（ＸＹＺ軸方向で同一）
Ｚ方向変倍率Ｓｃｚ＝Ｒｘｙ／Ｒｚ
Ｚ方向変倍後のＺ方向サイズＳｚ’＝Ｓｚ・Ｓｃｚ
【0195】
制御部１１は、上記した回転パラメータ行列Ｒに対して、初期値として単位行列を設定する。すなわち、Ｒ１１＝１、Ｒ１２＝０、Ｒ１３＝０、Ｒ１４＝０、Ｒ２１＝０、Ｒ２２＝１、Ｒ２３＝０、Ｒ２４＝０、Ｒ３１＝０、Ｒ３２＝０、Ｒ３３＝１、Ｒ４１＝０、Ｒ４２＝０、Ｒ４３＝０、Ｒ４４＝１と設定する。
そして、ＧＵＩの指示に従い、Ｘ軸中心回転Ｒｘ、Ｙ軸中心回転Ｒｙ、Ｚ軸中心回転Ｒｚ（角度単位：ラジアン）のいずれかを逐次指定し、以下のように、各々回転行列Ａを生成して回転パラメータ行列Ｒに左から乗算して、回転パラメータ行列Ｒを更新する。
そして、ＧＵＩの指示に従い、Ｘ軸中心回転Ｒｘ、Ｙ軸中心回転Ｒｙ、Ｚ軸中心回転Ｒｚ（角度単位：ラジアン）のいずれかを逐次指定し、以下のように、各々回転行列Ａを生成して回転パラメータ行列Ｒに左から乗算して、回転パラメータ行列Ｒを更新する。
【0196】
回転行列Ａを
Ａ＝［Ａ１１ Ａ１２ Ａ１３ Ａ１４；
Ａ２１ Ａ２２ Ａ２３ Ａ２４；
Ａ３１ Ａ３２ Ａ３３ Ａ３４；
Ａ４１ Ａ４２ Ａ４３ Ａ４４］
とすると、
Ｘ軸中心回転Ｒｘの場合の回転行列Ａの各要素は、
Ａ１１＝１、Ａ１２＝０、Ａ１３＝０、Ａ１４＝０
Ａ２１＝０、Ａ２２＝ｃｏｓＲｘ、Ａ２３＝ｓｉｎＲｘ、Ａ２４＝０
Ａ３１＝０、Ａ３２＝ｓｉｎＲｘ、Ａ３３＝ｃｏｓＲｘ、Ａ３４＝０
Ａ４１＝０、Ａ４２＝０、Ａ４３＝０、Ａ４４＝１
Ｙ軸中心回転Ｒｙの場合の回転行列Ａの各要素は、
Ａ１１＝ｃｏｓＲｙ、Ａ１２＝０、Ａ１３＝ｓｉｎＲｙ、Ａ１４＝０
Ａ２１＝０、Ａ２２＝１、Ａ２３＝０、Ａ２４＝０
Ａ３１＝－ｓｉｎＲｙ、Ａ３２＝０、Ａ３３＝ｃｏｓＲｙ、Ａ３４＝０
Ａ４１＝０、Ａ４２＝０、Ａ４３＝０、Ａ４４＝１
Ｚ軸中心回転Ｒｚの場合の回転行列Ａの各要素は、
Ａ１１＝ｃｏｓＲｚ、Ａ１２＝ｓｉｎＲｚ、Ａ１３＝０、Ａ１４＝０
Ａ２１＝－ｓｉｎＲｚ、Ａ２２＝ｃｏｓＲｚ、Ａ２３＝０、Ａ２４＝０
Ａ３１＝０、Ａ３２＝０、Ａ３３＝１、Ａ３４＝０
Ａ４１＝０、Ａ４２＝０、Ａ４３＝０、Ａ４４＝１
となる。
回転パラメータ行列Ｒは、Ｒ←Ａ×Ｒと更新される。
Ａ＝［Ａ１１ Ａ１２ Ａ１３ Ａ１４；
Ａ２１ Ａ２２ Ａ２３ Ａ２４；
Ａ３１ Ａ３２ Ａ３３ Ａ３４；
Ａ４１ Ａ４２ Ａ４３ Ａ４４］
とすると、
Ｘ軸中心回転Ｒｘの場合の回転行列Ａの各要素は、
Ａ１１＝１、Ａ１２＝０、Ａ１３＝０、Ａ１４＝０
Ａ２１＝０、Ａ２２＝ｃｏｓＲｘ、Ａ２３＝ｓｉｎＲｘ、Ａ２４＝０
Ａ３１＝０、Ａ３２＝ｓｉｎＲｘ、Ａ３３＝ｃｏｓＲｘ、Ａ３４＝０
Ａ４１＝０、Ａ４２＝０、Ａ４３＝０、Ａ４４＝１
Ｙ軸中心回転Ｒｙの場合の回転行列Ａの各要素は、
Ａ１１＝ｃｏｓＲｙ、Ａ１２＝０、Ａ１３＝ｓｉｎＲｙ、Ａ１４＝０
Ａ２１＝０、Ａ２２＝１、Ａ２３＝０、Ａ２４＝０
Ａ３１＝－ｓｉｎＲｙ、Ａ３２＝０、Ａ３３＝ｃｏｓＲｙ、Ａ３４＝０
Ａ４１＝０、Ａ４２＝０、Ａ４３＝０、Ａ４４＝１
Ｚ軸中心回転Ｒｚの場合の回転行列Ａの各要素は、
Ａ１１＝ｃｏｓＲｚ、Ａ１２＝ｓｉｎＲｚ、Ａ１３＝０、Ａ１４＝０
Ａ２１＝－ｓｉｎＲｚ、Ａ２２＝ｃｏｓＲｚ、Ａ２３＝０、Ａ２４＝０
Ａ３１＝０、Ａ３２＝０、Ａ３３＝１、Ａ３４＝０
Ａ４１＝０、Ａ４２＝０、Ａ４３＝０、Ａ４４＝１
となる。
回転パラメータ行列Ｒは、Ｒ←Ａ×Ｒと更新される。
【0197】
図２３のフローチャートの説明に戻る。
制御部１１は、ＯｐｅｎＧＬによるレンダリング処理を実行する（ステップＳ９４）。
図２５は、ＯｐｅｎＧＬによるレンダリング処理を示すフローチャートである。図２５に示すように、制御部１１は、拡大縮小倍率Ｓｃａｌｅ（ＸＹＺ軸方向で同一）、及びＺ方向変倍率Ｓｃｚに基づいて３Ｄテクスチャマップに対してスケーリング、及びＺ方向変倍処理を施し（ステップＳ１１１）、回転パラメータ行列Ｒに基づいて３Ｄテスチャマップに対して回転処理を施し（ステップＳ１１２）、Ｘ軸方向のオフセットＸｏｆｆ、Ｙ軸方向のオフセットＹｏｆｆ、及びＺ軸方向のオフセットＺｏｆｆに基づいて３Ｄテスチャマップに対してオフセット処理を施すことで（ステップＳ１１３）、３Ｄテクスチャマップに対して所定の座標変換を行う。
制御部１１は、ＯｐｅｎＧＬによるレンダリング処理を実行する（ステップＳ９４）。
図２５は、ＯｐｅｎＧＬによるレンダリング処理を示すフローチャートである。図２５に示すように、制御部１１は、拡大縮小倍率Ｓｃａｌｅ（ＸＹＺ軸方向で同一）、及びＺ方向変倍率Ｓｃｚに基づいて３Ｄテクスチャマップに対してスケーリング、及びＺ方向変倍処理を施し（ステップＳ１１１）、回転パラメータ行列Ｒに基づいて３Ｄテスチャマップに対して回転処理を施し（ステップＳ１１２）、Ｘ軸方向のオフセットＸｏｆｆ、Ｙ軸方向のオフセットＹｏｆｆ、及びＺ軸方向のオフセットＺｏｆｆに基づいて３Ｄテスチャマップに対してオフセット処理を施すことで（ステップＳ１１３）、３Ｄテクスチャマップに対して所定の座標変換を行う。
【0198】
続いて、制御部１１は、３次元空間のＸＹ座標面上の四角形をＺ軸方向に並べた積層四角形をＯｐｅｎＧＬで設定し、３Ｄテスクチャマップの対応付けを行う（ステップＳ１１４、図２１参照）。そして、制御部１１は、ＯｐｅｎＧＬによるレンダリング処理を実行し、積層四角形に３Ｄテクスチャマップを貼り付けながらスキャンコンバージョンし、フレームメモリ上で奥からＺ方向（視点方向）にアルファブレンディング処理を行う（ステップＳ１１５）。すなわち、所定の視点からＺ軸方向に平行な視線上の四角形のＸＹ座標に対応する座標変換後の３Ｄテスクチャマップのボクセルの色値（ＲＧＢ値）を、ボクセルの不透明度（α値）に基づいて視点から遠い四角形の順にアルファブレンディングして取得し、ボリュームレンダリング像の画素値として与える。
【0199】
【0200】
制御部１１は、ＣＰＵ（第１レンダリング部２８ａ）により実行されるレイキャスティング法によるレンダリング処理（図１６）、又はＧＰＵ（第２レンダリング部２８ｂ）により実行される３Ｄテクスチャマッピング法によるレンダリング処理（図２３）によりボリュームレンダリング像を生成し、生成したボリュームレンダリング像を表示部１６に表示する。
【0201】
図６のフローチャートの説明に戻る。
制御部１１（データ出力部２９）は、ユーザからデータの保存操作を受け付けることによって、各種データを出力して記憶部１２に保存することができる（図６のステップＳ１０）。例えば、制御部１１は、ステップＳ９において生成し表示されたボリュームレンダリング像（３次元画像）を出力して保存する。
制御部１１（データ出力部２９）は、ユーザからデータの保存操作を受け付けることによって、各種データを出力して記憶部１２に保存することができる（図６のステップＳ１０）。例えば、制御部１１は、ステップＳ９において生成し表示されたボリュームレンダリング像（３次元画像）を出力して保存する。
【0202】
また、制御部１１は、ボクセル構造体Ｖ（ボクセルデータ）を作成する際に参照可能なデータ（参照データ）を出力しても良い。これにより、同じ断層画像群Ｄｏに基づいて再度レンダリングを行う際、いくつかの処理を省略することができる。参照データとしては、例えば、以下のものが挙げられる。
【0203】
＜参照データ１＞
ボクセル構造体Ｖ（Ｖｃ、Ｖα）の作成（（式１２）、（式１３）参照）に使用されたカラーマップＣｍａｐ（ｖ、ｎ）（０≦ｎ≦３）のデータ
ボクセル構造体Ｖ（Ｖｃ、Ｖα）の作成（（式１２）、（式１３）参照）に使用されたカラーマップＣｍａｐ（ｖ、ｎ）（０≦ｎ≦３）のデータ
【0204】
＜参照データ２＞
補正テーブルＳα（（式１１）参照）の作成に使用された楕円パラメータＰ（ｚ）（（式９）、（式１０）参照）のデータ
補正テーブルＳα（（式１１）参照）の作成に使用された楕円パラメータＰ（ｚ）（（式９）、（式１０）参照）のデータ
【0205】
＜参照データ３＞
ボクセル構造体Ｖαの作成（（式１３）参照）に使用された補正テーブルＳα（（式１１）参照）のデータ
ボクセル構造体Ｖαの作成（（式１３）参照）に使用された補正テーブルＳα（（式１１）参照）のデータ
【0206】
＜参照データ４＞
断層画像群Ｄｏの各画素（ｘ、ｙ、ｚ）に適用される補正後の不透明度を格納した不透明度テーブルＴα（ｘ、ｙ、ｚ）（３次元α値分布）
断層画像群Ｄｏの各画素（ｘ、ｙ、ｚ）に適用される補正後の不透明度を格納した不透明度テーブルＴα（ｘ、ｙ、ｚ）（３次元α値分布）
【0207】
ここで、不透明度テーブルＴαは、以下のように定義される。
【0208】
（式１５）
Ｔα（ｘ、ｙ、ｚ）＝Ｃｍａｐ（Ｄｏ（ｘ、ｙ、ｚ）、３）・Ｓα（ｘ、ｙ、ｚ）
（０≦Ｔα（ｘ、ｙ、ｚ）≦２５５、０≦ｘ≦Ｓｘ－１、０≦ｙ≦Ｓｙ－１、０≦ｚ≦Ｓｚ－１；解像度：Ｒｘｙ、Ｒｚ）
Ｔα（ｘ、ｙ、ｚ）＝Ｃｍａｐ（Ｄｏ（ｘ、ｙ、ｚ）、３）・Ｓα（ｘ、ｙ、ｚ）
（０≦Ｔα（ｘ、ｙ、ｚ）≦２５５、０≦ｘ≦Ｓｘ－１、０≦ｙ≦Ｓｙ－１、０≦ｚ≦Ｓｚ－１；解像度：Ｒｘｙ、Ｒｚ）
【0209】
＜参照データ５＞
ボクセル構造体Ｖ（Ｖｃ、Ｖα）の作成（（式１２）、（式１３）参照）に使用されたカラーマップＣｍａｐ（ｖ、ｎ）（０≦ｎ≦３）と、補正テーブルＳα（式１１参照）の作成に使用された楕円パラメータＰ（ｚ）（（式９）、（式１０）参照）を含むデータ
ボクセル構造体Ｖ（Ｖｃ、Ｖα）の作成（（式１２）、（式１３）参照）に使用されたカラーマップＣｍａｐ（ｖ、ｎ）（０≦ｎ≦３）と、補正テーブルＳα（式１１参照）の作成に使用された楕円パラメータＰ（ｚ）（（式９）、（式１０）参照）を含むデータ
【0210】
＜参照データ６＞
ボクセル構造体Ｖ（Ｖｃ、Ｖα）の作成（（式１２）、（式１３）参照）に使用されたカラーマップＣｍａｐ（ｖ、ｎ）（０≦ｎ≦３）と、ボクセル構造体Ｖαの作成（（式１３）参照）に使用された補正テーブルＳα（（式１１）参照）を含むデータ
ボクセル構造体Ｖ（Ｖｃ、Ｖα）の作成（（式１２）、（式１３）参照）に使用されたカラーマップＣｍａｐ（ｖ、ｎ）（０≦ｎ≦３）と、ボクセル構造体Ｖαの作成（（式１３）参照）に使用された補正テーブルＳα（（式１１）参照）を含むデータ
【0211】
＜参照データ７＞
ボクセル構造体Ｖｃの作成（（式１２）参照）に使用されたカラーマップＣｍａｐ（ｖ、ｎ）（０≦ｎ≦２）（カラーパレット）と、不透明度テーブルＴα（ｘ、ｙ、ｚ）を含むデータ
ボクセル構造体Ｖｃの作成（（式１２）参照）に使用されたカラーマップＣｍａｐ（ｖ、ｎ）（０≦ｎ≦２）（カラーパレット）と、不透明度テーブルＴα（ｘ、ｙ、ｚ）を含むデータ
【0212】
以上のような参照データを保存しておくことで、同じ断層画像群Ｄｏに基づいて再度レンダリングを行う際に、いくつかの処理を省略でき、レンダリング像を得るまでの処理時間を短縮できる。
【0213】
例えば、参照データ１（カラーマップＣｍａｐ）を保存しておけば、図６のステップＳ５～Ｓ７（カラーマップの取得・調整処理）を省略できる。
また、参照データ２（楕円パラメータＰ（ｚ））を保存しておけば、図８のステップＳ２１、Ｓ２２（補正テーブル作成処理の一部）を省略できる。
また、参照データ２（楕円パラメータＰ（ｚ））を保存しておけば、図８のステップＳ２１、Ｓ２２（補正テーブル作成処理の一部）を省略できる。
【0214】
また、参照データ３（補正テーブルＳα）を保存しておけば、図６のステップＳ４（補正テーブル作成処理）を省略できる。
また、参照データ４（不透明度テーブルＴα）を保存しておけば、図６のステップＳ４（補正テーブル作成処理）を省略できるとともに、（式１３）の不透明度の補正処理（Ｖα（ｘ、ｙ、ｚ）＝Ｖα（ｘ、ｙ、ｚ）・Ｓα（ｘ、ｙ、ｚ））を省略できる。すなわち、不透明度テーブルＴαに格納されている不透明度をボクセル構造体Ｖαに設定するだけでよい（Ｖα（ｘ、ｙ、ｚ）＝Ｔα（ｘ、ｙ、ｚ））。
また、参照データ４（不透明度テーブルＴα）を保存しておけば、図６のステップＳ４（補正テーブル作成処理）を省略できるとともに、（式１３）の不透明度の補正処理（Ｖα（ｘ、ｙ、ｚ）＝Ｖα（ｘ、ｙ、ｚ）・Ｓα（ｘ、ｙ、ｚ））を省略できる。すなわち、不透明度テーブルＴαに格納されている不透明度をボクセル構造体Ｖαに設定するだけでよい（Ｖα（ｘ、ｙ、ｚ）＝Ｔα（ｘ、ｙ、ｚ））。
【0215】
また、参照データ５（カラーマップＣｍａｐ＋楕円パラメータＰ（ｚ））を保存しておけば、図６のステップＳ５～Ｓ７（カラーマップの取得・調整処理）を省略できるとともに、図８のステップＳ２１、Ｓ２２（補正テーブル作成処理の一部）を省略できる。
【0216】
また、参照データ６（カラーマップＣｍａｐ＋補正テーブルＳα）を保存しておけば、図６のステップＳ４～Ｓ７（補正テーブルの作成処理、及びカラーマップの取得・調整処理）を省略できる。
【0217】
また、参照データ７（カラーマップＣｍａｐ（カラーパレット）＋不透明度テーブルＴα）を保存しておけば、図６のステップＳ４～Ｓ７（補正テーブルの作成処理、及びカラーマップの取得・調整処理）を省略できるとともに、（式１３）の不透明度の補正処理（Ｖα（ｘ、ｙ、ｚ）＝Ｖα（ｘ、ｙ、ｚ）・Ｓα（ｘ、ｙ、ｚ））を省略できる。
【0218】
【0219】
[実施例]
最後に、従来手法と提案手法により生成されるレンダリング像の比較を行う。従来手法とは、図６のステップＳ４において取得したカラーマップＣｍａｐを断層画像群Ｄｏにそのまま適用してボクセル構造体Ｖ（Ｖｃ、Ｖα）を作成し、レンダリング像を生成する手法である。すなわち、従来手法では、以下の（式１６）に基づいてボクセル構造体Ｖ（Ｖｃ、Ｖα）を作成する。
最後に、従来手法と提案手法により生成されるレンダリング像の比較を行う。従来手法とは、図６のステップＳ４において取得したカラーマップＣｍａｐを断層画像群Ｄｏにそのまま適用してボクセル構造体Ｖ（Ｖｃ、Ｖα）を作成し、レンダリング像を生成する手法である。すなわち、従来手法では、以下の（式１６）に基づいてボクセル構造体Ｖ（Ｖｃ、Ｖα）を作成する。
【0220】
（式１６）
（１）Ｘｓ＜ｘ＜Ｘｅ、Ｙｓ＜ｙ＜Ｙｅ、及びＺｓ＜ｚ＜Ｚｅを全て満たす場合（関心領域ＲＯＩ内の場合）
Ｖｃ（ｘ、ｙ、ｚ、ｎ）＝Ｃｍａｐ（Ｄｏ（ｘ、ｙ、ｚ）、ｎ）（０≦ｎ≦２）
Ｖα（ｘ、ｙ、ｚ）＝Ｃｍａｐ（Ｄｏ（ｘ、ｙ、ｚ）、３）
（２）ｘ≦Ｘｓ、ｘ≧Ｘｅ、ｙ≦Ｙｓ、ｙ≧Ｙｅ、ｚ≦Ｚｓ、又はｚ≧Ｚｅのいずれかを満たす場合（関心領域ＲＯＩ外の場合）
Ｖｃ（ｘ、ｙ、ｚ、ｎ）＝０ （０≦ｎ≦２）
Ｖα（ｘ、ｙ、ｚ）＝０
（１）Ｘｓ＜ｘ＜Ｘｅ、Ｙｓ＜ｙ＜Ｙｅ、及びＺｓ＜ｚ＜Ｚｅを全て満たす場合（関心領域ＲＯＩ内の場合）
Ｖｃ（ｘ、ｙ、ｚ、ｎ）＝Ｃｍａｐ（Ｄｏ（ｘ、ｙ、ｚ）、ｎ）（０≦ｎ≦２）
Ｖα（ｘ、ｙ、ｚ）＝Ｃｍａｐ（Ｄｏ（ｘ、ｙ、ｚ）、３）
（２）ｘ≦Ｘｓ、ｘ≧Ｘｅ、ｙ≦Ｙｓ、ｙ≧Ｙｅ、ｚ≦Ｚｓ、又はｚ≧Ｚｅのいずれかを満たす場合（関心領域ＲＯＩ外の場合）
Ｖｃ（ｘ、ｙ、ｚ、ｎ）＝０ （０≦ｎ≦２）
Ｖα（ｘ、ｙ、ｚ）＝０
【0221】
提案手法は、本実施の形態に基づいてレンダリング像を生成する方法（図６の処理によりレンダリング像を生成する方法）である。ここで、補正倍率Ｓα（（式１１）参照）を算出する際のパラメータを、ｋ＝１．１（減衰係数）、ａｍｐ＝１．０（振幅倍率）、ｘｏｆｆｓｅｔ＝２（Ｘ方向のオフセット）、ｙｏｆｆｓｅｔ＝５（Ｙ方向のオフセット）、ａ＝０．５８（横方向のサイズの補正係数）、ｂ＝０．５８（縦方向のサイズの補正係数）と設定した。
また、提案手法、従来手法の双方ともカラーマップＣｍａｐの調整（図６のステップＳ６の処理）は行わないものとし、双方とも第１レンダリング部２８ａを用いてレンダリングを実行した（第２レンダリング部２８ｂを用いても殆ど同様なレンダリング像が得られる）。
また、提案手法、従来手法の双方ともカラーマップＣｍａｐの調整（図６のステップＳ６の処理）は行わないものとし、双方とも第１レンダリング部２８ａを用いてレンダリングを実行した（第２レンダリング部２８ｂを用いても殆ど同様なレンダリング像が得られる）。
【0222】
図２６は、従来手法により生成したレンダリング像を正面から観察した画像であり、図２７は、提案手法により生成したレンダリング像を正面から観察した画像である。図２６（従来手法）では、肋骨及び胸骨の透明化が不十分であり、内部の臓器等を観察することが困難であるのに対し、図２７（提案手法）では、胸骨は多少残るものの肋骨全体が透明化され、内部の臓器等を良好に観察することができる。図２８、図２９は、それぞれ提案手法により生成したレンダリング像を背面、側面から観察した画像である。図２８、図２９に示すように、椎骨は残るものの、図２７と同様に肋骨が良好に透明化されていることが分かる。
【0223】
以上、本開示の第１の実施の形態について説明した。第１の実施の形態によれば、カラーマップＣｍａｐを参照して断層画像群Ｄｏの各画素（ｘ、ｙ、ｚ）の信号値ｖを色値及び不透明度に変換してボクセル構造体Ｖ（ボクセルデータ）を作成する際、カラーマップＣｍａｐ（ｖ、３）（オパシティカーブ）に定義された不透明度をそのまま用いるのではなく、被写体領域の外郭を近似した楕円の情報（楕円パラメータＰ（ｚ））に基づいて算出された補正倍率Ｓα（ｘ、ｙ、ｚ）を乗算することで不透明度を補正する（（式１３）参照）。そして、補正された不透明度を用いてボクセル構造体（ボクセルデータ）を作成し、レンダリング像を生成する。これにより、従来のように３次元マスクを作成することなく、特定の領域を鮮明に可視化させたレンダリング像を得ることができる。従来の３次元マスクの代わりに、本開示では不透明度の補正倍率Ｓα（ｘ、ｙ、ｚ）に関する補正テーブルの作成を必要とするが、これらは前述のように事前に補正テーブルを作成する方法をとらず、ボクセル作成時に各ボクセルごとに（式１１）で都度算出する方法をとることもでき、３次元マスクのようにメモリ上に確保する必要はない。
【0224】
また、補正倍率Ｓαは、被写体領域の外郭を近似した楕円の中心から遠ざかるにつれその倍率が小さくなるので、被写体領域の中心から外側にいくにつれ不透明度が減衰し、骨領域（胸骨や頭蓋骨）や皮膚、外部の寝台、固定用治具等が透明化され、被写体領域の中心付近にある臓器等が鮮明に可視化される。
【0225】
また、補正倍率Ｓαは、被写体領域の外郭の形状を少数の変数で表現可能な楕円パラメータＰ（ｚ）（中心座標、横方向のサイズ、縦方向のサイズ）に基づいて解析的に算出されるので、複雑な処理を伴わずに、不透明度の制御を行える。
【0226】
また、楕円の近似精度が良好でない断層画像Ｄｏ（ｚ）の楕円パラメータＰ（ｚ）を、近似精度が良好な他の断層画像Ｄｏ（ｚ）の楕円パラメータＰ（ｚ）に基づいて補正することで、撮影条件等によって楕円近似が良好に行えない場合でも、楕円パラメータＰ（ｚ）を良好に算出することができる。
【0227】
また、カラーマップＣｍａｐ（ｖ、３）（オパシティカーブ）を信号値に応じて調整することで、骨領域をより効果的に減衰させることができる（図７参照）。
【0228】
また、従来のマスク処理のように、各ボクセルの可視／不可視の制御（オン／オフの制御）では、単一のボクセルに可視の内臓領域と不可視の骨領域が混在する境界領域を表現することができず、自然なレンダリング像が得られない場合があった。これに対して、本実施の形態では、マスク処理のようなオン／オフの制御でなく、被写体領域の外部に向かって不透明度を徐々に減衰させ、オンとオフの中間状態をもたせながら滑らかな制御を行うことで、自然なレンダリング像を得ることができる。
【0229】
［第２の実施の形態］
頭部は球形状に近いため、胸部に比べ、被写体領域に占める骨領域の面積比がスライス方向（Ｚ軸方向）に大きく変化する。特にスライス末端近くでは被写体領域に占める骨領域（頭頂骨や顎骨）の面積比が増大するため、骨領域（頭頂骨や顎骨）を透過させて脳の頭頂領域や脳幹領域を精度よく可視化することが困難となる場合がある。例えば、頭頂部では断層像の中心領域も骨領域であるため、同心楕円のパターンでは骨領域を透明にすることはできない。そこで、第２の実施の形態では、断層画像のスライス位置を更に考慮して補正倍率Ｓαを算出することで、特に頭部ＣＴにおけるレンダリングの精度向上を図る。
頭部は球形状に近いため、胸部に比べ、被写体領域に占める骨領域の面積比がスライス方向（Ｚ軸方向）に大きく変化する。特にスライス末端近くでは被写体領域に占める骨領域（頭頂骨や顎骨）の面積比が増大するため、骨領域（頭頂骨や顎骨）を透過させて脳の頭頂領域や脳幹領域を精度よく可視化することが困難となる場合がある。例えば、頭頂部では断層像の中心領域も骨領域であるため、同心楕円のパターンでは骨領域を透明にすることはできない。そこで、第２の実施の形態では、断層画像のスライス位置を更に考慮して補正倍率Ｓαを算出することで、特に頭部ＣＴにおけるレンダリングの精度向上を図る。
【0230】
第２の実施の形態では、第１の実施形態における図８のステップＳ２３の処理（補正テーブルＳαを作成する処理）が次のように変わる。すなわち、図８のステップＳ２３において、制御部１１（補正倍率算出部２６ｃ）は、断層画像Ｄｏ（ｚ）毎に、断層画像Ｄｏ（ｚ）に対応する幾何パラメータＰ（ｚ）（楕円パラメータＰ（ｚ））及び断層画像Ｄｏ（ｚ）のスライス位置の情報（以下、「スライスパラメータＰｓ（ｚ）」と表記）に基づいて、断層画像Ｄｏ（ｚ）の各画素（ｘ、ｙ）に対応する不透明度の補正倍率Ｓα（ｘ、ｙ、ｚ）を算出し、算出した補正倍率Ｓα（ｘ、ｙ、ｚ）が格納された３次元の補正テーブルＳαを作成する。
例えば、ｚ番目の断層画像Ｄｏ（ｚ）の各画素（ｘ、ｙ）の補正倍率Ｓα（ｘ、ｙ、ｚ）は、以下のように算出される。
例えば、ｚ番目の断層画像Ｄｏ（ｚ）の各画素（ｘ、ｙ）の補正倍率Ｓα（ｘ、ｙ、ｚ）は、以下のように算出される。
【0231】
まず、図８のステップＳ２３の（式１１）と同様に、制御部１１は、断層画像Ｄｏ（ｚ）に対応する幾何パラメータＰ（ｚ）（楕円パラメータＰ（ｚ））に基づいて、補正倍率Ｓα（ｘ、ｙ、ｚ）を、次のように算出する。
【0232】
（式１７）
Ｓα（ｘ、ｙ、ｚ）＝１－ｋ・ｒ（ｘ、ｙ、ｚ）・ａｍｐ
但し、０≦Ｓα（ｘ、ｙ、ｚ）≦１
ｒ（ｘ、ｙ、ｚ）＝｛（ｘ－Ｃｘ（ｚ）－ｘｏｆｆｓｅｔ）／（ａ・Ｗ（ｚ）／２）｝２＋｛（ｙ－Ｃｙ（ｚ）－ｙｏｆｆｓｅｔ）／（ｂ・Ｈ（ｚ）／２）｝２ｓ
Ｓα（ｘ、ｙ、ｚ）＝１－ｋ・ｒ（ｘ、ｙ、ｚ）・ａｍｐ
但し、０≦Ｓα（ｘ、ｙ、ｚ）≦１
ｒ（ｘ、ｙ、ｚ）＝｛（ｘ－Ｃｘ（ｚ）－ｘｏｆｆｓｅｔ）／（ａ・Ｗ（ｚ）／２）｝２＋｛（ｙ－Ｃｙ（ｚ）－ｙｏｆｆｓｅｔ）／（ｂ・Ｈ（ｚ）／２）｝２ｓ
【0233】
更に第２の実施の形態では、制御部１１は、次のように、スライスパラメータＰｓ（ｚ）に基づいて、Ｚ方向の補正倍率Ｓαｚ（ｚ）を算出し、Ｓα（ｘ、ｙ、ｚ）に乗算する。ここで、スライスパラメータＰｓ（ｚ）とは、スライス総数Ｓｚ及びスライス順位ｚのことを言う。
【0234】
（式１８）
Ｓαｚ（ｚ）＝１．０－ｋｚ・｛（ｚ－Ｓｚ／２－ｚｏｆｆｓｅｔ）／Ｓｚ・ｃ／２｝｝２・ａｍｐｚ
Ｓα（ｘ、ｙ、ｚ）＝Ｓα（ｘ、ｙ、ｚ）・Ｓαｚ（ｚ）
但し、０≦Ｓα（ｘ、ｙ、ｚ）、Ｓαｚ（ｚ）≦１
Ｓαｚ（ｚ）＝１．０－ｋｚ・｛（ｚ－Ｓｚ／２－ｚｏｆｆｓｅｔ）／Ｓｚ・ｃ／２｝｝２・ａｍｐｚ
Ｓα（ｘ、ｙ、ｚ）＝Ｓα（ｘ、ｙ、ｚ）・Ｓαｚ（ｚ）
但し、０≦Ｓα（ｘ、ｙ、ｚ）、Ｓαｚ（ｚ）≦１
【0235】
（式１７）（式１８）において、ｋはＸＹ平面、ｋｚはＺ方向の減衰係数である。初期値はｋ＝１．０、ｋｚ＝０．０であり、この場合、骨及び寝台が不可視となる。頭部ＣＴの場合はｋｚ＞０．０に設定する。
ａ、ｂはそれぞれ横方向のサイズＷ（ｚ）、縦方向のサイズＨ（ｚ）の補正係数であり、ｃはＺ方向幅Ｓｚ／２の補正係数（スライス総数倍率）である。補正しない場合はａ＝ｂ＝ｃ＝１．０に設定する。
ｘｏｆｆｓｅｔ、ｙｏｆｆｓｅｔはＸＹ面、ｚｏｆｆｓｅｔはＺ方向のオフセット値であり、初期値はｘｏｆｆｓｅｔ＝ｙｏｆｆｓｅｔ＝ｚｏｆｆｓｅｔ＝０（単位：ボクセル）である。
ａ、ｂはそれぞれ横方向のサイズＷ（ｚ）、縦方向のサイズＨ（ｚ）の補正係数であり、ｃはＺ方向幅Ｓｚ／２の補正係数（スライス総数倍率）である。補正しない場合はａ＝ｂ＝ｃ＝１．０に設定する。
ｘｏｆｆｓｅｔ、ｙｏｆｆｓｅｔはＸＹ面、ｚｏｆｆｓｅｔはＺ方向のオフセット値であり、初期値はｘｏｆｆｓｅｔ＝ｙｏｆｆｓｅｔ＝ｚｏｆｆｓｅｔ＝０（単位：ボクセル）である。
【0236】
（式１８）によって、断層画像Ｄｏ（ｚ）のスライス位置が中央から末端に位置するにつれ（スライス総数Ｓｚ／２とスライス順位ｚとの差が大きくなるにつれ）、補正倍率Ｓα（ｘ、ｙ、ｚ）が一律に減衰し、楕円の中心であってもＳα（ｘ、ｙ、ｚ）は１未満の値をとる。これによって、スライス末端近くの骨領域（頭頂骨や顎骨）を効果的に透過させることができる。
【0237】
図３０は、補正倍率Ｓα（ｘ、ｙ、ｚ）の分布パターンを示す図である。
図３０の上図は、（式１７）で算出される楕円パラメータＰ（ｚ）のみに基づいた補正倍率Ｓα（ｘ、ｙ、ｚ）の分布パターン（同心楕円分布）を示し、図３０の下図は、（式１８）で算出される楕円パラメータＰ（ｚ）及びスライスパラメータＰｓ（ｚ）に基づいた補正倍率Ｓα（ｘ、ｙ、ｚ）の分布パターン（同心楕円分布）を示す。図３０に示すように、Ｚ方向の補正倍率Ｓαｚ（ｚ）が乗算されることによって、補正倍率Ｓα（ｘ、ｙ、ｚ）がスライス位置に応じて一律に減衰する。すなわち、補正倍率Ｓαの分布パターンが、第１の実施の形態では２次元の同心楕円分布であったのに対し、第２の実施の形態では３次元の同心楕円体分布に拡張される。
図３０の上図は、（式１７）で算出される楕円パラメータＰ（ｚ）のみに基づいた補正倍率Ｓα（ｘ、ｙ、ｚ）の分布パターン（同心楕円分布）を示し、図３０の下図は、（式１８）で算出される楕円パラメータＰ（ｚ）及びスライスパラメータＰｓ（ｚ）に基づいた補正倍率Ｓα（ｘ、ｙ、ｚ）の分布パターン（同心楕円分布）を示す。図３０に示すように、Ｚ方向の補正倍率Ｓαｚ（ｚ）が乗算されることによって、補正倍率Ｓα（ｘ、ｙ、ｚ）がスライス位置に応じて一律に減衰する。すなわち、補正倍率Ｓαの分布パターンが、第１の実施の形態では２次元の同心楕円分布であったのに対し、第２の実施の形態では３次元の同心楕円体分布に拡張される。
【0238】
尚、（式１７）のａｍｐはＸＹ面、（式１８）のａｍｐｚはＺ方向の振幅倍率であり、初期値はａｍｐ＝ａｍｐｚ＝１．０である。図３１に示すように、被写体中心部の補正倍率を平坦にし、飽和領域を広くしたい場合は、ａｍｐ、ａｍｐｚ＞１．０に設定する（但し、陰影感が弱くなる）。
【0239】
制御部１１は、全ての断層画像Ｄｏ（ｚ）について、（式１７）（式１８）によりＺ方向の補正倍率Ｓαｚ（ｚ）を乗算した補正倍率Ｓα（ｘ、ｙ、ｚ）（０≦ｘ≦Ｓｘ－１、０≦ｙ≦Ｓｙ－１、０≦ｚ≦Ｓｚ－１）を算出し、算出された補正倍率Ｓα（ｘ、ｙ、ｚ）を格納する３次元の補正テーブルＳαを作成する。
【0240】
以上のように第２の実施の形態によれば、断層画像のスライス位置を考慮して補正倍率Ｓαを算出するので、更に柔軟な不透明度の制御を行うことができる。特に、頭部ＣＴ画像をレンダリングする際に有効である。頭部は球形状に近いため、胸部に比べ、被写体領域に占める骨領域の面積比がスライス方向（Ｚ軸方向）に大きく変化し、特に骨領域の面積比が大きくなるスライス末端近くの頭頂骨や顎骨の領域を透過させて脳の頭頂領域や脳幹領域を精度よく可視化することが困難となる場合がある。第２の実施の形態によれば、断層画像のスライス位置が中央から末端に位置するにつれ、当該断層画像の全画素に対応する補正倍率Ｓαを一律に小さくするので、スライス末端近くで被写体領域に占める骨領域（頭頂骨や顎骨）が増大する頭部においても、骨領域を効果的に透過させることができる。
【0241】
［第３の実施の形態］
第３の実施の形態は、ＧＵＩ上で、補正テーブルＳαを作成する際の各種パラメータをユーザが調整できるようにしたものである。
第３の実施の形態は、ＧＵＩ上で、補正テーブルＳαを作成する際の各種パラメータをユーザが調整できるようにしたものである。
【0242】
図３２は、第３の実施の形態における医用画像処理装置１ａの機能構成を示す図である。図３２に示すように、第３の実施形態では、第１の実施の形態における医用画像処理装置１の機能構成（図３参照）に加え、パラメータ調整部３０を更に備える。
【0243】
パラメータ調整部３０は、パラメータ調整画面４０を表示し、補正倍率Ｓαを算出する際の各種パラメータをユーザに調整させる。具体的には、補正倍率Ｓαが（式１１）に示したように楕円パラメータＰ（ｚ）に基づいて算出される場合（第１の実施の形態の場合）には、楕円の横方向のサイズの補正係数ａ、縦方向のサイズの補正係数ｂ、Ｘ方向のオフトセットｘｏｆｆｓｅｔ、Ｙ方向のオフセットｙｏｆｆｓｅｔ、ＸＹ平面の減衰係数ｋを調整させる。
【0244】
一方、補正倍率Ｓαが（式１７）（式１８）に示したように楕円パラメータＰ（ｚ）及びスライスパラメータＰｓ（ｚ）に基づいて算出される場合（第２の実施の形態の場合）には、上記パラメータに加え、スライス総数倍率ｃ、及びＺ方向のオフセットｚｏｆｆｓｅｔ、Ｚ方向の減衰係数ｋｚを更に調整可能とする。
【0245】
図３３は、第３の実施の形態における医用画像処理装置１ａの全体動作を示すフローチャートである。ボリュームレンダリング像を生成するまでの処理（ステップＳ２０１～ステップＳ２０９）は、図６の第１の実施の形態のステップＳ１～ステップＳ９と同様である。第３の実施の形態では、ステップＳ２０９においてボリュームレンダリング像を生成し表示部１６に表示した後、ユーザがパラメータの調整が必要と判断した場合（ステップＳ２１０；「Ｙｅｓ」）、パラメータ調整画面４０においてパラメータをユーザに調整させる（ステップＳ２１１）。
【0246】
図３４は、パラメータ調整画面４０の例を示す。図に示すように、パラメータ調整画面４０において、ＸＹ平面の減衰係数ｋ（図のＸＹ方向減衰［％］）、Ｚ方向の減衰係数ｋｚ（図のＺ方向減衰［％］）、楕円の横方向のサイズの補正係数ａ（図の楕円倍率Ｘ）、楕円の縦方向のサイズの補正係数ｂ（図の楕円倍率Ｙ）、スライス総数倍率ｃ（図の楕円倍率Ｚ）、Ｘ方向のオフトセットｘｏｆｆｓｅｔ（図の中心オフセットＸ）、Ｙ方向のオフセットｙｏｆｆｓｅｔ（図の中心オフセットＹ）、Ｚ方向のオフセットｚｏｆｆｓｅｔ（図の中心オフセットＺ）を調整することができる。
【0247】
パラメータの調整が終了すると（パラメータ調整画面４０の「ＯＫ」ボタンを押下すると）、ステップＳ２０７に戻り、制御部１１は、ステップＳ２０７～ステップＳ２０９の処理を再度実行する。すなわち、制御部１１は、調整されたパラメータに基づいて、再度、補正テーブルＳαを作成し直した上で（ステップＳ２０７）、ボクセル構造体Ｖ（ボクセルデータ）を作成し（ステップＳ２０８）、レンダリング処理を行う（ステップＳ２０９）。ステップＳ２０７～ステップＳ２０９の処理は、ステップＳ２０９において生成し表示されるレンダリング像を確認しながら、ユーザがパラメータの調整が十分と判断するまで（ステップＳ２１０；「Ｎｏ」）、繰り返し実行する。良好なレンダリング結果が得られ、ユーザがパラメータの調整が十分と判断すると（ステップＳ２１０；「Ｎｏ」）、制御部１１は、必要に応じてレンダリング結果等のデータを出力し（ステップＳ２１２）、処理を終了する。ステップＳ２１２のデータの出力処理は、図６の第１の実施の形態のステップＳ１０と同様である。
【0248】
以上、第３の実施の形態によれば、補正テーブルＳαを作成する際の各種パラメータをＧＵＩ上でユーザが調整できる。特に、被写体領域に撮影治具（寝台、固定具等）が映りこむと、被写体領域の外郭を近似する楕円の径が実際の被写体幅（頭幅や胸幅）より大きめに算出されたり、楕円の中心位置がズレたりし、楕円の近似精度が悪化する場合がある。第３の実施の形態によれば、ＧＵＩ上で楕円の径や中心位置等を適宜調整することができるので、上記のような場合でも、良好なレンダリング結果を得ることができる。
【0249】
尚、本実施の形態では、レンダリング結果を見ながらパラメータの調整が行えるよう、レンダリング後にパラメータを調整可能とする例を説明したが、パラメータはレンダリング後だけでなく、補正テーブルＳαを作成する前の任意の段階で調整可能である。また、ステップＳ２０７における調整されたパラメータに基づく補正テーブルＳαの再作成を行わずに、ステップＳ２０８におけるボクセル構造体Ｖ（ボクセルデータ）の再作成の段階で、（式１７）（式１８）に基づいて補正倍率をボクセルごとに都度算出する方法をとることもできる。
【0250】
［第４の実施の形態］
第４の実施の形態は、更にＭＰＲ（Multi Plane
Reconstruction）像を生成し表示することを特徴とする。また、第３の実施形態におけるパラメータ調整を、ＭＰＲ像を確認しながら行うことができる。
第４の実施の形態は、更にＭＰＲ（Multi Plane
Reconstruction）像を生成し表示することを特徴とする。また、第３の実施形態におけるパラメータ調整を、ＭＰＲ像を確認しながら行うことができる。
【0251】
図３５は、第４の実施の形態における医用画像処理装置１ｂの機能構成を示す図である。図３５に示すように、第４の実施の形態では、第３の実施の形態における医用画像処理装置１ａの機能構成（図３２参照）に加え、ＭＰＲ生成部３１、ＭＰＲ表示部３２を更に備える。
【0252】
ＭＰＲ生成部３１は、断層画像群Ｄｏに基づいて、任意断面の画像を生成し表示する。例えば、図３６に示すように、ＸＹ平面に平行な体軸断面（アキシャル（Ａｘｉａｌ）断面）、ＸＺ平面に平行な冠状断面（コロナル（Ｃｏｒｏｎａｌ）断面）、ＺＹ平面に平行な矢状断面（サジタル（Ｓａｇｉｔｔａｌ）断面）のＭＰＲ像を生成する。また、ＸＹＺ空間を斜めに切断した斜断面（オブリーク（Ｏｂｌｉｑｕｅ）断面）のＭＰＲ像を生成しても良い。
【0253】
特に、ＭＰＲ生成部３１は、断層画像群Ｄｏの各画素（ｘ、ｙ、ｚ）の信号値ｖに対して、各画素（ｘ、ｙ、ｚ）に対応する補正倍率Ｓα（ｘ、ｙ、ｚ）を乗算することで信号値ｖを補正し、補正された信号値ｖに基づいてＭＰＲ像を生成する。これにより、被写体形状（幾何形状）に応じて信号値ｖを減衰させたＭＰＲ像が得られる。
【0254】
ＭＰＲ表示部３２は、ＭＰＲ生成部３１により生成されたＭＰＲ像を表示部１６に表示する。また、信号値補正前と補正後の結果を比較するために、補正前の断層画像群Ｄｏ（原画像）に基づいて生成したＭＰＲ像を併せて表示しても良い。
【0255】
図３７は、第４の実施の形態における医用画像処理装置１ｂの全体動作を示すフローチャートである。補正テーブルＳαを作成するまでの処理（ステップＳ３０１～ステップＳ３０７）は、図３３の第３の実施の形態のステップＳ２０１～ステップＳ２０７（図６の第１の実施の形態のステップＳ１～ステップＳ７）と同様である。
【0256】
第４の実施の形態では、制御部１１は、ステップＳ３０７において補正テーブルＳαを作成した後、ＭＰＲ像を生成し表示部１６に表示する（ステップＳ３０８）。具体的には、以下のように、断層画像群Ｄｏの信号値ｖ（＝Ｄｏ（ｘ、ｙ、ｚ））を補正したうえで、ＭＰＲ像を生成し表示する。
【0257】
断層画像群Ｄｏの信号値が１６ビットの場合（（式１）参照）には、制御部１１は、次の手順で断層画像群Ｄｏの信号値ｖ（＝Ｄｏ（ｘ、ｙ、ｚ））を補正する。
【0258】
（式１９）
（１）断層画像群Ｄｏの中で、Ｓｚ／２番目の中間スライスＤｏ（ｚ／２）における全ての画素の最小値Ｄｍｉｎ、最大値Ｄｍａｘを算出する。
（２）下限値Ｌｍｉｎ＝（Ｄｍａｘ－Ｄｍｉｎ）・γ＋Ｄｍｉｎを設定する。ここで、γは階調圧縮画像のコントラスト調整幅で、０に近いほどコントラストが増大するが輝度が小さくなる。通常はγ＝０．１に設定する。
（３）次のように、断層画像群Ｄｏの信号値ｖ（＝Ｄｏ（ｘ、ｙ、ｚ））を補正する。
Ｄｏ（ｘ、ｙ、ｚ）≧Ｌｍｉｎの場合：
Ｄｏ（ｘ、ｙ、ｚ）＝（Ｄｏ（ｘ、ｙ、ｚ）－Ｌｍｉｎ）・Ｓα（ｘ、ｙ、ｚ）＋Ｌｍｉｎ
Ｄｏ（ｘ、ｙ、ｚ）≦Ｌｍｉｎの場合：
Ｄｏ（ｘ、ｙ、ｚ）＝Ｌｍｉｎ
（１）断層画像群Ｄｏの中で、Ｓｚ／２番目の中間スライスＤｏ（ｚ／２）における全ての画素の最小値Ｄｍｉｎ、最大値Ｄｍａｘを算出する。
（２）下限値Ｌｍｉｎ＝（Ｄｍａｘ－Ｄｍｉｎ）・γ＋Ｄｍｉｎを設定する。ここで、γは階調圧縮画像のコントラスト調整幅で、０に近いほどコントラストが増大するが輝度が小さくなる。通常はγ＝０．１に設定する。
（３）次のように、断層画像群Ｄｏの信号値ｖ（＝Ｄｏ（ｘ、ｙ、ｚ））を補正する。
Ｄｏ（ｘ、ｙ、ｚ）≧Ｌｍｉｎの場合：
Ｄｏ（ｘ、ｙ、ｚ）＝（Ｄｏ（ｘ、ｙ、ｚ）－Ｌｍｉｎ）・Ｓα（ｘ、ｙ、ｚ）＋Ｌｍｉｎ
Ｄｏ（ｘ、ｙ、ｚ）≦Ｌｍｉｎの場合：
Ｄｏ（ｘ、ｙ、ｚ）＝Ｌｍｉｎ
【0259】
一方、断層画像群Ｄｏの信号値が８ビットの場合（（式２）参照）には、制御部１１は、次のように信号値ｖ（Ｄｏ（ｘ、ｙ、ｚ））を補正する。
【0260】
（式２０）
Ｄｏ（ｘ、ｙ、ｚ）＝Ｄｏ（ｘ、ｙ、ｚ）・Ｓα（ｘ、ｙ、ｚ）
Ｄｏ（ｘ、ｙ、ｚ）＝Ｄｏ（ｘ、ｙ、ｚ）・Ｓα（ｘ、ｙ、ｚ）
【0261】
図３８、図３９は、頭部ＣＴ画像に基づいて生成されたＭＰＲ像の表示例を示す図である。図３８、図３９は、補正前の断層画像群Ｄｏ（原画像）に基づいて生成されたＭＰＲ像と、補正後の断層画像Ｄｏに基づいて生成されたＭＰＲ像と、を体軸断面、冠状断面、矢状断面の各断面について比較表示した例である。図に示すように、補正倍率Ｓαにより信号値ｖを補正することによって、頭蓋領域の信号値が減衰していることが分かる。このように、補正倍率Ｓαの適用効果をＭＰＲ像の表示により確認することができる。
【0262】
制御部１１は、ＭＰＲ像を生成し表示した後、ユーザがパラメータの調整が必要と判断した場合（ステップＳ３０９；「Ｙｅｓ」）、パラメータ調整画面４０（図３４参照）においてユーザにパラメータを調整させる（ステップＳ３１０）。
【0263】
パラメータの調整が終了すると、ステップＳ３０７に戻り、制御部１１は、ステップＳ３０７～ステップＳ３０８の処理を再度実行する。すなわち、制御部１１は、調整されたパラメータに基づいて、再度、補正テーブルＳαを作成し直した上で（ステップＳ３０８）、ＭＰＲ像を生成し表示部１６に表示する（ステップＳ３０８）。ステップＳ３０７～ステップＳ３０８の処理は、ステップＳ３０８において生成し表示されるＭＰＲ像を観察しながら、ユーザがパラメータの調整が十分と判断するまで（ステップＳ３０９；「Ｎｏ」）、繰り返し実行する。
【0264】
良好なＭＰＲ像が得られ、ユーザがパラメータの調整が十分と判断すると（ステップＳ３０９；「Ｎｏ」）、制御部１１は、ボクセル構造体Ｖ（ボクセルデータ）を作成し（ステップＳ３１１）、レンダリング処理を行い（ステップＳ３１２）、必要であればレンダリング結果等のデータを出力し（ステップＳ３１２）、処理を終了する。ステップＳ３１１～ステップＳ３１３の処理は、図６の第１の実施の形態のステップＳ８～ステップＳ１０と同様である。
【0265】
以上、第４の実施の形態によれば、断層画像群Ｄｏに基づいてＭＰＲ像を生成し表示する。特に、断層画像群Ｄｏの各画素の信号値ｖに対して、補正倍率Ｓαを乗算して信号値ｖを補正し、補正された信号値ｖに基づいてＭＰＲ像を生成する。これにより、３次元マスクを作成することなく、特定の領域を鮮明に可視化させたＭＰＲ像を得ることができる。従来の３次元マスクの代わりに、本開示では不透明度の補正倍率Ｓα（ｘ、ｙ、ｚ）に関する補正テーブルの作成を必要とするが、これらは前述のように事前にステップ３０７で補正テーブルを作成する方法をとらず、ステップ３０８においてＭＰＲ像を生成する際に各画素ごとに（式１７）（式１８）で都度算出する方法をとることもでき、３次元マスクのようにメモリ上に確保する必要はない。
【0266】
また、補正テーブルＳαの適用効果を、ボリュームレンダリング像を生成する前に、ＭＰＲ像にて確認することができ、また、ＭＰＲ像の結果を見ながら補正テーブルＳαのパラメータ調整を行うことができる。これにより、計算量が多いレンダリング処理の実行前に、補正テーブルＳαのパラメータ調整を行うことができるので、所望のレンダリング像を得るまでの全体時間を短縮することができる。
【0267】
尚、本実施の形態では、ＭＰＲ像の結果を見ながらパラメータの調整が行えるよう、ＭＰＲ表示後にパラメータを調整可能とする例を説明したが、パラメータはＭＰＲ表示後だけでなく、補正テーブルＳαを作成する前の任意の段階で調整可能である。また、第３の実施の形態のように、レンダリング後にパラメータの調整が行えるようにしても良い。
【0268】
以上、添付図面を参照しながら、本開示に係る医用画像処理装置等の好適な実施形態について説明したが、本開示はかかる例に限定されない。当業者であれば、本願で開示した技術的思想の範疇内において、各種の変更例又は修正例に想到し得ることは明らかであり、それらについても当然に本開示の技術的範囲に属するものと了解される。
【符号の説明】
【0269】
１、１ａ、１ｂ：医用画像処理装置
２１：画像取得部
２２：階調圧縮部
２３：領域指定部
２４：カラーマップ取得部
２５：カラーマップ調整部
２６：補正テーブル作成部
２６ａ：幾何情報取得部
２６ｂ：幾何情報補正部
２６ｃ：補正倍率算出部
２７：ボクセル作成部
２８：レンダリング部
２８ａ：第１レンダリング部
２８ｂ：第２レンダリング部
２９：データ出力部
３０：パラメータ調整部
３１：ＭＰＲ生成部
３２：ＭＰＲ表示部
４０：パラメータ調整画面
Ｓα：補正倍率、補正テーブル
２１：画像取得部
２２：階調圧縮部
２３：領域指定部
２４：カラーマップ取得部
２５：カラーマップ調整部
２６：補正テーブル作成部
２６ａ：幾何情報取得部
２６ｂ：幾何情報補正部
２６ｃ：補正倍率算出部
２７：ボクセル作成部
２８：レンダリング部
２８ａ：第１レンダリング部
２８ｂ：第２レンダリング部
２９：データ出力部
３０：パラメータ調整部
３１：ＭＰＲ生成部
３２：ＭＰＲ表示部
４０：パラメータ調整画面
Ｓα：補正倍率、補正テーブル
【図1】
【図2】
【図3】
【図4】
【図5】
【図6】
【図7】
【図8】
【図9】
【図10】
【図11】
【図12】
【図13】
【図14】
【図15】
【図16】
【図17】
【図18】
【図19】
【図20】
【図21】
【図22】
【図23】
【図24】
【図25】
【図26】
【図27】
【図28】
【図29】
【図30】
【図31】
【図32】
【図33】
【図34】
【図35】
【図36】
【図37】
【図38】
【図39】
