vtkBSplineInterpolateImageFunction.h

Go to the documentation of this file.

/*=auto=========================================================================

  Portions (c) Copyright 2005 Brigham and Women's Hospital (BWH) All Rights Reserved.

  See Doc/copyright/copyright.txt
  or http://www.slicer.org/copyright/copyright.txt for details.

  Program:   3D Slicer
  Module:    $RCSfile: vtkBSplineInterpolateImageFunction.h,v $
  Date:      $Date: 2006/05/26 20:02:42 $
  Version:   $Revision: 1.10 $

=========================================================================auto=*/


#ifndef __vtkBSplineInterpolateImageFunction_h
#define __vtkBSplineInterpolateImageFunction_h

#include "vtkImplicitFunction.h"
#include "vtkTeemConfigure.h"

#include <vtkstd/vector>  

#define VTK_INTEGRATE_MAJOR_EIGENVECTOR  0
#define VTK_INTEGRATE_MEDIUM_EIGENVECTOR 1
#define VTK_INTEGRATE_MINOR_EIGENVECTOR  2

#define ImageDimension 3

class vtkImageData;
class  VTK_Teem_EXPORT vtkBSplineInterpolateImageFunction : public vtkImplicitFunction
{
 public:
  static vtkBSplineInterpolateImageFunction *New();
  vtkTypeRevisionMacro(vtkBSplineInterpolateImageFunction, vtkImplicitFunction );
  virtual void PrintSelf(ostream& os, vtkIndent indent);

  virtual vtkFloatingPointType EvaluateFunction (vtkFloatingPointType x[ImageDimension]);

  virtual void EvaluateGradient (vtkFloatingPointType x[ImageDimension], 
                                 vtkFloatingPointType g[ImageDimension]);
  void SetInput(vtkImageData* dataset);

  unsigned int GetSplineOrder() { return this->SplineOrder; }
  void SetSplineOrder(unsigned int order);

protected:
  vtkBSplineInterpolateImageFunction() {
    this->Initialized = 0;
    this->SplineOrder = 0; 
    this->SetSplineOrder(3);
  }
  ~vtkBSplineInterpolateImageFunction() {}
  int                                 DataLength[ImageDimension];  
  unsigned int                        SplineOrder;    

  vtkImageData *                      Coefficients; 
  vtkFloatingPointType * Origin;
  vtkFloatingPointType * Spacing;
  int * Extent;
private:
  //BTX
  int Initialized;
  vtkBSplineInterpolateImageFunction(const vtkBSplineInterpolateImageFunction&);  
  void operator=(const vtkBSplineInterpolateImageFunction&);  

  void SetInterpolationWeights( vtkFloatingPointType *x, long *evaluateIndex[ImageDimension],
    double *weights[ImageDimension],unsigned int splineOrder ) const;

  void SetDerivativeWeights( vtkFloatingPointType *x, long *evaluateIndex[ImageDimension],
    double *weights[ImageDimension], unsigned int splineOrder ) const;

  void GeneratePointsToIndex();

  void DetermineRegionOfSupport( long *evaluateIndex[ImageDimension],
    vtkFloatingPointType x[], unsigned int splineOrder ) const;

  void ApplyMirrorBoundaryConditions(long *evaluateIndex[ImageDimension], 
    unsigned int splineOrder) const;


  unsigned int             MaxNumberInterpolationPoints; 
  std::vector<int>    PointsToIndex[ImageDimension];  
  //ETX
};


#endif