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
67
68
69
70
71
72
73
74
75
76
77
78
79
|
/*=========================================================================
Program: Insight Segmentation & Registration Toolkit
Module: $RCSfile: itkGaussianDensityFunctionTest.cxx,v $
Language: C++
Date: $Date: 2007-03-07 14:05:39 $
Version: $Revision: 1.5 $
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 "itkGaussianDensityFunction.h"
#include "itkVector.h"
int itkGaussianDensityFunctionTest(int, char* [] )
{
std::cout << "itkGaussianDensityFunctionTest Test \n \n";
typedef itk::Vector< double, 1 > MeasurementVectorType;
typedef itk::Statistics::GaussianDensityFunction<
MeasurementVectorType
> DensityFunctionType;
typedef DensityFunctionType::MeasurementVectorType MeasurementVectorType;
typedef DensityFunctionType::CovarianceType CovarianceType;
typedef DensityFunctionType::MeanType MeanType;
MeanType mean(1); // size = 1 because this is a scalar case.
CovarianceType covariance;
covariance.SetSize( 1, 1 ); // because this is a scalar case
const double Sigma = 7.0;
mean[0] = 19.0;
covariance[0][0] = Sigma * Sigma;
DensityFunctionType::Pointer densityFunction = DensityFunctionType::New();
densityFunction->SetMean( &mean );
densityFunction->SetCovariance( &covariance );
MeasurementVectorType inputValue;
inputValue[0] = mean[0];
const double value1 = densityFunction->Evaluate( inputValue );
const double nPI = 4.0 * atan( 1.0 );
const double gaussianNorm = 1.0 / ( sqrt( 2 * nPI ) * Sigma );
if( fabs( gaussianNorm - value1 ) > 1e-6 )
{
std::cerr << "ERROR in computation of the Gaussian " << std::endl;
return EXIT_FAILURE;
}
std::cout << "Evaluate at the mean = " << value1 << std::endl;
std::cout << "Expected at the mean = " << gaussianNorm << std::endl;
std::cout << "Test passed." << std::endl;
return EXIT_SUCCESS;
}
|
