// This program gives various examples of the usage of the IPRunLengthC and
// IPLinearAverageC filters. 


#include "amma/NumImage.hh"
#include "amma/Image/IPLinearAverage.hh"
#include "amma/Image/IPRunLength.hh"
#include "amma/DP/ComposeSE.hh"


int main()
{
  // create image and initialise it: all 1s except 1 x 7
  NumImageC<RealT> in(7,12); in.Fill(1); in[2][4]=7;
  NumImageC<RealT> out;

  // this line, + the associated header file, are requires to implement the
  // following "operator>>" style of applying image processing operators
  using namespace DPComposeSE;

  // filter with 2x3 mask run-length; default behaviour is: preserve original
  // image size, pad with black, shift 1/2 pixel up/left for even mask
  in >> IPRunLengthC<RealT>(2,3) >> out; 
  cout <<in.Rectangle() << "\n" << out << "\nrunlength orig\n\n";

  // exactly the same filter, but using the "Apply()" method (i.e. don't need
  // the "DPComposeSE" stuff)
  IPRunLengthC<RealT> f(2,3);
  out = f.Apply(in);
  cout <<in.Rectangle() << "\n" << out << "\nrunlength orig2\n\n";
  
  // same filter, but eroding to avoid border area problems
  in >> IPRunLengthC<RealT>(2,3,Erode) >> out; 
  cout <<in.Rectangle() << "\n" << out << "\nrunlength erode\n\n";

  // filter with 2x3 mask averaging; default behaviour is: preserve size, pad
  // with black & rescale border area, shift 1/2 pixel up/left for even mask
  // size (here we have an even no. of rows)
  in >> IPLinearAverageC<RealT>(2,3) >> out;
  cout <<in.Rectangle() << "\n" << out << "\nav orig TL\n\n";
  
  // same basic filter, but shift 1/2 pixel down/right for even mask size
  in >> IPLinearAverageC<RealT>(2,3,Original,PBRescale,ShiftDR) >> out;
  cout <<in.Rectangle() << "\n" << out << "\nav orig BR\n\n";
  
}