comhsb.m

来自「这是一个关于matlab提取图象中心位置的原代码」· M 代码 · 共 27 行

%% 计算得到 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