function [day_rgI, plant_rgI, day]=loci

load('dxc_930p.dat');
sensor=dxc_930p(:,2:4);

load('plant_107.dat');
plant=plant_107;

fid=fopen('cie.dat');
daylight=fscanf(fid, '%g', [58 107]);
status=fclose(fid);

fold=daylight*sensor;

nor=fold(18,:);
balance=[nor(2)/nor(1), 1, nor(2)/nor(3)];

for i=1:58
  fold(i,:)=fold(i,:).*balance;
end

div=fold(:,1)+fold(:,2)+fold(:,3);
d=100*fold./[div, div, div];
day=[d(:,3), d(:,2), d(:,1)];

I=zeros(58,3);

for i=1:58
  if (d(i,1)>d(i,2) & d(i,1)>d(i,3))
    I(i,:)=d(i,:)./[d(i,1),d(i,1),d(i,1)];
  elseif (d(i,2)>d(i,1) & d(i,2)>d(i,3)) 
    I(i,:)=d(i,:)./[d(i,2),d(i,2),d(i,2)];
  elseif (d(i,3)>d(i,1) & d(i,3)>d(i,2)) 
    I(i,:)=d(i,:)./[d(i,3),d(i,3),d(i,3)];
  else
    I(i,:)=[1,1,1];
  end
end

I=255*sqrt(sum((I.*I)'));

day_rgI=[d(:,3), d(:,2), I'];

fold=daylight*(sensor.*[plant, plant, plant]);


for i=1:58
  fold(i,:)=fold(i,:).*balance;
end

div=fold(:,1)+fold(:,2)+fold(:,3);
d=100*fold./[div, div, div];

I=zeros(58,1);

plant_rgI=[d(:,3), d(:,2), I];