homed-service-common/sun.cpp at master · u236/homed-service-common · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
#include "sun.h"

void Sun::updateSunrise(void)
{
    double t = julianCentury(m_julianDay), n = julianCentury(julianDay(t) + (720 - (m_longitude + degree(sunriseHourAngle(m_latitude, sunDeclination(t), ANGLE))) * 4 - timeEquation(t)) / 1440.0), m = round(720 - (m_longitude + degree(sunriseHourAngle(m_latitude, sunDeclination(n), ANGLE))) * 4 - timeEquation(n)) + m_offset;
    m_sunrise = QTime(static_cast <int> (m / 60), static_cast <int> (m) % 60);
}

void Sun::updateSunset(void)
{
    double t = julianCentury(m_julianDay), n = julianCentury(julianDay(t) + (720 - (m_longitude + degree(sunsetHourAngle(m_latitude, sunDeclination(t), ANGLE))) * 4 - timeEquation(t)) / 1440.0), m = round(720 - (m_longitude + degree(sunsetHourAngle(m_latitude, sunDeclination(n), ANGLE))) * 4 - timeEquation(n)) + m_offset;
    m_sunset = QTime(static_cast <int> (m / 60), static_cast <int> (m) % 60);
}

QTime Sun::fromString(const QString &string)
{
    QList <QString> itemList = string.split(QRegExp("[(\\-|\\+)]")), valueList = {"sunrise", "sunset"};
    QString value = itemList.value(0).toLower().trimmed();
    qint32 offset = itemList.value(1).toInt();

    if (string.mid(itemList.value(0).length(), 1) == "-")
        offset *= -1;

    switch (valueList.indexOf(value))
    {
        case 0:  return m_sunrise.addSecs(offset * 60);
        case 1:  return m_sunset.addSecs(offset * 60);
        default: return QTime::fromString(value, "h:mm");
    }
}

double Sun::position(void)
{
    double sunrise = m_sunrise.msecsSinceStartOfDay(), sunset = m_sunset.msecsSinceStartOfDay();
    return 1 - sin(M_PI * (QDateTime::currentDateTime().time().msecsSinceStartOfDay() - sunrise) / (sunset - sunrise));
}

double Sun::julianDay(double century)
{
    return century * 36525 + 2451545;
}

double Sun::julianCentury(double day)
{
    return (day - 2451545) / 36525;
}

double Sun::meanSunAnomaly(double t)
{
    return 357.52911 + t * (35999.05029 - t * 0.0001537);
}

double Sun::meanSunLongitude(double t)
{
    return fmod(t * (36000.76983 + t * 0.0003032) + 280.46646, 360);
}

double Sun::sunDeclination(double t)
{
    double r = radian(meanSunAnomaly(t)), l = meanSunLongitude(t) + sin(r) * (1.914602 - t * (0.004817 + t * 0.000014)) + sin(r * 2) * (0.019993 - t * 0.000101) + sin(r * 3) * 0.000289 - sin(radian(125.04 - 1934.136 * t)) * 0.00478 - 0.00569;
    return degree(asin(sin(radian(obliquityCorrection(t))) * sin(radian(l))));
}

double Sun::obliquityCorrection(double t)
{
    return cos(radian(125.04 - t * 1934.136)) * 0.00256 + (((21.448 - t * (46.8150 + t * (0.00059 - t * 0.001813))) / 60) + 26) / 60 + 23;
}

double Sun::timeEquation(double t)
{
    double e = 0.016708634 - t * (0.000042037 + t * 0.0000001267), y = pow(tan(radian(obliquityCorrection(t)) / 2), 2), a = meanSunAnomaly(t), l = meanSunLongitude(t), s = sin(radian(a));
    return degree(y * sin(radian(l) * 2) - e * s * 2 + e * y * s * cos(radian(l) * 2) * 4 - pow(y, 2) * sin(radian(l) * 4) / 2 - pow(e, 2) * sin(radian(a) * 2) * 1.25) * 4;
}

double Sun::sunriseHourAngle(double latitude, double declination, double offset)
{
    double l = radian(latitude), s = radian(declination);
    return acos(cos(radian(offset)) / (cos(l) * cos(s)) - tan(l) * tan(s));
}

double Sun::sunsetHourAngle(double latitude, double declination, double offset)
{
    return sunriseHourAngle(latitude, declination, offset) * -1;
}