From 9ec59b26d49bb693aac1a33cd87399bd3d1a0c99 Mon Sep 17 00:00:00 2001
From: Andreas <6977712+AndreasK79@users.noreply.github.com>
Date: Wed, 14 Feb 2024 10:50:01 +0100
Subject: [PATCH] [LBA] Add night shadow layer
---
assets/js/leaflet/L.Terminator.js | 148 ++++++++++++++++++++++++++
assets/js/sections/logbookadvanced.js | 16 ++-
2 files changed, 161 insertions(+), 3 deletions(-)
create mode 100644 assets/js/leaflet/L.Terminator.js
diff --git a/assets/js/leaflet/L.Terminator.js b/assets/js/leaflet/L.Terminator.js
new file mode 100644
index 000000000..8e6021061
--- /dev/null
+++ b/assets/js/leaflet/L.Terminator.js
@@ -0,0 +1,148 @@
+(function (global, factory) {
+ typeof exports === 'object' && typeof module !== 'undefined' ? module.exports = factory(require('leaflet')) :
+ typeof define === 'function' && define.amd ? define(['leaflet'], factory) :
+ (global.L = global.L || {}, global.L.terminator = factory(global.L));
+}(this, (function (L) { 'use strict';
+
+ L = L && L.hasOwnProperty('default') ? L['default'] : L;
+
+ /* Terminator.js -- Overlay day/night region on a Leaflet map */
+
+ function julian(date) {
+ /* Calculate the present UTC Julian Date. Function is valid after
+ * the beginning of the UNIX epoch 1970-01-01 and ignores leap
+ * seconds. */
+ return (date / 86400000) + 2440587.5;
+ }
+
+ function GMST(julianDay) {
+ /* Calculate Greenwich Mean Sidereal Time according to
+ http://aa.usno.navy.mil/faq/docs/GAST.php */
+ var d = julianDay - 2451545.0;
+ // Low precision equation is good enough for our purposes.
+ return (18.697374558 + 24.06570982441908 * d) % 24;
+ }
+
+ var Terminator = L.Polygon.extend({
+ options: {
+ color: '#00',
+ opacity: 0.5,
+ fillColor: '#00',
+ fillOpacity: 0.5,
+ resolution: 2
+ },
+
+ initialize: function (options) {
+ this.version = '0.1.0';
+ this._R2D = 180 / Math.PI;
+ this._D2R = Math.PI / 180;
+ L.Util.setOptions(this, options);
+ var latLng = this._compute(this.options.time);
+ this.setLatLngs(latLng);
+ },
+
+ setTime: function (date) {
+ this.options.time = date;
+ var latLng = this._compute(date);
+ this.setLatLngs(latLng);
+ },
+
+ _sunEclipticPosition: function (julianDay) {
+ /* Compute the position of the Sun in ecliptic coordinates at
+ julianDay. Following
+ http://en.wikipedia.org/wiki/Position_of_the_Sun */
+ // Days since start of J2000.0
+ var n = julianDay - 2451545.0;
+ // mean longitude of the Sun
+ var L$$1 = 280.460 + 0.9856474 * n;
+ L$$1 %= 360;
+ // mean anomaly of the Sun
+ var g = 357.528 + 0.9856003 * n;
+ g %= 360;
+ // ecliptic longitude of Sun
+ var lambda = L$$1 + 1.915 * Math.sin(g * this._D2R) +
+ 0.02 * Math.sin(2 * g * this._D2R);
+ // distance from Sun in AU
+ var R = 1.00014 - 0.01671 * Math.cos(g * this._D2R) -
+ 0.0014 * Math.cos(2 * g * this._D2R);
+ return {lambda: lambda, R: R};
+ },
+
+ _eclipticObliquity: function (julianDay) {
+ // Following the short term expression in
+ // http://en.wikipedia.org/wiki/Axial_tilt#Obliquity_of_the_ecliptic_.28Earth.27s_axial_tilt.29
+ var n = julianDay - 2451545.0;
+ // Julian centuries since J2000.0
+ var T = n / 36525;
+ var epsilon = 23.43929111 -
+ T * (46.836769 / 3600
+ - T * (0.0001831 / 3600
+ + T * (0.00200340 / 3600
+ - T * (0.576e-6 / 3600
+ - T * 4.34e-8 / 3600))));
+ return epsilon;
+ },
+
+ _sunEquatorialPosition: function (sunEclLng, eclObliq) {
+ /* Compute the Sun's equatorial position from its ecliptic
+ * position. Inputs are expected in degrees. Outputs are in
+ * degrees as well. */
+ var alpha = Math.atan(Math.cos(eclObliq * this._D2R)
+ * Math.tan(sunEclLng * this._D2R)) * this._R2D;
+ var delta = Math.asin(Math.sin(eclObliq * this._D2R)
+ * Math.sin(sunEclLng * this._D2R)) * this._R2D;
+
+ var lQuadrant = Math.floor(sunEclLng / 90) * 90;
+ var raQuadrant = Math.floor(alpha / 90) * 90;
+ alpha = alpha + (lQuadrant - raQuadrant);
+
+ return {alpha: alpha, delta: delta};
+ },
+
+ _hourAngle: function (lng, sunPos, gst) {
+ /* Compute the hour angle of the sun for a longitude on
+ * Earth. Return the hour angle in degrees. */
+ var lst = gst + lng / 15;
+ return lst * 15 - sunPos.alpha;
+ },
+
+ _latitude: function (ha, sunPos) {
+ /* For a given hour angle and sun position, compute the
+ * latitude of the terminator in degrees. */
+ var lat = Math.atan(-Math.cos(ha * this._D2R) /
+ Math.tan(sunPos.delta * this._D2R)) * this._R2D;
+ return lat;
+ },
+
+ _compute: function (time) {
+ var today = time ? new Date(time) : new Date();
+ var julianDay = julian(today);
+ var gst = GMST(julianDay);
+ var latLng = [];
+
+ var sunEclPos = this._sunEclipticPosition(julianDay);
+ var eclObliq = this._eclipticObliquity(julianDay);
+ var sunEqPos = this._sunEquatorialPosition(sunEclPos.lambda, eclObliq);
+ for (var i = 0; i <= 720 * this.options.resolution; i++) {
+ var lng = -360 + i / this.options.resolution;
+ var ha = this._hourAngle(lng, sunEqPos, gst);
+ latLng[i + 1] = [this._latitude(ha, sunEqPos), lng];
+ }
+ if (sunEqPos.delta < 0) {
+ latLng[0] = [90, -360];
+ latLng[latLng.length] = [90, 360];
+ } else {
+ latLng[0] = [-90, -360];
+ latLng[latLng.length] = [-90, 360];
+ }
+ return latLng;
+ }
+ });
+
+ function terminator(options) {
+ return new Terminator(options);
+ }
+
+ return terminator;
+
+})));
diff --git a/assets/js/sections/logbookadvanced.js b/assets/js/sections/logbookadvanced.js
index 113c872af..b346af9d1 100644
--- a/assets/js/sections/logbookadvanced.js
+++ b/assets/js/sections/logbookadvanced.js
@@ -8,6 +8,7 @@ var geojson;
var itugeojson;
var zonemarkers = [];
var ituzonemarkers = [];
+var nightlayer;
var defaultlinecolor = 'blue';
@@ -1065,7 +1066,8 @@ function loadMap(data, iconsList) {
div.innerHTML += ' Path lines
';
div.innerHTML += ' Gridsquares
';
div.innerHTML += ' CQ Zones
';
- div.innerHTML += ' ITU Zones';
+ div.innerHTML += ' ITU Zones
';
+ div.innerHTML += ' Show night shadow';
return div;
};
@@ -1086,9 +1088,9 @@ function loadMap(data, iconsList) {
filename: 'Wavelog',
exportOnly: true,
hideControlContainer: true
- }).addTo(map);
+ }).addTo(map);
- map.on('mousemove', onMapMove);
+ map.on('mousemove', onMapMove);
}
function createContentMessage(qso) {
@@ -1389,6 +1391,14 @@ function loadMap(data, iconsList) {
}
}
+ function toggleNightShadow(bool) {
+ if(!bool) {
+ map.removeLayer(nightlayer);
+ } else {
+ nightlayer = L.terminator().addTo(map);
+ }
+ }
+
function style(feature) {
var bordercolor = "black";
if (isDarkModeTheme()) {