劉子凡，許苑晶，柳毅浩，劉非，張敏，王文豪，鄧遷，顧越興，戴尅戎，王金武，△

(1.西南交通大學(xué)醫(yī)學(xué)院，成都 610031；2.上海交通大學(xué)生物醫(yī)學(xué)工程學(xué)院，上海 200030；3.上海市骨科內(nèi)植物重點(diǎn)實(shí)驗室，上海交通大學(xué)醫(yī)學(xué)院附屬第九人民醫(yī)院骨科，上海 200011；4.上海昕健醫(yī)療技術(shù)有限公司，上海 201800)

1 引 言

醫(yī)學(xué)影像可為臨床醫(yī)生診斷病情提供診療依據(jù)，傳統(tǒng)的醫(yī)學(xué)影像技術(shù)包括X線、CT、MRI等,輸出的為二維圖像，不便于立體觀察[1]。醫(yī)學(xué)圖像三維重構(gòu)技術(shù)通過組合醫(yī)學(xué)影像設(shè)備輸出的切片數(shù)據(jù)，重構(gòu)出三維圖像模型，被廣泛應(yīng)用于醫(yī)學(xué)教學(xué)、診斷、手術(shù)、仿真模擬等領(lǐng)域[2-3]，為醫(yī)生提供更逼真的顯示手段和定量分析工具[4]。三維重構(gòu)技術(shù)主要分為體繪制和面繪制[5]，體繪制以體素為單元，重構(gòu)時遍歷數(shù)據(jù)集、計算量大、耗時長，面繪制通過幾何單元拼接，重構(gòu)速度快、計算時間短[6]。

重構(gòu)精度是評價醫(yī)學(xué)圖像處理軟件的核心指標(biāo)，受閾值選擇影響產(chǎn)生的“啞鈴效應(yīng)(Dumb-bell Effect)”是造成重構(gòu)誤差的主要原因[7]。體模作為CT和MRI設(shè)備成像質(zhì)量的標(biāo)準(zhǔn)檢測工具，廣泛用于醫(yī)學(xué)影像設(shè)備的質(zhì)量研究[8-9]，其中Catphan 500型體模被國內(nèi)的計量、質(zhì)控部門廣泛使用[10]，其材料為有機(jī)玻璃，具有機(jī)械強(qiáng)度較高、耐高溫、耐腐蝕、易加工等優(yōu)點(diǎn)。因具備較高的準(zhǔn)確性和穩(wěn)定性而成為評價圖像設(shè)備成像質(zhì)量的標(biāo)準(zhǔn)驗證設(shè)備[11-12]，同時也應(yīng)用于CT成像算法、參數(shù)選擇的研究[13]。本研究通過比較體模實(shí)體的尺寸和使用Arigin3D Pro系統(tǒng)重構(gòu)的模型的測量尺寸，定量分析Arigin3D Pro(昕健醫(yī)療，中國)的重構(gòu)精度和穩(wěn)定性。

2 材料與方法

2.1 體模的準(zhǔn)備

2.1.1體模加工與制作 重構(gòu)軟件適用于CT和MRI兩種醫(yī)療設(shè)備生成的醫(yī)學(xué)影像，分別制作兩種成像設(shè)施的測試體模。由于CT和MRI的成像特性不同，CT體模分別設(shè)計了方形、圓柱形和塔形三種形狀，見圖1；MRI體模形狀設(shè)計為球形。體模的材料為有機(jī)玻璃，采用銑削和車削的工藝進(jìn)行加工，加工精度為±0.05 mm。

2.1.2掃描體模 CT和MRI如掃描參數(shù)不同，圖像質(zhì)量差別很大，為更好地對重構(gòu)結(jié)果進(jìn)行分析，掃描數(shù)據(jù)設(shè)置了統(tǒng)一的協(xié)議技術(shù)參數(shù)。

(1) CT數(shù)據(jù)掃描

在CT掃描過程中，為了驗證軟件處理數(shù)據(jù)的通用性，選取了Siemens SOMATOM Definition Flash(西門子，德國)、PhilipsBrilliance 64排螺旋CT(飛利浦，荷蘭)和聯(lián)影uCT 530(聯(lián)影，中國)三種品牌的CT，每臺機(jī)型的具體掃描參數(shù)，見表1。分別對方形、圓柱形和塔形3種體模進(jìn)行掃描，共得到9組掃描數(shù)據(jù)，以DICOM格式保存。

圖1 不同形狀的體模實(shí)物圖

表1 3種品牌的CT掃描協(xié)議技術(shù)參數(shù)Table 1 Technical parameters of three brands CT scanning protocols

(2)MRI數(shù)據(jù)掃描

選取聯(lián)影uMR 780核磁設(shè)備(聯(lián)影，中國)對球形體模進(jìn)行掃描，其掃描參數(shù)見表2，得到1組數(shù)據(jù)，以DICOM格式保存。

表2 聯(lián)影MRI掃描協(xié)議技術(shù)參數(shù)

