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
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
|
/*=========================================================================
Program: Insight Segmentation & Registration Toolkit
Module: itkGaussianImageSource.h
Language: C++
Date: $Date$
Version: $Revision$
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 __itkGaussianImageSource_h
#define __itkGaussianImageSource_h
#include "itkImageSource.h"
#include "itkFixedArray.h"
#include "itkSize.h"
namespace itk
{
/** \class GaussianImageSource
* \brief Generate an n-dimensional image of a Gaussian.
*
* GaussianImageSource generates an image of a Gaussian.
* m_Normalized determines whether or not the Gaussian is normalized
* (whether or not the sum over infinite space is 1.0)
* When creating an image, it is preferable to _not_ normalize the Gaussian
* m_Scale scales the output of the Gaussian to span a range
* larger than 0->1, and is typically set to the maximum value
* of the output data type (for instance, 255 for uchars)
*
* The output image may be of any dimension.
*
* \ingroup DataSources
*/
template <typename TOutputImage>
class ITK_EXPORT GaussianImageSource : public ImageSource<TOutputImage>
{
public:
/** Standard class typedefs. */
typedef GaussianImageSource Self;
typedef ImageSource<TOutputImage> Superclass;
typedef SmartPointer<Self> Pointer;
typedef SmartPointer<const Self> ConstPointer;
/** Typedef for the output image PixelType. */
typedef typename TOutputImage::PixelType OutputImagePixelType;
/** Typedef to describe the output image region type. */
typedef typename TOutputImage::RegionType OutputImageRegionType;
/** Spacing typedef support. Spacing holds the size of a pixel. The
* spacing is the geometric distance between image samples. */
typedef typename TOutputImage::SpacingType SpacingType;
/** Origin typedef support. The origin is the geometric coordinates
* of the index (0,0). */
typedef typename TOutputImage::PointType PointType;
/** Direction typedef support. The direction is the direction
* cosines of the image. */
typedef typename TOutputImage::DirectionType DirectionType;
/** Dimensionality of the output image */
itkStaticConstMacro(NDimensions, unsigned int, TOutputImage::ImageDimension);
/** Type used to store gaussian parameters. */
typedef FixedArray<double, itkGetStaticConstMacro(NDimensions)> ArrayType;
/** Size type matches that used for images */
typedef typename TOutputImage::SizeType SizeType;
typedef typename TOutputImage::SizeValueType SizeValueType;
/** Run-time type information (and related methods). */
itkTypeMacro(GaussianImageSource,ImageSource);
/** Method for creation through the object factory. */
itkNewMacro(Self);
/** Specify the size of the output image. */
virtual void SetSize( const SizeValueType * values);
/** Specify the size of the output image. */
virtual void SetSize( const SizeType values);
/** Get the size of the output image. */
itkGetVectorMacro(Size,const SizeValueType,NDimensions);
/** Specify the spacing of the output image. */
itkSetMacro(Spacing, SpacingType);
virtual void SetSpacing( const float* values);
virtual void SetSpacing( const double* values);
/** Get the spacing of the output image. */
itkGetConstReferenceMacro(Spacing,SpacingType);
/** Specify the origin of the output image. */
itkSetMacro(Origin, PointType);
virtual void SetOrigin( const float* values);
virtual void SetOrigin( const double* values);
/** Get the origin of the output image. */
itkGetConstReferenceMacro(Origin,PointType);
/** Specify the direction of the output image. */
itkSetMacro(Direction, DirectionType);
itkGetConstReferenceMacro(Direction, DirectionType);
/** Gets and sets for gaussian parameters */
itkSetMacro(Scale, double);
itkGetConstReferenceMacro(Scale, double);
itkSetMacro(Normalized, bool);
itkGetConstReferenceMacro(Normalized, bool);
itkSetMacro(Sigma, ArrayType);
itkGetConstReferenceMacro(Sigma, ArrayType);
itkSetMacro(Mean, ArrayType);
itkGetConstReferenceMacro(Mean, ArrayType);
protected:
GaussianImageSource();
~GaussianImageSource();
void PrintSelf(std::ostream& os, Indent indent) const;
void GenerateData();
virtual void GenerateOutputInformation();
private:
GaussianImageSource(const GaussianImageSource&); //purposely not implemented
void operator=(const GaussianImageSource&); //purposely not implemented
SizeValueType m_Size[NDimensions]; //size of the output image
SpacingType m_Spacing; //spacing
PointType m_Origin; //origin
DirectionType m_Direction; // direction
/** Parameters for the Gaussian. */
/** The standard deviation in each direction. */
ArrayType m_Sigma;
/** The mean in each direction. */
ArrayType m_Mean;
/** A scale factor multiplied by the true value of the Gaussian. */
double m_Scale;
/** Whether or not to normalize the Gaussian. */
bool m_Normalized;
};
} // end namespace itk
#ifndef ITK_MANUAL_INSTANTIATION
#include "itkGaussianImageSource.txx"
#endif
#endif
|
