何麗君，王曉強，云 健，包書哲

(大連民族學院計算機科學與工程學院，遼寧大連 116605)

三維數據場體繪制研究進展

何麗君，王曉強，云 健，包書哲

(大連民族學院計算機科學與工程學院，遼寧大連 116605)

介紹了體數據場按物理特征和幾何特征的分類，綜述了三維數據場體體視化的兩類繪制算法——面繪制和直接體繪制，描述了各種體繪制算法的基本思想，分析比較了相應的關鍵技術，詳細討論了體繪制技術的最新研究進展。

體數據場;體繪制;空間域;變換域

三維體數據處理及可視化的研究工作及應用實驗十分活躍，其技術水平正在從后處理轉向實時交互控制發(fā)展，并且將其與超級計算機、光纖高速網、高性能圖形工作站及虛擬現實等新技術結合起來，成為這一領域技術發(fā)展的重要方向。三維數據場可視化［1，4］主要是指運用交互計算機圖形學和圖像處理及計算機視覺技術，將大量的數據轉換為直觀的圖形和圖像信息在計算機屏幕上顯示出來，供人觀察研究，并進行交互處理的理論、方法和技術，包括體數據的表示、操作和繪制［5］。體數據描述的是三維(可能是時變的)實體及其內部性質，是對某個物理場的離散采樣，沒有明確的表面和邊界。體數據可以通過采樣、仿真或模型技術獲取，例如:由磁共振成像(magnetic resonance imaging，MRI)或計算機斷層攝影(Computed Tomography，CT)得到一系列二維切片數據構成三維體數據。

幾何圖元的繪制方法已經比較完善，早期的體繪制方法都是逼近體數據的一個等值面，然后使用幾何圖元繪制方法對等值面進行繪制的辦法。當體數據的繪制使用面繪制技術時，體數據的內部信息就會丟失，只能繪制表面信息。近年來大量技術應用于體繪制，產生了直接體繪制技術描述全部的三維信息，直接體繪制技術比面繪制描述的信息多，但繪制時間長，計算復雜。為了提高繪制效率，出現了大量優(yōu)化方法和加速技術。

1 體數據場及其繪制方法分類

體數據場是由一些稱為體素(voxel)的采樣點數據(x，y，z，v)構成的集合 V，這些采樣點數據(x，y，z，v)表示在三維位置(x，y，z)的數據一些性質的值為v。根據采樣點取值特點和采樣點位置分布特點可以對體數據場進行分類。

根據采樣點數據的物理特征，可將體數據場分為二值數據場、標量數據場、向量數據場和時變數據場四類。在每個三維空間位置，采樣點v的取值可以是二值、標量、向量或者時變值。若v的取值只取兩個值:0和i，其中0表示背景，i表示對象I，此時稱為二值數據場;若v的可取多個值，用來描述可度量標量性質，如灰度、密度、溫度、壓力等，只有大小而無方向，稱為標量數據場;若v是向量，比如用來描述每個位置的速度，不僅有數值的大小，還有方向的變化，稱為向量數據場;若v是時變的，此時數據場變?yōu)樗木S采樣點數據(x，y，z，t，v)的集合，稱為時變數據場。

根據采樣點空間分布幾何特征［6］，可將體數據場分為結構化數據場和非結構化數據場兩類。非結構化的又稱為散亂的數據場。若數據場在邏輯上是較規(guī)整的，則稱為結構化的。結構化數據場又可進一步分為規(guī)則結構化和非規(guī)則結構化數據。非規(guī)則結構化網格數據也稱曲線(Curvilinear)型網格數據，這類數據中，每一體素是邏輯上的六面體，相對的面并不要求平行，且每一面的四個定點可以不共面。規(guī)則結構化網格數據包含笛卡兒型、等距型、矩形型三種類型，笛卡兒型(均勻網格型):建立在笛卡爾坐標系上的體數據。每一個體素大小相同、也各維比例也完全相同，按照坐標軸方向均勻排列成正方體形狀。等距型的所有體素大小相同，按坐標軸方向排列成長方體。矩形型:沿每一坐標軸，體素間距各不相同，但體素仍是沿坐標軸排列的長方體。一般地，非規(guī)則結構化網格數據和非結構化網格數據可以通過計算幾何方法轉化成規(guī)則結構化體數據，所以只需考慮規(guī)則結構化體數據。

