基于多像素值提取優(yōu)化的重映射方法構建肝脾破裂傷三維體繪制模型

孫 堅，譚 珂，高 原，蔡守旺，王鵬飛

1 解放軍醫(yī)學院，北京 100853；2 解放軍總醫(yī)院研究生院 教育技術中心，北京 100853；3 解放軍總醫(yī)院第一醫(yī)學中心 肝膽胰外科學部，北京 100853

利用CT等醫(yī)療探查設備獲取數(shù)據(jù)信息并進行三維可視化建模，一直以來是醫(yī)學圖像領域的研究熱點。在肝膽外科手術研究領域，對肝腫瘤的形態(tài)和位置進行自動或手動標注分割，其結果可方便開展三維可視化建模分析，進而確定手術方案并實現(xiàn)術中規(guī)劃[1-3]。目前可視化建模的基本方法主要有移動立方體算法(marching cubes，MC)[4]、體繪制法[5]、多平面重建法[6-7]、最大密度投影法[8]、容積重建法[9]、有限元法[10]。Custodio等[11]提出了一種擴展的行進數(shù)據(jù)集算法，將MC中像素值的0、1標記改進為+、?和=，解決了拓撲結構劃分和網(wǎng)格質量提升的沖突。De Sales Tsuzuki等[12]提出了一種基于傳播的MC算法，采用開放邊界條件進行循環(huán)以確定體元的迭代處理順序，在高數(shù)據(jù)量的條件下保持了較高的運算速度。Gu和Kim[13]使用以CPU為核心的Bunyk算法優(yōu)化了光線投射算法，使用逐像素的交集表以保證非凸網(wǎng)格的準確性。上述研究在成像算法方面取得了建設性成果，但在實際應用中存在一定的欠缺，如描述器官的受傷情況。在日常生活中，車禍、墜落、鈍擊等意外沖擊造成脾損傷的發(fā)生率在腹部創(chuàng)傷中可達40%～50%，在腹部閉合性損傷中，脾破裂占20%～40%[14]。由于沖擊會造成肝脾包膜破裂、實質裂傷，甚至門靜脈、肝靜脈受損，肝脾葉血腫、離斷、異位、損毀，腹腔壓力變化造成血腫，整體解剖結構發(fā)生改變，相較于正常肝脾或肝腫瘤等結構，其CT圖像視野混亂，數(shù)據(jù)的分類、識別、標注易在實質損毀區(qū)域混淆，上述研究成果不能完全解決此類問題。為解決問題并提供思路，本文在體繪制流程中設計了多像素值的自動提取，并將體繪制中空間數(shù)據(jù)與增強CT的DICOM數(shù)據(jù)建立重映射，繪制的模型可清楚觀察肝脾破裂傷的解剖結構與病理信息，為后續(xù)相關研究提供了一定基礎。

資料和方法

1資料 以2013年1月-2018年11月在解放軍總醫(yī)院第一醫(yī)學中心肝膽胰外科學部15例肝脾破裂傷患者增強CT的DICOM數(shù)據(jù)作為基礎，其創(chuàng)傷來源為車禍沖擊、高空墜落、重物鈍擊等事故。依據(jù)人民衛(wèi)生出版社出版第8版《外科學》中關于肝外傷和脾外傷的分級以及肝膽胰外科學部診療實際情況，篩選出肝外傷Ⅱ～Ⅳ級共7例，脾外傷Ⅲ級4例(1例肝脾均破裂)，總計11例病例作為肝脾破裂傷可視化數(shù)據(jù)集，其中9例男性，1例女性，年齡21～47歲。

2設備 本實驗圖形工作站配置為：Intel Core i7-7700HQ@2.80 gHz四核；顯卡Nvidia GeForce GTX 1060；內(nèi)存16 gB；Windows 10 64位操作系統(tǒng)(DirectX 12)。在進行可視化模型構建與演示時系統(tǒng)流暢無卡頓。

3多像素值提取的CT數(shù)據(jù)體繪制 針對肝脾破裂傷的CT圖像解剖結構異常等特征，需對CT數(shù)據(jù)先進行體繪制。為更好體現(xiàn)繪制模型區(qū)分器官、血管與骨骼等不同部分，本文對傳統(tǒng)光線投射體繪制算法進行了多像素值提取改進，其實現(xiàn)基于視覺化工具函式庫(Visualization Tool Kit，VTK)，該函式庫是一種廣泛應用于計算機圖形、圖像處理、可視化的開源軟件系統(tǒng)。方法如下：1)利用DICOM Image Reader類庫讀取CT DICOM數(shù)據(jù)；2)利用Get Scalar Pointer (i,j,k)函數(shù)遍歷DICOM數(shù)據(jù)中所有體素，由于CT獲取的人體組織HU值為0~600[15]，所以創(chuàng)建容量為600的數(shù)列，歸類統(tǒng)計CT數(shù)據(jù)中大于0、小于600的HU值，并使用Bubble Sort排序方法對HU值的出現(xiàn)次數(shù)進行由大到小排序，出現(xiàn)頻率最高的3種HU值即對應骨骼、血管與器官的像素值信息；3)利用Opacity Transform(透明度變換)與Color Transform Function(像素值變換)方法確定多種透明度值和像素值的映射方式。傳統(tǒng)方法將像素值參數(shù)設為定值，本文優(yōu)化為將上述3個HU值作為參數(shù)輸入，以確保后續(xù)對易混淆區(qū)域有效繪制可視化模型；4)利用Volume Ray Cast Composite Function計算光照效應，將空間平行光線與三維影像數(shù)據(jù)求交，并將相交部分的光線線段等值劃分以確定采樣點；5)利用Set Interpolation Type To Linear方法，根據(jù)采樣點在空間坐標系下的位置確定與其距離最近的8個體素，將各體素值輸入三線性插值公式，計算出采樣點的像素值和透明度值；6)利用Volume Ray Cast Mapper模擬光線通過整個三維數(shù)據(jù)矩陣到達觀察點的過程，將相交部分光線線段中采樣點值由后至前累加；7)為體繪制過程完整渲染，需利用Render Window與Render Window Interacto定義繪制過程中顯示器鏡頭位置與光照模型。

