JP2012221448A - ボリュームデータセット内の表面様構造を可視化するための方法及び装置 - Google Patents
ボリュームデータセット内の表面様構造を可視化するための方法及び装置 Download PDFInfo
- Publication number
- JP2012221448A JP2012221448A JP2011089858A JP2011089858A JP2012221448A JP 2012221448 A JP2012221448 A JP 2012221448A JP 2011089858 A JP2011089858 A JP 2011089858A JP 2011089858 A JP2011089858 A JP 2011089858A JP 2012221448 A JP2012221448 A JP 2012221448A
- Authority
- JP
- Japan
- Prior art keywords
- coordinate system
- calculated
- predetermined position
- local coordinate
- lcs
- 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.)
- Withdrawn
Links
- 238000000034 method Methods 0.000 title claims abstract description 38
- 239000013598 vector Substances 0.000 claims abstract description 90
- 238000005286 illumination Methods 0.000 claims abstract description 13
- 230000001131 transforming effect Effects 0.000 claims abstract description 4
- 238000004364 calculation method Methods 0.000 claims description 20
- 238000012545 processing Methods 0.000 claims description 17
- 238000009877 rendering Methods 0.000 claims description 15
- 238000002604 ultrasonography Methods 0.000 claims description 11
- 238000006243 chemical reaction Methods 0.000 claims description 10
- 238000004590 computer program Methods 0.000 claims description 5
- 238000012800 visualization Methods 0.000 claims description 5
- 239000000523 sample Substances 0.000 description 88
- 230000008859 change Effects 0.000 description 5
- 238000013507 mapping Methods 0.000 description 5
- 230000009466 transformation Effects 0.000 description 4
- SYUPLLHVMCLXEM-UHFFFAOYSA-N 4,5-dichloro-2-(3-methylphenyl)pyridazin-3-one Chemical compound CC1=CC=CC(N2C(C(Cl)=C(Cl)C=N2)=O)=C1 SYUPLLHVMCLXEM-UHFFFAOYSA-N 0.000 description 3
- 230000008901 benefit Effects 0.000 description 3
- 238000005070 sampling Methods 0.000 description 3
- 238000002059 diagnostic imaging Methods 0.000 description 2
- 238000001914 filtration Methods 0.000 description 2
- 239000011159 matrix material Substances 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 238000012285 ultrasound imaging Methods 0.000 description 2
- 238000013459 approach Methods 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 238000005266 casting Methods 0.000 description 1
- 238000002591 computed tomography Methods 0.000 description 1
- 238000005094 computer simulation Methods 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000004993 emission spectroscopy Methods 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 238000007429 general method Methods 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 238000002595 magnetic resonance imaging Methods 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 238000013139 quantization Methods 0.000 description 1
Images
Landscapes
- Apparatus For Radiation Diagnosis (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
- Image Generation (AREA)
Abstract
【課題】ボリュームデータセット内の表面様構造を可視化するための方法及び装置を提供する。
【解決手段】ボリュームデータセットのサンプル点Pにおいて局所座標系LCSiを定義する段階と、外部パラメータを全体座標系GCSから局所座標系LCSiに変換する段階と、サンプル点Piの局所座標系LCSi内の勾配ベクトル成分Gai,Gbi,Gciを計算する段階と、勾配ベクトル成分Gai,Gbi,Gciを使用してボリュームデータセットの所定位置Qにおける表面法線Nを計算する段階とを含む。所定位置Qにおける全体座標系GCSからの外部パラメータを、サンプル点Piの局所座標系LCSiの変換済み外部パラメータを使用することによって計算し、所定位置Qにおける計算済み表面法線Nと所定位置Qにおける計算済み外部パラメータを使用して、所定位置Qにおけるシェーディング又は照明を行う。
【選択図】図2
【解決手段】ボリュームデータセットのサンプル点Pにおいて局所座標系LCSiを定義する段階と、外部パラメータを全体座標系GCSから局所座標系LCSiに変換する段階と、サンプル点Piの局所座標系LCSi内の勾配ベクトル成分Gai,Gbi,Gciを計算する段階と、勾配ベクトル成分Gai,Gbi,Gciを使用してボリュームデータセットの所定位置Qにおける表面法線Nを計算する段階とを含む。所定位置Qにおける全体座標系GCSからの外部パラメータを、サンプル点Piの局所座標系LCSiの変換済み外部パラメータを使用することによって計算し、所定位置Qにおける計算済み表面法線Nと所定位置Qにおける計算済み外部パラメータを使用して、所定位置Qにおけるシェーディング又は照明を行う。
【選択図】図2
Description
本発明は、ボリュームレンダリング等のボリュームデータセット内の表面様構造を可視化し、特に超音波データ取得によるサンプルデータの実時間レンダリングのためにそのような表面様構造を照明するための方法及び装置に関する。本発明は、具体的には、三次元データボリュームから二次元画像をボリュームレンダリングする三次元医用画像化に関する。更に、本発明は、独立クレーム1及び11に記載された方法及び装置を実施することができるコンピュータプログラム製品に関する。
三次元超音波イメージングだけなくCT又はMRイメージングも、一般にデカルト座標系でない座標系でフォーマットされた三次元データを生成する。そのような非デカルト格子は、例えば人間の心臓の像を取得する例えば超音波プローブの音響格子(acoustic grid)である。体積を表わすそのようなデータは、通常、ボリュームレンダリングされる前にデカルト座標系に沿って走査変換される。走査された物体を表示するときは、物体の表面を表示する必要があり、そのような表面は、「できるだけ自然」に表示されなければならない。それを行う一方法は、仮想外部光源を使用し、そのような外部光源の寄与分をボリュームレンダリング積分に組み込むことである。そのような三次元データセットをレンダリングするためにはそのような三次元シェーディングが使用され、三次元シェーディングにはOpenGLやDirectX等のアプリケーションプログラミングインタフェースが使用される。
Blinn−Phongシェーディングモデル等の照明モデルの場合は、通常、選択された位置におけるシェーディングモデルを評価するために表面法線を定義しなければならない。Blinn−Phongシェーディングモデルは、Phong反射モデルの改良版であり、表示すべき表面の各点に対して実行される。これや他のシェーディングモデルの場合、通常、拡散反射を計算するために光ベクトルとそのような面の特定位置の表面法線を定義し、また正反射を計算するために更に他のベクトル(視線ベクトル又はカメラベクトル等)を定義しなければならない。
特定位置の表面法線を近似する一般的な方法は、前方差分法、後方差分法、中心差分法等の有限差分法(finite differencing scheme)を使用することである。ボリュームデータの場合、そのような表面法線は、次のような空間勾配ベクトルとも呼ばれる3つ全ての次元の差分(空間導関数)を使用して計算される。
スカラ場の勾配は、最も大きい変化の方向であり、従って、一定位置における等値面(isosurface)の表面法線は、次のような正規化勾配と等しい。
サンプルデータボリュームの勾配を計算する安価で迅速な方法は、中心差分を計算することであり、隣接格子値(neighbouring grid value)が使用される。これは、本出願の図1に示されている。
図1は、人間の心臓等の走査される物体の二次元断面画像を得るために超音波ビーム等のビーム1を送出している音響プローブの先端Aを示す。選択されたサンプル点Pにおける勾配Gを計算するために(ここでは、本発明をよく理解できるように二次元計算だけを示す)、2個の勾配成分Ga及びGbを計算しなければならない。なお、aとbは局所(有限)座標系の軸である。
Ga=(f(P+x)−f(P-x))/2
Gb=(f(P-y)−f(P+y))/2
式3
Gb=(f(P-y)−f(P+y))/2
式3
三次元ボリュームデータセットでは、6個の位置において、三次元ボリュームデータをサンプリングするfの値を求めなければならないことは明らかである。シェーディング作業に使用される中心サンプルf(P)を含め、各陰影位置(例えば、三次元ボクセル空間内の各ボクセル)について7個のfの値を求めなければならない。
先行技術では、各離散位置における勾配ベクトルを事前に計算することによってそのような計算を高速化することは既知である。三次元デカルト座標系において各離散ボクセル位置で差が計算され、元のスカラ値と一緒に4要素ベクトル(Ga,Gb,Gc,f)で記憶される。性能の点では、そのような三次元ボクセル空間内の面をシェーディングするときの参照数が減少する。シェーディングは、勾配ベクトル成分とスカラ値だけを使用することで行われ、事前に計算された4要素ベクトルを使用することによってそれらの位置を照明できる。しかしながら、主な欠点は、データサイズの増大である。ボリュームデータは少なくとも4倍に増え、これは、少ないメモリ量のマシンでは問題になりやすい。
非デカルトデータの場合、空間導関数(勾配)とその結果得られる表面法線の数値計算はもっと複雑になり、照明計算は、全域デカルト座標系に基づき、光源、観察者位置(カメラ位置)等は、そのようなデカルト座標系で与えられる。Blinn−PhongやPhong等の幾つかのシェーディングモデルは全て、シェーディング計算にデカルト座標系を使用する。
例えばいわゆる音響座標で与えられる超音波データ等の非デカルトデータの場合、Blinn−PhongモデルやPhong反射モデル等のモデルをそのまま使用することはできず、このために前述のように所定のサンプル位置での勾配ベクトルが使用され、前述の差分法のいずれか使用される。これは座標系がデカルト座標系ではないためである。これは図1にも示されている。軸aと軸bを有する局所座標系は、値f(P+x)及びf(P-x)がそのような空間内にないため、方向aの正しい勾配ベクトル成分を提供しない。従って、勾配ベクトルGを計算するために、最初に三次元データセットを再サンプリングしてデカルト座標系のデータにして勾配成分を計算し、これにより、勾配成分は全域デカルト座標系で提供される。
別の手法が特許文献1に示されており、この特許では、最初に、音響ドメイン内のデータから勾配が決定され、次にデカルト座標又は表示スクリーンドメインに変換される。変換マトリックスは、サンプリングされた三次元音響データ内の各位置ごとに指定される。前述の解決策は、高価なデカルト座標配置(geometry)へのサンプリング段階を必要とする。
従って、本発明の目的は、デカルト座標系で提供されないボリュームデータセット内の表面様構造を可視化し、それにより既知の照明又はシェーディングモデルの使用を可能にする方法及び装置を提供することである。具体的には、本発明は、実時間の画像処理ユニット(GPU)を使用した音響格子におけるボリュームレンダリングを行うための高速勾配計算及び表面照明のソリューションを提供するものとする。
本発明は、これらの目的を独立クレーム1及び11の特徴によって解決する。好ましい利点は、対応するそれぞれの従属クレーム内に定義され請求されている。
本発明の主な利点は、データをデカルト座標配置に再サンプリングしなくてもよいことであり、本発明は、また、計算をより簡単に高速で実行し、それにより既知のシェーディング及び照明モデルの使用を可能にすることである。
本発明の方法は、
a)前記ボリュームデータセットのサンプル点(P)において局所座標系(LCSi)を定義する段階と、
b)前記サンプル点(Pi)の前記局所座標系(LCSi)内の勾配ベクトル成分(Gai,Gbi,Gci)を計算する段階と、
c)前記外部パラメータを全体座標系(OCS)から選択されたサンプル点又は基準格子の点における前記局所座標系(LCSi)に変換する段階と、
d)前記勾配ベクトル成分(Gai,Gbi,Gci)を使用して、前記ボリュームデータセットの所定位置(Q)における表面法線(N)を計算する段階と、
e)前記選択されたサンプル点又は基準格子の点における前記局所座標系(LCSi)の変換済み外部パラメータを使用することによって、前記所定位置(Q)における前記全体座標系(GCS)からの外部パラメータを計算する段階とを含む。
a)前記ボリュームデータセットのサンプル点(P)において局所座標系(LCSi)を定義する段階と、
b)前記サンプル点(Pi)の前記局所座標系(LCSi)内の勾配ベクトル成分(Gai,Gbi,Gci)を計算する段階と、
c)前記外部パラメータを全体座標系(OCS)から選択されたサンプル点又は基準格子の点における前記局所座標系(LCSi)に変換する段階と、
d)前記勾配ベクトル成分(Gai,Gbi,Gci)を使用して、前記ボリュームデータセットの所定位置(Q)における表面法線(N)を計算する段階と、
e)前記選択されたサンプル点又は基準格子の点における前記局所座標系(LCSi)の変換済み外部パラメータを使用することによって、前記所定位置(Q)における前記全体座標系(GCS)からの外部パラメータを計算する段階とを含む。
外部パラメータは、データセットの可視化(ボリュームレンダリング)を決定する位置及び/又は方向をそれぞれ定義する位置座標及び/又はベクトルを含むことが好ましい。外部パラメータは、例えば、1個又は幾つかの光源の位置、観察者位置、レンダリングに使用される光ベクトル、及び/又は観察者のベクトル/方向でよい。
例えば様々な観察者方向から表面構造を可視化するには、これらの外部パラメータを実時間で変更できることが望ましいので、本発明によって提供されるような高速アルゴリズムは重要である。
一実施形態によれば、ボリュームデータセットは、超音波画像診断、コンピュータ断層撮影法、磁気共鳴映像法、ポジトロン放射分光法等の医用画像技術によって取得された。サンプル点は、データセットを構成する点である。
局所座標系は、各サンプル点において、ボリュームデータセットの座標系(例えば、音響座標系)と整合された座標系であることが好ましい。従って、局所座標系は、各サンプル点ごとに僅かに異なり、局所座標系は、サンプル点間で滑らかに変化することが好ましい。
全体座標系は、ボリュームデータセットの位置と向き、並びに1個又は幾つかの光源の位置、観察者の位置、レンダリングに使用される光ベクトル、及び/又は観察者の方向等のボリュームレンダリング段階で必要とされる外部パラメータが全て同一の座標系で定義される座標系である。これは、一般にデカルト座標系であるが、必要条件ではない。
サンプル点の前記局所座標系内の勾配ベクトル成分を計算することによって、サンプル点(及び、対応する勾配ベクトル成分)を全域デカルト座標系に変換する必要はない。実際に、光源又はカメラ若しくは観察者の位置等の外部パラメータは、全体座標系から各サンプル点又は選択されたサンプル点の局所座標系に変換される。外部パラメータのそのような変換は、シェーディングすべき位置でそのような局所座標系を定義することによって可能である。
可能な局所座標系は、接ベクトル空間であり、例えば局所的バンプマッピング計算に使用される。ここで、レンダリングされる物体の表面法線の摂動は、各サンプル点におけるテクスチャマップで参照され、照明計算が行われる前に適用される。それにより、より豊かで細やかな曲面表現が得られる(例えば、Phongシェーディングを参照)。ノーマルマッピングとパララックスマッピングは、最も一般に使用されるバンプマッピング技術である。この場合、局所座標系を定義するために接ベクトル空間がしばしば使用される。
表面の接ベクトル空間を作成するには、3本の垂直軸(即ち、T,B,N)を計算しなければならない。T(接線ベクトル)は、パラメトリック曲面上でS又はTの増大する方向と平行である。N(法線ベクトル)は、局所表面と垂直である。B(従法線)は、NとTの両方に垂直であり、Tと同様に表面上にある。これは、表面に付けられ且つシェーディングされるサンプル点と共に移動する「移動座標系」である。
音響座標の場合には、3本の直交軸が、3方向のアジマス、エレベーション及びレンジと位置合わせされることが好ましい。例えば、Nはレンジに対応し、Tはエレベーションに対応し、Sはアジマスに対応する。各サンプル点における局所座標系は、同じように音響座標と整合されることが好ましい。
本発明は、各サンプル点においてそのようなサンプル点の局所座標系内で勾配ベクトル成分が計算される技術を使用する。前記サンプル点は、特殊な座標系である音響格子の格子点であることが好ましく、その理由は、例えば、音響格子の方向の直交軸(例えば、レンジ、エレベーション及びアジマス)を有する少なくとも三次元座標系超音波画像が走査され、次にそのようなサンプル点で局所座標系が定義されるからである。前記局所座標系内の勾配ベクトル成分の計算は、例えば、好ましくは、それぞれの隣接するサンプル点の離散値を用い、サンプル点の有限差分を求めることによって全ての次元の空間導関数を計算することで行うことができる。
前記ボリュームデータセットの所定位置(サンプル点又はそのようなサンプル点間のボクセル)を照明すべき場合は、サンプル点の勾配ベクトル成分を使用することによってそのような所定位置における表面法線を計算することができる。所定位置は、ボリュームレンダリングされた画像を作成するために光反射を計算することが望まれる位置である。可視化又はボリュームレンダリングプロセスが、レイキャスティング法によって行われる場合は、所定位置は、例えばそのような光線上にあり、スライスに基づく方法が使用される場合は、所定位置は、そのようなスライス上の画素位置になる。所定位置の全てではないにしても殆どは、通常、サンプル点と同一でなく、選択されたサンプル点とも基準格子上の点とも同一でない。
本発明の好ましい実施形態によれば、所定位置における前記表面法線は、最初に隣接するサンプル点の非正規化勾配ベクトル成分を補間し、次にその位置の補間された勾配ベクトルを正規化することによって計算される。
次に所定位置における表面法線を計算しなければならないとき、本発明の好ましい実施形態により、例えば8個の隣接するサンプル点の非正規化勾配ベクトル成分を補間することができ、サンプル点は、そのような三次元(ボクセル)空間内のそのような所定位置に隣接している。それにより、計算がサンプル点だけに減少するが、全ての他の(ボクセル)位置は、それぞれの隣接するサンプル点の勾配ベクトル成分を補間することにより計算される。
本発明の重要な様相によれば、所定位置での対応する勾配ベクトルを得るために、これらのサンプル点の正規化勾配ベクトル又は法線ベクトルを補間するのではなく、非正規化勾配ベクトル成分を補間し、後でそのような所定位置で正規化を行なう。
一実施形態によれば、隣接するサンプル点の勾配ベクトル成分間のそのような補間は、勾配ベクトル成分が異なる局所座標系に存在するという事実を考慮せずに行われる。換言すると、補間は、寄与する全てのサンプル点からの勾配ベクトル成分が同一の局所座標系内にあるかのように行われる。この近似は、少なくとも音響座標では座標系の幾何学形状に鋭い変化がなく、即ち隣接するサンプル点の局所座標系の向きは僅かしか相違しないために行うことができる。この状況は、互いに離れたサンプル点を比較しなければならない場合には異なるが、そうでないときは、この手順から生じる不正確さは結果を損なわない。
本発明の別の様相によれば、光源(及び他の外部パラメータ)は、三次元データ格子の幾つかの点における局所座標系に変換される。これらの点は、選択されたサンプル点でもよく、音響格子内で定義されるが音響格子のサンプル点と一致しなくてもよい基準格子の点でもよい。本発明の一実施形態によれば、外部パラメータの局所座標系へ転換/変換は、極めて計算集約的になるので、全てのサンプル点で行われるわけではなく、データセットの幾つかの選択された点(例えば、各次元で3番目〜50番目ごとサンプル点)又は事前に定義された基準格子の点でのみ行われ、その場合、基準格子は、サンプル点の分解能より低い分解能(例えば、5〜30分の1)を有する。
外部パラメータのサンプル点又は基準格子上の点の局所座標系への変換は、例えば、全体座標系とそれぞれの局所座標系との間の変換マトリックスを計算し、外部座標をサンプル点又は基準点における局所座標系にマッピングすることによって行われる。
次に、本発明の好ましい実施形態によれば、前記所定位置における前記全体座標系からの外部パラメータを、選択されたサンプル点又は基準格子の点、好ましくは隣接する選択されたサンプル点又は基準格子の隣接する点の前記局所座標系の変換済み外部パラメータを使用することによって計算することもできる。これにより、外部パラメータの各所定位置における局所座標系への変換を計算する必要がなく、前記三次元(音響)データ格子内の選択されたサンプル点に関してそのような計算を行うだけでよいので計算時間が更に短縮される。前述のように、これらのサンプル点は、米国特許出願番号2005/0253841 A1に記載されたような事前に定義された基準格子の格子点であってもよい。従って、一実施形態によれば、デカルト(全体座標系)座標配置から音響(局所座標系)座標配置への変換は、全てのサンプル点に関して行われるわけではなく、より粗いレベルで行われる。
一実施形態によれば、所定位置における外部パラメータは、次に、外部パラメータが局所座標系に変換された最も近いサンプル点又は基準格子の点から補間される。これは近似であるが、この場合も、幾何学形状に鋭い変化がないので、補間を使用して十分に滑らかな結果が得られる。
本発明の別の好ましい様相によれば、所定位置のシェーディングは、従来の照明モデルによって行われ、前記所定位置における表面法線は、前記所定位置における前記全体座標系からの計算済み外部パラメータと共に使用される。
本発明によれば、サンプリングされた生の超音波データから導出されたボリュームデータセットを使用することができ、局所座標系内の前記勾配ベクトル成分は、生の超音波データの各サンプル点から計算され、4要素ベクトル(Ga,Gb,Gc,F)を含む新しいデータセットにサンプル点のスカラ値と一緒に記憶される。任意の所定位置のために、これらの勾配ベクトル成分と隣接するサンプル点のそのようなスカラ値が、次に、所定位置における表面法線を計算するために使用される。その後で、前記サンプル点の局所座標系の変換された外部パラメータは、所定位置における外部パラメータを計算するために使用され、所定位置は、次に、計算済み外部パラメータと計算済み表面法線を使用することによって従来の照明モデルによって照明される。表面法線、光ベクトル、観察者ベクトル(viewers vector)、半正規、反射ベクトル等の関連要素が、各サンプル点の局所座標系内にあり、各所定位置ごとに補間されるので、従来のBlinn−Phong照明モデルを使用することができる。
本発明によれば、サンプリングされた生の超音波データは、最初に、グラフィック処理装置(GPU)に送られ、次にそこでシェーダー(shader)を使用して計算され、全体座標系から外部パラメータの局所座標系への変換は、例えば、超音波コンピュータシステムのグラフィック処理ユニット(GPU)によって実行される。次に、局所座標系内の勾配ベクトル成分が、中心差分、ソーベル演算子又は他の推定量を使用することにより計算される。
本発明は、また、ボリュームデータセット内の表面様構造を可視化するための装置であって、
a)前記ボリュームデータセットのサンプル点において局所座標系を定義し、且つ
b)外部パラメータを全体座標系から前記局所座標系に変換するための第1の演算処理装置(PU1)と、
c)前記サンプル点の前記局所座標系内の勾配ベクトル成分を計算し、且つ
d)前記勾配ベクトル成分を使用して前記ボリュームデータセットの所定位置における表面法線を計算するための第2の演算処理装置(PU2)とにより可視化するための装置に関する。
a)前記ボリュームデータセットのサンプル点において局所座標系を定義し、且つ
b)外部パラメータを全体座標系から前記局所座標系に変換するための第1の演算処理装置(PU1)と、
c)前記サンプル点の前記局所座標系内の勾配ベクトル成分を計算し、且つ
d)前記勾配ベクトル成分を使用して前記ボリュームデータセットの所定位置における表面法線を計算するための第2の演算処理装置(PU2)とにより可視化するための装置に関する。
両方の演算処理装置(PU1,PU2)は、CPUが照明又はシェーディングに必要な処理時間を最小にするのを支援するグラフィック演算処理装置であることが好ましい。グラフィック処理装置(GPU)は、一般に表示スクリーンドメインと同一のデカルト座標系内の計算済み外部パラメータと計算済み表面法線を使用することによって、従来の照明モデルによって照明された所定位置を表示する。この場合、グラフィック処理装置は、一般に、また、所定位置における全体座標系からの外部パラメータを、隣接するサンプル点の局所座標系の変換済み外部パラメータを補間してそのような局所格子点における変換され補間された関連する外部パラメータを得ることによって計算しており、そのような点は、音響格子の点(又は、これらの格子点間の点)、又は米国出願番号2005/0253841 A1に記載されたような基準格子の点である。
米国特許出願番号2005/0253841 A1は、基準格子(座標系)の変形を三次元データセット内で定義された切断平面の関数として開示している。米国特許出願番号2005/0253841 A1は、走査ボリュームを表わす代理幾何学形状(proxy geometry)と切断平面の交点を使用することを提案し、グラフィック処理ユニット(GPU)の頂点プロセッサは、次に、基準格子を変形し、デカルト座標とテクスチャ座標を基準格子のそれぞれの格子点に関して決定する。しかしながら、本発明では、格子点をデカルト座標に変換しなくてもよい。
本発明は、また、コンピュータ可読媒体に記憶されたプログラムコードを含むコンピュータプログラム製品であって、そのようなコンピュータプログラム製品がコンピュータ上で実行された場合に前述の方法を実行するコンピュータプログラム製品に関する。
本発明の好ましい実施形態は、添付図面について説明することにより更に詳しく説明される。
図1は、頂点A、ビーム1及びロー2を有する従来の超音波(音響)格子を示し、ビーム1とロー2は、二次元表現の音響座標系(例えば、それぞれレンジとアジマス)を表わす。本発明をよりよく理解するために、図1には第3次元(例えば、球座標系のエレベーション)は示されていない。
前述のように、局所勾配ベクトルGは、サンプル点Pの局所勾配ベクトル成分Ga及びGbによって定義され、これらの成分は、前述のような差分方法によって計算される。
図2は、光源LS等の外部パラメータの全体座標系GCSから局所座標系(即ち、第1の局所座標系LCS1と第2の局所座標系LCS2)への変換を概略的に示し、これらの局所座標系は、図1に示されたような音響格子では第1のサンプル点P1と第2のサンプル点P2で定義される。
各サンプル点Piは、局所勾配ベクトル成分Gai及びGbiを有する局所勾配ベクトルGiを有する。第1のサンプル点P1は、第1の勾配ベクトル成分Ga1と第2の勾配ベクトル成分Gb1を有し、第1の局所座標系LCS1の軸b1は、ビーム1と平行であり、一方第2の局所座標系軸a1は、ロー2の接線ベクトルであり且つ第1の座標系軸b1に垂直である。
第1の局所勾配ベクトルG1は、隣接するサンプル点(サンプル点P2や他のサンプル点等)を使用することにより定義される。この際、中心差分法(式1)等の周知の差分法、ソーベルフィルタ等の離散的畳み込みフィルタリング、又は立方体Bスプラインやその派生物等の連続的畳み込みフィルタリング等が使用される(隣接する4x4x4)。
図2は、また、軸x及びyを有する全体座標系GCSと、そのような全体座標系GCSで与えられた光ベクトルLGを示す。サンプル点P1,P2を全体座標系GCSに再サンプリングする代わりに、光源LSが局所座標系LCS1及びLCS2に変換され、それにより第1の局所光ベクトルLL1と第2の局所光ベクトルLL2が得られる。一実施形態によれば、光源等の外部パラメータは、全てのサンプル点の局所座標系に変換されるわけではなく、サンプル点のより粗い格子の局所座標系又は基準格子の点の局所座標系だけに変換される。
図2は、また、所定位置Qを示し、そのような所定位置Qの局所光ベクトルLQは、必ずしも光源LSを所定位置Qの局所座標系LCSQに変換することによって計算されるわけではなく、そのような光ベクトルLQは、サンプル点P1及びP2の第1の局所光ベクトルLL1及び第2の局所光ベクトルLL2と、必要に応じて、更に他のサンプル点の更に他の局所光ベクトルを使用することにより計算されることが好ましい。従って、外部パラメータは、音響格子又は前記基準格子の選択された隣接するサンプル点から、既に変換された外部パラメータを補間することによって、所定位置Qの局所座標系LCSQに変換される。
同様に、所定位置Qの法線ベクトルNは、第1のサンプル点P1の勾配ベクトル成分Ga1,Gb1及びGC1並びに第2のサンプル点P2(図示せず)の勾配ベクトル成分Ga2,Gb2及びGc2を使用することによって計算され、従って、座標系又は外部パラメータを変換することにより各所定位置における照明モデルに必要なベクトルを計算しなくてもよく、離散的サンプル点における外部パラメータの変換値を計算し、それぞれのサンプル点P1の非正規化勾配ベクトル成分を補間することにより所定位置Qにおける表面法線N(及び他のベクトル)を得るだけでよい。この場合も、外部パラメータは、音響格子又は前記基準格子の選択された隣接するサンプル点から、既に変換された外部パラメータを補間することによって、所定位置Qの局所座標系LCSQに変換される。
そのようにして、局所空間内だけで音響格子又は基準格子の格子点でだけの光源LSを計算し、また中間の点又はピクセルについての直線補間を使用することによって、計算コストを削減することができる。それにより、変換計算コストの大部分を削減することができる。何故なら、格子が粗いためである。
本発明は、各所定位置(三次元デカルトボクセル空間内の各ボクセル)における値は補間によって求められることを受け入れることにより、格子(音響格子、基準格子、又はこれらの格子の一部分若しくは組み合わせ)の内の幾つか又は全ての点だけにおいて外部パラメータを変換することによって計算モデルを縮小できる。そのような補間結果は照明計算だけに使用されるので、これは、理にかなった近似であると考えられる。人間の目は、照明の高周波数変化に極めて敏感である(例えば、鋭い縁や硬い切り口は見えやすい)。しかしながら、幾何学形状に鋭い変化がないので、本発明の補間技術を使用することによりかなり滑らかな結果を得ることができる。人間の目は、通常、照明結果の違いに気付かない。
図3は、生命体三次元超音波走査等の実時間用途の本発明の別の好ましい実施形態を示す。実時間用途が必要なとき勾配計算がしばしば障害になる。現在、処理前の(pre-process)勾配がCPUで計算され、後でGPUに転送される。本発明によれば、レンダリング段階で元の生データを転送して予備計算を実行することが可能となった。最新ハードウェアは、三次元ボリュームテクスチャの個別スライスへのそのようなレンダリングを支援する。
図3の段階Iに示されたような生データは、最初に、その固有の形式(例えば、11サンプル当たり8ビット)でGPUに転送される。第2の段階IIで、各サンプルが繰り返され、勾配が計算される。中心差分、ソーベル演算子又は他の推定量を使用することができる。その後で、非正規化勾配成分(Ga,Gb,Gc,F)が、図3に概略的に示されたような新しい三次元ボリュームで記憶される。ここで、いかなる量子化アーティファクトも回避するために高精度(例えば、16又は32ビット浮動)を使用することが好ましい。
本発明によれば、事前計算された勾配をこの段階で正規化するのではなく、その代わりに正規化されていない勾配ベクトル成分を記憶することが有利である。
図3に段階IIIに示されたように、所定位置Qにおける表面法線Nは、第1のサンプル点P1の勾配G1、第2のサンプル点P2の勾配G2、第3のサンプル点P3の勾配G3、及び第4のサンプル点P4の勾配G4の勾配ベクトル成分を補間することにより計算される。各勾配は、それぞれの勾配ベクトル成分を有し、これらの勾配ベクトル成分は全て所定位置Qにおける勾配ベクトルを定義するために使用され、この勾配ベクトルを正規化するだけで表面法線Nを得ることができる。
本発明は、使用する方法が基準格子の採用に上手く適合するという利点を有する。何故なら、ピクセルシェーダ段での作業負荷が少なく高性能な実時間勾配レンダリングが得られるためである。更に、全ての照明計算は局所座標系で行われる。音響空間内の勾配ボリュームは、照明が迅速化されるように事前に計算し記憶することが可能である。
Claims (14)
- ボリュームデータセット内の表面様構造を可視化する方法であって、可視化は、可視化を定義する位置座標及び/又は方向を含む1組の外部パラメータによって決定され、
a)前記ボリュームデータセットのサンプル点(P)において局所座標系(LCSi)を定義する段階と、
b)前記サンプル点(Pi)の前記局所座標系(LCSi)内の勾配ベクトル成分(Gai,Gbi,Gci)を計算する段階と、
c)前記外部パラメータを全体座標系(GCS)から選択されたサンプル点又は基準格子の点における前記局所座標系(LCSi)に変換する段階と、
d)前記勾配ベクトル成分(Gai,Gbi,Gci)を使用して前記ボリュームデータセットの所定位置(Q)における表面法線(N)を計算する段階と、
e)前記選択されたサンプル点又は基準格子の点における前記局所座標系(LCSi)の変形済み外部パラメータを使用することによって、前記所定位置(Q)における前記全体座標系(GCS)からの外部パラメータを計算する段階とを含む方法。 - 前記所定位置(Q)における前記表面法線(N)が、最初に隣接するサンプル点(P)の非正規化勾配ベクトル成分(Gai,Gbi,Gci)を補間し、次にそのような位置(Q)における前記補間された勾配ベクトル(G)を正規化することによって計算される、請求項1に記載の方法。
- 前記サンプル点(P)が、音響格子の格子点であり、前記局所座標系(LCSi)が少なくとも3本の直交軸を有する、請求項1又は2に記載の方法。
- 前記所定位置(Q)が、サンプル点(P)間の位置であり、前記所定位置(Q)における前記表面法線(N)が、少なくとも8個の隣接するサンプル点(P)における非正規化勾配ベクトル成分(Gai,Gbi,Gci)を補間することにより計算される、請求項1〜3の何れかに記載の方法。
- 空間導関数を好ましくは有限差分法によって全ての次元で計算することによって、前記サンプル点(Pi)の前記局所座標系(LCSi)内の前記勾配ベクトル成分(Gai,Gbi,Gci)を計算する、請求項1〜4の何れかに記載の方法。
- 前記所定位置(Q)における前記全体座標系(GCS)からの外部パラメータが、最も近い選択されたサンプル点又は基準格子の最も近い点からの変換済み外部パラメータ間で補間することにより計算される、請求項の1〜5の何れかに記載の方法。
- g)前記所定位置(Q)における前記表面法線(N)と前記所定位置(Q)における前記全体座標系(GCS)からの前記計算済み外部パラメータを使用して前記所定位置(Q)を従来の照明モデルによってシェーディングする追加の段階を含む、請求項の1〜6の何れかに記載の方法。
- 前記ボリュームデータセットが、サンプリングされた生の超音波データから導出され、前記局所座標系(LCSi)内の前記勾配ベクトル成分(Gai,Gbi,Gci)が、前記生の超音波データの各サンプル点(Pi)ごとに計算され、前記サンプル点(Pi)のスカラ値(F)と一緒に新しいデータセット(Gai,Gbi,Gci,F)に記憶され、
任意の所定位置(Q)に関して、前記勾配ベクトル成分(Gai,Gbi,Gci)と隣接するサンプル点(P1)の前記スカラ値(F)が、前記所定位置(Q)における表面法線(N)を計算するために使用され、
前記選択されたサンプル点又は基準格子の点の前記局所座標系(LCSi)の前記変換済み外部パラメータが、前記所定位置(Q)における外部パラメータを計算するために使用され、前記所定位置(Q)が、前記計算済み外部パラメータと前記計算済み表面法線(N)を使用することによって従来の照明モデルによって照明される、請求項の1〜7の何れかに記載の方法。 - サンプリングされた生の超音波データが、最初にグラフィック処理装置(GPU)に送られ、そこで頂点シェーダ(vertex shader)を使用して計算が行われ、全体座標系(GCS)からの外部パラメータの前記局所座標系(LCSi)への前記変換が超音波コンピュータシステムの演算処理装置によって実行される、請求項8に記載の方法。
- 前記局所座標系(LCSi)内の前記勾配ベクトル成分(Gai,Gbi,Gci)が、中心差分、ソーベル演算子又は他の推定量を使用して計算される、請求項1〜9のいずれか1項に記載の方法。
- ボリュームデータセットの表面様構造を可視化するための装置であって、
a)前記ボリュームデータセットのサンプル点(P)において局所座標系(LCSi)を定義し、且つ
b)外部パラメータを全体座標系(GCS)から前記局所座標系(LCSi)に変換するための第1の演算処理装置(PU1)と、
c)前記サンプル点(Pi)の前記局所座標系(LCSi)内の勾配ベクトル成分(Gai,Gbi,Gci)を計算し、且つ、
d)前記勾配ベクトル成分(Gai,Gbi,Gci)を使用して、前記ボリュームデータセットの所定位置(Q)における表面法線(N)を計算するための第2の演算処理装置(PU2)とを含む装置。 - 前記グラフィック処理装置(GPU)が、前記所定位置(Q)を表示し、前記所定位置(Q)が、前記計算済み外部パラメータと前記計算済み表面法線(N)を使用することによって、三次元ボリュームレンダリング中に、デカルト座標系で従来の照明モデルによって照明される、請求項11に記載の装置。
- 前記グラフィック処理装置(GPU)が、また、前記所定位置(Q)における前記全体座標系(GCS)から外部パラメータを、最も近い選択されたサンプル点又は基準格子の最も近い点の変換された外部パラメータ間を補間することによって計算する、請求項11又は12に記載の装置。
- コンピュータ可読媒体に記憶され且つコンピュータ上で実行された場合に請求項1〜10の何れかに記載された方法を実行するプログラムコードを含むコンピュータプログラム製品。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2011089858A JP2012221448A (ja) | 2011-04-14 | 2011-04-14 | ボリュームデータセット内の表面様構造を可視化するための方法及び装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2011089858A JP2012221448A (ja) | 2011-04-14 | 2011-04-14 | ボリュームデータセット内の表面様構造を可視化するための方法及び装置 |
Publications (1)
Publication Number | Publication Date |
---|---|
JP2012221448A true JP2012221448A (ja) | 2012-11-12 |
Family
ID=47272829
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2011089858A Withdrawn JP2012221448A (ja) | 2011-04-14 | 2011-04-14 | ボリュームデータセット内の表面様構造を可視化するための方法及び装置 |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP2012221448A (ja) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2016174905A (ja) * | 2015-03-18 | 2016-10-06 | 株式会社東芝 | 超音波診断装置、医用画像処理装置及び医用画像処理プログラム |
WO2018159709A1 (ja) * | 2017-02-28 | 2018-09-07 | 富士フイルム株式会社 | 3次元画像処理装置および方法並びにプログラム |
CN110060347A (zh) * | 2019-03-01 | 2019-07-26 | 新疆远山矿产资源勘查有限公司 | 一种球面坐标***下生成等值线的方法 |
-
2011
- 2011-04-14 JP JP2011089858A patent/JP2012221448A/ja not_active Withdrawn
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2016174905A (ja) * | 2015-03-18 | 2016-10-06 | 株式会社東芝 | 超音波診断装置、医用画像処理装置及び医用画像処理プログラム |
WO2018159709A1 (ja) * | 2017-02-28 | 2018-09-07 | 富士フイルム株式会社 | 3次元画像処理装置および方法並びにプログラム |
JPWO2018159709A1 (ja) * | 2017-02-28 | 2019-08-08 | 富士フイルム株式会社 | 3次元画像処理装置および方法並びにプログラム |
US10832459B2 (en) | 2017-02-28 | 2020-11-10 | Fujifilm Corporation | Three-dimensional image display apparatus for displaying target to which color is assigned |
CN110060347A (zh) * | 2019-03-01 | 2019-07-26 | 新疆远山矿产资源勘查有限公司 | 一种球面坐标***下生成等值线的方法 |
CN110060347B (zh) * | 2019-03-01 | 2023-12-26 | 新疆远山矿产资源勘查有限公司 | 一种球面坐标***下生成等值线的方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US8736609B2 (en) | Method and device for visualizing surface-like structures in volumetric data sets | |
JP5120926B2 (ja) | 画像処理装置、画像処理方法およびプログラム | |
US5831623A (en) | Volume rendering apparatus and method | |
US7239319B2 (en) | Rendering outline fonts | |
US9177416B2 (en) | Space skipping for multi-dimensional image rendering | |
JP5288795B2 (ja) | ボリュームデータの画像処理 | |
US20060262112A1 (en) | System and method for three-dimensional shape generation from partial and incomplete views, and interactive design system using same | |
US20050143654A1 (en) | Systems and methods for segmented volume rendering using a programmable graphics pipeline | |
US7893940B2 (en) | Super resolution contextual close-up visualization of volumetric data | |
US20050237336A1 (en) | Method and system for multi-object volumetric data visualization | |
US6563501B2 (en) | Bicubic surface rendering | |
US20060227131A1 (en) | Flat texture volume rendering | |
US20040138560A1 (en) | Real-time scan conversion and rendering of ultrasound data | |
US7586501B2 (en) | Simultaneous projection of multi-branched vessels and their context on a single image | |
US7602180B2 (en) | System and method for fast texture-based tensor field visualization for DT-MRI | |
CN104933749B (zh) | 图形图元的裁剪 | |
US7893938B2 (en) | Rendering anatomical structures with their nearby surrounding area | |
Jense | Voxel-based methods for CAD | |
CN105825471A (zh) | 一种基于Unity 3D的三维体表面重构与渲染方法 | |
JP4885042B2 (ja) | 画像処理方法および装置ならびにプログラム | |
Sakas | Interactive volume rendering of large fields | |
JP2006000127A (ja) | 画像処理方法および装置並びにプログラム | |
JP2012221448A (ja) | ボリュームデータセット内の表面様構造を可視化するための方法及び装置 | |
US20070188492A1 (en) | Architecture for real-time texture look-up's for volume rendering | |
US20040179010A1 (en) | Apparatus and method for semi-automatic classification of volume data |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A300 | Application deemed to be withdrawn because no request for examination was validly filed |
Free format text: JAPANESE INTERMEDIATE CODE: A300 Effective date: 20140701 |