/vol/vipdata/irtk/packages/transformation/include/irtkLinearFreeFormTransformation.h

/*=========================================================================

  Library   : Image Registration Toolkit (IRTK)
  Module    : $Id: irtkLinearFreeFormTransformation.h 2 2008-12-23 12:40:14Z dr $
  Copyright : Imperial College, Department of Computing
              Visual Information Processing (VIP), 2008 onwards
  Date      : $Date: 2008-12-23 12:40:14 +0000 (Tue, 23 Dec 2008) $
  Version   : $Revision: 2 $
  Changes   : $Author: dr $

=========================================================================*/

#ifndef _IRTKLINEARFREEFORMTRANSFORMATION_H

#define _IRTKLINEARFREEFORMTRANSFORMATION_H

#include <irtkGeometry.h>

class irtkLinearFreeFormTransformation : public irtkFreeFormTransformation3D
{

public:

  irtkLinearFreeFormTransformation();

  irtkLinearFreeFormTransformation(irtkBaseImage &, double = 1, double = 1, double = 1);

  irtkLinearFreeFormTransformation(double x1, double y1, double z1,
                                   double x2, double y2, double z2,
                                   double dx, double dy, double dz,
                                   double* xaxis, double* yaxis, double* zaxis);

  irtkLinearFreeFormTransformation(const class irtkLinearFreeFormTransformation &);

  irtkLinearFreeFormTransformation(const class irtkBSplineFreeFormTransformation3D &);

  virtual ~irtkLinearFreeFormTransformation();

  virtual double Approximate(double *, double *, double *,
                             double *, double *, double *, int);

  virtual void Interpolate(double* dxs, double* dys, double* dzs);

  virtual void Subdivide();

  virtual void FFD1(double &, double &, double &) const;

  virtual void FFD2(double &, double &, double &) const;

  virtual void Transform(double &, double &, double &, double = 0);

  virtual void Transform2(double &, double &, double &, double = 0);

  virtual void GlobalTransform(double &, double &, double &, double = 0);

  virtual void LocalTransform (double &, double &, double &, double = 0);

  virtual void GlobalDisplacement(double &, double &, double &, double = 0);

  virtual void LocalDisplacement(double &, double &, double &, double = 0);

  virtual void Jacobian(irtkMatrix &, double, double, double, double = 0);

  virtual void LocalJacobian(irtkMatrix &, double, double, double, double = 0);

  virtual void GlobalJacobian(irtkMatrix &, double, double, double, double = 0);

  virtual double Bending(double x, double y, double z);

  virtual void BoundingBox(int, irtkPoint &, irtkPoint &, double = 1) const;

  virtual void BoundingBox(int, double &, double &, double &,
                           double &, double &, double &, double = 1) const;

  virtual void BoundingBox(irtkGreyImage *, int, int &, int &, int &,
                           int &, int &, int &, double = 1) const;

  virtual void Print();

  static int CheckHeader(char *);

  virtual const char *NameOfClass();

  virtual irtkCifstream& Read(irtkCifstream&);

  virtual irtkCofstream& Write(irtkCofstream&);

  virtual istream& Import(istream&);

  virtual ostream& Export(ostream&);

};

inline void irtkLinearFreeFormTransformation::Transform(double &x, double &y, double &z, double)
{
  double u, v, w;

  u = x;
  v = y;
  w = z;

  // Convert world coordinates in to FFD coordinates
  this->WorldToLattice(u, v, w);

  // Calculate FFD
  this->FFD1(u, v, w);

  // Add FFD to world coordinates
  x += u;
  y += v;
  z += w;

}

inline void irtkLinearFreeFormTransformation::Transform2(double &x, double &y, double &z, double)
{
  double u, v, w;

  u = x;
  v = y;
  w = z;

  // Convert world coordinates in to FFD coordinates
  this->WorldToLattice(u, v, w);

  // Calculate FFD
  this->FFD2(u, v, w);

  // Add FFD to world coordinates
  x += u;
  y += v;
  z += w;
}

inline void irtkLinearFreeFormTransformation::GlobalTransform(double &x, double &y, double &z, double)
{}

inline void irtkLinearFreeFormTransformation::LocalTransform(double &x, double &y, double &z, double t)
{
  double u, v, w;

  u = x;
  v = y;
  w = z;

  // calculate displacement
  this->LocalDisplacement(u, v, w, t);

  // Add displacement
  x += u;
  y += v;
  z += w;
}

inline void irtkLinearFreeFormTransformation::GlobalDisplacement(double &x, double &y, double &z, double)
{
  x = 0;
  y = 0;
  z = 0;
}

inline void irtkLinearFreeFormTransformation::LocalDisplacement(double &x, double &y, double &z, double)
{
  // Convert world coordinates in to FFD coordinates
  this->WorldToLattice(x, y, z);

  // Calculate FFD
  this->FFD1(x, y, z);
}

inline const char *irtkLinearFreeFormTransformation::NameOfClass()
{
  return "irtkLinearFreeFormTransformation";
}

#endif

Generated on Mon Mar 23 16:55:53 2009 for IRTK by Doxygen 1.5.4