(57)【特許請求の範囲】
【請求項１】
外観比Ａを有するボクセルを有している体積測定データの表面を補間なしに表示する方法であって、ボクセルの値の各々は、特定された位置及びインデックスにおいて存在しており、
（ａ） 中心ボクセルを選択する工程と、
（ｂ） 該中心ボクセルを所定の閾値に対して試験する工程と、
（ｃ） 該所定の閾値を上回る各々の中心ボクセルに対して、
（ｉ） 前記中心ボクセルと異なる平面上に存在する複数の対角線ボクセルあって、前記中心ボクセルに対し対角の位置に存在し、相互に隣接しない複数の対角線ボクセル及び該ボクセルの値を決定する工程と、
（ｉｉ） 各々の対角線ボクセルの値を前記所定の閾値と比較する工程と、
（ｉｉｉ） 少なくとも１つの対角線ボクセルの値が前記閾値を下回る場合に、前記中心ボクセルを表面ボクセルとして分類する工程と、
（ｉｖ） 前記中心ボクセルを跨っているボクセルの値の間の差分から各々の表面ボクセルに対する表面法線ベクトルを算出する工程と、
（ｄ） 前記中心ボクセルとしての複数の異なるボクセルに対して工程（ａ）〜工程（ｃ）を繰り返して、複数の表面ボクセル及び表面法線ベクトルを得る工程と、
（ｅ） 該表面ボクセルの所望のビューのビューイング角度を受け取る工程と、
（ｆ） 所望のビューイング角度（方位角、仰角γ）での該表面ボクセル及び前記表面法線ベクトルをスクリーン・ピクセル・データにレンダリングする工程であって、該スクリーン・ピクセル・データの各々は、前記外観比Ａを考慮して修正されたビューイング角度から見えるように、スクリーン位置（ｕ，ｖ）、ピクセル・カラー及び強度を有している、レンダリングする工程と、
（ｇ） 前記ピクセル・データを画像としてオペレータに表示する工程とを備えた体積測定データの表面を補間なしに表示する方法。
【請求項２】
前記表面法線を算出する工程は、
（ａ） 前記中心ボクセルの左方に隣り合ったボクセルを左方ボクセルとして選択する工程と、
（ｂ） 前記中心ボクセルの右方に隣り合ったボクセルを右方ボクセルとして選択する工程と、
（ｃ） 前記左方ボクセルを前記右方ボクセルから減算して、「ｘ」の差分を決定する工程と、
（ｄ） 前記中心ボクセルの前方に隣り合ったボクセルを前方ボクセルとして選択する工程と、
（ｅ） 前記中心ボクセルの後方に隣り合ったボクセルを後方ボクセルとして選択する工程と、
（ｆ） 前記前方ボクセルを前記後方ボクセルから減算して、「ｚ」の差分を決定する工程と、
（ｇ） 前記中心ボクセルの下方に隣り合ったボクセルを下方ボクセルとして選択する工程と、
（ｈ） 前記中心ボクセルの上方に隣り合ったボクセルを上方ボクセルとして選択する工程と、
（ｉ） 前記下方ボクセルを前記上方ボクセルから減算して、「ｙ」の差分を決定する工程と、
（ｊ） ｇ＝｛（「ｘ」の差分）2＋（「ｙ」の差分）2＋（「ｚ」の差分）2｝1/2
Ｎｘ＝（「ｘ」の差分）／ｇ
Ｎｙ＝（「ｙ」の差分）／ｇ
Ｎｚ＝（「ｚ」の差分）／（Ａ×ｇ）
に従って、「ｘ」、「ｙ」及び「ｚ」の差分をベクトル成分Ｎｘ、Ｎｙ及びＮｚに規格化する工程と、
（ｋ） ベクトル成分Ｎｘ、Ｎｙ及びＮｚを表面法線ベクトルとして、関連するボクセル・インデックス又は位置と共に記憶する工程と
を含んでいる請求項１に記載の体積測定データの表面を補間なしに表示する方法。
【請求項３】
前記表示する工程の前に、
Ｅs ＝（（Ａｃｏｓ（γ））2＋ｓｉｎ2（γ））1/2
（ｕ，ｖ’）＝（ｕ，Ｅs ×ｖ）
に従って、画像をスケーリングするように新しい位置（ｕ，ｖ’）に対して前記ピクセル・データの位置（ｕ，ｖ）を補う工程を更に含んでおり、ここで、γは、前記所望のビューイング角度のうちの仰角である請求項１に記載の体積測定データの表面を補間なしに表示する方法。
【請求項４】
前記レンダリングする工程の前に、
Ψ＝ｔａｎ-1（（１／Ａ）ｔａｎ（γ））
に従って、前記仰角γを有効仰角Ψに調節する工程を更に含んでおり、該有効仰角Ψは、前記表面ボクセルをレンダリングするために仰角γの代わりに用いられる請求項１に記載の体積測定データの表面を補間なしに表示する方法。
【請求項５】
特定された位置及びインデックスにおける各々のボクセル値を有している体積測定データの表面を補間なしに表示するグラフィックス・システムであって、
（ａ） 前記体積測定データを記憶するように構成されている体積測定メモリと、
（ｂ） それぞれが該体積測定メモリに接続されており、中心ボクセルのインデックスを読み取り、該中心ボクセルの値を所定の閾値と比較し、該所定の閾値を上回る各々の中心ボクセルに対して、前記中心ボクセルと異なる平面上に存在する複数の対角線ボクセルあって、前記中心ボクセルし対角の位置に存在し、相互に隣接しない複数の対角線ボクセルを決定すると共に、該対角線ボクセルの値を所定の閾値と比較する複数の比較器と、
（ｃ） 表面ボクセルのインデックスと、各々の表面ボクセルに対する表面に垂直なベクトルとを記憶することが可能な表面／法線メモリと、
（ｄ） 前記比較器に接続されおり、各々の中心ボクセルに対するインデックスを選択すると共に前記体積メモリ内の中心ボクセルのインデックスを記憶し、前記比較器の出力を受け取り、当該中心ボクセルの少なくとも１つの対角線ボクセルの値が中心ボクセルの試験とは閾値の相容れない側にある場合に中心ボクセルを表面ボクセルとして分類して、表面ボクセルのインデックスを前記表面／法線メモリに記憶する論理装置と、
（ｅ） 前記中心ボクセルを「ｘ」方向に跨っているボクセルのインデックスを算出し、該ボクセルの値を減算して、「ｘ」の差分を得る「ｘ」減算装置と、
（ｆ） 前記中心ボクセルを「ｙ」方向に跨っているボクセルのインデックスを算出し、該ボクセルの値を減算して、「ｙ」の差分を得る「ｙ」減算装置と、
（ｇ） 前記中心ボクセルを「ｚ」方向に跨っているボクセルのインデックスを算出し、該ボクセルの値を減算して、「ｚ」の差分を得る「ｚ」減算装置と、
（ｈ） 前記減算装置に接続されており、「ｘ」、「ｙ」及び「ｚ」の差分を受け取り、表面ベクトル成分（Ｎx ，Ｎy ，Ｎz ）に規格化して、該ベクトル成分を前記表面／法線メモリに記憶する規格化器と、
（ｉ） 前記表面／法線メモリに接続されており、異方性に対するビューイング角度を受け取ると共に調節し、表面ボクセル内のインデックス及び法線ベクトルを読み出し、スクリーン位置（ｕ，ｖ）及びピクセル・カラーであるピクセル・データを形成し、所望のビューイング角度から見られる表面をレンダリングして、前記ピクセル・データをビデオ・メモリに記憶するグラフィックス・エンジンと、
（ｊ） 前記ビデオ・メモリに接続されており、前記ピクセル・データを画像として表示する表示装置とを備えたグラフィックス・システム。
【請求項６】
前記ビデオ・メモリに接続されており、前記ピクセル・データ及び前記ビューイング角度を受け取ると共に、
Ｅｓ＝（（Ａｃｏｓ（γ））2＋ｓｉｎ2（γ））1/2
（ｕ，ｖ’）＝（ｕ，Ｅｓ×ｖ）
に従って、前記画像をスケーリングするように新しい位置（ｕ，ｖ’）に対して前記ピクセル・データの位置（ｕ，ｖ）を補うストレッチング装置を更に含んでおり、ここで、Ａは、前記体積測定データのボクセルの外観比であり、γは、前記画像が異方性により発生された歪みを補正するようストレッチされるように、所望のビューイング角度のうちの仰角である請求項５に記載のグラフィックス・システム。
【発明の詳細な説明】
【０００１】
【産業上の利用分野】
本発明は、コンピュータ・グラフィック・モデルの表面を迅速に表示するシステムに関する。
【０００２】
【従来の技術】
体積測定データから表面を形成するシステムは、１９８９年４月１１日に発行されたCline、Ludke、Lorensenの米国特許第４，８２１，２１３号「固体内の２つ又はそれ以上の内部表面の同時表示のためのシステム（マーチング・キューブ法）（"System For The Simultanuous Display Of Two Or More Internal Surfaces Within A Solid Object"（'Marching Cubes Method'））」に記載されている。このシステムは、表面を効率的に発生するが、体積測定データ点の数が多くなるので、１９８８年１月１２日に発行されたCline、Ludke、Lorensenの米国特許第４，７１９，５８５号「剛体の内部領域内に含まれた表面構造の表示のための分割キューブ・システム及び方法（分割キューブ法）（"Dividing Cubes System And Method For The Display Of Surface Structures Contained Within The Interior Region Of A Solid Body"（'Dividing Cubes Method'））」に記載された他の方法に対してより適するようになっている。「分割キューブ法」は、３変数線形（trilinear）補間を用いて各々のボクセルを点及び垂直線に小分割する。これらの点及び垂直線は、レンダリング（描写）された表面を表している。
【０００３】
３変数線形補間は、特に補間を素早く実行することのできないより能力の劣ったデスクトップ・コンピュータでは、計算上の負担となる。
従って、体積データの表面を迅速に表示する単純なコンピュータ・グラフィック・システムであって、計算上の負担がより少ないコンピュータ・グラフィック・システムに対する必要性が存在する。
【０００４】
【発明の要約】
グラフィックス・システムが、複数のキューブ（立方体）から成っているキューブ／ベクトル・モデルのような表面モデルを表示する。中心ボクセル（体積要素）は、試験されるために選択されている。中心ボクセルは、提供されてもよく又は予め決定されていてもよい閾値に対して試験される。中心ボクセルが閾値を下回るならば、そのボクセルは物体の外部であると考えられ、新たな中心ボクセルが選択される。中心ボクセルが閾値を上回るならば、各々が中心ボクセルの頂点を共有している８つの対角線ボクセルが試験されて、それらのボクセルが閾値を下回るかを決定する。この比較は、速度を増すために並列して行われてもよい。対角線ボクセル値のうちの少なくとも１つが閾値を下回り、且つ中心ボクセルが閾値を上回るならば、中心ボクセルは、表面ボクセルとして分類される。
【０００５】
上の記載は、所望の物体が骨のような高いピクセル値を有していることを仮定している。本発明は又、肺のような低いピクセル値を有する構造に対しても作用する。選択された閾値は、より低い値を有しており、表面ボクセルは、少なくとも１つの対角線ボクセルが閾値を上回っている状態で、閾値を下回る中心ボクセルとして定義されている。
【０００６】
閾範囲が用いられてもよい。この場合、対角線ボクセルがその範囲の外側にあるかを決定するために試験されて、中心ボクセルはその閾範囲内にあるかを決定するために試験される。これは、イメージされるべき所望の物体が小さな範囲で変動するピクセル値を有している場合に用いられる。
次に、上方、下方、左方、右方、前方及び後方に隣接したボクセルの値が、ベクトルを決定するために用いられる。このベクトルは、ボクセルの異方性質が存在するならば、その異方性質に対して補正される。プロセスの並行性によって再び、比較が「ｘ」、「ｙ」及び「ｚ」次元において同時に実行される。
【０００７】
このプロセスは、複数の異なる中心ボクセルに対して繰り返されて、表面法線（surface normals）を有する複数の表面ボクセルが得られる。
所望の表面ボクセルを見る角度が、典型的にはオペレータによって与えられるが、システムによって予め決定されても又は選択されてもよい。このビューイング角度は次いで、有効角度を形成するために異方性に対して調整される。異方性表面ボクセルは次いで、有効ビューイング角度から見ているものとして表示される。
【０００８】
用いられた最初のボクセル・データが異方性であるならば、２次元画像は、スクリーン座標（ｕ，ｖ）に歪みのない画像を発生するために、非一様にスケーリングされる。スケーリング・ファクタＥs が計算され、スクリーン画像の「ｖ」座標の各々によって掛けられて、歪みのない画像の新たなスクリーン座標（ｕ，ｖ’）が得られる。
【０００９】
これによって、以前から可能であった相当能力の劣ったグラフィック装置に３次元表面画像を迅速に発生させることのできるシステムが得られる。
【００１０】
【本発明の目的】
本発明の目的は、体積測定データの表面を迅速に表示するコンピュータ・グラフィック・システムを提供することにある。
本発明のもう１つの目的は、相当能力の劣った設備を用いた従来のシステムと同様の速度で体積測定データの表面を表示する表面レンダリングの方法を提供することにある。
【００１１】
本発明の新規な特徴は、特許請求の範囲に具体的に記載されているが、本発明の構成及び要旨は、図面に関連した以下の詳細な記載から本発明の他の目的及び特徴と共に理解されるであろう。
【００１２】
【実施例】
デスクトップ・パーソナル・コンピュータ（「ＰＣ」）のようなロー・エンドのグラフィックス・システムには、３次元多角形モデルの表面をレンダリングするために用いられる３次元グラフィックス・アクセラレータ・カードが存在しなくてもよい。グラフィックス・カードなしにＰＣで動作する上で参照した「分割キューブ法」は、リアルタイム画像を発生するのに要するよりもはるかに低い速度でモデルを描写する。このことは、医師が対話的なリアルタイムのフィードバックを必要とするかもしれない臨床環境において用いるためにはシステムを非現実的なものにする可能性がある。
【００１３】
他の問題は、「分割キューブ法」がその方法に関連して用いられる表示リストを記憶するための大きなメモリを必要とすることである。このことは、メモリの要求が大き過ぎてＰＣで動作できなくする可能性がある。
サンプルされたデータの３次元モデルが補間なしに描写され得るならば、効率が大きく増大する。しかしながら、計算機式軸方向断層撮影法（ＣＡＴ）からのデータ・セットのような取得される多くのデータ・セットは異方性を有する。データは、磁気共鳴（ＭＲ）、超音波又は他のイメージング装置のような他の手段を用いて取得されるものであってもよい。
【００１４】
等方性体積データ１２及び異方性データ・セット１１が図１に示されている。体積ピクセル（「ボクセル」）を有している体積データ・セットの各々は、その中心で測定された物理的性質を示すデータ値を有している。１つの次元におけるデータの「スライス」の間の間隔Ｔは、体積データ・セット１１の他の次元におけるデータ点の間の間隔Ｓに等しくなく、これが異方性を生じさせている。イメージング装置は又、隣接したデータ点の間及びデータ・セット内のデータのスライスの間の間隔を記憶してもよい。補間なしに異方性ボクセルの表面を描写しようとすれば、結果として得られる画像は歪むであろう。
【００１５】
データ・セット１１に対する座標系は、「ｚ」軸に沿って並んでいるボクセルの長い軸で示されている。
体積データ・セット１２は、すべての次元において等しい間隔Ｓを有している等方性を成し、各々のボクセルは、立方体形状を有している。
表面ボクセルの決定
体積データ・セット内の表面がイメージングされることが望ましい。体積測定データ・セットのボクセルのサブセットが図２に示されている。ボクセルの層の各々の図が概略形式で図３に示されている。
【００１６】
本発明は、ボクセルを中心ボクセル５１として選択することにより作用する。中心ボクセル５１は次いで、オペレータによって提供されてもよいし又は予め決定されていてもよい閾値に対して試験される。中心ボクセル５１が閾値を下回っていれば、その中心ボクセルは物体の外部であると考えられ、新たな中心ボクセルが選択される。中心ボクセルが閾値を上回るならば、各々が中心ボクセルの頂点を共有している対角線ボクセルが試験されて、それらのボクセルが閾値を下回るかを決定する。
【００１７】
対角線ボクセル５３、５５、５７、５９、６３、６５、６７及び６９の各々は、中心ボクセル５１の単一のコーナ（角）を共有している。少なくとも１つの対角線ボクセルが閾値を上回るならば、中心ボクセル５１は、表面ボクセルとして分類される。これらの比較は、速度を増すために並列して行われてもよい。
上の記載は、所望の物体が骨のような高いピクセル値を有していることを仮定している。本発明は又、肺のような低いピクセル値を有する所望の構造に対しても作用する。この場合には、より低い閾値が選択される。表面ボクセルは、少なくとも１つの対角線ボクセルが閾値を上回っている状態で、閾値を下回る値を有している中心ボクセルとして定義されている。
【００１８】
閾範囲が用いられてもよい。この場合、対角線ボクセルの少なくとも１つがその範囲の外側にあるかを決定するために試験されて、中心ボクセルはその閾範囲内にあるかを決定するために試験される。これは、イメージされるべき所望の物体が小さな範囲で変動するピクセル値を有している場合に用いられる。
中心ボクセルは他のボクセル、典型的には隣り合っているボクセルに増大（インクレメント）させられる。次に、新たな３×３×３のグループのボクセルが選択されて、このプロセスが繰り返される。これは、すべてのボクセルが試験されるまで継続する。
【００１９】
もう１つの実施例では、中心ボクセルの後に、中間層ボクセル５２、５６、５４及び５８のような追加のボクセルが試験されてもよい。これは、速度とのトレードオフとして追加的な正確さをもたらす。
歪んだ画像の表示
次元（Ｓ，Ｔ）を有している異方性ボクセルの代わりに次元（Ｓ，Ｓ）を有している等方性ボクセルを用いると共にその差分を補正しようとすると、補間に対する必要がないので、イメージングはずっと早くなる。
【００２０】
典型的には、表面ボクセルが見えるビューイング角度（仰角、方位角、ロール（roll）角）をオペレータが与える。ビューイング角度は又、予め決定されていてもよいし又は予め設定された範囲にわたって走査されてもよい。ビューイング角度は、「ｘ」、「ｙ」及び「ｚ」軸の周りの回転として、その順でそれぞれ定義される。
【００２１】
イメージングでは、典型的には方位角及び仰角のみが用いられる。例えば、向かい合って直立している被検体に関して、その被検体の頭部の画像を撮る場合を考える。「ｘ」軸は被検体の耳を通過しており、「ｙ」軸は頭の頂部を通って首の中心を通っている。「ｚ」軸は鼻の頭部を通過して頭蓋骨の後ろから出ている。「ｘ」軸の周りの回転は、「はい」と合図しているように頭を上下に回転させる。「ｙ」軸の周りの回転は、「いいえ」と合図しているように頭を左右に回転させる。「ｚ」軸の周りの回転は単に、画像を逆さまに回転させる。従って、３つの方向角のうちの最後の角、ロール角は、時計方向回り又は反時計方向回りの様式でスクリーン画像を回転させるのと等価であるので、画像に殆ど影響を及ぼさない。イメージングの目的のためには、方位角及び仰角のみが検討される。
【００２２】
先ず実行されるのは方位角の回転であり、体積データ・セットは、方位角回転の軸が異方性ボクセルの長い軸と平行になるように整列させられ、その結果は、等方性ボクセルを用いたのと同じになる。
従って、補正するために残された唯一のビューイング角度は、仰角である。１９９３年７月６日に発行されたCline、Ludke、Dumoulin、Souzaの米国特許第５，２２６，１１３号「反対回転の放射線鋳物を用いた体積測定投影レンダリングのための方法及び装置（"Method and Apparatus for Volumetric Projection Rendering Using Reverse Ray Casting"）」は、本願発明の譲受人に譲渡され、ここに参照されるべきものであるが、本特許には、結果として得られる画像を正確な投影が与えられるようにストレッチする（引き伸ばす）ことにより画像が形成された後に、投影画像における異方性が補正された。
【００２３】
図４は、仰角に対する異方性の効果を示す側面図において異方性体積データ・ボクセル１１及び等方性体積データ・ボクセル１２を示している。表面ボクセルを仰角γで見ようとすると、等方性データを用いて実際には角度Ψになる。γはΨに次式によって変換される。
Ψ＝ａｒｃｔａｎ（ｔａｎ（γ）／Ａ）
ここで、Ａ＝Ｔ／Ｓ
表面法線
各々の表面ボクセルについて、表面に対する法線ベクトルが形成される。図３における中心ボクセル５１として示される各々の表面ボクセルに対して、６つの隣り合うボクセルが、表面法線ベクトルを決定するために用いられる。中心ボクセル５１の真下の下方ボクセル７３が、「ｙ」方向における差分を決定するために上方ボクセル７１から減算された値を有している。同様に、左方ボクセル７５が、「ｘ」方向における結果を得るために右方ボクセル７７から減算される。又、前方ボクセル８１の値が、「ｚ」方向における結果を得るために後方ボクセル８３から減算される。「ｘ」、「ｙ」及び「ｚ」の差分は、次式に従って、大きさｇを決定するために用いられる。
【００２４】
ｇ＝｛（「ｘ」の差分）2 ＋（「ｙ」の差分）2
＋（「ｚ」の差分）2 ｝1/2
ここで、（「ｘ」の差分）＝（ボクセル７７）−（ボクセル７５）
（「ｙ」の差分）＝（ボクセル７１）−（ボクセル７３）
（「ｚ」の差分）＝（ボクセル８３）−（ボクセル８１）
各々は次式に従って、規格化されている。
【００２５】
Ｎx ＝（「ｘ」の差分）／ｇ
Ｎy ＝（「ｙ」の差分）／ｇ
Ｎz ＝（「ｚ」の差分）／（Ａ×ｇ）
Ａ＝Ｔ／Ｓ
ここで、Ａは、ボクセルの長い辺Ｔのより短いか等しい辺Ｓに対する比である外観比（aspect Ratio）である。
【００２６】
同時に、（Ｎx ，Ｎy ，Ｎz ）は、中心ボクセル５１に位置している表面に垂直なベクトル８５を定義している。これは、中心ボクセル５１のような表面ボクセルの各々に関して繰り返される。
画像のストレッチ
ビューイング角度Ψに対する補正を要求することの他に、異方性ボクセルを用いることは又、スクリーン高さＨ’及びＨに対する補正を要求する。
【００２７】
発生された画像の座標に掛け算して１つのスクリーン方向における画像を補うために、スケーリング・ファクタＥs が用いられる。表面ボクセルは次いで、方位角Ψから見られるように表示されて、その結果、各々がスクリーン位置（ｕ，ｖ）を有しているピクセルが得られる。座標系が、異方性ボクセルの長い軸に垂直な軸から測定される仰角に対して設定されているならば、垂直スクリーン座標、「ｖ」座標は、下記に定義されるように、スケーリング・ファクタＥs によって掛け算される。
【００２８】
Ｅs ＝（（Ａｃｏｓ（γ））2 ＋ｓｉｎ2 （γ））1/2
（ｕ，ｖ’）＝（ｕ，ｖ×Ｅs ）
図５は、本発明の簡略されたブロック図を示している。イメージング装置１が体積測定データ・セットを形成している。上述した装置の他に、体積測定データを発生するＣＡＤ、ＲＡＤＡＲ又は他の測定装置を採用してもよい。
【００２９】
３次元体積測定データ・セットは、より大きな共有メモリ３０の一部であってもよい体積メモリ３１に記憶される。この体積メモリは又、外観比Ａ、即ち異方性ボクセルの寸法を記憶してもよい。
論理装置が、試験されるべき中心ボクセルを決定する。この論理装置は、中心ボクセル５１の位置（インデックス）を記憶し、又は中心ボクセル５１の位置を比較器４１に対して直接に通過させる。比較器４１は又、等価表面（isosurface）を定義している閾値を受け取る。この閾値は、オペレータによって手動で供給されてもよいし、予め決定されていてもよいし、又は他の装置から計算されて受け取られてもよい。上述したように、閾範囲は、試験されるボクセルが閾範囲内にあるか又は閾範囲外にあるかを決定するために用いられ得る。比較器４１は次いで、体積メモリ３１から中心ボクセル値を引き出す。この比較器は次いで、中心ボクセル値を閾値に対して試験する。中心ボクセルが閾値を下回れば、比較器はこのことを論理装置４９に指示し、論理装置は他の中心ボクセル・インデックスを選択して、このプロセスは、（所望の物体が上述したような高いボクセル値を有していると仮定して）繰り返される。
【００３０】
中心ボクセルの値が閾値を上回れば、複数の比較器４３、４５及び４７の各々が、閾値、及び体積メモリ３１からの中心ボクセルのインデックスを受け取ると共に、図２及び図３の参照番号５３、５５、５７、５９、６３、６５、６７及び６９のような対角線ボクセルを体積メモリ３１から引き出す。対角線ボクセルはすべて、中心ボクセル５１に隣り合うボクセルである。各々の比較器は、対角線ボクセルの値を、受け取った閾値と比較する。
【００３１】
比較の結果は、論理装置４９に通過させられて、論理装置は、表面ボクセルを示している、閾値を下回る少なくとも１つの対角線ボクセルの値が存在するかを決定する。
論理装置４９は、表面ボクセルの位置（又はデータ・セット・インデックス）を共有メモリ３０の一部であってもよい表面／法線メモリ３３に記憶する。表面ボクセルに遭遇するごとに、論理装置は、表面ボクセルの位置（又はインデックス）を直接に又は共有メモリを介して「ｘ」、「ｙ」及び「ｚ」減算装置２１、２３及び２５にそれぞれ示す。
【００３２】
「ｘ」減算装置２１は、左方ボクセル７５を得るために、表面ボクセル・インデックスの「ｘ」インデックスを減少させる。「ｘ」減算装置は次いで、右方ボクセル７７の値を得るために、表面ボクセル・インデックスの「ｘ」インデックスを増大させる。「ｘ」減算装置は次いで、右方ボクセルから左方ボクセルを減算して、「ｘ」の差分を決定する。
【００３３】
同様に、「ｙ」減算装置２３は、下方及び上方ボクセル７３及び７１を得るために、表面ボクセルの「ｙ」インデックスを減少及び増大させて、「ｙ」の差分を決定する。
最後に、「ｚ」減算装置２５は、前方及び後方ボクセル８１及び８３を得るために、表面ボクセルの「ｚ」インデックスを減少及び増大させて、「ｚ」の差分を決定する。
【００３４】
規格化器（ノーマライザ）２７が「ｘ」、「ｙ」及び「ｚ」の差分を受け取り、次式に従って、大きさｇを算出する。
ｇ＝｛（「ｘ」の差分）2 ＋（「ｙ」の差分）2
＋（「ｚ」の差分）2 ｝1/2
規格化器２７は、外観比Ａを体積メモリ３１から読み取り、次式に従って、ベクトル成分Ｎx 、Ｎy 及びＮz を計算する。
【００３５】
Ｎx ＝（「ｘ」の差分）／ｇ
Ｎy ＝（「ｙ」の差分）／ｇ
Ｎz ＝（「ｚ」の差分）／（Ａ×ｇ）
Ｎx 、Ｎy 及びＮz は、中心ボクセル５１の位置における表面に垂直なベクトル８５を定義しており、それが関連している表面ボクセルの位置（インデックス）と共に表面／法線メモリ３３に記憶される。
【００３６】
これは、レンダリングしようとするすべてのボクセルが試験されるまで繰り返される。グラフィックス・エンジン５０が、表面／法線メモリ３３の記憶内容を読み取ると共に、オペレータから手動で入力されてもよいし、予め決定されていてもよいし又は他の装置によって供給されてもよい所望のビューイング角度（方位角、仰角（γ）、ロール角）を受け取る。グラフィックス・エンジン５０は、所望の仰角γを上述したような異方性データに対する有効仰角Ψに変換し、次いで、方位角、有効仰角Ψから見えるように所望の表面をレンダリングする。グラフィックス・エンジン５０は、その出力を２次元スクリーン・ピクセル・データとしてビデオ・メモリ３５に記憶する。等方性データの場合には、２次元スクリーン・ピクセル・データはその後、表示装置６０に表示されてもよい。
【００３７】
異方性データ・セットが用いられるときに、ストレッチング装置５５がビデオ・メモリ３５に結合され、２次元スクリーン・ピクセル・データを読み取り、仰ビューイング角度を受け取って、スクリーン・ピクセル位置の各々の垂直位置「ｖ」に次式に従ってスケーリング・ファクタＥs を掛ける。
Ｅs ＝（（Ａｃｏｓ（γ））2 ＋ｓｉｎ2 （γ））1/2
（ｕ，ｖ’）＝（ｕ，Ｅs ×ｖ）
これは、異方性によって招かれる歪みとは反対の方向に画像を事実上歪ませるボクセルの間隔及びビューイング角度の非一様性の相対的割合に従って２次元スクリーン・ピクセル・データのデータ値の位置を補う（オフセットする）。
【００３８】
ビデオ・メモリ３５は、正確にレンダリングされた表面を示している表示装置６０に表示されてもよい。
図６は、被検体の頭部の計算機式断層撮影（ＣＴ）体積測定データのシングル・スライスである。これは、体積測定データ・セットにおける９３スライスのうちの１つのスライスであり、このスライスは、１．５ｍｍに等しいピクセル寸法Ｓを有している。
【００３９】
図６の画像を形成するために用いられた体積測定データ・セットは、ＰＣを採用して本発明に従って処理された。図７の３次元画像は、２２４のハンスフィールド（Hounsfield）単位の閾値を用いて本発明により発生された。
新規な本発明についての現在好ましい様々な実施例をここに詳細に記述してきたが、当業者には多くの改変及び変更が明らかになるであろう。従って、本発明は、例示の方法によって明らかにされた特定な詳細によってではなく、特許請求の範囲の要旨によってのみ限定されるものとする。
【図面の簡単な説明】
【図１】等方性及び異方性の体積測定データから成るボクセルを例示する図である。
【図２】本発明に従って表面ボクセルを決定するために用いられる体積測定データのボクセルを同定することを例示する図である。
【図３】本発明に従って表面ボクセルを決定するために用いられる体積測定データのボクセルを同定することを例示する図である。
【図４】仰角の異方性の効果を例示する図である。
【図５】本発明によるコンピュータ・グラフィック・モデル化システムの簡略されたブロック図である。
【図６】被検体の頭部から取得された計算機式断層撮影体積測定データ・セットのシングル・スライスの画像である。
【図７】図６を形成するのに用いられた体積測定データ・セットの骨表面の３次元画像である。
【符号の説明】
１  イメージング装置
１１  異方性体積データ・ボクセル
１２  等方性体積データ・ボクセル
２１  「ｘ」減算装置
２３  「ｙ」減算装置
２５  「ｚ」減算装置
２７  規格化器
３０  共有メモリ
３１  体積メモリ
３３  表面／法線メモリ
３５  ビデオ・メモリ
４１、４３、４５、４７  比較器
４９  論理装置
５０  グラフィックス・エンジン
５１  中心ボクセル
５２、５６、５４、５８  中間層ボクセル
５３、５５、５７、５９、６３、６５、６７、６９  対角線ボクセル
５５  ストレッチング装置
６０  表示装置
７１  上方ボクセル
７３  下方ボクセル
７５  左方ボクセル
７７  右方ボクセル
８１  前方ボクセル
８３  後方ボクセル
８５  垂直なベクトル
【図１】
0004063380.tif 000002
【図２】
0004063380.tif 000003
【図３】
0004063380.tif 000004
【図４】
0004063380.tif 000005
【図５】
0004063380.tif 000006
【図６】
0004063380.tif 000007
【図７】
0004063380.tif 000008