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

itkStandardDeviationProjectionImageFilter.h

Go to the documentation of this file.
00001 /*=========================================================================
00002 
00003   Program:   Insight Segmentation & Registration Toolkit
00004   Module:    $RCSfile: itkStandardDeviationProjectionImageFilter.h,v $
00005   Language:  C++
00006   Date:      $Date: 2009-04-01 14:36:35 $
00007   Version:   $Revision: 1.3 $
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 __itkStandardDeviationProjectionImageFilter_h
00018 #define __itkStandardDeviationProjectionImageFilter_h
00019 
00020 #include "itkProjectionImageFilter.h"
00021 #include "itkNumericTraits.h"
00022 #include "itkConceptChecking.h"
00023 
00024 namespace itk {
00047 namespace Function {
00048 template <class TInputPixel, class TAccumulate>
00049 class StandardDeviationAccumulator
00050 {
00051 public:
00052   typedef typename NumericTraits<TInputPixel>::RealType RealType;
00053 
00054   StandardDeviationAccumulator( unsigned long size )
00055     {
00056     m_Size = size;
00057     m_Values.reserve( size );
00058     }
00059   ~StandardDeviationAccumulator(){}
00060 
00061   inline void Initialize()
00062     {
00063     m_Sum = NumericTraits< TAccumulate >::Zero;
00064     m_Values.clear();
00065     }
00066 
00067   inline void operator()( const TInputPixel &input ) 
00068     {
00069     m_Sum = m_Sum + input;
00070     m_Values.push_back( input );
00071     }
00072 
00073   inline RealType GetValue()
00074     {
00075     // to avoid division by zero
00076     if( m_Size <= 1 )
00077       return NumericTraits<RealType>::Zero;
00078 
00079     typename NumericTraits<TInputPixel>::RealType mean =
00080       ((RealType) m_Sum) / m_Size;
00081     typename std::vector<TInputPixel>::iterator it;
00082     RealType squaredSum = NumericTraits<RealType>::Zero;
00083     for( it = m_Values.begin(); it != m_Values.end(); it++ )
00084       {
00085       squaredSum += vnl_math_sqr(*it - mean);
00086       }
00087     return vcl_sqrt( squaredSum / ( m_Size - 1) );
00088     }
00089 
00090   TAccumulate              m_Sum;
00091   unsigned long            m_Size;
00092   std::vector<TInputPixel> m_Values;
00093 };
00094 } // end namespace Function
00095 
00096 
00097 template <class TInputImage, 
00098           class TOutputImage, 
00099           class TAccumulate= ITK_TYPENAME 
00100           NumericTraits< ITK_TYPENAME TOutputImage::PixelType >
00101                            ::AccumulateType >
00102 class ITK_EXPORT StandardDeviationProjectionImageFilter :
00103     public
00104     ProjectionImageFilter<TInputImage, TOutputImage,
00105       Function::StandardDeviationAccumulator< ITK_TYPENAME 
00106                          TInputImage::PixelType, TAccumulate > >
00107 {
00108 public:
00109   typedef StandardDeviationProjectionImageFilter Self;
00110 
00111   typedef ProjectionImageFilter<TInputImage, TOutputImage, 
00112     Function::StandardDeviationAccumulator< typename 
00113                   TInputImage::PixelType, TAccumulate > > Superclass;
00114 
00115   typedef TInputImage                        InputImageType;
00116   typedef typename InputImageType::PixelType InputPixelType; 
00117 
00118   typedef SmartPointer<Self>        Pointer;
00119   typedef SmartPointer<const Self>  ConstPointer;
00120 
00122   itkTypeMacro(StandardDeviationProjectionImageFilter, ProjectionImageFilter);
00123 
00125   itkNewMacro(Self);
00126 
00127 #ifdef ITK_USE_CONCEPT_CHECKING
00128 
00129   itkConceptMacro(InputPixelToOutputPixelTypeGreaterAdditiveOperatorCheck,
00130     (Concept::AdditiveOperators<TAccumulate,
00131                                 InputPixelType,
00132                                 TAccumulate>));
00133   itkConceptMacro(InputHasNumericTraitsCheck,
00134     (Concept::HasNumericTraits<InputPixelType>));
00136 
00137   itkConceptMacro(AccumulateHasNumericTraitsCheck,
00138     (Concept::HasNumericTraits<TAccumulate>));
00139  
00141 #endif
00142 
00143 
00144 protected:
00145   StandardDeviationProjectionImageFilter() {}
00146   virtual ~StandardDeviationProjectionImageFilter() {}
00147 
00148 private:
00149   //purposely not implemented
00150   StandardDeviationProjectionImageFilter(const Self&); 
00151 
00152   void operator=(const Self&); //purposely not implemented
00153 
00154 }; // end StandardDeviationProjectionImageFilter
00155 
00156 } //end namespace itk
00157 
00158 #endif
00159 

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