EarthOrientation.hpp

Go to the documentation of this file.

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

//============================================================================
//
//  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
//  
//  Copyright 2004, The University of Texas at Austin
//
//============================================================================

//============================================================================
//
//This software developed by Applied Research Laboratories at 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. 
//
//Pursuant to DoD Directive 523024 
//
// DISTRIBUTION STATEMENT A: This software has been approved for public 
//                           release, distribution is unlimited.
//
//=============================================================================

#ifndef CLASS_EARTHORIENT_INCLUDE
#define CLASS_EARTHORIENT_INCLUDE

//------------------------------------------------------------------------------------
// system includes
#include <iostream>
#include <iomanip>
#include <string>
#include <map>
// GPSTk
#include "Exception.hpp"
#include "DayTime.hpp"

//------------------------------------------------------------------------------------
namespace gpstk {

   class EarthOrientation
   {
   public:
      double xp;           // arcseconds
      double yp;           // arcseconds
      double UT1mUTC;      // seconds

      EarthOrientation() : xp(0),yp(0),UT1mUTC(0) {}

      friend std::ostream& operator<<(std::ostream& s, const EarthOrientation&);
   };

   //---------------------------------------------------------------------------------
   class EOPPrediction
   {
      double tv;
           double ta,tb;
           double A,B,C1,C2,D1,D2,E,F,G1,G2,H1,H2,I,J,K1,K2,K3,K4,L1,L2,L3,L4;
           double P1,P2,Q1,Q2,R1,R2,R3,R4;

   public:
           int TAIUTC;
      int SerialNo;
      std::string Info;
   
      int getValidTime(void) const
         throw() { return int(tv+0.5); }

      int loadFile(std::string filename)
         throw(FileMissingException);
   
      static int getSerialNumber(DayTime& t)
         throw(DayTime::DayTimeException);
   
      static int getSerialNumber(int mjd)
         throw(DayTime::DayTimeException)
         {
            DayTime t;
            t.setMJD(double(mjd));
            return getSerialNumber(t);
         }

      EarthOrientation computeEOP(int& mjd) const
         throw(DayTime::DayTimeException);

      EarthOrientation computeEOP(DayTime& t) const
         throw();
   
      friend std::ostream& operator<<(std::ostream& s, const EOPPrediction&);
   };

   //---------------------------------------------------------------------------------
   class EOPStore
   {
      std::map<int,EarthOrientation> mapMJD_EOP;

      int begMJD,endMJD;
   
   public:
      EOPStore() : begMJD(-1), endMJD(-1) {}

      void addEOP(int MJD,
                  EarthOrientation& eop)
         throw();

      int addEOP(int MJD,
                 EOPPrediction& eopp)
         throw(DayTime::DayTimeException);

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

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

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

      void edit(int mjdmin,
                int mjdmax)
         throw();

      void edit(const DayTime& tmin,
                const DayTime& tmax)
         throw()
         { edit(int(tmin.MJD()+0.5),int(tmax.MJD()+1.5)); }

      int size(void)
         throw() { return mapMJD_EOP.size(); }

      void clear(void)
         throw() { mapMJD_EOP.clear(); begMJD=endMJD=-1; }

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

      int getFirstMJD(void) throw()
      { return begMJD; }

      int getLastMJD(void) throw()
      { return endMJD; }

      EarthOrientation getEOP(DayTime& t) const
         throw(InvalidRequest);

   };

}  // end namespace gpstk

#endif // nothing below this

---