Meteo1DInterpolator.h

Go to the documentation of this file.

// SPDX-License-Identifier: LGPL-3.0-or-later
/***********************************************************************************/
/*  Copyright 2009 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 METEO1DINTERPOLATOR_H
#define METEO1DINTERPOLATOR_H
 
#include <meteoio/dataClasses/MeteoData.h>
#include <meteoio/Config.h>
#include <meteoio/meteoResampling/ResamplingAlgorithms.h>
#include <meteoio/meteoFilters/ProcessingBlock.h> //for ProcessingProperties
 
#include <string>
#include <vector>
#include <map>
#include <memory>
 
namespace mio {
class ResamplingStack {
    public: 
                ResamplingStack();
 
        void addAlgorithm(std::shared_ptr<ResamplingAlgorithms> algo, const double& max_gap_size);
        std::vector<std::shared_ptr<ResamplingAlgorithms>> buildStack(const ResamplingAlgorithms::gap_info& gap) const;
 
                void resetResampling();
                void resample(const std::string &stationHash, const size_t &index, const ResamplingAlgorithms::ResamplingPosition elementpos, const size_t &ii, const std::vector<MeteoData> &vecM, MeteoData &md, const double& i_window_size) const;
                std::string getStackStr() const;
                bool empty() const;
 
    private:
        std::vector<double> max_gap_sizes;
        std::vector<std::shared_ptr<ResamplingAlgorithms>> stack;
};
 
class Meteo1DInterpolator {
        public:
 
                Meteo1DInterpolator(const Config& in_cfg, const char& rank=1, const IOUtils::OperationMode &mode=IOUtils::STD);
                Meteo1DInterpolator(const Meteo1DInterpolator& org) = default;
 
                bool resampleData(const Date& date, const std::string& stationHash, const std::vector<MeteoData>& vecM, MeteoData& md);
                
                void resetResampling();
 
                void getWindowSize(ProcessingProperties& o_properties) const;
 
                Meteo1DInterpolator& operator=(const Meteo1DInterpolator&); 
                const std::string toString() const;
 
        private:
                std::vector< std::pair<std::string, std::string> > getArgumentsForAlgorithm(const std::string& parname, const std::string& algorithm) const;
                std::string getAlgorithmsForParameter(const std::string& parname) const;
 
                void processAlgorithms(bool first_time, const std::string& parname, const std::vector<std::pair<std::string, std::string>>& vecAlgos, const IOUtils::OperationMode& mode=IOUtils::STD, const char& rank=1);
                void createResamplingStacks(const IOUtils::OperationMode& mode, const char& rank);
 
                std::map< std::string, ResamplingStack > mapAlgorithms; //per parameter interpolation algorithms
                const Config& cfg;
                double window_size; 
                bool enable_resampling, data_qa_logs; 
        
        public:
                static const std::string interpol_section;
                static const std::string interpol_pattern;
                static const std::string arg_pattern;
};
 
 
} //end namespace
 
#endif