Sophie

Sophie

distrib > Mandriva > 2010.0 > i586 > media > contrib-release > by-pkgid > d5e74628f0e673bb8680aebce32b2c04 > files > 16965

itk-doc-3.12.0-1mdv2010.0.i586.rpm

<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.0 Transitional//EN">
<html><head><meta name="robots" content="noindex">
<meta http-equiv="Content-Type" content="text/html;charset=iso-8859-1">
<title>ITK: itkQuadEdgeMeshDiscreteGaussianCurvatureEstimator.h Source File</title>
<link href="DoxygenStyle.css" rel="stylesheet" type="text/css">
</head><body bgcolor="#ffffff">


<!--  Section customized for INSIGHT : Tue Jul 17 01:02:45 2001 -->
<center>
<a href="index.html" class="qindex">Main Page</a>&nbsp;&nbsp; 
<a href="modules.html" class="qindex">Groups</a>&nbsp;&nbsp;
<a href="namespaces.html" class="qindex">Namespace List</a>&nbsp;&nbsp;
<a href="hierarchy.html" class="qindex">Class Hierarchy</a>&nbsp;&nbsp;
<a href="classes.html" class="qindex">Alphabetical List</a>&nbsp;&nbsp;
<a href="annotated.html" class="qindex">Compound List</a>&nbsp;&nbsp; 
<a href="files.html" class="qindex">File
List</a>&nbsp;&nbsp; 
<a href="namespacemembers.html" class="qindex">Namespace Members</a>&nbsp;&nbsp; 
<a href="functions.html" class="qindex">Compound Members</a>&nbsp;&nbsp; 
<a href="globals.html" class="qindex">File Members</a>&nbsp;&nbsp;
<a href="pages.html" class="qindex">Concepts</a></center>


<!-- Generated by Doxygen 1.5.9 -->
  <div class="navpath"><a class="el" href="dir_02e23d841c198393a6e961e9a0678fff.html">Code</a>&nbsp;&raquo&nbsp;<a class="el" href="dir_dd352b1cbb04f9ff09fa5f33693fed9f.html">Review</a>
  </div>
