颜色12个特征值数据新.m

clear all;clc;
A=[];
for i=1:52 %如果文件夹有7张图，可更改
ii=num2str(i);
I=imread(['E:\新建文件夹\A\',ii,'.jpg']); %读取当前文件夹，' '可以变成其他文件
R=I(:,:,1); %红色分量
G=I(:,:,2); %绿色分量
B=I(:,:,3); %蓝色分量
R=double(R); G=double(G); B=double(B); %利用double()函数将变量类型转为double型
Ravg1=mean2(R); %红色分量均值
Gavg1=mean2(G); %绿色分量均值
Bavg1=mean2(B); %蓝色分量均值
Rstd1=std(std(R)); %红色分量的方差
Gstd1= std(std(G)); %绿色分量的方差
Bstd1=std(std(B)); %蓝色分量的方差
%计算I分量
I=(R+G+B);
%计算S分量
tmp1=min(min(R,G),B);
tmp2=R+G+B;
tmp2(tmp2==0)=eps;
S=1-3.*tmp1./tmp2;
%计算H分量
tmp1=0.5*((R-G)+(R-B));
tmp2=sqrt((R-G).^2+(R-B).*(G-B));
theta=acos(tmp1./(tmp2+eps));
H=theta;
H(B>G)=2*pi-H(B>G);
H=H/(2*pi);
H(S==0)=0;
hsi=cat(3,H,S,I);
Havg1=mean2(H); %红色分量均值
Savg1=mean2(S); %绿色分量均值
Iavg1=mean2(I); %蓝色分量均值
Hstd1=std(std(H)); %红色分量的方差
Sstd1= std(std(S)); %绿色分量的方差
Istd1=std(std(I)); %蓝色分量的方差

a=[Ravg1,Gavg1,Bavg1,Rstd1,Gstd1,Bstd1,Havg1,Savg1,Iavg1,Hstd1,Sstd1,Istd1];
A=[A;a];
xlswrite('E:\新建文件夹\A\color颜色均值12个特征值.xls',A);
end