LoadFromBST.h

Go to the documentation of this file.

#pragma once

#include "VolumeSize.h"
#include "VoxelSize.h"

#include <Eigen/Core>
#include <ImageStack/ImageStack.h>
#include <ImageStack/ImageStackLoaderBST.h>
#include <ImageStack/ResolutionDecorator.h>
#include <gsl/gsl>

namespace CortidQCT {

namespace IO {

template <class F> void loadFromBST(std::string const &filename, F &&f) {
  using Decorator = ImageStack::ResolutionDecorator;
  using Img = ImageStack::ImageStack<float, ImageStack::HostStorage, Decorator>;
  using Loader = ImageStack::ImageStackLoaderBST<Img>;

  auto const img = Img(Loader(filename));

  auto const mapping = img.map();
  auto const size = img.size().template cast<std::size_t>().eval();
  auto const resol = img.resolution.template cast<float>().eval();
  VoxelSize voxelSize{resol.x(), resol.y(), resol.z()};
  VolumeSize volumeSize{size.x(), size.y(), size.z()};

  f(mapping.data(), volumeSize, voxelSize);
}

} // namespace IO

} // namespace CortidQCT