text: /*=========================================================================
   *
   *  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 itkBorderQuadEdgeMeshFilter_h
  #define itkBorderQuadEdgeMeshFilter_h
  
  #include "itkAutoPointer.h"
  #include "itkQuadEdgeMesh.h"
  #include "itkQuadEdgeMeshToQuadEdgeMeshFilter.h"
  #include "itkQuadEdgeMeshBoundaryEdgesMeshFunction.h"
  
  namespace itk
  {
  /**
   * \class BorderQuadEdgeMeshFilter
   * \brief Transform one border of a QuadEdgeMesh into either a circle
   * (conformal) or a square (arclength-wise).
   *
   * This class is one important step when computing a planar parameterization
   * of one mesh.
   *
   * If the input mesh has several boundaries, one can choose
   * the one which would be transformed via the variable m_BorderPick.
   *
   * \li <tt>m_BorderPick == Self::LONGEST</tt> refers to the boundary
   * \f$ b \f$ which satisfies:
   * \f[ b = \arg \max_{b^k} \sum_{i=1}^{N^k} \left\| x_{i}^k - x_{i+1}^k \right\| \f]
   *
   * \li <tt>m_BorderPick == Self::LARGEST</tt> refers to the boundary
   * \f$ b \f$ which satisfies:
   * \f[ b = \arg \max_{b^k} N^k \f]
   *
   * \sa ParameterizationQuadEdgeMeshFilter
   * \ingroup ITKQuadEdgeMeshFiltering
   */
  template< typename TInputMesh, typename TOutputMesh=TInputMesh >
  class ITK_TEMPLATE_EXPORT BorderQuadEdgeMeshFilter:
    public QuadEdgeMeshToQuadEdgeMeshFilter< TInputMesh, TOutputMesh >
  {
  public:
    /** Basic types. */
    typedef BorderQuadEdgeMeshFilter    Self;
    typedef QuadEdgeMeshToQuadEdgeMeshFilter< TInputMesh,
                                              TOutputMesh >
                                        Superclass;
    typedef SmartPointer< Self >        Pointer;
    typedef SmartPointer< const Self >  ConstPointer;
  
    typedef TInputMesh                                  InputMeshType;
    typedef typename InputMeshType::ConstPointer        InputMeshConstPointer;
    typedef typename InputMeshType::CoordRepType        InputCoordRepType;
    typedef typename InputMeshType::PointType           InputPointType;
    typedef typename InputMeshType::Traits              InputTraits;
    typedef typename InputMeshType::PointIdentifier     InputPointIdentifier;
    typedef typename InputMeshType::QEType              InputQEType;
    typedef typename InputQEType::IteratorGeom          InputIteratorGeom;
    typedef typename InputMeshType::VectorType          InputVectorType;
    typedef typename InputMeshType::EdgeListType        InputEdgeListType;
    typedef AutoPointer< InputEdgeListType >            InputEdgeListPointerType;
    typedef typename InputEdgeListType::iterator        InputEdgeListIterator;
    typedef typename InputMeshType::EdgeCellType        InputEdgeCellType;
    typedef typename InputMeshType::PolygonCellType     InputPolygonCellType;
    typedef typename InputMeshType::PointIdList         InputPointIdList;
    typedef typename InputMeshType::PointsContainer     InputPointsContainer;
    typedef typename InputMeshType::PointsContainerConstIterator
    InputPointsContainerConstIterator;
    typedef typename InputMeshType::CellsContainerConstIterator
    InputCellsContainerConstIterator;
  
    typedef TOutputMesh                              OutputMeshType;
    typedef typename OutputMeshType::Pointer         OutputMeshPointer;
    typedef typename OutputMeshType::CoordRepType    OutputCoordRepType;
    typedef typename OutputMeshType::PointType       OutputPointType;
    typedef typename OutputMeshType::Traits          OutputTraits;
    typedef typename OutputMeshType::PointIdentifier OutputPointIdentifier;
    typedef typename OutputMeshType::QEType          OutputQEType;
    typedef typename OutputMeshType::VectorType      OutputVectorType;
    typedef typename OutputMeshType::EdgeListType    OutputEdgeListType;
    typedef typename OutputMeshType::EdgeCellType    OutputEdgeCellType;
    typedef typename OutputMeshType::PolygonCellType OutputPolygonCellType;
    typedef typename OutputMeshType::PointIdList     OutputPointIdList;
    typedef typename OutputMeshType::PointsContainer OutputPointsContainer;
    typedef typename OutputMeshType::PointsContainerConstIterator
    OutputPointsContainerConstIterator;
    typedef typename OutputMeshType::CellsContainerConstIterator
    OutputCellsContainerConstIterator;
  
    itkNewMacro(Self);
    itkTypeMacro(BorderQuadEdgeMeshFilter, QuadEdgeMeshToQuadEdgeMeshFilter);
    itkStaticConstMacro(PointDimension, unsigned int,
                        InputTraits::PointDimension);
  
    typedef std::vector< InputPointType >                           InputVectorPointType;
    typedef std::map< InputPointIdentifier, OutputPointIdentifier > MapPointIdentifier;
    typedef typename MapPointIdentifier::iterator                   MapPointIdentifierIterator;
  
    typedef QuadEdgeMeshBoundaryEdgesMeshFunction< InputMeshType > BoundaryRepresentativeEdgesType;
    typedef typename BoundaryRepresentativeEdgesType::Pointer      BoundaryRepresentativeEdgesPointer;
  
    enum BorderTransformType {
      SQUARE_BORDER_TRANSFORM = 0,
      DISK_BORDER_TRANSFORM
      };
  
    enum BorderPickType {
      LONGEST = 0,
      LARGEST
      };
  
    itkSetMacro(TransformType, BorderTransformType);
    itkGetConstMacro(TransformType, BorderTransformType);
  
    itkSetMacro( BorderPick, BorderPickType );
    itkGetConstMacro( BorderPick, BorderPickType );
  
    itkSetMacro(Radius, InputCoordRepType);
    itkGetConstMacro(Radius, InputCoordRepType);
  
    void ComputeTransform();
  
    MapPointIdentifier GetBoundaryPtMap();
  
    InputVectorPointType GetBorder();
  
  protected:
    BorderQuadEdgeMeshFilter();
  
    virtual ~BorderQuadEdgeMeshFilter() ITK_OVERRIDE {}
  
    void PrintSelf(std::ostream & os, Indent indent) const ITK_OVERRIDE;
  
    BorderTransformType m_TransformType;
    BorderPickType      m_BorderPick;
  
    InputCoordRepType m_Radius;
  
    InputVectorPointType m_Border;
  
    MapPointIdentifier m_BoundaryPtMap;
  
    void GenerateData() ITK_OVERRIDE;
  
    void ComputeBoundary();
  
    InputQEType* ComputeLongestBorder();
  
    InputQEType* ComputeLargestBorder();
  
    void DiskTransform();
  
    InputPointType GetMeshBarycentre();
  
    InputCoordRepType RadiusMaxSquare();
  
    void ArcLengthSquareTransform();
  
  private:
    ITK_DISALLOW_COPY_AND_ASSIGN(BorderQuadEdgeMeshFilter);
  };
  } // end namespace itk
  
  #ifndef ITK_MANUAL_INSTANTIATION
  #include "itkBorderQuadEdgeMeshFilter.hxx"
  #endif
  
  #endif
  