dataclipper.h

Go to the documentation of this file.

#pragma once

/*+
________________________________________________________________________

 (C) dGB Beheer B.V.; (LICENSE) http://opendtect.org/OpendTect_license.txt
 Author:        Kristofer Tingdahl
 Date:          09-02-2002
________________________________________________________________________


-*/

#include "algomod.h"
#include "typeset.h"
#include "ranges.h"
template <class T> class ValueSeries;
template <class T> class ArrayND;

mExpClass(Algo) DataClipper
{
public:
                                DataClipper();
    bool                        operator ==(const DataClipper&) const;
    inline bool                 isEmpty() const { return samples_.isEmpty(); }

    void                        setApproxNrValues(od_int64 nrsamples,
                                                  int statsize=2000);
    void                        putData(float);
    void                        putData(const float*,od_int64 sz);
    void                        putData(const ValueSeries<float>&,od_int64 sz);
    void                        putData(const ArrayND<float>&);

    bool                        calculateRange(float cliprate,Interval<float>&);
    bool                        calculateRange(float lowcliprate,
                                               float highcliprate,
                                               Interval<float>&);
    static bool                 calculateRange(float* vals, od_int64 nrvals,
                                               float lowcliprate,
                                               float highcliprate,
                                               Interval<float>&);
    bool                        fullSort();
    bool                        getRange(float cliprate,Interval<float>&) const;
    bool                        getRange(float lowcliprate,float highcliprate,
                                         Interval<float>&) const;
    bool                        getSymmetricRange(float cliprate,float midval,
                                         Interval<float>&) const;
    void                        reset();

    const LargeValVec<float>&   statPts() const { return samples_; }

protected:

    int                         approxstatsize_;
    float                       sampleprob_;
    bool                        subselect_;
    LargeValVec<float>          samples_;
    Interval<float>             absoluterg_;
};


mExpClass(Algo) DataClipSampler
{
public:
                        DataClipSampler(int bufsz=10000);
                        ~DataClipSampler()      { delete [] vals_; }

    void                reset()                 { count_ = 0; }

    void                add(float);
    void                add(const float*,od_int64);
    void                finish() const;

    od_int64            nrVals() const;
    const float*        vals() const            { return vals_; }

    void                report(IOPar&) const;
    Interval<float>     getRange(float clipratio) const;


protected:

    float*              vals_;
    const int           maxnrvals_;
    od_int64            count_;
    Interval<float>     rg_;
    bool                finished_;

    void                doAdd(float);
    BufferString        getClipRgStr(float) const;
};