背景技術：

本發(fā)明涉及基于醫(yī)學圖像數(shù)據(jù)的冠狀動脈中的血液動力學計算，并且更特別地涉及基于醫(yī)學圖像數(shù)據(jù)模擬冠狀動脈中的血流以用于冠狀動脈狹窄的非侵入式功能評價。

心血管疾?。–VD）在世界范圍是死亡的主要原因。在各種CVD之中，冠狀動脈疾?。–AD）占據(jù)那些死亡中的接近百分之五十。盡管有醫(yī)學成像和其它診斷模態(tài)中的顯著改進，針對CAD患者的過早發(fā)病率和死亡率中的增加仍舊非常高。用于冠狀狹窄的診斷和管理的當前臨床實踐涉及視覺上或通過定量冠狀動脈造影術（QCA）對患病血管的評價。這樣的評價為臨床醫(yī)生提供狹窄區(qū)段和父血管（parent vessel）的解剖學概覽，包括面積減小、病變（lesion）長度和最小內(nèi)腔直徑，但是不提供病變對通過血管的血流的影響的功能評價。通過將壓力導絲插入到狹窄的血管中來測量血流儲備分數(shù)（FFR）已經(jīng)示出是用于引導血管再生決定的更好選項，因為FFR在識別導致病變的局部缺血方面相比于侵入式造影術而言更加有效。QCA僅評估狹窄的形態(tài)學明顯度并且具有數(shù)個其它限制。基于壓力導絲的FFR測量涉及與將壓力導絲插入到血管中所必要的介入相關聯(lián)的風險，并且對于非常窄的狹窄，壓力導絲可能引起附加的血壓下降。

技術實現(xiàn)要素：

本發(fā)明提供用于從醫(yī)學圖像數(shù)據(jù)計算冠狀動脈中的血流的方法和系統(tǒng)。本發(fā)明的實施例利用醫(yī)學圖像數(shù)據(jù)中的冠狀動脈的解剖學測量來以患者特定的方式估計靜息血流速率。本發(fā)明的實施例使用患者的所估計的靜息血流速率來執(zhí)行血流模擬并且計算冠狀動脈中的導出的血液動力學度量。

在本發(fā)明的一個實施例中，從患者的醫(yī)學圖像數(shù)據(jù)提取冠狀動脈樹的患者特定解剖學測量。從所述冠狀動脈樹的患者特定解剖學測量估計針對所述冠狀動脈樹中的多個分支中的每一個的參考半徑?；卺槍λ龉跔顒用}樹的多個分支中的每一個的參考半徑計算流動速率。計算針對所述冠狀動脈樹的多個總流動速率估計，其中所述多個總流動速率估計中的每一個從來自多個分支世代中的相應一個的分支的流動速率計算?；谒龆鄠€總流動速率估計計算所述冠狀動脈樹的總流動速率。

通過參照以下詳細描述和附圖，本發(fā)明的這些和其它優(yōu)點將對本領域普通技術人員是顯而易見的。

附圖說明

圖1圖示了根據(jù)本發(fā)明的實施例的從醫(yī)學圖像數(shù)據(jù)計算冠狀動脈中的血流的方法；

圖2圖示了用于生成冠狀血管樹的患者特定解剖學模型的示例性結果；

圖3圖示了具有分配到每一個分支的世代號碼的冠狀動脈樹的示例；以及

圖4是能夠實現(xiàn)本發(fā)明的計算機的高級框圖。

具體實施方式

本發(fā)明涉及用于從醫(yī)學圖像數(shù)據(jù)計算冠狀動脈中的血流的方法和系統(tǒng)。在本文中描述本發(fā)明的實施例以給出用于計算冠狀動脈中的血流的方法的視覺理解。數(shù)字圖像通常包括一個或多個對象（或形狀）的數(shù)字表示。對象的數(shù)字表示在本文中通常在識別和操縱對象方面來描述。這樣的操縱是在計算機系統(tǒng)的存儲器或其它電路/硬件中實現(xiàn)的虛擬操縱。相應地，要理解的是，本發(fā)明的實施例可以在計算機系統(tǒng)內(nèi)使用存儲在計算機系統(tǒng)內(nèi)的數(shù)據(jù)來執(zhí)行。

