JP4709177B2 - 三次元画像処理装置および方法並びにプログラム - Google Patents

Description

本発明は医用画像のボリュームレンダリング、詳しくは被写体の形状を表す形態画像と被写体の厚さや特性等の機能を表したマッピング画像を用いてボリュームレンダリングを実行する画像処理装置および方法並びにプログラムに関する。

ＣＴ装置、ＭＲＩ装置、超音波診断装置（エコー）等の発展により、詳細なレベルでの医用断面画像の取得が可能になったが、それと同時に、一回の検査で撮影される画像数も膨大になったため、撮影された断面画像を一枚ずつ観察することによって診断する従来の方法では、膨大な時間がかかるようになった。また、断面画像に基づいて被写体の三次元的構造を把握するには、十分な経験を要することから、読影者の経験の差によって、読影精度にも差が生じてしまっていた。

そこで、三次元コンピュータグラフィックスの技術を応用して、多数の二次元断面画像に基づいて、二次元平面に立体的に可視化した擬似三次元画像を生成する画像処理が行われるようになった。

このような擬似三次元画像を生成する技術としては、ボリュームレンダリング法が知られている。これは、三次元画像データを構成する各画素（ボクセルデータ）に対して設定された不透明度（オパシティ）と輝度値とに基づき、視線に沿った各探索点におけるこれらの値をサンプリングし、加算していくことによって投影画素の出力画素値を求め、半透明な疑似三次元画像を生成する処理である。

一方、心筋シンチグラフィ（ＳＰＥＣＴ）によって直接撮影されるものやＣＴ装置、ＭＲＩ装置等から得られた解析結果により得られた被写体の各部の厚さや機能等を表した機能画像が知られている。

特許文献１に開示された手法においては、このような機能画像から注目したい一部を抽出し、X線ＣＴ装置、ＭＲＩ装置または超音波診断装置等により撮像された被写体の形状を表す形態画像と合成したボリュームデータを作成し、この合成ボリュームデータを対象にボリュームレンダリング手法を行う装置が開示されている。
特開２００６−２８８４９５号公報

しかしながら、被写体を表す三次元画像データの被写体の一部における形態画像と機能画像を合成し、ボリュームレンダリングを実行すると、被写体の一部の形態を表す形態画像の形状を保持できずに、機能画像の形状に依存した疑似三次元画像が生成されるというという問題があった。

本発明は、上記事情に鑑み、形態画像に対応する位置に存在する機能画像等の情報をマッピングすることが可能な三次元画像処理装置および方法並びにプログラムを提供することを目的とするものである。

本発明による三次元画像処理装置は、被写体を表す三次元画像を構成するボクセルデータの各画素に少なくとも不透明度を割り当てられてなる第１のボクセルデータを用いて、被写体の一部を表すマッピング画像を生成するマッピング画像生成手段と、被写体を表す三次元画像を構成する第２のボクセルデータのうち少なくとも被写体の一部を含む対象の形態を表した形態画像を生成する形態画像生成手段と、マッピング画像と形態画像の各位置を対応させて、各位置に対応するマッピング画像の不透明度に基づきボリュームレンダリングを実行することにより疑似三次元画像を生成する画像生成手段とを備えたものであることを特徴とするものである。

本発明による三次元画像処理方法は、被写体を表す三次元画像を構成するボクセルデータの各画素に少なくとも不透明度を割り当てられてなる第１のボクセルデータを用いて、被写体の一部を表すマッピング画像を生成する工程と、被写体を表す三次元画像を構成する第２のボクセルデータのうち少なくとも被写体の一部を含む対象の形態を表した形態画像を生成する工程と、マッピング画像と形態画像の各位置を対応させて、各位置に対応するマッピング画像の不透明度に基づきボリュームレンダリングを実行することにより疑似三次元画像を生成する工程とを備えることを特徴とするものである。

