Astrometry.h

Go to the documentation of this file.

#ifndef SZA_UTIL_ASTROMETRY_H
#define SZA_UTIL_ASTROMETRY_H

#include "carma/szautil/Angle.h"
#include "carma/szautil/HourAngle.h"
#include "carma/szautil/DecAngle.h"

namespace sza {
  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);

      static void meanToApparentPlace(HourAngle& meanRa, DecAngle& meanDec, TimeVal& date, double equinox,
                                                  HourAngle& apparentRa, DecAngle& apparentDec);

      static void j2000ToApparentPlace(HourAngle& meanRa, DecAngle& meanDec, TimeVal& date, 
                                                   HourAngle& apparentRa, DecAngle& apparentDec);

      static void apparentToMeanPlace(HourAngle& apparentRa, DecAngle& apparentDec, TimeVal& date, 
                                                  double equinox,
                                                  HourAngle& meanRa, DecAngle& meanDec);

      static void apparentToJ2000Place(HourAngle& apparentRa, DecAngle& apparentDec, TimeVal& date, 
                                                   HourAngle& meanRa, DecAngle& meanDec);

      static void b1950ToJ2000(HourAngle& raB1950, DecAngle& decB1950, HourAngle& raJ2000, DecAngle& decJ2000);

      static void j2000ToB1950(HourAngle& raJ2000, DecAngle& decJ2000, HourAngle& raB1950, DecAngle& decB1950);

      static Angle angularSeparation(HourAngle ra1, DecAngle dec1, HourAngle ra2, DecAngle dec2);

    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 sza



#endif // End #ifndef SZA_UTIL_ASTROMETRY_H