SRT2DReconstruction.h

Go to the documentation of this file.

#ifndef __stir_analytic_SRT2D_SRT2DReconstruction_H__
#define __stir_analytic_SRT2D_SRT2DReconstruction_H__
/*
    Copyright (C) 2012-2016, 2023-2024, Dimitra Kyriakopoulou
    Copyright (C) 2024, University College London
    This file is part of STIR.
 
    SPDX-License-Identifier: Apache-2.0
 
    See STIR/LICENSE.txt for details
*/
 
#include "stir/recon_buildblock/AnalyticReconstruction.h"
#include "stir/RegisteredParsingObject.h"
#include <string>
#include <vector>
#include "stir/shared_ptr.h"
 
START_NAMESPACE_STIR
 
template <int num_dimensions, typename elemT>
class DiscretisedDensity;
class Succeeded;
class ProjData;

 
class SRT2DReconstruction
    : public RegisteredParsingObject<SRT2DReconstruction, Reconstruction<DiscretisedDensity<3, float>>, AnalyticReconstruction>
{
  // typedef AnalyticReconstruction base_type;
  typedef RegisteredParsingObject<SRT2DReconstruction, Reconstruction<DiscretisedDensity<3, float>>, AnalyticReconstruction>
      base_type;
#ifdef SWIG
  // work-around swig problem. It gets confused when using a private (or protected)
  // typedef in a definition of a public typedef/member
 public:
#else
private:
#endif
  typedef DiscretisedDensity<3, float> TargetT;
 
public:
  static constexpr const char* const registered_name = "SRT2D";

  SRT2DReconstruction();
  explicit SRT2DReconstruction(const std::string& parameter_filename);
 
  SRT2DReconstruction(const shared_ptr<ProjData>& proj_data_ptr_v, const int num_segments_to_combine = -1);
  virtual std::string method_info() const;
 
  virtual void ask_parameters();
 
  virtual Succeeded set_up(shared_ptr<TargetT> const& target_data_sptr);
 
protected: // make parameters protected such that doc shows always up in doxygen
  // parameters used for parsing


  int num_segments_to_combine;
 
private:
  Succeeded actual_reconstruct(shared_ptr<DiscretisedDensity<3, float>> const& target_image_ptr);
 
  virtual void set_defaults();
  virtual void initialise_keymap();
  virtual bool post_processing();

  void spline(const std::vector<float>& x, const std::vector<float>& y, int n, std::vector<float>& y2) const;

  float hilbert_der(float x,
                    const std::vector<float>& f,
                    const std::vector<float>& ddf,
                    const std::vector<float>& p,
                    int sp,
                    const std::vector<float>& lg,
                    float termC) const;

  float integ(float dist, int max, const std::vector<float>& ff) const;
};
 
END_NAMESPACE_STIR
 
#endif