最近学习图像处理,发现一个特别酷的操作,如标题所述。
原图
切换灰度图
img2gray()函数就可以实现
对灰度图进行闭运算消除噪声
seDisk=strel('disk',5,0); %构造圆盘
imgClose=imclose(imgGray,seDisk);
大津阈值分割
level=graythresh(imgClose); %大津阈值法
imgBw=imbinarize(imgClose,level);
获得异常细胞
需要查看细胞核和细胞的面积比,然后找到不正常的比例,就能确定为异常细胞,这边处理的代码有点长,我就假设第11个细胞异常是已知的(其实是通过代码计算出来的)。
获得异常细胞
Lc=bwlabel(1-imgBw);
ccc=bwconncomp(1-imgBw);
areaAllc=regionprops(Lc,'Area');
imgCell=regionprops(Lc,'Image');
figure;
imshow(imgCell(11).Image);
在原图中标记细胞
之前都是常规操作,这部分才是我想写这博文的关键。
先看结果,标记出了原图中细胞的位置(虽然可以通过中心坐标值来简单实现,但是我没想到击中击不中的方法真的能找到这个点!!)
imgJi=bwhitmiss(1-imgBw,imgCell(11).Image,1-imgCell(11).Image); %击中击不中
se=strel('disk',10);
imgPen=imdilate(imgJi,se);%膨胀算法
注意击中击不中处理,默认背景是白色,前景是黑色,所以代码里面才会出现1-imgBw。
代码中间跳了不少[笑哭],有问题可以留言讨论,我尽量及时回答。