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