JP4435329B2 - 物体レンダリングシステム及び方法 - Google Patents
物体レンダリングシステム及び方法 Download PDFInfo
- Publication number
- JP4435329B2 JP4435329B2 JP13836199A JP13836199A JP4435329B2 JP 4435329 B2 JP4435329 B2 JP 4435329B2 JP 13836199 A JP13836199 A JP 13836199A JP 13836199 A JP13836199 A JP 13836199A JP 4435329 B2 JP4435329 B2 JP 4435329B2
- Authority
- JP
- Japan
- Prior art keywords
- simulated
- image
- ray beam
- ray
- shell
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
- 238000009877 rendering Methods 0.000 title claims abstract description 45
- 238000000034 method Methods 0.000 title claims description 21
- 239000007787 solid Substances 0.000 claims abstract description 69
- 239000000463 material Substances 0.000 claims abstract description 16
- 239000013598 vector Substances 0.000 claims description 26
- 230000002238 attenuated effect Effects 0.000 claims description 6
- 230000004044 response Effects 0.000 claims description 5
- 239000011800 void material Substances 0.000 claims description 5
- 238000002591 computed tomography Methods 0.000 claims description 2
- 238000009659 non-destructive testing Methods 0.000 abstract description 4
- 230000008569 process Effects 0.000 description 10
- 230000006870 function Effects 0.000 description 9
- 230000008901 benefit Effects 0.000 description 4
- 238000004364 calculation method Methods 0.000 description 3
- 230000015572 biosynthetic process Effects 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 239000002131 composite material Substances 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000003384 imaging method Methods 0.000 description 2
- 238000003908 quality control method Methods 0.000 description 2
- 230000005855 radiation Effects 0.000 description 2
- 238000010561 standard procedure Methods 0.000 description 2
- 238000010521 absorption reaction Methods 0.000 description 1
- 238000010009 beating Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000000872 buffer Substances 0.000 description 1
- 238000013016 damping Methods 0.000 description 1
- 230000002950 deficient Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 238000012804 iterative process Methods 0.000 description 1
- 210000004165 myocardium Anatomy 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T15/00—3D [Three Dimensional] image rendering
Landscapes
- Engineering & Computer Science (AREA)
- Computer Graphics (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Image Generation (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Processing Or Creating Images (AREA)
- Analysing Materials By The Use Of Radiation (AREA)
Description
【発明の属する技術分野】
この発明は、物体レンダリングシステム及び方法に関し、より詳細には合成X線イメージをレンダリングする非ボクセルベースのシステム及び方法に関する。
【0002】
【従来の技術】
合成X線イメージをレンダリングする標準的な方法では、ボクセルをベースとしたレンダリングを用いる。この方法論において、物体はそれぞれボクセルと呼ばれる小さな立体要素からなる。ボクセルは、光線追跡されてレンダリングされたイメージを生成する。
【0003】
【発明が解決しようとする課題】
しかしながら、このプロセスには2つの制限がある。
第一に、物体に近似させるために必要なボクセルの数が非常に多くなる可能性がある。これにより、そのような大量データはある用途に十分速く処理できない場合もあるという点において、大部分のコンピュータ・グラフィックス・システムに課題を課す。
【0004】
第二に、ボクセルは形を変更することができないという点において、フレキシブルでない物体を画定する。物体の全体的な形状が変化した場合、「ボクセレーション(voxelation)」すなわちボクセルを新たに集成して新たな全体形状を描写することは、計算に費用がかかり、アルゴリズムが複雑である場合がある。
【0005】
この発明は上述した点に鑑みてなされたもので、物体をレンダリングしてX線イメージの形成をシミュレートするより単純でより非計算集約的な物体レンダリングシステム及び方法を提供するものである。
【0006】
【課題を解決するための手段】
この発明に係る物体レンダリングシステムは、物体の材質及び幾何記述を模擬X線イメージとしてレンダリングする物体レンダリングシステムであって、上記物体を、表面及び所定の質量減衰係数を有する、入れ子構造をした複数の立体の3Dモデルに分解する手段と、上記3Dモデルを、模擬X線ビームのソース及び上記複数の立体に関連するシーングラフを含む3D座標系に登録する手段と、上記複数の立体を通して上記模擬X線ビームを通過させる手段と、上記模擬X線ビームが上記複数の立体の特定の1つに入射するかまたは出射するかを判断する手段と、これに応答して上記模擬X線ビームの経路に沿った空所を位置決めする手段と、上記材質及び幾何記述と上記空所に従って上記模擬X線ビームを減衰させる手段と、上記模擬X線ビームの減衰に応答して、上記模擬X線イメージとして上記複数の立体を独立してレンダリングする手段とを含むものである。
また、この発明に係る物体レンダリング方法は、物体の材質及び幾何記述を上記物体の模擬X線イメージにレンダリングする物体レンダリング方法であって、上記物体の入れ子構造をした複数の立体の3Dモデルを提供するステップと、上記立体をデジタル化して、上記物体のそれぞれに対応する表面を提供するステップと、質量減衰係数を、それぞれ対応する表面により境界づけられたそれぞれの立体に割り当てると共に、対応する質量減衰特性を上記対応する表面に割り当てるステップであって、上記表面はシェルを形成し、各シェルは対応する所定の質量減衰特性を有する、ステップと、上記シェルをシーングラフに位置登録して、上記シェルの表面ジオメトリを関連させるための登録されたシェルジオメトリを提供するステップと、上記質量減衰係数及び上記質量減衰特性に従って模擬X線ビームを減衰させながら、上記物体の上記材質及び幾何記述を通して上記模擬X線ビームを通過させるステップと、上記シェルのそれぞれの深さイメージ画像を提供し、それにより、上記対応する立体を介して通過する上記模擬X線ビームの所定の光線によって横断される距離を示すように、上記シェルをレンダリングすることにより、上記模擬X線ビームが通過する上記シェルにより境界づけられたそれぞれの立体の経路長を得るステップと、上記対応する質量減衰係数を使用して各シェルに対する減衰イメージを計算することによって、上記深さシェルイメージのそれぞれを関連する質量減衰係数と組み合わせ、それによってその密度を示すために影の付けられたシェル減衰イメージを提供するステップと、
上記影の付けられた減衰イメージを組み合わせることにより、そこで表現された立体を有する組み合わせられた模擬X線イメージを提供するステップとを含むものである。
【0018】
【発明の実施の形態】
この発明おけるシステムは、物体をX線が通過するようにシミュレートし、減衰をシミュレートし、そして例えば境界内の境界、表面内の表面、あるいはシェル内のシェルといった、物質内に埋め込まれた物質の一つとなる物体を考慮することによって、非ボクセル物体データをX線フォーマットで表示可能なデータに変換する。
【0019】
レンダリングを達成するために、物体は、一方を他方の内側でネストしたと考えられる多数の立体に分解され得ることがわかってきた。これらの立体は、一般に、その立体中および外を通過する光線を利用することにより検出可能な表面によって境界が定められる。これらの表面は、フォトンビームが通過する際、立体によって減衰されるフォトンビームのシミュレートされた通過について位置決めすることができる。
【0020】
したがって、立体の表面は、一実施の形態において所与のフォトン群の方向に対応する方向および減衰していないフォトン数に対応する大きさを有するグラフィックス空間においてベクトルを利用する従来のコンピュータグラフィックスの描写によって得ることができる。
【0021】
したがって、本システムは立体を取り上げ、それをネストされた多数の立体に分解する。このネスティングの目的は、上記のネスティングによって規定される内部位置での立体の所定の特徴を得ることを可能にすることである。ネストされた表面あるいは境界を決定することは、ボクセルによって立体の内部要素それぞれを規定するよりもはるかに簡略である。
【0022】
立体の表面は、標準的な3Dビジュアルモデル化パッケージによってモデル化されて、ネストされた境界を作成する。各立体内部には、その実世界値(real world value)に基づいた質量減衰係数(mass-attenuation coefficient)およびその境界を取り巻く立体の質量減衰係数が割り当てられる。次いで、これらのモデルは、フォトンエミッタおよび結像面(imaging surface)を含む3D座標系に、シーン中の全ての物体に関連するシーングラフによって記録される。
【0023】
その後、シーングラフ中のデータは、X線が物体を通過するときに得られるイメージに近似させるようにレンダリングされる。より詳細には、本システムは、物体の3Dモデル、物体を通過するシミュレートされたフォトン放射、およびそのネストされた立体の全てを使用して、そしてイメージを出力するシミュレートされたX線カメラによって検出されるX線状イメージを提供する。
【0024】
一実施の形態において、3Dモデルは、モデル化される立体についての専門家の知識に基づいて手動で行われる。たとえば、外科医は心筋層の組織密度および各種心室の境界を特定することによって、心臓の内層を特定したいと望む場合がある。あるいは、その3DモデルのデータセットをCTスキャンから得て、その結果としてボクセル化されたデータを得ることもできる。これは、類似し近接するボクセルの多くを一つの立体に変えることができると仮定した場合、有益である。その結果、ネストされた層または表面すべての立体イメージを生成するレンダリングプロセスを、かなり少ないデータおよび同様に低減されたプロセスで達成することができる。
【0025】
上記により、新たなセットのボクセルを計算するのに対して、表面形状を変えることによって立体/ボリュームのリアルタイムアニメーションが可能となる。その結果、本システムにより、X線イメージを見慣れている医療実習生が容易に理解するフォーマットで拍動する心臓をシミュレートすることが可能となる。
【0026】
本技術の別の主要な用途は、非破壊的な試験の分野にある。ここでは、欠陥のある機械加工パーツをモデル化し、その結果得られる3Dモデルを利用して、評価すべきこれらのイメージをX線フォーマットで品質管理の専門家に提示して教えるトレーニングシステムを工夫することができる。この用途において、品質管理の専門家がX線検査システムを利用することが普通であり、したがって、本システムにより価値のあるトレーニングツールが提供される。
【0027】
より詳細には、本システムは、立体の物体を通してフォトンビームが減衰するという興味ある数学的特性を利用するアルゴリズムを採用している。フォトンビームの減衰プロセスを考慮して、各辺の長さがSである正方形の断面を有する長さLの立体の棒が与えられている場合、フォトンビームがその棒の一面に投射されると何が起こるのかを考察することができる。簡単に言えば、一または大部分のフォトンは、散乱するかあるいは棒によって吸収される。この吸収を説明する数学的な関係は以下の数式によって与えられる。
【0028】
I(E)=IO(E)exp(−mu(E)×L) (1)
ここで、
IO(E)=フォトンエネルギー(特定期間中の棒の表面に衝突する特定エネルギーでのフォトン数)の関数としての入射フォトンストリームの強度
I(E)=フォトンエネルギーの関数としての、反対の端において棒から離れたフォトン数
mu(E)=フォトンエネルギーの関数としてのこの物体についての減衰係数
exp( )=オイラーの定数を累乗する関数
L=棒の長さ
E=入射ビーム中のフォトンあたりのエネルギー
【0029】
特定のマテリアルが別のマテリアル内に埋め込まれた場合に何が起こるのかを考察することができる。この減衰は以下のように計算することができる。
【0030】
I=IO(E)×exp(−mu1(E)×(F2−F1))
×exp(−mu2(E)×exp(−mu1(E)×(B1−B2))
=IO(E)×exp(−mu2(E)−mu1(E))×(B2−F2))×exp(−mu1(E)×(B1−F1))
【0031】
レンダリングアルゴリズムは、完全な物体をレンダリングすることは、複数の単純物体をレンダリングすることに相当するという事実を利用している。一般に、n個の物体が互いにおいてネストされる場合、これらの物体をn個の離散物体に分解することができ、物体それぞれの減衰係数はその当初の減衰係数とそれを取り囲む物体の当初の係数との間の差である。
【0032】
特に、レンダリングアルゴリズムは以下のように実行することができる。
初期のX線レンダリングアルゴリズムでは、正確な結果を生成するために、「ネストされた」凸状物体を使用する必要があった。その後、技術が発見され、それによって、より一般的な性質、すなわち凸状または非凸状の「ネストされた」物体を用いることが可能となった。以下説明することは、初期のレンダリングアルゴリズムにおける変化を再現してこの新しいレベルの一般原理を理解するよう意図されている。
【0033】
オリジナルアルゴリズムは以下のように説明される。
I.X線イメージを定数値IOに設定する。
II.レンダリングされる全ての物体について、
A.最大範囲イメージを小さな数に設定する。
B.最小範囲イメージを大きな数に設定する。
C.最大および最小範囲値を最大および最小範囲イメージそれぞれに保存し、物体をレンダリングする。
D.最大と最小範囲イメージとの間の絶対差を取ることによって、各ピクセルの経路の長さを計算する。
E.式(1)に記載されている減衰関数を使用して減衰イメージを計算する。
F.X線イメージに減衰イメージで倍加することによって減衰イメージをX線イメージに合成し、その結果をX線イメージに記憶する。
なお、範囲イメージは、シミュレートされたフォトン源から対応する立体表面までの範囲を知って生成されたイメージを指す。
【0034】
新しいアルゴリズムはステップII.A.〜II.D.に代わり、その結果以下のようになる。
I.X線イメージを定数値IOに設定する。
II.レンダリングされる全ての物体について、
A.瞥見ベクトル(eye beam vector)および物体上の表面のパッチをこのベクトルに沿って定義する。
B.その物体上の表面のパッチすべてについて、
1.その表面のパッチに対応するすべてのイメージピクセルについて、a.表面法線ベクトルとともに視線ベクトル(eye vector)のドット積を計算する。
b.ドット積が0ではない場合、
1.視点から表面ポイントまでの距離を計算する。
2.上記距離の符号をドット積の反対符号に設定する。
3.上記距離を範囲イメージに加える。
注: この時点で、範囲イメージには、その物体が影響を及ぼすイメージの各ピクセルについて物体までの経路の長さが含まれる。
C.前ステップで計算したイメージとともに、数式1に記載した減衰関数を使用して減衰イメージを計算する。
D.X線イメージに減衰イメージを倍加することによって減衰イメージをX線イメージ中に合成し、その結果をX線イメージに記憶する。
【0035】
非凸状物体を用いることに加えて、このアルゴリズムにより最小および最大範囲バッファの必要性がなくなる。
【0036】
本システムの主要な利点は、レンダリング可能な物体を表示するために、単に固定ボクセルだけではなく、任意の形状を使用できることである。これによって、物体を表示するために必要なデータ量が減少し、高データ帯域幅の問題がなくなる。さらに、物体の形状が柔軟であるため、物体の形状を変えることがしばしば境界の表示を簡単に変えられることを意味し、新しい形状についての計算の複雑さを減少させる。
【0037】
要約すると、専門化されたアルゴリズムを有する物体レンダリングシステムが、医療での用途および産業上の非破壊的試験で使用するため、物体の物質的及び幾何学的描写をその物体のX線イメージに変換するために提供され、それによって、X線の結果ではなく物体表現をX線フォーマットで行うことが可能となる。本システムは、各物体がネストされた幾何学的立体からなり、それぞれが囲み面によって画定され、ネストされた立体それぞれのマテリアルの特性はX線の質量減衰係数によって定義されることを考慮することによって機能する。囲み面それぞれの質量減衰係数を調整する手段は、各立体を個別にレンダリングすることを可能にするため、提供される。
【0038】
本システムでは、物体の形状はしばしば「境界表現」と呼ばれる立体形状で閉曲面によって描写される立体であり、そのアルゴリズムによって物体が、境界をネストすることによって互いの内側でネストすることができる。
【0039】
上記のアルゴリズムへの別の入力は、その物体を通過するシミュレートされたX線ビームの数学的な描写と、視点、焦点距離およびシミュレートされたX線ビームに対するフィルタの応答を含むカメラの形状である。
【0040】
このようなレンダラー(renderer)の利点は、レンダラーがX線結像装置の使用をシミュレートするシステムでの構成要素として機能することができるという事実である。これらのシミュレータは、実際のX線生成装置は、動作中に、装置を操作している人または付近にいる人に有害である可能性があることから、周囲を訓練させる際に望ましい。レンダリングアルゴリズムは、標準的なボクセルレンダリング技術に適さない物体の幾何学的描写を取り扱う。
【0041】
この発明のこれらおよび別の特徴は、図面と共に詳細な説明に従ってよりよく理解されるであろう。
図1および図2を参照して、ネストされた立体(nested solids)14、16および18を有する3Dモデル12が入力として与えられるシステム10を示す。これらのネストされた立体の外側は、それぞれ20、22および24においてデジタル化される。これは、ネストされた立体の内側表面、中間表面および外側表面を効果的にデジタル化する。さらに、26、28および30において、質量減衰が内側、中間および外側の立体それぞれに割り当てられる。その結果、所定の質量減衰特性を有する均一な物質(material)で満たされた、これら表面に対応するシェル(shell)32、34および36のセットが得られる。
【0042】
38において、これらのシェルは表面形状を互いに関連させるため、記録されたシェル形状40を生成するようネストされたまま記録される。これらの表面形状は、各レンダラー42、44、および46においてレンダリングされる。
【0043】
レンダラーへの他の入力は、X線ビーム描写50およびカメラ描写52である。各レンダラーは、それぞれのシェル画像を提供し、それにより内側シェルの深さイメージ54が内側シェルについて、中間シェルの深さイメージ56が中間シェルについて、および外側シェルの深さイメージ58が作成される。シェルの深さイメージを生成する目的は、所与のフォトン光線が立体を通過した距離を示すことである。
【0044】
このことは、イメージにおいて結果として得られるピクセル値はフォトン光線がその物体まで動いた部分的な距離の関数であるため、重要である。そして、深さシェルイメージは質量減衰係数と組み合わされて、60、62、および64において減衰イメージの計算が可能となる。これにより、それらの密度を示すために影が付けられた内側、中間、および外側シェルの減衰イメージ66、68および70が提供される。減衰イメージは、倍加プロセス72においてピクセル毎に組み合わされ、その結果、例示するように、シミュレートされたX線イメージ74がネストされた立体14’、16’、および18’と組み合わされる。
【0045】
その結果、内側の表面を描写するイメージが形成され、選択された表面を容易に識別でき、相対的な密度も示すように、当初の立体に一貫性を持たせる。これによって、問題となっている物体の内側の特性について訓練生を指導するための優れたトレーニングツールが提供される。
【0046】
図3を参照して、膨らんだ上部領域82、より大きく膨らんだ下部領域84、およびくびれ部分86を有する非凸状物体80を示している。図3の目的は、非凸状物体を通る正確な経路の長さを得ることが可能であることを示すことである。
【0047】
本システムは、ネストされたシェルまでの経路の長さの描写が単純であるため、ネストされた凸状シェルにしか適応できないと当初考えられていた。しかしながら、視線ベクトルおよび表面法線ベクトルのドット積の比較を利用することによって、立体の非凸状部分を適切に描写できることが判明している。
【0048】
説明のため、シミュレートされたX線ビームはいわゆる視線ベクトルに沿って物体まで動くものであるといえる。ここで、視線ベクトルは、従来のようにベクトルが目89を指すように、矢印88で示され、視線ベクトルが通過する物体の多様な表面が90、92、94および96において示され、境界面のパッチと呼ばれる。視線ベクトルの入射あるいは出射ポイントでの表面のパッチに対する法線は、98、100、102、および104において示される。
【0049】
これらの法線は続いて立体を出射する方向を示す。上記のドット積を使用する目的は、視線ベクトルに沿う立体の部分間の空間を判断できるようにすることである。この視線ベクトルに沿って隙間すなわちオープンスペースを判断するため、視線ベクトルの入射ポイントおよび出射ポイントを判断するする必要がある。
【0050】
これは、一実施の形態において、視線ベクトルのドット積を与えることにより、また表面法線により、そしてこのドット積が正であるかあるいは負であるか、またはゼロであるかどうかを判断することにより達成される。
【0051】
このドット積がゼロである場合、光線は物体をそれて、入射または出射ポイントとしては考慮されない。ドット積が負である場合、交点は入射ポイントであると考えられる。ドット積が正である場合、交点は出射ポイントであると考えられる。目あるいは光線源から交点までの距離である深さ値の符号は、ドット積の符号とは反対に設定される。
【0052】
次に、補正された深さ値は、イメージに集積される。集積が完了すると、物体を通過するフォトン光線の距離が集積されたイメージ中に含められる。
出射ポイントおよび続く入射ポイントを知ることによって、出射ポイントと次の入射ポイントとの間に隙間があることが明らかになる。
【0053】
図4及び図5を参照して、本システムのフローチャートを示す。
ここで、プロセスは、図1および図2において、40で示される記録されたシェルの形状で開始され、レンダラーに進み、イメージ減衰計算、減衰イメージの生成、上記イメージを一緒にピクセル毎に倍加し、74で示すように最終的なX線状イメージをもって終了する。
【0054】
理解されるように、図3を参照して、120で開始されると、122で示すようにイメージのビーム強度を得る必要があり、この目的は、X線放射を特徴付け、イメージ化される立体が減衰しないときを検出するためであり、これによりベースラインが生成される。したがって、ベースラインイメージの強度はIOで示される。
【0055】
これは、図1および図2のシステムの第一通過に対応し、ここでは、ピクセル毎にイメージを倍加するユニット72が、レンダリングされたイメージが合成される初期イメージを必要とする。
【0056】
その後、124で示すように、集積イメージでのピクセル値がユニット122から導出されるようにIOの値に設定される。次に、126で示すようにiで示されるレンダリングすべき立体それぞれについてのステップがブロック128、130、および132に示される。
【0057】
説明するように、ループが境界のパッチjそれぞれに実行され、ここで、「境界のパッチ」という用語は、法線ベクトルがこのポイントおよびこのパッチについて確立できるように、光線が通過する立体の境界のわずかな部分を指す。これらの法線がそのパッチに投射された各ピクセルについて確立されることは重要である。立体の境界を単に数学的に説明された小さなパッチに分解することは、コンピュータグラフィックスでの標準的な技法である。
【0058】
X線ビームの描写とカメラの描写がパッチ毎に各レンダラーに適用されるのは、この継ぎ目においてである。立体の境界はパッチの集成であるため、ブロック126に示すループの目的は、パッチのリストとして表現されるネストされた立体それぞれについて各種境界を得ることである。
【0059】
このループプロセスは、各立体について実行される。図示するように、各立体について、範囲イメージRでのピクセル値を128で示すように0に設定し、130で示すようにその物体iの境界表現を得るとともに、物体iの質量減衰係数(mu)を得なければならず、その結果、その立体の完全な特徴付けがなされる。
【0060】
その後、134で示すように、各パッチ(j)はこのループで処理されて、その立体の境界を合成するパッチに対する基準を確立する。その後、136で示すように別のループが実行され、境界のパッチに投射されるイメージピクセルが同様に基準とされる。
【0061】
この基準化はブロック138、140、142、144、146および148によって利用される。先に説明したように、その立体までの経路の長さを計算できるように、境界のパッチ上の各ポイントについて法線ベクトルを得、視線ベクトルとのドット積を得ることは重要である。この経路の長さは、光線の入射ポイントと出射ポイントとの距離を定義するため、レンダリングにおいて重要である。この距離を確立してから、隙間があるかどうかを判断するため、ドット積の符号は、入射ポイントまたは出射ポイントのいずれかを示す。
【0062】
理解されるように、まず、境界のパッチj上にあるポイントkの法線ベクトル(nv)を得る。これによりパッチでのポイントと最終的なレンダリングイメージに形成されるイメージピクセルとの間の対応(correspondence)が確立される。
【0063】
その後、140に示すように、視線ベクトルと法線ベクトルのドット積を得て、ドット積dpをもたらす。ドット積が0に等しい場合、判断ブロック142に示すように、プロセスは150で示すように継続し、この境界のパッチでの別のピクセルを処理するためブロック136にループバックする。
【0064】
ドット積が0ではない場合、144で示すように、本システムは視点から境界のパッチjのポイントkまでの間の距離を計算する。この距離はブロック146に与えられ、このブロックにおいて、距離dの符号がドット積の符号とは反対に設定される。距離の符号をドット積の符号の反対に設定した上で、その距離を集積イメージに加え、それにより位置kでのピクセル値Rがdだけ増分される。増分した結果、立体の物体を通る入射光線の経路の長さが得られる。
【0065】
すべての境界のパッチが処理されると、152に示すように、ループにより結果として得られたイメージが合成される。これは、ピクセル位置mでのイメージRからピクセル値を得ることによってなされる。この値は、154に示すようにxと呼ばれる。
【0066】
図1および図2のブロック60、62および64で示すように、減衰計算は、出力として光線の減衰率を与えるため、上記の式(1)に従って156に示すように実行される。言い換えると、経路の長さと質量減衰値が得た後に、ここでイメージのレンダリングを完了することが可能となる。これは、初めはIOを保持する減衰値を倍加して集積イメージにすることによって達成され、ブロック158で達成される。これは、160に示すように、すべての立体がメモリにレンダリングされるまで、各ピクセルについて続けられる反復プロセスであることに留意する。
【0067】
そして、162に示すように、集積されたイメージが表示され、164においてプロセスが終了する。
次に、X線状イメージのレンダリングに関してのプログラムリストを提示する。
【0068】
ここで、この発明の2、3の実施の形態を説明してきたが、前述したものは単に例示的であり、制限するものではなく、例としてのみのために提示されていることは、当業者には明白であるはずである。多くの変形および別の実施の形態が当業者の意図にあり、添付する特許請求の範囲およびそれと同等のものによって定義されるこの発明の範囲内にあるものとして考慮される。
【0069】
【発明の効果】
以上のように、この発明によれば、物体をレンダリングしてX線イメージの形成をシミュレートすることでより単純でより非計算集約的な物体レンダリングシステム及び方法を得ることができ、物体の物質的及び幾何学的描写をその物体のX線イメージに変換するために提供され、それによって、X線の結果ではなく物体表現をX線フォーマットで行うことが可能となり、専門化されたアルゴリズムにより医療での用途および産業上の非破壊的な試験に使用することができる。
【図面の簡単な説明】
【図1】 この発明のシステムの組み合わされたデータの流れ図であり、外面および質量減衰を割り当てるためデジタル化された内部の立体を有する初期の3D物体、およびX線状のイメージを生成するよう組み合わされたシェル減衰イメージ(shell attenuation image)を与える本アルゴリズムの応用例を示す。
【図2】 図1に続く流れ図である。
【図3】 ベクトルがシミュレートされたX線ビームの減衰および方向について得られた非凸状立体を交差するシミュレートされたX線を示すものであり、この発明のシステムが非凸状面形状を有する立体を処理できることを示している説明図である。
【図4】 ネストされた立体の3DモデルからX線状イメージの生成するためのこの発明に係るアルゴリズムのフローチャートである。
【図5】 図4に続くフローチャートである。
【符号の説明】
10 この発明のシステム、12 3Dモデル、14,16,18 ネストされた立体、20 内側の表面をデジタル化、22 中間の表面をデジタル化、
24 外側の表面をデジタル化、26,28,30 質量減衰の割り当て、
38 ネストとして記録されたシェル、40 記録されたシェルの形状、
42、44、46 レンダラー、50 X線ビーム描写、52 カメラ描写、
54 内側シェルの深さイメージ、56 中間シェルの深さイメージ、58 外側シェルの深さイメージ、60、62、64 減衰イメージの計算、66 内側シェルの減衰イメージ、68 中間シェルの減衰イメージ、70 外側シェルの減衰イメージ、72 ピクセル毎にイメージを倍加。
Claims (7)
- 物体の材質及び幾何記述を模擬X線イメージとしてレンダリングする物体レンダリングシステムであって、
上記物体を、表面及び所定の質量減衰係数を有する、入れ子構造をした複数の立体の3Dモデルに分解する手段と、
上記3Dモデルを、模擬X線ビームのソース及び上記複数の立体に関連するシーングラフを含む3D座標系に登録する手段と、
上記複数の立体を通して上記模擬X線ビームを通過させる手段と、
上記模擬X線ビームが上記複数の立体の特定の1つに入射するかまたは出射するかを判断する手段と、
これに応答して上記模擬X線ビームの経路に沿った空所を位置決めする手段と、
上記材質及び幾何記述と上記空所に従って上記模擬X線ビームを減衰させる手段と、
上記模擬X線ビームの減衰に応答して、上記模擬X線イメージとして上記複数の立体を独立してレンダリングする手段と
を含む物体レンダリングシステム。 - 上記模擬X線ビームが上記複数の立体の特定の1つに入射するかまたは出射するかを判断する手段は、上記模擬X線ビームが上記特定の1つに入射するかまたは出射するポイントにおいて、上記模擬X線ビームに関連する減衰されない模擬フォトンの数に対応する大きさと、上記模擬X線ビームから上記ポイントにおいて散乱した上記模擬フォトンの方向に対応する方向とを有する、表面法線ベクトルを確立する手段と、上記模擬X線ビームの方向における一方向と、上記模擬X線ビームのソースと上記模擬X線ビームが上記特定の1つに入射するかまたは出射するポイントとの間の距離に等しい大きさを有する、視線ベクトルを確立する手段とを含むことを特徴とする請求項1記載の物体レンダリングシステム。
- 上記空所を位置決めする手段は、上記模擬X線ビームが入射または出射される上記特定の入り子の立体の表面上のポイントに対応するイメージピクセル位置を特定し、上記模擬X線ビームが上記表面に入射するかまたは出射するかを確立するために上記表面法線ベクトルと上記視線ベクトルとのドット積を計算する手段を含む請求項2記載の物体レンダリングシステム。
- 上記3Dモデルは、手動で形成される請求項1記載の物体レンダリングシステム。
- 上記3Dモデルは、実際の物理的な物体のCTスキャンから得られる請求項1記載の物体レンダリングシステム。
- 物体の材質及び幾何記述を上記物体の模擬X線イメージにレンダリングする物体レンダリング方法であって、
上記物体の入れ子構造をした複数の立体の3Dモデルを提供するステップと、
上記立体をデジタル化して、上記物体のそれぞれに対応する表面を提供するステップと、
質量減衰係数を、それぞれ対応する表面により境界づけられたそれぞれの立体に割り当てると共に、対応する質量減衰特性を上記対応する表面に割り当てるステップであって、上記表面はシェルを形成し、各シェルは対応する所定の質量減衰特性を有する、ステップと、
上記シェルをシーングラフに位置登録して、上記シェルの表面ジオメトリを関連させるための登録されたシェルジオメトリを提供するステップと、
上記質量減衰係数及び上記質量減衰特性に従って模擬X線ビームを減衰させながら、上記物体の上記材質及び幾何記述を通して上記模擬X線ビームを通過させるステップと、
上記シェルのそれぞれの深さイメージ画像を提供し、それにより、上記対応する立体を介して通過する上記模擬X線ビームの所定の光線によって横断される距離を示すように、上記シェルをレンダリングすることにより、上記模擬X線ビームが通過する上記シェルにより境界づけられたそれぞれの立体の経路長を得るステップと、
上記対応する質量減衰係数を使用して各シェルに対する減衰イメージを計算することによって、上記深さシェルイメージのそれぞれを関連する質量減衰係数と組み合わせ、それによってその密度を示すために影の付けられたシェル減衰イメージを提供するステップと、
上記影の付けられた減衰イメージを組み合わせることにより、そこで表現された立体を有する組み合わせられた模擬X線イメージを提供するステップと
を含む物体レンダリング方法。 - 上記シェルをレンダリングするステップは、所定の模擬X線ビーム記述および所定のカメラ記述を提供するステップと、上記模擬X線ビーム記述および上記カメラ記述を使用してレンダリングされた経路長を計算するステップとを含む請求項6記載の物体レンダリング方法。
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US09/102260 | 1998-06-22 | ||
US09/102,260 US6061469A (en) | 1998-06-22 | 1998-06-22 | Object rendering system to produce X-ray like images |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2000035406A JP2000035406A (ja) | 2000-02-02 |
JP4435329B2 true JP4435329B2 (ja) | 2010-03-17 |
Family
ID=22288964
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP13836199A Expired - Fee Related JP4435329B2 (ja) | 1998-06-22 | 1999-05-19 | 物体レンダリングシステム及び方法 |
Country Status (5)
Country | Link |
---|---|
US (1) | US6061469A (ja) |
EP (1) | EP0967578B1 (ja) |
JP (1) | JP4435329B2 (ja) |
AT (1) | ATE362628T1 (ja) |
DE (1) | DE69936078T2 (ja) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR20150141889A (ko) * | 2014-06-10 | 2015-12-21 | 한국건설기술연구원 | 3D X-ray CT 촬영을 이용한 재료의 이방성 측정방법 |
Families Citing this family (20)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6793496B2 (en) * | 1999-04-15 | 2004-09-21 | General Electric Company | Mathematical model and a method and apparatus for utilizing the model |
KR100769779B1 (ko) * | 2000-01-13 | 2007-10-24 | 닛토덴코 가부시키가이샤 | 광학 필름 및 액정 디스플레이 장치 |
DE10037491A1 (de) * | 2000-08-01 | 2002-02-14 | Stryker Leibinger Gmbh & Co Kg | Verfahren zum dreidimensionalen Visualisieren von Strukturen im Körperinneren |
US20040236550A1 (en) * | 2002-02-28 | 2004-11-25 | Edic Peter Michael | Mathematical model and a method and apparatus for utilizing the model |
EP1579202B1 (en) * | 2002-12-10 | 2007-04-18 | Commonwealth Scientific And Industrial Research Organisation | Radiographic equipment |
US7366278B2 (en) * | 2004-06-30 | 2008-04-29 | Accuray, Inc. | DRR generation using a non-linear attenuation model |
US9024949B2 (en) * | 2004-10-13 | 2015-05-05 | Sony Corporation | Object representation using distance functions |
US7974461B2 (en) * | 2005-02-11 | 2011-07-05 | Deltasphere, Inc. | Method and apparatus for displaying a calculated geometric entity within one or more 3D rangefinder data sets |
US7893938B2 (en) * | 2005-05-04 | 2011-02-22 | Siemens Medical Solutions Usa, Inc. | Rendering anatomical structures with their nearby surrounding area |
DE102005021068B4 (de) * | 2005-05-06 | 2010-09-16 | Siemens Ag | Verfahren zur Voreinstellung der Aufnahmeparameter bei der Erstellung von zweidimensionalen Durchlicht-Röntgenaufnahmen |
DE102005032686A1 (de) * | 2005-07-06 | 2007-01-11 | Carl Zeiss Industrielle Messtechnik Gmbh | Verfahren und Anordnung zum Untersuchen eines Messobjekts mittels invasiver Strahlung |
CN101743567A (zh) * | 2007-05-18 | 2010-06-16 | Uab研究基金会 | 虚拟互动式出席系统及方法 |
DE102009008115B4 (de) * | 2009-02-09 | 2012-05-24 | Siemens Aktiengesellschaft | Verfahren und Vorrichtung zur Erzeugung eines digital rekonstruierten 2-D-Bildes aus einem 3-D-Bilddatensatz |
US8831331B2 (en) * | 2009-02-10 | 2014-09-09 | Optosecurity Inc. | Method and system for performing X-ray inspection of a product at a security checkpoint using simulation |
JP2013231700A (ja) * | 2012-05-01 | 2013-11-14 | Tokyo Electron Ltd | X線検査方法及びx線検査装置 |
JP2014016239A (ja) * | 2012-07-09 | 2014-01-30 | Tokyo Electron Ltd | X線検査方法及びx線検査装置 |
US9091628B2 (en) | 2012-12-21 | 2015-07-28 | L-3 Communications Security And Detection Systems, Inc. | 3D mapping with two orthogonal imaging views |
CN105160705B (zh) * | 2015-09-10 | 2018-02-23 | 李晓强 | 一种基于三维引擎对移动终端三维图像数据的处理方法 |
EP3279865B2 (en) * | 2016-08-01 | 2022-02-16 | 3mensio Medical Imaging B.V. | Method, device and system for simulating shadow images |
WO2020217348A1 (ja) * | 2019-04-24 | 2020-10-29 | 三菱電機エンジニアリング株式会社 | 画像投影システム |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB2212267B (en) * | 1987-11-11 | 1992-07-29 | Circulation Res Ltd | Methods and apparatus for the examination and treatment of internal organs |
FR2624634B1 (fr) * | 1987-12-09 | 1994-04-29 | Schlumberger Etienne | Procede et dispositif permettant de reconstituer la forme et la position d'objets dans l'espace |
US4920491A (en) * | 1988-05-16 | 1990-04-24 | General Electric Company | Enhancement of image quality by utilization of a priori information |
US5226113A (en) * | 1989-10-30 | 1993-07-06 | General Electric Company | Method and apparatus for volumetric projection rendering using reverse ray casting |
-
1998
- 1998-06-22 US US09/102,260 patent/US6061469A/en not_active Expired - Fee Related
-
1999
- 1999-03-23 AT AT99105844T patent/ATE362628T1/de not_active IP Right Cessation
- 1999-03-23 DE DE69936078T patent/DE69936078T2/de not_active Expired - Lifetime
- 1999-03-23 EP EP99105844A patent/EP0967578B1/en not_active Expired - Lifetime
- 1999-05-19 JP JP13836199A patent/JP4435329B2/ja not_active Expired - Fee Related
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR20150141889A (ko) * | 2014-06-10 | 2015-12-21 | 한국건설기술연구원 | 3D X-ray CT 촬영을 이용한 재료의 이방성 측정방법 |
Also Published As
Publication number | Publication date |
---|---|
US6061469A (en) | 2000-05-09 |
DE69936078T2 (de) | 2008-04-30 |
DE69936078D1 (de) | 2007-06-28 |
EP0967578A3 (en) | 2001-04-04 |
ATE362628T1 (de) | 2007-06-15 |
EP0967578B1 (en) | 2007-05-16 |
JP2000035406A (ja) | 2000-02-02 |
EP0967578A2 (en) | 1999-12-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP4435329B2 (ja) | 物体レンダリングシステム及び方法 | |
JP4828422B2 (ja) | 適応直接ボリュームレンダリングの方法及びシステム | |
US9384325B2 (en) | System and method for generating enhanced density distribution in a three dimensional model of a structure for use in skeletal assessment using a limited number of two dimensional views | |
US8730240B2 (en) | Modeling and rendering of heterogeneous translucent materals using the diffusion equation | |
JP3483929B2 (ja) | 3次元画像生成方法 | |
US9001124B2 (en) | Efficient determination of lighting effects in volume rendering | |
CN109584349A (zh) | 用于渲染材料属性的方法和设备 | |
EP0262304A2 (en) | Method and system for solid modelling | |
US20060262112A1 (en) | System and method for three-dimensional shape generation from partial and incomplete views, and interactive design system using same | |
CA1315902C (en) | Minimization of directed points generated in three-dimensional dividing cubes method | |
US8149237B2 (en) | Information processing apparatus and program | |
Bellon et al. | aRTist–analytical RT inspection simulation tool | |
CN115828642A (zh) | 基于Unity的GPU加速X射线数字成像仿真方法 | |
Van Eyndhoven et al. | Region-based iterative reconstruction of structurally changing objects in CT | |
Vidal et al. | Simulation of X-ray Attenuation on the GPU. | |
CN107316282A (zh) | 一种图像降噪的系统和方法 | |
JPH0242585A (ja) | 空間内の物体の形状及び位置を再構成する方法及び装置 | |
Ito et al. | Orientation optimization and jig construction for xray ct scanning | |
US20110075899A1 (en) | Efficient correction of polychromy effects in image reconstruction | |
JP7178621B2 (ja) | 画像処理装置及び画像処理方法 | |
JP7178620B2 (ja) | 検査装置及び検査方法 | |
Inanc et al. | Human body radiography simulations: development of a virtual radiography environment | |
Albiol et al. | Densitometric radiographic imaging with contour sensors | |
JPH0728976A (ja) | 画像表示装置 | |
JP2809955B2 (ja) | グラフィックス・ワークステーションおよび三次元グラフィックス画像の生成方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20060330 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20080728 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20081209 |
|
A601 | Written request for extension of time |
Free format text: JAPANESE INTERMEDIATE CODE: A601 Effective date: 20090309 |
|
A602 | Written permission of extension of time |
Free format text: JAPANESE INTERMEDIATE CODE: A602 Effective date: 20090312 |
|
A601 | Written request for extension of time |
Free format text: JAPANESE INTERMEDIATE CODE: A601 Effective date: 20090409 |
|
A602 | Written permission of extension of time |
Free format text: JAPANESE INTERMEDIATE CODE: A602 Effective date: 20090414 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20090424 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20090630 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20090928 |
|
TRDD | Decision of grant or rejection written | ||
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20091208 |
|
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20091224 |
|
R150 | Certificate of patent or registration of utility model |
Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130108 Year of fee payment: 3 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130108 Year of fee payment: 3 |
|
S111 | Request for change of ownership or part of ownership |
Free format text: JAPANESE INTERMEDIATE CODE: R313113 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130108 Year of fee payment: 3 |
|
R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
LAPS | Cancellation because of no payment of annual fees |