本発明のプログラムは、コンピュータを、被写体を表す三次元画像を構成するボクセルデータの各画素に少なくとも不透明度を割り当てられてなる第１のボクセルデータを用いて、被写体の一部を表すマッピング画像を生成するマッピング画像生成手段と、被写体を表す三次元画像を構成する第２のボクセルデータのうち少なくとも被写体の一部を含む対象の形態を表した形態画像を生成する形態画像生成手段と、マッピング画像と形態画像の各位置を対応させて、各位置に対応するマッピング画像の不透明度に基づきボリュームレンダリングを実行することにより疑似三次元画像を生成する画像生成手段として機能させることを特徴とするものである。

「マッピング画像」とは、三次元コンピュータグラフィックスにおいて、ボリュームレンダリング等の際に、物体の表面に模様や質感などの様々な効果を貼る際に用いられる被写体の一部の特性等を表す画像をいう。具体例としては、機能画像が考えられる。

「機能画像」とは、被写体の一部における各位置の機能を視認できるように表した画像を意味する。

「ボリュームレンダリング」とは、三次元医用画像を構成する各画素（ボクセルデータ）に対して設定された不透明度（オパシティ）と輝度値等に基づき、視線に沿った各探索点におけるこれらの値をサンプリングし、加算していくことによって投影画素を生成する処理を意味する。具体例としては、レイキャスト法が考えられる。

「輝度値」とは三次元画像を構成する画素値、三次元画像を構成する画素値に対応した色情報、探索点における画像勾配と光源の関係から求まる照度等から計算される探索点における輝度の値を意味する。

「レイキャスト法」とは、物体に対して投影面から仮想的な光線を照射し、ボクセル値に対応した不透明度や輝度値等に基づいて、物体内部からの仮想的な反射光による三次元画像を作成することにより、投影面に物体内部の三次元構造を透視する投影画像を生成する手法を意味する。

本発明の三次元画像処理装置および方法並びにプログラムによれば、マッピング画像と形態画像の各位置を対応させて、各位置に対応するマッピング画像の不透明度に基づきボリュームレンダリングを実行するので、形態画像に対応する位置に存在するマッピング画像の情報をマッピングし、マッピング画像の形状に依存することなく形態画像の形状通りに疑似三次元画像の作成が可能となる効果を奏する。

以下、図面を参照して本発明の三次元画像処理装置の実施形態を詳細に説明する。図１は本発明の三次元画像処理装置の概略構成図である。図１のような三次元画像処理装置１の構成は、補助記憶装置に読み込まれた三次元画像処理プログラムをコンピュータ上で実行することにより実現される。このとき、この三次元画像処理プログラムは、ＣＤ−ＲＯＭ等の記憶媒体に記憶され、もしくはインターネット等のネットワークを介して配布され、コンピュータにインストールされる。

本発明の三次元画像処理装置１は、被写体を撮影して得た三次元医用画像を構成する第１のボクセルデータ１００を記憶する第１の記憶手段１０と、第１のボクセルデータ１００を用いて心臓の各位置における機能を表した機能画像２００を生成する機能画像生成手段２０と、前記被写体を撮影して得た三次元画像を構成する第２のボクセルデータ１１０を記憶する第２の記憶手段３０と、被写体を撮影して得た三次元画像を構成する第２のボクセルデータ１１０のうちまた、被写体を表す画像の一部である心臓の領域であるボクセルデータを用いて、前記心臓の形態を描画した形態画像を生成する形態画像生成手段４０と、第１のボクセルデータ１００と第２のボクセルデータ１１０を位置合せする位置合せ手段５０と、形態画像２１０と機能画像２００とを用いてボリュームレンダリング手法を実行する画像生成手段６０と、画像表示手段６０で作成された疑似三次元画像を表示する表示手段２と、更にオペレータによりデータ入力を可能とする入力手段３を備える。

