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__