EOP.hpp

浏览该文件的文档.

 
 
#pragma once
 
#include "AstGlobal.h"
#include <vector>
 
AST_NAMESPACE_BEGIN
 
 
class AST_CORE_API EOP
{
public:
    struct Entry
    {
        int mjd;                
        int dat;                
        double x, y;            
        double ut1_utc;         
        double lod;             
        double dpsi, deps;      
        double dx, dy;          
    };
public:
    EOP() = default;
    ~EOP() = default;
 
    errc_t loadDefault();
 
    errc_t load(StringView filepath);
 
    void unload();
 
    static errc_t load(StringView filepath, std::vector<Entry>& data);
 
    size_t size() const { return m_data.size(); }
 
    const Entry* getEntry(int mjd) const;
 
 
    errc_t setEntry(int mjd, const Entry& entry);
 
    //----------------
    // ut1 - utc
    //----------------
 
    double getUT1MinusUTC(const TimePoint& tp) const;
 
    double getUT1MinusUTC_UTC(const JulianDate& jdUTC) const;
 
    double getUT1MinusUTC_UTCMJD(double mjdUTC) const;
 
 
    //----------------
    // pole motion
    //----------------
 
 
    void getPoleMotion(const TimePoint& tp, double& x, double& y) const;
 
    void getPoleMotionUTC(const JulianDate& jdUTC, double& x, double& y) const;
 
 
    void getPoleMotionUTCMJD(double mjdUTC, double& x, double& y) const;
    
 
    //----------------
    // lod
    //----------------
 
    double getLOD(const TimePoint& tp) const;
 
    double getLOD_UTC(const JulianDate& jdUTC) const;
 
    double getLOD_UTCMJD(double mjdUTC) const;
 
 
    //----------------
    // xy correction
    //----------------
 
    void getXYCorrection(const TimePoint& tp, array2d& xyCorrection) const;
 
    void getXYCorrectionUTC(const JulianDate& jdUTC, array2d& xyCorrection) const;
 
    void getXYCorrectionUTCMJD(double mjdUTC, array2d& xyCorrection) const;
    
protected:
    template<double Entry::* Member>
    A_ALWAYS_INLINE double getValue(size_t index, double frac) const
    {
        return m_data[index].*Member + frac * (m_data[index+1].*Member - m_data[index].*Member);
    }
 
    void findEntryIndex(double mjdUTC, int& index, double& frac) const;
protected:
    std::vector<Entry> m_data;
    int m_startMJD = 0;
    int m_endMJD = 0;
};
 
 
AST_NAMESPACE_END