00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
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 bool operator!=( const Sqrt & ) const
00040 {
00041 return false;
00042 }
00043 bool operator==( const Sqrt & other ) const
00044 {
00045 return !(*this != other);
00046 }
00047 inline TOutput operator()( const TInput & A )
00048 {
00049 return (TOutput)vcl_sqrt((double)A);
00050 }
00051 };
00052 }
00053 template <class TInputImage, class TOutputImage>
00054 class ITK_EXPORT SqrtImageFilter :
00055 public
00056 UnaryFunctorImageFilter<TInputImage,TOutputImage,
00057 Function::Sqrt< typename TInputImage::PixelType,
00058 typename TOutputImage::PixelType> >
00059 {
00060 public:
00062 typedef SqrtImageFilter Self;
00063 typedef UnaryFunctorImageFilter<TInputImage,TOutputImage,
00064 Function::Sqrt< typename TInputImage::PixelType,
00065 typename TOutputImage::PixelType> > Superclass;
00066 typedef SmartPointer<Self> Pointer;
00067 typedef SmartPointer<const Self> ConstPointer;
00068
00070 itkNewMacro(Self);
00071
00073 itkTypeMacro(SqrtImageFilter,
00074 UnaryFunctorImageFilter);
00075
00076 #ifdef ITK_USE_CONCEPT_CHECKING
00077
00078 itkConceptMacro(InputConvertibleToDoubleCheck,
00079 (Concept::Convertible<typename TInputImage::PixelType, double>));
00080 itkConceptMacro(DoubleConvertibleToOutputCheck,
00081 (Concept::Convertible<double, typename TOutputImage::PixelType>));
00082
00084 #endif
00085
00086 protected:
00087 SqrtImageFilter() {}
00088 virtual ~SqrtImageFilter() {}
00089
00090 private:
00091 SqrtImageFilter(const Self&);
00092 void operator=(const Self&);
00093
00094 };
00095
00096 }
00097
00098
00099 #endif
00100