spacecruft
/
satnogs-network
1
0
Fork 0
Code Releases Activity

Refactor polar plot in observation view

environments/stage/deployments/136
Fabian P. Schmidt 2018-05-10 15:29:58 +02:00
parent c1f56da5af
commit 1384a447b3
2 changed files with 86 additions and 51 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

61
network/static/js/observation_view.js
View File

@ -1,4 +1,4 @@
/* global WaveSurfer URI satellite moment */
/* global WaveSurfer URI calcPolarPlotSVG */
$(document).ready(function() {
'use strict';
@ -117,59 +117,18 @@ $(document).ready(function() {
end: new Date($('svg#polar').data('timeframe-end'))
};
const pi = Math.PI;
const deg2rad = pi / 180.0;
const rad2deg = 180 / pi;
// Get the observer at lat/lon in RADIANS, altitude in km above ground (NOTE: BUG, should be elevation?)
var observerGd = {
longitude: $('svg#polar').data('groundstation-lon') * deg2rad,
latitude: $('svg#polar').data('groundstation-lat') * deg2rad,
height: $('svg#polar').data('groundstation-alt') / 1000
var groundstation = {
lon: $('svg#polar').data('groundstation-lon'),
lat: $('svg#polar').data('groundstation-lat'),
alt: $('svg#polar').data('groundstation-alt')
};
var polarGetXY = function(az, el) {
var ret = new Object();
ret.x = (90 - el) * Math.sin(az * deg2rad);
ret.y = (el - 90) * Math.cos(az * deg2rad);
return ret;
};
const polarPlotSVG = calcPolarPlotSVG(timeframe,
groundstation,
tleLine1,
tleLine2);
var polarOrbit = document.createElementNS('http://www.w3.org/2000/svg', 'path');
polarOrbit.setAttributeNS(null, 'fill', 'none');
polarOrbit.setAttributeNS(null, 'stroke', 'blue');
polarOrbit.setAttributeNS(null, 'stroke-opacity', '1.0');
polarOrbit.setAttributeNS(null, 'stroke-width', '2');
// Initialize the satellite record
var satrec = satellite.twoline2satrec(tleLine1, tleLine2);
// Draw the orbit pass on the polar az/el plot
var g = '';
for (var t = moment(timeframe.start); t < moment(timeframe.end); t.add(20, 's')) {
var positionAndVelocity = satellite.propagate(satrec, t.toDate());
var gmst = satellite.gstime(t.toDate());
var positionEci = positionAndVelocity.position;
var positionEcf = satellite.eciToEcf(positionEci, gmst);
var lookAngles = satellite.ecfToLookAngles(observerGd, positionEcf);
var azimuth = lookAngles.azimuth * rad2deg,
elevation = lookAngles.elevation * rad2deg;
var coord = polarGetXY(azimuth, elevation);
if (g == '') {
// Start of line
g += 'M';
} else {
// Continue line
g += ' L';
}
g += coord.x + ' ' + coord.y;
}
polarOrbit.setAttribute('d', g);
$('svg#polar').append(polarOrbit);
$('svg#polar').append(polarPlotSVG);
// Hotkeys bindings
$(document).bind('keyup', function(event){

76
network/static/js/polar_svg.js 100644
View File

@ -0,0 +1,76 @@
/* global satellite moment */
/* exported calcPolarPlotSVG */
/**
* @typedef {Object} timeframe
* @property {Date} start - Start of the observation.
* @property {Date} end - End of the observation.
*/
/**
* Returns a polar plot of a pass at the given groundstation as SVG in a string.
*
* @param {timeframe} Timeframe of the oberservation.
* @param {groundstation} The observing groundstation.
* @param {tleLine1} TLE line 1 of the observed satellite.
* @param {tleLine2} TLE line 2 of the observed satellite.
*/
function calcPolarPlotSVG(timeframe, groundstation, tleLine1, tleLine2) {
'use strict';
const pi = Math.PI;
const deg2rad = pi / 180.0;
const rad2deg = 180 / pi;
// Get the observer at lat/lon in RADIANS, altitude in km above ground (NOTE: BUG, should be elevation?)
var observerGd = {
longitude: groundstation.lon * deg2rad,
latitude: groundstation.lat * deg2rad,
height: groundstation.alt / 1000
};
var polarGetXY = function(az, el) {
var ret = new Object();
ret.x = (90 - el) * Math.sin(az * deg2rad);
ret.y = (el - 90) * Math.cos(az * deg2rad);
return ret;
};
var svg_namespace = 'http://www.w3.org/2000/svg';
var polarOrbit = document.createElementNS(svg_namespace, 'path');
polarOrbit.setAttributeNS(null, 'fill', 'none');
polarOrbit.setAttributeNS(null, 'stroke', 'blue');
polarOrbit.setAttributeNS(null, 'stroke-opacity', '1.0');
polarOrbit.setAttributeNS(null, 'stroke-width', '2');
// Initialize the satellite record
var satrec = satellite.twoline2satrec(tleLine1, tleLine2);
// Draw the orbit pass on the polar az/el plot
var g = '';
for (var t = moment(timeframe.start); t < moment(timeframe.end); t.add(20, 's')) {
var positionAndVelocity = satellite.propagate(satrec, t.toDate());
var gmst = satellite.gstime(t.toDate());
var positionEci = positionAndVelocity.position;
var positionEcf = satellite.eciToEcf(positionEci, gmst);
var lookAngles = satellite.ecfToLookAngles(observerGd, positionEcf);
var azimuth = lookAngles.azimuth * rad2deg,
elevation = lookAngles.elevation * rad2deg;
var coord = polarGetXY(azimuth, elevation);
if (g == '') {
// Start of line
g += 'M';
} else {
// Continue line
g += ' L';
}
g += coord.x + ' ' + coord.y;
}
polarOrbit.setAttribute('d', g);
return polarOrbit:
}