itkComplexBSplineInterpolateImageFunction.h

Go to the documentation of this file.

/*=========================================================================
 *
 *  Copyright NumFOCUS
 *
 *  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_MOVE(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 SizeType = typename Superclass::SizeType;
 
  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;
 
  SizeType
  GetRadius() const override
  {
    return SizeType::Filled(m_SplineOrder + 1);
  }
 
protected:
  ComplexBSplineInterpolateImageFunction();
  ~ComplexBSplineInterpolateImageFunction() override = default;
 
  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 itk
 
#ifndef ITK_MANUAL_INSTANTIATION
#  include "itkComplexBSplineInterpolateImageFunction.hxx"
#endif
 
#endif