00001 /*========================================================================= 00002 00003 Program: Insight Segmentation & Registration Toolkit 00004 Module: $RCSfile: itkSqrtImageFilter.h,v $ 00005 Language: C++ 00006 Date: $Date: 2002/01/15 19:36:14 $ 00007 Version: $Revision: 1.14 $ 00008 00009 Copyright (c) 2002 Insight 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 __itkSqrtImageFilter_h 00018 #define __itkSqrtImageFilter_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 Sqrt 00035 { 00036 public: 00037 Sqrt() {} 00038 ~Sqrt() {} 00039 inline TOutput operator()( const TInput & A ) 00040 { 00041 return (TOutput)sqrt((double)A); 00042 } 00043 }; 00044 } 00045 template <class TInputImage, class TOutputImage> 00046 class ITK_EXPORT SqrtImageFilter : 00047 public 00048 UnaryFunctorImageFilter<TInputImage,TOutputImage, 00049 Function::Sqrt< typename TInputImage::PixelType, 00050 typename TOutputImage::PixelType> > 00051 { 00052 public: 00054 typedef SqrtImageFilter Self; 00055 typedef UnaryFunctorImageFilter<TInputImage,TOutputImage, 00056 Function::Sqrt< typename TInputImage::PixelType, 00057 typename TOutputImage::PixelType> > Superclass; 00058 typedef SmartPointer<Self> Pointer; 00059 typedef SmartPointer<const Self> ConstPointer; 00060 00062 itkNewMacro(Self); 00063 00064 protected: 00065 SqrtImageFilter() {} 00066 virtual ~SqrtImageFilter() {} 00067 00068 private: 00069 SqrtImageFilter(const Self&); //purposely not implemented 00070 void operator=(const Self&); //purposely not implemented 00071 00072 }; 00073 00074 } // end namespace itk 00075 00076 00077 #endif