AlignmentDatabase.h

#ifndef ALIGNMENTDATABASE_H
#define ALIGNMENTDATABASE_H

#include <vector>
#include "Wm5ConvexHull3.h"
#include "Wm5Triangle3.h"
#include "Wm5Matrix3.h"

class CAlignmentDatabase
{
public:
    typedef enum { HOURS_DEGREES, DEGREES_DEGREES, RADIANS_RADIANS} CoordinateSystemUnits;

    CAlignmentDatabase();
    virtual ~CAlignmentDatabase();

    void SetLowerCoordinateSystemUnits(CoordinateSystemUnits Units) { LowerCoordinateSystemUnits = Units; }

    bool AddAlignmentPoint(const double ReferenceTime, const double RightAscension, const double Declination, const double Axis1, const double Axis2, bool Recalculate = false);

    bool RecalculateDatabase();

    bool TransformUpperToLower(const double RightAscension, const double Declination, double& Axis1, double& Axis2);

    bool TransformLowerToUpper(const double Axis1, const double Axis2, double& RightAscension, double& Declination);

    bool Save(const std::string& Filename);

    bool Load(const std::string& Filename);

    bool Append(const std::string& Filename);

protected:
private:
    static const double Epsilon = 0.0000001; // An arbitrary guess
    struct AlignmentPoint
    {
        AlignmentPoint(double ReferenceTime, double RightAscension, double Declination, double Axis1, double Axis2) :
                        ReferenceTime(ReferenceTime),
                        RightAscension(RightAscension),
                        Declination(Declination),
                        Axis1(Axis1),
                        Axis2(Axis2) {}
        // All angles are in radians
        double ReferenceTime;
        double RightAscension;
        double Declination;
        double Axis1;
        double Axis2;
    };

    std::vector<AlignmentPoint> AlignmentPoints;

    struct TinFacet
    {
        Wm5::Triangle3d Facet;
        Wm5::Matrix3d TransformationMatrix;
    };
    typedef std::vector<TinFacet> TinFacets;

    TinFacets UpperLayerTinFacets;
    TinFacets LowerLayerTinFacets;

    CoordinateSystemUnits LowerCoordinateSystemUnits;

    double RightAscensionToDummyHourAngle(const double RightAscension, const double ReferenceTime);
    double DummyHourAngleToRightAscension(const double DummyHourAngle, const double ReferenceTime);
    Wm5::Vector3d MakeUpperDirectionCosine(const double RightAscension, const double Declination, const double ReferenceTime);
    Wm5::Vector3d MakeLowerDirectionCosine(const double Axis1, const double Axis2, const double ReferenceTime);
    const double UTCNow();

    typedef enum { UPPER_TO_LOWER, LOWER_TO_UPPER} Direction;
    const Wm5::Matrix3d FindTransformationmatrix(const Wm5::Vector3d DirectionCosine, Direction Direction);
};


#endif // ALIGNMENTDATABASE_H