陳更新 劉崇偉 王錦輝 殷 焱 皋 源 聞大翔

·論著·

利用3D打印技術(shù)制作靜脈血管模型的研究

陳更新 劉崇偉 王錦輝*殷 焱 皋 源 聞大翔

3D打印技術(shù)是目前科研和商業(yè)領域的熱點之一，可直接由CAD 數(shù)據(jù)快速精確的制造出復雜幾何形狀的產(chǎn)品。在醫(yī)學上主要應用于醫(yī)學模型、骨骼、器管和細胞打印。目前CT技術(shù)雖然可以通過軟件實現(xiàn)三維成像，但難以用于手術(shù)模擬。3D打印技術(shù)的迅猛發(fā)展使得實體醫(yī)學模型的獲得成為可能。本文提出了基于CT掃描圖片通過3D打印技術(shù)制作透明靜脈血管模型的方法。該方法主要依據(jù)患者的CT切片圖，通過MATLAB、CAD等軟件建立血管的三維點云以及立體模型，并成功用3D打印方法制造出兩種不同精密度的血管模型，還原真實的鎖骨下靜脈、頸內(nèi)靜脈以及頭臂靜脈的人體血管結(jié)構(gòu)。此模型建立為鎖骨下靜脈穿刺置管術(shù)的手術(shù)模擬和改進提供了便利，同時該方法亦適用其他人體血管3D模型的制作。

3D打印 三維設計 血管模型

3D打印（3D printing）是一種以數(shù)字模型文件為基礎，運用粉末狀金屬或塑料等可粘合材料，通過逐層打印的方式來構(gòu)造物體的技術(shù)?？捎蒀AD數(shù)據(jù)快速精確的制造出復雜幾何形狀的產(chǎn)品，真正實現(xiàn)無模制造，具有成型速度快、成型精度高、表面質(zhì)量好，處理簡單省時等優(yōu)點。目前大量應用至航天、建筑、汽車、醫(yī)療等很多領域［1］。在醫(yī)學上主要應用于醫(yī)學模型、骨骼、器管和細胞打?。?-5］。本文提出的3D打印血管的方法是通用的，無論是何種血管，只要可以獲取一組清晰的增強CT圖，即可以按照該方法打印出仿真血管模型，對于疾病診斷與治療有很大的幫助。

1 數(shù)據(jù)的采集

圖1 CT血管切片圖

自患者的肘靜脈或手背靜脈注入含碘造影劑后進行頸胸部增強CT掃描，從而使頸內(nèi)靜脈、鎖骨下靜脈及頭臂靜脈增強顯影，從而得到增強CT血管掃描切片164張，同時利用CT機自帶的后處理工作站軟件，除去了血管周圍的組織，得到了獨立的血管切片圖（見圖1）。其中白色為血管，血管周圍的黑色密度結(jié)構(gòu)為脂肪組織。此CT圖片采集自仁濟醫(yī)院放射 科2015年10月數(shù)據(jù)（患者女，38歲，身高165cm，體重68kg）。另外通過CT測量得到了167例患者的靜脈血管數(shù)據(jù)，根據(jù)不同年齡段分組計算平均值，數(shù)據(jù)見表1。

表1 不同年齡段靜脈血管參數(shù)

2 3D模型的建立

2.1 簡化3D打印模型 參考表1數(shù)據(jù)作者利用UniGraphics建模軟件［6］制作了一個正常靜脈角模型（靜脈角75°）和與其對比的一個極端靜脈角模型（靜脈角110°）。正常靜脈角模型和極端靜脈角模型中右鎖骨下靜脈直徑為10mm，右頸內(nèi)靜脈直徑為14mm，右頭臂靜脈直徑為12.5mm。模型建立步驟：依照表中兩組數(shù)據(jù)建立血管中心線的方程，由于血管曲率較小，此處以二次方程模擬。中心線方程為三部分，代表3條血管：鎖骨下靜脈、頸內(nèi)靜脈和頭臂靜脈。根據(jù)靜脈血管的形狀，取三條血管在同一個平面上，開口分別向下、向左和向右；并將血管稍作旋轉(zhuǎn)，使得頸內(nèi)靜脈和頭臂靜脈交點處兩條曲線共切線。在三中心線交點處填補圓柱之間的空隙，得到血管模型。打印出來的普通患者和肥胖患者鎖骨下靜脈穿刺血管實物模型見圖2。