基于計算流體動力學（CFD）的血流模擬可以用于估計血液動力學度量（諸如血流儲備分數(shù)（FFR））以用于冠狀動脈中的狹窄嚴重性的非侵入式功能評價。這樣的血流模擬典型地要求關于由每一個冠狀動脈供應的心肌的體積質量的一些信息以提供每一個冠狀動脈中的血流的初始估計。本發(fā)明的實施例不依賴于心肌的這樣的特性。本發(fā)明的實施例僅利用醫(yī)學圖像數(shù)據(jù)中的冠狀動脈的解剖學測量以估計冠狀動脈中的患者特定流動靜息血流速率。此外，本發(fā)明的實施例可以在左冠狀樹和右冠狀樹中獨立地應用，從而消除如在其它方法中的準備全面的冠狀動脈樹模型的需要。

圖1圖示了根據(jù)本發(fā)明的實施例的用于從醫(yī)學圖像數(shù)據(jù)計算冠狀動脈中的血流的方法。圖1的方法對表示患者的冠狀區(qū)的醫(yī)學圖像數(shù)據(jù)進行變換以生成患者的冠狀動脈的患者特定解剖學模型并且模擬患者的冠狀動脈中的血流。在步驟102處，接收患者的醫(yī)學圖像數(shù)據(jù)?？梢越邮諄碜砸粋€或多個成像模態(tài)的醫(yī)學圖像數(shù)據(jù)。例如，醫(yī)學圖像數(shù)據(jù)可以包括計算機斷層掃描術（CT）、Dyna CT、磁共振（MR）、造影術、超聲、單光子發(fā)射計算機斷層掃描術（SPECT）和任何其它類型的醫(yī)學成像模態(tài)。醫(yī)學圖像數(shù)據(jù)可以是2D、3D或4D（3D+時間）醫(yī)學圖像數(shù)據(jù)。醫(yī)學圖像數(shù)據(jù)可以直接從諸如CT掃描儀、MR掃描儀、造影術掃描儀、超聲設備等之類的一個或多個圖像獲取設備接收，或者醫(yī)學圖像數(shù)據(jù)可以通過加載針對患者的之前存儲的醫(yī)學圖像數(shù)據(jù)來接收。

在有利實施例中，在CT掃描儀上獲取3D冠狀CT造影術（CTA）圖像。CTA圖像確保通過使用注射到患者中的造影劑來對包括包含狹窄的（多個）血管的冠狀脈管系統(tǒng)適當?shù)爻上瘛Ｔ谠撾A段，臨床醫(yī)生可以被提供有通過在圖像上交互式地查看感興趣的病變（狹窄）而識別它們的選項。該步驟還可以在從圖像數(shù)據(jù)提取的患者特定解剖學模型上執(zhí)行（步驟104）?？商鎿Q地，可以使用用于冠狀動脈狹窄的自動檢測的算法在圖像數(shù)據(jù)中自動檢測狹窄，諸如在美國公開專利申請?zhí)?011/0224542中描述的用于冠狀動脈狹窄的動態(tài)檢測的方法，其通過引用并入本文。除了醫(yī)學圖像數(shù)據(jù)之外，還可以獲取其它非侵入式臨床測量，諸如患者的心率以及心臟收縮和心臟舒張血壓。

