JP3762482B2

JP3762482B2 - ボリュームレンダリング法による陰影づけ方法及び装置

Info

Publication number: JP3762482B2
Application number: JP15835796A
Authority: JP
Inventors: 良洋後藤
Original assignee: Hitachi Medical Corp
Current assignee: Hitachi Healthcare Manufacturing Ltd
Priority date: 1996-06-19
Filing date: 1996-06-19
Publication date: 2006-04-05
Anticipated expiration: 2016-06-19
Also published as: JPH1011604A

Description

【０００１】
【発明の属する技術分野】
本発明は、ボリュームデータを含む複数の断層像から得る三次元画像の陰影づけに係り、特にボリュームレンダリング法による陰影づけ方法及び装置に関する。
【０００２】
【従来の技術】
ボリュームレンダリング法による陰影づけを行うと、淡い表示や半透明表示ができることが知られている。それは、光が物体内部まで減衰・反射しながら通過し、光量が一定の割合に減衰するまでに通過したところを全て画像化するという方法をとっているためである。
【０００３】
【発明が解決しようとする課題】
ボリュームレンダリング法による半透明表示の場合、いくつも見える臓器のうち、特定の臓器を明瞭（コントラストをつけて）にして、ほかの臓器は参考程度にうっすらと（コントラストをあまりつけず）表示したいときがある。例えば人間の頭に関して、頭は皮膚と脳と骨とで構成されるが皮膚と脳をうっすらと描画して、骨を明瞭に描画する如き例がある。
【０００４】
かかる従来のボリュームレンダリング法の一例（周藤安造著．「医学における三次元画像処理」の１０４ページ〜１１０ページ。１９９５年３月１０日．株式会社コロナ社発行）を以下に示す。
ボリュームレンダリングによる三次元画像の画素値を求めるには、先ず傾斜を考慮した画素値を下記の式に従って求める。
【数１】

ここで、Ｇｒａｄは傾斜から得られる画素値を示すためのシンボル、ｆ（ｉ、ｊ、ｋ）は、座標（ｉ、ｊ、ｋ）における画素値、Ｇｒａｄｆ（ｉ、ｊ、ｋ）は傾斜から求まる座標（ｉ、ｊ、ｋ）での画素値を示す。更に右辺での（ｆ(ｉ＋１、ｊ、ｋ)−ｆ(ｉ−１、ｊ、ｋ)）／２とは、座標ｉの前後の座標ｉ＋１、ｉ−１での画素値の差分を求め、これを２で除算（即ち１つの座標ｉでの傾斜に割り当てるための処理）するとの意味をなす。他の座標ｊ、ｋでも同様である。
【０００５】
簡単のため、数１に関して、座標（ｉ、ｊ、ｋ）をｘ_iとおく。そして勾配を考慮した任意の位置ｘ_iでのボクセルの不透明度α（ｘ_i）を、次式で表す。ここで、α_a、α_b、α_cは典型的な組織（例えば、ａ、ｂ、ｃの３つの組織例）毎に固有な不透明度を示す定数であり、テーブルに登録しておく。ｆ_a、ｆ_b、ｆ_cは複数の異なる生体組織（ａ、ｂ、ｃの３つの組織例）の境界におけるスカラー値（一般にボクセル値）であり、例えばｆ_a≦ｆ_b≦ｆ_cである。従って、ｆ_aとα_a、ｆ_bとα_b、ｆ_cとα_cとはそれぞれ複数の異なる生体組織ａ、ｂ、ｃの境界におけるスカラ値とその値に対する不透明度を示すことになる。ｆ（ｘ_i）と、ｆ_a、ｆ_b、ｆ_cとの大小の対比で不透明度α（ｘ_i）は以下の４つの形態となる。
【０００６】
【数２】

