Signal1D.hpp

Go to the documentation of this file.

/* ImLib3D
 * Copyright (c) 2001, ULP-IPB Strasbourg.
 *
 * 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.
 * 
 * You should have received a copy of the GNU General Public License
 * along with this program; if not, write to the Free Software
 * Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA.
 */
#ifndef _Signal1D_hpp
#define _Signal1D_hpp

#include<ImLib3D/CppTools.hpp>
#include<ImLib3D/ImLib3DFile.hpp>
#include<math.h>


class Signal1Df  : public Streamable
{
    float x0,x1;
    vector<float> values;
public:
    void ChangeSupport(float _x0,float _x1){x0=_x0;x1=_x1;}
    pair<float,float> GetSupport(){return pair<float,float>(x0,x1);}
    int NbSambles() const {return values.size();}
    int Size()      const {return values.size();}
    void GnuPlot()  const 
    {
        vector<float> coords(values.size());
        for(size_t i=0;i<values.size();i++){coords[i]=ISampleToX(i);}
        ::GnuPlot(coords,values);
    }
    float XToFSample(float x) const {return          (((x-x0)/(x1-x0))*values.size());}
    int   XToISample(float x) const {return (int)rint(((x-x0)/(x1-x0))*values.size());}
    float ISampleToX(int   i) const {return (i*(x1-x0)/values.size())+x0;}
    float &operator[](int i)        {return(values[i]);}
    float  operator[](int i) const  {return(values[i]);}
    float  operator()(float x) const{return InterpolatedValue(x);}
    float &SafeSample(size_t i)  
    {
        if(i<0){return values[0];}
        else
        if(i>=values.size()){return values.back();}
        else
        {return values[i];}
    }
    float SafeSample(size_t i)  const
    {
        if(i<0){return values[0];}
        else
        if(i>=values.size()){return values.back();}
        else
        {return values[i];}
    }
    float InterpolatedValue(float x) const
    {
        x=XToFSample(x);
        int i0=(int)x;
        float dx=x-i0;
        return(SafeSample(i0)*(1-dx)+SafeSample(i0+1)*dx);
    }
    float &Value(float rx){return SafeSample(XToISample(rx));}
    void Stats2ndOrder(float &avg,float &sigma)  const
    {
        float tot=0;
        float tot2=0;
        for(uint i=0;i<values.size();i++)
        {
            tot +=values[i];
            tot2+=values[i]*values[i];
        }
        avg=tot/values.size();
        sigma=sqrt(tot2/values.size()-avg*avg);
    }
    void Smooth(int fradius=3)  
    {
        // smooth this a bit
        vector<float> v0=values;
        for(size_t x=0;x<values.size();x++)
        {
            double s=0;
            int ct=0;
            for(int xx=(int)x-fradius;xx<=(int)x+fradius;xx++)
            {
                if(xx<0 || xx>=(int)values.size()){continue;}
                s+=v0[xx];
                ct++;
            }
            values[x]=(ct>0 ? s/ct :0);
        }
    }

    void AddSampleAtBegining(float val)
    {
        x0-=(x1-x0)/values.size();
        values.insert(values.begin(),val);
    }
    void AddSampleAtEnd(float val)
    {
        x1+=(x1-x0)/values.size();
        values.push_back(val);
    }
public:
    bool operator==(const Signal1Df &other)
    {
        return values==other.values && x0==other.x0 && x1==other.x1;
    }
    bool operator!=(const Signal1Df &other){return !operator==(other);}
        
private:
    void Read( ImLib3DFile *file,xmlpp::Element *parentNode=NULL,xmlpp::Element *node=NULL);
    void Write(ImLib3DFile *file,xmlpp::Element *parentNode=NULL,xmlpp::Element *node=NULL) const;
public:
    virtual ~Signal1Df(){}
    Signal1Df(int _size,float _x0=0,float _x1=1):
        x0(_x0),x1(_x1)
    {
        values.assign(_size,0);
    }
    Signal1Df(const string &filename)
    {
        ReadFromFile(filename);
    }

    Signal1Df():
        x0(0),x1(0)
    {
    }
};
void LinearCompleteLowConfidence(Signal1Df &fct,vector<double> &confidence,float minconf);

#endif// _Signal1D_hpp

---