clear;
d11=imread('football.jpg');
%d1=double(i);
[m n r]=size(d11);
z=max(max(max(d11)));
d1=d11;
min1=double(min(min(min(d1))));
%mx1=max(max(d1(:,:,1)));
%mx2=max(max(d1(:,:,2)));
%mx3=max(max(d1(:,:,3)));
d(:,:,1)=double(d1(:,:,1));
d(:,:,2)=double(d1(:,:,2));
d(:,:,3)=double(d1(:,:,3));
I=(d(:,:,1)+d(:,:,2)+d(:,:,3))/3;
s=1-min1./I;
rmg= d(:,:,1)-d(:,:,2);
rmb= d(:,:,1)-d(:,:,3);
gmb= d(:,:,2)-d(:,:,3);
r1=rmg+rmb;
r2=rmg.*rmg+rmb.*gmb;
for t=1:m
for t1=1:n
if(r2(t,t1)==0)
r2(t,t1)=1;
end
end
end
r3=sqrt(r2);
a=(0.5*r1./r3);
a1=a/(max(max(a)));
%for t=1:m
% for t1=1:n
% a=0.5*r1./r3;
% if(a>1.0)
% a=1.0;
% end
% end
%end
tht=acosd(a1);
for t=1:m
for t1=1:n
if( d(t,t1,3)>d(t,t1,2))
H=360-tht;
else
H=tht;
end
end
end
fil(:,:,1)=H;
fil(:,:,2)=s;
fil(:,:,3)=I;
imshow(uint8(H));
figure;
imshow(uint8(I));
figure;
imshow((s));
figure;
imshow(fil);