圖2 靜脈血管3D打印模型。A靜脈角為75°；B靜脈角為110°

2.2 仿真血管模型 為更加逼真地反應血管的形狀，作者根據(jù)增強CT切片圖構(gòu)建了仿真血管模型。模型建立步驟：先將每一張CT圖中血管進行一系列處理，提取邊界，轉(zhuǎn)化為一系列的散點；根據(jù)每張CT切片的間距將每組散點組合放置于同一坐標系中，在空間中形成點云［7］；再利用計算機輔助設計軟件將點云重新組合成一個新的柱面；將柱面加厚，得到仿真血管模型。從切片圖到散點是圖片數(shù)字化處理過程，要取出有效的點坐標，也就是點云；從點云到立體圖是一個逆向工程。從CT圖片到散點的過程中經(jīng)歷步驟見圖3。以其中的一張CT片（如圖1）為例，相關變化見圖4。原始CT片（見圖1）周圍有相關數(shù)據(jù)記錄，所以先將其剪切掉排除干擾；此時血管為白色（圖4A）。黑白互換得到反相圖（圖4B），因為黑色圖像處理更加方便；此時血管為黑色。純化目標是遍歷每一個點，將其中非黑色的點全部變成白色，使黑色的血管更加突出（見圖4C）。然后利用MATLAB函數(shù)庫中的算法［8］，獲取黑色血管的邊界（見圖4D），此時圖片自動被黑白反相所以需要再次反相得到黑色的血管輪廓（見圖4E）。從散點至點云過程利用了MATLAB軟件的scatter3d函數(shù)，建立點云見圖5。獲得點云后，利用中望CAD+計算機輔助設計軟件［9］，根據(jù)點云數(shù)據(jù)重新構(gòu)造相應血管的完整曲面。利用“蒙皮”功能，其原理是利用一個較大的柱面不斷收縮，直至與散點接觸，從而得到了散點的外圍形狀。獲得面之后，利用Powershape軟件對面進行“潤化”處理，即使模型的面與面之間過度更加平滑。之后使用“加厚”功能，此處認為血管的厚度是均勻的，將柱面向外平行擴展，將無厚度的柱面轉(zhuǎn)化為有厚度的空心柱體，最終的三維血管模型見圖6。得到血管模型以后，將文件以.stl格式輸出至3D打印機，即可打印出3D打印實物，見圖7。打印模型強度較好，光透明性很好，經(jīng)人工拋光之后可以達到完全透明。

圖3 CT圖片處理流程圖

圖4 從CT圖片到散點設計過程：A剪取CT片中間部分；B反相圖（黑白互換）；C純化目標；D獲取邊界；E再次黑白反相

圖5 點云的三維顯示

圖6 靜脈血管仿真模型

圖7 3D打印血管模型。A正視圖、B側(cè)視圖、C俯視圖

3 討論

3D打印主要分為三個步驟：三維設計、切片處理和打印。三維設計是指通過計算機建模軟件建立三維模型。3D打印模型的標準文件格式是STL文件，使用三角面來近似模擬物體表面。三角面越小，其生成的分辨率越高。切片處理和打印兩個步驟是利用現(xiàn)有的3D打印機完成的。

以本文所用激光SPS600B快速成型機為例，其所用打印材料為透明光敏樹脂。首先打印機讀入STL文件，對模型進行從下向上切片，將三維模型分解成多個二維截面。根據(jù)二維截面輪廓要求，在計算機控制下，沿樹脂槽中盛滿液態(tài)光敏樹脂的液面進行掃描，并用紫外激光快速固化掃描區(qū)域，從而得到該截面輪廓的樹脂薄片，然后工作臺下降一層切片的厚度，已經(jīng)固化的樹脂薄片被一層新的液態(tài)樹脂所覆蓋，再進行掃描固化，新的固化層與前一層粘結(jié)在一起。這樣的過程重復進行直至樣品完整成型。最后取出樣品進行清洗，二次固化及表面處理即可。

