prestackattrib.h

Go to the documentation of this file.

#ifndef prestackattrib_h
#define prestackattrib_h

/*+
________________________________________________________________________

 (C) dGB Beheer B.V.; (LICENSE) http://opendtect.org/OpendTect_license.txt
 Author:        B.Bril & H.Huck
 Date:          14-01-2008
 RCS:           $Id$
________________________________________________________________________

-*/


#include "attributesmod.h"
#include "attribprovider.h"
#include "prestackprop.h"
#include "multiid.h"
#include "datapack.h"

class SeisPSReader;
class IOObj;

namespace PreStack { class AngleComputer; class Gather;
                     class GatherSetDataPack; class ProcessManager; }


namespace Attrib
{

mExpClass(Attributes) PSAttrib : public Provider
{ mODTextTranslationClass(PSAttrib)
public:

    static void         initClass();

                        PSAttrib(Desc&);

    static const char*  attribName()            { return "PreStack"; }
    static const char*  offStartStr()           { return "offstart"; }
    static const char*  offStopStr()            { return "offstop"; }
    static const char*  preProcessStr()         { return "preprocessor"; }
    static const char*  calctypeStr()           { return "calctype"; }
    static const char*  stattypeStr()           { return "stattype"; }
    static const char*  lsqtypeStr()            { return "lsqtype"; }
    static const char*  offsaxisStr()           { return "offsaxis"; }
    static const char*  valaxisStr()            { return "valaxis"; }
    static const char*  componentStr()          { return "comp"; }
    static const char*  apertureStr()           { return "aperture"; }
    static const char*  velocityIDStr()         { return "velocityid"; }
    static const char*  useangleStr()           { return "useangle"; }
    static const char*  angleStartStr()         { return "anglestart"; }
    static const char*  angleStopStr()          { return "anglestop"; }
    static const char*  rayTracerParamStr()     { return "raytracerparam"; }
    static const char*  gathertypeStr()         { return "gathertype"; }
    static const char*  xaxisunitStr()          { return "xaxisunit"; }
    static const char*  angleDPIDStr()          { return "angleid"; }
    static const char*  angleSmoothType()       { return "anglesmoothtype"; }
    static const char*  angleFiltFunction()     { return "anglefiltfunc"; }
    static const char*  angleFiltValue()        { return "anglefiltval"; }
    static const char*  angleFiltLength()       { return "anglefiltlen"; }
    static const char*  angleFFTF3Freq()        { return "anglefftf3freq"; }
    static const char*  angleFFTF4Freq()        { return "anglefftf4freq"; }

    enum GatherType     { Off, Ang };
                        mDeclareEnumUtils(GatherType)
    enum XaxisUnit      { Deg, Rad };
                        mDeclareEnumUtils(XaxisUnit)

    const PreStack::PropCalc::Setup&    setup() const   { return setup_; }
    const MultiID&                      psID() const    { return psid_; }
    const MultiID&                      preProcID() const { return preprocid_; }
    const MultiID&                      velocityID() const
                                        { return velocityid_; }

    void                updateCSIfNeeded(TrcKeyZSampling&) const;
    void                setAngleComp(PreStack::AngleComputer*);
    void                setAngleData(DataPack::ID);

protected:

                        ~PSAttrib();
    static Provider*    createInstance(Desc&);
    static void         updateDesc(Desc&);

    bool                allowParallelComputation() const        { return true;}
    bool                getInputOutput(int input,TypeSet<int>& res) const;
    bool                getInputData(const BinID&, int idx);
    bool                getAngleInputData();
    bool                computeData(const DataHolder&,const BinID& relpos,
                                    int t0,int nrsamples,int threadid) const;
    void                prepPriorToBoundsCalc();
    void                setSmootheningPar();
    float               getXscaler(bool isoffset, bool isindegrees) const;
    DataPack::ID        getPreProcessedID(const BinID& relbid);
    bool                getGatherData(const BinID& bid,DataPack::ID& gatherid,
                                      DataPack::ID& angglegid);

    MultiID                     psid_;
    IOObj*                      psioobj_;
    SeisPSReader*               psrdr_;
    int                         component_;
    PreStack::ProcessManager*   preprocessor_;
    PreStack::PropCalc*         propcalc_;
    PreStack::PropCalc::Setup   setup_;
    PreStack::AngleComputer*    anglecomp_;
    int                         anglegsdpid_;

    MultiID                     preprocid_;
    int                         dataidx_;
    const DataHolder*           inputdata_;
    MultiID                     velocityid_;

    ObjectSet<PreStack::Gather>    gatherset_;

private:

    PreStack::GatherSetDataPack*        getMemoryGatherSetDP() const;
    void                        setGatherIsAngle(PreStack::Gather&);
};

}; // namespace Attrib

#endif