以下總體上涉及醫(yī)學成像，并且利用針對管腔結(jié)構(gòu)的分割的具體應用進行了描述。

背景技術(shù)：

文獻指出，冠狀動脈疾病是全球唯一的最大死亡原因。冠狀動脈疾病是經(jīng)由血液向心肌供氧的冠狀動脈的管腔變窄。冠狀動脈的管腔直徑不等。例如，健康成年個體的左主冠狀動脈管腔的直徑通常為4.5±0.5mm，遠端左前降支冠狀動脈管腔通常為1.9±0.4mm。

評估冠狀管腔變窄(例如，狹窄)的一種無創(chuàng)技術(shù)是心臟計算機斷層攝影血管造影(ccta)。在ccta中，使用冠狀動脈的一幅或多幅三維(3d)圖像(例如，計算機斷層攝影(ct)圖像)來識別和量化管腔變窄。另一定量測量使用基于管腔的3d模型的血流儲備分數(shù)(ffr)模擬的壓差。

可視成像分辨率或精確度(例如，利用ccta成像)通常為約1.5mm，例如，能夠準確分辨圖像中的差異的一個空間尺寸。使用圖像來確定小于可見分辨率的三倍到四倍的管腔的變窄對以下造成影響：首先影響是否檢測到變窄；其次，影響受到變窄影響的、對管腔的結(jié)構(gòu)進行分割以量化測量結(jié)果的準確度。重建算法通常利用低通濾波器來使圖像平滑，這會使得管腔邊界模糊。在管腔直徑小于低通濾波器支持能力的情況下，平滑處理可能導致分割算法過高估計管腔直徑。動脈中的鈣沉積使動脈變窄，這可能導致成像中的暈開(blooming)偽影和對比潛在的真實管腔更大的管腔進行模擬的典型的分割算法中的暈開偽影。

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

在本文中描述的各方面解決了上述問題和其他問題。

以下描述了使用圖形最小切割算法利用改善的圖形來檢測血管管腔的亞分辨率變窄并分割血管管腔的系統(tǒng)和方法。

在一個方面中，一種成像系統(tǒng)包括：亞分辨率管腔變窄檢測器，其通過中心線輪廓分析來檢測圖像體積中的血管管腔的亞分辨率變窄，并且通過利用所檢測的亞分辨率變窄修正近似的可見管腔直徑來計算亞分辨率確定的直徑。

在另一方面中，一種分割管腔血管的方法包括：利用對圖像體積中的血管管腔中的近似的可見管腔直徑的中心線輪廓分析來檢測亞分辨率變窄；通過基于所檢測的亞分辨率變窄修正近似的可見管腔直徑來計算亞分辨率確定的直徑。

在另一方面中，一種成像系統(tǒng)包括：圖形最小切割變化分割器，其利用在由節(jié)點表示的圖像體積的經(jīng)分割的區(qū)域中的每個體素來構(gòu)建圖形。每個節(jié)點連接到管腔節(jié)點和背景節(jié)點。將每個節(jié)點連接到所述管腔節(jié)點的邊利用可能性li與可能性ld的乘積進行加權(quán)，其中，所述可能性li是體素根據(jù)最接近的中心線體素的強度以及經(jīng)訓練的模型而是經(jīng)分割的血管管腔的部分的可能性，所述可能性ld是體素基于所述體素距最接近的中心線點的空間距離以及亞分辨率確定的直徑而是所述經(jīng)分割的血管管腔的部分的可能性。將每個節(jié)點連接到所述背景節(jié)點的邊通過1-li*ld進行加權(quán)。鄰近體素之間的邊通過平滑度懲罰進行加權(quán)。所述圖形最小切割變化分割器并對所述圖形進行切割，所述的對所述圖形進行切割基于最小成本將所述圖形分成連接到所述管腔節(jié)點的第一部分和連接到所述背景節(jié)點的第二部分，并且所述第一部分中的所述節(jié)點生成經(jīng)分割的血管管腔的體積。

附圖說明

本發(fā)明可以采取各種部件和部件的布置，以及各個步驟和步驟的安排的形式。附圖僅出于圖示優(yōu)選實施例的目的，并且不得被解釋為對本發(fā)明的限制。

圖1示意性地圖示了與成像系統(tǒng)相結(jié)合的示例性亞分辨率管腔變窄檢測和分割系統(tǒng)。

圖2圖示了使用半峰全寬(fwhm)規(guī)則的近似的管腔直徑的范例。