<div class="contents">
<h1>itkQuadEdgeMeshDiscreteGaussianCurvatureEstimator.h</h1><a href="itkQuadEdgeMeshDiscreteGaussianCurvatureEstimator_8h.html">Go to the documentation of this file.</a><div class="fragment"><pre class="fragment"><a name="l00001"></a>00001 <span class="comment">/*=========================================================================</span>
<a name="l00002"></a>00002 <span class="comment"></span>
<a name="l00003"></a>00003 <span class="comment">  Program:   Insight Segmentation &amp; Registration Toolkit</span>
<a name="l00004"></a>00004 <span class="comment">  Module:    $RCSfile: itkQuadEdgeMeshDiscreteGaussianCurvatureEstimator.h,v $</span>
<a name="l00005"></a>00005 <span class="comment">  Language:  C++</span>
<a name="l00006"></a>00006 <span class="comment">  Date:      $Date: 2008-10-01 21:40:25 $</span>
<a name="l00007"></a>00007 <span class="comment">  Version:   $Revision: 1.3 $</span>
<a name="l00008"></a>00008 <span class="comment"></span>
<a name="l00009"></a>00009 <span class="comment">  Copyright (c) Insight Software Consortium. All rights reserved.</span>
<a name="l00010"></a>00010 <span class="comment">  See ITKCopyright.txt or http://www.itk.org/HTML/Copyright.htm for details.</span>
<a name="l00011"></a>00011 <span class="comment"></span>
<a name="l00012"></a>00012 <span class="comment">     This software is distributed WITHOUT ANY WARRANTY; without even</span>
<a name="l00013"></a>00013 <span class="comment">     the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR</span>
<a name="l00014"></a>00014 <span class="comment">     PURPOSE.  See the above copyright notices for more information.</span>
<a name="l00015"></a>00015 <span class="comment"></span>
<a name="l00016"></a>00016 <span class="comment">=========================================================================*/</span>
<a name="l00017"></a>00017 
<a name="l00018"></a>00018 <span class="preprocessor">#ifndef __itkQuadEdgeMeshDiscreteGaussianCurvatureEstimator_h</span>
<a name="l00019"></a>00019 <span class="preprocessor"></span><span class="preprocessor">#define __itkQuadEdgeMeshDiscreteGaussianCurvatureEstimator_h</span>
<a name="l00020"></a>00020 <span class="preprocessor"></span>
<a name="l00021"></a>00021 <span class="preprocessor">#include "<a class="code" href="itkQuadEdgeMeshDiscreteCurvatureEstimator_8h.html">itkQuadEdgeMeshDiscreteCurvatureEstimator.h</a>"</span>
<a name="l00022"></a>00022 <span class="preprocessor">#include &lt;vnl/vnl_math.h&gt;</span>
<a name="l00023"></a>00023 
<a name="l00024"></a>00024 <span class="comment">//#include "itkTriangleHelper.h"</span>
<a name="l00025"></a>00025 
<a name="l00026"></a>00026 <span class="keyword">namespace </span>itk
<a name="l00027"></a>00027 {
<a name="l00037"></a>00037 <span class="keyword">template</span>&lt; <span class="keyword">class</span> TInputMesh, <span class="keyword">class</span> TOutputMesh &gt;
<a name="l00038"></a><a class="code" href="classitk_1_1QuadEdgeMeshDiscreteGaussianCurvatureEstimator.html">00038</a> <span class="keyword">class </span><a class="code" href="classitk_1_1QuadEdgeMeshDiscreteGaussianCurvatureEstimator.html" title="see the following paper title: Discrete Differential-Geometry Operators for Triangulated...">QuadEdgeMeshDiscreteGaussianCurvatureEstimator</a> :
<a name="l00039"></a>00039   <span class="keyword">public</span> <a class="code" href="classitk_1_1QuadEdgeMeshDiscreteCurvatureEstimator.html" title="FIXME.">QuadEdgeMeshDiscreteCurvatureEstimator</a>&lt; TInputMesh, TOutputMesh &gt;
<a name="l00040"></a>00040 {
<a name="l00041"></a>00041 <span class="keyword">public</span>:
<a name="l00042"></a><a class="code" href="classitk_1_1QuadEdgeMeshDiscreteGaussianCurvatureEstimator.html#f2e185e98615f304485ff52fdc6c9058">00042</a>   <span class="keyword">typedef</span> <a class="code" href="classitk_1_1QuadEdgeMeshDiscreteGaussianCurvatureEstimator.html" title="see the following paper title: Discrete Differential-Geometry Operators for Triangulated...">QuadEdgeMeshDiscreteGaussianCurvatureEstimator</a>        <a class="code" href="classitk_1_1LightObject.html" title="Light weight base class for most itk classes.">Self</a>;
<a name="l00043"></a><a class="code" href="classitk_1_1QuadEdgeMeshDiscreteGaussianCurvatureEstimator.html#5f79efb978d0b87574e61bb3c2816901">00043</a>   <span class="keyword">typedef</span> <a class="code" href="classitk_1_1SmartPointer.html">SmartPointer&lt; Self &gt;</a>                                  <a class="code" href="classitk_1_1SmartPointer.html">Pointer</a>;
<a name="l00044"></a><a class="code" href="classitk_1_1QuadEdgeMeshDiscreteGaussianCurvatureEstimator.html#a305b700503b486be6b3ab188c1115df">00044</a>   <span class="keyword">typedef</span> <a class="code" href="classitk_1_1SmartPointer.html">SmartPointer&lt; const Self &gt;</a>                            <a class="code" href="classitk_1_1SmartPointer.html">ConstPointer</a>;
<a name="l00045"></a>00045   <span class="keyword">typedef</span> <a class="code" href="classitk_1_1QuadEdgeMeshDiscreteCurvatureEstimator.html" title="FIXME.">QuadEdgeMeshDiscreteCurvatureEstimator</a>&lt; 
<a name="l00046"></a><a class="code" href="classitk_1_1QuadEdgeMeshDiscreteGaussianCurvatureEstimator.html#958b46a02a2519b5064e4be460de6d3d">00046</a>     TInputMesh, TOutputMesh &gt;                                   <a class="code" href="classitk_1_1LightObject.html" title="Light weight base class for most itk classes.">Superclass</a>;
<a name="l00047"></a>00047 
<a name="l00048"></a><a class="code" href="classitk_1_1QuadEdgeMeshDiscreteGaussianCurvatureEstimator.html#127c150868f29269a379532f3c80abe5">00048</a>   <span class="keyword">typedef</span> <span class="keyword">typename</span> <a class="code" href="classitk_1_1QuadEdgeMeshDiscreteCurvatureEstimator.html#91442a1ff621ee0cea1e72ad060251e9">Superclass::InputMeshType</a>                    <a class="code" href="classitk_1_1QuadEdgeMeshDiscreteGaussianCurvatureEstimator.html#127c150868f29269a379532f3c80abe5">InputMeshType</a>;
<a name="l00049"></a><a class="code" href="classitk_1_1QuadEdgeMeshDiscreteGaussianCurvatureEstimator.html#836ace6e2e5a4cab80aa1f3643a73e52">00049</a>   <span class="keyword">typedef</span> <span class="keyword">typename</span> <a class="code" href="classitk_1_1QuadEdgeMeshDiscreteCurvatureEstimator.html#ab5f323e579ba0f817e8cc23df22fb66">Superclass::InputMeshPointer</a>                 <a class="code" href="classitk_1_1QuadEdgeMeshDiscreteGaussianCurvatureEstimator.html#836ace6e2e5a4cab80aa1f3643a73e52">InputMeshPointer</a>;
<a name="l00050"></a>00050 
<a name="l00051"></a><a class="code" href="classitk_1_1QuadEdgeMeshDiscreteGaussianCurvatureEstimator.html#5671c72f6201cd24eb3d0f6ed7d6e77f">00051</a>   <span class="keyword">typedef</span> <span class="keyword">typename</span> <a class="code" href="classitk_1_1QuadEdgeMeshDiscreteCurvatureEstimator.html#06e0df9db398b91c45a64783648fcc7c">Superclass::OutputMeshType</a>                   <a class="code" href="classitk_1_1QuadEdgeMeshDiscreteGaussianCurvatureEstimator.html#5671c72f6201cd24eb3d0f6ed7d6e77f">OutputMeshType</a>;
<a name="l00052"></a><a class="code" href="classitk_1_1QuadEdgeMeshDiscreteGaussianCurvatureEstimator.html#1d4acc802a850164f71d5fc48d19234f">00052</a>   <span class="keyword">typedef</span> <span class="keyword">typename</span> <a class="code" href="classitk_1_1QuadEdgeMeshDiscreteCurvatureEstimator.html#039918d090d8f4d684e9f49b040a5a97">Superclass::OutputMeshPointer</a>                <a class="code" href="classitk_1_1QuadEdgeMeshDiscreteGaussianCurvatureEstimator.html#1d4acc802a850164f71d5fc48d19234f">OutputMeshPointer</a>;
<a name="l00053"></a><a class="code" href="classitk_1_1QuadEdgeMeshDiscreteGaussianCurvatureEstimator.html#c045d30330ae477d5b0b3a9c63631766">00053</a>   <span class="keyword">typedef</span> <span class="keyword">typename</span> <a class="code" href="classitk_1_1QuadEdgeMeshDiscreteCurvatureEstimator.html#65226e1cda75c903a516261a4143d7ff">Superclass::OutputPointsContainerPointer</a>     <a class="code" href="classitk_1_1QuadEdgeMeshDiscreteGaussianCurvatureEstimator.html#c045d30330ae477d5b0b3a9c63631766">OutputPointsContainerPointer</a>;
<a name="l00054"></a><a class="code" href="classitk_1_1QuadEdgeMeshDiscreteGaussianCurvatureEstimator.html#7c22d80b6e308619b56e8788d35c5311">00054</a>   <span class="keyword">typedef</span> <span class="keyword">typename</span> <a class="code" href="classitk_1_1QuadEdgeMeshDiscreteCurvatureEstimator.html#723bb495bca3a65338c46b962d3ff9a2">Superclass::OutputPointsContainerIterator</a>    <a class="code" href="classitk_1_1QuadEdgeMeshDiscreteGaussianCurvatureEstimator.html#7c22d80b6e308619b56e8788d35c5311">OutputPointsContainerIterator</a>;
<a name="l00055"></a><a class="code" href="classitk_1_1QuadEdgeMeshDiscreteGaussianCurvatureEstimator.html#b11a090c2d9d95058620c4ed93b53340">00055</a>   <span class="keyword">typedef</span> <span class="keyword">typename</span> <a class="code" href="classitk_1_1QuadEdgeMeshDiscreteCurvatureEstimator.html#fe038b4fc856a0aa78984238bdfce8e3">Superclass::OutputPointType</a>                  <a class="code" href="classitk_1_1QuadEdgeMeshDiscreteGaussianCurvatureEstimator.html#b11a090c2d9d95058620c4ed93b53340">OutputPointType</a>;
<a name="l00056"></a><a class="code" href="classitk_1_1QuadEdgeMeshDiscreteGaussianCurvatureEstimator.html#e89b3d10df30a826c1fc160e8e3b7640">00056</a>   <span class="keyword">typedef</span> <span class="keyword">typename</span> <a class="code" href="classitk_1_1QuadEdgeMeshToQuadEdgeMeshFilter.html#4cd8abe396063fd821a1b8c0f35abf14">Superclass::OutputVectorType</a>                 <a class="code" href="classitk_1_1QuadEdgeMeshDiscreteGaussianCurvatureEstimator.html#e89b3d10df30a826c1fc160e8e3b7640">OutputVectorType</a>;
<a name="l00057"></a><a class="code" href="classitk_1_1QuadEdgeMeshDiscreteGaussianCurvatureEstimator.html#842b1a97df009518ea5e6ad88b59ff99">00057</a>   <span class="keyword">typedef</span> <span class="keyword">typename</span> <a class="code" href="classitk_1_1QuadEdgeMeshDiscreteCurvatureEstimator.html#46eabf0b5dc961091b4e36fb3754de7a">Superclass::OutputCoordType</a>                  <a class="code" href="classitk_1_1QuadEdgeMeshDiscreteGaussianCurvatureEstimator.html#842b1a97df009518ea5e6ad88b59ff99">OutputCoordType</a>;
<a name="l00058"></a><a class="code" href="classitk_1_1QuadEdgeMeshDiscreteGaussianCurvatureEstimator.html#2dbf0b4e18a719b44a343f60c753db2d">00058</a>   <span class="keyword">typedef</span> <span class="keyword">typename</span> <a class="code" href="classitk_1_1QuadEdgeMeshDiscreteCurvatureEstimator.html#13f35915a451eccfdca9960f07f7deb2">Superclass::OutputPointIdentifier</a>            <a class="code" href="classitk_1_1QuadEdgeMeshDiscreteGaussianCurvatureEstimator.html#2dbf0b4e18a719b44a343f60c753db2d">OutputPointIdentifier</a>;
<a name="l00059"></a><a class="code" href="classitk_1_1QuadEdgeMeshDiscreteGaussianCurvatureEstimator.html#00cdba8a4a143c303bfd85a0c6e50bcb">00059</a>   <span class="keyword">typedef</span> <span class="keyword">typename</span> <a class="code" href="classitk_1_1QuadEdgeMeshDiscreteCurvatureEstimator.html#be1ead64ad8b0be5ecc43ee6fc3a0c0d">Superclass::OutputCellIdentifier</a>             <a class="code" href="classitk_1_1QuadEdgeMeshDiscreteGaussianCurvatureEstimator.html#00cdba8a4a143c303bfd85a0c6e50bcb">OutputCellIdentifier</a>;
<a name="l00060"></a><a class="code" href="classitk_1_1QuadEdgeMeshDiscreteGaussianCurvatureEstimator.html#1d801f96ea699bf99cbc92e5131022c9">00060</a>   <span class="keyword">typedef</span> <span class="keyword">typename</span> <a class="code" href="classitk_1_1QuadEdgeMeshDiscreteCurvatureEstimator.html#8c37516d77b4acfad00daac1817e4e96">Superclass::OutputQEType</a>                     <a class="code" href="classitk_1_1QuadEdgeMeshDiscreteGaussianCurvatureEstimator.html#1d801f96ea699bf99cbc92e5131022c9">OutputQEType</a>;
<a name="l00061"></a><a class="code" href="classitk_1_1QuadEdgeMeshDiscreteGaussianCurvatureEstimator.html#3c07284bbd2bd17107e56263f097d67c">00061</a>   <span class="keyword">typedef</span> <span class="keyword">typename</span> <a class="code" href="classitk_1_1QuadEdgeMeshDiscreteCurvatureEstimator.html#d4895f665cd810d9f5dbb8fab76ecf35">Superclass::OutputMeshTraits</a>                 <a class="code" href="classitk_1_1QuadEdgeMeshDiscreteGaussianCurvatureEstimator.html#3c07284bbd2bd17107e56263f097d67c">OutputMeshTraits</a>;
<a name="l00062"></a><a class="code" href="classitk_1_1QuadEdgeMeshDiscreteGaussianCurvatureEstimator.html#1fe6258e373e333d23003685ca6efffa">00062</a>   <span class="keyword">typedef</span> <span class="keyword">typename</span> <a class="code" href="classitk_1_1QuadEdgeMeshDiscreteCurvatureEstimator.html#aee25f84b9a0b54e8a48d9cc1e64213e">Superclass::OutputCurvatureType</a>              <a class="code" href="classitk_1_1QuadEdgeMeshDiscreteGaussianCurvatureEstimator.html#1fe6258e373e333d23003685ca6efffa">OutputCurvatureType</a>;
<a name="l00063"></a><a class="code" href="classitk_1_1QuadEdgeMeshDiscreteGaussianCurvatureEstimator.html#f24e8f0b7f81ba35b3838abd0c6f7214">00063</a>   <span class="keyword">typedef</span> <span class="keyword">typename</span> <a class="code" href="classitk_1_1QuadEdgeMeshDiscreteCurvatureEstimator.html#014ea5654478d537811c5623c153b3ad">Superclass::TriangleType</a>                     <a class="code" href="classitk_1_1TriangleHelper.html" title="Convenient class for various triangles elements computation in 2D or 3D.">TriangleType</a>;
<a name="l00064"></a>00064   
<a name="l00066"></a>00066   <a class="code" href="itkMacro_8h.html#574c1f86661d409dd32c742307a0e96a">itkTypeMacro</a>( <a class="code" href="classitk_1_1QuadEdgeMeshDiscreteGaussianCurvatureEstimator.html" title="see the following paper title: Discrete Differential-Geometry Operators for Triangulated...">QuadEdgeMeshDiscreteGaussianCurvatureEstimator</a>, <a class="code" href="classitk_1_1QuadEdgeMeshDiscreteCurvatureEstimator.html" title="FIXME.">QuadEdgeMeshDiscreteCurvatureEstimator</a> );
<a name="l00067"></a>00067 
<a name="l00069"></a>00069   <a class="code" href="itkMacro_8h.html#0fdc15685afa271857c57e7cd482f30c">itkNewMacro</a>( <a class="code" href="classitk_1_1LightObject.html" title="Light weight base class for most itk classes.">Self</a> );
<a name="l00070"></a>00070 
<a name="l00071"></a>00071 <span class="keyword">protected</span>:
<a name="l00072"></a><a class="code" href="classitk_1_1QuadEdgeMeshDiscreteGaussianCurvatureEstimator.html#09175946f17433a7aaa248563eca97a4">00072</a>   <a class="code" href="classitk_1_1QuadEdgeMeshDiscreteGaussianCurvatureEstimator.html#09175946f17433a7aaa248563eca97a4">QuadEdgeMeshDiscreteGaussianCurvatureEstimator</a>() {}
<a name="l00073"></a><a class="code" href="classitk_1_1QuadEdgeMeshDiscreteGaussianCurvatureEstimator.html#cb8e530212b77b71f83c97f11374a877">00073</a>   <a class="code" href="classitk_1_1QuadEdgeMeshDiscreteGaussianCurvatureEstimator.html#cb8e530212b77b71f83c97f11374a877">~QuadEdgeMeshDiscreteGaussianCurvatureEstimator</a>() {}
<a name="l00074"></a>00074 
<a name="l00075"></a><a class="code" href="classitk_1_1QuadEdgeMeshDiscreteGaussianCurvatureEstimator.html#09169990a8a7551a6dbe45b05d47292e">00075</a>   <span class="keyword">virtual</span> <a class="code" href="classitk_1_1QuadEdgeMeshDiscreteGaussianCurvatureEstimator.html#1fe6258e373e333d23003685ca6efffa">OutputCurvatureType</a> <a class="code" href="classitk_1_1QuadEdgeMeshDiscreteGaussianCurvatureEstimator.html#09169990a8a7551a6dbe45b05d47292e">EstimateCurvature</a>( <span class="keyword">const</span> <a class="code" href="classitk_1_1QuadEdgeMeshDiscreteGaussianCurvatureEstimator.html#b11a090c2d9d95058620c4ed93b53340">OutputPointType</a>&amp; iP )
<a name="l00076"></a>00076     {
<a name="l00077"></a>00077     <a class="code" href="classitk_1_1QuadEdgeMeshDiscreteGaussianCurvatureEstimator.html#1d4acc802a850164f71d5fc48d19234f">OutputMeshPointer</a> output = this-&gt;<a class="code" href="classitk_1_1MeshSource.html#591e055a2e79ba11022b7cbc6e117bd4">GetOutput</a>();
<a name="l00078"></a>00078 
<a name="l00079"></a>00079     <a class="code" href="classitk_1_1QuadEdgeMeshDiscreteGaussianCurvatureEstimator.html#1d801f96ea699bf99cbc92e5131022c9">OutputQEType</a>* qe = iP.GetEdge( );
<a name="l00080"></a>00080 
<a name="l00081"></a>00081     <span class="keywordflow">if</span>( qe != 0 )
<a name="l00082"></a>00082       {
<a name="l00083"></a>00083       <a class="code" href="classitk_1_1QuadEdgeMeshDiscreteGaussianCurvatureEstimator.html#1d801f96ea699bf99cbc92e5131022c9">OutputQEType</a>* qe_it = qe;
<a name="l00084"></a>00084       <a class="code" href="classitk_1_1QuadEdgeMeshDiscreteGaussianCurvatureEstimator.html#1d801f96ea699bf99cbc92e5131022c9">OutputQEType</a>* qe_it2;
<a name="l00085"></a>00085 
<a name="l00086"></a>00086       <a class="code" href="classitk_1_1QuadEdgeMeshDiscreteGaussianCurvatureEstimator.html#b11a090c2d9d95058620c4ed93b53340">OutputPointType</a> q0, q1;
<a name="l00087"></a>00087 
<a name="l00088"></a>00088       <a class="code" href="classitk_1_1QuadEdgeMeshDiscreteGaussianCurvatureEstimator.html#1fe6258e373e333d23003685ca6efffa">OutputCurvatureType</a> sum_theta = 0.;
<a name="l00089"></a>00089       <a class="code" href="classitk_1_1QuadEdgeMeshDiscreteGaussianCurvatureEstimator.html#1fe6258e373e333d23003685ca6efffa">OutputCurvatureType</a> area = 0.;
<a name="l00090"></a>00090 
<a name="l00091"></a>00091       <span class="keywordflow">do</span>
<a name="l00092"></a>00092         {
<a name="l00093"></a>00093         <span class="comment">// cell_id = qe_it-&gt;GetLeft();</span>
<a name="l00094"></a>00094         qe_it2 = qe_it-&gt;GetOnext();
<a name="l00095"></a>00095         q0 = output-&gt;GetPoint( qe_it-&gt;GetDestination() );
<a name="l00096"></a>00096         q1 = output-&gt;GetPoint( qe_it2-&gt;GetDestination() );
<a name="l00097"></a>00097 
<a name="l00098"></a>00098         <span class="comment">// Compute Angle;</span>
<a name="l00099"></a>00099         sum_theta += <span class="keyword">static_cast&lt;</span> <a class="code" href="classitk_1_1QuadEdgeMeshDiscreteGaussianCurvatureEstimator.html#1fe6258e373e333d23003685ca6efffa">OutputCurvatureType</a> <span class="keyword">&gt;</span>(
<a name="l00100"></a>00100               <a class="code" href="classitk_1_1TriangleHelper.html#cbb906e3a4742f84afb5e66e5fc7ed8c">TriangleType::ComputeAngle</a>( q0, iP, q1 ) );
<a name="l00101"></a>00101         area += <a class="code" href="classitk_1_1QuadEdgeMeshDiscreteCurvatureEstimator.html#745ce1ca40f600aa72d6fa926e4a4599">ComputeMixedArea</a>( qe_it, qe_it2 );
<a name="l00102"></a>00102         qe_it = qe_it2;
<a name="l00103"></a>00103         } <span class="keywordflow">while</span>( qe_it != qe );
<a name="l00104"></a>00104 
<a name="l00105"></a>00105       <span class="keywordflow">return</span> ( 2.0 * vnl_math::pi - sum_theta ) / area;
<a name="l00106"></a>00106 
<a name="l00107"></a>00107       }
<a name="l00108"></a>00108 
<a name="l00109"></a>00109     <span class="keywordflow">return</span> 0.;
<a name="l00110"></a>00110   
<a name="l00111"></a>00111     }
<a name="l00112"></a>00112 
<a name="l00113"></a>00113 <span class="keyword">private</span>:
<a name="l00114"></a>00114   <a class="code" href="classitk_1_1QuadEdgeMeshDiscreteGaussianCurvatureEstimator.html#09175946f17433a7aaa248563eca97a4">QuadEdgeMeshDiscreteGaussianCurvatureEstimator</a>( <span class="keyword">const</span> <a class="code" href="classitk_1_1QuadEdgeMeshDiscreteGaussianCurvatureEstimator.html#f2e185e98615f304485ff52fdc6c9058">Self</a>&amp; ); <span class="comment">// purposely not implemented</span>
<a name="l00115"></a>00115   <span class="keywordtype">void</span> operator = ( <span class="keyword">const</span> <a class="code" href="classitk_1_1QuadEdgeMeshDiscreteGaussianCurvatureEstimator.html#f2e185e98615f304485ff52fdc6c9058">Self</a>&amp; ); <span class="comment">// purposely not implemented</span>
<a name="l00116"></a>00116 
<a name="l00117"></a>00117 };
<a name="l00118"></a>00118 
<a name="l00119"></a>00119 }
<a name="l00120"></a>00120 
<a name="l00121"></a>00121 <span class="preprocessor">#endif</span>
<a name="l00122"></a>00122 <span class="preprocessor"></span>
</pre></div></div>
<hr><address><small>
Generated at Thu May 7 23:55:02 2009 for ITK  by <a href="http://www.stack.nl/~dimitri/doxygen/index.html"> <img 
src="http://www.stack.nl/~dimitri/doxygen/doxygen.png" alt="doxygen"
align="middle" border=0 width=110 height=53>
</a> 1.5.9 written by <a href="mailto:dimitri@stack.nl">Dimitri van Heesch</a>,
 &copy;&nbsp;1997-2000</small></address>
</body>
</html>