abstractdynamicvolumedata.h

#ifndef CTL_ABSTRACTDYNAMICVOLUMEDATA_H
#define CTL_ABSTRACTDYNAMICVOLUMEDATA_H

#include "spectralvolumedata.h"
#include "models/xydataseries.h"
#include "processing/range.h" // for Range<T>

/*
 * NOTE: This is header only.
 */

namespace CTL {

class AbstractDynamicVolumeData : public SpectralVolumeData
{
    // abstract interface
    protected:virtual void updateVolume() = 0;
    public: SpectralVolumeData* clone() const override = 0;

public:
    AbstractDynamicVolumeData(const SpectralVolumeData& other);

    virtual XYDataSeries timeCurve(uint x, uint y, uint z, const std::vector<float>& timePoints);

    void setTime(double seconds);
    double time() const;
    XYDataSeries timeCurve(uint x, uint y, uint z, float tStart, float tEnd, uint nbSamples);
    XYDataSeries timeCurve(uint x, uint y, uint z, SamplingRange timeRange, uint nbSamples);

protected:
    AbstractDynamicVolumeData(const AbstractDynamicVolumeData&) = default;
    AbstractDynamicVolumeData(AbstractDynamicVolumeData&&) = default;
    AbstractDynamicVolumeData& operator=(const AbstractDynamicVolumeData&) = default;
    AbstractDynamicVolumeData& operator=(AbstractDynamicVolumeData&&) = default;

private:
    double _time = 0.0; 
};

inline AbstractDynamicVolumeData::AbstractDynamicVolumeData(const SpectralVolumeData& other)
    : SpectralVolumeData(other)
{
}

inline XYDataSeries AbstractDynamicVolumeData::timeCurve(uint x, uint y, uint z,
                                                         const std::vector<float>& timePoints)
{
    XYDataSeries ret;
    const auto timePtCache = time();

    for(const auto& smp : timePoints)
    {
        setTime(smp);
        ret.append(smp, this->operator()(x,y,z));
    }

    // reset time point to initial value
    setTime(timePtCache);

    return ret;
}

inline void AbstractDynamicVolumeData::setTime(double seconds)
{
    _time = seconds;
    updateVolume();
}

inline double AbstractDynamicVolumeData::time() const { return _time; }

inline XYDataSeries AbstractDynamicVolumeData::timeCurve(uint x, uint y, uint z,
                                                         float tStart, float tEnd, uint nbSamples)
{
    return timeCurve(x, y, z, SamplingRange(tStart, tEnd).linspace(nbSamples));
}

inline XYDataSeries AbstractDynamicVolumeData::timeCurve(uint x, uint y, uint z,
                                                         SamplingRange timeRange, uint nbSamples)
{
    return timeCurve(x, y, z, timeRange.linspace(nbSamples));
}

} // namespace CTL

#endif // CTL_ABSTRACTDYNAMICVOLUMEDATA_H