1.1 體數據場的插值

在體繪制中，要繪制的對象本身最原始的時候都應該是連續(xù)的信號(數據)。而體繪制中最初的輸入都是通過采樣得到的三維離散數據場。體繪制算法中為了真實、準確地反映原始數據的特征，在繪制前，首先要將這些離散的數據重構成連續(xù)的三維數據場，然后根據繪制分辨率對重構的連續(xù)信號進行重新采樣，最后選擇某種具體的體繪制算法進行渲染。通常通過插值來完成三維數據場的重構，體繪制圖像的質量很大程度上取決于所用的重構(插值)函數的精確度，另外重采樣過程是制約體繪制運算速度的一個關鍵因素。

1.2 體數據場繪制方法分類

根據數據描述方法和繪制方法的不同，體視化技術可以分為兩大類:面繪制方法(Surface Based Rendering，SBR)［7－8］和直接體繪制方法(DirectVolumeRendering，DVR，簡 稱 體 繪制)［9－16］。隨著研究的不斷深入，每一種技術又包含很多具體方法，如圖1是體視化技術分類圖。面繪制是指體表面的重建，它可以有效地繪制三維體數據的某個表面，但缺乏體數據內部信息的表示;體繪制也稱為直接體繪制，直接由體數據生成三維體的圖像，以體素作為基本單元，能夠表示體數據的內部信息，但相應的計算量大，包括體數據的采樣、重構、重采樣、繪制等操作。

圖1 體視化技術分類圖

2 面繪制

在面繪制中，首先由三維空間數據場構造出中間幾何圖元(多邊形或三角形面片)來逼近等值面，即逼近數據場中滿足某一特性的物質，然后利用已有的算法對等值面進行光照計算和渲染。面繪制方法可以利用現有的圖形硬件對繪制進行加速，使生成圖像及其變換的速度加快，并且可以產生的等值面圖像較清晰。然而面繪制方法生成的圖像只能反映原始三維數據場的局部信息，不能反映其全局信息，生成的圖像也不夠精細，而且該方法需要對數據進行二值分割來生成等值面，分割精度要求較高。目前常見的面繪制方法有以下四種:輪廓線連接法、移動立方體法、移動四面體法和剖分立方體法。

2.1 輪廓線連接法

輪廓線連接法是首先提取每層圖像的輪廓線，然后選取不同圖像層間輪廓線上的頂點，連接這些頂點形成三角面片，最后通過光照計算進行渲染。

輪廓線連接法存在對應問題、鑲嵌問題和分支問題，因此只適用于層間等值面變化較小或者具有幾何形狀上比較相似，且繪制精度要求較低的場合，這些缺點限制了輪廓線連接法的進一步應用。

2.2 移動立方體法

移動立方體方法是最有影響的面繪制重建方法，是進行規(guī)則體數據場等值面生成和繪制的經典算法。由兩層間8個相鄰體素構成的體素內部構造代表等值面的三角面片，首先根據二值分割確定的等值面，然后在每一個體素內逐個比較其頂點值，判斷其位于等值面內或等值面之外。如果體素的一條邊兩個頂點一個在等值面內，一個在等值面外，則表示等值面與該邊相交，求出其交點坐標，將這些交點作為三角面片的頂點。

雖然移動立方體方法的基本算法存在面片過多和二義性(Ambiguity)的缺陷，但經過改進后，現已經廣泛用于各種體數據的可視化。

2.3 移動四面體法

移動四面體法的主要目的是為了解決移動立方體法二義性問題，首先將立方體體素剖分成四面體，然后在四面體內部構造等值面。

由于四面體是構成任何類型的多面體的基本元素，是最簡單的多面體，因此移動四面體法比移動立方體法應用面更廣。此外，移動四面體法比移動立方體法構建的等值面精度更高。

2.4 剖分立方體法

當離散三維數據場的密度很高，接近或者超過計算機屏幕的顯示分辨率時，采用移動立方體法構造的三角面片與屏幕像素大小差不多，甚至更小，此時不必通過插值來計算小三角面片，剖分立方體法正好改進移動立方體法的這個缺點。

