OSSPSReconstruction.h

Go to the documentation of this file.

 
//
//
/*
    Copyright (C) 2002- 2009, Hammersmith Imanet Ltd
    This file is part of STIR.
 
    SPDX-License-Identifier: Apache-2.0
 
    See STIR/LICENSE.txt for details
*/
 
#ifndef __stir_OSSPS_OSSPSReconstruction_h__
#  define __stir_OSSPS_OSSPSReconstruction_h__
 
#  include "stir/recon_buildblock/IterativeReconstruction.h"
#  include "stir/recon_buildblock/BinNormalisation.h"
#  include "stir/ProjData.h"
#  include "stir/RegisteredParsingObject.h"
 
START_NAMESPACE_STIR

template <class TargetT>
class OSSPSReconstruction
    : public RegisteredParsingObject<OSSPSReconstruction<TargetT>, Reconstruction<TargetT>, IterativeReconstruction<TargetT>>
// public IterativeReconstruction<TargetT>
{
private:
  typedef RegisteredParsingObject<OSSPSReconstruction<TargetT>, Reconstruction<TargetT>, IterativeReconstruction<TargetT>>
      base_type;
 
public:
  static const char* const registered_name;

  OSSPSReconstruction();

  explicit OSSPSReconstruction(const std::string& parameter_filename);

  void ask_parameters() override;
  OSSPSReconstruction& get_parameters() { return *this; }

  const OSSPSReconstruction& get_parameters() const { return *this; }

  std::string method_info() const override;


  Succeeded precompute_denominator_of_conditioner_without_penalty();

  Succeeded set_up(shared_ptr<TargetT> const& target_image_ptr) override;

  void update_estimate(TargetT& current_image_estimate) override;
 
protected: // could be private, but this way the doxygen comments are always listed
  int enforce_initial_positivity;

  double upper_bound;

  int write_update_image;


  std::string precomputed_denominator_filename;
 
#  if 0
  bool do_line_search;
#  endif

  float relaxation_parameter;
  float relaxation_gamma;
 
  void set_defaults() override;
  void initialise_keymap() override;
  bool post_processing() override;
 
private:
  shared_ptr<TargetT> precomputed_denominator_ptr;


  shared_ptr<ProjData> fwd_ones_sptr;
 
  GeneralisedPrior<TargetT>* get_prior_ptr() { return this->get_objective_function().get_prior_ptr(); }
 
#  if 0
  float
    line_search(const TargetT& current_estimate, const TargetT& additive_update);
#  endif
};
 
END_NAMESPACE_STIR
 
#endif
 
// __OSSPSReconstruction_h__