1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
|
/*=========================================================================
Program: Insight Segmentation & Registration Toolkit
Module: $RCSfile: itkSpatialFunctionTest.cxx,v $
Language: C++
Date: $Date: 2003-09-10 14:30:10 $
Version: $Revision: 1.12 $
Copyright (c) Insight Software Consortium. All rights reserved.
See ITKCopyright.txt or http://www.itk.org/HTML/Copyright.htm for details.
This software is distributed WITHOUT ANY WARRANTY; without even
the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR
PURPOSE. See the above copyright notices for more information.
=========================================================================*/
#if defined(_MSC_VER)
#pragma warning ( disable : 4786 )
#endif
#include <stdio.h>
// Spatial function stuff
#include "itkSphereSpatialFunction.h"
int itkSpatialFunctionTest(int, char* [] )
{
// Change this parameter (and the positions, below) to work in higher or lower dimensions
const unsigned int dim = 3;
//---------Create and initialize a spatial function-----------
typedef itk::SphereSpatialFunction<dim> TFunctionType;
typedef TFunctionType::InputType TFunctionPositionType;
// Create and initialize a new sphere function
TFunctionType::Pointer spatialFunc = TFunctionType::New();
spatialFunc->SetRadius( 5 );
TFunctionPositionType center;
center[0]=10;
center[1]=10;
center[2]=10;
spatialFunc->SetCenter(center);
// Test the Get macros as well
spatialFunc->GetCenter();
spatialFunc->GetRadius();
std::cout << "Sphere spatial function created\n";
//----------------Test evaluation of funtion------------------
// We're going to evaluate it at the center of the sphere (10,10,10)
bool funcVal = spatialFunc->Evaluate(center);
printf("Sphere function value is %i\n", funcVal);
// The function should have returned a value of 1, since the center is inside
// the sphere
if(funcVal == 1)
return EXIT_SUCCESS;
else
return EXIT_FAILURE;
}
|