鎖骨下靜脈置管術(shù)能夠直接快速有效的輸入大量液體進入血液循環(huán)，尤其在手術(shù)室危重患者的輸液、搶救及中心靜脈壓檢測等方面起到了極其關鍵的作用，因而目前得到廣泛的應用［10-11］。由于右鎖骨下靜脈置管術(shù)手術(shù)風險相比左鎖骨下靜脈置管術(shù)較小，使用更多［12］。但是穿刺異位率較高［13］。作者提出的磁鐵引導法鎖骨下靜脈置管術(shù)可以明顯降低穿刺異位率［14］。利用3D打印制作透明可視的仿真鎖骨下穿刺靜脈血管模型對探究磁引導鎖骨下靜脈穿刺術(shù)提供了較大便利，可以直觀地觀察導絲的走向。另外個性化3D打印的靜脈模型亦可用于鎖骨下靜脈穿刺置管術(shù)的手術(shù)模擬［15］。

鎖骨下靜脈穿刺血管模型設計分為兩種，分別是不同靜脈角的簡單血管模型設計和仿真血管模型設計。簡單血管模型設計是將血管當作光滑的圓柱體，模型制作方便，獲得了血管的基本數(shù)據(jù)后，使用普通三維設計軟件就可以完成，可以根據(jù)需要設定靜脈角；而仿真設計是依據(jù)164張CT掃描圖片，還原度極高，可以真實地表現(xiàn)出血管的外觀形態(tài)如血管表面的凸起、凹陷和細微的內(nèi)徑變化等。

從表1 中可以看出不同年齡段靜脈血管參數(shù)有所不同。右鎖骨下靜脈直徑平均值介于9.4~10.6 mm之間；右頸內(nèi)靜脈直徑平均值介于12.5~14.2mm之間；而右頭臂靜脈直徑平均值介于12.3~12.9mm之間。靜脈角平均值介于73.3°~79.8°之間。表1結(jié)果與靜脈角及附近腔內(nèi)結(jié)構(gòu)的解剖數(shù)據(jù)基本一致［9］。

圖2為靜脈角分別為75°的“正常”靜脈角實物模型和110°極端靜脈角實物模型。作者將此實物模型右鎖骨下靜脈處和透明人體模型膠結(jié)在一起，血管放置在人體模型上表面下方3.5cm處，水平距離鎖骨1.5cm。鎖骨下靜脈在鎖骨中點對應位置灼穿一個小孔，作為穿刺點。人體模型在穿刺點上方對應位置也灼穿小孔。小孔邊緣打磨光滑，導絲伸入小孔內(nèi)運動不會受到阻礙。將導絲經(jīng)鎖骨下靜脈緩慢插入，可以清楚地看到導絲在靜脈匯合處的走向。對于靜脈角為75°的正常靜脈角實物模型，導絲都基本向下進入了頭臂靜脈，而對于靜脈角為110°的極端靜脈角模型，導絲向上進入頸內(nèi)靜脈，發(fā)生穿刺異位。因此3D打印的靜脈角模型可以深入地研究靜脈穿刺情況。這對手術(shù)模擬或?qū)︽i骨下靜脈置管術(shù)的培訓均有所幫助。圖7為精密度更高的靜脈角實物模型。從三視圖上可以清晰地觀察到鎖骨下靜脈、頸內(nèi)靜脈和頭臂靜脈的表面形態(tài)，及其尺寸的變化。

綜上所述，作者通過CT測量及增強CT切片圖得到頸內(nèi)靜脈，鎖骨下靜脈及頭臂靜脈血管的相關數(shù)據(jù)，利用三維建模軟件設計了靜脈角分別為75°和110°的光滑血管模型以及完全根據(jù)CT片切片圖得到的仿真血管模型，并利用3D打印機制作了透明的血管實物模型。