剖分立方體法和移動立方體法一樣對數據場中的體素逐一進行處理。當某一個體素8個頂點的體素值中有的在值面內，有的在等值面外，且此體素在屏幕上的投影大于屏幕像素時，則將此體素沿x，y，z三個方向進行剖分，直到體素都在等值面之內或都在等值面之外或者其投影等于或小于像素，將投影等于或小于像素的小體素投影到屏幕上，形成所需要的等值面圖像。

3 直接體繪制

直接體繪制方法不構造中間幾何元素，直接將三維數據場投射到二維屏幕上，能生成三維數據場的整體圖像，包括體數據內部信息，繪制的圖像質量高，能得到數據深度信息。但其計算量巨大，且很難利用傳統(tǒng)的圖形硬件進行加速，導致計算時間較長。然而，直接體繪制近年來發(fā)展迅速。

3.1 空間域體繪制

空間域體繪制方法的實質可概括為三維數據場的重采樣和顏色合成兩大步驟。根據繪制次序，空間域體繪制算法可以分為三類:以圖像空間為序(簡稱像序)的體繪制算法、以對象空間為序(簡稱物序)的體繪制算法、及圖像和對象空間混合序(混合序)的體繪制算法。

(1)像序體繪制算法由屏幕的像素點發(fā)出光線穿過數據場決定像素顏色值。其具有下述兩個共同特點:一是從屏幕上的每一個像素點出發(fā)，根據視點位置，發(fā)射穿過三維體數據場的光線，沿這些光線進行離散化采樣，按照一定原則選取若干重采樣點;二是通過插值近似計算這些重采樣點間的顏色和不透明度融合計算，計算對應屏幕像素的顏色。常見的像序體繪制方法根據采樣方式和模式的不同包括X線繪制、最大強度投影法和等值面繪制三種模式，其中X線繪制算法直接將插值后的重采樣點直接和值作為像素顏色，最大強度投影法則直接使用插值后的重采樣點的最大值作為像素顏色。

(2)物序體繪制算法將體數據映射到二維的圖像屏幕，將對二維屏幕像素的有影響的每一個體素的貢獻計算出來，然后將其合成為像素的顏色來繪制圖像，生成圖像的步驟如下:首先將每一個體素點投影到屏幕空間的坐標，計算對屏幕的影響范圍和貢獻，Splatting算法是一種典型的物序體繪制算法。

(3)混合序體繪制算法結合了像序體繪制和物序體繪制算法兩種算法的優(yōu)點。它通過將三維離散數據場的投影變換分解為Shear變換和Warp變換，實現將三維空間的重采樣簡化為二維平面的重采樣，減少了大量計算量，在不顯著降低圖像質量的前提下，可以在微機上接近實時速度繪制體數據。但由于Shear－Warp算法對二維圖像空間的采樣不可能得到完整、正確的三維信息，是以犧牲圖像質量和準確性為代價來加速繪制的。

3.2 變換域體繪制

變換域體繪制首先將空間體數據變換到變換域，然后在變換域內直接生成投影或借助變換域的信息生成投影，從而進行顯示。常見的變換域包括壓縮變換域、頻域和小波域，因此，對應的變換域體繪制方法有壓縮域體繪制、頻域體繪制和小波域體繪制。

(1)壓縮域體繪制是不需要將壓縮后的數據進行解壓縮，直接在壓縮域進行繪制，不需要裝入全部未壓縮的體數據，繪制結果可以立即顯示。因此對計算機的內存、計算和傳輸要求降低，加快了繪制速度，是將繪制與壓縮進行有機結合的好方法。首先在空間域對體數據用向量量化技術進行壓縮，然后直接對量化塊用一般的空間體繪制方法進行繪制。

(2)頻域體繪制算法是根據Fourier切片定理提出來的，是在頻域中進行體繪制過程的技術。根據Fourier切片定理，在三維數據場經過Fourier變換后相對應的頻域場中，按視線方向抽取一個經過原點的截面，再對這個截面做Fourier逆變換，得到的恰好就是在空域中的圖像平面的投影。頻域體繪制先對三維數據場進行三維快速傅里葉變換，得到數據場的頻域表示，然后對任意投影方向(即觀察方向)，在頻域內通過原點且垂直于投影方向的平面內進行二維切片重采樣，最后對二維采樣得到的數據場進行二維Fourier逆變換，即得到空間域表示的三維數據場沿投影方向的二維平行投影圖像。

