00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017 #ifndef __itkTernaryMagnitudeImageFilter_h
00018 #define __itkTernaryMagnitudeImageFilter_h
00019
00020 #include "itkTernaryFunctorImageFilter.h"
00021
00022 namespace itk
00023 {
00024
00034 namespace Function {
00035
00036 template< class TInput1, class TInput2, class TInput3, class TOutput>
00037 class Modulus3
00038 {
00039 public:
00040 Modulus3() {}
00041 ~Modulus3() {}
00042 bool operator!=( const Modulus3 & ) const
00043 {
00044 return false;
00045 }
00046 bool operator==( const Modulus3 & other ) const
00047 {
00048 return !(*this != other);
00049 }
00050 inline TOutput operator()( const TInput1 & A,
00051 const TInput2 & B,
00052 const TInput3 & C)
00053 { return (TOutput) vcl_sqrt((double)(A*A + B*B + C*C) ); }
00054 };
00055 }
00056
00057 template <class TInputImage1, class TInputImage2,
00058 class TInputImage3, class TOutputImage>
00059 class ITK_EXPORT TernaryMagnitudeImageFilter :
00060 public
00061 TernaryFunctorImageFilter<TInputImage1,TInputImage2,
00062 TInputImage3,TOutputImage,
00063 Function::Modulus3<
00064 typename TInputImage1::PixelType,
00065 typename TInputImage2::PixelType,
00066 typename TInputImage3::PixelType,
00067 typename TOutputImage::PixelType> >
00068 {
00069 public:
00071 typedef TernaryMagnitudeImageFilter Self;
00072 typedef TernaryFunctorImageFilter<TInputImage1,TInputImage2,
00073 TInputImage3,TOutputImage,
00074 Function::Modulus3<
00075 typename TInputImage1::PixelType,
00076 typename TInputImage2::PixelType,
00077 typename TInputImage3::PixelType,
00078 typename TOutputImage::PixelType>
00079 > Superclass;
00080 typedef SmartPointer<Self> Pointer;
00081 typedef SmartPointer<const Self> ConstPointer;
00082
00084 itkNewMacro(Self);
00085
00087 itkTypeMacro(TernaryMagnitudeImageFilter,
00088 TernaryFunctorImageFilter);
00089
00090 protected:
00091 TernaryMagnitudeImageFilter() {}
00092 virtual ~TernaryMagnitudeImageFilter() {}
00093
00094 private:
00095 TernaryMagnitudeImageFilter(const Self&);
00096 void operator=(const Self&);
00097
00098
00099 };
00100
00101 }
00102
00103
00104 #endif
00105