clear;clc;close?all;
%%%%%%%%%%%%%%%%%?讀取原圖，并處理成二值圖像?%%%%%%%%%%%%%%%%%%%%%%%
%I=imread（‘梨2.jpg‘）;
%I=imread（‘桃子2.jpg‘）;
%I=imread（‘蘋果2.jpg‘）;
%I=imread（‘香蕉.jpg‘）;
%I=imread（‘青椒.jpg‘）;
I=imread（‘test4.jpg‘）;

I2=rgb2gray（I）;
BW=im2bw（I20.9）;
figure（1）subplot（131）imshow（I）title（‘原始圖像‘）;
subplot（132）imshow（I2）title（‘灰度圖像‘）;
subplot（133）imshow（BW）title（‘二值圖像‘）;

%%%%%%%%%%%%%%%%?進行邊緣檢測?得到不連續的圖形邊界?%%%%%%%%%%%%%%%%%%%%%%
%得到各個圖形的連續邊界
SE=strel（‘rectangle‘[40?30]）;??%?結構定義
J2=imopen（BWSE）;????????????%?開啟運算
figure（2）imshow（J2）title（‘對二值圖像進行開運算‘）;

SE=strel（‘square‘5）;?%?采用方形結構元素進行腐蝕
J=imerode（~J2SE）;
BW2=（~J2）-J;????????%?檢測邊緣
figure（3）imshow（BW2）title（‘邊緣檢測‘）;

%填充了已有的檢測的連續形狀邊界
B?=?imfill（BW2‘holes‘）;%圖像填充
B?=?bwmorph（B‘remove‘）;%移除內部像素
figure（4）imshow（B）title（‘提取出的邊界圖像‘）;

%將不同的圖形進行分別標記，num表示連接的圖形對象的個數
[Labelnum]?=?bwlabel（B8）;

%%%%%%%%%%%%%%%%%%%計算各個圖形單元邊界像素點數%%%%%%%%%%%%%%%%%%%%%%%%%%%
????for?i?=?1?:?num
????????Premeter（i）?=?0;
????end

????[rowcol]?=?size（Label）;
????for?i?=?1?:?row
????????for?j?=?1?:?col
????????????if（Label（ij）?>?0）
????????????????Premeter（Label（ij））?=?Premeter（Label（ij））?+?1;%計算標記后的各塊圖形邊界中像素的個數的總數
????????????end
????????end
????end

%%%%%%%%%%%%%%%%%%%?計算各個圖形單元的面積%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
FilledLabel?=?imfill（Label‘holes‘）;??%填充標記過的邊界線中間圍成的圖形區域
figureimshow（FilledLabel）title（‘標記過并被填充的結果‘）;
for?i?=?1?:?num
????Area（i）?=?0;
end

[rowcol]?=?size（FilledLabel）;
for?i?=?1?:?row
????for?j?=?1?:?col
????????if（FilledLabel（ij）?>?0）
????????????Area（FilledLabel（ij））?=?Area（FilledLabel（ij））?+?1;???%通過統計像素點個數的方式來求各形狀的面積
????????end
????end
end

%%%%%%%%%%%%%%%%%%?計算各個圖形單元的圓度?%%%%%%%%%%%%%%%%%%%%%%%%%%
for?i?=?1?:?num?????
????Ecllipseratio（i）?=?4*pi*Area（i）/Premeter（i）^2;
end

%%%%%%%%%%%%%%%%%%%%?計算各個圖像的色度?%%%%%%%%%%%%%%%%%%%%%

HSV?=?rgb2hsv（I）;???%轉換為HSV，為后面的色度元素的提取做準備

[rowcol]?=?size（FilledLabel）;???%統計填充后的圖形中各塊圖形所含像素的個數的多少
MeanHue?=?zeros（1num）;
????for?i?=?1?:?num
????????Hue?=?zeros（Area（i）1）;
????????nPoint?=?0;
????????for?j?=?1?:?row
????????????for?k?=?1?:?col
????????????????if（FilledLabel（jk）?==?i）
????????????????????nPoint?=?nPoint?+?1;
????????????????????Hue（nPoint1）?=?HSV（jk1）;
????????????????end
????????????end
????????end
????????
????????Hue（:i）?=?sort（Hue（:1））;
????????for?j?=?floor（nPoint*0.1）?:?floor（nPoint*0.9）
????????????MeanHue（i）?=?MeanHue（i）?+?Hue（j1）;
????????end
????????MeanHue（i）?=?MeanHue（i）?/?（0.8*nPoint）;??%計算出平均的色度值
????end
????
%%%%%%%%%%%%%%%%%%%%%%%%?各種水果識別?%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%識別梨，在二維特征空間對各個圖像進行類別區分
pear=0;result1=0;
for?i=1:num
????if（MeanHue（i）<0.125）??%判斷各個圖形中平均色度值小于0.125的為梨
????????pear=i;
????????result1=1;
????end
end
%對分出來的梨構建相應的圖像掩膜，并用對原圖的亮度圖像進行掩膜操作
pearHSV=HSV;
????????for?j?=?1?:?row
????????????for?k?=?1?:?col
????????????????if（FilledLabel（jk）?~=pear）
???????????????????pearHSV（jk3）=0;
????????????

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2014-05-04?11:04??大設計\
?????文件????????5191??2014-05-03?21:10??大設計\final.m
?????文件???????23733??2014-05-03?10:27??大設計\test1.jpg
?????文件???????20824??2014-05-02?22:48??大設計\test2.jpg
?????文件???????18189??2014-05-03?09:11??大設計\test3.jpg
?????文件???????81877??2014-05-03?10:35??大設計\test4.jpg
?????文件?????1401421??2014-05-04?10:54??大設計\大設計報告.docx
?????文件??????923905??2014-05-04?10:56??大設計\大設計報告.pdf
?????文件????????5776??2014-05-02?14:52??大設計\無.png
?????文件???????17213??2014-05-03?10:10??大設計\桃子1.jpg
?????文件???????14903??2014-05-03?09:37??大設計\桃子2.jpg
?????文件????????6996??2012-10-25?08:59??大設計\梨1.jpg
?????文件???????22336??2014-05-03?10:15??大設計\梨2.jpg
?????文件????????9356??2014-05-02?16:41??大設計\蘋果1.jpg
?????文件???????19645??2014-05-02?16:44??大設計\蘋果2.jpg
?????文件???????12359??2014-05-02?21:00??大設計\青椒.jpg
?????文件???????35508??2012-10-24?18:25??大設計\香蕉.jpg