TinyGPSPlus based LST calculator

Author: MichelDequick

Software/Arduino code/Sidereal.cpp

@@ -0,0 +1,30 @@
+#include "Sidereal.hpp"
+
+static DayTime Sidereal::calculateByGPS(TinyGPSPlus* gps){
+    DayTime timeUTC = DayTime(gps->time.hour(), gps->time.minute(), gps->time.second());
+    int deltaJd = calculateDeltaJd(gps->date.year(), gps->date.month(), gps->date.day());
+    float deltaJ = (float) deltaJd + ((timeUTC.getTotalHours()) / 24.0);
+    return DayTime((float) (calculateTheta(deltaJ, (float) gps->location.lng(), timeUTC.getTotalHours()) / 15.0));
+}
+
+static const float Sidereal::calculateTheta(float deltaJ, float longitude, float timeUTC){
+    float theta = C1;
+    theta += C2 * deltaJ;
+    theta += C3 * timeUTC;
+    theta += longitude;
+    return fmod(theta, 360.0);
+}
+
+static const int Sidereal::calculateDeltaJd(int year, int month, int day){
+    const int daysInMonth[12] = {31, 28, 31, 30, 31, 30, 31, 31, 30, 31, 30, 31};
+
+    // Calculating days without leapdays
+    long int deltaJd = (year - REFERENCE_JEAR) * 365 + day; 
+    for (int i=0; i < month - 1; i++){deltaJd += daysInMonth[i];}
+
+    // Add leapdays
+    if (month <= 2){year--;}
+    deltaJd += (year / 4 - year / 100 + year / 400);
+    deltaJd -= (REFERENCE_JEAR / 4 - REFERENCE_JEAR / 100 + REFERENCE_JEAR / 400);
+    return deltaJd; 
+}

Software/Arduino code/Sidereal.hpp

@@ -0,0 +1,24 @@
+#pragma once
+
+#include <math.h>
+#include <TinyGPS++.h>
+
+#include "DayTime.hpp"
+
+// Constants for sidereal calculation       
+// Source: http://www.stargazing.net/kepler/altaz.html
+#define C1                 100.46
+#define C2                 0.985647
+#define C3                 15
+#define REFERENCE_JEAR     2000
+
+
+class Sidereal
+{
+ public:
+    static DayTime calculateByGPS(TinyGPSPlus* gps);
+
+ private:
+    static const float calculateTheta(float deltaJ, float longitude, float timeUTC);
+    static const int calculateDeltaJd(int year, int month, int day);
+};