ITK/Examples/Meshes/AreaAndVolume: Difference between revisions

From KitwarePublic
< ITK‎ | Examples
Jump to navigationJump to search
No edit summary
(Deprecated content that is moved to sphinx)
 
(2 intermediate revisions by 2 users not shown)
Line 1: Line 1:
This example demonstrates how to calculate the area and volume of an itk::Mesh by (a) ensuring that it is triangular and (b) converting it to an itk::SimplexMesh.
{{warning|1=The media wiki content on this page is no longer maintainedThe examples presented on the https://itk.org/Wiki/*  pages likely require ITK version 4.13 or earlier releases.   In many cases, the examples on this page no longer conform to the best practices for modern ITK versions.}}
 
Contributed by: Davis Vigneault
 
==AreaAndVolume.cxx==
 
<source lang="cpp">
 
#include <itkSimplexMesh.h>
#include <itkRegularSphereMeshSource.h>
#include <itkTriangleMeshToSimplexMeshFilter.h>
#include <itkSimplexMeshVolumeCalculator.h>
 
typedef itk::Mesh< float, 3 >        TMesh;
typedef itk::SimplexMesh< float, 3 > TSimplex;
typedef itk::RegularSphereMeshSource< TMesh >                  TSphere;
typedef itk::TriangleMeshToSimplexMeshFilter< TMesh, TSimplex > TConvert;
typedef itk::SimplexMeshVolumeCalculator< TSimplex >            TVolume;
 
int main()
{
 
  // Create a spherical mesh with known radius and resolution.
  TSphere::Pointer source = TSphere::New();
  TSphere::VectorType scale;
  scale.Fill( 5.0 );
  source->SetScale( scale );
  source->SetResolution( 5 );
  source->Update();
 
  // Ensure that all cells of the mesh are triangles.
  for (TMesh::CellsContainerIterator it = source->GetOutput()->GetCells()->Begin();
      it != source->GetOutput()->GetCells()->End();
      ++it)
    {  
    TMesh::CellAutoPointer cell;
    source->GetOutput()->GetCell(it->Index(), cell);
    if (3 != cell->GetNumberOfPoints())
      { 
      std::cerr << "ERROR: All cells must be trianglar." << std::endl;  
      return EXIT_FAILURE;
      } 
    } 
 
  // Convert the triangle mesh to a simplex mesh.
  TConvert::Pointer convert = TConvert::New();
  convert->SetInput( source->GetOutput() );
  convert->Update();
 
  // Calculate the volume and area of the simplex mesh.
  TVolume::Pointer volume = TVolume::New();
  volume->SetSimplexMesh( convert->GetOutput() );
  volume->Compute();
 
  // Compare with the volume and area of an ideal sphere.
  std::cout << "Ideal Volume: "      << 4.0/3.0*M_PI*pow(5.0,3) << std::endl;
  std::cout << "Mesh Volume: "        << volume->GetVolume()    << std::endl;
  std::cout << "Ideal Surface Area: " << 4.0*M_PI*pow(5.0,2)    << std::endl;
  std::cout << "Mesh Surface Area: "  << volume->GetArea()      << std::endl;
 
  return EXIT_SUCCESS;
 
}
 
</source>
{{ITKCMakeLists|{{SUBPAGENAME}}}}

Latest revision as of 16:51, 7 June 2019

Warning: The media wiki content on this page is no longer maintained. The examples presented on the https://itk.org/Wiki/* pages likely require ITK version 4.13 or earlier releases. In many cases, the examples on this page no longer conform to the best practices for modern ITK versions.