在步驟104處，從醫(yī)學圖像數(shù)據(jù)提取冠狀動脈樹的患者特定解剖學測量。在示例性實施例中，醫(yī)學圖像數(shù)據(jù)在靜息狀態(tài)獲取，并且從在靜息狀態(tài)獲取的圖像數(shù)據(jù)提取冠狀動脈的測量。冠狀動脈樹的患者特定解剖學測量可以是針對患者的完全冠狀動脈樹的患者特定解剖學測量，或者針對小于患者的完全冠狀動脈樹的任何部分的患者特定解剖學測量。在可能的實現(xiàn)方式中，冠狀動脈樹的患者特定解剖學測量可以是僅左冠狀動脈（LCA）樹或右冠狀動脈（RCA）樹的患者特定解剖學測量。在有利實施例中，冠狀動脈的測量通過生成從醫(yī)學圖像數(shù)據(jù)生成的冠狀動脈樹的患者特定解剖學模型來提取，但是本發(fā)明不限于此?；颊咛囟ń馄蕦W模型可以是患者的完全冠狀動脈樹的任何部分的患者特定解剖學模型。為了生成冠狀動脈樹的患者特定解剖學模型，可以使用自動化冠狀動脈中心線提取算法在3D醫(yī)學圖像數(shù)據(jù)中對冠狀動脈進行分段。例如，可以使用美國公開專利申請?zhí)?010/0067760描述的方法在CT體積中對冠狀動脈進行分段，其通過引用并入本文。一旦冠狀動脈中心線被提取，可以在中心線樹的每一個點處生成截面輪廓。每一個中心線點處的截面輪廓給出在冠狀動脈中的該點處的對應截面區(qū)域測量。然后生成針對經(jīng)分段的冠狀動脈的幾何表面模型。例如，在美國專利號7,860,290和美國專利號7,953,266中描述用于冠狀動脈的解剖學建模的方法，其二者通過引用并入本文。除了冠狀之外，患者特定解剖學模型可以包括主動脈根連同主動脈的近端部分。還使用類似算法來提取每一個狹窄的詳細3D模型，其包括近端血管直徑和面積、末端血管直徑和面積、最小內(nèi)腔直徑和面積以及狹窄的長度的量化。

圖2圖示了用于生成冠狀血管樹的患者特定解剖學模型的示例性結果。圖2的圖像200示出冠狀CTA數(shù)據(jù)。圖像210示出從CTA數(shù)據(jù)提取的中心線樹212。圖像220示出在中心線樹212的每一個點處提取的截面輪廓222。圖像230示出冠狀動脈、主動脈根和主動脈的近端部分的2D表面網(wǎng)格232。要理解的是，患者的冠狀樹的解剖學模型可以被輸出和顯示，例如在計算機系統(tǒng)的顯示屏幕上。

返回到圖1，在步驟106處，基于患者特定解剖學測量針對冠狀動脈樹的每一個分支估計參考半徑和靜息流動速率。特別地，針對每一個分支基于針對該分支的患者特定解剖學測量估計參考半徑，并且基于針對該分支的參考半徑計算針對每一個分支的靜息流動速率。已經(jīng)在動脈循環(huán)中的流動速率與血管參考半徑之間建立了以下冪定律關系，并且特別地對于冠狀循環(huán)：

其中q 是血管中的平均流動速率，r_ref是血管的參考半徑，k 是比例常數(shù)，并且n 是冪系數(shù)，其取2（對于大動脈）和3（對于小動脈）之間的值。相應地，特定冠狀動脈分支中的血流速率可以表述為針對該分支的參考半徑和冪系數(shù)的函數(shù)：

。

冪系數(shù)k 的值可以根據(jù)文獻數(shù)據(jù)設置，以便獲得典型靜息狀態(tài)流動速率值?？商鎿Q地，可以通過匹配侵入式測量與大患者數(shù)據(jù)庫中的計算值來設置k的值。在該情況下，針對患者計算的值可以直接是指流動速率，或者其它量，比如速度或血液動力學指數(shù)（FFR、CFR、iFR等）?？赡艿氖?，相同冪系數(shù)值可以用于所有患者。然而，還可能的是，可以基于諸如年齡、性別等之類的患者特性而針對不同患者來設置不同冪系數(shù)值。還可以針對不同血管（LAD、LCX、RCA、側分支等）來設置不同冪系數(shù)值。在可能的實施例中，k的值可以被調(diào)節(jié)以計及增加的冠狀流的條件，例如當患者處于充血狀態(tài)時。在該情況下，k的值從用于估計靜息流動速率的值增加成使得所估計的流動速率對應于患者的充血狀態(tài)處的血流。

