close?all;?
clear?all;
clc;??
[FileNamePathName]?=?uigetfile（‘*.*‘）;%彈出對話框得到打開的路徑
Image=imread（strcat（PathNameFileName））;%打開得到的路徑下的文件?
[MNO]?=?size（Image）;?
[hsv]?=?rgb2hsv（Image）;?
H?=?h;?S?=?s;?V?=?v;?
h?=?h*360;???
%將hsv空間非等間隔量化：?
%??h量化成8級；
%??s量化成3級；
%??v量化成3級；?
for?i?=?1:M?
??for?j?=?1:N??
??????if?h（ij）<=20||h（ij）>315?????????????
??????????H（ij）?=?0;?????????
??????end??????????
??????if?h（ij）<=40&&h（ij）>20?????????????
??????????H（ij）?=?1;?????????
??????end
??????if?h（ij）<=75&&h（ij）>40?????????????
??????????H（ij）?=?2;?????????
??????end
??????if?h（ij）<=155&&h（ij）>75?????????????
??????????H（ij）?=?3;?????????
??????end
??????if?h（ij）<=190&&h（ij）>155?????????????
??????????H（ij）?=?4;????????
??????end
??????if?h（ij）<=270&&h（ij）>190?????????????
??????????H（ij）?=?5;?????????
??????end
??????if?h（ij）<=295&&h（ij）>270?????????????
??????????H（ij）?=?6;?????????
??????end
??????if?h（ij）<=315&&h（ij）>295?????????????
??????????H（ij）?=?7;?????????
??????end
??end
end
for?i?=?1:M?????
????for?j?=?1:N??????????
????????if?s（ij）<=0.2&&s（ij）>0?????????????
????????????S（ij）?=?0;?????????
????????end
????????if?s（ij）<=0.7&&s（ij）>0.2?????????????
????????????S（ij）?=?1;?????????
????????end
????????if?s（ij）<=1&&s（ij）>0.7?????????????
????????????S（ij）?=?2;?????????
????????end
????end
end
for?i?=?1:M?????
????for?j?=?1:N??????????
????????if?v（ij）<=0.2&&v（ij）>0????????????
????????????V（ij）?=?0;????????????
????????end
????????if?v（ij）<=0.7&&v（ij）>0.2????????????
????????????V（ij）?=?1;?????????
????????end
????????if?v（ij）<=1&&v（ij）>0.7????????????
????????????V（ij）?=?2;?????????
????????end
????end
end
%將三個顏色分量合成為一維特征向量：L?=?H*Qs*Qv+S*Qv+V；QsQv分別是?S和V的量化級數?L取值范圍[071]??
%取Qs?=?3;?Qv?=?3?
for??i?=?1:M?????
????for?j?=?1:N?????????
????????L（ij）?=?H（ij）*9+S（ij）*3+V（ij）;?????
????end
end
%計算L的直方圖?
for?i?=?0:71??????
????Hist（i+1）?=?size（find（L==i）1）;?
end
T0=Hist;???
%?循環讀入圖像?
A=dir（‘D:\graduate?design\test_pictures\*.jpg‘）;?
for?k=1:size（A）?????
????B=strcat（‘D:\graduate?design\test_pictures\‘A（k）.name）;????
????Image=imread（B）;???
[MNO]?=?size（Image）;???
%?計算每一幅圖像的顏色直方圖???
[hsv]?=?rgb2hsv（Image）;?
H?=?h;?S?=?s;?V?=?v;
h?=?h*360;???
%將hsv空間非等間隔量化：?
%??h量化成8級；
%??s量化成3級；
%??v量化成3級；?
for?i?=?1:M
????for?j?=?1:N