itkImageToImageFilterDetail.h

Go to the documentation of this file.

/*=========================================================================

  Program:   Insight Segmentation & Registration Toolkit
  Module:    $RCSfile: itkImageToImageFilterDetail.h,v $
  Language:  C++
  Date:      $Date: 2005/05/13 19:44:21 $
  Version:   $Revision: 1.8 $

  Copyright (c) Insight Software Consortium. All rights reserved.
  See ITKCopyright.txt or http://www.itk.org/HTML/Copyright.htm for details.

  Portions of this code are covered under the VTK copyright.
  See VTKCopyright.txt or http://www.kitware.com/VTKCopyright.htm for details.

     This software is distributed WITHOUT ANY WARRANTY; without even 
     the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR 
     PURPOSE.  See the above copyright notices for more information.

=========================================================================*/
#ifndef __itkImageToImageFilterDetail_h
#define __itkImageToImageFilterDetail_h

#include "itkImageRegion.h"
#include "itkSmartPointer.h"

namespace itk
{

namespace ImageToImageFilterDetail
{

  struct DispatchBase {};

  template <bool>
  struct BooleanDispatch {};

  template <int>
  struct IntDispatch : public DispatchBase {};

  template <unsigned int>
  struct UnsignedIntDispatch : public DispatchBase {};

  template <bool B1, bool B2>
  struct BinaryBooleanDispatch
  {

    typedef BooleanDispatch<B1> FirstType;
    typedef BooleanDispatch<B2> SecondType;
  };

  template <int D1, int D2>
  struct BinaryIntDispatch
  {

    typedef IntDispatch<D1> FirstType;
    typedef IntDispatch<D2> SecondType;
  };

  template <unsigned int D1, unsigned int D2>
  struct BinaryUnsignedIntDispatch : public DispatchBase
  {

    typedef UnsignedIntDispatch<D1> FirstType;
    typedef UnsignedIntDispatch<D2> SecondType;

    typedef IntDispatch<(D1 > D2) - (D1 < D2)> ComparisonType;
    typedef IntDispatch<0>  FirstEqualsSecondType;  
    typedef IntDispatch<1>  FirstGreaterThanSecondType;  
    typedef IntDispatch<-1> FirstLessThanSecondType;  
  };

  template <unsigned int D1, unsigned int D2>
  void ImageToImageFilterDefaultCopyRegion(const typename
                  BinaryUnsignedIntDispatch<D1, D2>::FirstEqualsSecondType &,
                  ImageRegion<D1> &destRegion,
                  const ImageRegion<D2> &srcRegion)
  {
    destRegion = srcRegion;
  }

  template <unsigned int D1, unsigned int D2>
  void ImageToImageFilterDefaultCopyRegion(const typename
                  BinaryUnsignedIntDispatch<D1, D2>::FirstLessThanSecondType &,
                  ImageRegion<D1> &destRegion,
                  const ImageRegion<D2> &srcRegion)
  {
    // Source dimension is greater than the destination dimension, copy the
    // first part of the source into the destination
    unsigned int dim;
    Index<D1> destIndex;
    Size<D1>  destSize;
    const Index<D2> &srcIndex = srcRegion.GetIndex();
    const Size<D2> &srcSize = srcRegion.GetSize();

    // copy what we can
    for (dim=0; dim < D1; ++dim)
      {
      destIndex[dim] = srcIndex[dim];
      destSize[dim] = srcSize[dim];
      }
    
    destRegion.SetIndex(destIndex);
    destRegion.SetSize(destSize);
  }
  
  template <unsigned int D1, unsigned int D2>
  void ImageToImageFilterDefaultCopyRegion(const typename
               BinaryUnsignedIntDispatch<D1, D2>::FirstGreaterThanSecondType &,
               ImageRegion<D1> &destRegion,
               const ImageRegion<D2> &srcRegion)
  {
    // Source dimension is less than the destination dimension, copy source
    // into the first part of the destination and set zeros elsewhere.
    unsigned int dim;
    Index<D1> destIndex;
    Size<D1>  destSize;
    const Index<D2> &srcIndex = srcRegion.GetIndex();
    const Size<D2> &srcSize = srcRegion.GetSize();

    // copy what we can
    for (dim=0; dim < D2; ++dim)
      {
      destIndex[dim] = srcIndex[dim];
      destSize[dim] = srcSize[dim];
      }
    // fill in the rest of the dimensions with zero/one
    for (; dim < D1; ++dim)
      {
      destIndex[dim] = 0;
      destSize[dim] = 1;
      }
    
    destRegion.SetIndex(destIndex);
    destRegion.SetSize(destSize);
  }


  template <unsigned int D1, unsigned int D2>
    class ITK_EXPORT ImageRegionCopier 
  {
  public:
    virtual void operator() (ImageRegion<D1> &destRegion,
                            const ImageRegion<D2> &srcRegion) const
    {
      typedef typename BinaryUnsignedIntDispatch<D1, D2>::ComparisonType ComparisonType;
      ImageToImageFilterDefaultCopyRegion<D1, D2>(
                      ComparisonType(),
                      destRegion, srcRegion);
    }
    virtual ~ImageRegionCopier() {}
  };


  template<unsigned int D1, unsigned int D2>
  std::ostream & operator<<(std::ostream &os, const ImageRegionCopier<D1, D2> &
                            copier)
  {
    os << "ImageRegionCopier: "
       << typeid(ImageRegionCopier<D1, D2>).name() << std::endl;
    return os;
  }

  template<unsigned int D1, unsigned int D2>
  bool operator!=(const ImageRegionCopier<D1, D2> &c1,
                  const ImageRegionCopier<D1, D2> &c2)
  {
    return &c1 != &c2;
  }

} // end of namespace ImageToImageFilterDetail
    
} // end namespace itk

//#ifndef ITK_MANUAL_INSTANTIATION
//#include "itkImageToImageFilterDetail.txx"
//#endif

#endif

---