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();
 
    JplDe(const JplDe&) = delete;
    JplDe& operator=(const JplDe&) = delete;
 
    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,
        Euler& ang,
        Euler& angRate
    );
    
 
    errc_t       getLibration(
        const TimePoint& time,
        Euler& ang
    );
    
    uint32_t getEphemVersion() const{ return ephemVersion_; }
 
 
    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        isSameEndian_;        // 二进制文件与系统大小端是否一致
    uint32_t    ephemVersion_{ 0 };    // De星历版本
    uint32_t    numConstants_{};        // 文件的常量个数
    uint32_t    numDataBlock_{0};     // 星历总的数据块个数      
    uint32_t    numCoeff_{};            // 星历每个数据块的系数数量
    uint32_t    ipt_[15][3]{};        // 数据指针[数据起始位置；每个坐标分量的数据量；子块个数：如4代表了32天要分为四块，每个子块包含了8天的数据]
    double      ephemStart_{0};       // 星历起始时间(JED)
    double      ephemEnd_{0};         // 星历结束时间(JED)
    double      ephemStep_{0};        // 星历每个数据块的天数
    double      au_{};                  // 光速 km
    double      emMassRatio_{};         // 地球月球质量比
    FILE*       deFile_{ NULL };      // De二进制文件
    double**    dataBlocks_{NULL};    // 星历数据块的内存缓存
    std::mutex  dataBlockMutex_;      // 数据块缓存互斥锁
};
 
 
 
AST_NAMESPACE_END