skypoint.h

Go to the documentation of this file.

/***************************************************************************
                          skypoint.h  -  K Desktop Planetarium
                             -------------------
    begin                : Sun Feb 11 2001
    copyright            : (C) 2001-2005 by Jason Harris
    email                : jharris@30doradus.org
    copyright            : (C) 2004-2005 by Pablo de Vicente
    email                : p.devicente@wanadoo.es
 ***************************************************************************/

/***************************************************************************
 *                                                                         *
 *   This program is free software; you can redistribute it and/or modify  *
 *   it under the terms of the GNU General Public License as published by  *
 *   the Free Software Foundation; either version 2 of the License, or     *
 *   (at your option) any later version.                                   *
 *                                                                         *
 ***************************************************************************/

#ifndef SKYPOINT_H_
#define SKYPOINT_H_


#include <QList>

#include "kstars/dms.h"
#include "kstars/kstarsdatetime.h"

class KSNumbers;
class KSSun;

class SkyPoint {
public:
    SkyPoint( const dms& r, const dms& d ) :
    lastPrecessJD( J2000 ), RA0(r), Dec0(d),
            RA(r),  Dec(d)
    {}

    
    //FIXME: this (*15.0) thing is somewhat hacky.
    explicit SkyPoint( double r, double d ) :
    lastPrecessJD( J2000 ), RA0(r*15.0), Dec0(d), RA(r*15.0),  Dec(d)
    {}
    
    SkyPoint();

    virtual ~SkyPoint();


    void set( const dms& r, const dms& d );

    inline void setRA0( dms r ) { RA0 = r; }

    inline void setRA0( double r ) { RA0.setH( r ); }

    inline void setDec0( dms d ) { Dec0 = d; }

    inline void setDec0( double d ) { Dec0.setD( d ); }

    inline void setRA( dms r ) { RA = r; }

    inline void setRA( double r ) { RA.setH( r ); }

    inline void setDec( dms d ) { Dec = d; }

    inline void setDec( double d ) { Dec.setD( d ); }

    inline void setAlt( dms alt ) { Alt = alt; }

    inline void setAlt( double alt ) { Alt.setD( alt ); }

    inline void setAz( dms az ) { Az = az; }

    inline void setAz( double az ) { Az.setD( az ); }


    inline const dms& ra0() const { return RA0; }

    inline const dms& dec0() const { return Dec0; }

    inline const dms& ra() const { return RA; }

    inline const dms& dec() const { return Dec; }

    inline const dms& az() const { return Az; }

    inline const dms& alt() const { return Alt; }

    dms altRefracted() const;


    void EquatorialToHorizontal( const dms* LST, const dms* lat );

    void HorizontalToEquatorial( const dms* LST, const dms* lat );

    void findEcliptic( const dms *Obliquity, dms &EcLong, dms &EcLat );

    void setFromEcliptic( const dms *Obliquity, const dms& EcLong, const dms& EcLat );

    void Equatorial1950ToGalactic(dms &galLong, dms &galLat);

    void GalacticToEquatorial1950(const dms* galLong, const dms* galLat);


    virtual void updateCoords( KSNumbers *num, bool includePlanets=true, const dms *lat=0, const dms *LST=0, bool forceRecompute = false );

    void apparentCoord(long double jd0, long double jdf);

    void nutate(const KSNumbers *num);

    bool checkBendLight();

    bool bendlight();

    SkyPoint deprecess( const KSNumbers *num, long double epoch=J2000 );

    void aberrate(const KSNumbers *num);

    void precessFromAnyEpoch(long double jd0, long double jdf);

    SkyPoint Eterms(void);

    void B1950ToJ2000(void);

    void J2000ToB1950(void);

    void addEterms(void);

    void subtractEterms(void);

    dms angularDistanceTo(const SkyPoint *sp, double * const positionAngle = 0) const;

    inline bool operator == ( SkyPoint &p ) { return ( ra() == p.ra() && dec() == p.dec() ); }

    double vRSun(long double jd);

    double vHeliocentric(double vlsr, long double jd);

    double vHelioToVlsr(double vhelio, long double jd);

    double vREarth(long double jd0);

    double vGeocentric(double vhelio, long double jd);

    double vGeoToVHelio(double vgeo, long double jd);

    double vRSite(double vsite[3]);

    double vTopocentric(double vgeo, double vsite[3]);

    double vTopoToVGeo(double vtopo, double vsite[3]);

    SkyPoint moveAway( const SkyPoint &from, double dist );

    bool checkCircumpolar( const dms *gLat );

    static double refractionCorr( double alt );

    static double refract(const double alt);

    static double unrefract(const double alt);

    static inline dms refract(const dms alt) { return dms( refract( alt.Degrees() ) ); }

    static inline dms unrefract(const dms alt) { return dms( unrefract( alt.Degrees() ) ); }

    static const double altCrit;

protected:
    void precess(const KSNumbers *num);

    long double   lastPrecessJD; // JD at which the last coordinate update (see updateCoords) for this SkyPoint was done

private:
    dms RA0, Dec0; //catalog coordinates
    dms RA, Dec; //current true sky coordinates
    dms Alt, Az;
    static KSSun *m_Sun;
};

#endif