對于患者特定幾何形狀，分支的半徑沿其中心線連續(xù)改變。另外，分支可以包含輕微到嚴重的狹窄，其進而可以是病灶性的（focal）、彌散性的（diffuse）等。因此，根據(jù)有利實施例，為了計算分支的參考半徑，向縱向變化的半徑應用運算符：

其中f₂是運算符，r 是分支的半徑，并且x 是沿分支的中心線的位置。在可能的實現(xiàn)方式中，運算符f₂可以計算沿分支的整個長度或分支的一部分的健康半徑的平均值。當在本文中使用時，“健康半徑”是指分支的健康（非狹窄）部分的半徑。如以上描述的，狹窄區(qū)可以被自動檢測或手動識別并且標記在冠狀動脈樹的患者特定解剖學模型上。相應地，可以從冠狀動脈樹的患者特定解剖學模型中自動識別每一個分支的健康部分。在另一可能的實現(xiàn)方式中，運算符f₂可以計算沿分支的整個長度或分支的一部分、排除健康半徑值的最大x%和最小y%的健康半徑的平均值。在另一可能的實現(xiàn)方式中，運算符f₂可以計算沿分支的整個長度或分支的一部分的健康半徑的最大或最小值。要理解的是，運算符f₂不一定限于這些運算并且其它可能的計算也可以用于估計分支的參考半徑。

在患者特定幾何形狀中，如果分支非常短（例如分支的長度小于閾值）或分支沿其整個長度表現(xiàn)出彌散性疾病（例如患病的分支的百分數(shù)大于閾值），針對該分支的參考半徑可以作為父代、同代和/或子代分支的參考半徑值的函數(shù)而計算：

其中是針對在冠狀動脈樹中分支從其裂開的父代分支的參考半徑值，是針對同代分支（即從相同父代分支分裂的另一分支）的參考半徑值，并且和是針對冠狀動脈樹中從該分支分裂的子代分支的參考半徑值。例如，分支的參考半徑可以基于父代、同代和子代分支的參考半徑值來確定為：

或。

一旦針對冠狀動脈樹中的每一個分支估計出參考半徑值，就使用以上描述的冪定律關系基于針對冠狀動脈樹中的每一個分支的參考半徑值來計算靜息流動速率?；趨⒖及霃结槍γ恳粋€分支計算的靜息流動速率是針對該分支中的血流的平均流動速率或整體（bulk）流動速率的估計。

在步驟108處，從多個分支世代計算總靜息流動速率估計。圖3圖示了具有分配到每一個分支的世代號碼的冠狀動脈樹300的示例。如圖3中所示，冠狀動脈樹300的根分支302具有世代號碼0，并且在每一個分叉處世代號碼增加一?？梢愿鶕?jù)每一個世代號碼的分支估計針對冠狀動脈樹的總流動速率的單獨的估計。使用具有世代號碼g 的分支計算的針對冠狀樹的總流動速率估計計算如下。在估計總流動之前，向每一個分支分配權重w_i，其表示用于針對該分支的估計的參考半徑的正確性的置信度值?？梢曰诜种У拈L度和/或患病分支的百分數(shù)來對分支進行加權。相應地，短分支（諸如在圖3中具有等于1的世代號碼的底部分支304）或整個患病的分支（諸如在圖3中具有等于2的世代號碼的彌散性患病的分支306）被分配低置信度值，而沒有半徑不規(guī)則性的長分支被分配大置信度值。權重（置信度值）可以被歸一化到0（最小置信度）和1（最大置信度）之間的范圍。權重還可以從在醫(yī)學圖像中觀察到的特性（諸如圖像質量、鈣化水平等）導出。例如，具有相對高鈣化的分支可以被分配低置信度，而具有低或零鈣化的分支可以被分配非常高的置信度值。如果在用于重構分支的圖像的部分中存在偽像（例如運動偽像、金屬偽像、射束硬化偽像等），權重還可以是該感興趣的區(qū)的局部圖像質量的函數(shù)。然后使用針對分支計算的流動速率和分配到分支的權重基于來自世代g 的分支來估計針對整個冠狀動脈樹的總流動速率：

