ComputeWindUp.hpp

Go to the documentation of this file.

#pragma ident "$Id: ComputeWindUp.hpp 1662 2009-02-04 09:00:34Z architest $"

#ifndef GPSTK_COMPUTEWINDUP_HPP
#define GPSTK_COMPUTEWINDUP_HPP

//============================================================================
//
//  This file is part of GPSTk, the GPS Toolkit.
//
//  The GPSTk is free software; you can redistribute it and/or modify
//  it under the terms of the GNU Lesser General Public License as published
//  by the Free Software Foundation; either version 2.1 of the License, or
//  any later version.
//
//  The GPSTk is distributed in the hope that it will be useful,
//  but WITHOUT ANY WARRANTY; without even the implied warranty of
//  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
//  GNU Lesser General Public License for more details.
//
//  You should have received a copy of the GNU Lesser General Public
//  License along with GPSTk; if not, write to the Free Software Foundation,
//  Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307  USA
//
//  Dagoberto Salazar - gAGE ( http://www.gage.es ). 2007, 2008, 2009
//
//============================================================================



#include <string>
#include "ProcessingClass.hpp"
#include "Triple.hpp"
#include "Position.hpp"
#include "SunPosition.hpp"
#include "XvtStore.hpp"
#include "SatDataReader.hpp"
#include "geometry.hpp"



namespace gpstk
{



   class ComputeWindUp : public ProcessingClass
   {
   public:

      ComputeWindUp()
         : pEphemeris(NULL), nominalPos(0.0, 0.0, 0.0),
           satData("PRN_GPS"), fileData("PRN_GPS")
      { setIndex(); };


      ComputeWindUp( XvtStore<SatID>& ephem,
                     const Position& stapos,
                     string filename="PRN_GPS" )
         : pEphemeris(&ephem), nominalPos(stapos), satData(filename),
           fileData(filename)
      { setIndex(); };


      virtual satTypeValueMap& Process( const DayTime& time,
                                        satTypeValueMap& gData )
         throw(ProcessingException);


      virtual gnssSatTypeValue& Process(gnssSatTypeValue& gData)
         throw(ProcessingException)
      { Process(gData.header.epoch, gData.body); return gData; };


      virtual gnssRinex& Process(gnssRinex& gData)
         throw(ProcessingException)
      { Process(gData.header.epoch, gData.body); return gData; };


      virtual string getFilename(void) const
      { return fileData; };


      virtual ComputeWindUp& setFilename(const string& name);


      virtual Position getNominalPosition(void) const
      { return nominalPos; };


      virtual ComputeWindUp& setNominalPosition(const Position& stapos)
        { nominalPos = stapos; return (*this); };


      virtual XvtStore<SatID> *getEphemeris(void) const
      { return pEphemeris; };


      virtual ComputeWindUp& setEphemeris(XvtStore<SatID>& ephem)
      { pEphemeris = &ephem; return (*this); };


      virtual int getIndex(void) const;


      virtual std::string getClassName(void) const;


      virtual ~ComputeWindUp() {};


   private:


      XvtStore<SatID> *pEphemeris;


      Position nominalPos;


      SatDataReader satData;


      string fileData;


      struct phaseData
      {
            // Default constructor initializing the data in the structure
         phaseData() : previousPhase(0.0) {};

         double previousPhase;      
      };


      map<SatID, phaseData> phase_station;


      map<SatID, phaseData> phase_satellite;


      map<SatID, double> satArcMap;


      virtual double getWindUp( const SatID& sat,
                                const DayTime& time,
                                const Triple& satpos,
                                const Triple& sunpos );


      static int classIndex;

      int index;

      void setIndex(void)
      { index = classIndex++; };


   }; // End of class 'ComputeWindUp'




}  // End of namespace gpstk

#endif // GPSTK_COMPUTEWINDUP_HPP

---