Astrometry.h

Go to the documentation of this file.

#ifndef GCP_UTIL_ASTROMETRY_H
#define GCP_UTIL_ASTROMETRY_H

#include "gcp/util/common/Angle.h"
#include "gcp/util/common/HourAngle.h"

namespace gcp {
  namespace util {
    
    class QuadraticInterpolator;
    
    class Astrometry {
    public:
      
      struct Date {
    int year_;
    int month_;
    int day_;
    int hour_;
    int min_;
    int sec_;
    int nsec_;
      };

      Astrometry();
      
      virtual ~Astrometry();
      
      void extendUt1Utc(double mjd, double ut1utc);
      
      void extendEqnEqx(double tt, double eqneqx);
      
      double getUt1Utc(double mjd);
      
      double getEqnEqx(double tt);
      
      bool canBracket(double mjdUtc);
      
      static HourAngle mjdUtcToLst(double mjdUtc, Angle longitude, 
                           double ut1Utc, double eqnEqx);
      
      HourAngle mjdUtcToLst(double mjdUtc, Angle longitude);
      
      static double mjdUtcToMjdTt(double mjdUtc);
      static Date mjdUtcToCalendarDate(double mjdUtc);
      
    private:
      
      static const double secondsPerDay_;
      static const double pi_;
      static const double twopi_;
      
      QuadraticInterpolator* ut1Utc_;
      
      QuadraticInterpolator* eqnEqx_;
      
    }; // End class Astrometry
    
  } // End namespace util
} // End namespace gcp




#endif // End #ifndef GCP_UTIL_ASTROMETRY_H