mirror of
https://github.com/f4exb/sdrangel.git
synced 2024-08-26 07:46:33 -04:00
124 lines
4.1 KiB
C
124 lines
4.1 KiB
C
|
///////////////////////////////////////////////////////////////////////////////////
|
||
|
// Copyright (C) 2019 Edouard Griffiths, F4EXB //
|
||
|
// Copyright (C) 2021 Jon Beniston, M7RCE //
|
||
|
// //
|
||
|
// This program is free software; you can redistribute it and/or modify //
|
||
|
// it under the terms of the GNU General Public License as published by //
|
||
|
// the Free Software Foundation as version 3 of the License, or //
|
||
|
// (at your option) any later version. //
|
||
|
// //
|
||
|
// This program is distributed in the hope that it will be useful, //
|
||
|
// but WITHOUT ANY WARRANTY; without even the implied warranty of //
|
||
|
// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the //
|
||
|
// GNU General Public License V3 for more details. //
|
||
|
// //
|
||
|
// You should have received a copy of the GNU General Public License //
|
||
|
// along with this program. If not, see <http://www.gnu.org/licenses/>. //
|
||
|
///////////////////////////////////////////////////////////////////////////////////
|
||
|
|
||
|
#ifndef INCLUDE_RADIOASTRONOMYSINK_H
|
||
|
#define INCLUDE_RADIOASTRONOMYSINK_H
|
||
|
|
||
|
#include "dsp/channelsamplesink.h"
|
||
|
#include "dsp/nco.h"
|
||
|
#include "dsp/interpolator.h"
|
||
|
#include "dsp/fftfactory.h"
|
||
|
#include "dsp/fftengine.h"
|
||
|
#include "dsp/fftwindow.h"
|
||
|
#include "util/movingaverage.h"
|
||
|
#include "util/messagequeue.h"
|
||
|
|
||
|
#include "radioastronomysettings.h"
|
||
|
|
||
|
class ChannelAPI;
|
||
|
class RadioAstronomy;
|
||
|
|
||
|
class RadioAstronomySink : public ChannelSampleSink {
|
||
|
public:
|
||
|
RadioAstronomySink(RadioAstronomy *aisDemod);
|
||
|
~RadioAstronomySink();
|
||
|
|
||
|
virtual void feed(const SampleVector::const_iterator& begin, const SampleVector::const_iterator& end);
|
||
|
|
||
|
void applyChannelSettings(int channelSampleRate, int channelFrequencyOffset, bool force = false);
|
||
|
void applySettings(const RadioAstronomySettings& settings, bool force = false);
|
||
|
void setMessageQueueToChannel(MessageQueue *messageQueue) { m_messageQueueToChannel = messageQueue; }
|
||
|
void setChannel(ChannelAPI *channel) { m_channel = channel; }
|
||
|
void startMeasurements();
|
||
|
void stopMeasurements();
|
||
|
void startCal(bool hot);
|
||
|
void clearCal();
|
||
|
|
||
|
double getMagSq() const { return m_magsq; }
|
||
|
|
||
|
void getMagSqLevels(double& avg, double& peak, int& nbSamples)
|
||
|
{
|
||
|
if (m_magsqCount > 0)
|
||
|
{
|
||
|
m_magsq = m_magsqSum / m_magsqCount;
|
||
|
m_magSqLevelStore.m_magsq = m_magsq;
|
||
|
m_magSqLevelStore.m_magsqPeak = m_magsqPeak;
|
||
|
}
|
||
|
|
||
|
avg = m_magSqLevelStore.m_magsq;
|
||
|
peak = m_magSqLevelStore.m_magsqPeak;
|
||
|
nbSamples = m_magsqCount == 0 ? 1 : m_magsqCount;
|
||
|
|
||
|
m_magsqSum = 0.0f;
|
||
|
m_magsqPeak = 0.0f;
|
||
|
m_magsqCount = 0;
|
||
|
}
|
||
|
|
||
|
|
||
|
private:
|
||
|
struct MagSqLevelsStore
|
||
|
{
|
||
|
MagSqLevelsStore() :
|
||
|
m_magsq(1e-12),
|
||
|
m_magsqPeak(1e-12)
|
||
|
{}
|
||
|
double m_magsq;
|
||
|
double m_magsqPeak;
|
||
|
};
|
||
|
|
||
|
RadioAstronomy *m_radioAstronomy;
|
||
|
RadioAstronomySettings m_settings;
|
||
|
ChannelAPI *m_channel;
|
||
|
int m_channelSampleRate;
|
||
|
int m_channelFrequencyOffset;
|
||
|
|
||
|
int m_fftSequence;
|
||
|
FFTEngine *m_fft;
|
||
|
FFTWindow m_fftWindow;
|
||
|
int m_fftCounter;
|
||
|
QList<int> m_filterBins;
|
||
|
|
||
|
Real *m_fftSum;
|
||
|
Real *m_fftTemp;
|
||
|
int m_fftSumCount;
|
||
|
|
||
|
bool m_enabled;
|
||
|
bool m_cal;
|
||
|
bool m_hot;
|
||
|
|
||
|
NCO m_nco;
|
||
|
Interpolator m_interpolator;
|
||
|
Real m_interpolatorDistance;
|
||
|
Real m_interpolatorDistanceRemain;
|
||
|
|
||
|
double m_magsq;
|
||
|
double m_magsqSum;
|
||
|
double m_magsqPeak;
|
||
|
int m_magsqCount;
|
||
|
MagSqLevelsStore m_magSqLevelStore;
|
||
|
|
||
|
MessageQueue *m_messageQueueToChannel;
|
||
|
|
||
|
MovingAverageUtil<Real, double, 16> m_movingAverage;
|
||
|
|
||
|
void processOneSample(Complex &ci);
|
||
|
MessageQueue *getMessageQueueToChannel() { return m_messageQueueToChannel; }
|
||
|
};
|
||
|
|
||
|
#endif // INCLUDE_RADIOASTRONOMYSINK_H
|