itkSimplexMesh.h

Go to the documentation of this file.

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

Program:   Insight Segmentation & Registration Toolkit
Module:    $RCSfile: itkSimplexMesh.h,v $
Language:  C++
Date:      $Date: 2007-02-01 18:08:15 $
Version:   $Revision: 1.11 $

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 __itkSimplexMesh_h
#define __itkSimplexMesh_h

#if defined(_MSC_VER)
#pragma warning ( disable : 4786 )
#endif

#include "itkMesh.h"
#include "itkSimplexMeshGeometry.h"
#include "itkVertexCell.h"
#include "itkTriangleCell.h"
#include "itkCellInterface.h"
#include "itkMapContainer.h"
#include "itkFixedArray.h"
#include "itkNumericTraits.h"
#include <itkCovariantVector.h>
#include <vector>
#include <algorithm>
#include <set>

namespace itk
  {
  template <typename TPixelType, unsigned int VDimension = 3,
    typename TMeshTraits = DefaultStaticMeshTraits< TPixelType , VDimension, VDimension, TPixelType , TPixelType , TPixelType >
  >
class SimplexMesh : public Mesh<TPixelType, VDimension, TMeshTraits>
  {
  public:
    typedef SimplexMesh                Self;

    typedef Mesh<TPixelType, VDimension, TMeshTraits>  Superclass;

    typedef SmartPointer<Self>  Pointer;

    typedef SmartPointer<const Self>  ConstPointer;

    typedef typename SimplexMeshGeometry::IndexArray            IndexArray;

    typedef std::set<unsigned long>                             NeighborSetType;

    typedef typename NeighborSetType::iterator                  NeighborSetIterator;

    typedef std::vector<unsigned long>                          NeighborListType;

    typedef typename TMeshTraits::PointType                     PointType;

    typedef typename TMeshTraits::PointIdentifier               PointIdentifier;

    typedef typename PointType::VectorType VectorType;

    typedef CovariantVector<typename VectorType::ValueType, 3>  CovariantVectorType;


    typedef typename Superclass::CellType                       CellType;

    typedef typename CellType::CellAutoPointer                  CellAutoPointer;

    typedef itk::LineCell<CellType>                             LineType;

    typedef itk::MapContainer<unsigned long, SimplexMeshGeometry *>   GeometryMapType;

    typedef typename GeometryMapType::Pointer GeometryMapPointer;

    typedef typename GeometryMapType::Iterator GeometryMapIterator;


    itkNewMacro(Self);

    itkTypeMacro(SimplexMesh, Mesh);

    typedef TMeshTraits                                          MeshTraits;
    typedef typename MeshTraits::PixelType                       PixelType;  
    typedef typename MeshTraits::PointsContainer                 PointsContainer;
    typedef typename Superclass::PointsContainerPointer          PointsContainerPointer;
    typedef typename Superclass::PointsContainer::Iterator       PointsContainerIterator;
    typedef typename Superclass::PointsContainerConstIterator    PointsContainerConstIterator;
    typedef typename Superclass::CellsContainerPointer           CellsContainerPointer;
    typedef typename Superclass::CellsContainerIterator          CellsContainerIterator;

    itkSetMacro(GeometryData, GeometryMapPointer );

    itkGetConstReferenceMacro(GeometryData, GeometryMapPointer );

    itkSetMacro(LastCellId, unsigned long);

    itkGetMacro(LastCellId, unsigned long);

    virtual void CopyInformation(const DataObject *data);


    unsigned long AddEdge(unsigned long startPointId, unsigned long endPointId);


    unsigned long AddFace(CellAutoPointer &cellPointer);

    unsigned long ReplaceFace(unsigned long replaceIndex, CellAutoPointer &cellPointer);

    IndexArray GetNeighbors(unsigned long pointId) const;

    NeighborListType* GetNeighbors(unsigned long pointId, unsigned int radius, NeighborListType* list = NULL) const;

    void AddNeighbor(unsigned long pointId, unsigned long neighborId);

    void ReplaceNeighbor(unsigned long pointId, unsigned long oldNeighborId,unsigned long newNeighborIdx);

    void SwapNeighbors(unsigned long pointId, unsigned long firstNeighborId,unsigned long secondNeighborId);

    void SetGeometryData(unsigned long pointId, SimplexMeshGeometry*);

    void SetBarycentricCoordinates(unsigned long idx, PointType values);

    PointType GetBarycentricCoordinates(unsigned long idx) const;

    void SetReferenceMetrics(unsigned long idx, PointType values);

    PointType GetReferenceMetrics(unsigned long idx) const;

    void SetPhi(unsigned long idx, double values);

    double GetPhi(unsigned long idx) const;

    void SetMeanCurvature(unsigned long idx, double values);

    double GetMeanCurvature(unsigned long idx) const;

    void SetRadius(unsigned long idx, double values);

    double GetRadius(unsigned long idx) const;

    void SetDistance(unsigned long idx, double values);

    double GetDistance(unsigned long idx) const;

    CovariantVectorType ComputeNormal(unsigned long idx ) const;

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

    GeometryMapPointer m_GeometryData;

    /*  */
    unsigned long m_LastCellId;


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



  }; // End Class:  SimplexMesh

  } // end namespace itk

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

#endif