matlab 分水岭法,matlab中的分水岭算法

这是使用分水岭分割图像的一种方法.你可以做更多的事情(例如,如果他们尚未完成胞质分裂,则融合两个细胞核的细胞),但下面的步骤应该给你一个第一个想法.

(1)确定细胞 – 背景阈值,细胞核阈值

%# read image

img = imread('http://i.stack.imgur.com/nFDkX.png');

%# normalize to 0...1

imgN = double(img-min(img(:)))/(max(img(:)-min(img(:))));

th1=graythresh(imgN);

th2 = graythresh(imgN(imgN>th1));

cellMsk = imgN>th1;

nucMsk = imgN>th2;

figure,imshow(cellMsk+nucMsk,[])

(2)平滑原始图像(以避免过度分割)并将原子核强加为最小值

[xx,yy]=ndgrid(-5:5,-5:5);

gf = exp((-xx.^2-yy.^2)/20);

filtImg = conv2(imgN,gf,'same');

figure,imshow(filtImg,[])

filtImgM = imimposemin(-filtImg,nucMsk);

(3)分水岭,掩模细胞和显示器

ws = watershed(filtImgM);

ws(~cellMsk) = 0;

lblImg = bwlabel(ws);

figure,imshow(label2rgb(lblImg,'jet','k','shuffle'));

(4)现在,您可以在带标签的图像上使用REGIONPROPS来提取所需的统计数据.