第１、第２のボクセルデータ１００、１１０は、ＣＴ（Ｃｏｍｐｕｔｅｄ Ｔｏｍｏｇｒａｐｈｙ）装置やＭＲＩ（Ｍａｇｎｅｔｉｃ Ｒｅｓｏｎａｎｃｅ Ｉｍａｇｉｎｇ）装置などにより撮影されて得られた三次元画像である。第１のボクセルデータ１００は、機能画像２００を生成するために用いられ、第２のボクセルデータ１１０は、形態画像２１０を生成するために用いられる。

第１、第２の記憶手段１０、３０は、大容量の記憶装置であり、ハードディスクや画像サーバなどである。被写体をＣＴ装置やＭＲＩ装置などを用いて撮影した第１、第２のボクセルデータ１００、１１０が記憶される。

マッピング画像とは、三次元コンピュータグラフィックスにおいて、ボリュームレンダリング等の際に、物体の表面に模様や質感などの様々な効果を貼る（マッピング）際に用いられる被写体の一部の特性等を表す画像をいう。マッピング画像の具体例として、機能画像が考えられる。以下、機能画像２００について説明する。

機能画像２００は、心臓の動きや心室の内径や心筋の厚さなどに応じて、心臓の機能の評価を表す評価値の分布（例えば、色や数値等）を表す画像である。具体的には、あるフェーズにおける心室の径を表した心室内径画像、拡張期における心室の径を表した拡張末期心室内径の画像、収縮期における心室の径を表した収縮末期心室内径の画像、区画された領域毎の駆出率を表した局所駆出率の画像、あるフェーズでの心筋の厚みを表した壁厚の画像、拡張期における心筋の厚さを表した拡張末期壁厚の画像、収縮期における心筋の厚さを表した収縮末期壁厚の画像、拡張期の心筋の厚さから収縮期の心筋の厚さの差を表した壁厚変化量の画像、拡張期の心筋の厚さをA、収縮期の心筋の厚さをBとした時、B-A/Aで計算された値を表した壁厚増加率の画像、拡張期の心室の径から収縮期の心室の径の差を表した壁運動量の画像、心筋シンチグラフィの画像などがある。

心筋シンチグラフィを表す機能画像２００は、心筋に集まる薬を被検者の腕から注射して、その薬の分布を体の外から測定したデータを画像として表したもので、薬を変えることで、心筋の血流状態、心筋組織の代謝、神経の働きなどを表すことができる。

また、機能画像２００で心臓の動きの評価を表す場合には、動いている心臓を複数フェーズ撮影して画像間の違いから心臓の機能を評価する値を得て機能画像２００を作成する。このように複数フェーズ撮影する場合には、被曝することがないＭＲＩ装置で複数フェーズ撮影した画像から機能画像２００を生成するのが望ましい。

一方、形態画像２１０は心臓の形態を表す画像であるが、各臓器の構造が明瞭に表れるＣＴ装置で撮影した三次元画像から生成するのが好ましい。

また、形態画像２１０は機能画像２００の被写体の一部を含む対象の形態を表すことが好ましい。

また、形態画像２１０の作成に用いる第２のボクセルデータには、形態を詳細に表すために、ＣＴ装置やＭＲＩ装置などで細かいピッチ（例えば、１〜３ｍｍ間隔）で撮影した断層画像から生成したボクセルデータを用いるのが好ましい。一方、機能画像２００の作成するために心臓の動きに応じて複数のフェーズで撮影する必要がある場合は、細かいピッチで撮影を行ったのでは、複数フェーズを撮影することができない場合があるので、第１のボクセルデータは、ＭＲＩ装置などで広いピッチ（例えば、５〜１０ｍｍ間隔）で撮影した断層画像から生成したボクセルデータを用いることが多い。

