00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017 #ifndef __itkGradientMagnitudeRecursiveGaussianImageFilter_h
00018 #define __itkGradientMagnitudeRecursiveGaussianImageFilter_h
00019
00020 #include "itkNthElementImageAdaptor.h"
00021 #include "itkImage.h"
00022 #include "itkPixelTraits.h"
00023 #include "itkRecursiveGaussianImageFilter.h"
00024
00025 namespace itk
00026 {
00027
00039
00040
00041
00042 template <typename TInputImage,
00043 typename TOutputImage= TInputImage >
00044 class ITK_EXPORT GradientMagnitudeRecursiveGaussianImageFilter:
00045 public ImageToImageFilter<TInputImage,TOutputImage>
00046 {
00047 public:
00049 typedef GradientMagnitudeRecursiveGaussianImageFilter Self;
00050 typedef ImageToImageFilter<TInputImage,TOutputImage> Superclass;
00051 typedef SmartPointer<Self> Pointer;
00052 typedef SmartPointer<const Self> ConstPointer;
00053
00054
00056 typedef TInputImage InputImageType;
00057 typedef typename InputImageType::PixelType PixelType;
00058
00059
00061 itkStaticConstMacro(ImageDimension, unsigned int,
00062 TInputImage::ImageDimension);
00063
00064 typedef typename NumericTraits<PixelType>::RealType RealType;
00065
00070 typedef float InternalRealType;
00071 typedef Image<InternalRealType,
00072 itkGetStaticConstMacro(ImageDimension) > RealImageType;
00073
00075 typedef RecursiveGaussianImageFilter<
00076 RealImageType,
00077 RealImageType
00078 > GaussianFilterType;
00079
00081 typedef RecursiveGaussianImageFilter<
00082 InputImageType,
00083 RealImageType
00084 > DerivativeFilterType;
00085
00087 typedef typename GaussianFilterType::Pointer GaussianFilterPointer;
00088
00090 typedef typename DerivativeFilterType::Pointer DerivativeFilterPointer;
00091
00093 typedef typename TOutputImage::Pointer OutputImagePointer;
00094
00096 typedef TOutputImage OutputImageType;
00097 typedef typename OutputImageType::PixelType OutputPixelType;
00098
00100 typedef Image< InternalRealType,
00101 itkGetStaticConstMacro(ImageDimension) > CumulativeImageType;
00102 typedef typename CumulativeImageType::Pointer CumulativeImagePointer;
00103
00105 itkNewMacro(Self);
00106
00108 void SetSigma( RealType sigma );
00109
00111 void SetNormalizeAcrossScale( bool normalizeInScaleSpace );
00112 itkGetMacro( NormalizeAcrossScale, bool );
00113
00119 virtual void GenerateInputRequestedRegion() throw(InvalidRequestedRegionError);
00120
00121 protected:
00122 GradientMagnitudeRecursiveGaussianImageFilter();
00123 virtual ~GradientMagnitudeRecursiveGaussianImageFilter() {};
00124 void PrintSelf(std::ostream& os, Indent indent) const;
00125
00127 void GenerateData( void );
00128
00129
00130 void EnlargeOutputRequestedRegion(DataObject *output);
00131
00132 private:
00133 GradientMagnitudeRecursiveGaussianImageFilter(const Self&);
00134 void operator=(const Self&);
00135
00136 GaussianFilterPointer m_SmoothingFilters[ImageDimension-1];
00137 DerivativeFilterPointer m_DerivativeFilter;
00138
00139 CumulativeImagePointer m_CumulativeImage;
00140
00142 bool m_NormalizeAcrossScale;
00143
00144
00145 };
00146
00147 }
00148
00149 #ifndef ITK_MANUAL_INSTANTIATION
00150 #include "itkGradientMagnitudeRecursiveGaussianImageFilter.txx"
00151 #endif
00152
00153 #endif
00154
00155
00156
00157