其中索引i是指來自世代g的所有分支和具有小于g的世代號碼的所有終端分支。終端分支是不具有任何子代分支（即不具有到具有更高世代號碼的另外的分支的任何分叉）的分支。在示例性實現(xiàn)方式中，基于來自世代g的分支的總流動速率估計可以計算為：

。

通過基于來自g_min和g_max之間的每一個世代g 的分支計算相應總流動速率估計來計算多個總流動速率估計。在有利實現(xiàn)方式中，最小世代級g_min可以是0，但是還可以大于0，如果根節(jié)點非常短的話。最大世代級可以被設置以確定在計算冠狀動脈樹的總靜息流動速率中使用多少世代。在有利實現(xiàn)方式中，用于最大世代級g_max的值可以設置成3或4。更高世代的分支變得越來越小，其使得使用更高世代分支的參考半徑和對應流動速率的準確估計更加困難。另外，當解剖學模型是從醫(yī)學圖像重構的時，可能在模型中未計及小的側分支。因而，世代號碼越高，未考慮到的側分支的數(shù)目將越高，從而造成流動速率估計中的更大的誤差。

返回到圖1，在步驟110處，基于從多個分支世代計算的總流動速率估計而計算冠狀動脈樹的總靜息流動速率。為了改進針對冠狀動脈樹的總流動速率估計的準確性，從多個不同分支世代計算的流動速率估計被用于計算最終總流動速率值。特別地，從g_min和g_max之間的每一個世代g 計算的總流動速率估計被用于計算針對冠狀動脈樹的最終總流動速率值。如以上描述的，最小世代級g_min可以為0，但是還可以大于0，如果根節(jié)點非常短的話，并且最大世代級可以被設置以確定在計算冠狀動脈樹的總靜息流動速率中使用多少世代。在估計針對冠狀動脈樹的總流動速率之前，權重v_i被分配到每一個世代號碼，其表示針對從具有對應世代號碼的分支計算的總流動速率估計的正確性的置信度值。對于冠狀動脈樹的患者特定解剖學模型，低世代號碼可以被分配大權重，而大世代號碼可以被分配低權重，因為隨著世代號碼增加，可能錯過較小側分支，這可能造成流動速率值的低估。例如，權重可以具有與世代號碼的逆向關系。權重（置信度值）可以歸一化到0（最小置信度）和1（最大置信度）之間的范圍。最終總流動速率值被估計為針對各種世代的總流動速率估計和分配到世代的對應權重的函數(shù)：

其中索引j 是指g_min和g_max之間的世代。在示例性實現(xiàn)方式中，針對冠狀動脈樹的總流動速率可以被計算為：

。

在步驟112處，使用總靜息流動速率作為入口邊界條件來模擬冠狀動脈樹中的血流。在有利實現(xiàn)方式中，一個或多個計算流體動力學（CFD）模擬被用于模擬冠狀動脈樹中的血流。針對冠狀動脈樹計算的總靜息流動速率被用作用于CFD模擬的患者特定入口邊界條件。在CFD模擬中，血液被建模為牛頓流體，并且速度場通過在剛性壁假設之下對離散納維-斯托克斯方程（連續(xù)性和動量方程）進行數(shù)值求解來獲得。離散納維-斯托克斯方程用于隨時間推移遞增地模擬冠狀動脈內(nèi)的壓力和血流的速度。冠狀動脈的患者特定解剖模型還被輸入到CFD建模以便約束基于患者特定解剖模型的血流模擬。CFD模擬導致隨時間推移針對通過冠狀動脈的血流的速度和壓力的模擬值。

