CelestialBody.hpp

浏览该文件的文档.

 
 
#pragma once
 
#include "AstGlobal.h"
#include "AstCore/Object.hpp"
#include "AstCore/GravityField.hpp"
#include "AstCore/JplDe.hpp"
#include "AstCore/BodyEphemeris.hpp"
#include "AstCore/BodyOrientation.hpp"
#include "AstCore/AxesBodyInertial.hpp"
#include "AstCore/AxesBodyFixed.hpp"
#include "AstCore/AxesBodyMOD.hpp"
#include "AstCore/AxesBodyTOD.hpp"
#include "AstCore/Point.hpp"
#include "AstCore/Frame.hpp"
#include "AstCore/AxesAPI.hpp"
#include "AstCore/BodyShape.hpp"
#include "AstUtil/SharedPtr.hpp"
#include "AstUtil/ScopedPtr.hpp"
 
AST_NAMESPACE_BEGIN
 
class CelestialBody;
using Body = CelestialBody;
 
 
enum class EEphemerisSource
{
    eBodyEphemeris,      
    eJplDE,              
    eJplSpice,           
    eJplSpiceBarycenter, 
};
 
 
 
AST_CORE_CAPI errc_t aGetGravityParameter(const Body& body, StringView gravityModel, double& gm);
 
 
class AST_CORE_API CelestialBody : public Point
{
public:
    AST_OBJECT(CelestialBody)
    AST_PROPERT(Shape)
    AST_PROPERT(Ephemeris)
    AST_PROPERT(Orientation)
    static CelestialBody* Resolve(StringView value);
    
    CelestialBody();
    explicit CelestialBody(SolarSystem* solarSystem);
    explicit CelestialBody(CelestialBody* parentBody);
    ~CelestialBody();
 
    int getJplSpiceId() const { return jplSpiceId_; }
    void setJplSpiceId(int id) { jplSpiceId_ = id; }
 
    int getJplIndex() const { return jplIndex_; }
    void setJplIndex(int index);
 
    std::string getDirpath() const;
 
    SolarSystem* getSolarSystem() const;
    
    CelestialBody* getParent() const { return parent_.get(); }
public:
 
    const std::string& getGravityModel() const{ return gravityField_.getModelName(); }
 
    errc_t load(StringView filepath);
 
    double getRadius() const { return radius_; }
 
    double getGM() const { return gm_; }
 
    double getSystemGM() const { return systemGM_; }
 
    errc_t setGravityModel(StringView model);
 
    double getGravityRefDistance() const { return gravityField_.refDistance(); }
    
    double getJ2() const { return getJn(2); }
 
    double getJ3() const { return getJn(3); }
    
    double getJ4() const { return getJn(4); }
    
    double getJ5() const { return getJn(5); }
    
    double getJ6() const { return getJn(6); }
 
    bool isEarth() const { return jplIndex_ == JplDe::eEarth; }
 
    bool isLuna() const { return jplIndex_ == JplDe::eMoon; }
 
    errc_t getPosICRF(const TimePoint& tp, Vector3d& pos) const;
 
 
    errc_t getPosVelICRF(const TimePoint& tp, Vector3d& pos, Vector3d& vel) const;
 
public: // 从Point继承重写的函数
 
    Frame* getFrame() const final;
    errc_t getPos(const TimePoint& tp, Vector3d& pos) const final;
    errc_t getPosVel(const TimePoint& tp, Vector3d& pos, Vector3d& vel) const final;
 
PROPERTIES: // 天体的形状、重力场、星历、姿态
 
    BodyShape* getShape() const { return shape_.get(); }
    BodyShape* shape() const {return shape_.get(); }
 
    BodyEphemeris* getEphemeris() const { return ephemeris_.get(); }
    BodyEphemeris* ephemeris() const {return ephemeris_.get(); }
 
    BodyOrientation* getOrientation() const { return orientation_.get(); }
    BodyOrientation* orientation() const {return orientation_.get(); }
public:
    BodyEphemeris* getEphemeris(EEphemerisSource ephemerisSource) const;
public:
    const GravityField& getGravityField() const { return gravityField_; }
    const GravityField& gravityField() const {return gravityField_; }
 
public:
 
    Point* getPointCenter() const { return const_cast<CelestialBody*>(this); }
    Point* pointCenter() const {return const_cast<CelestialBody*>(this); }
 
