Raster.h

Go to the documentation of this file.

 
#pragma once
 
#include <array>
#include <utility>
 
#include "MathLib/Point3d.h"
 
namespace GeoLib {
 
struct RasterHeader
{
    std::size_t n_cols; // width
    std::size_t n_rows; // height
    std::size_t n_depth; // depth (for 3d image)
    MathLib::Point3d origin; // lower left corner
    double cell_size; // edge length of each pixel
    double no_data; // no data value
};
 
class Raster {
public:
    using const_iterator = const double*;
    using iterator = double*;
 
    template <typename InputIterator>
    Raster(RasterHeader header, InputIterator begin, InputIterator end)
        : _header(std::move(header)),
          _raster_data(new double[_header.n_cols * _header.n_rows])
    {
        std::copy(begin, end, _raster_data);
    }
 
    Raster(Raster const&) = delete;
    Raster(Raster&&) = delete;
    Raster& operator=(Raster const&) = delete;
    Raster& operator=(Raster&&) = delete;
 
    RasterHeader const& getHeader() const { return _header; }
 
    void refineRaster(std::size_t scaling);
 
    const_iterator begin() const { return _raster_data; }
    const_iterator end() const { return _raster_data + _header.n_rows*_header.n_cols; }
 
    double getValueAtPoint(const MathLib::Point3d &pnt) const;
 
    double interpolateValueAtPoint(const MathLib::Point3d &pnt) const;
 
    bool isPntOnRaster(MathLib::Point3d const& pnt) const;
 
    ~Raster();
 
private:
    void setCellSize(double cell_size);
    void setNoDataVal (double no_data_val);
 
    GeoLib::RasterHeader _header;
    double* _raster_data;
};
 
}  // namespace GeoLib