welltrack.h

Go to the documentation of this file.

#pragma once
 
/*+
________________________________________________________________________
 
 (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(); }
 
public:
 
    bool                extendIfNecessary(const Interval<float>& dahrg);
 
};
 
 
}; // namespace Well