01用Distance Mapping對Clustered細胞計數

Introduction:
血球計數女蜈。如果血球細胞數量異常滔悉,可能說明你患有疾病型凳。人工計數是通過顯微鏡觀察血細胞切片沼沈,然后數數流酬。自動計數要面對的一個問題是:分割或者說如果處理clustered cells 和 粘連細胞币厕。clustered cells 和其相鄰的細胞接觸到一起,導致只有一小部分不相連的細胞可以看得清楚芽腾。兩個相鄰的細胞旦装,他們之間的對比度是很小的。因此想要將粘連細胞分離開比分離細胞和背景難得多摊滔。(在二值圖里阴绢,背景是黑色的,細胞是白色的)

Discuss Methods
A艰躺、 Edge Detection
Canny or Sobel 算子的邊緣檢測都是利用圖像intensity(亮度) 梯度強度來檢測物體的邊緣呻袭。但是在兩個粘連細胞里,他們之間的對比度很小腺兴,所以很難分別畫出他們內部的輪廓左电。

B、 Hough Transform (不是很懂)
霍夫變換是用來檢測一些形狀含长。通常券腔,霍夫變換需要邊緣檢測作為前處理。然后霍夫變換estimates the parameters of the shape from its edge拘泞,然后在參數空間進行voting纷纫。他可以用來檢測部分粘連的圓形細胞。對于圖片要求比較高陪腌。

C辱魁、 Model Based Contour Tracing
利用細胞形狀的先驗知識來促進輪廓描繪。它使用一個二階多項式模型來自動檢測圖片诗鸭。但是這個方法很垃圾基本只能檢測獨立的細胞染簇。下面是原話。

this method can trace the contours of fully visible or minimally occluded cells only

D强岸、 Morphological Image Processing
用一個預先定義好的(也是根據先驗知識確定的)幾何形狀 (structuring elements)锻弓,來進行形態(tài)學操作。一般是用disk元素蝌箍。這個方法的缺點是當粘連細胞很嚴重時青灼,腐蝕就不是很準確了。

E妓盲、 Template-Based Matching
利用一個模板圖片在圖片里去搜索想要的目標杂拨。在識別相似的目標時,模板會與圖片的不同區(qū)域相匹配悯衬。通常我們是使用模板和圖片之間的互相關cross-correlation來進行匹配的弹沽。因為cross-correlation具有抗噪性,所以它的魯棒性還可以。另一方面策橘,模板匹配在處理object occlusion時不是那么好炸渡,因此對于高聚類的細胞效果不好。

F役纹、 Distance Transform
一張二值圖片的距離變換標記了每個目標像素到離它最近的背景像素的距離偶摔。因為目標中心centroid是距離背景像素最遠的,因此距離變換被用來可以檢測細胞中心點促脉。為了分割二值圖像中的overlapping目標,距離變換和分水嶺變換通常結合起來去分割overlapping目標策州。首先瘸味,計算出二值圖像的距離變換。接著够挂,距離變化圖distance map 被complemented so that 細胞中心點就對應了局部最小值旁仿。最后,分水嶺變換被用到complemented distance map上孽糖,把細胞看作disjoint集水盆地枯冈。但是,分水嶺算法會產生過分割办悟,因為圖片里有許多局部最小值尘奏,就會產生許多的盆地。

First, the distance transform of the binary thresholded microscopy image is computed. Next, this distance map is complemented so that the cell centroids now correspond to the local minima. Finally, the watershed transform is applied to the complemented distance map in order to identify the cells as disjoint catchment basins. However, the watershed method produces a severe oversegmentation of the image, i.e., many small catchment basins are produced due to many local minima in the input image.

(a)overlapping disks (b)edge detection (c)erosion (d)Hough circle transform(detected circles are overlaid on the parameter space (e)template matching) (f)distance transform

