itkComplexBSplineInterpolateImageFunction.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 itkComplexBSplineInterpolateImageFunction_h
#define itkComplexBSplineInterpolateImageFunction_h

#include "itkBSplineInterpolateImageFunction.h"
#include "itkComplexToRealImageFilter.h"
#include "itkComplexToImaginaryImageFilter.h"

namespace itk
{
template< typename TImageType, typename TCoordRep = double, typename TCoefficientType = double >
class ITK_TEMPLATE_EXPORT ComplexBSplineInterpolateImageFunction:
  public InterpolateImageFunction< TImageType, TCoordRep >
{
public:
  ITK_DISALLOW_COPY_AND_ASSIGN(ComplexBSplineInterpolateImageFunction);

  using Self = ComplexBSplineInterpolateImageFunction;

  using Superclass = InterpolateImageFunction< TImageType, TCoordRep >;

  using Pointer = SmartPointer< Self >;

  using ConstPointer = SmartPointer< const Self >;

  itkTypeMacro(ComplexBSplineInterpolateImageFunction, InterpolateImageFunction);

  itkNewMacro(Self);

  static constexpr unsigned int ImageDimension = Superclass::ImageDimension;

  using OutputType = typename Superclass::OutputType;

  using InputImageType = typename Superclass::InputImageType;

  using IndexType = typename Superclass::IndexType;

  using ContinuousIndexType = typename Superclass::ContinuousIndexType;

  using PointType = typename Superclass::PointType;

  using InternalImageType = Image< double, Self::ImageDimension >;

  using RealFilterType = ComplexToRealImageFilter< InputImageType, InternalImageType >;
  using ImaginaryFilterType = ComplexToImaginaryImageFilter< InputImageType, InternalImageType >;

  using InterpolatorType = BSplineInterpolateImageFunction< InternalImageType, TCoordRep, TCoefficientType >;

  OutputType EvaluateAtContinuousIndex(const ContinuousIndexType & index) const override;

/*  using CovariantVectorType = CovariantVector< OutputType, Self::ImageDimension  >;

  CovariantVectorType EvaluateDerivative( const PointType & point ) const
  {
   ContinuousIndexType index;
   this->GetInputImage()->TransformPhysicalPointToContinuousIndex( point, index );
   return ( this->EvaluateDerivativeAtContinuousIndex( index ) );
  }

  CovariantVectorType EvaluateDerivativeAtContinuousIndex( const ContinuousIndexType & x ) const;
*/

  void SetSplineOrder(unsigned int SplineOrder);

  itkGetConstMacro(SplineOrder, int);

  void SetInputImage(const TImageType *inputData) override;

protected:
  ComplexBSplineInterpolateImageFunction();
  ~ComplexBSplineInterpolateImageFunction() override {}

  void PrintSelf(std::ostream & os, Indent indent) const override;

private:
  unsigned int m_SplineOrder;

  typename InterpolatorType::Pointer m_RealInterpolator;
  typename InterpolatorType::Pointer m_ImaginaryInterpolator;

  typename RealFilterType::Pointer m_RealFilter;

  typename ImaginaryFilterType::Pointer m_ImaginaryFilter;
}; // class
} // namespace

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

#endif