邵玶等

[摘要]目的：研究根據(jù)不同的患者，不同的牙列情況，建立基于螺旋CT 掃描圖像的上頜牙列的三維重建模型和三維有限元模型，為上頜牙列的正畸矯治提供一定的理論指導(dǎo)。方法：選擇男女兩名自愿者作為建模素材，通過螺旋CT掃描技術(shù)得到頭部二維CT圖像，應(yīng)用三維影像技術(shù)取出含上頜牙列部分圖像，對截取圖像重建三維模型，最后在ANSYS中建立上頜牙列的三維有限元模型。結(jié)果：獲得了細(xì)致逼真的上頜牙列三維重建影像和三維有限元模型。結(jié)論：綜合運(yùn)用CT掃描、三維建模等技術(shù)可以獲得不同個(gè)體的上頜牙列三維有限元模型，為臨床上頜牙列的正畸矯治提供一定的理論依據(jù)。

[關(guān)鍵詞]CT圖像； 三維重建；實(shí)體模型；有限元模型

[中圖分類號]R783.5[文獻(xiàn)標(biāo)識碼]A[文章編號]1008-6455（2012）10-0-0

對人體器官的螺旋CT圖像的實(shí)體建模屬于三維重建技術(shù)領(lǐng)域，是當(dāng)前醫(yī)學(xué)研究的一個(gè)熱點(diǎn)課題。國內(nèi)外對CT圖像的三維重建進(jìn)行了大量的研究，很多學(xué)者采用不同的方法對頭顱[1]、肱骨、下頜骨、足跖骨等進(jìn)行三維重建，用于醫(yī)學(xué)研究。1979年，Artzy和Herman等首先利用CT掃描的二維圖像進(jìn)行了三維重建工作，克服了以往只能從CT掃描的二維圖像中重疊相加進(jìn)行思維綜合的缺點(diǎn)。因此，利用高分辨率CT掃描圖像在計(jì)算機(jī)上進(jìn)行三維重建迅速得到了學(xué)者們的廣泛關(guān)注。尤其到了20世紀(jì)90年代，螺旋CT作為一種新的工具出現(xiàn)后，因具有快速掃描成像、容積數(shù)據(jù)采集、優(yōu)良的多軸面和三維重建圖像等優(yōu)點(diǎn)逐漸取代了傳統(tǒng)CT，使基于螺旋CT掃描圖像的計(jì)算機(jī)三維重建更加精確[2]。

對三維物體的描述方法目前應(yīng)用較多的有：表面重建法、體素重建法以及截面重建法等。表面重建法采用多邊形對物體的表面進(jìn)行擬合，忽略了物體內(nèi)部的信息，不利于進(jìn)行一些分析處理或者實(shí)體制造。體素重建法不需要確立表面的幾何顯示，直接在體數(shù)據(jù)的基礎(chǔ)上進(jìn)行顯示，速度慢，設(shè)備要求比較高，結(jié)果數(shù)據(jù)文件比較大，不利于遠(yuǎn)程傳輸。而截面重建法的重建結(jié)果圖像只具備深度信息，不具備形狀信息，旋轉(zhuǎn)顯示的可視側(cè)面效果差，立體感不明顯，重建時(shí)不能任意旋轉(zhuǎn)不同側(cè)面，重建靈活性差。

本文擬利用螺旋CT圖像，研究面向?qū)嶓w模型的上頜牙列的三維模型建模方法，該方法的基本思想是：選取2個(gè)正在進(jìn)行正畸治療的患者，對二者頭部拍攝CT圖像，對CT斷層圖像進(jìn)行填充，對填充后圖像進(jìn)行濾波，修整，然后三維重建實(shí)體模型，并根據(jù)重建的實(shí)體模型建立三維有限元模型。為后續(xù)的應(yīng)力分析提供實(shí)驗(yàn)基礎(chǔ)。

1材料和方法

1.1 建模素材與實(shí)驗(yàn)設(shè)備：從臨床病例中選取男女兩位成年志愿者，男性志愿者A是牙合平面偏斜患者，治療前狀態(tài)，女性志愿者B是上頜前突患者，上頜左、右第四顆牙齒已拔除，治療中狀態(tài)，牙列有托槽。并且無明顯牙周疾患及牙槽骨吸收，無牙體缺損及明顯磨耗。

硬件64排螺旋CT機(jī)（美國GE公司），軟件 Mimics10.1（比利時(shí)Materalise公司），Ansys12.1 （美國Ansys公司）。

1.2 螺旋CT二維圖像的采集[3-4]：患者取仰臥位，下頜處于息止頜位，避免上、下牙列接觸重疊。采用美國GE公司64排LightSpeed VCT螺旋掃描CT機(jī)水平位薄層掃描被檢者頭部，掃描層厚 0.625mm，數(shù)字矩陣512×512，窗寬：2000，窗位：100。

