/usr/include/InsightToolkit/Review/itkQuadEdgeMeshDiscreteGaussianCurvatureEstimator.h is in libinsighttoolkit3-dev 3.20.1-1.
This file is owned by root:root, with mode 0o644.
The actual contents of the file can be viewed below.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 | /*=========================================================================
Program: Insight Segmentation & Registration Toolkit
Module: itkQuadEdgeMeshDiscreteGaussianCurvatureEstimator.h
Language: C++
Date: $Date$
Version: $Revision$
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 __itkQuadEdgeMeshDiscreteGaussianCurvatureEstimator_h
#define __itkQuadEdgeMeshDiscreteGaussianCurvatureEstimator_h
#include "itkQuadEdgeMeshDiscreteCurvatureEstimator.h"
#include <vnl/vnl_math.h>
//#include "itkTriangleHelper.h"
namespace itk
{
/**
* \class QuadEdgeMeshDiscreteGaussianCurvatureEstimator
* \brief see the following paper
* title: Discrete Differential-Geometry Operators for Triangulated 2-Manifolds
* authors: Mark Meyer, Mathieu Desbrun, Peter Schroder, Alan H. Barr
* conference: VisMath '02
* location: Berlin (Germany)
* \author: Arnaud Gelas, Alexandre Gouaillard
*/
template< class TInputMesh, class TOutputMesh >
class QuadEdgeMeshDiscreteGaussianCurvatureEstimator :
public QuadEdgeMeshDiscreteCurvatureEstimator< TInputMesh, TOutputMesh >
{
public:
typedef QuadEdgeMeshDiscreteGaussianCurvatureEstimator Self;
typedef SmartPointer< Self > Pointer;
typedef SmartPointer< const Self > ConstPointer;
typedef QuadEdgeMeshDiscreteCurvatureEstimator<
TInputMesh, TOutputMesh > Superclass;
typedef typename Superclass::InputMeshType InputMeshType;
typedef typename Superclass::InputMeshPointer InputMeshPointer;
typedef typename Superclass::OutputMeshType OutputMeshType;
typedef typename Superclass::OutputMeshPointer OutputMeshPointer;
typedef typename Superclass::OutputPointsContainerPointer OutputPointsContainerPointer;
typedef typename Superclass::OutputPointsContainerIterator OutputPointsContainerIterator;
typedef typename Superclass::OutputPointType OutputPointType;
typedef typename Superclass::OutputVectorType OutputVectorType;
typedef typename Superclass::OutputCoordType OutputCoordType;
typedef typename Superclass::OutputPointIdentifier OutputPointIdentifier;
typedef typename Superclass::OutputCellIdentifier OutputCellIdentifier;
typedef typename Superclass::OutputQEType OutputQEType;
typedef typename Superclass::OutputMeshTraits OutputMeshTraits;
typedef typename Superclass::OutputCurvatureType OutputCurvatureType;
typedef typename Superclass::TriangleType TriangleType;
/** Run-time type information (and related methods). */
itkTypeMacro( QuadEdgeMeshDiscreteGaussianCurvatureEstimator, QuadEdgeMeshDiscreteCurvatureEstimator );
/** New macro for creation of through a Smart Pointer */
itkNewMacro( Self );
protected:
QuadEdgeMeshDiscreteGaussianCurvatureEstimator() {}
~QuadEdgeMeshDiscreteGaussianCurvatureEstimator() {}
virtual OutputCurvatureType EstimateCurvature( const OutputPointType& iP )
{
OutputMeshPointer output = this->GetOutput();
OutputQEType* qe = iP.GetEdge( );
if( qe != 0 )
{
OutputQEType* qe_it = qe;
OutputQEType* qe_it2;
OutputPointType q0, q1;
OutputCurvatureType sum_theta = 0.;
OutputCurvatureType area = 0.;
do
{
// cell_id = qe_it->GetLeft();
qe_it2 = qe_it->GetOnext();
q0 = output->GetPoint( qe_it->GetDestination() );
q1 = output->GetPoint( qe_it2->GetDestination() );
// Compute Angle;
sum_theta += static_cast< OutputCurvatureType >(
TriangleType::ComputeAngle( q0, iP, q1 ) );
area += ComputeMixedArea( qe_it, qe_it2 );
qe_it = qe_it2;
} while( qe_it != qe );
return ( 2.0 * vnl_math::pi - sum_theta ) / area;
}
return 0.;
}
private:
QuadEdgeMeshDiscreteGaussianCurvatureEstimator( const Self& ); // purposely not implemented
void operator = ( const Self& ); // purposely not implemented
};
}
#endif
|