三次元画像のためのレンダリング処理は次のようになる。投影線毎にその投影線に沿って奥行き方向に並んだボクセルについて、前方の光と自らが発する光、あるいは光源に照らされて反射する光（前記不透明度α（ｘ_i）よって定まる）とを合成して、次々とボクセルに渡していき、最後のボクセルから出る光がレンダリング画像のビクセル値（画素値）になる。なお、反射係数はあらかじめテーブルで与えておく。
こうしたボリュームレンダリング法では、前述の如き骨を明瞭に、皮膚と脳とをうっすらと描く如き描き分けは不可能であった。
【０００７】
本発明の目的は、こうした課題の解消をはかり、必要とする部位の画素値を任意に設定可能な陰影づけ方法及び装置を提供するものである。
【０００８】
【課題を解決するための手段】
本発明は、ボリュームデータを含む複数の断層像から三次元画像を構成し、その際に前記断層像の画素値に対応する情報を含ませて陰影を施すようにする陰影づけ装置において、ボリュームレンダリング法の要素である陰影づけの勾配を画素点の位置と投影面との距離に関係づけた値で冪乗する手段と、この冪乗で得た勾配に従って各ボクセルの不透明度を算出する手段と、かくして得た不透明度に従って、各投影線上でのボリュームレンダリング画素値を算出する手段と、該算出したボリュームレンダリング画素値を投影面上の画素値として表示する手段と、より成るボリュームレンダリング法による陰影づけ装置を開示する。
【０００９】
更に本発明は、ボリュームデータを含む複数の断層像から三次元画像を構成し、その際に前記断層像の画素値に対応する情報を含ませて陰影を施すようにする陰影づけ方法において、第１、第２、第３の算出手段と表示手段とを具えて、第１の算出手段によってボリュームレンダリング法の要素である陰影づけの勾配を画素点の位置と投影面との距離に関係づけた値で冪乗するステップと、第２の算出手段によって、この冪乗で得た勾配に従って各ボクセルの不透明度を算出するステップと、第３の算出手段によってこの不透明度に従って、各投影線上でのボリュームレンダリング画素値を算出するステップと、表示手段によって該算出したボリュームレンダリング画素値を投影面上の画素値として表示するステップと、より成るボリュームレンダリング法による陰影づけ方法を開示する。
【００１１】
【発明の実施の形態】
本実施の形態では、前記（数１）はそのまま使用し、前記（数２）を下記の如き数式とする。
【数３】

