HUToMuImageProcessor.h

Go to the documentation of this file.

/*
    Copyright (C) 2020, UCL
    See STIR/LICENSE.txt for details
*/

#ifndef __stir_HUToMuImageProcessor_H__
#define __stir_HUToMuImageProcessor_H__

#include "stir/RegisteredParsingObject.h"
#include "stir/DataProcessor.h"
#include "stir/DiscretisedDensity.h"

START_NAMESPACE_STIR

template <typename TargetT>
class HUToMuImageProcessor
    : public RegisteredParsingObject<HUToMuImageProcessor<TargetT>, DataProcessor<TargetT>, DataProcessor<TargetT>>
{
private:
  typedef RegisteredParsingObject<HUToMuImageProcessor<TargetT>, DataProcessor<TargetT>, DataProcessor<TargetT>> base_type;

public:
  static const char* const registered_name;

  HUToMuImageProcessor();

  void set_slope_filename(const std::string& filename);
  void set_manufacturer_name(const std::string& name);
  void set_kilovoltage_peak(const float kVp);
  void set_target_photon_energy(const float gamma_energy);

  void apply_scaling_to_HU(TargetT& output_image, const TargetT& input_image) const;

  void set_slope(float a1, float a2, float b1, float b2, float breakPoint);

protected:
  // parsing functions

  void set_defaults() override;
  void initialise_keymap() override;
  bool post_processing() override;


  Succeeded virtual_set_up(const TargetT& image) override;

  void virtual_apply(TargetT& out_density, const TargetT& in_density) const override;
  void virtual_apply(TargetT& density) const override;

private:
  std::string filename;
  std::string manufacturer_name;
  float kilovoltage_peak;
  float target_photon_energy;

  // parameters for piecewise linear curve
  float a1;
  float b1;

  float a2;
  float b2;
  float breakPoint;

#ifdef HAVE_JSON
  void get_record_from_json();
#endif
};

END_NAMESPACE_STIR

#endif