#ifdef _USE_ITK_

#include "bbitkImageRegion.h"
#include "bbitkPackage.h"

namespace bbitk
{

  BBTK_BLACK_BOX_IMPLEMENTATION(ImageRegionCreator,
				bbtk::AtomicBlackBox);

  BBTK_ADD_BLACK_BOX_TO_PACKAGE(itk,ImageRegionCreator);


  void ImageRegionCreator::DoIt()
  {
    const std::vector<long>& index = bbGetInputIndex();
    const std::vector<long>& size = bbGetInputSize();
    //   	   std::cout<< "isize="<<index.size()<<std::endl;
    // 	   std::cout<< "ssize="<<size.size()<<std::endl;
    unsigned long maxs = index.size() > size.size() ? index.size() : size.size();
    switch (maxs) 
      {
      case 2 : DoIt<2>(); break;
      case 3 : DoIt<3>(); break;
      case 4 : DoIt<4>(); break;
      default : bbtkError("ImageRegionCreator : cannot build a region of dimension "<<maxs);
      }
     

  }

  template <unsigned int Dimension>
  void ImageRegionCreator::DoIt()
  {
    std::vector<long> index = bbGetInputIndex();
    std::vector<long> size = bbGetInputSize();
    //    	   std::cout<< "isize="<<index.size()<<std::endl;
    //    	   std::cout<< "ssize="<<size.size()<<std::endl;
 
    int ds = index.size() - size.size();
    if (ds<0) 
      {
	for (int i=0;i<-ds;++i) 
 index.push_back(0);
      }
    else 
      {
	for (int i=0;i<ds;++i) size.push_back(0);
      }

    typename itk::ImageRegion<Dimension>::IndexType I;
    typename itk::ImageRegion<Dimension>::SizeType S;

    std::vector<long>::const_iterator ii,si;
    int d(0);
    for (ii=index.begin(), si=size.begin(); 
	 ii!=index.end(); ++ii, ++si, ++d)
      {
	I[d] = *ii;
	S[d] = *si;
	//  std::cout << "o="<<*ii<<" d="<<*si<<std::endl;
      }


    bbSetOutputOut ( itk::ImageRegion<Dimension>(I,S) );
  }

}
// eo namespace

#endif