圖3示意性地圖示了在可見成像分辨率下對變化的血管直徑的可視化的范例。

圖4示出了減小的強度管腔變窄的圖形范例。

圖5示出了鈣暈開的管腔變窄的范例。

圖6示出了在管腔橫截面處的最小閾值強度和最大閾值強度與擬合的線性模型的圖形范例。

圖7示出了在管腔橫截面的方差閾值強度與擬合的線性模型的圖形范例。

圖8用流程圖圖示了亞分辨率管腔變窄檢測和分割的示例性方法。

圖9用流程圖圖示了對可見管腔直徑進行近似的示例性方法。

圖10用流程圖圖示了檢測亞分辨率管腔變窄的示例性方法。

具體實施方式

首先參考圖1，圖示了示例性亞分辨率管腔變窄檢測和分割系統(tǒng)100。諸如計算機斷層攝影(ct)成像設備、超聲(us)、磁共振成像(mri)、正電子發(fā)射斷層攝影(pet)、單光子發(fā)射計算機斷層攝影(spect)等的成像設備102生成成像數(shù)據(jù)，例如，在成像設備102的成像膛104中的對象的感興趣區(qū)域的投影數(shù)據(jù)、原始數(shù)據(jù)等，所述感興趣區(qū)域包括具有在成像設備的分辨率內(nèi)的管腔的結(jié)構(gòu)，例如，脈管結(jié)構(gòu)。重建器106接收成像數(shù)據(jù)，并且根據(jù)成像數(shù)據(jù)來重建一幅或多幅圖像，例如，(一個或多個)重建的圖像體積。例如，ct成像設備生成心臟ct血管造影(ccta)圖像。能夠從成像設備102和/或從存儲設備108接收所重建的圖像，所述存儲設備108例如為圖片存檔與通信系統(tǒng)(pacs)、放射學信息系統(tǒng)(ris)、醫(yī)院信息系統(tǒng)(his)、電子病歷(emr)等。

血管中心線生成器110生成通過選擇的、要被分割的血管的管腔的中心線。血管中心線生成器110能夠使用本領域已知的技術(shù)和/或其它技術(shù)來生成指示的脈管結(jié)構(gòu)(例如，冠狀動脈)的中心線。血管中心線生成器110接收所重建的圖像，并且從經(jīng)自動或手動指導的所接收的圖像中根據(jù)圖像體積來確定血管管腔的中心線。

亞分辨率(sub-res.)管腔變窄檢測器112基于圖像體積和所確定的中心線來對可見管腔直徑進行近似。亞分辨率管腔變窄檢測器112利用中心線強度輪廓分析來檢測到亞分辨率變窄，例如，精確度大于所重建的圖像的可見分辨率的腔變窄的測量結(jié)果。中心線強度輪廓分析包括根據(jù)強度輪廓的變化和指示鈣暈開的增加的強度來檢測亞分辨率變窄。基于所檢測的變窄，亞分辨率管腔變窄檢測器112將近似的可見管腔直徑修正為亞分辨率確定的直徑。

在一個非限制性實例中，圖形最小切割變化分割器114基于以逐體素強度為基礎的可能性、以逐體素直徑為基礎的可能性、平滑度懲罰和優(yōu)化來在亞分辨率下分割管腔。

針對每個體素以逐體素強度為基礎的可能性是使用到中心線的橫截面中的體素強度的特性的經(jīng)訓練的模型來計算的。分割區(qū)域(例如，血管的區(qū)域)中的每個體素被分配一個值li，該值表示體素根據(jù)最接近的中心線體素的強度以及經(jīng)訓練的模型而是管腔的部分的可能性。經(jīng)訓練的模型使用最接近的中心線體素的強度ictr來計算最大強度閾值fmax(ictr)和最小強度閾值fmin(ictr)。體素是經(jīng)分割的管腔的部分的可能性能夠包括以下函數(shù)：li＝1/(wshu*((ivoxel-ictr)/((fmax(ictr)-fminictr))/2))，其中，wshu是權(quán)重，ivoxel是體素強度。

針對每個體素以逐體素直徑為基礎的可能性是使用到中心線的橫截面中的體素強度的特性的經(jīng)訓練的模型來計算的。分割區(qū)域中的每個體素被分配一個值ld，該值表示體素基于該體素距最接近的中心線點的空間距離以及亞分辨率確定的直徑而是管腔的部分的可能性。體素基于距離而是經(jīng)分割的管腔的部分的可能性能夠包括以下函數(shù)：其中，wd是權(quán)重，dvoxel是體素與最接近的中心線點之間的空間距離，d是中心線點處的亞分辨率確定的直徑。

