00001 /*========================================================================= 00002 00003 Program: Insight Segmentation & Registration Toolkit 00004 Module: $RCSfile: itkComplexToModulusImageFilter.h,v $ 00005 Language: C++ 00006 Date: $Date: 2006/04/03 14:33:37 $ 00007 Version: $Revision: 1.7 $ 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 __itkComplexToModulusImageFilter_h 00018 #define __itkComplexToModulusImageFilter_h 00019 00020 #include "itkUnaryFunctorImageFilter.h" 00021 #include "vnl/vnl_math.h" 00022 00023 namespace itk 00024 { 00025 00031 namespace Function { 00032 00033 template< class TInput, class TOutput> 00034 class ComplexToModulus 00035 { 00036 public: 00037 ComplexToModulus() {} 00038 ~ComplexToModulus() {} 00039 bool operator!=( const ComplexToModulus & ) const 00040 { 00041 return false; 00042 } 00043 bool operator==( const ComplexToModulus & other ) const 00044 { 00045 return !(*this != other); 00046 } 00047 inline TOutput operator()( const TInput & A ) 00048 { return (TOutput)( vcl_sqrt(A.real() * A.real() + 00049 A.imag() * A.imag() ) ); } 00050 }; 00051 } 00052 00053 template <class TInputImage, class TOutputImage> 00054 class ITK_EXPORT ComplexToModulusImageFilter : 00055 public 00056 UnaryFunctorImageFilter<TInputImage,TOutputImage, 00057 Function::ComplexToModulus< 00058 typename TInputImage::PixelType, 00059 typename TOutputImage::PixelType> > 00060 { 00061 public: 00063 typedef ComplexToModulusImageFilter Self; 00064 typedef UnaryFunctorImageFilter<TInputImage,TOutputImage, 00065 Function::ComplexToModulus< typename TInputImage::PixelType, 00066 typename TOutputImage::PixelType> > Superclass; 00067 typedef SmartPointer<Self> Pointer; 00068 typedef SmartPointer<const Self> ConstPointer; 00069 00071 itkNewMacro(Self); 00072 00073 typedef typename TInputImage::PixelType InputPixelType; 00074 typedef typename TOutputImage::PixelType OutputPixelType; 00075 typedef typename NumericTraits< InputPixelType >::ValueType InputPixelValueType; 00076 00077 00078 #ifdef ITK_USE_CONCEPT_CHECKING 00079 00080 itkConceptMacro(InputMultiplyOperatorCheck, 00081 (Concept::MultiplyOperator<InputPixelValueType>)); 00082 00084 #endif 00085 00086 protected: 00087 ComplexToModulusImageFilter() {} 00088 virtual ~ComplexToModulusImageFilter() {} 00089 00090 private: 00091 ComplexToModulusImageFilter(const Self&); //purposely not implemented 00092 void operator=(const Self&); //purposely not implemented 00093 00094 }; 00095 00096 } // end namespace itk 00097 00098 00099 #endif 00100