ITK  4.9.0
Insight Segmentation and Registration Toolkit
WikiExamples/Meshes/QuadEdgeMeshParameterizationFilter.cxx
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;
// ** 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 **
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;
}