AdaptiveDespiker.h

Go to the documentation of this file.

// -*- mode: C++; indent-tabs-mode: nil; c-basic-offset: 4; tab-width: 4; -*-
// vim: set shiftwidth=4 softtabstop=4 expandtab:
/*
 ********************************************************************
 ** NIDAS: NCAR In-situ Data Acquistion Software
 **
 ** 2006, Copyright University Corporation for Atmospheric Research
 **
 ** This program is free software; you can redistribute it and/or modify
 ** it under the terms of the GNU General Public License as published by
 ** the Free Software Foundation; either version 2 of the License, or
 ** (at your option) any later version.
 **
 ** This program 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 General Public License for more details.
 **
 ** The LICENSE.txt file accompanying this software contains
 ** a copy of the GNU General Public License. If it is not found,
 ** write to the Free Software Foundation, Inc.,
 ** 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA.
 **
 ********************************************************************
*/
 
#ifndef NIDAS_CORE_ADAPTIVE_DESPIKER_H
#define NIDAS_CORE_ADAPTIVE_DESPIKER_H
 
#include <stdexcept>
 
namespace nidas { namespace core {
 
class AdaptiveDespiker
{
 
public:
    AdaptiveDespiker();
 
    void setOutlierProbability(float val);
 
    float getOutlierProbability() const { return _prob; }
 
    void setDiscLevelMultiplier(float val);
 
    float getDiscLevelMultiplier() const { return _levelMultiplier; }
 
    float getDiscLevel() const { return _level; }
 
    float despike(float u, bool* spike);
 
    float forecast() const
    {
    return _u1 * _corr + (1. - _corr) * _mean2;
    }
 
    void reset();
 
    static float discrLevel(float prob);
 
    static float adjustLevel(float corr);
 
private:
 
    static bool staticInit();
 
    static const size_t STATISTICS_SIZE = 100;
 
    static const int ADJUST_TABLE_SIZE = 100;
 
    static const int LEN_ERFC_ARRAY = 140;
 
    static float _adj[ADJUST_TABLE_SIZE][2];
 
    static bool staticInitDone;
 
    static void spline(float* x,float* y,int n,double yp1,double ypn,float* y2);
 
    static double splint (float* xa,float* ya,float* y2a,int n, double x);
 
    void initStatistics(float u);
 
    void incrementStatistics(float u);
 
    void updateStatistics(float u);
 
 
    float _prob;
 
    float _levelMultiplier;
 
    float _maxMissingFreq;
 
    float _u1;  
 
    double _mean1;
 
    double _mean2;
 
    double _var1;
 
    double _var2;
 
 
    double _corr;
 
    double _initLevel;
 
    double _level;
 
    float _missfreq;
 
    int _msize;
 
    size_t _npts;
};
 
}}  // namespace nidas namespace core
 
#endif