itkImageAlgorithm.h

Go to the documentation of this file.

/*=========================================================================
 *
 *  Copyright Insight Software Consortium
 *
 *  Licensed under the Apache License, Version 2.0 (the "License");
 *  you may not use this file except in compliance with the License.
 *  You may obtain a copy of the License at
 *
 *         http://www.apache.org/licenses/LICENSE-2.0.txt
 *
 *  Unless required by applicable law or agreed to in writing, software
 *  distributed under the License is distributed on an "AS IS" BASIS,
 *  WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
 *  See the License for the specific language governing permissions and
 *  limitations under the License.
 *
 *=========================================================================*/
#ifndef itkImageAlgorithm_h
#define itkImageAlgorithm_h

#include "itkImageRegionIterator.h"

#ifdef ITK_HAS_STLTR1_TYPE_TRAITS
#  include <type_traits>
#elif defined ITK_HAS_STLTR1_TR1_TYPE_TRAITS
#  include <tr1/type_traits>
#else
#  include "itkIsSame.h"
#endif

#ifdef ITK_HAS_CPP11_TYPETRAITS
#  define ITK_STD_TR1_NAMESPACE std
#else
#  define ITK_STD_TR1_NAMESPACE std::tr1
#endif

namespace itk
{

template <typename TPixelType, unsigned int VImageDimension > class VectorImage;


struct ImageAlgorithm
{

#if defined(ITK_HAS_STLTR1_TR1_TYPE_TRAITS) || defined(ITK_HAS_STLTR1_TYPE_TRAITS)
    typedef ITK_STD_TR1_NAMESPACE::true_type  TrueType;
    typedef ITK_STD_TR1_NAMESPACE::false_type FalseType;
#else
    typedef itk::TrueType  TrueType;
    typedef itk::FalseType FalseType;
#endif

  template<typename InputImageType, typename OutputImageType>
  static void Copy( const InputImageType *inImage, OutputImageType *outImage,
                    const typename InputImageType::RegionType &inRegion,
                    const typename OutputImageType::RegionType &outRegion )
  {
    ImageAlgorithm::DispatchedCopy( inImage, outImage, inRegion, outRegion );
  }

  template<typename TPixel1, typename TPixel2, unsigned int VImageDimension>
  static void Copy( const Image<TPixel1, VImageDimension> * inImage,
                               Image<TPixel2, VImageDimension> * outImage,
                               const typename Image<TPixel1, VImageDimension>::RegionType &inRegion,
                               const typename Image<TPixel2, VImageDimension>::RegionType &outRegion )
  {
    typedef Image<TPixel1, VImageDimension> _ImageType1;
    typedef Image<TPixel2, VImageDimension> _ImageType2;
    ImageAlgorithm::DispatchedCopy( inImage, outImage, inRegion, outRegion
#if defined(ITK_HAS_STLTR1_TR1_TYPE_TRAITS) || defined(ITK_HAS_STLTR1_TYPE_TRAITS)
                                   , ITK_STD_TR1_NAMESPACE::is_convertible<typename _ImageType1::PixelType,
                                   typename _ImageType2::PixelType>()
#else
                                    // note the above trait is
                                    // primarily used to get a better
                                    // error message
                                    , TrueType()
#endif
      );
  }

  template<typename TPixel1, typename TPixel2, unsigned int VImageDimension>
  static void Copy( const VectorImage<TPixel1, VImageDimension> * inImage,
                               VectorImage<TPixel2, VImageDimension> * outImage,
                               const typename VectorImage<TPixel1, VImageDimension>::RegionType &inRegion,
                               const typename VectorImage<TPixel2, VImageDimension>::RegionType &outRegion )
  {
    typedef VectorImage<TPixel1, VImageDimension> _ImageType1;
    typedef VectorImage<TPixel2, VImageDimension> _ImageType2;
    ImageAlgorithm::DispatchedCopy( inImage, outImage, inRegion, outRegion
#if defined(ITK_HAS_STLTR1_TR1_TYPE_TRAITS) || defined(ITK_HAS_STLTR1_TYPE_TRAITS)
                                   , ITK_STD_TR1_NAMESPACE::is_convertible<typename _ImageType1::PixelType,
                                   typename _ImageType2::PixelType>()
#else
                                    , TrueType()
#endif
      );
  }

  template<typename InputImageType, typename OutputImageType>
  static typename OutputImageType::RegionType
  EnlargeRegionOverBox(const typename InputImageType::RegionType & inputRegion,
                       const InputImageType* inputImage,
                       const OutputImageType* outputImage);

private:

  template<typename InputImageType, typename OutputImageType>
  static void DispatchedCopy( const InputImageType *inImage, OutputImageType *outImage,
                              const typename InputImageType::RegionType &inRegion,
                              const typename OutputImageType::RegionType &outRegion, TrueType isSpecialized );

  template<typename InputImageType, typename OutputImageType>
  static void DispatchedCopy( const InputImageType *inImage, OutputImageType *outImage,
                              const typename InputImageType::RegionType &inRegion,
                              const typename OutputImageType::RegionType &outRegion, FalseType isSpecialized = FalseType() );


  template <typename TImageType>
  struct PixelSize
  {
    static size_t Get( const TImageType *)
    {
      return 1;
    }
  };

  template <typename TPixelType, unsigned int VImageDimension>
  struct PixelSize< VectorImage<TPixelType, VImageDimension> >
  {
    typedef VectorImage<TPixelType, VImageDimension> ImageType;
    static size_t Get( const  ImageType * i )
    {
      const size_t vectorLength = ImageType::AccessorFunctorType::GetVectorLength(i);
      return vectorLength;
    }
  };

  template<typename TInputType, typename TOutputType>
  struct StaticCast
    : public std::unary_function<TInputType,TOutputType>
  {
    TOutputType operator() (const TInputType i) { return static_cast<TOutputType>(i); }
  };


  template<typename TType>
  static TType* CopyHelper(const TType *first, const TType *last, TType *result)
    {
      // Note: On some MS compilers the following may generate a
      // warning. Please include itkMacro.h before <algorithm> or
      // another stl header to avoid.
      return std::copy(first, last, result);
    }

  template<typename TInputType, typename TOutputType>
  static TOutputType* CopyHelper(const TInputType *first, const TInputType *last, TOutputType *result)
    {
      return std::transform(first, last, result, StaticCast<TInputType,TOutputType>());
    }

};
} // end namespace itk


#ifndef ITK_MANUAL_INSTANTIATION
#include "itkImageAlgorithm.hxx"
#endif


#endif //itkImageAlgorithm_h