Main Page   Groups   Namespace List   Class Hierarchy   Alphabetical List   Compound List   File List   Namespace Members   Compound Members   File Members   Concepts

itkGaussianBlurImageFunction.h

Go to the documentation of this file.
00001 /*=========================================================================
00002 
00003   Program:   Insight Segmentation & Registration Toolkit
00004   Module:    $RCSfile: itkGaussianBlurImageFunction.h,v $
00005   Language:  C++
00006   Date:      $Date: 2009-05-14 02:11:46 $
00007   Version:   $Revision: 1.17 $
00008 
00009   Copyright (c) Insight Software Consortium. All rights reserved.
00010   See ITKCopyright.txt or http://www.itk.org/HTML/Copyright.htm for details.
00011 
00012      This software is distributed WITHOUT ANY WARRANTY; without even
00013      the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR
00014      PURPOSE.  See the above copyright notices for more information.
00015 
00016 =========================================================================*/
00017 #ifndef __itkGaussianBlurImageFunction_h
00018 #define __itkGaussianBlurImageFunction_h
00019 
00020 #include "itkNeighborhoodOperatorImageFunction.h"
00021 #include "itkGaussianOperator.h"
00022 #include "itkGaussianSpatialFunction.h"
00023 #include "itkImageFunction.h"
00024 #include "itkCastImageFilter.h"
00025 
00026 namespace itk
00027 {
00028 
00038 template <class TInputImage,class TOutput=double>
00039 class ITK_EXPORT GaussianBlurImageFunction :
00040   public ImageFunction< TInputImage, TOutput >
00041 {
00042 public:
00043 
00045   typedef GaussianBlurImageFunction Self;
00046 
00048   typedef ImageFunction<TInputImage, TOutput> Superclass;
00049 
00051   typedef SmartPointer<Self>          Pointer;
00052   typedef SmartPointer<const Self>    ConstPointer;
00053 
00055   itkNewMacro(Self);
00056 
00058   itkTypeMacro( GaussianBlurImageFunction, ImageFunction );
00059 
00061   typedef TInputImage                                 InputImageType;
00062   typedef typename InputImageType::PixelType          InputPixelType;
00063   typedef typename Superclass::IndexType              IndexType;
00064   typedef typename Superclass::ContinuousIndexType    ContinuousIndexType;
00065 
00066 
00068   itkStaticConstMacro(ImageDimension, unsigned int, InputImageType::ImageDimension);
00069 
00070   typedef GaussianOperator<
00071     TOutput, itkGetStaticConstMacro(ImageDimension)>                            GaussianOperatorType;
00072   typedef Neighborhood<TOutput, itkGetStaticConstMacro(ImageDimension)>         NeighborhoodType;
00073   typedef FixedArray<NeighborhoodType,itkGetStaticConstMacro(ImageDimension)>   OperatorArrayType;
00074 
00075   typedef GaussianSpatialFunction<TOutput,1>                                    GaussianFunctionType;
00076   typedef typename GaussianFunctionType::Pointer                                GaussianFunctionPointer;
00077   typedef itk::Image<double,itkGetStaticConstMacro(ImageDimension)>             InternalImageType;
00078   typedef typename InternalImageType::Pointer                                   InternalImagePointer;
00079 
00080   typedef NeighborhoodOperatorImageFunction< InternalImageType, TOutput>        OperatorImageFunctionType;
00081   typedef typename OperatorImageFunctionType::Pointer                           OperatorImageFunctionPointer;
00082 
00083   typedef itk::CastImageFilter<InputImageType,InternalImageType>                CastImageFilterType;
00084   typedef typename CastImageFilterType::Pointer                                 CastImageFilterPointer;
00085 
00086   typedef itk::FixedArray< double, itkGetStaticConstMacro(ImageDimension) >     ErrorArrayType;
00087   typedef itk::FixedArray< double, itkGetStaticConstMacro(ImageDimension) >     ExtentArrayType;
00088   typedef itk::FixedArray< double, itkGetStaticConstMacro(ImageDimension) >     SigmaArrayType;
00089 
00091   typedef typename Superclass::PointType PointType;
00092 
00094    virtual TOutput Evaluate(const PointType& point) const;
00095 
00096 
00098   virtual TOutput EvaluateAtIndex( const IndexType & index ) const;
00099 
00101   virtual TOutput EvaluateAtContinuousIndex(
00102     const ContinuousIndexType & index ) const;
00103 
00110   void SetSigma( const double* sigma);
00111   void SetSigma( const float* sigma);
00112   void SetSigma( const double sigma);
00113   itkSetMacro( Sigma, SigmaArrayType );
00114   itkGetConstReferenceMacro( Sigma, SigmaArrayType );
00116 
00121   virtual void SetInputImage( const InputImageType * ptr );
00122 
00126   itkSetMacro( Extent, ExtentArrayType );
00127   itkGetConstReferenceMacro( Extent, ExtentArrayType );
00128   void SetExtent( const double* extent);
00129   void SetExtent( const double extent);
00131 
00135   itkSetMacro( MaximumError, ErrorArrayType );
00136   itkGetConstReferenceMacro( MaximumError, ErrorArrayType );
00138 
00143   itkSetMacro( MaximumKernelWidth, int );
00144   itkGetConstMacro( MaximumKernelWidth, int );
00146 
00152   itkSetMacro( UseImageSpacing, bool );
00153   itkGetConstMacro( UseImageSpacing, bool );
00154   itkBooleanMacro( UseImageSpacing );
00156 
00157 
00158 protected:
00159   GaussianBlurImageFunction();
00160   GaussianBlurImageFunction( const Self& ){};
00161 
00162   ~GaussianBlurImageFunction(){};
00163 
00164   void operator=( const Self& ){};
00165   void PrintSelf(std::ostream& os, Indent indent) const;
00166 
00167   void RecomputeGaussianKernel();
00168   void RecomputeContinuousGaussianKernel(const double* offset) const;
00169 
00170 private:
00171 
00172   virtual TOutput EvaluateAtIndex( 
00173     const IndexType & index, const OperatorArrayType & operatorArray ) const;
00174 
00175   SigmaArrayType                m_Sigma;
00176   OperatorImageFunctionPointer  m_OperatorImageFunction;
00177   mutable OperatorArrayType     m_OperatorArray;
00178   mutable OperatorArrayType     m_ContinuousOperatorArray;
00179   InternalImagePointer          m_InternalImage;
00180   CastImageFilterPointer        m_Caster;
00181 
00185   ErrorArrayType                m_MaximumError;
00186   ExtentArrayType               m_Extent;
00187 
00190   int m_MaximumKernelWidth;
00191 
00193   unsigned int m_FilterDimensionality;
00194 
00196   bool m_UseImageSpacing;
00197 
00199   GaussianFunctionPointer m_GaussianFunction;
00200 };
00201 
00202 } // end namespace itk
00203 
00204 // Define instantiation macro for this template.
00205 #define ITK_TEMPLATE_GaussianBlurImageFunction(_, EXPORT, x, y) namespace itk { \
00206   _(2(class EXPORT GaussianBlurImageFunction< ITK_TEMPLATE_2 x >)) \
00207   namespace Templates { typedef GaussianBlurImageFunction< ITK_TEMPLATE_2 x > \
00208                                                   GaussianBlurImageFunction##y; } \
00209   }
00210 
00211 #if ITK_TEMPLATE_EXPLICIT
00212 # include "Templates/itkGaussianBlurImageFunction+-.h"
00213 #endif
00214 
00215 #if ITK_TEMPLATE_TXX
00216 # include "itkGaussianBlurImageFunction.txx"
00217 #endif
00218 
00219 #endif
00220 

Generated at Tue Sep 15 03:01:14 2009 for ITK by doxygen 1.5.8 written by Dimitri van Heesch, © 1997-2000