i=imread('pic8.jpg');
imshow(i);
subplot(4,4,1)
imshow(i)
%%
iycbcr=rgb2ycbcr(i);
subplot(4,4,2)
imshow(iycbcr)
y=iycbcr(:,:,1);
cb=iycbcr(:,:,2);
cr=iycbcr(:,:,3);
ccb=cb.^2;
subplot(4,4,3)
imshow(ccb)
ccr=(255-cr).^2;
subplot(4,4,4)
imshow(ccr)
cbcr=ccb./cr;
subplot(4,4,5)
imshow(cbcr)
%%
igray=rgb2gray(i);
subplot(4,4,6)
imshow(igray)
g=1./3;
l=g*ccb;
m=g*ccr;
n=g*cbcr;
%%
eyemapchr=l+m+n;
subplot(4,4,7)
imshow(eyemapchr)
J=histeq(eyemapchr);
subplot(4,4,8)
imshow(J)
SE=strel('disk',15,8);
o=imdilate(igray,SE);
p=1+imerode(igray,SE);
eyemaplum=o./p;
subplot(4,4,9)
imshow(eyemaplum)
cc=and(J,eyemaplum);
subplot(4,4,10)
imshow(cc)