SRT2DSPECTReconstruction.h

Go to the documentation of this file.

//
//
#ifndef __stir_analytic_SRT2DSPECT_SRT2DSPECTReconstruction_H__
#define __stir_analytic_SRT2DSPECT_SRT2DSPECTReconstruction_H__
/*
    Copyright (C) 2014-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 "stir/shared_ptr.h"
#include <string>
#include <vector>
 
START_NAMESPACE_STIR
 
template <int num_dimensions, typename elemT>
class DiscretisedDensity;
class Succeeded;
class ProjData;

 
class SRT2DSPECTReconstruction : public RegisteredParsingObject<SRT2DSPECTReconstruction,
                                                                Reconstruction<DiscretisedDensity<3, float>>,
                                                                AnalyticReconstruction>
{
  // typedef AnalyticReconstruction base_type;
  typedef RegisteredParsingObject<SRT2DSPECTReconstruction, 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 = "SRT2DSPECT";

  SRT2DSPECTReconstruction();
  explicit SRT2DSPECTReconstruction(const std::string& parameter_filename);
 
  SRT2DSPECTReconstruction(const shared_ptr<ProjData>& proj_data_ptr_v);
 
  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


  std::string attenuation_projection_filename;
  float thres_restr_bound;
  std::vector<double> thres_restr_bound_vector;
  shared_ptr<ProjData> atten_data_ptr;
 
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();

  float hilbert_node(
      float x, const std::vector<float>& f, const std::vector<float>& ddf, const std::vector<float>& p, int sp, float fn) const;

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

  void hilbert_der_double(float x,
                          const std::vector<float>& f,
                          const std::vector<float>& ddf,
                          const std::vector<float>& f1,
                          const std::vector<float>& ddf1,
                          const std::vector<float>& p,
                          int sp,
                          float* dhp,
                          float* dh1p,
                          const std::vector<float>& lg) const;

  float splint(const std::vector<float>& xa, const std::vector<float>& ya, const std::vector<float>& y2a, int n, float x) const;

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

  float integ(float dist, int max, float ff[]) const;
};
 
END_NAMESPACE_STIR
 
#endif