基于針對鄰近體素p與鄰近物q的局部強度差、表面方向和局部強度差方向配對的乘積以及經(jīng)訓練的模型來計算平滑度懲罰。針對每對鄰近體素的局部強度差能夠由下式給出：其中，ip和iq是以亨氏單位(hu)為單位的鄰近體素的強度，dp,q是體素p與q之間的空間距離，fvar(ictr,p)是根據(jù)到p的最接近的中心線點處的以hu為單位的強度的方差的經(jīng)訓練的模型的函數(shù)。表面方向能夠由下式給出：其中，abs是絕對值函數(shù)，(p-q)是體素p與q之間的向量，dir(ctr,p)是到p的最接近的中心線點處的中心線的方向，θ是常數(shù)，例如在0.1至0.5的范圍內(nèi)，例如為0.3和/或其他值，局部強度差方向由下式給出：如果ip–iq>0或p(q)>m，則是1，其中m是預定閾值最大強度的常數(shù)，例如400-800hu的范圍內(nèi)的值，例如為600hu；否則為0。

基于以最小代價對構(gòu)建的圖形的切割(其將圖形分成管腔和背景)，優(yōu)化向經(jīng)分割的區(qū)域中的每個體素分配標簽，作為表示體素是否處于經(jīng)分割的管腔中的值，例如，背景。該圖形包括由節(jié)點表示的每個體素，并且兩個額外的節(jié)點表示管腔和背景，例如，非管腔。每個體素節(jié)點通過邊連接到管腔節(jié)點和背景節(jié)點。每個體素與管腔節(jié)點之間的邊的權(quán)重由li*ld給出，其中，li和ld是基于上述強度和距離而分配給體素的可能性值。每個體素與背景節(jié)點之間的邊的權(quán)重由1-li*ld給出。每個體素的節(jié)點由通過平滑度懲罰(其為局部強度差、表面方向和局部強度差方向的乘積)進行加權(quán)的邊連接到鄰近體素的節(jié)點，例如，預定的一組鄰近體素(例如，與體素鄰接的那些體素)。使用諸如ford-fulkerson、edmonds-karp、push-relabel、boykov-kolmogorov等的最小切割算法來執(zhí)行對包括節(jié)點和加權(quán)邊的構(gòu)建的圖形的切割。

圖形最小切割變化分割器114基于根據(jù)最小切割算法的經(jīng)分配的標簽在亞分辨率下分割管腔。連接到管腔的切割部分中的體素被定義為亞分辨率分割的管腔116。連接到背景的切割部分中的體素并不包括在經(jīng)分割的管腔中。亞分辨率分割的管腔能夠被顯示在顯示設備118上，被疊加在圖像體積上和/或與圖像數(shù)據(jù)108一起存儲。

血流儲備分數(shù)(frac.)分析器120能夠基于亞分辨率分割的管腔116來計算血流儲備分數(shù)。血流儲備分數(shù)分析器120使用更準確分割的結(jié)構(gòu)(例如，更精確分割的結(jié)構(gòu))來確定血流儲備分數(shù)，其包括根據(jù)減少的強度和/或鈣暈開的亞分辨率變窄。所計算的血流儲備分數(shù)能夠例如與圖像數(shù)據(jù)一起被顯示在顯示設備118上和/或被存儲在計算機存儲器中。

重建器106、血管中心線生成器110、亞分辨率管腔變窄檢測器、圖形最小切割變化分割器114和血流儲備分數(shù)分析器120由一個或多個配置的數(shù)據(jù)處理器122(例如，電子數(shù)據(jù)處理器、光學數(shù)據(jù)處理器、微處理器、數(shù)字處理器等)來適當?shù)貙嵤?。所配置的?shù)據(jù)處理器122運行被存儲在非瞬態(tài)計算機可讀存儲介質(zhì)(“存儲器”)中的計算機可讀指令。所配置的數(shù)據(jù)處理器122也可以運行由載波、信號或其他瞬態(tài)介質(zhì)承載的一個或多個計算機可讀指令。數(shù)據(jù)處理器122能夠包括諸如工作站、膝上型計算機、移動計算設備、智能電話、服務器、身體佩戴式計算設備等的計算設備124。數(shù)據(jù)處理器122能夠是分布式計算設備，例如，點對點計算設備、客戶端-服務器計算設備等。

