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

%***********************平滑***********************************%
h?=?fspecial（‘average‘3）;?%3*3均值濾波模板
I?=?imfilter（I?h?‘replicate‘）;?%重復填充邊界
%figure（2）;imshow（I）;

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

%**********************技巧性處理******************************%
%該處理主要是針對圖片四周的紅細胞的處理，便于下一步的填充處理，避免將邊緣的紅細胞去除
[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‘）;%對處理后的紅細胞進行填充處理
%figure（3）subplot（132）imshow（C）;
D?=?D?-?C;
D?=?imfill（D8‘holes‘）;
%figure（3）subplot（133）imshow（D）;

%*********************標記區域*********************************%
[LN]?=?bwlabel（D8）;?????%連通標記
s?=?regionprops（L‘Area‘）;%將連通標記結果對應連通區域計算面積
F?=?bwareaopen（D1308）;??%根據上面統計的面積，結合圖片，將面積小于130個像素點的區域刪除。注意，這里存在一個大氣泡和白細胞無法去除
[LN]?=?bwlabel（F8）;?????%對處理后的圖像重新進行連通標記
%figureimshow（F）;
B?=?bwboundaries（F）;??????%疑似紅細胞邊界提取

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

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

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

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

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

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

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

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