由于增強CT體素中器官、血管和器官實質的HU值均高于出血時血液的HU值，故繪制的可視化模型不含有出血血液，可以清晰準確地顯示患者受傷肝、脾的解剖結構。具體實現(xiàn)流程見圖1。

圖1 多像素值提取的體繪制流程圖Fig.1 Flowchart of body drawing for multi-pixel value extraction

4體繪制模型重映射與多層插值構建體繪制模型 前面利用CT數(shù)據(jù)繪制出的體模型是關于整個腹腔的可視化模型，為提取肝脾破裂傷的臟器實質和脈管部分，進行體繪制模型重映射與多層插值，具體步驟如下：1)利用自動勾選和手動劃分相結合的方式，去除體繪制模型中大部分無用區(qū)域。2)利用矩陣變換，將剩余三維體繪制模型所包含的體素重映射回二維CT數(shù)據(jù)[16]。3)使用多層插值法刪除模型冗余細節(jié)。首先在CT數(shù)據(jù)中選擇若干張圖片并勾選不需要的部分HU值，然后對其進行距離變換以查找多余HU值的主干部分，最后以各層主干部分為基礎，在整個空間進行插值并分割刪除[17]。4)更新體繪制模型的剩余數(shù)據(jù)并重新繪制，得到可觀察的最終模型。其具體實現(xiàn)演示如圖2。

圖2 體繪制模型重建步驟A：繪制全DICOM數(shù)據(jù)模型；B：首先去掉骨骼，然后勾選去除肝以外區(qū)域；C：重映射與多層插值法去除冗余部分；D：更新模型，完成肝破裂可視化模型構建Fig.2 Reconstruction steps of volume renderingA: Draw the full DICOM data model; B: First remove the bones, then remove the area outside the liver; C: Remove the redundant parts by using remap and multi-layer interpolation; D: Update to complete the liver rupture visualization model

結 果

1方法應用與案例分析 本文設計方法于2020年8-9月在解放軍總醫(yī)院第一醫(yī)學中心肝膽胰外科學部針對肝脾破裂傷增強CT數(shù)據(jù)進行體繪制模型構建應用，模型構建時間為4.01～5.83s。除工程師操作外，對學部副主任醫(yī)師、主治醫(yī)師和醫(yī)師進行培訓，經(jīng)1.5～2 h的練習后，可以獨立構建肝脾破裂傷模型，也能通過模型進行診斷與確定治療方案。

本文所構建患者的肝脾破裂可視化模型，能夠在三維空間多角度、多層次觀察，清楚分辨器官實質的受損部位和受損情況，并實際以此為依據(jù)分析患者病情、輔助決定治療方法與手段。如圖3A所示，患者為25歲男性，因重物砸傷腹部，其體繪制的模型可見門脈右后支五、六、七段均破裂；顯示肝右靜脈、右后下靜脈均在右肝破裂范圍內(nèi)，因肝右靜脈、右后下靜脈破裂，單純切除肝右后葉會造成肝右前葉血液回流受限，所以合理的手術方案應為右半肝切除術。如圖3B所示，患者為27歲男性，因肝破裂術后再次出血而復查，可見右肝表面包裹性壞死，門脈與靜脈系統(tǒng)完好，考慮給予對癥治療，定期進行肝CT掃描與血紅蛋白監(jiān)測，并服用保肝、止血、抗炎類藥物。圖4所示，患者為38歲女性，因車禍沖擊致使脾破裂，模型可見脾門處有3處裂傷，脾膈面上極有長破裂口，因脾質地較脆、破裂口較多，遂考慮行脾切除術。

圖3 肝可視化模型圖A：門脈右后支五、六、七段均破裂，肝右靜脈、右后下靜脈均在右肝破裂范圍內(nèi)；B：右肝表面包裹性壞死，門脈與靜脈系統(tǒng)完好Fig.3 Visualization model of the liverA: The right posterior branch of the portal vein in segments V, VI and VII are ruptured, and the right hepatic vein and right posterior inferior vein are within the ruptured right liver; B: Encapsulated necrosis can be found on the surface of the right liver, and portal and venous system are intact

