welld2tmodel.h

Go to the documentation of this file.

#ifndef welld2tmodel_h
#define welld2tmodel_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 "uistring.h"

class TimeDepthModel;
class uiString;

namespace Well
{

class Data;
class Track;
class Info;

mExpClass(Well) D2TModel : public DahObj
{ mODTextTranslationClass(D2TModel);
public:

                        D2TModel( const char* nm= 0 )
                        : DahObj(nm)    {}
                        D2TModel( const D2TModel& d2t )
                        : DahObj("")    { *this = d2t; }
    D2TModel&           operator =(const D2TModel&);
    bool                operator ==(const D2TModel&) const;
    bool                operator !=(const D2TModel&) const;

    float               getTime(float d_ah, const Track&) const;
    float               getDepth(float time, const Track&) const;
    float               getDah(float time, const Track&) const;
    double              getVelocityForDah(float d_ah,const Track&) const;
    double              getVelocityForDepth(float dpt,const Track&) const;
    double              getVelocityForTwt(float twt,const Track&) const;
    bool                getTimeDepthModel(const Well::Data&,
                                          TimeDepthModel&) const;

    inline float        t( int idx ) const      { return t_[idx]; }
    float               value( int idx ) const  { return t(idx); }
    float*              valArr()                { return t_.arr(); }
    const float*        valArr() const          { return t_.arr(); }

    BufferString        desc;
    BufferString        datasource;

    static const char*  sKeyTimeWell(); 
    static const char*  sKeyDataSrc();

    void                add( float d_ah, float tm )
                                                { dah_ += d_ah; t_ += tm; }
    bool                insertAtDah(float d_ah,float t);

    void                makeFromTrack(const Track&, float cstvel,
                                      float replvel);
    bool                ensureValid(const Well::Data&,uiString& errmsg,
                                    TypeSet<double>* zvals=0,
                                    TypeSet<double>* tvals=0);

protected:

    TypeSet<float>      t_;

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

    bool                getVelocityBoundsForDah(float d_ah,const Track&,
                                          Interval<double>& depths,
                                          Interval<float>& times) const;
    bool                getVelocityBoundsForTwt(float twt,const Track&,
                                          Interval<double>& depths,
                                          Interval<float>& times) const;
    int                 getVelocityIdx(float pos,const Track&,
                                       bool posisdah=true) const;

protected:

    inline float        getDepth( float time ) const { return mUdf(float); }
    bool                getOldVelocityBoundsForDah(float d_ah,const Track&,
                                             Interval<double>& depths,
                                             Interval<float>& times) const;
    bool                getOldVelocityBoundsForTwt(float twt,const Track&,
                                             Interval<double>& depths,
                                             Interval<float>& times) const;
    static bool getTVDD2TModel( Well::D2TModel& d2t, const Well::Data& wll,
                          TypeSet<double>& zvals, TypeSet<double>& tvals,
                          uiString& errmsg, uiString& warnmsg );
    static void convertDepthsToMD( const Well::Track& track,
                                  const TypeSet<double>& zvals,
                                  TypeSet<float>& dahs );
    static void shiftTimesIfNecessary( TypeSet<double>& tvals, double wllheadz,
                                 double vrepl, double origintwtinfile,
                                  uiString& msg );
    static void checkReplacementVelocity( Well::Info& info,
                                         double vreplinfile,
                                         uiString& msg );
};

mGlobal(Well) float     getDefaultVelocity();


}; // namespace Well

#endif