00001 
00002 
00003 
00004 
00005 
00006 
00007 
00008 
00009 
00010 
00011 
00012 
00013 
00014 
00015 
00016 
00017 
00018 
00019 
00020 
00021 
00022 
00023 
00024 
00025 
00026 
00027 
00028 
00029 
00030 
00031 
00032 
00033 
00034 
00035 
00036 
00037 
00038 
00054 #ifndef __vtkImageGaussianSmooth_h
00055 #define __vtkImageGaussianSmooth_h
00056 
00057 
00058 #include "vtkImageToImageFilter.h"
00059 
00060 class VTK_EXPORT vtkImageGaussianSmooth : public vtkImageToImageFilter
00061 {
00062 public:
00063   vtkTypeMacro(vtkImageGaussianSmooth,vtkImageToImageFilter);
00064   void PrintSelf(ostream& os, vtkIndent indent);
00065 
00069   static vtkImageGaussianSmooth *New();
00070   
00071   
00073   vtkSetVector3Macro(StandardDeviations, float);
00074   void SetStandardDeviation(float std)
00075         {this->SetStandardDeviations(std,std,std);}
00076   void SetStandardDeviations(float a,float b)
00077         {this->SetStandardDeviations(a,b,0.0);}
00078   vtkGetVector3Macro(StandardDeviations, float);
00079 
00082   void SetStandardDeviation(float a,float b)
00083         {this->SetStandardDeviations(a,b,0.0);}
00084   void SetStandardDeviation(float a,float b,float c) 
00085         {this->SetStandardDeviations(a,b,c);}
00086 
00090   vtkSetVector3Macro(RadiusFactors, float);
00091   void SetRadiusFactors(float f, float f2) {this->SetRadiusFactors(f,f2,1.5);}
00092   void SetRadiusFactor(float f) {this->SetRadiusFactors(f, f, f);}
00093   vtkGetVector3Macro(RadiusFactors, float);
00094 
00097   vtkSetMacro(Dimensionality, int);
00098   vtkGetMacro(Dimensionality, int);
00099 
00100 protected:
00101   vtkImageGaussianSmooth();
00102   ~vtkImageGaussianSmooth();
00103   vtkImageGaussianSmooth(const vtkImageGaussianSmooth&) {};
00104   void operator=(const vtkImageGaussianSmooth&) {};
00105 
00106   int Dimensionality;
00107   float StandardDeviations[3];
00108   float RadiusFactors[3];
00109   
00110   void ComputeKernel(double *kernel, int min, int max, double std);
00111   void ComputeInputUpdateExtent(int inExt[6], int outExt[6]);
00112   void ExecuteAxis(int axis, vtkImageData *inData, int inExt[6],
00113          vtkImageData *outData, int outExt[6],
00114          int *pcycle, int target, int *pcount, int total);
00115   void ThreadedExecute(vtkImageData *inData, 
00116              vtkImageData *outData, int outExt[6], int id);
00117   
00118 };
00119 
00120 #endif
00121 
00122 
00123 
00124 
00125 
00126 
00127 
00128 
00129 
00130