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

itkRescaleIntensityImageFilter.h

Go to the documentation of this file.
00001 /*=========================================================================
00002 
00003   Program:   Insight Segmentation & Registration Toolkit
00004   Module:    $RCSfile: itkRescaleIntensityImageFilter.h,v $
00005   Language:  C++
00006   Date:      $Date: 2007/09/27 11:36:41 $
00007   Version:   $Revision: 1.15 $
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 __itkRescaleIntensityImageFilter_h
00018 #define __itkRescaleIntensityImageFilter_h
00019 
00020 #include "itkUnaryFunctorImageFilter.h"
00021 
00022 namespace itk
00023 {
00024 
00025 // This functor class applies a linear transformation A.x + B 
00026 // to input values.
00027 namespace Functor {  
00028  
00029 template< typename TInput, typename  TOutput>
00030 class IntensityLinearTransform
00031 {
00032 public:
00033   typedef typename NumericTraits< TInput >::RealType RealType;
00034   IntensityLinearTransform()
00035   {
00036     m_Factor = 1.0;
00037     m_Offset = 0.0;
00038     m_Minimum = NumericTraits<TOutput>::NonpositiveMin();
00039     m_Maximum = NumericTraits<TOutput>::max();
00040   }
00041   ~IntensityLinearTransform() {}
00042   void SetFactor( RealType a ) { m_Factor = a; }
00043   void SetOffset( RealType b ) { m_Offset = b; }
00044   void SetMinimum( TOutput min ) { m_Minimum = min; }
00045   void SetMaximum( TOutput max ) { m_Maximum = max; }
00046   bool operator!=( const IntensityLinearTransform & other ) const
00047   {
00048     if( m_Factor != other.m_Factor ||
00049         m_Offset != other.m_Offset ||
00050         m_Maximum != other.m_Maximum    ||
00051         m_Minimum != other.m_Minimum )
00052       {
00053       return true;
00054       }
00055     return false;
00056    }
00057   bool operator==( const IntensityLinearTransform & other ) const
00058   {
00059     return !(*this != other);
00060   }
00061   inline TOutput operator()( const TInput & x )
00062   {
00063     RealType value  = static_cast<RealType>(x) * m_Factor + m_Offset;
00064     TOutput  result = static_cast<TOutput>( value );
00065     result = ( result > m_Maximum ) ? m_Maximum : result;
00066     result = ( result < m_Minimum ) ? m_Minimum : result;
00067     return result;
00068   }
00069 private:
00070   RealType m_Factor;
00071   RealType m_Offset;
00072   TOutput  m_Maximum;
00073   TOutput  m_Minimum;
00074 }; 
00075 
00076 }  // end namespace functor
00077 
00078 
00104 template <typename  TInputImage, typename  TOutputImage=TInputImage>
00105 class ITK_EXPORT RescaleIntensityImageFilter :
00106     public
00107 UnaryFunctorImageFilter<TInputImage,TOutputImage, 
00108                         Functor::IntensityLinearTransform< 
00109   typename TInputImage::PixelType, 
00110   typename TOutputImage::PixelType>   >
00111 {
00112 public:
00114   typedef RescaleIntensityImageFilter  Self;
00115   typedef UnaryFunctorImageFilter<TInputImage,TOutputImage, 
00116                                   Functor::IntensityLinearTransform< 
00117     typename TInputImage::PixelType, 
00118     typename TOutputImage::PixelType> >  Superclass;
00119   typedef SmartPointer<Self>   Pointer;
00120   typedef SmartPointer<const Self>  ConstPointer;
00121 
00122   typedef typename TOutputImage::PixelType OutputPixelType;
00123   typedef typename TInputImage::PixelType  InputPixelType;
00124   typedef typename NumericTraits<InputPixelType>::RealType RealType;
00125 
00127   itkNewMacro(Self);
00128 
00130   itkTypeMacro(RescaleIntensityImageFilter, 
00131                UnaryFunctorImageFilter);
00132 
00133   itkSetMacro( OutputMinimum, OutputPixelType );
00134   itkSetMacro( OutputMaximum, OutputPixelType );
00135   itkGetConstReferenceMacro( OutputMinimum, OutputPixelType );
00136   itkGetConstReferenceMacro( OutputMaximum, OutputPixelType );
00137 
00141   itkGetConstReferenceMacro( Scale, RealType );
00142   itkGetConstReferenceMacro( Shift, RealType );
00144 
00147   itkGetConstReferenceMacro( InputMinimum, InputPixelType );
00148   itkGetConstReferenceMacro( InputMaximum, InputPixelType );
00150 
00152   void BeforeThreadedGenerateData(void);
00153 
00155   void PrintSelf(std::ostream& os, Indent indent) const;
00156 
00157 #ifdef ITK_USE_CONCEPT_CHECKING
00158 
00159   itkConceptMacro(InputHasNumericTraitsCheck,
00160                   (Concept::HasNumericTraits<InputPixelType>));
00161   itkConceptMacro(OutputHasNumericTraitsCheck,
00162                   (Concept::HasNumericTraits<OutputPixelType>));
00163   itkConceptMacro(RealTypeMultiplyOperatorCheck,
00164                   (Concept::MultiplyOperator<RealType>));
00165   itkConceptMacro(RealTypeAdditiveOperatorsCheck,
00166                   (Concept::AdditiveOperators<RealType>));
00167 
00169 #endif
00170 
00171 protected:
00172   RescaleIntensityImageFilter();
00173   virtual ~RescaleIntensityImageFilter() {};
00174 
00175 private:
00176   RescaleIntensityImageFilter(const Self&); //purposely not implemented
00177   void operator=(const Self&); //purposely not implemented
00178 
00179   RealType m_Scale;
00180   RealType m_Shift;
00181 
00182   InputPixelType        m_InputMinimum;
00183   InputPixelType        m_InputMaximum;
00184 
00185   OutputPixelType       m_OutputMinimum;
00186   OutputPixelType       m_OutputMaximum;
00187 
00188 };
00189 
00190 
00191   
00192 } // end namespace itk
00193   
00194 #ifndef ITK_MANUAL_INSTANTIATION
00195 #include "itkRescaleIntensityImageFilter.txx"
00196 #endif
00197   
00198 #endif
00199 

Generated at Tue Jul 29 22:03:54 2008 for ITK by doxygen 1.5.1 written by Dimitri van Heesch, © 1997-2000