2015-11-18 21:27:37 -05:00
|
|
|
// Copyright (C) 2012 maintech GmbH, Otto-Hahn-Str. 15, 97204 Hoechberg, Germany //
|
|
|
|
// (C) 2015 John Greb //
|
|
|
|
// //
|
|
|
|
// 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 //
|
|
|
|
// //
|
|
|
|
// 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/>. //
|
|
|
|
///////////////////////////////////////////////////////////////////////////////////
|
|
|
|
|
|
|
|
#include "udpsrc.h"
|
|
|
|
|
2015-12-01 20:28:31 -05:00
|
|
|
#include <QUdpSocket>
|
2015-11-18 21:27:37 -05:00
|
|
|
#include <QThread>
|
|
|
|
#include "dsp/channelizer.h"
|
|
|
|
#include "udpsrcgui.h"
|
|
|
|
|
|
|
|
MESSAGE_CLASS_DEFINITION(UDPSrc::MsgUDPSrcConfigure, Message)
|
|
|
|
MESSAGE_CLASS_DEFINITION(UDPSrc::MsgUDPSrcSpectrum, Message)
|
|
|
|
|
2015-11-19 02:55:58 -05:00
|
|
|
UDPSrc::UDPSrc(MessageQueue* uiMessageQueue, UDPSrcGUI* udpSrcGUI, SampleSink* spectrum) :
|
2015-11-18 21:27:37 -05:00
|
|
|
m_settingsMutex(QMutex::Recursive)
|
|
|
|
{
|
2015-11-19 02:55:58 -05:00
|
|
|
setObjectName("UDPSrc");
|
2015-11-18 21:27:37 -05:00
|
|
|
|
2015-12-01 20:28:31 -05:00
|
|
|
m_socket = new QUdpSocket(this);
|
|
|
|
|
2015-11-18 21:27:37 -05:00
|
|
|
m_inputSampleRate = 96000;
|
|
|
|
m_sampleFormat = FormatSSB;
|
|
|
|
m_outputSampleRate = 48000;
|
|
|
|
m_rfBandwidth = 32000;
|
2015-11-19 02:55:58 -05:00
|
|
|
m_udpPort = 9999;
|
2015-11-18 21:27:37 -05:00
|
|
|
m_nco.setFreq(0, m_inputSampleRate);
|
|
|
|
m_interpolator.create(16, m_inputSampleRate, m_rfBandwidth / 2.0);
|
|
|
|
m_sampleDistanceRemain = m_inputSampleRate / m_outputSampleRate;
|
|
|
|
m_uiMessageQueue = uiMessageQueue;
|
2015-11-19 02:55:58 -05:00
|
|
|
m_udpSrcGUI = udpSrcGUI;
|
2015-11-18 21:27:37 -05:00
|
|
|
m_spectrum = spectrum;
|
|
|
|
m_spectrumEnabled = false;
|
|
|
|
m_nextSSBId = 0;
|
|
|
|
m_nextS16leId = 0;
|
|
|
|
|
|
|
|
m_last = 0;
|
|
|
|
m_this = 0;
|
|
|
|
m_scale = 0;
|
|
|
|
m_boost = 0;
|
|
|
|
m_magsq = 0;
|
2015-11-19 02:55:58 -05:00
|
|
|
m_sampleBufferSSB.resize(udpFftLen);
|
|
|
|
UDPFilter = new fftfilt(0.3 / 48.0, 16.0 / 48.0, udpFftLen);
|
2015-11-18 21:27:37 -05:00
|
|
|
// if (!TCPFilter) segfault;
|
|
|
|
}
|
|
|
|
|
|
|
|
UDPSrc::~UDPSrc()
|
|
|
|
{
|
2015-12-01 20:28:31 -05:00
|
|
|
delete m_socket;
|
2015-11-19 02:55:58 -05:00
|
|
|
if (UDPFilter) delete UDPFilter;
|
2015-11-18 21:27:37 -05:00
|
|
|
}
|
|
|
|
|
2015-12-01 20:28:31 -05:00
|
|
|
void UDPSrc::configure(MessageQueue* messageQueue, SampleFormat sampleFormat, Real outputSampleRate, Real rfBandwidth, QString& udpAddress, int udpPort, int boost)
|
2015-11-18 21:27:37 -05:00
|
|
|
{
|
2015-12-01 20:28:31 -05:00
|
|
|
Message* cmd = MsgUDPSrcConfigure::create(sampleFormat, outputSampleRate, rfBandwidth, udpAddress, udpPort, boost);
|
2015-11-18 21:27:37 -05:00
|
|
|
messageQueue->push(cmd);
|
|
|
|
}
|
|
|
|
|
|
|
|
void UDPSrc::setSpectrum(MessageQueue* messageQueue, bool enabled)
|
|
|
|
{
|
|
|
|
Message* cmd = MsgUDPSrcSpectrum::create(enabled);
|
|
|
|
messageQueue->push(cmd);
|
|
|
|
}
|
|
|
|
|
|
|
|
void UDPSrc::feed(const SampleVector::const_iterator& begin, const SampleVector::const_iterator& end, bool positiveOnly)
|
|
|
|
{
|
|
|
|
Complex ci;
|
|
|
|
fftfilt::cmplx* sideband;
|
|
|
|
Real l, r;
|
|
|
|
|
|
|
|
m_sampleBuffer.clear();
|
|
|
|
|
|
|
|
m_settingsMutex.lock();
|
|
|
|
|
|
|
|
// Rtl-Sdr uses full 16-bit scale; FCDPP does not
|
|
|
|
//int rescale = 32768 * (1 << m_boost);
|
|
|
|
int rescale = (1 << m_boost);
|
|
|
|
|
|
|
|
for(SampleVector::const_iterator it = begin; it < end; ++it) {
|
|
|
|
//Complex c(it->real() / 32768.0f, it->imag() / 32768.0f);
|
|
|
|
Complex c(it->real(), it->imag());
|
|
|
|
c *= m_nco.nextIQ();
|
|
|
|
|
|
|
|
if(m_interpolator.interpolate(&m_sampleDistanceRemain, c, &ci))
|
|
|
|
{
|
|
|
|
m_magsq = ((ci.real()*ci.real() + ci.imag()*ci.imag())*rescale*rescale) / (1<<30);
|
|
|
|
m_sampleBuffer.push_back(Sample(ci.real() * rescale, ci.imag() * rescale));
|
|
|
|
m_sampleDistanceRemain += m_inputSampleRate / m_outputSampleRate;
|
|
|
|
}
|
|
|
|
}
|
|
|
|
|
|
|
|
if((m_spectrum != 0) && (m_spectrumEnabled))
|
|
|
|
{
|
|
|
|
m_spectrum->feed(m_sampleBuffer.begin(), m_sampleBuffer.end(), positiveOnly);
|
|
|
|
}
|
|
|
|
|
2015-12-01 20:28:31 -05:00
|
|
|
if (m_sampleFormat == FormatSSB)
|
2015-11-18 21:27:37 -05:00
|
|
|
{
|
2015-12-01 20:28:31 -05:00
|
|
|
for(SampleVector::const_iterator it = m_sampleBuffer.begin(); it != m_sampleBuffer.end(); ++it)
|
|
|
|
{
|
2015-11-18 21:27:37 -05:00
|
|
|
//Complex cj(it->real() / 30000.0, it->imag() / 30000.0);
|
|
|
|
Complex cj(it->real(), it->imag());
|
2015-11-19 02:55:58 -05:00
|
|
|
int n_out = UDPFilter->runSSB(cj, &sideband, true);
|
2015-12-01 20:28:31 -05:00
|
|
|
|
|
|
|
if (n_out)
|
|
|
|
{
|
|
|
|
for (int i = 0; i < n_out; i+=2)
|
|
|
|
{
|
2015-11-18 21:27:37 -05:00
|
|
|
//l = (sideband[i].real() + sideband[i].imag()) * 0.7 * 32000.0;
|
|
|
|
//r = (sideband[i+1].real() + sideband[i+1].imag()) * 0.7 * 32000.0;
|
|
|
|
l = (sideband[i].real() + sideband[i].imag()) * 0.7;
|
|
|
|
r = (sideband[i+1].real() + sideband[i+1].imag()) * 0.7;
|
|
|
|
m_sampleBufferSSB.push_back(Sample(l, r));
|
|
|
|
}
|
2015-12-01 20:28:31 -05:00
|
|
|
|
|
|
|
m_socket->writeDatagram((const char*)&m_sampleBufferSSB[0], (qint64 ) (n_out * 2), m_udpAddress, m_udpPort);
|
2015-11-18 21:27:37 -05:00
|
|
|
m_sampleBufferSSB.clear();
|
|
|
|
}
|
|
|
|
}
|
|
|
|
}
|
2015-12-01 20:28:31 -05:00
|
|
|
else if (m_sampleFormat == FormatNFM)
|
|
|
|
{
|
|
|
|
for(SampleVector::const_iterator it = m_sampleBuffer.begin(); it != m_sampleBuffer.end(); ++it)
|
|
|
|
{
|
2015-11-18 21:27:37 -05:00
|
|
|
Complex cj(it->real() / 32768.0f, it->imag() / 32768.0f);
|
|
|
|
// An FFT filter here is overkill, but was already set up for SSB
|
2015-11-19 02:55:58 -05:00
|
|
|
int n_out = UDPFilter->runFilt(cj, &sideband);
|
2015-12-01 20:28:31 -05:00
|
|
|
|
|
|
|
if (n_out)
|
|
|
|
{
|
2015-11-18 21:27:37 -05:00
|
|
|
Real sum = 1.0;
|
2015-12-01 20:28:31 -05:00
|
|
|
for (int i = 0; i < n_out; i+=2)
|
|
|
|
{
|
2015-11-18 21:27:37 -05:00
|
|
|
l = m_this.real() * (m_last.imag() - sideband[i].imag())
|
|
|
|
- m_this.imag() * (m_last.real() - sideband[i].real());
|
|
|
|
m_last = sideband[i];
|
|
|
|
r = m_last.real() * (m_this.imag() - sideband[i+1].imag())
|
|
|
|
- m_last.imag() * (m_this.real() - sideband[i+1].real());
|
|
|
|
m_this = sideband[i+1];
|
|
|
|
m_sampleBufferSSB.push_back(Sample(l * m_scale, r * m_scale));
|
|
|
|
sum += m_this.real() * m_this.real() + m_this.imag() * m_this.imag();
|
|
|
|
}
|
|
|
|
// TODO: correct levels
|
2015-11-19 02:55:58 -05:00
|
|
|
m_scale = 24000 * udpFftLen / sum;
|
2015-12-01 20:28:31 -05:00
|
|
|
m_socket->writeDatagram((const char*)&m_sampleBufferSSB[0], (qint64 ) (n_out * 2), m_udpAddress, m_udpPort);
|
2015-11-18 21:27:37 -05:00
|
|
|
m_sampleBufferSSB.clear();
|
|
|
|
}
|
|
|
|
}
|
|
|
|
}
|
2015-12-01 20:28:31 -05:00
|
|
|
else
|
|
|
|
{
|
|
|
|
m_socket->writeDatagram((const char*)&m_sampleBuffer[0], (qint64 ) (m_sampleBuffer.size() * 4), m_udpAddress, m_udpPort);
|
|
|
|
}
|
2015-11-18 21:27:37 -05:00
|
|
|
|
|
|
|
m_settingsMutex.unlock();
|
|
|
|
}
|
|
|
|
|
|
|
|
void UDPSrc::start()
|
|
|
|
{
|
|
|
|
}
|
|
|
|
|
|
|
|
void UDPSrc::stop()
|
|
|
|
{
|
|
|
|
}
|
|
|
|
|
|
|
|
bool UDPSrc::handleMessage(const Message& cmd)
|
|
|
|
{
|
|
|
|
qDebug() << "UDPSrc::handleMessage";
|
|
|
|
|
|
|
|
if (Channelizer::MsgChannelizerNotification::match(cmd))
|
|
|
|
{
|
|
|
|
Channelizer::MsgChannelizerNotification& notif = (Channelizer::MsgChannelizerNotification&) cmd;
|
|
|
|
|
|
|
|
m_settingsMutex.lock();
|
|
|
|
|
|
|
|
m_inputSampleRate = notif.getSampleRate();
|
|
|
|
m_nco.setFreq(-notif.getFrequencyOffset(), m_inputSampleRate);
|
|
|
|
m_interpolator.create(16, m_inputSampleRate, m_rfBandwidth / 2.0);
|
|
|
|
m_sampleDistanceRemain = m_inputSampleRate / m_outputSampleRate;
|
|
|
|
|
|
|
|
m_settingsMutex.unlock();
|
|
|
|
|
2015-11-19 02:55:58 -05:00
|
|
|
qDebug() << "UDPSrc::handleMessage: MsgChannelizerNotification: m_inputSampleRate: " << m_inputSampleRate
|
2015-11-18 21:27:37 -05:00
|
|
|
<< " frequencyOffset: " << notif.getFrequencyOffset();
|
|
|
|
|
|
|
|
return true;
|
|
|
|
}
|
|
|
|
else if (MsgUDPSrcConfigure::match(cmd))
|
|
|
|
{
|
|
|
|
MsgUDPSrcConfigure& cfg = (MsgUDPSrcConfigure&) cmd;
|
|
|
|
|
|
|
|
m_settingsMutex.lock();
|
|
|
|
|
|
|
|
m_sampleFormat = cfg.getSampleFormat();
|
|
|
|
m_outputSampleRate = cfg.getOutputSampleRate();
|
|
|
|
m_rfBandwidth = cfg.getRFBandwidth();
|
|
|
|
|
2015-12-01 20:28:31 -05:00
|
|
|
if (cfg.getUDPAddress() != m_udpAddress.toString())
|
|
|
|
{
|
|
|
|
m_udpAddress.setAddress(cfg.getUDPAddress());
|
|
|
|
}
|
|
|
|
|
2015-11-19 02:55:58 -05:00
|
|
|
if (cfg.getUDPPort() != m_udpPort)
|
2015-11-18 21:27:37 -05:00
|
|
|
{
|
2015-11-19 02:55:58 -05:00
|
|
|
m_udpPort = cfg.getUDPPort();
|
2015-11-18 21:27:37 -05:00
|
|
|
}
|
|
|
|
|
|
|
|
m_boost = cfg.getBoost();
|
|
|
|
m_interpolator.create(16, m_inputSampleRate, m_rfBandwidth / 2.0);
|
|
|
|
m_sampleDistanceRemain = m_inputSampleRate / m_outputSampleRate;
|
|
|
|
|
|
|
|
if (m_sampleFormat == FormatSSB)
|
|
|
|
{
|
2015-11-19 02:55:58 -05:00
|
|
|
UDPFilter->create_filter(0.3 / 48.0, m_rfBandwidth / 2.0 / m_outputSampleRate);
|
2015-11-18 21:27:37 -05:00
|
|
|
}
|
|
|
|
else
|
|
|
|
{
|
2015-11-19 02:55:58 -05:00
|
|
|
UDPFilter->create_filter(0.0, m_rfBandwidth / 2.0 / m_outputSampleRate);
|
2015-11-18 21:27:37 -05:00
|
|
|
}
|
|
|
|
|
|
|
|
m_settingsMutex.unlock();
|
|
|
|
|
2015-11-19 02:55:58 -05:00
|
|
|
qDebug() << " - MsgUDPSrcConfigure: m_sampleFormat: " << m_sampleFormat
|
2015-11-18 21:27:37 -05:00
|
|
|
<< " m_outputSampleRate: " << m_outputSampleRate
|
|
|
|
<< " m_rfBandwidth: " << m_rfBandwidth
|
2015-12-01 20:28:31 -05:00
|
|
|
<< " m_boost: " << m_boost
|
|
|
|
<< " m_udpAddress: " << cfg.getUDPAddress()
|
|
|
|
<< " m_udpPort: " << m_udpPort;
|
2015-11-18 21:27:37 -05:00
|
|
|
|
|
|
|
return true;
|
|
|
|
}
|
|
|
|
else if (MsgUDPSrcSpectrum::match(cmd))
|
|
|
|
{
|
|
|
|
MsgUDPSrcSpectrum& spc = (MsgUDPSrcSpectrum&) cmd;
|
|
|
|
|
|
|
|
m_spectrumEnabled = spc.getEnabled();
|
|
|
|
|
2015-11-19 02:55:58 -05:00
|
|
|
qDebug() << " - MsgUDPSrcSpectrum: m_spectrumEnabled: " << m_spectrumEnabled;
|
2015-11-18 21:27:37 -05:00
|
|
|
|
|
|
|
return true;
|
|
|
|
}
|
|
|
|
else
|
|
|
|
{
|
|
|
|
if(m_spectrum != 0)
|
|
|
|
{
|
|
|
|
return m_spectrum->handleMessage(cmd);
|
|
|
|
}
|
|
|
|
else
|
|
|
|
{
|
|
|
|
return false;
|
|
|
|
}
|
|
|
|
}
|
|
|
|
}
|