在步驟114處，基于冠狀動脈樹中的模擬血流計算一個或多個血液動力學度量。例如，可以基于CFD血流模擬來計算患者特定血液動力學參數(shù)。CFD模擬導致隨時間推移針對通過冠狀動脈的血流的速度和壓力的模擬值。這些模擬值可以用于計算各種血液動力學參數(shù)，其可以用于評價冠狀動脈疾病。例如，流動速率和壓力下降可以用于評價狹窄的嚴重性，并且壁面剪應力可以用于評價斑塊（plaque）形成。諸如血流儲備分數(shù)（FFR）、靜息Pd/Pa、瞬態(tài)波-無波比（iFR）等之類的血液動力學度量可以從冠狀動脈中的模擬流動速率和壓力針對特定狹窄來計算。例如，靜息Pd/Pa是指在靜息狀態(tài)處狹窄上的壓力下降（遠離狹窄的壓力（Pd）和主動脈壓力（Pa）的比）并且可以使用總靜息流動速率作為邊界條件直接從CFD模擬中的模擬壓力計算。iFR是指心臟舒張中的無波周期期間的Pd/Pa之比并且還可以在靜息狀態(tài)處使用總靜息流動速率作為入口邊界條件從CFD模擬中的模擬壓力計算。

FFR被定義為狹窄血管中的最大血流與正常血管中的最大血流之比，并且用于表征狹窄的嚴重性?？梢葬槍跔顒用}樹中的狹窄通過在模擬充血狀態(tài)下模擬冠狀動脈樹的患者特定解剖模型中的流動和壓力并且計算遠離狹窄的時間平均的壓力（Pd）關于充血狀態(tài)處主動脈中的平均壓力（Pa）的比來近似FFR。在一個可能的實施例中，為了模擬充血處的血流和壓力，可以使用不同的運算符f₁來執(zhí)行步驟106以從參考半徑計算每一個的流動速率，使得針對每一個分支估計充血流動速率而不是靜息流動速率。步驟108和110如以上描述地那樣執(zhí)行，但是在該情況下導致總充血流動速率被針對冠狀動脈樹計算，并且在步驟112中，使用總充血流動速率作為入口邊界條件來執(zhí)行血流模擬，從而導致冠狀動脈樹中的充血血流和壓力的模擬。在另一可能的實施例中，針對冠狀動脈樹計算的總靜息流動速率可以被修改以反映充血狀態(tài)而不是靜息狀態(tài)。例如，在示例性實現(xiàn)方式中，冠狀動脈樹的總靜息流動速率可以如以上描述地那樣在步驟106-110中計算并且然后可以基于總靜息流動速率和基于針對患者非侵入式測量的心率、收縮血壓和舒張血壓計算的均值動脈壓力（MAP）來計算冠狀動脈樹的每一個分支的終端處的靜息微脈管阻力。每一個分支的終端處的充血微脈管阻力然后可以基于靜息微脈管阻力來計算，并且充血微脈管阻力可以用作用于冠狀動脈樹中的充血血流和壓力的模擬的邊界條件。關于計算靜息微脈管阻力MAP和充血微脈管阻力的附加細節(jié)描述在題為“Method and System for Non-Invasive Functional Assessment of Coronary Artery Stenosis”的美國專利公開號2013/0246034和題為“Method and System for Non-Invasive Functional Assessment of Coronary Artery Stenosis”的美國專利公開號2014/00058715中，其以其整體通過引用并入本文。在另一可能的實現(xiàn)方式中，針對冠狀動脈樹計算的總靜息流動速率可以使用描述在2015年1月19日提交的題為“System and Method for Mapping Patient Data from One Physiological State to Another Physiological State”的美國專利申請?zhí)?4/599,678中的方法而直接映射到總充血流動速率，其以其整體通過引用并入本文?？偝溲鲃铀俾嗜缓罂梢杂米饔糜诠跔顒用}樹中的充血血流和壓力的模擬的入口邊界條件。

