itkPSMRBFCorrespondenceInterpolator.h

/*=========================================================================
 *
 *  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 __itkPSMRBFCorrespondenceInterpolator_h
#define __itkPSMRBFCorrespondenceInterpolator_h

#include <vector>
#include <vnl/vnl_vector.h>
#include <vnl/vnl_matrix.h>
#include <vnl/algo/vnl_svd.h>

#include "itkFunctionBase.h"
#include "itkPoint.h"
#include "itkMacro.h"

namespace itk
{

template <unsigned int VDimension>
class ITK_EXPORT PSMRBFCorrespondenceInterpolator :
  public FunctionBase<Point<double,VDimension>, Point<double,VDimension> >
{
public:
  itkStaticConstMacro(Dimension, unsigned int, VDimension);

  typedef PSMRBFCorrespondenceInterpolator    Self;
  typedef Point<double,VDimension>            PointType;
  typedef FunctionBase <PointType, PointType> Superclass;
  typedef SmartPointer<Self>                  Pointer;
  typedef SmartPointer<const Self>            ConstPointer;

  itkNewMacro(Self);
  
  itkTypeMacro( PSMRBFCorrespondenceInterpolator, FunctionBase );

  void SetPointSetA(const std::vector<PointType> &v)
    {      
      m_PointSetA = v;
    }
  const std::vector<PointType> & GetPointSetA() const
  {
    return m_PointSetA;
  }

  void SetPointSetB(const std::vector<PointType> &v)
    {      
      m_PointSetB = v;
    }
  const std::vector<PointType> & GetPointSetB() const
  {
    return m_PointSetB;
  }

  void Initialize();

  itkGetMacro(Initialized, bool);

  virtual PointType Evaluate( const PointType& pt ) const;

protected:
  PSMRBFCorrespondenceInterpolator();
  virtual ~PSMRBFCorrespondenceInterpolator() {}
  void PrintSelf(std::ostream& os, Indent indent) const
  {
    Superclass::PrintSelf(os,indent);
  }

private:
  PSMRBFCorrespondenceInterpolator(const Self&); //purposely not implemented
  void operator=(const Self&); //purposely not implemented

  bool m_Initialized;

  std::vector<PointType> m_PointSetA;
  std::vector<PointType> m_PointSetB;

  vnl_matrix<double> m_P;
  vnl_matrix<double> m_C;
  
};

} // end namespace itk

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

#endif