itkAdvancedLinearInterpolateImageFunction.h

Go to the documentation of this file.

/*=========================================================================
 *
 *  Copyright UMC Utrecht and contributors
 *
 *  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 __itkAdvancedLinearInterpolateImageFunction_h
#define __itkAdvancedLinearInterpolateImageFunction_h
 
#include "itkLinearInterpolateImageFunction.h"
 
namespace itk
{
template< class TInputImage, class TCoordRep = double >
class AdvancedLinearInterpolateImageFunction :
  public LinearInterpolateImageFunction< TInputImage, TCoordRep >
{
public:
 
  typedef AdvancedLinearInterpolateImageFunction                   Self;
  typedef LinearInterpolateImageFunction< TInputImage, TCoordRep > Superclass;
  typedef SmartPointer< Self >                                     Pointer;
  typedef SmartPointer< const Self >                               ConstPointer;
 
  itkTypeMacro( AdvancedLinearInterpolateImageFunction, LinearInterpolateImageFunction );
 
  itkNewMacro( Self );
 
  typedef typename Superclass::OutputType OutputType;
 
  typedef typename Superclass::InputImageType  InputImageType;
  typedef typename InputImageType::SpacingType InputImageSpacingType;
 
  typedef typename Superclass::InputPixelType InputPixelType;
 
  typedef typename Superclass::RealType RealType;
 
  itkStaticConstMacro( ImageDimension, unsigned int, Superclass::ImageDimension );
 
  typedef typename Superclass::IndexType IndexType;
 
  typedef typename Superclass::ContinuousIndexType ContinuousIndexType;
  typedef typename ContinuousIndexType::ValueType  ContinuousIndexValueType;
 
  typedef CovariantVector< OutputType,
    itkGetStaticConstMacro( ImageDimension ) >        CovariantVectorType;
 
  CovariantVectorType EvaluateDerivativeAtContinuousIndex(
    const ContinuousIndexType & x ) const;
 
  void EvaluateValueAndDerivativeAtContinuousIndex(
    const ContinuousIndexType & x,
    OutputType & value,
    CovariantVectorType & deriv ) const
  {
    return this->EvaluateValueAndDerivativeOptimized(
      Dispatch< ImageDimension >(), x, value, deriv );
  }
 
 
protected:
 
  AdvancedLinearInterpolateImageFunction();
  ~AdvancedLinearInterpolateImageFunction() override{}
 
private:
 
  AdvancedLinearInterpolateImageFunction( const Self & ); // purposely not implemented
  void operator=( const Self & );                         // purposely not implemented
 
  struct DispatchBase {};
  template< unsigned int >
  struct Dispatch : public DispatchBase {};
 
  inline void EvaluateValueAndDerivativeOptimized(
    const Dispatch< 2 > &,
    const ContinuousIndexType & x,
    OutputType & value,
    CovariantVectorType & deriv ) const;
 
  inline void EvaluateValueAndDerivativeOptimized(
    const Dispatch< 3 > &,
    const ContinuousIndexType & x,
    OutputType & value,
    CovariantVectorType & deriv ) const;
 
  inline void EvaluateValueAndDerivativeOptimized(
    const DispatchBase &,
    const ContinuousIndexType & x,
    OutputType & value,
    CovariantVectorType & deriv ) const
  {
    return this->EvaluateValueAndDerivativeUnOptimized( x, value, deriv );
  }
 
 
  inline void EvaluateValueAndDerivativeUnOptimized(
    const ContinuousIndexType & x,
    OutputType & value,
    CovariantVectorType & deriv ) const
  {
    itkExceptionMacro( << "ERROR: EvaluateValueAndDerivativeAtContinuousIndex() "
                       << "is not implemented for this dimension ("
                       << ImageDimension << ")." );
  }
 
 
};
 
} // end namespace itk
 
#ifndef ITK_MANUAL_INSTANTIATION
#include "itkAdvancedLinearInterpolateImageFunction.hxx"
#endif
 
#endif