IERSConventions.hpp

Go to the documentation of this file.

#pragma ident "$Id: IERSConventions.hpp 2957 2011-10-30 18:15:34Z yanweignss $"

#ifndef GPSTK_IERSCONVENTIONS_HPP
#define GPSTK_IERSCONVENTIONS_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 <cmath>
#include <string>
#include <map>

#include "DayTime.hpp"
#include "Triple.hpp"
#include "Matrix.hpp"
#include "EOPDataStore.hpp"

namespace gpstk
{

   // IERS Data Handling
   //--------------------------------------------------------------------------
   
   void LoadIERSFile(const std::string& fileName);
   
   void LoadIGSFile(const std::string& fileName);

   void LoadSTKFile(const std::string& fileName);

   EOPDataStore::EOPData EOPData(const DayTime& UTC)
      throw(InvalidRequest);
   
   double PolarMotionX(const DayTime& UTC);

   double PolarMotionY(const DayTime& UTC);

   double UT1mUTC(const DayTime& UTC);

   double NutationDPsi(const DayTime& UTC);

   double NutationDEps(const DayTime& UTC);


   int TAImUTC(DayTime UTC)
      throw(InvalidRequest);


   double TTmTAI();


   double TAImGPST();


   // Time System Handling
   //--------------------------------------------------------------------------
   enum TimeSystem
   {
      TS_Unknown = 0,
      TS_UTC,
      TS_UT1,
      TS_GPST,
      TS_TAI,
      TS_TT
   };
   
   DayTime ConvertTimeSystem(const DayTime& time, TimeSystem from, TimeSystem to);

   DayTime GPST2UTC(const DayTime& GPST);
   DayTime UTC2GPST(const DayTime& UTC);

   DayTime UT12UTC(const DayTime& UT1);
   DayTime UTC2UT1(const DayTime& UTC);

   DayTime UT12UTC(const DayTime& UT1,double ut1mutc);
   DayTime UTC2UT1(const DayTime& UTC,double ut1mutc);

   DayTime TT2UTC(const DayTime& TT);
   DayTime UTC2TT(const DayTime& UTC);

   DayTime TAI2UTC(const DayTime& TAI);
   DayTime UTC2TAI(const DayTime& UTC);

   DayTime BDT2UTC(const DayTime& BDT);
   DayTime UTC2BDT(const DayTime& UTC);
   
   // Reference System Handling
   //--------------------------------------------------------------------------

   Matrix<double> Rx(const double& angle);

   Matrix<double> Ry(const double& angle);

   Matrix<double> Rz(const double& angle);

   Matrix<double> iauPmat76(const DayTime& TT);


   static double iauNut80Args(const DayTime& TT,double& eps, double& dpsi,double& deps)
      throw(Exception);


   // IAU1976/1980 model (IERS conventions 1996)
   void J2kToECEFMatrix(const DayTime& UTC, 
                        const EOPDataStore::EOPData& ERP,
                        Matrix<double>& POM, 
                        Matrix<double>& Theta, 
                        Matrix<double>& NP)
      throw(Exception);

   Triple J2kPosToECEF(const Triple& j2kPos, const DayTime& time, TimeSystem sys = TS_GPST);

   Triple ECEFPosToJ2k(const Triple& ecefPos, const DayTime& time, TimeSystem sys = TS_GPST);

   Matrix<double> J2kToECEFMatrix(const DayTime& UTC,const EOPDataStore::EOPData& ERP);
      

   Vector<double> J2kPosToECEF(const DayTime& UTC, const Vector<double>& j2kPos)
      throw(Exception);


   Vector<double> ECEFPosToJ2k(const DayTime& UTC, const Vector<double>& ecefPos)
      throw(Exception);   


   Vector<double> J2kPosVelToECEF(const DayTime& UTC, const Vector<double>& j2kPosVel)
      throw(Exception);


   Vector<double> ECEFPosVelToJ2k(const DayTime& UTC, const Vector<double>& ecefPosVel)
      throw(Exception);   
 
 
   Vector<double> sunJ2kPosition(const DayTime& TT);
      

   Vector<double> moonJ2kPosition(const DayTime& TT);

      

   double normalizeAngle(double a);
   

   void nutationAngles(DayTime TT, double& dpsi, double& deps);

   double meanObliquity(DayTime TT);

   double iauEqeq94(DayTime TT,double eps, double dPsi);

   double iauGmst82(DayTime UT1);


   Matrix<double> iauNmat(const double& epsa, 
                          const double& dpsi, 
                          const double& deps);

   double earthRotationAngle(DayTime UT1);

   double earthRotationAngleRate1(DayTime TT);


   
} // namespace gpstk

#endif  // GPSTK_IERSCONVENTIONS2003_HPP

---