ResamplingStack.h

Go to the documentation of this file.

// SPDX-License-Identifier: LGPL-3.0-or-later
/***********************************************************************************/
/*  Copyright 2026 WSL Institute for Snow and Avalanche Research    SLF-DAVOS      */
/***********************************************************************************/
/* This file is part of MeteoIO.
    MeteoIO is free software: you can redistribute it and/or modify
    it under the terms of the GNU Lesser General Public License as published by
    the Free Software Foundation, either version 3 of the License, or
    (at your option) any later version.
 
    MeteoIO is distributed in the hope that it will be useful,
    but WITHOUT ANY WARRANTY; without even the implied warranty of
    MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
    GNU Lesser General Public License for more details.
 
    You should have received a copy of the GNU Lesser General Public License
    along with MeteoIO.  If not, see <http://www.gnu.org/licenses/>.
*/
#ifndef RESAMPLINGSTACK_H
#define RESAMPLINGSTACK_H
 
#include <meteoio/dataClasses/MeteoData.h>
#include <meteoio/meteoResampling/ResamplingAlgorithms.h>
#include <meteoio/Config.h>
 
#include <string>
#include <vector>
#include <memory>
 
namespace mio {
 
class ResamplingStack {
    public:
                ResamplingStack(const Config& cfg, const std::string& i_gap_size_key, const std::string& parname, const IOUtils::OperationMode &mode=IOUtils::STD, const char& rank=1, const bool& isDefault=false, const double& i_fallback_gap_size=0.0);
 
                std::vector<std::shared_ptr<ResamplingAlgorithms>> buildStack(const ResamplingAlgorithms::gap_info& gap) const;
 
                void resetResampling();
                void resample(const std::string &stationHash, const size_t &pos, const ResamplingAlgorithms::ResamplingPosition elementpos, const size_t &paramindex, const std::vector<MeteoData> &vecM, MeteoData &md) const;
                std::string getStackStr() const;
                bool isEmpty() const {return stack.empty();}
 
        private:
                // Helper methods for stack construction
                double getAlgoMaxGapSize(const Config& cfg, const size_t& algo_index, const std::string& parname) const;
                double getStackMaxGapSize(const double& fallback) const;
                void processAlgorithm(const Config& cfg, const std::string& parname,
                                     const std::pair<std::string, std::string>& key_algo,
                                     const IOUtils::OperationMode &mode,
                                     const char& rank);
                static size_t extractAlgorithmIndex(const std::string& key);
                static void checkDeprecatedSyntax(const Config& cfg, const std::string& parname,
                                          const std::string& algo_name);
                void handleAccumulateAlgorithm(std::vector<std::pair<std::string, std::string>>& vecArgs,
                                              const std::string& algo_name,
                                              const Config& cfg,
                                              const IOUtils::OperationMode &mode,
                                              const char& rank);
 
    private:
                std::vector<std::shared_ptr<ResamplingAlgorithms>> stack;
                static const std::string interpol_section, interpol_pattern, arguments_ini_key;
                const std::string gap_size_key;
                const double max_gap_size_dflt, max_stack_gap_size; //both in days
};
 
} //end namespace
 
#endif