孟世洋，李晶

1.云南大學 物理科學技術學院，云南昆明 650032；2.昆明總醫(yī)院 醫(yī)學工程中心，云南 昆明 650032

一種眼眶容積測量的新方法

孟世洋1，李晶2

1.云南大學 物理科學技術學院，云南昆明 650032；2.昆明總醫(yī)院 醫(yī)學工程中心，云南 昆明 650032

目的 建立一種方便快速，重復性高，精度準的眼眶測量方法，有效指導眼眶手術的整復治療。方法 利用計算機三維影像技術和計算機圖形技術，對成人眼眶特征點進行標記，提取眼眶內表面，結合逆向工程軟件，實現(xiàn)對眼眶容積的測量。采用填砂測量10例正常成人20只眼眶容積作為對照組，驗證眼眶容積測量新方法的準確性。結果 參照對照組，測量結果誤差控制在-2.58%～2.33%，精確度0.15%～2.58%。對同一例成人左、右眼眶測量的方差分別為1.2%和1.4%。結論 這種眼眶容積測量的新方法具有方便快捷、操作具有較強可重復性、精準度高等特點，可以在臨床應用上進行推廣。

眼眶容積；三維影像技術；逆向工程；CT測量；填砂測量

眼眶是位于顱頂骨和顱面骨之間的兩個骨腔，容納并保護眼球及眶內重要結構，形態(tài)不規(guī)則，位置較深。眼眶爆裂骨折、眶內腫物等疾病，均能改變眼眶容積。測量眼眶容積，對于眶內腫物摘除術后眶內填充物填入量的估計以及由于爆裂骨折致眼眶容積擴大造成眼球內陷的修復都具有重要指導意義[1]。

長期以來，眼部病變多依賴普通X線片、超聲等影像學檢查方法，但是這此檢查手段只能對眼眶做二維測量，不能對眼眶容積進行測量[2]。隨著CT技術的發(fā)展成熟，特別是用于眼眶檢查之后，使眼眶疾病的診斷有了突破性進展。CT掃描具有組織密度分辨力高，解剖關系清楚，可清晰顯示眼眶的解剖結構及其與周圍結構的相互關系，能清晰顯示軟組織、骨骼和氣體，成像不受深度限制，CT重建技術可以指導矯形與整復時立體結構的精確修復，而且利用2D或3D CT重建圖像測量眼眶容積準確可行[3]。國內外對利用CT數(shù)據(jù)重建后測量的方法進行了大量的研究和報道[4-5]，目前測量的方法上主要有影像測量和實物填充測量。本文在以往研究的基礎上提出一種更加方便快捷的方法，改變繁瑣的測量方法，大大縮短了眼眶容積的測量時間，及時快速地為眼部的治療提供重要參考。

1 傳統(tǒng)眼眶容積的測量方法

傳統(tǒng)的眼眶容積的測量方法主要有影像測量和實物填充測量。

1.1 影像測量眼眶容積的方法

利用影像測量眼眶容積主要是利用CT對人體眼眶進行掃描，獲取像素矩陣的CT斷層圖像并行三維重建，而后利用編程程序對選定區(qū)域進行掃描計算出容積[6]。主要步驟如下：

（1）將掃描所得的CT數(shù)據(jù)進行三維重建。

（2）空間線段的測量，確定空間任意兩點間的距離。按公式 (1) 計算距離，其中D為A點(X1，Y1，Z1)到B點(X2，Y2，Z2)點的測量距離。

（3）面積和體積計算采用像素累加的方法完成。按公式 (2) 計算面積，其中Li為待測區(qū)域的周長，S為層厚。

眼眶體積按公式 (3)，其中Ai為每一層待測區(qū)域的面積，S為層厚。

同時為達到精確測量減小誤差的目的，通常采用取兩層的平均值進行累加計算，計算公式為：

