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: 2006/08/01 11:31:24 $
00007   Version:   $Revision: 1.14 $
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 
00129   itkSetMacro( OutputMinimum, OutputPixelType );
00130   itkSetMacro( OutputMaximum, OutputPixelType );
00131   itkGetConstReferenceMacro( OutputMinimum, OutputPixelType );
00132   itkGetConstReferenceMacro( OutputMaximum, OutputPixelType );
00133 
00137   itkGetConstReferenceMacro( Scale, RealType );
00138   itkGetConstReferenceMacro( Shift, RealType );
00140 
00143   itkGetConstReferenceMacro( InputMinimum, InputPixelType );
00144   itkGetConstReferenceMacro( InputMaximum, InputPixelType );
00146 
00148   void BeforeThreadedGenerateData(void);
00149 
00151   void PrintSelf(std::ostream& os, Indent indent) const;
00152 
00153 #ifdef ITK_USE_CONCEPT_CHECKING
00154 
00155   itkConceptMacro(InputHasNumericTraitsCheck,
00156                   (Concept::HasNumericTraits<InputPixelType>));
00157   itkConceptMacro(OutputHasNumericTraitsCheck,
00158                   (Concept::HasNumericTraits<OutputPixelType>));
00159   itkConceptMacro(RealTypeMultiplyOperatorCheck,
00160                   (Concept::MultiplyOperator<RealType>));
00161   itkConceptMacro(RealTypeAdditiveOperatorsCheck,
00162                   (Concept::AdditiveOperators<RealType>));
00163 
00165 #endif
00166 
00167 protected:
00168   RescaleIntensityImageFilter();
00169   virtual ~RescaleIntensityImageFilter() {};
00170 
00171 private:
00172   RescaleIntensityImageFilter(const Self&); //purposely not implemented
00173   void operator=(const Self&); //purposely not implemented
00174 
00175   RealType m_Scale;
00176   RealType m_Shift;
00177 
00178   InputPixelType        m_InputMinimum;
00179   InputPixelType        m_InputMaximum;
00180 
00181   OutputPixelType       m_OutputMinimum;
00182   OutputPixelType       m_OutputMaximum;
00183 
00184 };
00185 
00186 
00187   
00188 } // end namespace itk
00189   
00190 #ifndef ITK_MANUAL_INSTANTIATION
00191 #include "itkRescaleIntensityImageFilter.txx"
00192 #endif
00193   
00194 #endif
00195 

Generated at Mon Mar 12 02:44:57 2007 for ITK by doxygen 1.5.1 written by Dimitri van Heesch, © 1997-2000