計算設備124能夠包括顯示設備118，例如，陰極射線管(crt)、存儲管、雙穩(wěn)態(tài)顯示器、電子紙、矢量顯示器、平板顯示器、真空熒光顯示器(vf)、發(fā)光二極管(led)顯示器、電致發(fā)光顯示器(eld)、等離子體顯示面板(pdp)、液晶顯示器(lcd)、有機發(fā)光二極管顯示器(oled)、投影儀以及頭戴式顯示器等。計算設備124能夠包括一個或多個輸入設備126，例如，鍵盤、麥克風、鼠標、觸摸屏等。

參考圖2，圖示了使用半峰全寬(fwhm)規(guī)則確定的近似的可見管腔直徑的范例。在直的多平面視圖中示出左前降支冠狀管腔200。該視圖使用垂直軸上距所生成的中心線的距離(以毫米(mm)為單位)以及水平軸上距口的距離(以毫米為單位)。近似的可見管腔直徑202被指示為白線。

針對強度值對在以不同角度垂直于中心線的平面中的射線進行采樣。應用半峰全寬(fwhm)規(guī)則來確定針對每個采樣射線的近似的可見直徑。fwhm選擇直徑端點，其中，強度值是最大強度值的一半或峰值強度值的一半以及強度分布曲線的最大寬度。針對每個橫截面計算每個平面中的采樣射線上的直徑值(例如，中值)，以確定近似的可見管腔直徑。在一個實例中，針對直徑值計算樣本射線的平均值。濾波器能夠?qū)Υ怪庇谥行木€的平面上的直徑值進行平滑。

參考圖3，示意性地圖示了在可見成像分辨率下對變化的血管直徑的可視化的范例。基于所重建的圖像，血管直徑被過高估計。管腔輪廓300(例如，真實尺寸)被圖示有3mm302、1.5mm304、1.3mm306、1mm308和0.5mm310的不同輪廓或直徑。尺寸被示出在水平軸線上，中心線在4mm處，在垂直軸上強度為1.0。卷積312之后的管腔輪廓(例如如通過成像設備在分辨率下觀察的)被示為對應于真實的管腔輪廓中的每個。卷積后的管腔輪廓使用與真腔管輪廓相同的軸線。隨著直徑變窄，卷積之后的強度輪廓314的測量結(jié)果減小，例如由峰值、寬度、形狀和/或fwhm(定義如上)測量結(jié)果所指示的。例如，以大小減小的順序的管腔的峰值強度大約為1.1、1.0、0.9、0.8和0.4。半峰寬(例如卷積后輪廓的強度分別為0.55、0.5、0.45、0.4和0.2下的尺寸)過高估計了真實的管腔。

對應于不同直徑管腔中的每個圖示了管腔的二維(2d)可視化316。以白色指示的第一強度區(qū)域表示管腔318，由暗區(qū)指示的第二強度區(qū)域表示非管腔320。如所可視化的并且當管腔的直徑減小時，第一區(qū)域減小或變窄，并且類似地，在可視化中表示的管腔318的強度也減小。

參考圖4，示出了減小的強度管腔變窄的圖形范例。第一圖形400描繪沿著中心線的以hu為單位的強度值。垂直軸或y軸是以hu為單位的強度值，水平軸或x軸是沿著中心線距口的距離。圖形上的第一標繪線是沿著中心線402的強度值，第二標繪線是對強度值404的平滑。變窄406由在沿著中心線距口的距離為19mm與23mm之間處強度值的降低來指示，例如，低于最小閾值408。最小閾值408基于中心線強度值，例如，中心線強度中值。

亞分辨率管腔變窄檢測器112計算中心線輪廓的多尺度二階導數(shù)，例如，表示沿著中心線的強度值的經(jīng)平滑的線404的二階導數(shù)。在值(例如，線)低于閾值并且二階導數(shù)為正(例如，＞)的地方檢測到變窄406。

基于如參考圖3所描述的強度值的變化，第二圖形410將變窄412描繪為垂直軸上的百分比以及水平軸上沿著中心線距口的距離，與第一圖形400的標度相同。顯示對應的圖像414，其包括減小的強度區(qū)416，例如，變暗區(qū)或不太明亮的區(qū)，所述減小的強度區(qū)416對應于由亞分辨率管腔變窄檢測器112基于沿著中心線的強度值所檢測的變窄。對應的圖像414包括與第一圖形400和第二圖形410相同的水平軸。在一個實例中，亞分辨率管腔變窄檢測器112接收來自輸入設備126的輸入，所述輸入指示搜索變窄的區(qū)和/或呈現(xiàn)變窄的區(qū)。