(3)小波域體繪制首先對體數據進行三維離散小波變換，將體數據進行多分辨率表示，并進行數據壓縮;在繪制階段，不需要執(zhí)行任何解壓縮操作，直接將小波變換的小波系數代入體繪制方程生成二維圖像。小波變換具有在空間域和頻域同時具有局部性質，可以用少量小波系數表示體數據。

4 優(yōu)化網格數據的體繪制技術

優(yōu)化網格是指以截八面體為體素的體心立方(body－centered cubic)網格和以菱形十二面體為體素的面心立方(face－centered cubic)網格。優(yōu)化網格上的體素分布比立方網格更加緊湊和合理，在三維空間具有更好的采樣性質，有關優(yōu)化網格的直線生成［17］、Fourier變換［18］和立體顯示［19］的研究日益活躍。運用采樣性質更好的面心立方網格和體心立方網格代替?zhèn)鹘y(tǒng)的笛卡爾網格，在不影響質量的前提下，能減少大約30%的數據量，從而提高體繪制速度。因此，與立方網格體視化技術類似，根據數據描述方法和繪制方法的不同，優(yōu)化網格體視化技術可以分為兩大類:面繪制方法和直接體繪制方法，其中直接體繪制方法按處理數據域的不同可分為優(yōu)化網格空間域方法和優(yōu)化網格變換域方法，優(yōu)化網格空間域體繪制可以分為三類:優(yōu)化網格像序體繪制算法、優(yōu)化網格物序體繪制算法、優(yōu)化網格混合序的體繪制算法。

關于優(yōu)化網格面繪制方法，Carr等人［20］提出了體心立方網格上的移動八面體法、改進的移動八面體法和改進的移動六面體法，這些算法都生成體心立方網格體數據場的等值面。體心立方網格可視為由兩個交錯的立方網格構成，分別稱為主網格和次網格。構造采樣網格最簡單的方式為對角連接主網格和次網格點的Delaunay體元，這些Delaunay體元均為四面體，在主網格和次網格上具有公共邊的4個四面體可以組成一個八面體，這些八面體能鋪滿整個三維空間，與移動立方體算法類似，可以得到移動八面體算法，該算法同樣存在二義性的缺陷，但生成的三角面片數目比移動四面體算法少。移動六面體算法則是將具有公共對角邊的四面體組成六面體。Strand等人［21］則提出了面心立方網格上的類似于移動立方體法的等值面生成算法。

優(yōu)化網格像序體繪制算法中，Ibá?ez等人［22］將光線投影算法推廣到體心立方網格和面心立方網格，為了對穿過三維體數據場的光線進行采樣，提出了任何結構化網格上的直線向量遍歷算法逐個生成體素，并提出了四種策略(雙線性插值、三線性插值、重心插值和剪切三線性插值)進行重采樣。

優(yōu)化網格混合序體繪制算法將像序體繪制和物序體繪制算法有機結合起來。Sweeney等人［23］將Shear－Warp算法推廣到體心立方網格。由于體心立方網格可看成交錯的兩個立方網格構成，因此體心立方網格體數據可以非常方便地分解為兩個立方網格體數據，然后分別對這兩個立方網格體數據進行Shear變換和二維圖像的Warp兩步，最后再將其合成，值得注意的是在Z軸方向有半個網格距離的偏移量。

優(yōu)化網格的變換域體繪制算法在最近幾年被相繼提出，Dornhofer在他的博士畢業(yè)論文［24］中建立了一般規(guī)則網格上的離散Fourier變換公式并將其應用于體心立方網格的頻域體繪制。Alim等人［25］給出了面心立方網格和體心立方網格的一種快速Fourier變換算法。文獻［26］和［27］提出了應用體心立方網格Box樣條進行體繪制的算法。

5 體繪制的加速技術

由于體數據場的數據量巨大，體繪制計算復雜度較高，計算較為費時，導致體繪制速度與用戶要求實時、交互進行顯示處理之間的矛盾越來越嚴重。為解決這類矛盾，研究者提出了大量優(yōu)化體繪制的加速技術。這些加速技術主要包括軟件算法優(yōu)化和圖形硬件加速兩類方法。