［1］ 王康, 黃筱調(diào), 袁鴻. 3 D打印技術(shù)最新進展. 機械設計與制造工程, 2015, 44(10): 1-6.

［2］ 張云峰, 楊棟. 3D打印骨科模型臨床應用的初步探索. 中國臨床研究, 2014, 27(10): 1260-1261.

［3］ 許世兵, 單樂天, 金紅婷, 等. 3D打印技術(shù)在骨科的應用進展.中華創(chuàng)傷雜志, 2015, 31(1): 10-15.

［4］ 王鎵垠, 柴磊, 劉利彪. 人體器官3D打印的最新進展. 機械工程學報, 2014, 50(23): 119-127.

［5］ 周麗宏, 陳自強, 黃國友, 等. 細胞打印技術(shù)及應用. 中國生物工程雜志, 2010, 30(12): 95-104.

［6］ 羊玢, 萬茂松, 王健康, 等. 三維UniGraphics實體造型軟件在理工高校的教學實踐與探討. 機械制造與自動化, 2012, 41(6): 151-153.

［7］ 劉常青, 彭翔, 張宗華, 等. 基于逆反工程(RE)多數(shù)據(jù)點云的處理. 工程圖學學報, 2002, 23(4): 33-39.

［8］ 鄧薇. MATLAB函數(shù)速查手冊. 修訂版. 北京: 人民郵電出版社, 2010.

［9］ 李富德, 王月蓮. 靜脈角和附近靜脈腔內(nèi)結(jié)構(gòu)的觀測及臨床意義. 長治醫(yī)學院學報, 1998, 12(4): 246-247.

［10］ Gibson Thomas J. Subclavian vein catheterization. Operative techniques in general surgery, 2001, 3(4): 203-216.

［11］ 陳灝珠. 實用內(nèi)科學. 第12 版. 北京: 人民衛(wèi)生出版社, 2005: 809.

［12］ 王宏志, 任宏. 多種體位鎖骨下靜脈穿刺中心靜脈置管術(shù)的應用. 中國急救醫(yī)學, 2004, 24(2): 146-147

［13］ 羅輝遇, 黎勝駒, 謝長江, 等. 鎖骨下靜脈置管異位相關因素分析. 臨床誤診誤治, 2009, 22(3): 7-8.

［14］ 邢順鵬, 聞大翔, 皋源, 等. 磁鐵引導法在鎖骨下靜脈置管中的臨床應用. 臨床麻醉學雜志, 2013, 29(6): 573-575.

［15］ Peter Weinstock, Sanjay P. Prabhu. Optimizing cerebrovascular surgical and endovascular procedures in hildren via personalized 3D printing. J Neurosurg Pediatr, 16(5): 584-589.

3D printing is one of hotspots in scientif i c research and commercial area due to that any products with complex shapes can be rapidly and accurately manufactured. In medical domain，3D printing is used to fabricate medical model，bone，organ and cell. Nowadays three dimension pictures can be achieved with using Computed Tomography（CT）via professional software. However，it is diff i cult to simulate operation. In this paper we developed a method to prepare transparent vein models using 3D printing. The method was based on patients’ CT slice graphs，and a 3D point cloud model and 3D virtual model were fi rstly designed by software such as MATLAB and CAD，and then two kinds of vein models in different precision degrees were successfully fabricated by 3D printer. The vein models can clearly display the detailed structure of subclavian vein，internal jugular vein and brachiocephalic vein，and would help to simulate operation and promote the method of subclavian vein catheterization. Besides，the method can be easily employed to make other vein model.

3D printing 3D design Vein model

上海交通大學醫(yī)工（理）交叉基金資助項目（YG2013MS39），上海交通大學本科生研究計劃（PRP）資助項目（T072PRP28015）

200240 上海交通大學（陳更新 劉崇偉 王錦輝）200127上海交通大學醫(yī)學院附屬仁濟醫(yī)院（殷焱 皋源 聞大翔）

*通信作者