ComputeIonoModel.hpp

Go to the documentation of this file.

#pragma ident "$Id: ComputeIonoModel.hpp 2526 2011-03-09 17:46:16Z yanweignss $"

#ifndef GPSTK_COMPUTE_IONO_MODEL_HPP
#define GPSTK_COMPUTE_IONO_MODEL_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
//
//  Wei Yan - Chinese Academy of Sciences . 2009, 2010, 2011
//
//============================================================================



#include "ProcessingClass.hpp"
#include "IonexStore.hpp"
#include "IonoModelStore.hpp"


namespace gpstk
{


   class ComputeIonoModel : public ProcessingClass
   {
   public:
      enum IonoModelType
      {
         Zero = 0,     
         Klobuchar,    
         Ionex,        
         DualFreq      
      };

   public:

      ComputeIonoModel()
         : ionoType(Zero), nominalPos(0.0,0.0,0.0)
      { setIndex(); };


      ComputeIonoModel(const Position& stapos) 
         : ionoType(Zero), nominalPos(stapos)
      { 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 ComputeIonoModel& setZeroModel()
      { ionoType = Zero; return (*this); }

      virtual ComputeIonoModel& setKlobucharModel(const double a[4], 
                                                  const double b[4]);

      virtual ComputeIonoModel& setKlobucharModel(const IonoModel& im);


      virtual ComputeIonoModel& setklobucharModel(const std::string& brdcFile);


      virtual ComputeIonoModel& setIonosphereMap(const std::string& ionexFile);

      virtual ComputeIonoModel& setDualFreqModel()
      { ionoType=DualFreq; return (*this); }

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


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

      virtual int getIndex(void) const;


      virtual std::string getClassName(void) const;


      virtual ~ComputeIonoModel() {};


   private:

      IonoModelType  ionoType;

      IonoModelStore klbStore;

      IonexStore gridStore;

      Position nominalPos;


      static int classIndex;


      int index;


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


   }; // End of class 'ComputeIonoModel'


} // End of namespace gpstk

#endif   // GPSTK_COMPUTE_IONO_MODEL_HPP

---