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
00072 #ifdef ITK_USE_CONCEPT_CHECKING
00073
00074 itkConceptMacro(InputConvertibleToDoubleCheck,
00075 (Concept::Convertible<typename TInputImage::PixelType, double>));
00076 itkConceptMacro(DoubleConvertibleToOutputCheck,
00077 (Concept::Convertible<double, typename TOutputImage::PixelType>));
00078
00080 #endif
00081
00082 protected:
00083 SqrtImageFilter() {}
00084 virtual ~SqrtImageFilter() {}
00085
00086 private:
00087 SqrtImageFilter(const Self&);
00088 void operator=(const Self&);
00089
00090 };
00091
00092 }
00093
00094
00095 #endif
00096