function [Data,Data_gt_rgb,label]=HSI_read(HSI,Data_gt,ClassColor,Class_Numble)
%%
[row,colum,spec]=size(HSI);
feats{1}=(reshape(HSI,row*colum,spec))';
feats{1} = feats{1}./repmat(sqrt(sum(feats{1}.^2,1)),size(feats{1},1),1);%按列归一化
feats{2}=(reshape(Data_gt,row*colum,1))';
Data=feats{1};
%% Ground_truth map
label=cell(Class_Numble,1); % 各类样本点
label(1)={find(Data_gt==0)};
each_class_number=zeros(1,Class_Numble);% 各类样本数
Data_gt_rgb=zeros(3,row*colum);
for i=1:Class_Numble
index=find(Data_gt==i);
label(i+1)={index}; % 每类样本
each_class_number(i)=length(index);
Data_gt_rgb(:,index)=repmat(ClassColor(:,i),1,length(index));
end