IonexStore.hpp

Go to the documentation of this file.

#pragma ident "$Id: IonexStore.hpp 2741 2011-06-22 16:37:02Z nwu $"

#ifndef GPSTK_IONEXSTORE_HPP
#define GPSTK_IONEXSTORE_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
//
//  Octavian Andrei - FGI ( http://www.fgi.fi ). 2008, 2009
//
//============================================================================


#include <map>

#include "FileStore.hpp"
#include "IonexData.hpp"

#include "geometry.hpp"                   // DEG_TO_RAD
#include "icd_200_constants.hpp"          // LX_FREQ, with X = 1,2,5,6,7,8
#include "Triple.hpp"

namespace gpstk
{


   class IonexStore : public FileStore<IonexHeader>
   {
   public:


      IonexStore()
         throw()
         : initialTime(DayTime::END_OF_TIME),
           finalTime(DayTime::BEGINNING_OF_TIME)
      {};


      virtual ~IonexStore() {};


      virtual void loadFile(const std::string& filename)
         throw(FileMissingException);


      void addMap(const IonexData& iod)
         throw();


      void dump( std::ostream& s = std::cout,
                 short detail = 0 ) const
         throw();


      void clear() throw();


      Triple getIonexValue( const DayTime& t,
                            const Position& RX,
                            int strategy = 3 ) const
         throw(InvalidRequest);



   double getSTEC( const double& elevation,
                   const double& tecval,
                   const std::string& ionoMapType ) const
      throw (InvalidParameter);



      double getIono( const double& elevation,
                      const double& tecval,
                      const double& freq,
                      const std::string& ionoMapType ) const
         throw (InvalidParameter);


         // The next 6 functions define the interface for calculating
         // the ionospheric slant delay for a specific frequency

      double getIonoL1( const double& elevation,
                        const double& tecval,
                        const std::string& ionoMapType ) const
         throw (InvalidParameter)
      { return getIono(elevation, tecval, L1_FREQ, ionoMapType); };


      double getIonoL2( const double& elevation,
                        const double& tecval,
                        const std::string& ionoMapType ) const
         throw (InvalidParameter)
      { return getIono(elevation, tecval, L2_FREQ, ionoMapType); };


      double getIonoL5( const double& elevation,
                        const double& tecval,
                        const std::string& ionoMapType ) const
         throw (InvalidParameter)
      { return getIono(elevation, tecval, L5_FREQ, ionoMapType); };


      double getIonoL6( const double& elevation,
                        const double& tecval,
                        const std::string& ionoMapType ) const
         throw (InvalidParameter)
      { return getIono(elevation, tecval, L6_FREQ, ionoMapType); };


      double getIonoL7( const double& elevation,
                        const double& tecval,
                        const std::string& ionoMapType ) const
         throw (InvalidParameter)
      { return getIono(elevation, tecval, L7_FREQ, ionoMapType); };


      double getIonoL8( const double& elevation,
                        const double& tecval,
                        const std::string& ionoMapType ) const
         throw (InvalidParameter)
      { return getIono(elevation, tecval, L8_FREQ, ionoMapType); };


      double iono_mapping_function( const double& elevation,
                                    const std::string& ionoMapType ) const;


      DayTime getInitialTime() const
         throw(InvalidRequest)
      { return initialTime; }


      DayTime getFinalTime() const
         throw(InvalidRequest)
      { return finalTime; }


      double findDCB( const SatID sat,
                      const DayTime& time ) const
         throw(InvalidRequest);


   private:


      DayTime initialTime, finalTime;


      typedef std::map<IonexData::IonexValType, IonexData> IonexValTypeMap;


      typedef std::map<DayTime, IonexValTypeMap> IonexMap;


      IonexMap inxMaps;


      typedef std::map<DayTime, IonexHeader::SatDCBMap> IonexDCBMap;


      IonexDCBMap inxDCBMap;


   }; // End of class 'IonexStore'


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

---