PlanetEphemeris.hpp

Go to the documentation of this file.

#pragma ident "$Id: PlanetEphemeris.hpp 2950 2011-10-27 16:21:52Z yanweignss $"

#ifndef GPSTK_PLANETEPHEMERIS_HPP
#define GPSTK_PLANETEPHEMERIS_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 . 2011
//
//============================================================================

#include <iostream>
#include <fstream>
#include <string>
#include <vector>
#include <map>

#include "Exception.hpp"
#include "DayTime.hpp" 
#include "Position.hpp"           

namespace gpstk
{
   class PlanetEphemeris
   {
   public:
      enum Planet {
         // the following are relative to the solar system barycenter, except MOON
         None=0,                
         Mercury,               
         Venus,                 
         Earth,                 
         Mars,                  
         Jupiter,               
         Saturn,                
         Uranus,                
         Neptune,               
         Pluto,                 
         Moon,                  
         Sun,                   
         SolarSystemBarycenter, 
         EarthMoonBarycenter,   
         Nutations,             
         Librations             
      };

      PlanetEphemeris(void) throw() : EphemerisNumber(-1) {};

      void readASCIIheader(std::string filename) throw(gpstk::Exception);

      int readASCIIdata(std::vector<std::string>& filenames) throw(gpstk::Exception);

      int readASCIIdata(std::string filename) throw(gpstk::Exception);

      int writeASCIIheader(std::ostream& os) throw(gpstk::Exception);

      int writeASCIIdata(std::ostream& os) throw(gpstk::Exception);

      int writeBinaryFile(std::string filename) throw(gpstk::Exception);

      void clearStorage(void) throw() { store.clear(); }

      int readBinaryFile(std::string filename) throw(gpstk::Exception);

      int initializeWithBinaryFile(std::string filename) throw(gpstk::Exception);

      int computeState(double tt,
         Planet target,
         Planet center,
         double PV[6],
         bool kilometers = true)
         throw(gpstk::Exception);

      double AU(void) throw()
      { if(EphemerisNumber == -1) return -1.0; return constants["AU"]; }

      int JPLNumber(void) const throw()
      { return EphemerisNumber; }

      double getConstant(std::string name) throw() {
         if(EphemerisNumber == -1) return -1.0;
         if(constants.find(name) != constants.end()) return constants[name];
         return 0.0;
      }

      gpstk::DayTime startTime(void) const throw(gpstk::Exception)
      { gpstk::DayTime t; t.setMJD(startJD - gpstk::DayTime::JD_TO_MJD); return t; }

      gpstk::DayTime endTime(void) const throw(gpstk::Exception)
      { gpstk::DayTime t; t.setMJD(endJD - gpstk::DayTime::JD_TO_MJD); return t; }


   private:
      void writeBinary(std::ofstream& strm, const char *ptr, size_t size)
         throw(gpstk::Exception);

      void readBinary(char *ptr, size_t size) throw(gpstk::Exception);

      void readBinaryHeader(std::string filename) throw(gpstk::Exception);

      int readBinaryData(bool save) throw(gpstk::Exception);

      int readBinaryRecord(std::vector<double>& data_vector) throw(gpstk::Exception);

      enum computeID {
         // the following are relative to the solar system barycenter, except MOON
         NONE=-1,        
         MERCURY,        
         VENUS,          
         EMBARY,         
         MARS,           
         JUPITER,        
         SATURN,         
         URANUS,         
         NEPTUNE,        
         PLUTO,          
         MOON,           
         SUN,            
         NUTATIONS,      
         LIBRATIONS      
      };

      int seekToJD(double JD) throw(gpstk::Exception);

      void computeState(double tt, computeID which, double PV[6])
         throw(gpstk::Exception);

      // member data ---------------------------------------------------------

      std::ifstream istrm;  

      int EphemerisNumber;

      int Ncoeff;
      int Nconst;           
      std::string label[3]; 
      double startJD;       
      double endJD;         
      double interval;      
      int c_offset[13];     
      int c_ncoeff[13];     
      int c_nsets[13];      

      std::map<std::string,double> constants;

      std::map<double, std::vector<double> > store;

      std::map<double, long> fileposMap;

      std::vector<double> coefficients;

   }; // end class PlanetEphemeris
 
}   // End of namespace gpstk


#endif  //GPSTK_PLANETEPHEMERIS_HPP

---