即ち、傾斜画素値Ｇｒａｄｆ（ｘ_i）（絶対値）をδ乗化するようにした。ここで、値δとは、画像の性質（又は陰影化の仕方）に従って変化させるものとする。
【００１２】
図２には、画素点と投影面との距離ｒによって値δを特別の形で変化させた例である。形態Ａは、
ｒ＜ｒ₁では、δ＝δ₀₁、
ｒ₁≦ｒ≦ｒ₃では、δを三角形の形態で変更、
ｒ₃＜ｒでは、δ＝δ₀₁、
とした例である。
形態Ｂは、
ｒ＜ｒ₂では、δ＝δ₀₂、
ｒ₂≦ｒ≦ｒ₄では、δ＝３、
ｒ₄＜ｒでは、δ＝δ₀₂、
とした例である。形態Ａにするか形態Ｂにするかは、描きたい陰影化による。また、形態Ａ、Ｂ以外にも、距離ｒに応じて種々の形態を実現できる。形態ＡとＢとを比較すれば、
・形態Ａは形態Ｂよりも距離ｒの手前側でδを大きく設定している。
・形態Ａは三角形で且つピークが大きく、α₀₁＜α₀₂であるのに対し、形態Ｂは矩形でピークも小さく、α₀₂＞α₀₁である。
・従って、形態Ａを選べば投影面から近いところ（ｒ₁〜ｒ₂）で且つコントラストも三角形に従った集中形のものとなり、形態Ｂを選べば投影面から遠いところ（ｒ₂〜ｒ₄）で且つコントラストも矩形に従った多少の幅を持つものとなる。
また形態Ａに比べて形態Ｂはδ₀₂＞δ₀₁であることから、投影面から近い位置及び投影面から遠い位置では、明るく陰影化される。
こうした距離ｒと値δとの関係はメモリテーブルに格納しておき、必要に応じて読み出して使用する。
尚、図２での投影面とは、三次元画像の投影先の面のことであり、投影法には平行投影法や中心投影法等がある。投影面に投影する画像について陰影化をはかるが、かかる陰影化に際して、図２を利用することにしたのである。
【００１３】
図３は、ＣＴ値（画素値）及びＣＴ値の勾配に従って値δを変化させる例である。ＣＴ値をｘ軸、ＣＴ値の勾配をｙ軸とし、δをｚ軸とした。値δは、（ｘ₀、ｙ₀）でピーク値を呈し、その近傍でヒストグラム的な形状Ｃをなすように設定してある。ピーク値及びヒストグラムの形態及びその大きさは求めたい陰影化によって定める。
こうした図３の形態Ｃを採用することで、例えば皮膚のＣＴ値と脳のＣＴ値は互いに近いために区別しにくいが、境界での勾配は大きく異なる（皮膚は空気と接するので勾配は大きい）ので、皮膚と脳とを区別して表示できる利点を持つ。
【００１４】
図３のデータもメモリテーブルとして格納することで、読み出し自在に使用できる。テーブルに格納するに際しては、例えば、
ｘ＝ｘ₀、ｙ＝ｙ₀では、
δ＝Ｃ₂
ｘ₁≦ｘ≦ｘ₂、ｙ₁≦ｙ≦ｙ₂（但し、ｘ₁＜ｘ₀＜ｘ₂、ｙ₁＜ｙ₀＜ｙ₂とし、ｘ、ｙはｘ≠ｘ₀、ｙ≠ｙ₀とする）では、
δ＝Ｃ₀＋Ｃ₁／（ｘ²・ｙ²）
ｘ＜ｘ₁又はｘ＞ｘ₂、ｙ＜ｙ₁又はｙ＞ｙ₂では、
δ＝Ｃ₀
の如くする。ここで、Ｃ₀、Ｃ₁、Ｃ₂は定数である。Ｃ₀＜Ｃ₂である。
ヒストグラムとしたが、三角形や矩形等その形態は種々設定可能である。
【００１５】
図１には、図２又は図３を利用して得た画像の表示例を示す。表示画面上に、脳１と、皮膚３と骨２とを重ねて表示した。特に、骨２を強く表示し、脳１と皮膚３とを柔らかく表示させる陰影化をはかっている。
【００１６】
図４は、本実施の形態の装置例を示す。共通バス１８に、ＣＰＵ１０、主メモリ１１、磁気ディスク１２、表示メモリ１３、ＣＲＴ１４、コントローラ１５、マウス１６、キーボード１７を接続した構成である。磁気ディスク１２にボリュームレンダリング対象の画像（ＣＴ画像やＭＲ画像）が格納されており、これをキーボード１７やマウス１６の指示により主メモリ１１に送り、ＣＰＵ１０によってボリュームレンダリング処理を行う。この処理によって得られた陰影化された擬似三次元画像は表示メモリ１３に一時記憶されて、ＣＲＴ１４に表示される。その表示例が図１の如きものである。ここで、前記の参考文献に従うと、三次元画像のためのレンダリング処理は次のようになる。投影線に沿って並んだボクセルについて、前方の光と自らが発する光、あるいは光源に照らされて反射する光（前記不透明度によって定まる）とを合成して、次々とボクセルに渡していき、最後のボクセルから出る光がレンダリング画像のビクセル値（画素値）になる。尚、反射係数はあらかじめテーブルで与えておく。
尚、メモリテーブルに図２や図３の形態を格納させるとしたが、計算によって求めてもよく、又はキーボード１７からの入力によってもよい。
【００１７】
【発明の効果】
本発明によれば、半透明表示の場合、いくつも見える臓器のうち、特定の臓器を明瞭（コントラストをつけて）にして、他の臓器は参考程度にうっすらと（コントラストをあまりつけず）表示できるようになる。
【図面の簡単な説明】
【図１】本発明の陰影表示例図である。
【図２】本発明の陰影処理のためのδの設定例図である。
【図３】本発明の陰影処理のためのδの他の設定例図である。
【図４】本発明の陰影づけ装置の構成例図である。
【符号の説明】
１脳の画像
２骨の画像
３皮膚の画像

Claims

ボリュームデータを含む複数の断層像から三次元画像を構成し、その際に前記断層像の画素値に対応する情報を含ませて陰影を施すようにする陰影づけ装置において、ボリュームレンダリング法の要素である陰影づけの勾配を画素点の位置と投影面との距離に関係づけた値で冪乗する手段と、この冪乗で得た勾配に従って各ボクセルの不透明度を算出する手段と、かくして得た不透明度に従って、各投影線上でのボリュームレンダリング画素値を算出する手段と、該算出したボリュームレンダリング画素値を投影面上の画素値として表示する手段と、より成るボリュームレンダリング法による陰影づけ装置。
ボリュームデータを含む複数の断層像から三次元画像を構成し、その際に前記断層像の画素値に対応する情報を含ませて陰影を施すようにする陰影づけ方法において、第１、第２、第３の算出手段と表示手段とを具えて、第１の算出手段によってボリュームレンダリング法の要素である陰影づけの勾配を画素点の位置と投影面との距離に関係づけた値で冪乗するステップと、第２の算出手段によって、この冪乗で得た勾配に従って各ボクセルの不透明度を算出するステップと、第３の算出手段によってこの不透明度に従って、各投影線上でのボリュームレンダリング画素値を算出するステップと、表示手段によって該算出したボリュームレンダリング画素値を投影面上の画素値として表示するステップと、より成るボリュームレンダリング法による陰影づけ方法。