sc_star_catalog.h

Go to the documentation of this file.

//
//  sc_star_catalog.h
//  SCControl
//
//  Created by Stephanie Thomas on 11/17/15.
//
//
 
#ifndef __SCControl__sc_star_catalog__
#define __SCControl__sc_star_catalog__
 
#include <stdio.h>
#include <map>
#include <vector>
#include <string>
#include <algorithm>
using namespace std; // It may be a bad idea to put this in a header, but it works.
 
#ifdef AS_OS_WINDOWS
#include "matrixlib.h"
#else
#include <MatrixLib/MatrixLib.h>
#endif
 
class star
{
    private:
    int hipparcos_id;
    double right_ascension;
    double declination;
    double visual_magnitude;
    double parallax;
    double proper_motion_ascension;
    double proper_motion_declination;
    double radial_velocity;
    
public:
    star();
    star(int id, double VM, double right_ascension, double declination);
    void set_motion_terms(double parallax, double pm_ascension, double pm_declination, double radialV);
    ml_matrix CalculateVector( double jd_now, double jd_catalog ) const;
    ml_matrix UnitVector( void ) const;
    int get_id() const;
    double get_ascension() const;
    double get_declination() const;
    double get_magnitude() const;
};
 
class star_catalog
{
    friend class near_matrix;
private:
    map<int,star> star_map;
    double julianDate;
    string name;
    void LoadStarMatrix(ml_matrix starData);
    bool status;
public:
    star_catalog();
    star_catalog(ml_matrix starData,double date=2448347.75);
    star_catalog(const char* fileName,double date=2448347.75);
    void set_date(double newDate);
    void set_name(string newName);
    star* get_star(int id);
    const double get_date();
    const string get_name();
    const int size();
    const ml_matrix unit_vectors(double jd_now);
    const ml_matrix visual_magnitude();
    const ml_int_array hipparcos_ids();
    int remove_stars(double angleTol);
    bool is_valid() {return status;};
};
 
enum DotDirection { FORWARD=1, BACK=-1, NEUTRAL=0 };
typedef struct star_dot {
    double angle;
    DotDirection dir = NEUTRAL;
    int star_1;
    int star_2;
    int getFirst(void) {if(dir>=0) return star_1; else return star_2; };
    int getSecond(void) {if(dir>=0) return star_2; else return star_1; };
} StarDot;
// think I can write an operator for comparison in find() algorithm
bool operator== (const star_dot& i, const star_dot& j);
bool operator== (const star_dot& i, int j);
 
typedef struct star_triad {
    int star_1;
    int star_2;
    int star_3;
} StarTriad;
 
class near_matrix
{
private:
    double q;
    double slope;
    ml_int_array k_vector;
    vector<StarDot> dot_products;
public:
    near_matrix();
    near_matrix(star_catalog catalog, double field_of_view);
    near_matrix(ml_int_array ids, ml_matrix u_star, double field_of_view);
    vector<StarDot> get_angle_range(double angle, double tol);
    vector<StarDot> get_nearby_stars(int star_id);
    int size();
    
};
 
 
#endif /* defined(__SCControl__sc_star_catalog__) */