00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017 #ifndef __itkLaplacianRecursiveGaussianImageFilter_h
00018 #define __itkLaplacianRecursiveGaussianImageFilter_h
00019
00020 #include "itkRecursiveGaussianImageFilter.h"
00021 #include "itkNthElementImageAdaptor.h"
00022 #include "itkImage.h"
00023 #include "itkPixelTraits.h"
00024 #include "itkCommand.h"
00025
00026
00027 namespace itk
00028 {
00029
00039
00040
00041
00042 template <typename TInputImage,
00043 typename TOutputImage= TInputImage >
00044 class ITK_EXPORT LaplacianRecursiveGaussianImageFilter:
00045 public ImageToImageFilter<TInputImage,TOutputImage>
00046 {
00047 public:
00049 typedef LaplacianRecursiveGaussianImageFilter 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 typedef MemberCommand< Self > CommandType;
00106 typedef typename CommandType::Pointer CommandPointer;
00107
00109 itkNewMacro(Self);
00110
00112 void SetSigma( RealType sigma );
00113
00115 void SetNormalizeAcrossScale( bool normalizeInScaleSpace );
00116 itkGetMacro( NormalizeAcrossScale, bool );
00118
00124 virtual void GenerateInputRequestedRegion() throw(InvalidRequestedRegionError);
00125
00126 protected:
00127 LaplacianRecursiveGaussianImageFilter();
00128 virtual ~LaplacianRecursiveGaussianImageFilter() {};
00129 void PrintSelf(std::ostream& os, Indent indent) const;
00130
00132 void GenerateData( void );
00133
00134
00135 void EnlargeOutputRequestedRegion(DataObject *output);
00136
00138 void ReportProgress(const Object * object, const EventObject & event );
00139
00140 private:
00141 LaplacianRecursiveGaussianImageFilter(const Self&);
00142 void operator=(const Self&);
00143
00144 GaussianFilterPointer m_SmoothingFilters[ImageDimension-1];
00145 DerivativeFilterPointer m_DerivativeFilter;
00146
00147 CumulativeImagePointer m_CumulativeImage;
00148
00149 CommandPointer m_ProgressCommand;
00150 float m_Progress;
00151
00153 bool m_NormalizeAcrossScale;
00154
00155
00156 };
00157
00158 }
00159
00160 #ifndef ITK_MANUAL_INSTANTIATION
00161 #include "itkLaplacianRecursiveGaussianImageFilter.txx"
00162 #endif
00163
00164 #endif
00165
00166
00167
00168
00169