通過以上方法測量眼眶容積方便簡單，但其存在的主要存在以下不足：① 選擇區(qū)域體積的龐大致使數(shù)據(jù)量大，計算速度會變慢，這導致對硬件的要求較高，無疑增加了投入的成本；② 容積效應是產(chǎn)生測量誤差的另一原因，它使圖像中眼眶邊緣產(chǎn)生重疊。在進行掃描計算時，對眼眶邊緣的計算存在誤差，而且不適當?shù)拇敖M合，可掩蓋其重疊邊緣的密度差，使測量者判斷眼眶界限產(chǎn)生困難，窗寬過窄時測量容積趨小，窗寬過寬時測量容積趨大[7]。根據(jù)國內外對眼眶部位容積精度與窗位窗寬的關系的報道，在窗位200 HU、窗寬2000 HU時，可分清重疊的邊緣[5,8-9]。

1.2 實物填充測量方法

這種方法主要是利用重構好的眼眶并打印出眼眶實物或者在尸體上進行，應用培樂多泥封閉眶腔的眶上裂、眶下裂、視神經(jīng)孔和小的破裂孔。應用磨面細砂充填整個眶腔，并反復搖放使細砂均勻、結實。充填細砂的標準為細砂平面和眶緣平齊，用精密天平稱細砂的重量，每個眼眶充填多次，取平均值。將每個眼眶內充填的細砂重量除以細砂比重，可得到眶腔的容積[5,8]。

但這種方法的缺點如下：① 眼眶中存在較多的細孔，填充不便，尤其是對較深部位的填充，細沙容易漏出；②填充時以整個假體為模型，而我們面對的是實際的病人，即使可進行3D打印后填充，但耗時較多。

基于對傳統(tǒng)眼眶容積測量存在的一些不足，本文提出的方法不需要較高的硬件要求，一般的硬件即可完成，不需要設計算法掃描眼眶邊緣，有效避免了容積效應。可在活體CT掃描數(shù)據(jù)的基礎上進行測量，縮短了測量時間，真正到達為患者減輕痛苦。

2 材料與方法

2.1 眼眶CT數(shù)據(jù)的獲取

采用美國GE（MedicalSystems/LightSpeed16）公司生產(chǎn)的64排CT連續(xù)掃描，獲取原始CT數(shù)據(jù)。掃描條件：尸體取仰臥位，讓頭部與掃描平面垂直，對眼眶部位進行掃描，層厚0.562 mm，電壓120 kV，電流100 MA，矩陣為512×512，一次掃描時間約為15～20 s。CT測量人體器官體積的基礎，是通過對掃描圖像中興趣區(qū)面積的測定。而構成圖像的基本單位的像素大小，是由掃描視野、重建矩陣決定的。像素越小，測量結果越精確。體素是構成圖像的最小體積單位，即一定厚度的像素。厚度越大，體素代表的體積越大，測量產(chǎn)生的誤差越大[10-11]。

2.2 頭顱模型的重建

將獲取的CT數(shù)據(jù)保存為DICOM格式，導入Mimics10.01（交互式醫(yī)學影像控制系統(tǒng)）軟件[12]中，剛才已闡述過，采用窗位200 Hu、窗寬2000 Hu為眼眶最佳的顯示，可清晰的觀察病人眼眶部位，窗寬過寬會使圖像變得模糊不便于區(qū)分骨骼與軟組織。本文僅針對眼眶的重構，選用人體骨頭閾值進行分割，形成頭顱蒙板，采用0.1 mm重建層距，減小重建間距雖增加重建時間，但可以減小容積效應[13]。3D計算后實現(xiàn)眼眶部位的三維重建，見圖1。

圖1 三維重構結果

2.3 確定特征標記點

標記點的定位是關鍵的，能很好的幫助我們進行眼眶的提取，成人眼眶的表面有著清晰的特征，根據(jù)經(jīng)典人體頜面部解剖標志點，確定眼眶的標志點，見圖2：① 淚嵴點：上頜骨鼻突、額骨和淚骨的交匯點；② 外眶距點：顴額縫的外側緣點；③ 眶上緣點：眶上緣中點；④ 眶下緣點：眶下緣中點。

圖2 定位標記點

確定標記點后應用鼠標沿著骨性眶壁描記，確定提取的眼眶位置，并將文件存儲為STL格式，即利用三角網(wǎng)格來表現(xiàn)3D模型。

