• 沒有找到結果。

Cell compettition 的四種型式與 cost function

在新一代的 2D cell competition 演算法中,cell competition 的過程被視為是許多個 snakes 同 時形變的現象。一個 snake 包含了一個或數個 cells。且 cell competition 的型式也由過去的三種增 加為四種。Type I 的 cell competition 指得是一個僅包含一個 cell 的 snake 想要合併另一個同類型 的 snake。Type II 指得則是一個僅包含一個 cell 的 snake 想要自一個包含多個 cells 的 snake 中脫 離。Type III 則是指一個僅包含一個 cell 的 snake 與一個包含多個 cells 的 snake 間競爭一個 cell 的機制。Type IV 則是指兩個包含多個 cells 的 snake 間競爭一個 cell 的機制。此四種競爭型式的 圖例可參考附件二。

嶄新的 cost function 是此一新的 cell competition 的演算法主要特色。此一新的 cost function 包含兩個主要成份。一個區域一致性,另一個是邊緣強度。令ℜ 表示 cell competiton 第 i 個 iterationi 的區域一致性,它定義為

Θ

圖 V-B-1:第 i 個 iteration 中 snake 形變的一個剪影(snapshot),其中粗實線、細實線與細虛線 分別表示 snakes、cells 與 elementary cells 的邊緣。Snakes sipsqi 之間的邊緣上的小圈圈則是定 義為節點,它們乃是 snakes、cells 與 elementary cells 的邊緣在兩 snakes 間的共同邊界的交點。

Cost function 的第二個主要成份—邊緣強度—主要是量化所有相鄰 snakes 的共同邊界的相 似性。為了定義邊緣強度,圖 V-B-1 呈現第 i 個 iteration 中 snake 形變的一個剪影(snapshot),

其中粗實線、細實線與細虛線分別表示 snakes、cells 與 elementary cells 的邊緣。Snakes sipsqi 之間的邊緣上的小圈圈則是定義為結點,它們乃是 snakes、cells 與 elementary cells 的邊緣在兩 snakes 間的共同邊界的交點。定義兩 snakes 間的共同邊界上,兩個節點間的邊緣像素點集合為 一 boundary segment。令 表示兩個 elementary cells,

vertex

elementary-cell boundary cell boundary

snake boundary

i

sp

i

sq

圖 V-B-2:相鄰兩節點

υ

i2

υ

i3所定義的 boundary segment BS(ej1,ej2),其中粗實線為兩 snakes 間的共同邊界,而細虛線則為 elementary cells 的邊緣。

每一段 boundary segment 皆賦與一個 cost 以描述其所在的 edge profile 以及分享此一 boundary segment 之 elementary cells 的區域相似度。所為 edge profile 指得是此 boundary segment 附近的灰階值的變化。在新一代的 cell competition 演算法中,我們考慮了兩類的 edge profiles,

為方便起見,稱其為 classes I 和 II。假設圖 V-B-2 中的線 L 垂直於BS(ej1,ej2)且與BS(ej1,ej2)交 於像素點

υ

LBS(ej1,ej2)。圖 V-B-3 描繪了線 L 上,classes I 和 II 兩種 edge profiles。此兩種 edge profiles 的主要差異在於 class II 呈現出快速變化的 overshoot 或 undershoot。約略而言,class I edge profile 可見於一般的組織邊緣,而 class II edge profile 有時可見諸於兩個腫瘤的中隔。

每一段 boundary segment 的 cost 是以三個雙尾 Kolmogorov-Smirnov (KS) tests 的 p-values 定 義之。第一個 KS test(表示為 KS(e ,j1 e ))j2 ,檢定 elementary cells e 和j1 e 中灰階值分佈的相似j2 程度。其目的在於量化 class I edge 的強度。為了捕捉 class II edge 的 overshoot(或 undershoot)

的特徵,針對每一個

υ

LBS(ej1,ej2),令

ζ

L

υ

L在線 L 上最近NL個相鄰點中具有最大灰階值的

b 之三個k KS tests 的 p-values 中最小的一個,而

{

j

}

綜合而言,在新一代的 2D cell competition 演算法中,其第 i 個 iteration 的 total cost function(以 Ε 表示之)定義為: i

N 和c N 分別為第一個 iteration 開始時 cells 與 boundary segments 的數目。2D cell coompetitionbs 演算法採取 steepest gradient descent 的方式進行 cost 的最小化。在第 i 個 iteration 中,本演算法 選取會導致 cost 降低最多的 snake 形變進行之。

圖 V-B-3:延著線 L 之(a) Class I edge profile;(b) class II edge profile。

應用實例

在附件二中,我們對於新一代的 2D cell competition 演算法的效能有詳細的分析,在此我們 僅呈現在四張超音波乳房腫瘤影像應用的成果。圖 V-B-4 為原始影像,而圖 V-B-5 則為以本演 算法所得到的腫瘤邊緣。其結果顯示,所求得的腫瘤邊緣與視覺所定義的相當吻合,值得一提 的是第一張的腫瘤中隔也被順利偵測出。

(a) (b)

(c) (d)

圖 V-B-4: 四張原始超音波乳房腫瘤影像。

υL

L

NL

υL

L

(a)

(b) NL

(a) (b)

(c) (d)

圖 V-B-5:以新一代 2D cell competition 演算法所求得的乳房腫瘤邊緣。