itkManifoldParzenWindowsPointSetFunction.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 __itkManifoldParzenWindowsPointSetFunction_h
#define __itkManifoldParzenWindowsPointSetFunction_h

#include "itkPointSetFunction.h"

#include "itkGaussianMembershipFunction.h"
#include "itkMatrix.h"
#include "itkPointSet.h"
#include "itkPointsLocator.h"
#include "itkVector.h"

#include <vector>

namespace itk
{

template <typename TPointSet, typename TOutput = double, typename TCoordRep = double>
class ManifoldParzenWindowsPointSetFunction
: public PointSetFunction<TPointSet, TOutput, TCoordRep>
{
public:
  typedef ManifoldParzenWindowsPointSetFunction            Self;
  typedef PointSetFunction<TPointSet, TOutput, TCoordRep>  Superclass;
  typedef SmartPointer<Self>                               Pointer;
  typedef SmartPointer<const Self>                         ConstPointer;

  itkNewMacro(Self);

  itkStaticConstMacro( PointDimension, unsigned int, TPointSet::PointDimension );

  typedef typename Superclass::InputPointSetType   InputPointSetType;
  typedef typename Superclass::InputPointType      InputPointType;

  typedef TPointSet                                        PointSetType;
  typedef typename PointSetType::PointType                 PointType;
  typedef typename PointSetType::PointsContainer           PointsContainer;
  typedef typename PointsContainer::ElementIdentifier      PointIdentifier;

  typedef TOutput                                          RealType;
  typedef TOutput                                          OutputType;
  typedef TCoordRep                                        CoordRepType;

  typedef PointsLocator< PointsContainer>                       PointsLocatorType;
  typedef typename PointsLocatorType::NeighborsIdentifierType   NeighborsIdentifierType;

  typedef typename Statistics::
    GaussianMembershipFunction<PointType>                  GaussianType;
  typedef typename GaussianType::Pointer                   GaussianPointer;
  typedef typename GaussianType::ConstPointer              GaussianConstPointer;
  typedef std::vector<GaussianPointer>                     GaussianContainerType;
  typedef typename GaussianType::CovarianceMatrixType      CovarianceMatrixType;

  itkSetMacro( CovarianceKNeighborhood, unsigned int );

  itkGetConstMacro( CovarianceKNeighborhood, unsigned int );

  itkSetMacro( EvaluationKNeighborhood, unsigned int );

  itkGetConstMacro( EvaluationKNeighborhood, unsigned int );

  itkSetMacro( RegularizationSigma, RealType );

  itkGetConstMacro( RegularizationSigma, RealType );

  itkSetMacro( KernelSigma, RealType );

  itkGetConstMacro( KernelSigma, RealType );

  itkSetMacro( Normalize, bool );

  itkGetConstMacro( Normalize, bool );

  itkBooleanMacro( Normalize );

  itkSetMacro( UseAnisotropicCovariances, bool );

  itkGetConstMacro( UseAnisotropicCovariances, bool );

  itkBooleanMacro( UseAnisotropicCovariances );

  virtual void SetInputPointSet( const InputPointSetType * );

  virtual TOutput Evaluate( const InputPointType & ) const;

  GaussianConstPointer GetGaussian( PointIdentifier ) const;

  itkGetModifiableObjectMacro(PointsLocator, PointsLocatorType );

protected:
  ManifoldParzenWindowsPointSetFunction();
  virtual ~ManifoldParzenWindowsPointSetFunction();
  void PrintSelf( std::ostream& os, Indent indent ) const;

  void GenerateData();

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

  typename PointsLocatorType::Pointer           m_PointsLocator;

  unsigned int                                  m_CovarianceKNeighborhood;
  unsigned int                                  m_EvaluationKNeighborhood;
  RealType                                      m_RegularizationSigma;
  RealType                                      m_KernelSigma;

  GaussianContainerType                         m_Gaussians;
  bool                                          m_Normalize;
  bool                                          m_UseAnisotropicCovariances;
};

} // end namespace itk

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

#endif