seisblocksdata.h

Go to the documentation of this file.

#pragma once

/*
________________________________________________________________________

 (C) dGB Beheer B.V.; (LICENSE) http://opendtect.org/OpendTect_license.txt
 Author:        Bert
 Date:          Apr 2017
________________________________________________________________________

*/

#include "seiscommon.h"

class DataBuffer;
namespace Survey { class Geometry3D; }
template <class T> class DataInterpreter;


namespace Seis
{

namespace Blocks
{
    typedef short               IdxType;
    typedef unsigned short      SzType;
    typedef Survey::Geometry3D  SurvGeom;

#define mDefSeisBlockTripletClass(clss,typ) \
mExpClass(Seis) clss : public Triplet<typ> \
{ \
public: \
\
    inline      clss()                                  {} \
    inline      clss( short iidx, short xidx, short zidx ) \
                    : Triplet<typ>(iidx,xidx,zidx)      {} \
    inline bool operator ==( const clss& oth ) const \
                { return first == oth.first && second == oth.second \
                      && third == oth.third; } \
\
    inline typ  inl() const     { return first; } \
    inline typ& inl()           { return first; } \
    inline typ  crl() const     { return second; } \
    inline typ& crl()           { return second; } \
    inline typ  z() const       { return third; } \
    inline typ& z()             { return third; } \
}

mDefSeisBlockTripletClass(GlobIdx,IdxType);
mDefSeisBlockTripletClass(SampIdx,IdxType);
mDefSeisBlockTripletClass(Dimensions,SzType);


mExpClass(Seis) Block
{
public:

                        Block(GlobIdx,SampIdx start=SampIdx(),
                                Dimensions dims=defDims(),
                                OD::FPDataRepType fpr=OD::F32);
                        ~Block();

    const GlobIdx&      globIdx() const         { return globidx_; }
    const SampIdx&      start() const           { return start_; }
    const Dimensions&   dims() const            { return dims_; }
    const DataBuffer&   dataBuf() const         { return dbuf_; }

    void                zero();
    void                retire();
    bool                isRetired() const;

    float               value(const SampIdx&) const;
    void                getVert(SampIdx,float*,int sz) const;
    void                setValue(const SampIdx&,float);
    void                setVert(SampIdx,const float*,int sz);

    SampIdx             getSampIdx(const BinID&,const SurvGeom&) const;
    SampIdx             getSampIdx(const BinID&,float z,const SurvGeom&) const;
    IdxType             getSampZIdx(float,const SurvGeom&) const;

    static IdxType      globIdx4Inl(const SurvGeom&,int inl,SzType inldim);
    static IdxType      globIdx4Crl(const SurvGeom&,int crl,SzType crldim);
    static IdxType      globIdx4Z(const SurvGeom&,float z,SzType zdim);
    static IdxType      sampIdx4Inl(const SurvGeom&,int inl,SzType inldim);
    static IdxType      sampIdx4Crl(const SurvGeom&,int crl,SzType crldim);
    static IdxType      sampIdx4Z(const SurvGeom&,float z,SzType zdim);
    static int          inl4Idxs(const SurvGeom&,SzType inldim,IdxType globidx,
                                IdxType sampidx);
    static int          crl4Idxs(const SurvGeom&,SzType crldim,IdxType globidx,
                                IdxType sampidx);
    static float        z4Idxs(const SurvGeom&,SzType zdim,IdxType globidx,
                                IdxType sampidx);

    static Dimensions   defDims();

protected:

    const GlobIdx       globidx_;
    const SampIdx       start_;
    const Dimensions    dims_;
    DataBuffer&         dbuf_;
    const DataInterpreter<float>* interp_;

    int                 getBufIdx(const SampIdx&) const;
    inline int          nrSampsPerInl() const
                        { return ((int)dims_.crl()) * dims_.z(); }

};


mExpClass(Seis) DataStorage
{
public:

    virtual             ~DataStorage()          {}

    unsigned short      version() const         { return version_; }

    static BufferString fileNameFor(const GlobIdx&);

    static const char*  sKeyFileType()   { return "Column Cube"; }
    static const char*  sKeyDimensions() { return "Blocks.Max Dimensions"; }
    static const char*  sKeyGlobInlRg()  { return "Blocks.Inl ID Range"; }
    static const char*  sKeyGlobCrlRg()  { return "Blocks.Crl ID Range"; }
    static const char*  sKeyGlobZRg()    { return "Blocks.Z ID Range"; }
    static const char*  sKeyComponents() { return "Components"; }

protected:

                        DataStorage(const SurvGeom*);

    const SurvGeom&     survgeom_;
    const Dimensions    dims_;
    unsigned short      version_;

};


} // namespace Blocks

} // namespace Seis