在軟件算法優(yōu)化方面［28］，主要有兩大類技術可以有效提高體繪制算法速度:第一類是空間數據結構優(yōu)化，這類技術使用八叉樹、k－d樹、索引表、距離變換等數據結構對體數據進行優(yōu)化和編碼，剔除無效體素，達到體繪制加速的目的。第二類稱為提前不透明度截止，這類方法基本思想是:在應用由前向后融合算法的體繪制過程中，當累積的不透明度達到一個預先設置的閾值時，光線可以到達的任何其它體素可以不作處理，因為其將被遮擋。

在圖形硬件加速［29］方面，針對體繪制算法不能很好的借助傳統(tǒng)的圖形硬件加速的問題，研究者開發(fā)設計了專用體繪制體系結構。

針對大規(guī)模和超大規(guī)模體數據，并行和分布式體繪制研究［30－32］也取得了很大進展。根據各處理器間分配的數據類型不同，目前提出的并行體繪制主要可以分為兩大類:一類是基于圖像空間數據劃分的并行算法，比較適合于像序體繪制算法的并行計算。另外一類基于對象空間數據劃分的并行算法，比較適合于物序體繪制算法的并行計算。

6 體繪制研究展望

隨著數據場規(guī)模的不斷膨脹，數據量越來越巨大，體繪制技術不能滿足人們對數據場進行實時、精確、交互操作的要求。因此，提高算法效率是體繪制的發(fā)展方向，具體有以下幾類:一是設計和開發(fā)專用的體繪制相關的硬件加速設備;二是設計分布式并行體繪制算法，以滿足實時交互的要求;三是將硬件加速設備和分布并行算法相結合。

［1］ DEFANTI TA，BROWN MD.Brown.Visualization in Scientific Computing［J］.Advances in Computers，1991，33:247－307.

［2］DEFANTI TA，BROWN MD，MCCORMICK BH.Visualization:Expanding Scientific and Engineering Rendering Opportunities［J］.Computer，1989，22(8):12 －16.

［3］石教英，蔡文立.科學計算可視化算法與理論［M］.北京:科學出版社，1996.

［4］唐澤圣，孫延奎，鄧俊輝.科學計算可視化理論與應用研究進展［J］.清華大學學報:自然科學版，2001，41(4/5):199－202.

［5］MIN C，ARIE K，RONI Y.Volume Graphics，Springer［M］.London，February，2000.

［6］DON S，STEVE K.Volume probes:Interaetive data exploration on arbitrary grids［J］.ACM SIGGRAPH Computer Graphics，1990，24(5):5－12.

［7］ Marc L.Display of Surfaces from Volume Data，IEEE Computer Graphics＆Applications［J］.1988，8(3):29 －37.

［8］LORENSEN WE，CLINE HE.Marching cubes:A high resolution 3D surface construction algorithm［J］.ACM SIGGRAPH Computer Graphics，1987，21(4):163 －169.

［9］WESTOVER L.Footprint evaluation for volume rendering［J］.ACM SIGGRAPH Computer Graphics 1990，24(4):367－376.

［10］LACROUTE P，LEVOY M.Fast volume rendering using a shear－ warp factorization of the viewing transformation［C］.Proceedings of the 21st annual conference on Computer graphics and interactive techniques 1994:451－458.

［11］SWEENEY J，MUELLER K.Shear－ Warp Deluxe:The Shear－ Warp Algorithm Revisited［C］.Proceedings of the symposium on Data Visualisation 2002:95 －105.

［12］ BUNEMAN O.Conversion of FFT＇s to fast Hartley transforms［J］.SIAM Journal on Scientific and Statistical Computing，1986，7(2):624－638.

［13］MALZBENDER T.Fourier volume rendering［J］.ACM Transactions on Graphics，1993，12(3):233 －250.

［14］LEVOY M.Volume rendering using the Fourier projection－slice theorem［C］.In Proceedings of the Graphics Interface ＇92，1992:61 －69.

［15］TOTSUKA T，LEVOY M.Frequency domain volume rendering［J］.Computer Graphics SIGGRAPH ＇93 Proceedings，1993，27(4):271 －278.

［16］GROSS M H，LIPPERT L，DITTRICH R，et al.Two methods for wavelet－ based volume rendering［J］.Computers and Graphics，1997，21(2):237－252.

［17］何麗君，劉勇奎，孫世昶.三維面心立方網格上的直線生成算法［J］.計算機學報，2010，33(12):2407－2416.

