welltrack.h

Go to the documentation of this file.

#ifndef welltrack_h
#define welltrack_h

/*+
________________________________________________________________________

 (C) dGB Beheer B.V.; (LICENSE) http://opendtect.org/OpendTect_license.txt
 Author:        Bert Bril
 Date:          Aug 2003
 RCS:           $Id$
________________________________________________________________________


-*/

#include "wellmod.h"
#include "welldahobj.h"
#include "position.h"

namespace Well
{

class Data;

mExpClass(Well) Track : public DahObj
{
public:

                        Track( const char* nm=0 )
                        : DahObj(nm), zistime_(false) {}
                        Track( const Track& t )
                        : DahObj("")            { *this = t; }
    Track&              operator =(const Track&);

    bool                isEmpty() const;
    const Coord3&       pos( int idx ) const    { return pos_[idx]; }
    float               value( int idx ) const  { return (float) pos_[idx].z; }
    float               getKbElev() const;
    float               td() const
                        { return isEmpty() ? 0 : dah_.last(); }
    int                 size() const    { return pos_.size(); }
    bool                zIsTime() const         { return zistime_; }
    const Interval<double> zRangeD() const;
    const Interval<float> zRange() const;

    int                 insertPoint(const Coord3&);
    int                 insertPoint(const Coord&,float z);
    void                addPoint(const Coord3&,float dah=mUdf(float));
    void                addPoint(const Coord&,float z,float dah=mUdf(float));
    void                setPoint(int,const Coord3&);
    void                setPoint(int,const Coord&,float z);
    void                insertAfterIdx(int,const Coord3&);
    void                removePoint(int);

    Coord3              getPos(float d_ah) const;
    const TypeSet<Coord3>& getAllPos() const { return pos_; }

    float               getDahForTVD(double,float prevdah=mUdf(float)) const;
    float               getDahForTVD(float,float prevdah=mUdf(float)) const;
    float               nearestDah(const Coord3&) const;
                        // If zIsTime() z must be time

                        // If you know what you're doing:
    Coord3              coordAfterIdx(float d_ah,int) const;

    bool                insertAtDah(float dah,float zpos);

    bool                alwaysDownward() const;
    void                toTime(const Data&);

protected:


    TypeSet<Coord3>     pos_;
    bool                zistime_;

    void                removeAux( int idx )    { pos_.removeSingle(idx); }
    void                eraseAux()              { pos_.erase(); }

};


}; // namespace Well

#endif