function bmap(im)
  clf
  
  colormap([bone(256); autumn(256)])  
  d = double(im);
  image(d), axis image, axis off

  coords = getline('closed');
  [x,y] = meshgrid(1:size(im,1),1:size(im,2));
  z = inpolygon(x,y,coords(:,1), coords(:,2));
  d(z) = d(z) + 256;

  image(d), axis image, axis off
