/usr/include/InsightToolkit/SpatialObject/itkGaussianSpatialObject.h is in libinsighttoolkit3-dev 3.20.1-1.
This file is owned by root:root, with mode 0o644.
The actual contents of the file can be viewed below.
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 | /*=========================================================================
Program: Insight Segmentation & Registration Toolkit
Module: itkGaussianSpatialObject.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 __itkGaussianSpatialObject_h
#define __itkGaussianSpatialObject_h
#include "itkAffineTransform.h"
#include "itkFixedArray.h"
#include "itkSpatialObject.h"
#include "itkEllipseSpatialObject.h"
namespace itk
{
/** \class GaussianSpatialObject
*
* \brief Represents a multivariate Gaussian function.
*
* The Gaussian function G(x) is given by
* \f[
* G(\vec{x}) = m e^{-\|\S^{-1} \vec{x}\|^2 / 2},
* \f]
* where m is a scaling factor set by SetMaximum(), and \f$\S\f$ is the
* (invertible) matrix associated to the IndexToObjectTransform of the object
* multiplied by the Sigma parameter. If \f$\S\f$ is symmetric and positive
* definite, and m is chosen so that the integral of G(x) is 1, then G will
* denote a normal distribution with mean 0 and covariance matrix \f$\S \times
* Sigma\f$.
*/
template < unsigned int TDimension = 3 >
class ITK_EXPORT GaussianSpatialObject
: public SpatialObject< TDimension >
{
public:
typedef GaussianSpatialObject Self;
typedef double ScalarType;
typedef SmartPointer < Self > Pointer;
typedef SmartPointer < const Self > ConstPointer;
typedef SpatialObject< TDimension > Superclass;
typedef SmartPointer< Superclass > SuperclassPointer;
typedef typename Superclass::PointType PointType;
typedef typename Superclass::TransformType TransformType;
typedef typename Superclass::BoundingBoxType BoundingBoxType;
itkStaticConstMacro(NumberOfDimensions, unsigned int,
TDimension);
itkNewMacro( Self );
itkTypeMacro( GaussianSpatialObject, SpatialObject );
/** The Radius determines the bounding box, and which points are
* considered to be inside the SpatialObject. All points with
* z-score less than the radius are in the object. */
itkSetMacro(Radius,ScalarType);
itkGetConstReferenceMacro(Radius,ScalarType);
/** The Sigma parameter determines the fallout of the Gaussian inside of the
* region defined by the Radius parameter. */
itkSetMacro(Sigma,ScalarType);
itkGetConstReferenceMacro(Sigma,ScalarType);
/** The maximum value of the Gaussian (its value at the origin of
* the spatial object coordinate system). */
itkSetMacro(Maximum,ScalarType);
itkGetConstReferenceMacro(Maximum,ScalarType);
/** If the matrix S is returned by
* this->GetIndexToObjectTransform()->GetMatrix(), then SquaredZScore(x)
* returns |Sx| squared. */
ScalarType SquaredZScore( const PointType& point ) const;
/** Returns the value of the Gaussian at the given point. */
virtual bool ValueAt( const PointType & point, ScalarType & value,
unsigned int depth=0,
char * name=NULL) const;
/** Return true if the object provides a method to evaluate the value
* at the specified point, false otherwise. */
virtual bool IsEvaluableAt( const PointType & point,
unsigned int depth=0,
char * name=NULL) const;
/** Test whether a point is inside or outside the object */
virtual bool IsInside( const PointType & point,
unsigned int depth,
char * name) const;
/** Test whether a point is inside or outside the object
* For computational speed purposes, it is faster if the method does not
* check the name of the class and the current depth */
virtual bool IsInside( const PointType & point) const;
/** This function needs to be called every time one of the object's
* components is changed. */
virtual bool ComputeLocalBoundingBox() const;
/** Returns the sigma=m_Radius level set of the Gaussian function, as an
* EllipseSpatialObject. */
typename EllipseSpatialObject< TDimension >::Pointer GetEllipsoid() const;
protected:
GaussianSpatialObject(const Self&); //purposely not implemented
void operator=(const Self&); //purposely not implemented
GaussianSpatialObject( void );
~GaussianSpatialObject( void );
ScalarType m_Maximum;
ScalarType m_Radius;
ScalarType m_Sigma;
/** Print the object information in a stream. */
virtual void PrintSelf( std::ostream& os, Indent indent ) const;
};
} // end namespace itk
#ifndef ITK_MANUAL_INSTANTIATION
#include "itkGaussianSpatialObject.txx"
#endif
#endif // __itkGaussianSpatialObject_h
|