u3d_seis.hpp

/* u3d_seis.hpp */
/* $Id: u3d_seis.hpp,v 1.1 2008/05/13 12:19:18 sad Exp $ */
//------------------------------------------------------------------------------
#ifndef __u3d_seis_hpp
#define __u3d_seis_hpp

#include "u3d_os.hpp"
#include <mth/u3d_box.hpp>

enum U3D_OS_TYPE { OFFSET_VSP, WALKAWAY_VSP };

//==============================================================================
class U3dAbstractSeis 
{
public:
  
  U3dAbstractSeis();
  
  U3dAbstractSeis(U3dPointsAr sp_ar, U3dPointsAr rec_ar, U3dModAreaPos modelOrigin);
  
  virtual ~U3dAbstractSeis();
  
  int    getSPNum()  { return m_sp_ar.size(); }
  
  int    getRecNum() { return m_rec_ar.size(); }
  
  U3dObservSys* getOS(int spId);
  
  bool   getSP(int spId, U3dPoint &sp);
  
  bool   getRec(int recId, U3dPoint &recp);
  
  bool   getReady() { return m_ready; }

  
  virtual bool loadSPTraces(int spId) { return false; }
  
  double getAmpl(int recId, double time, double f, double az);
  
  double getAmplPure(int recId, double time, double f, double az);
  double getAmplParabolic(int recId, double time, double f, double az);
  double getAmplWindowSumm(int recId, double time, double f, double az);
  
  
  double getA(int recId, int sample_num, double f, double az);
  

  
  int getTraceLength() { return m_trLength; }
  
  double getTraceStep() { return m_trStep; }
  
  U3dBox getOSBox();
  
  /* DEBUG */
  float** getTraceXComp() { return m_seisDataXComp; }
  float** getTraceYComp() { return m_seisDataYComp; }
  float** getTraceZComp() { return m_seisDataZComp; }
  
  void    setReady(bool flag) { m_ready = flag; }
  
protected:
  
  void   sortRec(U3dPointsAr &recs);
  
  void   setSpsNRecs(U3dPointsAr sp_ar, U3dPointsAr recs_ar);
  
  bool       m_ready;             
  int        m_compNum;            
  int        m_trLength;          
  double     m_trStep;            
  int        m_currentSP;
  
  U3dModAreaPos m_modelOrigin;    /*< always in UNIVERS coordinate system */
  
  float     **m_seisDataXComp;  
  float     **m_seisDataYComp;  
  float     **m_seisDataZComp;  
  U3dPointsAr m_sp_ar;              
  U3dPointsAr m_rec_ar;             
};

//------------------------------------------------------------------------------
#endif /* u3d_seis.hpp */