ITK  5.2.0
Insight Toolkit
itkDiscretePrincipalCurvaturesQuadEdgeMeshFilter.h
Go to the documentation of this file.
1 /*=========================================================================
2  *
3  * Copyright NumFOCUS
4  *
5  * Licensed under the Apache License, Version 2.0 (the "License");
6  * you may not use this file except in compliance with the License.
7  * You may obtain a copy of the License at
8  *
9  * http://www.apache.org/licenses/LICENSE-2.0.txt
10  *
11  * Unless required by applicable law or agreed to in writing, software
12  * distributed under the License is distributed on an "AS IS" BASIS,
13  * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
14  * See the License for the specific language governing permissions and
15  * limitations under the License.
16  *
17  *=========================================================================*/
18 #ifndef itkDiscretePrincipalCurvaturesQuadEdgeMeshFilter_h
19 #define itkDiscretePrincipalCurvaturesQuadEdgeMeshFilter_h
20 
23 
24 namespace itk
25 {
33 template <typename TInputMesh, typename TOutputMesh = TInputMesh>
35  : public DiscreteCurvatureQuadEdgeMeshFilter<TInputMesh, TOutputMesh>
36 {
37 public:
38  ITK_DISALLOW_COPY_AND_MOVE(DiscretePrincipalCurvaturesQuadEdgeMeshFilter);
39 
44 
59 
61 
64 
66 
67 #ifdef ITK_USE_CONCEPT_CHECKING
68  // Begin concept checking
70  // End concept checking
71 #endif
72 
73 protected:
75  : m_Gaussian(0.0)
76  , m_Mean(0.0)
77  {}
79 
82 
83  void
85  {
86  OutputMeshPointer output = this->GetOutput();
87 
88  OutputQEType * qe = iP.GetEdge();
89 
90  m_Mean = 0.;
91  m_Gaussian = 0.;
92 
93  if (qe != nullptr)
94  {
95  OutputVectorType Laplace;
96  Laplace.Fill(0.);
97 
98  OutputQEType * qe_it = qe;
99 
100  OutputCurvatureType area(0.), sum_theta(0.);
101 
102  if (qe_it != qe_it->GetOnext())
103  {
104  qe_it = qe;
105  OutputQEType * qe_it2;
106 
107  OutputPointType q0, q1;
108  OutputVectorType face_normal;
109 
110  OutputVectorType normal;
111  normal.Fill(0.);
112 
113  OutputCurvatureType temp_area;
114  OutputCoordType temp_coeff;
115 
116  CoefficientType coefficent;
117 
118  do
119  {
120  qe_it2 = qe_it->GetOnext();
121  q0 = output->GetPoint(qe_it->GetDestination());
122  q1 = output->GetPoint(qe_it2->GetDestination());
123 
124  temp_coeff = coefficent(output, qe_it);
125  Laplace += temp_coeff * (iP - q0);
126 
127  // Compute Angle;
128  sum_theta += static_cast<OutputCurvatureType>(TriangleType::ComputeAngle(q0, iP, q1));
129 
130  temp_area = this->ComputeMixedArea(qe_it, qe_it2);
131  area += temp_area;
132 
133  face_normal = TriangleType::ComputeNormal(q0, iP, q1);
134  normal += face_normal;
135 
136  qe_it = qe_it2;
137  } while (qe_it != qe);
138 
139  if (area > 1e-10)
140  {
141  area = 1. / area;
142  Laplace *= 0.25 * area;
143  m_Mean = Laplace * normal;
144  m_Gaussian = (2. * itk::Math::pi - sum_theta) * area;
145  }
146  }
147  }
148  }
149 
150  virtual OutputCurvatureType
152  {
153  return std::max(static_cast<OutputCurvatureType>(0.), m_Mean * m_Mean - m_Gaussian);
154  }
155 };
156 } // namespace itk
157 
158 #endif
itk::DiscreteCurvatureQuadEdgeMeshFilter::ComputeMixedArea
OutputCurvatureType ComputeMixedArea(OutputQEType *iQE1, OutputQEType *iQE2)
Definition: itkDiscreteCurvatureQuadEdgeMeshFilter.h:81
itk::DiscreteCurvatureQuadEdgeMeshFilter::OutputMeshPointer
typename OutputMeshType::Pointer OutputMeshPointer
Definition: itkDiscreteCurvatureQuadEdgeMeshFilter.h:49
itk::DiscretePrincipalCurvaturesQuadEdgeMeshFilter
FIXME add documentation here.
Definition: itkDiscretePrincipalCurvaturesQuadEdgeMeshFilter.h:34
itk::DiscreteCurvatureQuadEdgeMeshFilter::InputMeshType
TInputMesh InputMeshType
Definition: itkDiscreteCurvatureQuadEdgeMeshFilter.h:45
itk::MeshSource::GetOutput
OutputMeshType * GetOutput()
itk::MeshSource::OutputMeshPointer
typename OutputMeshType::Pointer OutputMeshPointer
Definition: itkMeshSource.h:69
itk::DiscreteCurvatureQuadEdgeMeshFilter::OutputMeshType
TOutputMesh OutputMeshType
Definition: itkDiscreteCurvatureQuadEdgeMeshFilter.h:48
itk::QuadEdgeMeshToQuadEdgeMeshFilter::OutputPointsContainerPointer
typename OutputMeshType::PointsContainerPointer OutputPointsContainerPointer
Definition: itkQuadEdgeMeshToQuadEdgeMeshFilter.h:86
itk::MeshToMeshFilter::InputMeshPointer
typename InputMeshType::Pointer InputMeshPointer
Definition: itkMeshToMeshFilter.h:66
itk::DiscreteCurvatureQuadEdgeMeshFilter::InputMeshPointer
typename InputMeshType::Pointer InputMeshPointer
Definition: itkDiscreteCurvatureQuadEdgeMeshFilter.h:46
itkDiscreteCurvatureQuadEdgeMeshFilter.h
itk::SmartPointer< Self >
itk::DiscreteCurvatureQuadEdgeMeshFilter::OutputPointIdentifier
typename OutputMeshType::PointIdentifier OutputPointIdentifier
Definition: itkDiscreteCurvatureQuadEdgeMeshFilter.h:54
itk::DiscretePrincipalCurvaturesQuadEdgeMeshFilter::OutputCurvatureType
typename Superclass::OutputCurvatureType OutputCurvatureType
Definition: itkDiscretePrincipalCurvaturesQuadEdgeMeshFilter.h:58
itk::TriangleHelper
A convenience class for computation of various triangle elements in 2D or 3D.
Definition: itkTriangleHelper.h:31
itk::MeshToMeshFilter::InputMeshType
TInputMesh InputMeshType
Definition: itkMeshToMeshFilter.h:65
itk::DiscreteCurvatureQuadEdgeMeshFilter::TriangleType
TriangleHelper< OutputPointType > TriangleType
Definition: itkDiscreteCurvatureQuadEdgeMeshFilter.h:60
itk::DiscretePrincipalCurvaturesQuadEdgeMeshFilter::m_Gaussian
OutputCurvatureType m_Gaussian
Definition: itkDiscretePrincipalCurvaturesQuadEdgeMeshFilter.h:80
itk::LightObject
Light weight base class for most itk classes.
Definition: itkLightObject.h:59
itk::DiscreteCurvatureQuadEdgeMeshFilter::OutputPointsContainerPointer
typename OutputMeshType::PointsContainerPointer OutputPointsContainerPointer
Definition: itkDiscreteCurvatureQuadEdgeMeshFilter.h:50
itk::ConformalMatrixCoefficients
Compute a matrix filed by Conformal Coefficients of the edge wherever two vertices are connected by a...
Definition: itkQuadEdgeMeshParamMatrixCoefficients.h:127
itk::DiscreteCurvatureQuadEdgeMeshFilter::OutputCurvatureType
typename OutputMeshTraits::PixelType OutputCurvatureType
Definition: itkDiscreteCurvatureQuadEdgeMeshFilter.h:58
itk::QuadEdgeMeshToQuadEdgeMeshFilter::OutputVectorType
typename OutputMeshType::VectorType OutputVectorType
Definition: itkQuadEdgeMeshToQuadEdgeMeshFilter.h:83
itkQuadEdgeMeshParamMatrixCoefficients.h
itk::Concept::IsFloatingPoint
Definition: itkConceptChecking.h:946
itk::DiscreteCurvatureQuadEdgeMeshFilter::OutputMeshTraits
typename OutputMeshType::MeshTraits OutputMeshTraits
Definition: itkDiscreteCurvatureQuadEdgeMeshFilter.h:57
itk::DiscreteCurvatureQuadEdgeMeshFilter::OutputPointType
typename OutputMeshType::PointType OutputPointType
Definition: itkDiscreteCurvatureQuadEdgeMeshFilter.h:52
itk::DiscreteCurvatureQuadEdgeMeshFilter::OutputCoordType
typename OutputPointType::CoordRepType OutputCoordType
Definition: itkDiscreteCurvatureQuadEdgeMeshFilter.h:53
itk::DiscreteCurvatureQuadEdgeMeshFilter
FIXME.
Definition: itkDiscreteCurvatureQuadEdgeMeshFilter.h:35
itkConceptMacro
#define itkConceptMacro(name, concept)
Definition: itkConceptChecking.h:65
itk::DiscretePrincipalCurvaturesQuadEdgeMeshFilter::~DiscretePrincipalCurvaturesQuadEdgeMeshFilter
~DiscretePrincipalCurvaturesQuadEdgeMeshFilter() override=default
itk
The "itk" namespace contains all Insight Segmentation and Registration Toolkit (ITK) classes....
Definition: itkAnnulusOperator.h:24
itk::DiscreteCurvatureQuadEdgeMeshFilter::OutputPointsContainerIterator
typename OutputMeshType::PointsContainerIterator OutputPointsContainerIterator
Definition: itkDiscreteCurvatureQuadEdgeMeshFilter.h:51
itk::DiscretePrincipalCurvaturesQuadEdgeMeshFilter::ComputeDelta
virtual OutputCurvatureType ComputeDelta()
Definition: itkDiscretePrincipalCurvaturesQuadEdgeMeshFilter.h:151
itk::DiscreteCurvatureQuadEdgeMeshFilter::OutputQEType
typename OutputMeshType::QEType OutputQEType
Definition: itkDiscreteCurvatureQuadEdgeMeshFilter.h:56
itk::Math::e
static constexpr double e
Definition: itkMath.h:54
itk::MeshSource::OutputMeshType
TOutputMesh OutputMeshType
Definition: itkMeshSource.h:68
itk::DiscretePrincipalCurvaturesQuadEdgeMeshFilter::DiscretePrincipalCurvaturesQuadEdgeMeshFilter
DiscretePrincipalCurvaturesQuadEdgeMeshFilter()
Definition: itkDiscretePrincipalCurvaturesQuadEdgeMeshFilter.h:74
itk::DiscreteCurvatureQuadEdgeMeshFilter::OutputCellIdentifier
typename OutputMeshType::CellIdentifier OutputCellIdentifier
Definition: itkDiscreteCurvatureQuadEdgeMeshFilter.h:55
itk::Math::pi
static constexpr double pi
Definition: itkMath.h:64
itk::DiscretePrincipalCurvaturesQuadEdgeMeshFilter::ComputeMeanAndGaussianCurvatures
void ComputeMeanAndGaussianCurvatures(const OutputPointType &iP)
Definition: itkDiscretePrincipalCurvaturesQuadEdgeMeshFilter.h:84
itk::QuadEdgeMeshToQuadEdgeMeshFilter::OutputPointIdentifier
typename OutputMeshType::PointIdentifier OutputPointIdentifier
Definition: itkQuadEdgeMeshToQuadEdgeMeshFilter.h:81
itk::QuadEdgeMeshToQuadEdgeMeshFilter::OutputPointType
typename OutputMeshType::PointType OutputPointType
Definition: itkQuadEdgeMeshToQuadEdgeMeshFilter.h:80
itk::DiscretePrincipalCurvaturesQuadEdgeMeshFilter::m_Mean
OutputCurvatureType m_Mean
Definition: itkDiscretePrincipalCurvaturesQuadEdgeMeshFilter.h:81
itk::QuadEdgeMeshToQuadEdgeMeshFilter::OutputPointsContainerIterator
typename OutputMeshType::PointsContainerIterator OutputPointsContainerIterator
Definition: itkQuadEdgeMeshToQuadEdgeMeshFilter.h:85