位置合せ手段５０は、第１のボクセルデータ１００の心臓と第２のボクセルデータ１１０の心臓について位置合せを行う。例えば、第１のボクセルデータ１００がＭＲＩ装置で心臓を撮影した画像であり、第２のボクセルデータ１１０がＣＴ装置で心臓を撮影した画像である場合には、同じ被写体を撮影した場合であっても、全く同じ位置になるように撮影することは難しい。そこで、機能画像２００の対応する心臓の位置に心臓の周囲に存在する血管の形態を投影するためには、第１のボクセルデータ１００の心臓の形状と第２のボクセルデータ１１０の心臓の形状とを位置合せを行う。

この位置合せはレジストレーションと呼ばれる種々の手法を用いることができる。ここで、第１のボクセルデータ１００の画像Ａの点ｘ_Ａを第２のボクセルデータ１１０の画像Rの点ｘ_Ｒに変換する場合について具体的に考える。

変換関数をTとすると、

と表すことができる。

心臓の形状を剛体として考えた場合には、三次元空間内のX軸方向、Y軸方向、Z軸方向の各方向への移動（ｔ_ｘ、ｔ_ｙ、ｔ_ｚ）と、各軸の回転（α、β、γ）により変換関数が表され、下式のようになる。

心臓の形状は変化するため非剛体として考えた場合には、自由度が増えるためＴが多項式やスプライン関数などを含む複雑な関数で表される。

変換関数は、具体的には、心臓の特徴となる解剖学的な点をいくつか（例えば、三点）検出して、この点から求めるようにしてもよい。あるいは、一つの画像の輪郭を決定し、もう一つの画像の表面に相当する点列の距離が最も小さくなるように、フィティングを繰り返して変換関数を求めることもできる。あるいは、画像中の全てのボクセルの画素値を用いて類似度を調べることによって重ね合わせるようにして求めてもよい。このように画像中の全ての画素値を用いるようにすると、画素上の雑音成分が打ち消しあうので比較的安定した結果を得ることができる。（詳細は、日本放射線技術学会雑誌、第５３巻、第１号、Ｐ６０〜Ｐ６５、２００３年１月などを参照）。

位置合せは、上述のいずれのレジストレーション手法を用いてもよい。

このように座標変換関数が定義された上で第２のボクセルデータ１１０は独立したデータであるため、第２のボクセルデータ１１０の可視化条件はそのままで第１のボクセルデータ１００を変化させることで、同じ形状に対し異なった情報をマッピングすることが可能となる。

画像生成手段６０は、位置合せ手段５０によって算出された第１、第２のボクセルデータ１００，１１０の対応情報に基づき形態画像２１０と機能画像２００の各位置を対応させ、各位置に対応する機能画像の不透明度に基づき形態画像と機能画像のそれぞれの輝度値（ボクセルデータの各画素値に対し、割り当てられた色情報等）を用いてボリュームレンダリングを実行し、疑似三次元画像を生成するものである。

例えば、ボリュームレンダリングの手法として、具体的にはレイキャスト法を用いる。レイキャスト法は、濃度や密度で定義されるボリュームデータに対し、画面上から伸ばした視線上にある濃度／密度データを計算し内部情報も反映して画像を生成するものである。

被写体の三次元医用画像である形態画像２１０を取得する。この形態画像２１０の各ボクセルの位置は、被写体の左右方向をｘ軸、前後方向をｙ軸、上下方向をｚ軸とする三次元座標系で定義され、各ボクセルの画素値は、そのボクセルの位置の座標と関連づけられるものである。

次に、初期設定値を用いるか、またはオペレータによる入力手段３（例えばキーボードやマウス等）からの入力等により設定された視点Ｅ、光源Ｓ、投影面Ｆ（大きさ、位置、画素数）により、図２に示すように、その視点Ｅと投影面Ｆ上の投影画素の各々とを結ぶ複数の視線Ｅ_k（k＝１、２、…、Ｌ；Ｌは、視線の数）に沿って三次元医用画像Ｖを所定の間隔でサンプリングした複数の探査点ｖ_ｉ（ｉ＝１、２、…、ｎ）を設定することが可能となる。