    Axes* getAxesInertial() const { return axesInertial_.get(); }
    Axes* axesInertial() const {return axesInertial_.get(); }
    
    Axes* getAxesFixed() const { return axesFixed_.get(); }
    Axes* axesFixed() const {return axesFixed_.get(); }
 
    Axes* getAxesMOD() const { return axesMOD_.get(); }
    Axes* axesMOD() const {return axesMOD_.get(); }
 
    Axes* getAxesTOD() const { return axesTOD_.get(); }
    Axes* axesTOD() const {return axesTOD_.get(); }
 
    Axes* getAxes(StringView name) const;
    
#ifndef SWIG // 还没解决好swig封装智能指针的问题，暂时屏蔽
public:
    HAxes makeEpochAxes(Axes* sourceAxes, const TimePoint& tp, Axes* reference) const{return aMakeEpochAxes(sourceAxes, tp, reference);}
 
    HAxes makeEpochAxes(Axes* sourceAxes, EventTime* time, Axes* reference) const{return aMakeEpochAxes(sourceAxes, time, reference);}
 
 
    HFrame makeEpochFrame(Axes* sourceAxes, const TimePoint& tp, Axes* reference) const;
 
    HFrame makeEpochFrame(Axes* sourceAxes, EventTime* time, Axes* reference) const;
 
    HAxes makeAxesMOE(const TimePoint& tp) const;
 
    HAxes makeAxesMOE(EventTime* time) const;
 
    HAxes makeAxesTOE(const TimePoint& tp) const;
 
    HAxes makeAxesTOE(EventTime* time) const;
 
    HFrame makeFrameMOE(const TimePoint& tp) const;
 
    HFrame makeFrameMOE(EventTime* time) const;
 
    HFrame makeFrameTOE(const TimePoint& tp) const;
 
    HFrame makeFrameTOE(EventTime* time) const;
 
public:
    HFrame makeFrame(Axes* axes) const;
 
    HFrame makeFrameInertial() const;
 
    HFrame makeFrameFixed() const;
 
    HFrame makeFrameMOD() const;
 
    HFrame makeFrameTOD() const;
 
    HFrame makeFrameJ2000() const;
 
    HFrame makeFrameICRF() const;
public:
    Frame* getFrameInertial() const;
    Frame* frameInertial() const {return getFrameInertial();}
 
    Frame* getFrameFixed() const;
    Frame* frameFixed() const {return getFrameFixed();}
 
    Frame* getFrameICRF() const;
 
    Frame* getFrameJ2000() const;
 
#endif
protected:
 
    Axes* getEpochAxesReference() const;
protected:
    double getJn(int n) const { return gravityField_.getJn(n); }
 
    errc_t loadGravityModel(StringView model);
protected:
    friend class SolarSystem;
 
    errc_t loadAstroDefinition(BKVParser& parser);
 
    errc_t loadSpinData(BKVParser& parser);
 
    errc_t loadEphemerisData(BKVParser& parser);
 
    errc_t loadEarth(BKVParser& parser);   
    
    errc_t loadMoon(BKVParser& parser);   
 
    errc_t loadMeanEarthDefinition(BKVParser& parser);
 
    A_DISABLE_COPY(CelestialBody)
private:
    WeakPtr<SolarSystem>        solarSystem_{};                    
    SharedPtr<CelestialBody>    parent_{};                         
    double                      gm_{0.0};                        
    double                      systemGM_{0.0};                  
    double                      radius_{0.0};                    
    int                         jplSpiceId_{-1};                 
    int                         jplIndex_{-1};                   
    GravityField                gravityField_{};                   
    SharedPtr<BodyShape>        shape_{};                          
    SharedPtr<BodyOrientation>  orientation_{};                    
    SharedPtr<BodyEphemeris>    ephemeris_{};                      
 
    SharedPtr<AxesBodyInertial> axesInertial_{};                   
    SharedPtr<AxesBodyFixed>    axesFixed_{};                      
    SharedPtr<AxesBodyMOD>      axesMOD_{};                        
    SharedPtr<AxesBodyTOD>      axesTOD_{};                        
 
    mutable WeakPtr<Frame>      frameInertial_{};                  
    mutable WeakPtr<Frame>      frameFixed_{};                     
 