2.1.3重構(gòu)三維模型 將掃描的DICOM數(shù)據(jù)導(dǎo)入重構(gòu)軟件，并對數(shù)據(jù)進(jìn)行閾值分割、區(qū)域增長等預(yù)處理，進(jìn)行分割重構(gòu)，見圖2。重構(gòu)結(jié)束后，以.stl格式導(dǎo)出保存。

圖2 Arigin 3D Pro軟件重構(gòu)不同形狀體模數(shù)據(jù)界面圖

2.1.4模型測量 使用游標(biāo)卡尺測量不同形狀體模實(shí)物，使用Geomagic Wrap(3D system, 美國)自帶的尺寸測量功能測量重構(gòu)模型的固定參數(shù)，每個參數(shù)測量3次，記錄并計算平均值。

2.2 模型測量

本研究使用Arigin 3D Pro對不同體模進(jìn)行重構(gòu)。通過抽簽隨機(jī)選取重構(gòu)經(jīng)驗為初級(重構(gòu)經(jīng)驗≤1年)、一般(重構(gòu)經(jīng)驗為2～3年)和豐富(重構(gòu)經(jīng)驗≥4年)的操作者各1名，對9組CT和1組MRI掃描數(shù)據(jù)分別進(jìn)行重構(gòu)，每組數(shù)據(jù)均重構(gòu)3次，且每次重構(gòu)均使用Geomagic Wrap測量參數(shù)3次。操作者1、2、3的重構(gòu)經(jīng)驗分別為無、3年和5年。實(shí)體測量與重構(gòu)模型測量均由同一個人操作。

2.3 測量指標(biāo)

(1)方形體模：長(mm)、寬(mm)、高(mm)；

(2)圓柱形體模：直徑(mm)、高度(mm)；

(3)塔形體模：塔基直徑(mm)、塔基高度(mm)、第一層直徑(mm)、第一層高度(mm)、第二層直徑(mm)、第二層高度(mm)；

(4)球形體模：內(nèi)直徑(mm)。

2.4 統(tǒng)計學(xué)處理

應(yīng)用IBM SPSS 20.0統(tǒng)計學(xué)軟件和GraphPad Prism 8 數(shù)據(jù)處理軟件進(jìn)行統(tǒng)計分析并繪圖。使用Kolmogorov-Smirnov檢驗判斷計量資料是否為正態(tài)分布。其中多組之間的比較判斷是否同時滿足方差齊性，若滿足，則使用單因素方差分析進(jìn)行比較；若不符合正態(tài)分布或不滿足方差齊性，則采用Kruskal-Walls檢驗進(jìn)行多組之間的比較。P<0.05認(rèn)為差異有統(tǒng)計學(xué)意義。

3 結(jié)果

3.1 軟件重構(gòu)模型的精度

本研究通過分析實(shí)體模型與軟件重構(gòu)模型各參數(shù)之間的誤差，評價軟件重構(gòu)的精度。

表3 實(shí)體模型與重構(gòu)模型各測量參數(shù)的誤差分析

由表3可知，軟件重構(gòu)模型與實(shí)體模型均有一定的誤差， CT重構(gòu)相對誤差的絕對值范圍占實(shí)體模型的0.09%～1.28%，MRI重構(gòu)相對誤差為0.09%。

3.2 軟件處理CT掃描數(shù)據(jù)的通用性

本研究使用了飛利浦、聯(lián)影、西門子3種品牌的CT設(shè)備進(jìn)行數(shù)據(jù)掃描，分析軟件重構(gòu)不同品牌掃描的數(shù)據(jù)是否具有差異性。為避免由人為因素引起的誤差，本研究以塔形重構(gòu)模型為例，選擇任意1名操作者的數(shù)據(jù)進(jìn)行分析，結(jié)果見圖3。

圖3 軟件重構(gòu)不同CT品牌掃描體模數(shù)據(jù)的差異性

由圖3可知，不同品牌CT重構(gòu)的模型參數(shù)范圍中，最小為第二層高度參數(shù)誤差0.11 mm，最大為塔基直徑參數(shù)誤差1.21 mm。對不同CT品牌重構(gòu)的塔基直徑參數(shù)進(jìn)行組內(nèi)相關(guān)系數(shù)計算，ICC=0.930>0.75，有較高的一致性，因此，該軟件處理CT掃描數(shù)據(jù)時不受設(shè)備品牌的限制，具有通用性。

3.3 操作者多次重構(gòu)三維模型的組內(nèi)差異

本研究中3名操作者均對同一形狀的體模分別重構(gòu)3次，且每次重構(gòu)均測量3次，通過分析同一體模各重構(gòu)模型間參數(shù)的差異，評價軟件重構(gòu)模型的穩(wěn)定性。為避免因不同CT品牌掃描數(shù)據(jù)引起的誤差，本研究以方形重構(gòu)模型為例，并選擇任意一名操作者的數(shù)據(jù)進(jìn)行分析，結(jié)果見圖4。