各探査点ｖ_ｉの輝度値と不透明度とを決定した後、各視線Ｅ_kに沿って、その視線Ｅ_k上の各探索点v_ｉ（ｉ＝１、２、…、ｎ）の輝度値ｃ(v_ｉ)と不透明度α(v_ｉ)とに基づいて、その視線Ｅ_kが通る投影面上の投影画素の出力画素値Ｃを、次の式（３）により決定する。

このように、形態画像２１０の各画素値に対して不透明度を割り当てることで被写体の対象となる心臓を抽出するような画像の生成が可能となる。この不透明度が可視化対象となるデータを意味し、対象となるボリュームデータの形態を定義することになる。

これらの手法の考え方を応用して、本発明は以下の手法を提案するものである。

機能画像２００と形態画像２１０の少なくとも２つを入力とし、形態画像２１０を対象にレイキャスト法を実行する。可視化対象となる不透明な画素値を持つ場合、上述の位置合せ手段５０で対応する位置を求め、その機能画像２００から対応する位置での輝度値を取得して可視化をおこなう。

具体的には、図３に示すように、その視点Ｅと投影面Ｆ上の投影画素の各々とを結ぶ複数の視線に沿ってＶ１を所定の間隔でサンプリングした複数の探査点である形態画像２１０に関するｖ_ｉ（ｉ＝１、２、…、ｎ）とそれに対応するマッピング用の機能画像２００に関する探査点ｖ´_ｉ（ｉ＝１、２、…、ｎ）を設定する。

形態画像２１０での各探査点における不透明度と輝度値を取得する。また、形態画像２１０の探査点の位置に対応する機能画像２００の位置での不透明度と輝度値を取得する。各探査点で機能画像２００から取得した不透明度をアルファ値として形態画像２１０の輝度値と機能画像２００の輝度値とのアルファブレンディングすることで、形態画像２１０と機能画像２００をブレンドした出力画素値Ｃを式（４）により決定する。

また、表示手段２は画像生成手段６０を介して、図４に示すような機能画像Ａ１と形態画像Ｂ１を重ね合わせされた疑似三次元画像を表示する。図４に示すように、重ね合わせて表示することによって、血管の位置から心臓のどのあたりの機能が正常であるか異常であるかなどを画像表示手段２により観察することができる。

このように形態画像２１０に対応する位置に存在する機能画像２００の情報をマッピングするため、形態画像での形状はそのままに機能画像から得られる情報が付加された形での可視化が可能となる。

なお、オペレータにより入力手段３を介して、機能画像２００を他の機能画像に変更することも可能である。そのような場合、形態画像２１０と他の機能画像が重畳したような疑似三次元画像を表示手段２に表示することが可能となる。

次に、三次元画像処理装置１を用いて、機能画像２００上に形態画像２１０を重ねた画像を表示するときの流れについて図５のフローチャートに従って説明する。

まず、被写体の心臓をＭＲＩ装置で撮影し、撮影した画像を第１のボクセルデータ１００として第１の記憶手段１０に記憶する（Ｓ１００）。さらに、同じ被写体をＣＴ装置で撮影した画像を第２のボクセルデータ１１０として第２の記憶手段２０に記憶する（Ｓ１０１）。

機能画像生成手段２０で、第１のボクセルデータ１００より心臓の機能を表す機能画像２００を生成する（Ｓ１０２）。

次に、第２のボクセルデータ１１０を用いて、形態画像生成手段４０で形態画像２１０を生成する。（Ｓ１０３）。

位置合せ手段５０で、第１のボクセルデータ１００の心臓と第２のボクセルデータ１１０の心臓の形状を位置合せする。具体的には、例えば、心臓を剛体として考えて、ボクセルデータ１００，１１０の画素値を用いて相関係数を用いて評価し、機能画像の座標系から携帯画像の座標系へ変換するときの回転量、移動量などの変位量を求めることができる（Ｓ１０４）。