1.3 CT 圖像的傳輸：由于螺旋CT 掃描機(jī)不具備直接圖像文件數(shù)據(jù)的傳輸功能，因此，將螺旋CT掃描圖像以DICOM格式存儲在內(nèi)置計(jì)算機(jī)上，然后通過光盤刻錄的方式導(dǎo)出至后續(xù)處理用微型計(jì)算機(jī)中，圖1是其中一層的CT掃描橫斷面圖像。

1.4 CT 圖像的三維重建：CT圖像中，不同的組織，就顯現(xiàn)出不同的灰度值。由于骨密質(zhì)，骨松質(zhì)的存在，使得所得到的CT圖像邊緣極不連貫，并且中央圖像灰度值分布也很不規(guī)律，存在大量的空穴，在這種情況下，如果不進(jìn)行任何處理就進(jìn)行三維重建，所得的三維模型含有眾多不規(guī)則孔洞，對于醫(yī)學(xué)圖像的復(fù)雜性，如果要利用程序全自動的解決這個(gè)問題，將會是十分困難的。因此，本文采用了人工填充的方法。對經(jīng)過上述處理后的橫斷面圖像進(jìn)行牙列的三維面重建。圖2為人工填充后重建的2個(gè)患者上頜牙列的三維實(shí)體模型。

僅僅經(jīng)過填充就進(jìn)行重建的三維模型，并不能直接用以建立三維有限元模型。由于醫(yī)學(xué)圖像數(shù)據(jù)的獲取過程中，影響設(shè)備中各電子器件的隨機(jī)擾動不可避免地會產(chǎn)生噪聲，因此需要對獲得的橫斷面圖像進(jìn)行濾波平滑，抑制噪聲，增強(qiáng)圖像特征，提高信噪比。尤其對于牙列上粘有托槽的B患者，其CT圖像有強(qiáng)烈的偽影干擾，必須通過處理去除偽影。圖3是優(yōu)化處理后得到的三維實(shí)體模型。

1.5 三維有限元模型：由于所研究的對象為復(fù)雜形狀實(shí)體模型，為保證所建有限元模型準(zhǔn)確，選擇ANSYS有限元軟件進(jìn)行仿真分析。將上頜牙列三維實(shí)體模型導(dǎo)入ANSYS。選取合適的空間單元，并自由網(wǎng)格劃分，保存上頜牙列網(wǎng)格模型文件并導(dǎo)入Mimics 10.1。在Mimics10.1的FEA模塊中，利用軟件自帶的“CT 圖像灰度-組織力學(xué)材料性質(zhì)”關(guān)系公式，對三維上頜牙列網(wǎng)格模型的各單元進(jìn)行材質(zhì)分配。本研究將模型中各種材料和組織考慮為各向同性的線彈性材料，材質(zhì)分配完成后的體網(wǎng)格模型再次返回導(dǎo)入ANSYS12.1。至此得出完整的上頜牙列有限元網(wǎng)格模型。

2結(jié)果

獲得了兩個(gè)患者上頜牙列的螺旋CT三維重建實(shí)體模型，模型細(xì)致逼真、可被任意旋轉(zhuǎn)、剖開、透視、截取，從不同角度觀察重建影像形態(tài)還原性良好、重建效果滿意?；趯?shí)體模型有限元模型網(wǎng)格均勻，能夠滿足后續(xù)力學(xué)分析。模型圖見圖3、圖4。

3討論

3.1 實(shí)體模型及有限元模型：能夠獲得可用于有限元力學(xué)分析的三維重建模型，一直是致力于有限元理論分析的國內(nèi)外學(xué)者感興趣的課題。早先采用的切片、磨片法因組織丟失、誤差大且不適合活體逐漸被淘汰[5]。上世紀(jì)90年代，國外知名公司如西門子公司、美國通用公司等都相繼開發(fā)了自己的具有三維重建功能的螺旋CT機(jī)。正是基于螺旋CT 技術(shù)和計(jì)算機(jī)技術(shù)的發(fā)展，才使三維重建技術(shù)在活體有限元模型建立并獲得理想、精確的生物力學(xué)研究結(jié)果中起到了重要的作用。

目前，醫(yī)學(xué)圖像三維重建的基本步驟包括，從外設(shè)獲得斷層數(shù)字（灰度）圖像序列；對斷層圖像進(jìn)行分割，對斷層間進(jìn)行插值以生成立方體體素?cái)?shù)據(jù)；進(jìn)行立體感顯示[6-10]。隨著今后影像學(xué)和計(jì)算機(jī)技術(shù)的進(jìn)一步發(fā)展，三維重建技術(shù)必然在口腔醫(yī)學(xué)領(lǐng)域中具有更加廣闊的應(yīng)用前景。

