ajout

function ajout
global im
global xi
global yi


I=get(im,'CData');

[x,y]=ginput(1);

[m,ind]=min(sqrt((xi-x).*(xi-x)+(yi-y).*(yi-y)));
xi1=[xi(1:ind) ;x ;xi(ind+1:end)];
yi1=[yi(1:ind); y; yi(ind+1:end)];
xi2=[xi(1:ind-1) ;x ;xi(ind:end)];
yi2=[yi(1:ind-1); y; yi(ind:end)];

d1=sum(sqrt((xi1(1:end)-[xi1(2:end);xi1(1)]).^2+(yi1(1:end)-[yi1(2:end);yi1(1)]).^2));
d2=sum(sqrt((xi2(1:end)-[xi2(2:end);xi1(1)]).^2+(yi2(1:end)-[yi2(2:end);yi2(1)]).^2));
if d1>d2
    xi=xi2;yi=yi2;
else
    xi=xi1;yi=yi1;
end
Mask=roipoly(I,xi,yi);

I(:,:,3)=255*Mask;

set(im,'CData',I);