機能画像生成手段２０で生成された第１のボクセルデータ１００より心臓の機能を表示した機能画像２００に割り当てられた不透明度をアルファ値とし、アルファブレンド処理により、機能画像２００と形態画像２１０それぞれの輝度値等を用いて、画像生成手段６０で、ボリュームレンダリング処理を実行する（Ｓ１０５）。

画像生成手段６０は、上述のボリュームレンダリング処理を経て、機能画像２００と形態画像２１０とを重ね合うように疑似三次元画像を生成する（Ｓ１０６）。

また、第１のボクセルデータと第２のボクセルデータがＣＴ装置によって撮影されたものとＭＲＩ装置によって撮影されたものである場合や、共にＭＲＩ装置あるいはＣＴ装置によって撮影されたものであってもスライス間隔が異なる場合には、上述のレジストレーションの手法を用いて位置合せが行なわれるが、心筋シンチグラフィを対象とする場合は、画像の撮影目的が異なり形態画像と機能画像との関連が小さいため、画素値の相関を用いたレジストレーションでは上手く位置合せできないことがある。そこでこのような場合には、それぞれの画像から対応する特徴点を選択し、その点間でレジストレーションをおこない座標変換のための回転量、移動量を算出するようにしてもよい。

また、上述では機能画像２００の座標系を形態画像２１０の座標系に変換する場合について説明したが、心臓上の各位置が一致するように変換するのであれば、形態画像２１０の座標系を機能画像２００の座標系に変換しても、双方を変形するようにしてもよい。

あるいは、機能画像２００と形態画像２１０のいずれか一方の画像を生成し、位置合せの情報を用いて他方の画像を先に生成した画像の位置に対応させて描画するようにしてもよい。

以上、機能画像と形態画像の各位置を対応させて、各位置に対応する機能画像の不透明度に基づきボリュームレンダリングを実行することにより、形態画像に対応する位置に存在する機能画像の情報をマッピングし、機能画像の形状に依存することなく形態画像の形状通りに疑似三次元画像の作成が可能となる。

本発明は、形態画像の可視化に対してはボリュームレンダリング手法がベースとなっており、内部構造を保持したデータとして不透明な物体として可視化できるボリュームレンダリングの特徴を引き継いでいるため、機能画像からマッピングされる情報は形状の表面でなく内部にも行われるので、心臓の外観形状の表面のみにマッピングのみならず、より透明な不透明度を設定することで内部にもマッピングされた可視化が可能となる。

更にオペレータが入力手段３を介して、表示手段２に表示された画像データを可逆又は非可逆圧縮データ（例えば、ＪＰＥＧ形式）として出力するようにしてもよい。これにより大容量である医療画像を圧縮し、低容量の画像としてデータ転送やプレゼンテーション資料等への利用が容易になる。

本発明の疑似三次元画像生成装置の概略構成図 レイキャスト法における変換図 形態画像用座標系からマッピング画像用座標系への変換図 形態画像と機能画像を重ね合わせた画像の一例（その１） 本発明の疑似三次元画像を生成するまでを説明するためのフローチャート

符号の説明

１ 三次元画像生成装置
２ 表示装置
１０ 第１の記憶手段
２０ 機能画像生成手段
３０ 第２の記憶手段
４０ 形態画像生成手段
５０ 位置合せ手段
６０ 画像生成手段
１００ 第１のボクセルデータ
１１０ 第２のボクセルデータ
２００ 機能画像
２１０ 形態画像

Claims (8)