以往建立的醫(yī)學(xué)有限元模型由于軟件功能單一、網(wǎng)格劃分算法少、人體組織結(jié)構(gòu)形狀不規(guī)則等原因，需要進(jìn)行適當(dāng)?shù)暮喕图僭O(shè)，在精度方面有一定的缺陷，而且模型創(chuàng)建過程環(huán)節(jié)多、成本高、周期長，進(jìn)一步限制了其應(yīng)用。以Mimics10.1建立實(shí)體模型并結(jié)合ansys12.1建立三維有限元模型，簡化了建模中提取輪廓線的過程，縮短了建模時(shí)間；再進(jìn)一步的自動重劃功能中，能進(jìn)行更加專業(yè)地半自動或手動劃分；通過掃描數(shù)據(jù)的亨氏單位，對體網(wǎng)格進(jìn)行材質(zhì)分配，從而能顯著提高了有限元模型分析（FEA）結(jié)果的可靠性和精確性，在整個(gè)重建過程中最大限度地排除了人為的干擾，體現(xiàn)了便捷、準(zhǔn)確的建模優(yōu)勢。

3.2 真實(shí)牙列：實(shí)驗(yàn)采用某自愿者的真實(shí)牙列做實(shí)驗(yàn)，能夠盡可能還原患者矯治過程，體現(xiàn)臨床矯治特點(diǎn)。以往人們采用仿真牙列進(jìn)行實(shí)驗(yàn)時(shí)，所有性能參數(shù)都是固定的，重建的三維模型也是相同的，力學(xué)仿真試驗(yàn)結(jié)果有很大相似性。采用真實(shí)患者牙列實(shí)驗(yàn)時(shí)，由于每個(gè)患者的骨質(zhì)密度不同，牙根形狀不同，對正畸力的反應(yīng)不同，同一種矯治方案，對不同患者的矯治效果和矯治周期都會不同。因此，針對不同患者通過試驗(yàn)可以指定個(gè)性化矯治方案。

3.3 托槽影響：拍攝CT時(shí)，托槽能對正常CT圖像形成強(qiáng)烈的干擾，在重建實(shí)體模型之前必須經(jīng)過細(xì)致處理，否則重建的實(shí)體模型無法輸入ANSYS進(jìn)行有限元建模。

4結(jié)論

本文利用64排螺旋CT技術(shù)，結(jié)合Mimics和ANSYS軟件建立了上頜牙列的實(shí)體模型和三維有限元模型。實(shí)體模型形態(tài)細(xì)致逼真，模型形態(tài)與實(shí)體標(biāo)本一致，可以被任意旋轉(zhuǎn)并從不同角度觀察，能得到上頜骨及牙列骨骼鮮明、直觀、整體的印象。模型可編輯性強(qiáng)，可進(jìn)行任意分割、復(fù)制和存儲。上頜牙列的三維有限元模型具有非常精確的組織結(jié)構(gòu)信息，為下一步的模型應(yīng)力分析打下了良好的基礎(chǔ)。

[參考文獻(xiàn)]

[1]Heman GT，Liu HK.Three-dimensional Display of Human Organs from CT[J].Comtp Graph Image Processings，1979，（9）：1-21.

[2]Szolar DH，Groell R，Braun H，et al.UI trafast computed tomography and three-dimensional image processing of CT sialography in patients with parotid masses poorly defined by magnetic resonance imaging[J].Acta Otolaryngol，1996，116（1）：112-119.

[3]胡凱，榮起國，方競，等.人顳下頜關(guān)節(jié)CT 三維重建及其有限元實(shí)體建模[J].中國醫(yī)學(xué)影像學(xué)雜志， 1999，7（2）：137-139.

[4]楊星，馬彪，蘇勤，等.螺旋CT 三維重建的方法及技巧[J].中華放射學(xué)雜志，1999，33（7）：492-493.

[5]王橋，馮曉剛，陳文歸.下頜56 缺失兩基牙雙端固定橋修復(fù)的生物力學(xué)分析[J].中國醫(yī)學(xué)物理學(xué)雜志， 2008，25（2）：678-679.

[6]秦緒佳，歐宗瑛，侯建華. 醫(yī)學(xué)圖像重建模型的剖切與立體視窗剪裁[J]. 計(jì)算機(jī)輔助設(shè)計(jì)與圖形學(xué)學(xué)報(bào)， 2002，14（3）：275-279

[7]張智勇，周洪，劉明智，等.T形曲作用下上頜三維有限元應(yīng)力比較分析[J].中國美容醫(yī)學(xué)， 2007，16（5）： 677-679.

[8]王彬，談龍，葉湘玉，等.上頜切牙段三維有限元模型的建立[J].中國美容醫(yī)學(xué)，2003，12（2）： 138-139.

[9]張曉慧，丁寅，劉東，等.顱上頜復(fù)合體三維有限元模型的建立[J].中國美容醫(yī)學(xué)，2005，14（2）： 194-196.

[10]張晗，邵玶，閆偉軍.基于MEAW的牙合平面左右向偏斜治療的三維有限元研究[J].中國美容醫(yī)學(xué)，2011， 20（8）：1282-1284.

[收稿日期]2012-07-05 [修回日期]2012-08-30

編輯/張惠娟