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

itkHistogramToLogProbabilityImageFilter.h

Go to the documentation of this file.
00001 /*=========================================================================
00002 
00003   Program:   Insight Segmentation & Registration Toolkit
00004   Module:    $RCSfile: itkHistogramToLogProbabilityImageFilter.h,v $
00005   Language:  C++
00006   Date:      $Date: 2009-04-06 11:51:06 $
00007   Version:   $Revision: 1.11 $
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 __itkHistogramToLogProbabilityImageFilter_h
00018 #define __itkHistogramToLogProbabilityImageFilter_h
00019 
00020 #include "itkHistogramToImageFilter.h"
00021 
00022 namespace itk
00023 {
00024   
00046 namespace Function {  
00047 template< class TInput, class TOutput=double >
00048 class HistogramLogProbabilityFunction
00049 {
00050 public:
00051   
00052   //Probability function = Number of occurances in each bin /
00053   //   Total Number of occurances. 
00054   //
00055   // Returns pixels of float.. 
00056   typedef  TOutput  OutputPixelType;
00057   
00058   
00059   HistogramLogProbabilityFunction(): 
00060       m_TotalFrequency(1) {}
00061     
00062   ~HistogramLogProbabilityFunction() {};
00063   
00064   inline OutputPixelType operator()( const TInput & A ) const
00065     {
00066     if( A )
00067       {
00068       return static_cast<OutputPixelType>(vcl_log( static_cast<OutputPixelType>(A) / 
00069         static_cast<OutputPixelType>(m_TotalFrequency)) / vcl_log(2.0) );
00070       }
00071     else
00072       { // Check for Log 0. Always assume that the frequency is atleast 1.
00073       return static_cast<OutputPixelType>(vcl_log( static_cast<OutputPixelType>(A+1) / 
00074         static_cast<OutputPixelType>(m_TotalFrequency)) / vcl_log(2.0) );
00075       }
00076     }
00077 
00078   void SetTotalFrequency( unsigned long n ) 
00079     {
00080     m_TotalFrequency = n;
00081     }
00082   
00083   unsigned long GetTotalFrequency( ) const 
00084     {
00085     return m_TotalFrequency;
00086     }
00087 
00088 private:
00089   unsigned long  m_TotalFrequency;
00090 }; 
00091 }
00092 
00093 template <class THistogram, class TOutputPixel=double >
00094 class ITK_EXPORT HistogramToLogProbabilityImageFilter :
00095   public HistogramToImageFilter< THistogram, 
00096   Function::HistogramLogProbabilityFunction< unsigned long, TOutputPixel > > 
00097 {
00098 public:
00099   
00101   typedef HistogramToLogProbabilityImageFilter Self;
00102 
00104   typedef HistogramToImageFilter< THistogram, 
00105     Function::HistogramLogProbabilityFunction< unsigned long, TOutputPixel > >
00106                                                Superclass;
00107 
00108   typedef SmartPointer<Self>                   Pointer;
00109   typedef SmartPointer<const Self>             ConstPointer;
00110   
00112   itkTypeMacro( HistogramToLogProbabilityImageFilter, HistogramToImageFilter );
00113 
00115   itkNewMacro(Self);
00116 
00117 protected:
00118   HistogramToLogProbabilityImageFilter() {}
00119   virtual ~HistogramToLogProbabilityImageFilter() {}
00120   
00121 private:
00122   HistogramToLogProbabilityImageFilter(const Self&); //purposely not implemented
00123   void operator=(const Self&); //purposely not implemented
00124 
00125 };
00126 
00127 } // end namespace itk
00128 
00129 #endif
00130 

Generated at Thu May 28 10:12:15 2009 for ITK by doxygen 1.5.5 written by Dimitri van Heesch, © 1997-2000