SolarSystem.hpp

Go to the documentation of this file.

#pragma ident "$Id: SolarSystem.hpp 2293 2010-02-12 18:14:16Z btolman $"

// ======================================================================
// This software was developed by Applied Research Laboratories, The
// University of Texas at Austin under contract to an agency or agencies
// within the U.S. Department of Defense. The U.S. Government retains all
// rights to use, duplicate, distribute, disclose, or release this
// software.
// 
// Copyright 2008 The University of Texas at Austin
// ======================================================================

//------------------------------------------------------------------------------------
#ifndef SOLAR_SYSTEM_EPHEMERIS_INCLUDE
#define SOLAR_SYSTEM_EPHEMERIS_INCLUDE

//------------------------------------------------------------------------------------
// includes
// system
#include <iostream>
#include <fstream>
#include <string>
#include <vector>
#include <map>

// GPSTk
#include "Exception.hpp"
#include "DayTime.hpp"             // only for WGS84SolarSystemPosition()
#include "Position.hpp"            // only for WGS84SolarSystemPosition()
#include "EarthOrientation.hpp"    // only for WGS84SolarSystemPosition()

namespace gpstk {
//------------------------------------------------------------------------------------
class SolarSystem{
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             
   };

   SolarSystem(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; }

   gpstk::Position WGS84Position(Planet body,
                                 const gpstk::DayTime time,
                                 const gpstk::EarthOrientation& eo)
      throw(gpstk::Exception);

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 ---------------------------------------------------------

   // input stream, for use by readBinary...()
   std::ifstream istrm;  

   // header information
   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 SolarSystem

}  // end namespace gpstk

#endif // SOLAR_SYSTEM_EPHEMERIS_INCLUDE
// nothing below this

---