itkSimplexMesh.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 itkSimplexMesh_h
#define itkSimplexMesh_h

#include "itkMesh.h"
#include "itkSimplexMeshGeometry.h"
#include "itkVertexCell.h"
#include "itkTriangleCell.h"
#include "itkFixedArray.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< SizeValueType > NeighborSetType;

  typedef typename NeighborSetType::iterator NeighborSetIterator;

  typedef std::vector< SizeValueType > 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< SizeValueType, SimplexMeshGeometry * > GeometryMapType;

  typedef typename GeometryMapType::Pointer GeometryMapPointer;

  typedef typename GeometryMapType::Iterator        GeometryMapIterator;
  typedef typename GeometryMapType::ConstIterator   GeometryMapConstIterator;

  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::CellsContainerConstPointer   CellsContainerConstPointer;
  typedef typename Superclass::CellsContainerIterator       CellsContainerIterator;
  typedef typename Superclass::CellsContainerConstIterator  CellsContainerConstIterator;
  typedef typename Superclass::CellIdentifier               CellIdentifier;

  itkSetMacro(GeometryData, GeometryMapPointer);

  itkGetConstReferenceMacro(GeometryData, GeometryMapPointer);

  itkSetMacro(LastCellId, CellIdentifier);

  itkGetConstMacro(LastCellId, CellIdentifier);

  virtual void CopyInformation(const DataObject *data) ITK_OVERRIDE;

  CellIdentifier AddEdge(PointIdentifier startPointId, PointIdentifier endPointId);

  CellIdentifier AddFace(CellAutoPointer & cellPointer);

  CellIdentifier ReplaceFace(CellIdentifier replaceIndex, CellAutoPointer & cellPointer);

  IndexArray GetNeighbors(PointIdentifier pointId) const;

  NeighborListType * GetNeighbors(PointIdentifier pointId, unsigned int radius, NeighborListType *list = ITK_NULLPTR) const;

  void AddNeighbor(PointIdentifier pointId, PointIdentifier neighborId);

  void ReplaceNeighbor(PointIdentifier pointId, PointIdentifier oldNeighborId, PointIdentifier newNeighborIdx);

  void SwapNeighbors(PointIdentifier pointId, PointIdentifier firstNeighborId, PointIdentifier secondNeighborId);

  void SetGeometryData(PointIdentifier pointId, SimplexMeshGeometry *);

  void SetBarycentricCoordinates(PointIdentifier idx, PointType values);

  PointType GetBarycentricCoordinates(PointIdentifier idx) const;

  void SetReferenceMetrics(PointIdentifier idx, PointType values);

  PointType GetReferenceMetrics(PointIdentifier idx) const;

  void SetPhi(PointIdentifier idx, double values);

  double GetPhi(PointIdentifier idx) const;

  void SetMeanCurvature(PointIdentifier idx, double values);

  double GetMeanCurvature(PointIdentifier idx) const;

  void SetRadius(PointIdentifier idx, double values);

  double GetRadius(PointIdentifier idx) const;

  void SetDistance(PointIdentifier idx, double values);

  double GetDistance(PointIdentifier idx) const;

  CovariantVectorType ComputeNormal(PointIdentifier idx) const;

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

  GeometryMapPointer m_GeometryData;

  CellIdentifier 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.hxx"
#endif

#endif