|
|
| (3 intermediate revisions by 2 users not shown) |
| Line 1: |
Line 1: |
| This example reads a mesh from a vtk file. It then compute the planar parameterization
| | {{warning|1=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. |
| | }} |
|
| |
|
| The sample sphere file is [http://www.itk.org/Wiki/images/8/89/Sphere.vtk.zip here].
| | [https://itk.org/ITKExamples[ITK Sphinx Examples]] |
| | |
| ==QuadEdgeMeshNormalFilter.cxx==
| |
| <source lang="cpp">
| |
| | |
| #include "itkVTKPolyDataReader.h"
| |
| #include "itkVTKPolyDataWriter.h"
| |
| #include "VNLIterativeSparseSolverTraits.h"
| |
| #include "itkParameterizationQuadEdgeMeshFilter.h"
| |
| | |
| int main( int argc, char* argv[] )
| |
| {
| |
| // ** ERROR MESSAGE AND HELP ** //
| |
| if( argc < 5 )
| |
| {
| |
| std::cout <<"Requires 4 arguments: " << std::endl;
| |
| std::cout <<"1-Input file name " << std::endl;
| |
| std::cout <<"2-Border Type" << std::endl;
| |
| std::cout <<" * 0: SQUARE" << std::endl;
| |
| std::cout <<" * 1: DISK" << std::endl;
| |
| std::cout <<"3-CoefficientType Type" << std::endl;
| |
| std::cout <<" * 0: OnesMatrixCoefficients" << std::endl;
| |
| std::cout <<" * 1: InverseEuclideanDistanceMatrixCoefficients" << std::endl;
| |
| std::cout <<" * 2: ConformalMatrixCoefficients" << std::endl;
| |
| std::cout <<" * 3: AuthalicMatrixCoefficients" << std::endl;
| |
| std::cout <<" * 4: HarmonicMatrixCoefficients" << std::endl;
| |
| std::cout <<"4-Output file name " << std::endl;
| |
| | |
| return EXIT_FAILURE;
| |
| }
| |
| | |
| | |
| // ** TYPEDEF **
| |
| typedef double Coord;
| |
| | |
| typedef itk::QuadEdgeMesh< Coord, 3 > MeshType;
| |
| typedef itk::VTKPolyDataReader< MeshType > ReaderType;
| |
| typedef itk::VTKPolyDataWriter< MeshType > WriterType;
| |
| typedef itk::BorderQuadEdgeMeshFilter< MeshType, MeshType > BorderTransformType;
| |
| typedef VNLIterativeSparseSolverTraits< Coord > SolverTraits;
| |
| | |
| typedef itk::ParameterizationQuadEdgeMeshFilter< MeshType, MeshType, SolverTraits > ParametrizationType;
| |
| | |
| | |
| // ** READ THE FILE IN **
| |
| 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;
| |
| }
| |
| | |
| MeshType::Pointer mesh = reader->GetOutput( );
| |
| | |
| // ** CHOSE< COMPUTE AND SET BORDER TRANSFORM **
| |
| BorderTransformType::Pointer border_transform = BorderTransformType::New( );
| |
| border_transform->SetInput( mesh );
| |
| // two following line for coverage
| |
| border_transform->SetRadius( border_transform->GetRadius() );
| |
| | |
| int border;
| |
| std::stringstream ssout( argv[2] );
| |
| ssout >>border;
| |
| switch( border ) // choose border type
| |
| {
| |
| case 0: // square shaped domain
| |
| border_transform->SetTransformType( BorderTransformType::SQUARE_BORDER_TRANSFORM );
| |
| break;
| |
| case 1: // disk shaped domain
| |
| border_transform->SetTransformType( BorderTransformType::DISK_BORDER_TRANSFORM );
| |
| break;
| |
| default: // handle .... user ....
| |
| std::cerr << "2nd argument must be " << std::endl;
| |
| std::cerr << "0 for SQUARE BORDER TRANSFORM or "
| |
| << "1 for DISK BORDER TRANSFORM" << std::endl;
| |
| return EXIT_FAILURE;
| |
| }
| |
| std::cout << "Transform type is: " << border_transform->GetTransformType( );
| |
| std::cout << std::endl;
| |
| | |
| // ** CHOOSE AND SET BARYCENTRIC WEIGHTS **
| |
| itk::OnesMatrixCoefficients< MeshType > coeff0;
| |
| itk::InverseEuclideanDistanceMatrixCoefficients< MeshType > coeff1;
| |
| itk::ConformalMatrixCoefficients< MeshType > coeff2;
| |
| itk::AuthalicMatrixCoefficients< MeshType > coeff3;
| |
| itk::HarmonicMatrixCoefficients< MeshType > coeff4;
| |
| | |
| ParametrizationType::Pointer param = ParametrizationType::New( );
| |
| param->SetInput( mesh );
| |
| param->SetBorderTransform( border_transform );
| |
| | |
| int param_type;
| |
| std::stringstream ssout3( argv[3] );
| |
| ssout3 >> param_type;
| |
| | |
| switch( param_type )
| |
| {
| |
| case 0:
| |
| param->SetCoefficientsMethod( &coeff0 );
| |
| break;
| |
| case 1:
| |
| param->SetCoefficientsMethod( &coeff1 );
| |
| break;
| |
| case 2:
| |
| param->SetCoefficientsMethod( &coeff2 );
| |
| break;
| |
| case 3:
| |
| param->SetCoefficientsMethod( &coeff3 );
| |
| break;
| |
| case 4:
| |
| param->SetCoefficientsMethod( &coeff4 );
| |
| break;
| |
| default:
| |
| std::cerr << "3rd argument must be " << std::endl;
| |
| std::cerr << "0, 1, 2, 3 or 4" << std::endl;
| |
| std::cerr << "Here it is: " << param_type << std::endl;
| |
| return EXIT_FAILURE;
| |
| }
| |
| | |
| // ** PROCESS **
| |
| param->Update( );
| |
| MeshType::Pointer output = param->GetOutput( );
| |
| | |
| // ** WRITE OUTPUT **
| |
| WriterType::Pointer writer = WriterType::New( );
| |
| writer->SetInput( param->GetOutput( ) );
| |
| writer->SetFileName( argv[4] );
| |
| writer->Update( );
| |
| | |
| // ** PRINT **
| |
| std::cout << "BorderTransform: \n" << border_transform;
| |
| std::cout << "Parametrization: \n" << param;
| |
| | |
| return EXIT_SUCCESS;
| |
| }
| |
| </source>
| |
| {{ITKCMakeLists|QuadEdgeMeshParameterizationFilter}}
| |