［18］姚繼鋒，孫家昶.平行十二面體區(qū)域上的快速離散傅立葉變換及其并行實現［J］.數值計算與計算機應用，2004(4):303－314.

［19］李莉，楊忠，邢建芳，等.面向立體顯示的點采樣柵格優(yōu)化策略及其性能分析［J］.山東大學學報(工學版)，2008，38(3):1－6.

［20］CARR H，THEU?L T，M?LLER T.Isosurfaces on Optimal Regular Samples［C］.Proceedings of Eurographics Visualization Symposium，2003:39－48.

［21］ STRAND R，STELLDINGER P.Topology Preserving Marching Cubes－like Algorithms on the Face－Centered Cubic Grid［C］.Proceedings of the 14th International Conference on Image Analysis and Processing，2007:781－788.

［22］IB??EZ L，HAMITOUCHE C，ROUX C.A Vectorial Algorithm for Tracing Discrete Straight Lines in N － Dimensional Generalized Grids［J］.IEEE Transactions on Visualization and Computer Graphics，2001(7):97 －108.

［23］SWEENEY J，MUELLER K.Shear－ Warp Deluxe:The Shear－ Warp Algorithm Revisited［C］.Proceedings of the symposium on Data Visualisation，2002:95 －105.

［24］DORNHOFER A.A Discrete Fourier Transform Pair for Arbitrary Sampling Geometries with Applications to Frequency Domain Volume Rendering on the Body－Centered Cubic Lattice［D］.Vienna:Vienna University of Technology，2003.

［25］ALIM UR，M?LLER T.A fast Fourier transform with rectangular output on the BCC and FCC lattices［C］.Proceedings of the eighth international conference on sampling theory and applications(SampTA ＇09)，Marseille，2009.

［26］ ENTEZARI A，VAN DE VILLE D，M?LLER T.Practical box splines for volume rendering on the body centered cubic lattice［J］.IEEE Transactions on Visualization and Computer Graphics，2008，14(2):313 －328.

［27］BERNHARD F，ALIREZA E，DIMITRI V，et al.Efficient volume rendering on the body centered cubic lattice using box splines［J］.Computers＆ Graphics，2010，34(4):409－423.

［28］HERNELL F，LJUNG P，YNNERMAN A.Local Ambient Occlusion in Direct Volume Rendering［J］.IEEE Transactions on Visualization and Computer Graphics，2010，16(4):548 － 559.

［29］JENS S，MARTIN K，R?DIGER W.GPU －Based Euclidean Distance Transforms and Their Application to Volume Rendering［J］.Communications in Computer and Information Science，2010，68(4):215－228.

［30］KWAN －LIU M，JAMES S，CHARLES D，et al.Parallel volume rendering using binary－swap image composition［C］.ACM SIGGRAPH ASIA，2008.

［31］PETERKA T，HONGFENG Y，ROSS R，et al.Endto－End Study of Parallel Volume Rendering on the IBM Blue Gene/P［C］.International Conference on Parallel Processing，2009.Vienna，566 － 573.

［32］Jeff AS，Cheng－Kai C，Kwan－Liu M，etal.Multi－GPU volume rendering using MapReduce［C］.Proceedings of the 19th ACM International Symposium on High Performance Distributed Computing，2010.

Recent Advances in Volume Rendering for 3D Data Sets

HE Li－jun，WANG Xiao－qiang，YUN Jian，BAO Shu－zhe
(School of Computer Science＆ Engineering，Dalian Nationalities University，Dalian Liaoning 116605，China)

Volume rendering for 3D data sets，a hot research topic in computer graphics，has important significance in science and worthiness in practical application.Volume data sets are classfied by physical characteristics and geometrical characteristics.The basic ideas and the key techniques of the two kinds of methods of the 3D data sets called surface rendering and direct volume rendering，are introduced into this paper.The recent advances in volume rendering is discussed in detail.

volume data sets;volume rendering;spatial domain;transform domain

TP391

A

1009－315X(2012)05－0486－06

2011－11－16;最后

2012－02－28

遼寧省教育廳科學技術研究項目(L2010096);中央高?；究蒲袠I(yè)務費專項資金資助項目(DC10010114);大連民族學院青年基金(2007A204)。

何麗君(1974－)，女，湖南新邵人，副教授，博士研究生，主要從事計算機圖形學、計算幾何、數據庫技術研究。

(責任編輯 劉敏)