    mutable SharedPtr<BodyEphemeris> ephemerisDE_{};              
    mutable SharedPtr<BodyEphemeris> ephemerisSpice_{};           
    mutable SharedPtr<BodyEphemeris> ephemerisSpiceBarycenter_{}; 
};
 
 
enum ESpiceId: int;
 
AST_CORE_CAPI ESpiceId aGetPlanetBarycenterId(ESpiceId planetId);
 
 
AST_CORE_CAPI errc_t aBodyOrbitNormalIn
(
    Axes& axes,
    const Body& body,
    const TimePoint& tp,
    Vector3d& normal
);
 
 
AST_CORE_CAPI errc_t aBodyOrbitNormalInParentInertial
(
    const Body& body,
    const TimePoint& tp,
    Vector3d& normal
);
 
 
AST_CORE_CAPI errc_t aBodyOrbitNormalInICRF
(
    const Body& body,
    const TimePoint& tp,
    Vector3d& normal
);
 
 
 
enum ESpiceId: int
{
    eSolarSystemBarycenter  = 0,    
    eMercuryBarycenter      = 1,    
    eVenusBarycenter        = 2,    
    eEarthBarycenter        = 3,    
    eMarsBarycenter         = 4,    
    eJupiterBarycenter      = 5,    
    eSaturnBarycenter       = 6,    
    eUranusBarycenter       = 7,    
    eNeptuneBarycenter      = 8,    
    ePlutoBarycenter        = 9,    
 
    eSun                    = 10,   
 
    eMercury                = 199,  
 
    eVenus                  = 299,  
 
    eEarth                  = 399,  
    eMoon                   = 301,  
 
    eMars                   = 499,  
    ePhobos                 = 401,
    eDeimos                 = 402,
 
 
    eJupiter                = 599,  
    eIo                     = 501,
    eEuropa                 = 502,
    eGanymede               = 503,
    eCallisto               = 504,
    eAmalthea               = 505,
    eHimalia                = 506,
    eElara                  = 507,
    ePasiphae               = 508,
    eSinope                 = 509,
    eLysithea               = 510,
    eCarme                  = 511,
    eAnanke                 = 512,
    eLeda                   = 513,
    eThebe                  = 514,
    eAdrastea               = 515,
    eMetis                  = 516,
 
    eSaturn                 = 699,  
    eMimas                  = 601,
    eEnceladus              = 602,
    eTethys                 = 603,
    eDione                  = 604,
    eRhea                   = 605,
    eTitan                  = 606,
    eHyperion               = 607,
    eIapetus                = 608,
    ePhoebe                 = 609,
    eJanus                  = 610,
    eEpimetheus             = 611,
    eHelene                 = 612,
    eTelesto                = 613,
    eCalypso                = 614,
    eAtlas                  = 615,
    ePrometheus             = 616,
    ePandora                = 617,
    ePan                    = 618,
    eMethone                = 632,
    ePallene                = 633,
    ePolydeuces             = 634,
    eDaphnis                = 635,
    eAnthe                  = 649,
    eAegaeon                = 653,
 
 
    eUranus                 = 799,  
    eAriel                  = 701,
    eUmbriel                = 702,
    eTitania                = 703,
    eOberon                 = 704,
    eMiranda                = 705,
    eCordelia               = 706,
    eOphelia                = 707,
    eBianca                 = 708,
    eCressida               = 709,
    eDesdemona              = 710,
    eJuliet                 = 711,
    ePortia                 = 712,
    eRosalind               = 713,
    eBelinda                = 714,
    ePuck                   = 715,
 
 
    eNeptune                = 899,  
    eTriton                 = 801,
    eNereid                 = 802,
    eNaiad                  = 803,
    eThalassa               = 804,
    eDespina                = 805,
    eGalatea                = 806,
    eLarissa                = 807,
    eProteus                = 808,
 
    ePluto                  = 999,  
    eCharon                 = 901,
};
 
 
 
using HCelestialBody = SharedPtr<CelestialBody>;        
using WCelestialBody = WeakPtr<CelestialBody>;          
using HBody = HCelestialBody;                           
using WBody = WCelestialBody;                           
 
 
using PCelestialBody = CelestialBody*;                  
using PBody = PCelestialBody;                           
 
 
AST_NAMESPACE_END
 
 
#include "SolarSystem.hpp"