ITK/Examples/Meshes/QuadEdgeMeshNormalFilter: Difference between revisions

From KitwarePublic
< ITK‎ | Examples
Jump to navigationJump to search
No edit summary
Line 116: Line 116:
==CMakeLists.txt==
==CMakeLists.txt==
<source lang="cmake">
<source lang="cmake">
 
cmake_minimum_required(VERSION 2.6)
PROJECT(OtsuThresholdImageFilter)
FIND_PACKAGE(ITK REQUIRED)
INCLUDE(${ITK_USE_FILE})
ADD_EXECUTABLE(QuadEdgeMeshNormalFilter QuadEdgeMeshNormalFilter.cxx)
TARGET_LINK_LIBRARIES(QuadEdgeMeshNormalFilter ITKAlgorithms)
</source>
</source>

Revision as of 16:53, 13 December 2010

This example reads a mesh from a vtk file. It then compute normals for the mesh.

The sample sphere file is here.

QuadEdgeMeshNormalFilter.cxx

<source lang="cpp">

  1. include "itkVector.h"
  2. include "itkQuadEdgeMesh.h"
  3. include "itkVTKPolyDataReader.h"
  1. include "itkQuadEdgeMeshExtendedTraits.h"
  2. include "itkQuadEdgeMeshNormalFilter.h"
  3. include <stdlib.h>

int main( int argc, char* argv[] ) {

 if( argc < 2 )
   {
   std::cout <<"This program requires at least 1 argument" <<std::endl;
   std::cout <<"1- Input filename" <<std::endl;
   std::cout <<"2- Weight type" <<std::endl;
   std::cout <<"   * 0:  GOURAUD" <<std::endl;
   std::cout <<"   * 1:  THURMER" <<std::endl;
   std::cout <<"   * 2:  AREA" <<std::endl;
   return EXIT_FAILURE;
   }
 const unsigned int    Dimension = 3;
 typedef double        CoordType;
 typedef itk::QuadEdgeMesh< CoordType, Dimension > InputMeshType;
 typedef itk::Vector< CoordType, Dimension > VectorType;
 typedef itk::QuadEdgeMeshExtendedTraits <
   VectorType,
   Dimension,
   2,
   CoordType,
   CoordType,
   VectorType,
   bool,
   bool > Traits;
 typedef itk::QuadEdgeMesh < VectorType, Dimension, Traits > OutputMeshType;
 typedef itk::VTKPolyDataReader< InputMeshType > ReaderType;
 typedef itk::QuadEdgeMeshNormalFilter< InputMeshType, OutputMeshType > NormalFilterType;
 NormalFilterType::WeightType weight_type;
 int param = atoi( argv[2] );
 if( ( param < 0 ) || ( param > 2 ) )
   {
   std::cout <<"Weight type must be either: " <<std::endl;
   std::cout <<"   * 0:  GOURAUD" <<std::endl;
   std::cout <<"   * 1:  THURMER" <<std::endl;
   std::cout <<"   * 2:  AREA" <<std::endl;
   return EXIT_FAILURE;
   }
 else
   {
   switch( param )
     {
     default:
     case 0:
       weight_type = NormalFilterType::GOURAUD;
       break;
     case 1:
       weight_type = NormalFilterType::THURMER;
       break;
     case 2:
       weight_type = NormalFilterType::AREA;
       break;
     }
   }
 ReaderType::Pointer reader = ReaderType::New( );
 reader->SetFileName( argv[1] );
 try
   {
   reader->Update( );
   }
 catch( itk::ExceptionObject & excp )
   {
   std::cerr << "Exception thrown while reading the input file " << std::endl;
   std::cerr << excp << std::endl;
   return EXIT_FAILURE;
   }
 InputMeshType::Pointer mesh = reader->GetOutput( );
 NormalFilterType::Pointer normals = NormalFilterType::New( );
 normals->SetInput( mesh );
 normals->SetWeight( weight_type );
 try
   {
   normals->Update( );
   }
 catch( itk::ExceptionObject & excp )
   {
   std::cerr << excp << std::endl;
   return EXIT_FAILURE;
   }
 OutputMeshType::Pointer output = normals->GetOutput( );
 std::cout << normals << std::endl;
 return EXIT_SUCCESS;

} </source>

CMakeLists.txt

<source lang="cmake"> cmake_minimum_required(VERSION 2.6)

PROJECT(OtsuThresholdImageFilter)

FIND_PACKAGE(ITK REQUIRED) INCLUDE(${ITK_USE_FILE})

ADD_EXECUTABLE(QuadEdgeMeshNormalFilter QuadEdgeMeshNormalFilter.cxx) TARGET_LINK_LIBRARIES(QuadEdgeMeshNormalFilter ITKAlgorithms) </source>