function bmap(im)
  clf, zoom on
  
  colormap([bone(128); autumn(128)])  
  d = double(im)/2;
  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) + 128;

  image(d), axis image, axis off