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
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
|
/*=========================================================================
Program: Insight Segmentation & Registration Toolkit
Module: $RCSfile: itkListSampleBase.h,v $
Language: C++
Date: $Date: 2010-02-02 13:45:16 $
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.
=========================================================================*/
#ifndef __itkListSampleBase_h
#define __itkListSampleBase_h
#include "itkSample.h"
#include <vector>
namespace itk {
namespace Statistics {
/** \class ListSampleBase
* \brief This class is the base class for Samples that store measurements in a list
*
* ListSampleBase stores measurements in a list type structure (as
* opposed to a Histogram, etc.). ListSampleBase allows duplicate
* measurements. ListSampleBase is not sorted.
*
* ListSampleBase does not allow the user to specify the frequency of
* a measurement directly. The GetFrequency() methods returns 1 if
* the measurement exists in the list, 0 otherwise.
*
* Recent API changes:
* The static const macro to get the length of a measurement vector,
* 'MeasurementVectorSize' has been removed to allow the length of a measurement
* vector to be specified at run time. Please use the function
* GetMeasurementVectorSize() instead.
*
*\sa Sample, Histogram
*/
template< class TMeasurementVector >
class ITK_EXPORT ListSampleBase : public Sample< TMeasurementVector >
{
public:
/** Standard class typedef. */
typedef ListSampleBase Self;
typedef Sample< TMeasurementVector > Superclass;
typedef SmartPointer< Self > Pointer;
typedef SmartPointer<const Self> ConstPointer;
/** Standard macros */
itkTypeMacro(ListSampleBase, Sample);
/** Typedefs inherited from the superclass */
typedef typename Superclass::MeasurementVectorType MeasurementVectorType;
typedef typename Superclass::MeasurementVectorSizeType MeasurementVectorSizeType;
typedef typename Superclass::MeasurementType MeasurementType;
typedef typename Superclass::FrequencyType FrequencyType;
typedef typename Superclass::InstanceIdentifier InstanceIdentifier;
/** Vector of InstanceIdentifiers used for returning search
* results. */
typedef std::vector< InstanceIdentifier > SearchResultVectorType;
/** Search for measurements within the specified radius of a search
* point. A vector of InstanceIdentifiers is returned. */
inline void Search(MeasurementVectorType center, double radius,
SearchResultVectorType& result) const
{
if (radius == 0.0)
{
itkGenericExceptionMacro("Search radius should be greater than zero.");
}
unsigned int j;
double squaredRadius;
double distance;
double coordinateDistance;
MeasurementVectorType tempVector;
squaredRadius = radius * radius;
result.clear();
for ( InstanceIdentifier identifier = 0; identifier < this->Size(); ++identifier )
{
distance = 0.0;
tempVector = this->GetMeasurementVector( identifier );
for (j = 0; j < this->GetMeasurementVectorSize() && distance < squaredRadius; j++)
{
coordinateDistance = (double)tempVector[j] - center[j];
if (vnl_math_abs(coordinateDistance) > radius )
{
distance = squaredRadius;
}
}
for (j = 0; j < this->GetMeasurementVectorSize() && distance < squaredRadius; j++)
{
coordinateDistance = (double)tempVector[j] - center[j];
distance += coordinateDistance * coordinateDistance;
}
if (distance < squaredRadius)
{
result.push_back( identifier );
}
}
}
protected:
ListSampleBase() {}
virtual ~ListSampleBase() {}
void PrintSelf(std::ostream& os, Indent indent) const
{
Superclass::PrintSelf(os,indent);
}
private:
ListSampleBase(const Self&); //purposely not implemented
void operator=(const Self&); //purposely not implemented
};
} // end of namespace Statistics
} // end of namespace itk
#endif
|
