MeasurementModel.h

Go to the documentation of this file.

#pragma once

#include "DiscreteRange.h"
#include "Optional.h"

#include <chrono>
#include <set>
#include <string>
#include <unordered_map>
#include <vector>

namespace CortidQCT {

#pragma clang diagnostic push
#pragma clang diagnostic ignored "-Wpadded"

class MeasurementModel {
public:
  using Label = unsigned int;

private:
  struct VOIData {
    Label label;
    double scale;
    std::string name;
    std::vector<float> data;
  };
  using DataStorage = std::unordered_map<Label, VOIData>;

public:

  std::optional<std::string> name;
  std::optional<std::string> description;
  std::optional<std::string> author;
  std::optional<std::string> creationDate;

  float kernelSigma = 1.f;
  float sliceSpacing = 1.f;

  DiscreteRange<float> samplingRange{-2.f, 2.f, 0.1f};
  DiscreteRange<float> densityRange{-1000.f, 2000.f, 1.f};
  DiscreteRange<float> angleRange{0.f, 90.f, 1.f};



  inline MeasurementModel() noexcept(noexcept(DataStorage())) {}

  explicit MeasurementModel(std::string const &filename);

  inline static MeasurementModel fromFile(std::string const &filename) {
    return MeasurementModel(filename);
  }



  MeasurementModel &loadFromFile(std::string const &filename);



  inline bool isEmpty() const noexcept { return data_.empty(); }

  inline auto labels() const {
    std::set<Label> labels;
    for (auto &&entry : data_) { labels.insert(entry.first); }

    return labels;
  }

  inline std::size_t labelCount() const noexcept { return data_.size(); }

  inline double densityScale(Label const &label) const {
    if (auto const voiIt = data_.find(label); voiIt != data_.end()) {
      return voiIt->second.scale;
    }
    return std::numeric_limits<double>::quiet_NaN();
  }


  template <class F>
  inline auto withUnsafeDataPointer(Label label, F &&f) const {
    if (auto const store = data_.find(label); store != data_.end()) {
      return f(store->second.data.data());
    }

    throw std::invalid_argument("Label " + std::to_string(label) +
                                " not found");

    // Should never be called, but GCC complains if nothing is returned.
    return f(nullptr);
  }


private:

  void reorderData();

  DataStorage data_;
};
#pragma clang diagnostic pop

} // namespace CortidQCT