The Proposed Method : Distance Mapping
在距離變換不被細胞的overlapping和occlusion影響的情況下病蛉,基于距離變換的方法對clustered cells的分割炫加、計數效果不錯。然而铺然, 當噪聲存在或者distance map里有局部不規(guī)律時俗孝,距離變換會造成過分割,從而導致很多False positive魄健。這種過分割在面對有噪聲的high degree of cell overlapping的顯微圖像時赋铝,顯得更加明顯。
我們提出了一種Distance Mapping方法沽瘦,把基于模式匹配的模板用到distance transform map上革骨,來最小化過分割。這種模板匹配有助于檢測細胞中心點其垄,并且抵消了distance map中假的local maxima作用苛蒲。這些local maxima造成了過分割和計數中的False positive值。因此绿满,我們的模板匹配有助于減弱過分割現象臂外,(由于定冠詞一直是the,我也分不清楚到底是哪個模板匹配)。具體步驟見下漏健。

A嚎货、 Steps of Counting
1)Pre-Processing:
首先,原始original圖片Io用Contrast-Limited Adaptive Histogram Equal(CLAHE)處理從而得到增強enhanced圖片Ie蔫浆。CLAHE增強了圖片的對比度殖属,修正了光照的不均勻。
2)Global Threshold:
接下來瓦盛,把細胞和背景分割開來洗显。用Otsu閾值分割,得到二值圖像B原环。
3)Distance Transform:
然后挠唆,計算二值圖像的distance transform。這一步把二值圖像B變成了a distance map D嘱吗,D代表了每個cell pixel距離它最近的背景像素的distance玄组。

Template Generation:(a) a circular disk (b) the corresponding distance transform map. The distance map is used as the template image

4)Template Generation:
為了運用模板匹配,模板圖片T是由一個圓盤的距離變化產生的谒麦。圓盤的半徑是根據細胞的平均size來選擇的俄讹。
5)Template Matching:
模板圖片T和distance transform map D匹配來識別細胞的中心點。 模板匹配 is performed by T和D之間的normalized cross-correlation绕德。

S代表了相似矩陣患膛,矩陣里的值都是[0,1]。0代表了毫不相似迁匠,1代表了完美匹配剩瓶。Trow 和 Tcol 代表了模板圖片的行數與列數,Ed 和 Et 代表了 the energy of the distance map 和 the energy of the template image城丧。為了檢測邊界上的partial cells延曙,我們使用symmetric border handling during the cross-correlation computation。
Energy of the picture

6)Watershed Transform:
最后亡哄,相似矩陣S is complemented 并且所有的背景像素 are set to -Inf枝缔。然后把分水嶺變換用到這個 complemented 相似矩陣-S上。(蚊惯?愿卸?)這樣就把相似矩陣S分成了separate disjoint regions L。每個區(qū)域都被標記和計數截型,這樣就能得到細胞數目趴荸。
B、 The Algorithm
distance mapping 算法見table I宦焦。


distance mapping algorithm

結果
A发钝、 Cell-Background Separation
因為細胞和背景之間的對比度很高顿涣,所以Otsu閾值分割方法可以很準確地分割背景和目標。但是不能分開clustered cells酝豪。

B涛碑、 Cell Counting
下面是distance mapping 的示意圖。


(a) original pic (b) thresholded pic (c)similarity matrix from the template matching on the distance map (d) the watershed segmentation for cell counting
Note:high (60%) overlap probability

然后是不同方法的對比孵淘,紅色部分是false positive值蒲障。

image.png

特別是distance transform 和 distance mapping。
(a) distance transform (b) distance mapping

