ITK  4.1.0
Insight Segmentation and Registration Toolkit
itkSqrtImageFilter.h
Go to the documentation of this file.
00001 /*=========================================================================
00002  *
00003  *  Copyright Insight Software Consortium
00004  *
00005  *  Licensed under the Apache License, Version 2.0 (the "License");
00006  *  you may not use this file except in compliance with the License.
00007  *  You may obtain a copy of the License at
00008  *
00009  *         http://www.apache.org/licenses/LICENSE-2.0.txt
00010  *
00011  *  Unless required by applicable law or agreed to in writing, software
00012  *  distributed under the License is distributed on an "AS IS" BASIS,
00013  *  WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
00014  *  See the License for the specific language governing permissions and
00015  *  limitations under the License.
00016  *
00017  *=========================================================================*/
00018 #ifndef __itkSqrtImageFilter_h
00019 #define __itkSqrtImageFilter_h
00020 
00021 #include "itkUnaryFunctorImageFilter.h"
00022 #include "vnl/vnl_math.h"
00023 
00024 namespace itk
00025 {
00026 namespace Functor
00027 {
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   }
00044 
00045   bool operator==(const Sqrt & other) const
00046   {
00047     return !( *this != other );
00048   }
00049 
00050   inline TOutput operator()(const TInput & A) const
00051   {
00052     return (TOutput)vcl_sqrt( (double)A );
00053   }
00054 };
00055 }
00065 template< class TInputImage, class TOutputImage >
00066 class ITK_EXPORT SqrtImageFilter:
00067   public
00068   UnaryFunctorImageFilter< TInputImage, TOutputImage,
00069                            Functor::Sqrt< typename TInputImage::PixelType,
00070                                            typename TOutputImage::PixelType >   >
00071 {
00072 public:
00074   typedef SqrtImageFilter Self;
00075   typedef UnaryFunctorImageFilter<
00076     TInputImage, TOutputImage,
00077     Functor::Sqrt< typename TInputImage::PixelType,
00078                    typename TOutputImage::PixelType > >  Superclass;
00079 
00080   typedef SmartPointer< Self >       Pointer;
00081   typedef SmartPointer< const Self > ConstPointer;
00082 
00084   itkNewMacro(Self);
00085 
00087   itkTypeMacro(SqrtImageFilter,
00088                UnaryFunctorImageFilter);
00089 
00090 #ifdef ITK_USE_CONCEPT_CHECKING
00091 
00092   itkConceptMacro( InputConvertibleToDoubleCheck,
00093                    ( Concept::Convertible< typename TInputImage::PixelType, double > ) );
00094   itkConceptMacro( DoubleConvertibleToOutputCheck,
00095                    ( Concept::Convertible< double, typename TOutputImage::PixelType > ) );
00096 
00098 #endif
00099 protected:
00100   SqrtImageFilter() {}
00101   virtual ~SqrtImageFilter() {}
00102 private:
00103   SqrtImageFilter(const Self &); //purposely not implemented
00104   void operator=(const Self &);  //purposely not implemented
00105 };
00106 } // end namespace itk
00108 
00109 #endif
00110