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)
        throw(std::range_error);

    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