%% 計算得到 HSI三個分量
a=double(imread('D:\王浩\標準板1-1旋轉.bmp'));  %% 讀入圖像
R=a(:,:,1);  %% 得到R分量
G=a(:,:,2);  %% 得到G分量 
B=a(:,:,3);  %% 得到B分量
[M,N]=size(a); %% 得到圖像大小
N=N/3;  %% 得到圖像的列，因為24為圖則需要除以3
H=zeros(M,N);  %% 定義存儲H分量矩陣
S=zeros(M,N);  %% 定義存儲S分量矩陣
I=zeros(M,N);  %% 定義存儲I分量矩陣
temp=zeros(1,3); %% 定義臨時變量
I=(R+G+B)/(3*255); %% 計算得到I分量
H=acos((2*R-G-B)./(0.0000001+2*sqrt((R-G).*(R-G)+(R-B).*(G-B))))*180/pi;  %%初步計算H分量，這里將分母加0.0000001為了避免分母為零的情況
tic
%% 計算S 以及H
for i=1:M
    for j=1:N
        temp(1)=R(i,j);
        temp(2)=G(i,j);
        temp(3)=B(i,j);
        S(i,j)=1-3*min(temp)/(sum(temp)+0.0000001);
        if(G(i,j)<B(i,j))
            H(i,j)=360-H(i,j);
        end
    end
end