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