mirror of
https://github.com/f4exb/sdrangel.git
synced 2026-06-03 06:24:48 -04:00
Star Tracker updates.
Plot Sun and Moon on sky temperature chart. Plot markers on Galactic line-of-sight chart. Create animations from Galactic line-of-sight chart. Allow weather at antenna location to be downloaded from openweathermap.org Allow target to be entered as Galactic longitude / latitude. Add azimuth and elevation offsets to support scans around targets. Add S7, S8 and S9 targets. Refactor some code from GUI to main plugin, so computed values can be used in other plugins.
This commit is contained in:
Jon Beniston
@@ -25,14 +25,20 @@
|
||||
#include "SWGFeatureReport.h"
|
||||
#include "SWGFeatureActions.h"
|
||||
#include "SWGDeviceState.h"
|
||||
#include "SWGStarTrackerDisplaySettings.h"
|
||||
|
||||
#include "dsp/dspengine.h"
|
||||
#include "util/weather.h"
|
||||
#include "util/units.h"
|
||||
#include "maincore.h"
|
||||
|
||||
#include "startrackerreport.h"
|
||||
#include "startrackerworker.h"
|
||||
#include "startracker.h"
|
||||
|
||||
MESSAGE_CLASS_DEFINITION(StarTracker::MsgConfigureStarTracker, Message)
|
||||
MESSAGE_CLASS_DEFINITION(StarTracker::MsgStartStop, Message)
|
||||
MESSAGE_CLASS_DEFINITION(StarTracker::MsgSetSolarFlux, Message)
|
||||
|
||||
const char* const StarTracker::m_featureIdURI = "sdrangel.feature.startracker";
|
||||
const char* const StarTracker::m_featureId = "StarTracker";
|
||||
@@ -45,8 +51,17 @@ StarTracker::StarTracker(WebAPIAdapterInterface *webAPIAdapterInterface) :
|
||||
m_worker = new StarTrackerWorker(this, webAPIAdapterInterface);
|
||||
m_state = StIdle;
|
||||
m_errorMessage = "StarTracker error";
|
||||
connect(&m_updatePipesTimer, SIGNAL(timeout()), this, SLOT(updatePipes()));
|
||||
m_updatePipesTimer.start(1000);
|
||||
m_networkManager = new QNetworkAccessManager();
|
||||
connect(m_networkManager, SIGNAL(finished(QNetworkReply*)), this, SLOT(networkManagerFinished(QNetworkReply*)));
|
||||
m_weather = nullptr;
|
||||
m_solarFlux = 0.0f;
|
||||
// Unfortunately, can't seem to access resources in static global constructor
|
||||
m_temps.append(new FITS(":/startracker/startracker/150mhz_ra_dec.fits"));
|
||||
m_temps.append(new FITS(":/startracker/startracker/408mhz_ra_dec.fits"));
|
||||
m_temps.append(new FITS(":/startracker/startracker/1420mhz_ra_dec.fits"));
|
||||
m_spectralIndex = new FITS(":/startracker/startracker/408mhz_ra_dec_spectral_index.fits");
|
||||
}
|
||||
|
||||
StarTracker::~StarTracker()
|
||||
@@ -58,6 +73,13 @@ StarTracker::~StarTracker()
|
||||
}
|
||||
|
||||
delete m_worker;
|
||||
if (m_weather)
|
||||
{
|
||||
disconnect(m_weather, &Weather::weatherUpdated, this, &StarTracker::weatherUpdated);
|
||||
delete m_weather;
|
||||
}
|
||||
qDeleteAll(m_temps);
|
||||
delete m_spectralIndex;
|
||||
}
|
||||
|
||||
void StarTracker::start()
|
||||
@@ -71,8 +93,8 @@ void StarTracker::start()
|
||||
m_state = ok ? StRunning : StError;
|
||||
m_thread.start();
|
||||
|
||||
StarTrackerWorker::MsgConfigureStarTrackerWorker *msg = StarTrackerWorker::MsgConfigureStarTrackerWorker::create(m_settings, true);
|
||||
m_worker->getInputMessageQueue()->push(msg);
|
||||
m_worker->getInputMessageQueue()->push(StarTrackerWorker::MsgConfigureStarTrackerWorker::create(m_settings, true));
|
||||
m_worker->getInputMessageQueue()->push(MsgSetSolarFlux::create(m_solarFlux));
|
||||
}
|
||||
|
||||
void StarTracker::stop()
|
||||
@@ -107,12 +129,52 @@ bool StarTracker::handleMessage(const Message& cmd)
|
||||
|
||||
return true;
|
||||
}
|
||||
else if (MsgSetSolarFlux::match(cmd))
|
||||
{
|
||||
MsgSetSolarFlux& msg = (MsgSetSolarFlux&) cmd;
|
||||
m_solarFlux = msg.getFlux();
|
||||
m_worker->getInputMessageQueue()->push(new MsgSetSolarFlux(msg));
|
||||
return true;
|
||||
}
|
||||
else if (MainCore::MsgStarTrackerDisplaySettings::match(cmd))
|
||||
{
|
||||
MainCore::MsgStarTrackerDisplaySettings& settings = (MainCore::MsgStarTrackerDisplaySettings&) cmd;
|
||||
if (m_guiMessageQueue) {
|
||||
m_guiMessageQueue->push(new MainCore::MsgStarTrackerDisplaySettings(settings));
|
||||
}
|
||||
return true;
|
||||
}
|
||||
else if (MainCore::MsgStarTrackerDisplayLoSSettings::match(cmd))
|
||||
{
|
||||
MainCore::MsgStarTrackerDisplayLoSSettings& settings = (MainCore::MsgStarTrackerDisplayLoSSettings&) cmd;
|
||||
if (m_guiMessageQueue) {
|
||||
m_guiMessageQueue->push(new MainCore::MsgStarTrackerDisplayLoSSettings(settings));
|
||||
}
|
||||
return true;
|
||||
}
|
||||
else
|
||||
{
|
||||
return false;
|
||||
}
|
||||
}
|
||||
|
||||
void StarTracker::updatePipes()
|
||||
{
|
||||
QList<AvailablePipeSource> availablePipes = updateAvailablePipeSources("startracker.display", StarTrackerSettings::m_pipeTypes, StarTrackerSettings::m_pipeURIs, this);
|
||||
|
||||
if (availablePipes != m_availablePipes)
|
||||
{
|
||||
m_availablePipes = availablePipes;
|
||||
if (getMessageQueueToGUI())
|
||||
{
|
||||
MsgReportPipes *msgToGUI = MsgReportPipes::create();
|
||||
QList<AvailablePipeSource>& msgAvailablePipes = msgToGUI->getAvailablePipes();
|
||||
msgAvailablePipes.append(availablePipes);
|
||||
getMessageQueueToGUI()->push(msgToGUI);
|
||||
}
|
||||
}
|
||||
}
|
||||
|
||||
QByteArray StarTracker::serialize() const
|
||||
{
|
||||
return m_settings.serialize();
|
||||
@@ -143,6 +205,10 @@ void StarTracker::applySettings(const StarTrackerSettings& settings, bool force)
|
||||
<< " m_dec: " << settings.m_dec
|
||||
<< " m_az: " << settings.m_az
|
||||
<< " m_el: " << settings.m_el
|
||||
<< " m_l: " << settings.m_l
|
||||
<< " m_b: " << settings.m_b
|
||||
<< " m_azOffset: " << settings.m_azOffset
|
||||
<< " m_elOffset: " << settings.m_elOffset
|
||||
<< " m_latitude: " << settings.m_latitude
|
||||
<< " m_longitude: " << settings.m_longitude
|
||||
<< " m_serverPort: " << settings.m_serverPort
|
||||
@@ -215,6 +281,52 @@ void StarTracker::applySettings(const StarTrackerSettings& settings, bool force)
|
||||
if ((m_settings.m_rgbColor != settings.m_rgbColor) || force) {
|
||||
reverseAPIKeys.append("rgbColor");
|
||||
}
|
||||
if ((m_settings.m_az != settings.m_az) || force) {
|
||||
reverseAPIKeys.append("azimuth");
|
||||
}
|
||||
if ((m_settings.m_el != settings.m_el) || force) {
|
||||
reverseAPIKeys.append("elevation");
|
||||
}
|
||||
if ((m_settings.m_l != settings.m_l) || force) {
|
||||
reverseAPIKeys.append("l");
|
||||
}
|
||||
if ((m_settings.m_b != settings.m_b) || force) {
|
||||
reverseAPIKeys.append("b");
|
||||
}
|
||||
if ((m_settings.m_azOffset != settings.m_azOffset) || force) {
|
||||
reverseAPIKeys.append("azimuthOffset");
|
||||
}
|
||||
if ((m_settings.m_elOffset != settings.m_elOffset) || force) {
|
||||
reverseAPIKeys.append("elevationOffset");
|
||||
}
|
||||
|
||||
if ((m_settings.m_owmAPIKey != settings.m_owmAPIKey) || force)
|
||||
{
|
||||
if (m_weather)
|
||||
{
|
||||
disconnect(m_weather, &Weather::weatherUpdated, this, &StarTracker::weatherUpdated);
|
||||
delete m_weather;
|
||||
m_weather = nullptr;
|
||||
}
|
||||
if (!settings.m_owmAPIKey.isEmpty())
|
||||
{
|
||||
m_weather = Weather::create(settings.m_owmAPIKey);
|
||||
if (m_weather) {
|
||||
connect(m_weather, &Weather::weatherUpdated, this, &StarTracker::weatherUpdated);
|
||||
}
|
||||
}
|
||||
}
|
||||
|
||||
if ( (m_settings.m_owmAPIKey != settings.m_owmAPIKey)
|
||||
|| (m_settings.m_latitude != settings.m_latitude)
|
||||
|| (m_settings.m_longitude != settings.m_longitude)
|
||||
|| (m_settings.m_weatherUpdatePeriod != settings.m_weatherUpdatePeriod)
|
||||
|| force)
|
||||
{
|
||||
if (m_weather) {
|
||||
m_weather->getWeatherPeriodically(m_settings.m_latitude, m_settings.m_longitude, settings.m_weatherUpdatePeriod);
|
||||
}
|
||||
}
|
||||
|
||||
StarTrackerWorker::MsgConfigureStarTrackerWorker *msg = StarTrackerWorker::MsgConfigureStarTrackerWorker::create(
|
||||
settings, force
|
||||
@@ -320,6 +432,13 @@ void StarTracker::webapiFormatFeatureSettings(
|
||||
response.getStarTrackerSettings()->setReverseApiPort(settings.m_reverseAPIPort);
|
||||
response.getStarTrackerSettings()->setReverseApiFeatureSetIndex(settings.m_reverseAPIFeatureSetIndex);
|
||||
response.getStarTrackerSettings()->setReverseApiFeatureIndex(settings.m_reverseAPIFeatureIndex);
|
||||
|
||||
response.getStarTrackerSettings()->setAzimuth(settings.m_az);
|
||||
response.getStarTrackerSettings()->setElevation(settings.m_el);
|
||||
response.getStarTrackerSettings()->setL(settings.m_l);
|
||||
response.getStarTrackerSettings()->setB(settings.m_b);
|
||||
response.getStarTrackerSettings()->setAzimuthOffset(settings.m_azOffset);
|
||||
response.getStarTrackerSettings()->setElevationOffset(settings.m_elOffset);
|
||||
}
|
||||
|
||||
void StarTracker::webapiUpdateFeatureSettings(
|
||||
@@ -396,6 +515,24 @@ void StarTracker::webapiUpdateFeatureSettings(
|
||||
if (featureSettingsKeys.contains("reverseAPIFeatureIndex")) {
|
||||
settings.m_reverseAPIFeatureIndex = response.getStarTrackerSettings()->getReverseApiFeatureIndex();
|
||||
}
|
||||
if (featureSettingsKeys.contains("azimuth")) {
|
||||
settings.m_az = response.getStarTrackerSettings()->getAzimuth();
|
||||
}
|
||||
if (featureSettingsKeys.contains("elevation")) {
|
||||
settings.m_el = response.getStarTrackerSettings()->getElevation();
|
||||
}
|
||||
if (featureSettingsKeys.contains("l")) {
|
||||
settings.m_l = response.getStarTrackerSettings()->getL();
|
||||
}
|
||||
if (featureSettingsKeys.contains("b")) {
|
||||
settings.m_b = response.getStarTrackerSettings()->getB();
|
||||
}
|
||||
if (featureSettingsKeys.contains("azimuthOffset")) {
|
||||
settings.m_azOffset = response.getStarTrackerSettings()->getAzimuthOffset();
|
||||
}
|
||||
if (featureSettingsKeys.contains("elevationOffset")) {
|
||||
settings.m_elOffset = response.getStarTrackerSettings()->getElevationOffset();
|
||||
}
|
||||
}
|
||||
|
||||
void StarTracker::webapiReverseSendSettings(QList<QString>& featureSettingsKeys, const StarTrackerSettings& settings, bool force)
|
||||
@@ -463,6 +600,24 @@ void StarTracker::webapiReverseSendSettings(QList<QString>& featureSettingsKeys,
|
||||
if (featureSettingsKeys.contains("rgbColor") || force) {
|
||||
swgStarTrackerSettings->setRgbColor(settings.m_rgbColor);
|
||||
}
|
||||
if (featureSettingsKeys.contains("azimuth") || force) {
|
||||
swgStarTrackerSettings->setAzimuth(settings.m_az);
|
||||
}
|
||||
if (featureSettingsKeys.contains("elevation") || force) {
|
||||
swgStarTrackerSettings->setElevation(settings.m_el);
|
||||
}
|
||||
if (featureSettingsKeys.contains("l") || force) {
|
||||
swgStarTrackerSettings->setL(settings.m_l);
|
||||
}
|
||||
if (featureSettingsKeys.contains("b") || force) {
|
||||
swgStarTrackerSettings->setB(settings.m_b);
|
||||
}
|
||||
if (featureSettingsKeys.contains("azimuthOffset") || force) {
|
||||
swgStarTrackerSettings->setAzimuthOffset(settings.m_azOffset);
|
||||
}
|
||||
if (featureSettingsKeys.contains("elevationOffset") || force) {
|
||||
swgStarTrackerSettings->setElevationOffset(settings.m_elOffset);
|
||||
}
|
||||
|
||||
QString channelSettingsURL = QString("http://%1:%2/sdrangel/featureset/%3/feature/%4/settings")
|
||||
.arg(settings.m_reverseAPIAddress)
|
||||
@@ -504,3 +659,172 @@ void StarTracker::networkManagerFinished(QNetworkReply *reply)
|
||||
|
||||
reply->deleteLater();
|
||||
}
|
||||
|
||||
void StarTracker::weatherUpdated(float temperature, float pressure, float humidity)
|
||||
{
|
||||
if (!isnan(temperature)) {
|
||||
m_settings.m_temperature = temperature;
|
||||
}
|
||||
if (!isnan(pressure)) {
|
||||
m_settings.m_pressure = pressure;
|
||||
}
|
||||
if (!isnan(humidity)) {
|
||||
m_settings.m_humidity = humidity;
|
||||
}
|
||||
|
||||
m_worker->getInputMessageQueue()->push(StarTrackerWorker::MsgConfigureStarTrackerWorker::create(m_settings, false));
|
||||
if (m_guiMessageQueue) {
|
||||
m_guiMessageQueue->push(MsgConfigureStarTracker::create(m_settings, false));
|
||||
}
|
||||
}
|
||||
|
||||
double StarTracker::applyBeam(const FITS *fits, double beamwidth, double ra, double dec, int& imgX, int& imgY) const
|
||||
{
|
||||
const double halfBeamwidth = beamwidth/2.0;
|
||||
// Use cos^p(x) for approximation of radiation pattern
|
||||
// (Essentially the same as Gaussian of exp(-4*ln(theta^2/beamwidth^2))
|
||||
// (See a2 in https://arxiv.org/pdf/1812.10084.pdf for Elliptical equivalent))
|
||||
// We have gain of 0dB (1) at 0 degrees, and -3dB (~0.5) at half-beamwidth degrees
|
||||
// Find exponent that correponds to -3dB at that angle
|
||||
double minus3dBLinear = pow(10.0, -3.0/10.0);
|
||||
double p = log(minus3dBLinear)/log(cos(Units::degreesToRadians(halfBeamwidth)));
|
||||
// Create an matrix with gain as a function of angle
|
||||
double degreesPerPixelH = abs(fits->degreesPerPixelH());
|
||||
double degreesPerPixelV = abs(fits->degreesPerPixelV());
|
||||
int numberOfCoeffsH = ceil(beamwidth/degreesPerPixelH);
|
||||
int numberOfCoeffsV = ceil(beamwidth/degreesPerPixelV);
|
||||
if ((numberOfCoeffsH & 1) == 0) {
|
||||
numberOfCoeffsH++;
|
||||
}
|
||||
if ((numberOfCoeffsV & 1) == 0) {
|
||||
numberOfCoeffsV++;
|
||||
}
|
||||
double *beam = new double[numberOfCoeffsH*numberOfCoeffsV];
|
||||
double sum = 0.0;
|
||||
int y0 = numberOfCoeffsV/2;
|
||||
int x0 = numberOfCoeffsH/2;
|
||||
int nonZeroCount = 0;
|
||||
for (int y = 0; y < numberOfCoeffsV; y++)
|
||||
{
|
||||
for (int x = 0; x < numberOfCoeffsH; x++)
|
||||
{
|
||||
double xp = (x - x0) * degreesPerPixelH;
|
||||
double yp = (y - y0) * degreesPerPixelV;
|
||||
double r = sqrt(xp*xp+yp*yp);
|
||||
if (r < halfBeamwidth)
|
||||
{
|
||||
beam[y*numberOfCoeffsH+x] = pow(cos(Units::degreesToRadians(r)), p);
|
||||
sum += beam[y*numberOfCoeffsH+x];
|
||||
nonZeroCount++;
|
||||
}
|
||||
else
|
||||
{
|
||||
beam[y*numberOfCoeffsH+x] = 0.0;
|
||||
}
|
||||
}
|
||||
}
|
||||
|
||||
// Get centre pixel coordinates
|
||||
double centreX;
|
||||
if (ra <= 12.0) {
|
||||
centreX = (12.0 - ra) / 24.0;
|
||||
} else {
|
||||
centreX = (24 - ra + 12) / 24.0;
|
||||
}
|
||||
double centreY = (90.0-dec) / 180.0;
|
||||
imgX = centreX * fits->width();
|
||||
imgY = centreY * fits->height();
|
||||
|
||||
// Apply weighting to temperature data
|
||||
double weightedSum = 0.0;
|
||||
for (int y = 0; y < numberOfCoeffsV; y++)
|
||||
{
|
||||
for (int x = 0; x < numberOfCoeffsH; x++)
|
||||
{
|
||||
weightedSum += beam[y*numberOfCoeffsH+x] * fits->scaledWrappedValue(imgX + (x-x0), imgY + (y-y0));
|
||||
}
|
||||
}
|
||||
// From: https://www.cv.nrao.edu/~sransom/web/Ch3.html
|
||||
// The antenna temperature equals the source brightness temperature multiplied by the fraction of the beam solid angle filled by the source
|
||||
// So we scale the sum by the total number of non-zero pixels (i.e. beam area)
|
||||
// If we compare to some maps with different beamwidths here: https://www.cv.nrao.edu/~demerson/radiosky/sky_jun96.pdf
|
||||
// The values we've computed are a bit higher..
|
||||
double temp = weightedSum/nonZeroCount;
|
||||
|
||||
delete[] beam;
|
||||
|
||||
return temp;
|
||||
}
|
||||
|
||||
bool StarTracker::calcSkyTemperature(double frequency, double beamwidth, double ra, double dec, double& temp) const
|
||||
{
|
||||
const FITS *fits;
|
||||
int imgX, imgY;
|
||||
|
||||
if ((frequency >= 1.4e9) && (frequency <= 1.45e9))
|
||||
{
|
||||
// Adjust temperature from 1420MHz FITS file, just using beamwidth
|
||||
fits = getTempFITS(2);
|
||||
if (fits && fits->valid())
|
||||
{
|
||||
temp = applyBeam(fits, beamwidth, ra, dec, imgX, imgY);
|
||||
return true;
|
||||
}
|
||||
else
|
||||
{
|
||||
qDebug() << "StarTracker::calcSkyTemperature: 1420MHz FITS temperature file not valid";
|
||||
return false;
|
||||
}
|
||||
}
|
||||
else
|
||||
{
|
||||
// Adjust temperature from 408MHz FITS file, taking in to account
|
||||
// observation frequency and beamwidth
|
||||
fits = getTempFITS(1);
|
||||
if (fits && fits->valid())
|
||||
{
|
||||
double temp408 = applyBeam(fits, beamwidth, ra, dec, imgX, imgY);
|
||||
|
||||
// Scale according to frequency - CMB contribution constant
|
||||
// Power law at low frequencies, with slight variation in spectral index
|
||||
// See:
|
||||
// Global Sky Model: https://ascl.net/1011.010
|
||||
// An improved Model of Diffuse Galactic Radio Emission: https://arxiv.org/pdf/1605.04920.pdf
|
||||
// A high-resolution self-consistent whole sky foreground model: https://arxiv.org/abs/1812.10084
|
||||
// (De-striping:) Full sky study of diffuse Galactic emission at decimeter wavelength https://www.aanda.org/articles/aa/pdf/2003/42/aah4363.pdf
|
||||
// Data here: http://cdsarc.u-strasbg.fr/viz-bin/cat/J/A+A/410/847
|
||||
// LFmap: https://www.faculty.ece.vt.edu/swe/lwa/memo/lwa0111.pdf
|
||||
double iso408 = 50 * pow(150e6/408e6, 2.75); // Extra-galactic isotropic in reference map at 408MHz
|
||||
double isoT = 50 * pow(150e6/frequency, 2.75); // Extra-galactic isotropic at target frequency
|
||||
double cmbT = 2.725; // Cosmic microwave backgroud;
|
||||
double spectralIndex;
|
||||
const FITS *spectralIndexFITS = getSpectralIndexFITS();
|
||||
if (spectralIndexFITS && spectralIndexFITS->valid())
|
||||
{
|
||||
// See https://www.aanda.org/articles/aa/pdf/2003/42/aah4363.pdf
|
||||
spectralIndex = spectralIndexFITS->scaledValue(imgX, imgY);
|
||||
}
|
||||
else
|
||||
{
|
||||
// See https://arxiv.org/abs/1812.10084 fig 2
|
||||
if (frequency < 200e6) {
|
||||
spectralIndex = 2.55;
|
||||
} else if (frequency < 20e9) {
|
||||
spectralIndex = 2.695;
|
||||
} else {
|
||||
spectralIndex = 3.1;
|
||||
}
|
||||
}
|
||||
double galactic480 = temp408 - cmbT - iso408;
|
||||
double galacticT = galactic480 * pow(408e6/frequency, spectralIndex); // Scale galactic contribution by frequency
|
||||
temp = galacticT + cmbT + isoT; // Final temperature
|
||||
|
||||
return true;
|
||||
}
|
||||
else
|
||||
{
|
||||
qDebug() << "StarTracker::calcSkyTemperature: 408MHz FITS temperature file not valid";
|
||||
return false;
|
||||
}
|
||||
}
|
||||
}
|
||||
|
||||
Reference in New Issue
Block a user