JplDe.hpp

浏览该文件的文档.

 
 
#pragma once
 
#include "AstGlobal.h"
#include "AstCore/RunTime.hpp"
#include <stdint.h>      // for uint32_t
#include <cstdio>       // for FILE
#include <mutex>         // for std::mutex
 
AST_NAMESPACE_BEGIN
 
 
class AST_CORE_API JplDe
{
public:
    enum EDataCode
    {
        eNotAvailable = -1,     
        eMercury = 0,           
        eVenus = 1,             
        eEarth = 2,             
        eMars = 3,              
        eJupiter = 4,           
        eSaturn = 5,            
        eUranus = 6,            
        eNeptune = 7,           
        ePluto = 8,             
        eMoon = 9,              
        eSun = 10,              
        eSSBarycenter = 11,     
        eEMBarycenter = 12,     
    };
 
    JplDe();
    ~JplDe();
 
    errc_t open(const char* filepath);
 
 
    errc_t openDefault();
 
    bool isOpen() const;
 
 
    void close();
 
 
    errc_t getPosVelICRF(
        const TimePoint& time,
        EDataCode target,
        EDataCode referenceBody,
        Vector3d& pos,
        Vector3d& vel
    );
 
    errc_t getPosVelICRF(
        const TimePoint& time,
        EDataCode target,
        EDataCode referenceBody,
        Vector3d& pos,
        Vector3d* vel
    );
 
 
    errc_t getPosICRF(
        const TimePoint& time,
        EDataCode target,
        EDataCode referenceBody,
        Vector3d& pos
    );
    
 
    errc_t  getNutation(
        const TimePoint& time,
        double& nutLong,
        double& nutObl
    );
 
    
    errc_t      getLibration(
        const TimePoint& time,
        Vector3d& ang,
        Vector3d& angRate
    );
    
    errc_t       getLibration(
        const TimePoint& time,
        Vector3d& ang
    );
 
    errc_t       getLibration(
        const TimePoint& time,
        Euler& ang
    );
    
    uint32_t getEphemVersion() const{ return m_EphemVerion; }
 
 
    errc_t getInterval(TimeInterval& interval) const;
 
private:
    errc_t         readDataBlock(size_t idx);
    errc_t         getStateTDB(const JulianDate& jdTDB, int dataid, double pos[], double vel[]);
    errc_t         getStateTDB(const JulianDate& jdTDB, int datalist[11], double statelist[11][6]);
 
    errc_t         getState(const TimePoint& time, int dataid, double pos[], double vel[]);
    errc_t         getState(const TimePoint& time, int datalist[11], double statelist[11][6]);
protected:
    bool        m_IsSameEndian;        // 二进制文件与系统大小端是否一致
    uint32_t    m_EphemVerion{ 0 };    // De星历版本
    uint32_t    m_NumConstants{};        // 文件的常量个数
    uint32_t    m_NumDataBlock{0};     // 星历总的数据块个数      
    uint32_t    m_NumCoeff{};            // 星历每个数据块的系数数量
    uint32_t    m_ipt[15][3]{};        // 数据指针[数据起始位置；每个坐标分量的数据量；子块个数：如4代表了32天要分为四块，每个子块包含了8天的数据]
    double      m_EphemStart{0};       // 星历起始时间(JED)
    double      m_EphemEnd{0};         // 星历结束时间(JED)
    double      m_EphemStep{0};        // 星历每个数据块的天数
    double      m_AU{};                  // 光速 km
    double      m_EMMassRatio{};         // 地球月球质量比
    FILE*       m_DeFile{ NULL };      // De二进制文件
    double**    m_DataBlocks{NULL};    // 星历数据块的内存缓存
    std::mutex  m_DataBlockMutex;      // 数据块缓存互斥锁
};
 
 
 
AST_NAMESPACE_END