圖4 脾可視化模型圖Fig.4 Visualization model of the spleen

2 與傳統(tǒng)MC算法構建模型對比 模型對比如圖5所示。相較于MC算法構建的三維面模型，體繪制三維模型更容易實現(xiàn)從整體受傷區(qū)域觀察與其他器官組織的毗鄰關系。因為注射后的造影劑在影像數(shù)據(jù)中顯示為高信號值，在體繪制過程中可以更加清楚地觀察肝脾內(nèi)脈管系統(tǒng)的形態(tài)走向以及受損情況。不同于骨骼外傷等病情，肝脾等臟器受外傷會導致實質表面凹凸不平，易與血液凝塊混合，血腫時更會與其他臟器粘連，傳統(tǒng)MC重建出的可視化模型凹凸不平，難以將數(shù)據(jù)信息準確呈現(xiàn)給外科醫(yī)生，而過度的平滑處理則會遺失部分解剖信息。

圖5 相同一套CT數(shù)據(jù)繪制的面模型(A)、體模型(B)與本文體繪制模型(C)對比(針對肝破裂傷，傳統(tǒng)面繪制模型在繪制時難以描述受損區(qū)域，而傳統(tǒng)體繪制模型難以區(qū)分肝的脈管與實質)Fig.5 Comparison of the surface rendering model (A), volume rendering model (B) and our model (C) established by the same set of CT data(for liver rupture injury, the traditional surface rendering model is difficult to describe the damaged area, while the traditional volume rendering model is difficult to distinguish the vascular and parenchyma of the liver)

3 與傳統(tǒng)體繪制構建模型對比 傳統(tǒng)體繪制模型將DICOM數(shù)據(jù)的HU值與像素值建立線性映射，其最終繪制的肝脾可視化模型難以清晰區(qū)分實質與脈管部分[18]，易造成影像結果的誤判，而解決上述問題的傳統(tǒng)方法需提前對數(shù)據(jù)進行分類標記，對CT圖像讀取能力、分辨有用與無用信息的要求較高。本文自動提取HU值并對應像素值，進而建立顏色區(qū)間，最終的模型能準確描述受損區(qū)域的解剖結構。在進行傳統(tǒng)體繪制模型構建時，需手動一次性去除無用數(shù)據(jù)，操作要求高，模型與實際情況相差大，而本文方法將體模型包含的三維數(shù)據(jù)與CT中HU值建立映射，在全DICOM數(shù)據(jù)體繪制的前提下，交替進行三維空間分割與二維圖像分割，能夠方便直觀地構建肝脾破裂傷可視化模型。

討 論

利用CT等影像學設備對肝脾的位置、體積、動靜脈與膽道結構進行成像，一直以來是外科醫(yī)生了解患者情況并進行診斷和手術設計的必要手段。本文優(yōu)化的三維重建體繪制方法可以同時展示患者的器官、組織、血管等信息，便于直觀地將三維解剖結構準確地呈現(xiàn)給外科醫(yī)生。其成果首先可輔助術前準備時更細致地觀察患者受傷位置、受傷情況及各器官組織的空間結構變化等一系列解剖信息，輔助指導外科醫(yī)生進行臨床診斷、設計手術方案與并開展術前討論；其次可以與三維空間定位傳感技術相結合，在術中實時更新三維模型的空間注冊信息，從而為術中導航和輔助手術操作提供支持；且更便于對比評估術后患者恢復的情況。

上述方法對于提升手術質量與評估術后恢復效果有著十分重要的價值。其重建結果也可以應用于進一步構建三維空間網(wǎng)格模型，為個性化手術仿真模擬提供器官實質和血管的疊加物理仿真計算模型；將該模型在臨床教學中應用，可方便學生快速直觀地理解相關疾病的解剖結構，具有良好的教學應用前景。但當特殊情況導致器官內(nèi)血管變窄、信號變低時，其CT數(shù)據(jù)區(qū)分度差，如肝中央?yún)^(qū)損傷引起的大范圍血腫形成包塊，肝中靜脈受壓變細，造影劑信號值與肝實質信號值相似，難以構建出準確直觀的體繪制模型。

目前器官組織的三維重建在外科手術領域展現(xiàn)了良好的發(fā)展前景，其研究成果近幾年在三維數(shù)據(jù)重建、混合現(xiàn)實顯示、術中規(guī)劃導航、支持力反饋的虛擬手術模擬等方面有著長足的進步。但目前為止，其成果算法仍未完全平衡好快速性與準確性之間的關系，本文為相關問題發(fā)展與解決提供了改進思路。未來醫(yī)學三維重建會向著精確化、快速化、設備小型化、掃描自動化、配準實時化的方向發(fā)展，其研究成果會更加廣泛地運用于肝膽外科、顱腦外科、泌尿外科等一系列學科的臨床手術中，大幅提高外科手術的效率和精準性。