itkConformalFlatteningMeshFilter.h

Go to the documentation of this file.

/*=========================================================================

  Program:   Insight Segmentation & Registration Toolkit
  Module:    $RCSfile: itkConformalFlatteningMeshFilter.h,v $
  Language:  C++
  Date:      $Date: 2008-07-09 20:17:38 $
  Version:   $Revision: 1.12 $

  Copyright (c) Insight Software Consortium. All rights reserved.
  See ITKCopyright.txt or http://www.itk.org/HTML/Copyright.htm for details.

     This software is distributed WITHOUT ANY WARRANTY; without even
     the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR
     PURPOSE.  See the above copyright notices for more information.

=========================================================================*/
#ifndef __itkConformalFlatteningMeshFilter_h
#define __itkConformalFlatteningMeshFilter_h

#include "itkMesh.h"
#include "itkMeshToMeshFilter.h"

// vnl headers
#include <vnl/vnl_math.h>
#include <vnl/vnl_sparse_matrix.h>

// #include "itkConformalFlatteningFunction.h"

namespace itk
{

template <class TInputMesh, class TOutputMesh>
class ITK_EXPORT ConformalFlatteningMeshFilter :
    public MeshToMeshFilter<TInputMesh, TOutputMesh>
{
public:
  typedef ConformalFlatteningMeshFilter  Self;

  typedef TInputMesh   InputMeshType;
  typedef TOutputMesh  OutputMeshType;

  typedef MeshToMeshFilter<TInputMesh,TOutputMesh> Superclass;

  typedef SmartPointer<Self>        Pointer;
  typedef SmartPointer<const Self>  ConstPointer;

  typedef typename InputMeshType::Pointer      InputMeshPointer;
  typedef typename OutputMeshType::Pointer     OutputMeshPointer;
  typedef typename InputMeshType::PointType    InputPointType;
  typedef typename OutputMeshType::PointType   OutputPointType;

  //typedef typename TInputMesh::CoordRepType          CoordRepType;
  typedef double CoordRepType;

  itkNewMacro(Self);

  itkTypeMacro(ConformalFlatteningMeshFilter,MeshToMeshFilter);

  void SetInput(TInputMesh *input);

  itkStaticConstMacro(InputPointDimension, unsigned int,
     ::itk::GetMeshDimension< TInputMesh >::PointDimension );
  itkStaticConstMacro(OutputPointDimension, unsigned int,
     ::itk::GetMeshDimension< TOutputMesh >::PointDimension );

  typedef typename InputMeshType::PointsContainer    PointsContainer;
  typedef typename InputMeshType::CellsContainer     CellsContainer;
  typedef typename InputMeshType::PointIdentifier    PointIdentifier;
  typedef typename InputMeshType::CellIdentifier     CellIdentifier;
  typedef typename PointsContainer::ConstIterator    PointIterator;
  typedef typename CellsContainer::ConstIterator     CellIterator;
  typedef typename InputMeshType::CellType           CellType;
  typedef typename CellType::PointIdIterator         PointIdIterator;
  typedef typename CellType::CellAutoPointer         CellAutoPointer;

  void SetPolarCellIdentifier( CellIdentifier cellId );

  void SetScale( double );

  void MapToSphere( void );

  void MapToPlane( void );

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

  virtual void GenerateData( void );

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

  typedef vnl_vector< CoordRepType >       VectorCoordType;
  typedef vnl_sparse_matrix<CoordRepType>  SparseMatrixCoordType;

  unsigned int m_PolarCellIdentifier;

  bool    m_MapToSphere;

  double  m_MapScale;
};

} // end namespace itk

#ifndef ITK_MANUAL_INSTANTIATION
#include "itkConformalFlatteningMeshFilter.txx"
#endif

#endif