ComputeMOPSWeights.hpp

Go to the documentation of this file.

#pragma ident "$Id: ComputeMOPSWeights.hpp 2479 2010-10-26 08:02:55Z architest $"

#ifndef COMPUTEMOPSWEIGHTS_HPP
#define COMPUTEMOPSWEIGHTS_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 ). 2006, 2007, 2008, 2010
//
//============================================================================



#include "DayTime.hpp"
#include "DataStructures.hpp"
#include "WeightBase.hpp"
#include "EngEphemeris.hpp"
#include "TabularEphemerisStore.hpp"
#include "GPSEphemerisStore.hpp"
#include "ComputeIURAWeights.hpp"
#include "TropModel.hpp"
#include "geometry.hpp"             // DEG_TO_RAD


namespace gpstk
{



   class ComputeMOPSWeights : public ComputeIURAWeights
   {
   public:

      ComputeMOPSWeights() : receiverClass(2), defaultIono(TypeID::ionoL1)
      { pBCEphemeris = NULL; pTabEphemeris = NULL; setIndex(); }


      ComputeMOPSWeights( const Position& pos,
                          GPSEphemerisStore& bcephem,
                          int rxClass = 2 )
         : receiverClass(rxClass), nominalPos(pos), defaultIono(TypeID::ionoL1)
      { setDefaultEphemeris(bcephem); setIndex(); };


      ComputeMOPSWeights( const Position& pos,
                          TabularEphemerisStore& tabephem,
                          int rxClass = 2 )
         : receiverClass(rxClass), nominalPos(pos), defaultIono(TypeID::ionoL1)
      { setDefaultEphemeris(tabephem); setIndex(); };


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


      virtual gnssSatTypeValue& Process(gnssSatTypeValue& gData)
         throw(ProcessingException);


      virtual gnssRinex& Process(gnssRinex& gData)
         throw(ProcessingException);


      virtual ComputeMOPSWeights& setPosition(const Position& pos)
      { nominalPos = pos; return (*this); };


      virtual ComputeMOPSWeights& setDefaultEphemeris(XvtStore<SatID>& ephem);


      virtual ComputeMOPSWeights& setDefaultEphemeris(GPSEphemerisStore& ephem)
      { pBCEphemeris = &ephem; pTabEphemeris = NULL; return (*this); };


      virtual ComputeMOPSWeights& setDefaultEphemeris(
                                             TabularEphemerisStore& ephem )
      { pBCEphemeris = NULL; pTabEphemeris = &ephem; return (*this); };


      virtual TypeID getDefaultIono() const
      { return defaultIono; };


      virtual ComputeMOPSWeights& setDefaultIono(const TypeID& type)
      { defaultIono = type; return (*this); };


      virtual int getIndex(void) const;


      virtual std::string getClassName(void) const;


      virtual ~ComputeMOPSWeights() {};


   private:


      int receiverClass;


      Position nominalPos;


      TypeID defaultIono;


      virtual double getWeight( const SatID& sat,
                                typeValueMap& tvMap )
         throw(InvalidWeights);


         // Compute ionospheric sigma^2 according to Appendix J.2.3
         // and Appendix A.4.4.10.4 in MOPS-C
      double sigma2iono( const double& ionoCorrection,
                         const double& elevation,
                         const double& azimuth,
                         const Position& rxPosition )
         throw(InvalidWeights);


      static int classIndex;

      int index;

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


   }; // End of class 'ComputeMOPSWeights'


}  // End of namespace gpstk
#endif   // COMPUTEMOPSWEIGHTS_HPP

---