2.4 眼眶的提取

在前一步的基礎上，導入逆向工程軟件Geomagic Studio 12.0中[14]，利用軟件中的畫筆工具，沿著眼眶骨性標記進行描記，界定眼眶外側面積范圍，見圖3。然后利用選取工具中的套索工具選擇眼眶的整個內表面，見圖4。隨后翻轉選區(qū)，刪除多余部分，得到眼眶的內表面，見圖5。

圖3 沿標記點描記眼眶外緣

圖4 選取眼眶內表面

圖5 眼眶內表面

2.5 填充

填充是最重要的一個環(huán)節(jié)，首先要確定所填充孔的獨立性，就是孔與孔之間不存在相交，否則容易導致多個孔被填充，從而影響精度。在填充的過程中，要合理選擇填充的方式。Geomagic Studio 12.0軟件提供了填充完整孔、填充邊界孔和生成橋填充3種填充方式，其中完全填充孔適用于完全封閉邊界線形成的孔，填充邊界孔模式適用于填充部分孔，包括邊界缺口或圓周孔的一部分，生成橋填充模式通過生成跨越孔的橋，從而將長窄孔分割成多個孔，并分別填充。對孔的填充要充分結合集中填充模式進行。填充前、后的結果，見圖6。

2.6 容積的計算

將封裝好的眼眶投影到平面，計算出體積。沿著眼眶外緣建立投影平面，使眼眶整個軸心平行于全局坐標z軸，見圖7。

圖6 填充效果

圖7 計算體積

3 結果

3.1 重復性檢測

對1例32歲男性尸體標本的左、右眼眶CT數(shù)據(jù)進行重復操作5次，眼眶容積測量結果：左眼依次為23.64、23.83、23.74、23.51、23.70 cm3；右眼依次為23.88、23.79、23.61、23.91、23.75 cm3。由測量結果計算，左眼的平均值為23.68 cm3，右眼為23.79 cm3，方差分別為0.012和0.014。從計算結果看每次測量的數(shù)據(jù)波動較小，從而保證了測量結果的可重復性。

3.2 體積精度的驗證

實物填充方法雖然耗時較多，不能對病人直接填充，但目前眼眶容積測量方法中其精度較高，常作為其他研究方法的參考標準。為驗證我們提出的方法的正確性，選取10具成年人尸體，男性6具，女性4具，年齡21～45歲，平均年齡28.5歲。以實物填充20只眼眶測量的體積作為標準組，與文中提出的新方法進行比較，見表1。測量結果誤差在-2.58%～2.33%，精確度0.15%～2.58%。

表1 新方法測量誤差與精度值

4 討論

本文提出的新方法能有效方便地指導醫(yī)生進行術前的規(guī)劃，可重復操作。相對于實物填充，省去了細沙填充和眼眶實物的模型打印步驟。相對于影像測量方法，有效地避免了繁瑣的程序設計和不精準的掃描算法。并且可通過多次練習熟練掌握軟件的使用，來減少耗時。本文方法誤差來源主要是在對眼眶的填充步驟，能否保證眼眶的骨性特征是關鍵點。首先是要對眼眶的結構熟悉，才能在填充的過程中保證眼眶的外形，其次是要合理利用軟件提供的3種填充模式，做到填充精細化。綜上所述，本文提出的眼眶容積測量的新方法具有簡單快速、重復性高、精度準的特點，可以在臨床應用上進行推廣。

[1] Koppel DA,Foy RH,Mccaul JA,et al.The reliability of “Analyze”software in measuring orbital volume utilizing CT-derived data[J]. J Cranio Maxill Surg,2003,31(2):88-91.

[2] Bullock P,Dunaway D,Mcgurk L,et al.Integration of image guidance and rapid prototyping technology in craniofacial surgery[J].Int J Oral Max Surg,2013,42:970-973.

[3] Bettschart C,Kruse A,Matthews F,et al.Point-to-point registration with mandibulo-maxillary splint in openand closed jawposition. Evaluation of registration accuracy for computer-aided surgery of themandible[J].J Cranio Maxill Surg,2012,40:592-598.

