Image1 = zeros(lines,samples,3);
for i = 1:lines
for j = 1:samples
if(HSI(i,j,1)>=0 && HSI(i,j,1)<120)
Image1(i,j,3) = round(HSI(i,j,3)*(1-HSI(i,j,2)));
Image1(i,j,1) = round(HSI(i,j,3)*(1+HSI(i,j,2)*cos(pi*HSI(i,j,1)/180)/cos(pi*(60-HSI(i,j,1))/180)));
Image1(i,j,2) = round(3*HSI(i,j,3)- Image1(i,j,1)-Image1(i,j,3));
end
if(HSI(i,j,1)>=120 && HSI(i,j,1)<240)
Image1(i,j,1) = round(HSI(i,j,3)*(1-HSI(i,j,2