0001 function msk = getInteractiveMask( grayimg, markers )
0002 fh = figure;
0003 grayimg = grayimg - min(min(grayimg));
0004 grayimg = grayimg / max(max(grayimg));
0005 msk = false(size(grayimg));
0006 rgb = cat(3,grayimg/2,grayimg/3,grayimg);
0007
0008 while 1
0009 tmp =rgb;
0010 tmp(:,:,1)=max(max(tmp(:,:,1)))*msk;
0011 imagesc(tmp);
0012 hold on;
0013 plot(markers,'*y');
0014 nb = roipoly;
0015 if ~any(nb)
0016 break;
0017 end
0018 msk = msk | nb;
0019 end
0020
0021
0022
0023 close(fh);
0024 return