transform_3d_object.h

Go to the documentation of this file.

//
//
/*
    Copyright (C) 2003- 2011, Hammersmith Imanet Ltd
    For internal GE use only
*/
#ifndef __stir_motion_transformed_3d_object_H__
#define __stir_motion_transformed_3d_object_H__

#include "stir/common.h"

// for caching functions
#include <utility>
#include "stir/Array.h"
#include "stir/assign.h"
#include "stir/BasicCoordinate.h"

START_NAMESPACE_STIR

class RigidObject3DTransformation;
template <int num_dimensions, class elemT>
class ObjectTransformation;

template <int num_dimensions, typename elemT>
class DiscretisedDensity;
class ProjData;
class Succeeded;


template <class ObjectTransformationT, class PushInterpolatorT>
inline Succeeded transform_3d_object_push_interpolation(DiscretisedDensity<3, float>& out_density,
                                                        const DiscretisedDensity<3, float>& in_density,
                                                        const ObjectTransformationT& transformation_in_to_out,
                                                        const PushInterpolatorT& interpolator,
                                                        const bool do_jacobian);


template <class ObjectTransformationT, class PullInterpolatorT>
inline Succeeded transform_3d_object_pull_interpolation(DiscretisedDensity<3, float>& out_density,
                                                        const DiscretisedDensity<3, float>& in_density,
                                                        const ObjectTransformationT& transformation_out_to_in,
                                                        const PullInterpolatorT& interpolator,
                                                        const bool do_jacobian);


Succeeded transform_3d_object(DiscretisedDensity<3, float>& out_density,
                              const DiscretisedDensity<3, float>& in_density,
                              const RigidObject3DTransformation& transformation_in_to_out);
//                  const ObjectTransformation<3,float>& transformation_in_to_out);

// ugly functions for storing transformed points.
// TODO clean up at some point

Array<3, BasicCoordinate<3, float>>
find_grid_coords_of_transformed_centres(const DiscretisedDensity<3, float>& source_density,
                                        const DiscretisedDensity<3, float>& target_density,
                                        const ObjectTransformation<3, float>& transformation_source_to_target);

// TODO need this for now to get Array<std::pair<>> to work
template <class T2>
inline void
assign(std::pair<BasicCoordinate<3, float>, float>& x, T2 y)
{
  BasicCoordinate<3, float> tmp;
  assign(tmp, y);
  x = std::make_pair(tmp, float(y));
}

Array<3, std::pair<BasicCoordinate<3, float>, float>>
find_grid_coords_of_transformed_centres_and_jacobian(const DiscretisedDensity<3, float>& source_density,
                                                     const DiscretisedDensity<3, float>& target_density,
                                                     const ObjectTransformation<3, float>& transformation_source_to_target);


Succeeded transform_3d_object(ProjData& out_proj_data,
                              const ProjData& in_proj_data,
                              const RigidObject3DTransformation& object_transformation);


Succeeded transpose_of_transform_3d_object(DiscretisedDensity<3, float>& out_density,
                                           const DiscretisedDensity<3, float>& in_density,
                                           const RigidObject3DTransformation& transformation_in_to_out);
//                  const ObjectTransformation<3,float>& transformation_in_to_out);


Succeeded transform_3d_object(ProjData& out_proj_data,
                              const ProjData& in_proj_data,
                              const RigidObject3DTransformation& rigid_object_transformation,
                              const int min_in_segment_num_to_process,
                              const int max_in_segment_num_to_process);

END_NAMESPACE_STIR

#include "stir_experimental/motion/transform_3d_object.inl"

#endif