要理解的是，其他血液動力學度量可以使用處于靜息、充血或其它生理狀態(tài)的血流模擬來類似地計算。在可替換的實施例中，血液動力學度量可以直接從針對冠狀動脈樹的分支估計的流動速率計算而不執(zhí)行血流估計。

根據(jù)本發(fā)明的有利實施例，總靜息流動速率可以在不對LCA和RCA二者進行建模的情況下針對左冠狀動脈（LCA）或右冠狀動脈（RCA）獨立地計算。相應地，圖1的方法的優(yōu)點在于該方法可以針對LCA或RCA中的一個單獨地執(zhí)行以便確定入口邊界條件并且執(zhí)行LCA或RCA中的血流模擬，并且計算針對LCA或RCA中的狹窄的一個或多個血液動力學度量。例如，圖1的方法可以在LCA或RCA中的一個中單獨執(zhí)行以評價LCA或RCA中的目標狹窄，而不執(zhí)行不包含目標狹窄的LCA或RCA中的另一個的任何建模。

如以上描述的，圖1的方法估計冠狀動脈樹中的總靜息流動速率和所計算的血流。該方法可以類似地應用于其它動脈血管樹，諸如也用于腎循環(huán)和腦循環(huán)。此外，圖1的方法可以類似地應用于估計針對處于除靜息狀態(tài)之外的其它生理狀態(tài)（例如充血）的動脈樹的總流動速率，例如通過改變用于將針對每一個分支的參考半徑映射到靜息流動速率的運算符f₁。

用于計算冠狀動脈樹中的血流和冠狀動脈狹窄的非侵入式評價的以上描述的方法可以實現(xiàn)在使用公知計算機處理器、存儲器單元、存儲設備、計算機軟件和其它組件的計算機上。這樣的計算機的高級框圖在圖4中圖示。計算機402包含處理器404，其通過執(zhí)行定義這樣的操作的計算機程序指令來控制計算機402的總體操作。計算機程序指令可以存儲在存儲設備412（例如磁盤）中并且當期望執(zhí)行計算機程序指令時被加載到存儲器410中。因此，圖1的方法的步驟可以由存儲在存儲器410和/或儲存器412中的計算機程序指令定義并且由執(zhí)行計算機程序指令的處理器404控制。諸如CT掃描設備、MR掃描設備、超聲設備等之類的圖像獲取設備420可以連接到計算機402以向計算機402輸入圖像數(shù)據(jù)?？赡艿氖菍D像獲取設備420和計算機402實現(xiàn)為一個設備。還可能的是圖像獲取設備420和計算機402通過網(wǎng)絡無線地通信。在可能的實施例中，計算機402可以關于圖像獲取設備420遠程定位，其可以作為服務器或基于云的服務的部分而執(zhí)行方法步驟。計算機402還包括一個或多個網(wǎng)絡接口406以用于經(jīng)由網(wǎng)絡與其它設備通信。計算機402還包括使得能夠實現(xiàn)與計算機402的用戶交互的其它輸入/輸出設備408（例如顯示器、鍵盤、鼠標、揚聲器、按鈕等）。這樣的輸入/輸出設備408可以結合計算機程序集合來用作注釋工具以注釋從圖像獲取設備420接收的體積。本領域技術人員將認識到，實際計算機的實現(xiàn)方式也可以包含其它組件，并且圖4是用于說明性目的的這樣的計算機的一些組件的高級表示。

前述具體實施方式要理解為在每一個方面都是說明性和示例性的，但是不是限制性的，并且本文所公開的發(fā)明的范圍不從具體實施方式確定，而是從如根據(jù)專利法所準許的完整范圍解釋的權利要求確定。要理解的是，本文所示出和描述的實施例僅僅說明本發(fā)明的原理并且各種修改可以由本領域技術人員在不脫離本發(fā)明的范圍和精神的情況下實現(xiàn)。本領域技術人員可以在不脫離本發(fā)明的范圍和精神的情況下實現(xiàn)各種其它特征組合。