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
 *
 *         https://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 TCoordinate = double, typename TCoefficientType = double>
class ITK_TEMPLATE_EXPORT ComplexBSplineInterpolateImageFunction
  : public InterpolateImageFunction<TImageType, TCoordinate>
{
public:
  ITK_DISALLOW_COPY_AND_MOVE(ComplexBSplineInterpolateImageFunction);

  using Self = ComplexBSplineInterpolateImageFunction;

  using Superclass = InterpolateImageFunction<TImageType, TCoordinate>;

  using Pointer = SmartPointer<Self>;

  using ConstPointer = SmartPointer<const Self>;

  itkOverrideGetNameOfClassMacro(ComplexBSplineInterpolateImageFunction);

  itkNewMacro(Self);

  static constexpr unsigned int ImageDimension = Superclass::ImageDimension;

  using typename Superclass::OutputType;

  using typename Superclass::InputImageType;

  using typename Superclass::IndexType;

  using typename Superclass::SizeType;

  using typename Superclass::ContinuousIndexType;

  using typename Superclass::PointType;

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

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

  using InterpolatorType = BSplineInterpolateImageFunction<InternalImageType, TCoordinate, 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, unsigned 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