由圖4可知，方形體模的長、寬、高三次重構(gòu)之間均有一定的誤差，最大值為方形模型的長度參數(shù)誤差0.82 mm，最小值為寬度參數(shù)誤差0.56 mm，組內(nèi)差異占實(shí)體模型的0.56%～0.82%。

3.4 不同操作者重構(gòu)模型的組間差異

本研究選擇3名不同專業(yè)背景且對軟件的使用時間和熟練程度各異的操作者，通過分析不同操作者重構(gòu)結(jié)果的差異來驗證該軟件的普適性。

首先分析操作者1、操作者2、操作者3各自重構(gòu)出的模型與實(shí)體模型之間的相對誤差，見表4。

圖4 方形體模各重構(gòu)次數(shù)間參數(shù)的比較

表4 3名操作者的偏差分析結(jié)果

由表4可知，3種形狀的重構(gòu)模型與實(shí)體模型之間存在一定的誤差，操作者1的偏差范圍為0.34%～1.05%，操作者2的偏差范圍為0.07%～1.15%，操作者3的偏差范圍為0.18%～1.37%。進(jìn)一步分析3名操作者之間在重構(gòu)模型時是否具有差異，見表5。

表5 操作者1、操作者2、操作者3之間的組間比較

由表5可知，無論重構(gòu)哪種形狀的體模，3組之間均有參數(shù)具有統(tǒng)計學(xué)差異，為分析具體哪兩組之間不同，需進(jìn)一步做兩兩比較分析。本研究每個形狀的參數(shù)若符合正態(tài)分布，則采用單因素方差分析中的LSD法進(jìn)行多重比較分析；若參數(shù)若不符合正態(tài)分布，則采用Kruskal-Wallis檢驗中的兩兩比較法(Pairwise Comparisons)進(jìn)行分析。計算發(fā)現(xiàn)，方形體模的高、圓柱形體模的高和塔形體模的塔基直徑不滿足正態(tài)分布，限于篇幅，選擇這三個參數(shù)進(jìn)行分析，結(jié)果見圖5。兩組之間的具體比較值，見表6。

由圖5和表6可知，在方形、圓柱形和塔形體模的數(shù)據(jù)中，對操作者1與操作者3進(jìn)行比較時，差異均有統(tǒng)計學(xué)意義，而對操作者2與操作者3的數(shù)據(jù)進(jìn)行兩兩比較時，差異均無統(tǒng)計學(xué)意義；當(dāng)比較操作者1與操作者2時，差異在方形和圓柱形體模中有統(tǒng)計學(xué)意義。

表6 操作者1、操作者2和操作者3之間的兩兩比較

圖5 操作者1、操作者2和操作者3在不同形狀體模之間的兩兩比較

4 討論

本研究使用醫(yī)學(xué)重構(gòu)軟件Arigin 3D Pro進(jìn)行研究，根據(jù)體模的材質(zhì)使用固定閾值-155～1 000，避免了閾值變動引起的誤差，有利于結(jié)果比較。為了減少測量點(diǎn)的選取誤差，參考CT設(shè)備質(zhì)量體模標(biāo)準(zhǔn)檢測方法，采用Catphan體模參數(shù)測量，并將其作為標(biāo)準(zhǔn)實(shí)體模型用于對比。

通過Arigin 3D Pro重構(gòu)出四種形狀的三維模型，與相對應(yīng)的實(shí)體模型進(jìn)行測量參數(shù)的對比，排除客觀的機(jī)器掃描誤差后，重構(gòu)模型與實(shí)體模型存在一定的偏差，偏差不到實(shí)體模型的2%，相對模型的整體尺寸而言，該誤差可忽略，表明Arigin 3D Pro可重構(gòu)出比較精準(zhǔn)的三維模型。以往研究也表明，Arigin 3D Pro與醫(yī)療行業(yè)同類軟件相比，重構(gòu)精度相當(dāng)[14]。此外，本研究重構(gòu)了不同品牌的CT掃描數(shù)據(jù)，ICC顯示有較高的一致性，因此，該軟件適用于不同品牌掃描的數(shù)據(jù)，具有通用性。

本研究在驗證過程中任意選取一名操作者，研究發(fā)現(xiàn)不同重構(gòu)次數(shù)的結(jié)果存在差異，但對于整體模型而言，組內(nèi)差異極小，實(shí)際操作中可以忽略。此外，本研究發(fā)現(xiàn)不同操作者之間重構(gòu)的模型具有統(tǒng)計學(xué)差異，導(dǎo)致差異的因素可能與操作者的重構(gòu)資歷有關(guān)，操作者的重構(gòu)經(jīng)驗越豐富，重構(gòu)出的模型越一致。

綜上，本研究從定量評價角度出發(fā)對Arigin 3D Pro醫(yī)學(xué)重構(gòu)軟件的精度和穩(wěn)定性進(jìn)行評價，反映了Arigin 3D Pro軟件可重構(gòu)出精準(zhǔn)模型的特點(diǎn)，為后續(xù)在臨床中的應(yīng)用奠定了基礎(chǔ)。