最后編輯于
?著作權歸作者所有,轉載或內容合作請聯系作者
  • 序言:七十年代末瘫证,一起剝皮案震驚了整個濱河市揉阎,隨后出現的幾起案子,更是在濱河造成了極大的恐慌背捌,老刑警劉巖余黎,帶你破解...
    沈念sama閱讀 216,651評論 6 501
  • 序言:濱河連續(xù)發(fā)生了三起死亡事件,死亡現場離奇詭異载萌,居然都是意外死亡,警方通過查閱死者的電腦和手機巡扇,發(fā)現死者居然都...
    沈念sama閱讀 92,468評論 3 392
  • 文/潘曉璐 我一進店門扭仁,熙熙樓的掌柜王于貴愁眉苦臉地迎上來,“玉大人厅翔,你說我怎么就攤上這事乖坠。” “怎么了刀闷?”我有些...
    開封第一講書人閱讀 162,931評論 0 353
  • 文/不壞的土叔 我叫張陵熊泵,是天一觀的道長。 經常有香客問我甸昏,道長顽分,這世上最難降的妖魔是什么? 我笑而不...
    開封第一講書人閱讀 58,218評論 1 292
  • 正文 為了忘掉前任施蜜,我火速辦了婚禮卒蘸,結果婚禮上,老公的妹妹穿的比我還像新娘翻默。我一直安慰自己缸沃,他們只是感情好,可當我...
    茶點故事閱讀 67,234評論 6 388
  • 文/花漫 我一把揭開白布修械。 她就那樣靜靜地躺著趾牧,像睡著了一般。 火紅的嫁衣襯著肌膚如雪肯污。 梳的紋絲不亂的頭發(fā)上翘单,一...
    開封第一講書人閱讀 51,198評論 1 299
  • 那天吨枉,我揣著相機與錄音,去河邊找鬼县恕。 笑死东羹,一個胖子當著我的面吹牛,可吹牛的內容都是我干的忠烛。 我是一名探鬼主播属提,決...
    沈念sama閱讀 40,084評論 3 418
  • 文/蒼蘭香墨 我猛地睜開眼,長吁一口氣:“原來是場噩夢啊……” “哼美尸!你這毒婦竟也來了冤议?” 一聲冷哼從身側響起,我...
    開封第一講書人閱讀 38,926評論 0 274
  • 序言:老撾萬榮一對情侶失蹤师坎,失蹤者是張志新(化名)和其女友劉穎恕酸,沒想到半個月后,有當地人在樹林里發(fā)現了一具尸體胯陋,經...
    沈念sama閱讀 45,341評論 1 311
  • 正文 獨居荒郊野嶺守林人離奇死亡蕊温,尸身上長有42處帶血的膿包…… 初始之章·張勛 以下內容為張勛視角 年9月15日...
    茶點故事閱讀 37,563評論 2 333
  • 正文 我和宋清朗相戀三年,在試婚紗的時候發(fā)現自己被綠了遏乔。 大學時的朋友給我發(fā)了我未婚夫和他白月光在一起吃飯的照片义矛。...
    茶點故事閱讀 39,731評論 1 348
  • 序言:一個原本活蹦亂跳的男人離奇死亡,死狀恐怖盟萨,靈堂內的尸體忽然破棺而出凉翻,到底是詐尸還是另有隱情,我是刑警寧澤捻激,帶...
    沈念sama閱讀 35,430評論 5 343
  • 正文 年R本政府宣布制轰,位于F島的核電站,受9級特大地震影響胞谭,放射性物質發(fā)生泄漏垃杖。R本人自食惡果不足惜,卻給世界環(huán)境...
    茶點故事閱讀 41,036評論 3 326
  • 文/蒙蒙 一韭赘、第九天 我趴在偏房一處隱蔽的房頂上張望缩滨。 院中可真熱鬧,春花似錦泉瞻、人聲如沸脉漏。這莊子的主人今日做“春日...
    開封第一講書人閱讀 31,676評論 0 22
  • 文/蒼蘭香墨 我抬頭看了看天上的太陽侧巨。三九已至,卻和暖如春鞭达,著一層夾襖步出監(jiān)牢的瞬間司忱,已是汗流浹背皇忿。 一陣腳步聲響...
    開封第一講書人閱讀 32,829評論 1 269
  • 我被黑心中介騙來泰國打工, 沒想到剛下飛機就差點兒被人妖公主榨干…… 1. 我叫王不留坦仍,地道東北人鳍烁。 一個月前我還...
    沈念sama閱讀 47,743評論 2 368
  • 正文 我出身青樓,卻偏偏與公主長得像繁扎,于是被迫代替她去往敵國和親幔荒。 傳聞我的和親對象是個殘疾皇子,可洞房花燭夜當晚...
    茶點故事閱讀 44,629評論 2 354

推薦閱讀更多精彩內容