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