亞分辨率管腔變窄檢測器112基于所檢測的變窄根據(jù)近似的可見管腔直徑來計算亞分辨率確定的直徑。該計算能夠包括公式其中，dsub是亞管腔變窄的直徑，dapprox是近似的可見管腔直徑，ictr是中心線的強度值，μ0是參考管腔強度值，例如，主動脈的強度值，α、β和γ是常數(shù)，例如分別為0.5、0.16和0.14。

參考圖5，示出了隨著管腔變窄的鈣暈開的范例。亞分辨率管腔變窄檢測器112基于大于一閾值的強度值來檢測由于鈣暈開而引起的變窄。第一圖形500包括表示沿著中心線的強度值的線502，以及表示用濾波器平滑的強度值的第二線504。在垂直軸或y軸上表示強度值，在水平軸或x軸上表示沿著中心線的距離。超過閾值508的強度值506被指示為在距口9mm至12mm之間。

亞分辨率管腔變窄檢測器112檢測鈣暈開的變窄，其中，值(例如，線)高于最大閾值508并且二階導數(shù)為負(例如，＜0)。上限閾值的強度值能夠使用參考值，例如，主動脈的強度。基于所檢測的由于鈣暈開而引起的變窄，亞分辨率管腔變窄檢測器112將近似的可見分辨率直徑修正為亞分辨率確定的直徑。

第二圖形510指示變窄512作為垂直軸或y軸上的管腔的百分比以及在與第一圖形500相同的水平軸上沿著中心線距口的距離。在對應于由于鈣暈開而引起的所檢測的變窄的9至12mm距離處顯示對應的圖像514，所述對應的圖像514包括增加的強度區(qū)516，例如，較亮區(qū)。對應的圖像514包括與第一圖形500和第二圖形510公用的水平軸。

參考圖6，示出了在管狀橫截面處的最小閾值強度600與模型的擬合的線性函數(shù)602以及在管腔橫截面處的最大閾值強度604與模型的擬合的線性函數(shù)606的圖形范例。圖形最小切割變化分割器114基于沿著中心線的橫截面的強度特性來生成限定用于分割管腔的體素特異性強度極限的模型。

圖形最小切割變化分割器114對沿著中心線的橫截面中的體素的特性強度進行采樣。例如，在每個橫截面處，計算要分割的區(qū)域中的強度的最小值、最大值和標準偏差的特性。圖6是散點圖，在水平軸或x軸上是以hu為單位的中心線強度，在垂直軸或y軸上是以hu為單位的最大強度和最小強度。

圖形最小切割變化分割器114基于所采樣的特性來生成模型。該模型包括中心線強度的函數(shù)。最大強度606變暗，并且將第一函數(shù)606擬合到最大強度以導出最大強度閾值極限。第一函數(shù)606能夠包括形式fmax＝aictr+b，其中，fmax是垂直于由計算的亞分辨率確定的直徑界定的中心線的平面中的體素的最大強度閾值，ictr是中心線處的體素的強度，a和b是通過線性回歸分析計算的常數(shù)。最小強度600較亮，并且第二函數(shù)602被擬合到最小強度600以導出最小強度閾值。第二函數(shù)602能夠包括形式fmin＝cictr+d，其中，imin是垂直于中心線的平面中的體素的最小強度，ictr是中心線處的體素的強度，c和d是通過線性回歸分析計算的常數(shù)。該模型不同于全局值，例如從主動脈強度導出的值或從中心線的值導出的值(例如，中心線值的平均值)。

參考圖7，示出了管腔橫截面的標繪的變化強度700和擬合的線性函數(shù)702的圖形范例。該模型包括采樣的中心線橫截面強度的方差的擬合函數(shù)702。在水平軸或x軸上表示中心線強度，以hu為單位，在垂直軸或y軸上表示方差(例如，標準偏差)，以hu為單位。第三函數(shù)702被擬合到采樣的方差上。第三函數(shù)702能夠包括形式fvar＝eictr+f，其中，ivar是垂直于中心線的平面中體素的強度的方差，ictr是中心線處的體素的強度，e和f是通過線性回歸分析計算的常數(shù)。

