00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018 #ifndef __itkGaussianSpatialObject_h
00019 #define __itkGaussianSpatialObject_h
00020
00021 #include "itkAffineTransform.h"
00022 #include "itkFixedArray.h"
00023 #include "itkSpatialObject.h"
00024 #include "itkEllipseSpatialObject.h"
00025
00026 namespace itk
00027 {
00028
00045 template < unsigned int TDimension = 3 >
00046 class GaussianSpatialObject
00047 : public SpatialObject< TDimension >
00048 {
00049
00050 public:
00051
00052 typedef GaussianSpatialObject Self;
00053 typedef double ScalarType;
00054 typedef SmartPointer < Self > Pointer;
00055 typedef SmartPointer < const Self > ConstPointer;
00056 typedef SpatialObject< TDimension > Superclass;
00057 typedef SmartPointer< Superclass > SuperclassPointer;
00058 typedef typename Superclass::PointType PointType;
00059 typedef typename Superclass::TransformType TransformType;
00060 typedef typename Superclass::BoundingBoxType BoundingBoxType;
00061
00062 itkStaticConstMacro(NumberOfDimensions, unsigned int,
00063 TDimension);
00064
00065 itkNewMacro( Self );
00066 itkTypeMacro( GaussianSpatialObject, SpatialObject );
00067
00071 itkSetMacro(Radius,ScalarType);
00072 itkGetConstReferenceMacro(Radius,ScalarType);
00074
00077 itkSetMacro(Maximum,ScalarType);
00078 itkGetConstReferenceMacro(Maximum,ScalarType);
00080
00084 ScalarType SquaredZScore( const PointType& point ) const;
00085
00087 virtual bool ValueAt( const PointType & point, ScalarType & value,
00088 unsigned int depth=0,
00089 char * name=NULL) const;
00090
00092 virtual bool IsEvaluableAt( const PointType & , unsigned int,
00093 char *) const { return true; }
00094
00096 virtual bool IsInside( const PointType & point,
00097 unsigned int depth,
00098 char * name) const;
00099
00103 virtual bool IsInside( const PointType & point) const;
00104
00107 virtual bool ComputeLocalBoundingBox() const;
00108
00111 typename EllipseSpatialObject< TDimension >::Pointer GetEllipsoid() const;
00112
00113 protected:
00114
00115 GaussianSpatialObject( void );
00116 ~GaussianSpatialObject( void );
00117
00118 ScalarType m_Maximum;
00119 ScalarType m_Radius;
00120
00122 virtual void PrintSelf( std::ostream& os, Indent indent ) const;
00123
00124 };
00125
00126 }
00127
00128 #ifndef ITK_MANUAL_INSTANTIATION
00129 #include "itkGaussianSpatialObject.txx"
00130 #endif
00131
00132 #endif // __itkGaussianSpatialObject_h
00133