[4] Christoph Z,Gellrich NC,Metzger MC,et al.Computer-assisted Reconstruction of orbital foor based on cone beam tomography[J]. Brit J Oral Max Surg,2007,45:79-80.

[5] Molina DG,Cebrian JL,Carretero MP.Facial symmetry assessment following computer-aided 3D zygomatico-orbital Reconstruction[J].Magn Reson Imaging,2013,15:123-128.

[6] Seong-hakim S,Dae-seungkim D,Kyung-hoehuh K.Direct and continuous localization of anatomical landmarks for[J]. Compute Aided Surg,2013,116:402-410.

[7] Kolk A,Christoph P,Wiener E,et al.Isotropic proton-densityweighted high-resolution MRI for volume measurement of reconstructed orbital fractures-acomparison with multislice CT[J].Magn Reson Imaging,2008,26:1167-1174.

[8] Ko EC,Chen PT,Tai IH,et al.Fronto-facial monobloc distraction in syndromic craniosynostosis.Three- dimensional evaluation of treatmen toutcome and facial growth[J].Int J Oral Max Surg, 2012,41:20-27.

[9] Schmelzeisen R,Gellrich NC,Schoen R,et al.Navigation-aided reconstruction of medial orbital wall and f oor contour in craniomaxillofacial reconstruction[J].Injury-Int J Care Injured,2004, 35(2):955-962.

[10] Modabber A,R?sch M,Ghassemi M,et al.Non-invasive threedimensional evaluation of periorbital asymmetry in isolated unilateral orbital f oor fractures[J].J Oral Maxil Surg,2014,118: 128-152.

[11] Erik N,van Bezooijen JS,Koudstaal MJ,et al.Orbital change following Le Fort III advancement in syndromic craniosynostosis:Quantitative evaluation of orbital volume,infra-orbital rim and globe position[J]. J Craniomaxillofac Surg,2012,40(3):223–228.

[12] Adolphs N,Haberl EJ,Liu W,et al.Virtual planning for craniomaxillofacial surgery-7 Years of experience[J].J Craniomaxillofac Surg,2014,42(5):e289-e295.

[13] Sicurezza E,Giuseppe P,Leonardia R,et al.Three-dimensional computerized tomographic orbital volume and aperture width evaluation[J].Oral Surg Oral Med Oral Pathol Oral Radiol,2011, 111(4):503-507.

[14] Biao L,Lei Z,Hao S,et al.A novel method of computer aided orthognathic surgery using individual CAD/CAM templates:a combination of osteotomy and repositioning guides[J].Brit J Oral Maxillofac Surg,2013,51:239-244.

A New Method of Orbital Volume Measurement

MENG Shi-yang1, LI Jing2
1. College of Physics Science Technique, Yunnan University, Kunming Yunnan 650032, China; 2. Department of Medical Engineering, Kunming General Hospital, Kunming Yunnan 650032, China

Objective To establish a convenient, quick, repetitive and accurate method of orbital volume measurement so as to provide effective guidance for plastic and reconstructive treatment of orbital surgeries. Methods The 3D imaging and graphics technique was used to mark adult orbital characteristics and extract orbital surface so as to realize the measurement of orbital volume in combination with reverse engineering software. The sand filling measurement was taken to measure 20 orbital volumes of 10 formal adults as Control Group and to verify the accuracy of the new method. Results In contrast with Control Group, error control of the measurement results ranged from -2.58% to 2.33% and the degree of precision was 0.15% to 2.58%. In the same case of one adult, the error variance of the left and right eye was 1.2% and 1.4%, respectively. Conclusion The method of orbital volume measurement in this paper is fast and convenient with strong repetitive operation and high precision, which should be widely promoted clinically.

orbital volume; three-dimensional imaging technique; reverse engineering; CT measurement; sand f lling measurement

R777.5

A

10.3969/j.issn.1674-1633.2016.10.012

1674-1633(2016)10-0040-04

2014-11-12

2014-12-09

李晶，高級工程師。

通訊作者郵箱：li_jingkm@126.com