1. 被写体を表す三次元画像を構成するボクセルデータの各画素に少なくとも不透明度を割り当てられてなる第１のボクセルデータを用いて、前記被写体の一部を表すマッピング画像を生成するマッピング画像生成手段と、
   前記被写体を表す三次元画像を構成する第２のボクセルデータのうち少なくとも前記被写体の一部を含む対象の形態を表した形態画像を生成する形態画像生成手段と、
   前記マッピング画像と前記形態画像の各位置を対応させて、該各位置に対応する前記マッピング画像の不透明度に基づき、該マッピング画像の不透明度の値が大きいほど前記マッピング画像の輝度値の比率が高くなるように、前記形態画像の輝度値と前記マッピング画像の輝度値とをブレンドした輝度値を算出し、該ブレンドされた輝度値を用いてボリュームレンダリングを実行することにより疑似三次元画像を生成する疑似三次元画像生成手段とを備えるものであることを特徴とする画像処理装置。
2. 前記第２のボクセルデータは、各画素に不透明度が割り当てられたものであり、
   前記疑似三次元画像生成手段は、前記ブレンドされた輝度値と、前記第２のボクセルデータの各画素に割り当てられた不透明度との積算値とを用いて、前記第１のボクセルデータの各画素に割り当てられた不透明度を用いずに、前記ボリュームレンダリングを実行するものであることを特徴とする請求項１に記載の画像処理装置。
3. 前記疑似三次元画像生成手段が、レイキャスト法を用いて前記ボリュームレンダリングを実行するものであることを特徴とする請求項１または請求項２記載の画像処理装置。
4. 前記疑似三次元画像生成手段が、前記各位置に対応する前記マッピング画像の不透明度に基づきアルファ値を設定し、前記各位置に対応する前記マッピング画像と前記形態画像のそれぞれの輝度値に対してアルファブレンド処理を用いてボリュームレンダリングを実行するものであることを特徴とする請求項１から請求項３いずれかに記載の画像処理装置。
5. 前記疑似三次元画像生成手段は、前記マッピング画像と前記形態画像が重畳された前記疑似三次元画像データを作成することを特徴とする請求項１から請求項４いずれかに記載の画像処理装置。
6. 前記マッピング画像は、前記対象組織の機能を表した機能画像であることを特徴とする請求項１から請求項５いずれかに記載の画像処理装置。
7. 被写体を表す三次元画像を構成するボクセルデータの各画素に少なくとも不透明度を割り当てられてなる第１のボクセルデータを用いて、前記被写体の一部を表すマッピング画像を生成する工程と、
   前記被写体を表す三次元画像を構成する第２のボクセルデータのうち少なくとも前記被写体の一部を含む対象の形態を表した形態画像を生成する工程と、
   前記マッピング画像と前記形態画像の各位置を対応させて、該各位置に対応する前記マッピング画像の不透明度に基づき、該マッピング画像の不透明度の値が大きいほど前記マッピング画像の輝度値の比率が高くなるように、前記形態画像の輝度値と前記マッピング画像の輝度値とをブレンドした輝度値を算出し、該ブレンドされた輝度値を用いてボリュームレンダリングを実行することにより疑似三次元画像を生成する工程とを備えることを特徴とする画像処理方法。
8. コンピュータを、
   被写体を表す三次元画像を構成するボクセルデータの各画素に少なくとも不透明度を割り当てられてなる第１のボクセルデータを用いて、前記被写体の一部を表すマッピング画像を生成するマッピング画像生成手段と、
   前記被写体を表す三次元画像を構成する第２のボクセルデータのうち少なくとも前記被写体の一部を含む対象の形態を表した形態画像を生成する形態画像生成手段と、
   前記マッピング画像と前記形態画像の各位置を対応させて、該各位置に対応する前記マッピング画像の不透明度に基づき、該マッピング画像の不透明度の値が大きいほど前記マッピング画像の輝度値の比率が高くなるように、前記形態画像の輝度値と前記マッピング画像の輝度値とをブレンドした輝度値を算出し、該ブレンドされた輝度値を用いてボリュームレンダリングを実行することにより疑似三次元画像を生成する疑似三次元画像生成手段として機能させるプログラム。