GloriaAstrometry.h

/*
 * GloriaAstrometry.h
 *
 *  Created on: Sep 14, 2012
 *      Author: Lech Wiktor Piotrowski
 */

#ifndef GLORIAASTROMETRY_H_
#define GLORIAASTROMETRY_H_

#include "math.h"
#include "GloriaFitsTable.h"
#include "leblind.h"

namespace gloria
{

class GloriaAstrometry: public gloria::GloriaFitsTable
{
public:
        GloriaAstrometry();

        GloriaAstrometry(leblind_struct *lestr);

        GloriaAstrometry(GloriaFitsHeader *header);

        void InitFromLeblind(leblind_struct *lestr);

        void InitFromHeader(GloriaFitsHeader *header);

        void WriteAstrometry(GloriaFitsHeader *header);

        void DeleteAstrometry(GloriaFitsHeader *header);

        bool IsAstrometryPresent(GloriaFitsHeader *header);

        void Print(std::ostream& os);

        void GetRaDecFromXY(double x, double y, double *ra, double *dec);

        inline void GetRaDecFromXY(double x, double y, double *radec)
                                { GetRaDecFromXY( x, y, radec, radec+1); };

        void GetXYFromRaDec(double ra, double dec, double *x, double *y);

        inline void GetXYFromRaDec(double ra, double dec, double *xy)
                                { GetXYFromRaDec( ra, dec, xy, xy+1); };

        virtual ~GloriaAstrometry();

public:

        // WCS astrometry values
        int crpix[2];
        double crval[2];
        double cd[4];
        int sip_order_a;
        int sip_order_b;
        int sip_order_ap;
        int sip_order_bp;

        vector< vector<double> > sip_a;

        vector< vector<double> > sip_b;

        vector< vector<double> > sip_ap;

        vector< vector<double> > sip_bp;

 private:

        double ref_ra, ref_dec;

        void InitWCSVec();

        double wcs_ref_vec[3];
        double wcs_ra_vec[3];
        double wcs_dec_vec[3];

        double deg2rad;
        double rad2deg;
};

} /* namespace gloria */
#endif /* GLORIAASTROMETRY_H_ */