NonRigidObjectTransformationUsingBSplines.h

Go to the documentation of this file.

//
//
/*
    Copyright (C) 2005- 2009 , Hammersmith Imanet Ltd
    For internal GE use only
*/
#ifndef __stir_motion_NonRigidObjectTransformationUsingBSplines_H__
#define __stir_motion_NonRigidObjectTransformationUsingBSplines_H__

#include "stir_experimental/motion/ObjectTransformation.h"
#include "stir/CartesianCoordinate3D.h"
#include "stir/numerics/BSplines.h"
#include "stir/numerics/BSplinesRegularGrid.h"
#include "stir/RegisteredParsingObject.h"
#include "stir/shared_ptr.h"
//#include "stir/DiscretisedDensityOnCartesianGrid.h"

START_NAMESPACE_STIR
class Succeeded;

template <int num_dimensions, class elemT>
class DeformationFieldOnCartesianGrid : public BasicCoordinate<num_dimensions, Array<num_dimensions, elemT>>
//  public DiscretisedDensityOnCartesianGrid<num_dimensions, BasicCoordinate<num_dimensions, elemT> >
{
public:
  DeformationFieldOnCartesianGrid() {}
};

template <int num_dimensions, class elemT>
class NonRigidObjectTransformationUsingBSplines
    : public RegisteredParsingObject<NonRigidObjectTransformationUsingBSplines<num_dimensions, elemT>,
                                     ObjectTransformation<num_dimensions, elemT>,
                                     ObjectTransformation<num_dimensions, elemT>>
{
public:
  static const char* const registered_name;

  // Default constructor
  NonRigidObjectTransformationUsingBSplines();

  NonRigidObjectTransformationUsingBSplines(const std::string& filename_x,
                                            const std::string& filename_y,
                                            const std::string& filename_z,
                                            const int bspline_order);

  NonRigidObjectTransformationUsingBSplines(const std::string& filename, const int bspline_order);

  virtual BasicCoordinate<num_dimensions, elemT> transform_point(const BasicCoordinate<num_dimensions, elemT>& point) const;

  float jacobian(const BasicCoordinate<num_dimensions, elemT>& point) const;

  virtual void set_defaults();
  virtual void initialise_keymap();
  virtual bool post_processing();
  virtual void set_key_values();

private:
  BasicCoordinate<num_dimensions, BSpline::BSplinesRegularGrid<num_dimensions, elemT, elemT>> interpolator;
  // BasicCoordinate<num_dimensions,elemT> _grid_spacing;
  // BasicCoordinate<num_dimensions,elemT> _origin;
  CartesianCoordinate3D<elemT> _grid_spacing;
  CartesianCoordinate3D<elemT> _origin;
  BSpline::BSplineType _bspline_type;

  // use for parsing only
  shared_ptr<DeformationFieldOnCartesianGrid<num_dimensions, elemT>> deformation_field_sptr;
  int _bspline_order;
  // for NCAT only
  std::string _deformation_field_from_NCAT_file;
  CartesianCoordinate3D<int> _deformation_field_from_NCAT_size;
  // from file
  std::string _deformation_field_from_file_x;
  std::string _deformation_field_from_file_y;
  std::string _deformation_field_from_file_z;
  // multicomponent filename
  std::string _deformation_field_multicomponent_filename;
};
#if 0

std::ostream&
operator<<(std::ostream& out,
           const NonRigidObjectTransformationUsingBSplines& rigid_object_transformation);

std::istream&
operator>>(std::istream& ,
           NonRigidObjectTransformationUsingBSplines& rigid_object_transformation);
#endif

END_NAMESPACE_STIR

#endif