---

/home/xchen/xia/src/image_processing/demoDilation.cc

#include <Morphology.H>

using namespace xchen;

void dilate(BitMap const& image)
{
  vector<iE2> cntr_pnts;
  Morphology :: ComputeRegionCentralPoints(image, cntr_pnts);

  cout << "\n\nimage is\n" << image << endl;
  cout << "cntr_pnts are\n";
  copy(cntr_pnts.begin(), cntr_pnts.end(), ostream_iterator<iE2>(cout, " "));
  cout << endl;
}


int main()
{
  BitMap image(20,20);
  fill_n(image.Iterate(), 20*20, 0);

  for(int i=0;i<20; i++) image[0][i] = image[19][i] = image[i][0] = image[i][19] = 1;
  dilate(image);

  for(int i=0;i<20; i++)  image[i][i] = 1; 
  dilate(image);

  for(int i=0;i<20; i++) image[i][19-i] = 1;
  dilate(image);

  BitMap big_image(1000,1000);
  fill_n(big_image.Iterate(), 1000*1000, 0);
  for(int i=0;i<1000; i++)
  {
    big_image[0][i] = big_image[999][i] = big_image[i][0] = big_image[i][999] = 1;
  }
  dilate(big_image);

  return 0;
}

---