參考圖8，用流程圖圖示了亞分辨率管腔變窄檢測和分割的示例性方法。在800處，根據(jù)諸如ccta數(shù)據(jù)的圖像數(shù)據(jù)和血管中心線檢測管腔變窄，以計算亞分辨率確定的管腔直徑。使用fwhm規(guī)則和/或其他方法根據(jù)圖像數(shù)據(jù)和中心線對可見管腔直徑進行近似。在沿著中心線的橫截面處對強度輪廓進行采樣。如參考圖3、圖4和圖5所描述地分析中心線強度輪廓。如參考圖4和圖5所描述地檢測變窄，其包括由低于最小閾值的強度以及高于最大閾值的強度而檢測的變窄?；谒鶛z測的變窄和近似的可見管腔直徑，計算針對沿著中心線的橫截面的亞分辨率確定的直徑。

在806處，根據(jù)中心線強度輪廓和經(jīng)訓練的模型的函數(shù)導出逐體素強度極限。參考圖6和圖7描述經(jīng)訓練的模型的函數(shù)的范例。

在810處，基于以逐體素強度為基礎的可能性、以逐體素直徑為基礎的可能性、平滑度懲罰和優(yōu)化來構(gòu)建圖形。通過最小切割分割算法將圖形切成兩部分。一個部分限定經(jīng)分割的管腔812，一個部分是背景，例如，不經(jīng)分割的管腔。經(jīng)分割的管腔812被表示為三維空間模型，例如，冠狀動脈管腔。

參考圖9，用流程圖圖示了對可見管腔直徑進行近似的示例性方法。在900處，接收具有生成的血管的中心線的圖像體積。在902處，沿著中心線對橫截面進行采樣。每個橫截面以橫截面的平面的射線進行采樣，并且射線處于k個角度，其中，k是＞0的整數(shù)。

在904處，針對每個橫截面，使用fwhm規(guī)則來計算下限閾值。在906處，針對每個橫截面，使用常數(shù)(例如，2.3倍的下限閾值)來計算上限閾值。

在908處，從每個橫截面中消除低于下限閾值或上限閾值的異常值，例如處于一角度的射線。在910處，針對橫截面中的每個剩余角度估計直徑。

在912處，根據(jù)針對每個橫截面的剩余角度處的估計的直徑來計算近似的可見直徑，例如，中值直徑、平均值直徑等。

在914處，能夠?qū)ρ刂行木€的橫截面處的近似的可見直徑進行濾波。所述濾波能夠使用平滑濾波器，例如在橫截面上求平均，這對從橫截面到橫截面的近似的直徑的變化進行了平滑。

參考圖10，用流程圖圖示了檢測亞分辨率管腔變窄的示例性方法。在1000處，接收具有生成的血管中心線和估計的管腔的近似的可見直徑的圖像體積。

在1002處，確定強度上限閾值和強度下限閾值。所述上限閾值能夠根據(jù)主動脈管腔的強度來確定。所述下限閾值能夠根據(jù)中心線的強度來確定。

在1004處，計算中心線輪廓的多尺度二階導數(shù)。

在1006處，檢測變窄。大于最大閾值并且具有為負的二階導數(shù)(例如，局部最大值)的強度值識別由于沿著中心線的橫截面處的鈣暈開而引起的變窄。小于最小閾值并且具有為正的二階導數(shù)(例如，局部最小值)的強度值識別基于強度在沿著中心線的橫截面處所檢測的變窄。

在1008處，基于所檢測的變窄來修正近似的可見直徑，以計算每個橫截面處的亞分辨率確定的直徑，所檢測的變窄是基于根據(jù)所減少的強度的變窄或由于鈣暈開而引起的變窄的。

以上可以以計算機可讀指令的方式實施，所述計算機可讀指令被編碼或嵌入在計算機可讀存儲介質(zhì)上，所述計算機可讀指令當由(一個或多個)計算機處理器執(zhí)行時，令(一個或多個)處理器執(zhí)行所描述的動作。額外地或備選地，計算機可讀指令中的至少一個由信號、載波或其他瞬態(tài)介質(zhì)承載。

已經(jīng)參考優(yōu)選實施例描述了本發(fā)明。他人在閱讀和理解前面的具體描述的情況下可以想到修改和替代。本文旨在將本發(fā)明解釋為包括所有這樣的修改和替代，只要它們落入權(quán)利要求書及其等價方案的范圍內(nèi)。