%*************************************************************%
%作者：楊呂鵬
%時(shí)間：2018/7/21
%實(shí)現(xiàn)功能：對血液圖片進(jìn)行處理，提取圖片中的紅細(xì)胞的形態(tài)特征，實(shí)現(xiàn)紅細(xì)胞計(jì)數(shù)。正確率大于95%。
close?all;clear?all;clc;
%*************************************************************%
I?=?imread（‘blood.jpg‘）;%讀取血液圖像
I?=?rgb2gray（I）;????????%將彩圖轉(zhuǎn)為灰度圖
%figure（1）;imshow（I）;

%***********************平滑***********************************%
h?=?fspecial（‘a(chǎn)verage‘3）;?%3*3均值濾波模板
I?=?imfilter（I?h?‘replicate‘）;?%重復(fù)填充邊界
%figure（2）;imshow（I）;

%*********************對原圖像進(jìn)行二值化處理********************%
t?=?graythresh（I）;????????%計(jì)算閾值t
C?=?im2bw（It）;???????????%根據(jù)閾值二值化圖像
%C?=?not（C）;
%figure（3）;subplot（131）imshow（C）;

%**********************技巧性處理******************************%
%該處理主要是針對圖片四周的紅細(xì)胞的處理，便于下一步的填充處理，避免將邊緣的紅細(xì)胞去除
[rc]?=?size（C）;
C（1:）?=?1;
%C（r:）?=?1;
C（:1）?=?1;
C（:c）?=?1;?????????
C（2:）?=?0;
%C（r-1:）?=?0;
C（:2）?=?0;
C（:c-1）?=?0;?
C?=?not（C）;

%**********************填充處理********************************%
D=imfill（C8‘holes‘）;%對處理后的紅細(xì)胞進(jìn)行填充處理
%figure（3）subplot（132）imshow（C）;
D?=?D?-?C;
D?=?imfill（D8‘holes‘）;
%figure（3）subplot（133）imshow（D）;

%*********************標(biāo)記區(qū)域*********************************%
[LN]?=?bwlabel（D8）;?????%連通標(biāo)記
s?=?regionprops（L‘Area‘）;%將連通標(biāo)記結(jié)果對應(yīng)連通區(qū)域計(jì)算面積
F?=?bwareaopen（D1308）;??%根據(jù)上面統(tǒng)計(jì)的面積，結(jié)合圖片，將面積小于130個(gè)像素點(diǎn)的區(qū)域刪除。注意，這里存在一個(gè)大氣泡和白細(xì)胞無法去除
[LN]?=?bwlabel（F8）;?????%對處理后的圖像重新進(jìn)行連通標(biāo)記
%figureimshow（F）;
B?=?bwboundaries（F）;??????%疑似紅細(xì)胞邊界提取

%*******************目標(biāo)的坐標(biāo)提取*****************************%
for?g?=?1:length（B）
????boundary?=?B{g};??????%依次提取每個(gè)目標(biāo)的邊界
????Y?=?boundary（:1）;
????X?=?boundary（:2）;
????Ymax?=?max（Y）;????????%疑似紅細(xì)胞上邊界
????Ymin?=?min（Y）;????????%疑似紅細(xì)胞下邊界
????Xmax?=?max（X）;????????%疑似紅細(xì)胞左邊界
????Xmin?=?min（X）;????????%疑似紅細(xì)胞右邊界
????
????x（g）?=?floor（（Xmax?+?Xmin）/2）;%記錄疑似紅細(xì)胞的橫坐標(biāo)
????y（g）?=?floor（（Ymax?+?Ymin）/2）;%記錄疑似紅細(xì)胞的縱坐標(biāo)?
end

%***********************目標(biāo)提取******************************%
figureimshow（I）;
for?i?=?1:g?
????rectangle（‘Position‘[x（i）-21?y（i）-21?42?42]...
??????????????‘edgecolor‘‘r‘‘LineWidth‘2）;
end

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----

?????文件?????131438??2018-07-21?13:00??紅細(xì)胞提取\blood.jpg

?????文件???????2451??2018-07-21?17:27??紅細(xì)胞提取\blood.m

?????目錄??????????0??2018-07-21?17:27??紅細(xì)胞提取

-----------?---------??----------?-----??----

???????????????133889????????????????????3