1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
|
#ifndef WATERFALL_GLOBAL_DATA_CPP
#define WATERFALL_GLOBAL_DATA_CPP
#include <waterfallGlobalData.h>
#include <cstdio>
WaterfallData::WaterfallData(const double minimumFrequency,
const double maximumFrequency,
const uint64_t fftPoints,
const unsigned int historyExtent)
#if QWT_VERSION < 0x060000
: QwtRasterData(QwtDoubleRect(minimumFrequency /* X START */, 0 /* Y START */,
maximumFrequency - minimumFrequency /* WIDTH */,
static_cast<double>(historyExtent)/* HEIGHT */))
#else
: QwtRasterData()
#endif
{
_intensityRange = QwtDoubleInterval(-200.0, 0.0);
_fftPoints = fftPoints;
_historyLength = historyExtent;
_spectrumData = new double[_fftPoints * _historyLength];
#if QWT_VERSION >= 0x060000
setInterval(Qt::XAxis, QwtInterval(minimumFrequency, maximumFrequency));
setInterval(Qt::YAxis, QwtInterval(0, historyExtent));
setInterval(Qt::ZAxis, QwtInterval(-200, 0.0));
#endif
Reset();
}
WaterfallData::~WaterfallData()
{
delete[] _spectrumData;
}
void WaterfallData::Reset()
{
memset(_spectrumData, 0x0, _fftPoints*_historyLength*sizeof(double));
_numLinesToUpdate = -1;
}
void WaterfallData::Copy(const WaterfallData* rhs)
{
#if QWT_VERSION < 0x060000
if((_fftPoints != rhs->GetNumFFTPoints()) ||
(boundingRect() != rhs->boundingRect()) ){
_fftPoints = rhs->GetNumFFTPoints();
setBoundingRect(rhs->boundingRect());
delete[] _spectrumData;
_spectrumData = new double[_fftPoints * _historyLength];
}
#else
if(_fftPoints != rhs->GetNumFFTPoints()) {
_fftPoints = rhs->GetNumFFTPoints();
delete[] _spectrumData;
_spectrumData = new double[_fftPoints * _historyLength];
}
#endif
Reset();
SetSpectrumDataBuffer(rhs->GetSpectrumDataBuffer());
SetNumLinesToUpdate(rhs->GetNumLinesToUpdate());
#if QWT_VERSION < 0x060000
setRange(rhs->range());
#else
setInterval(Qt::XAxis, rhs->interval(Qt::XAxis));
setInterval(Qt::YAxis, rhs->interval(Qt::YAxis));
setInterval(Qt::ZAxis, rhs->interval(Qt::ZAxis));
#endif
}
void WaterfallData::ResizeData(const double startFreq,
const double stopFreq,
const uint64_t fftPoints)
{
#if QWT_VERSION < 0x060000
if((fftPoints != GetNumFFTPoints()) ||
(boundingRect().width() != (stopFreq - startFreq)) ||
(boundingRect().left() != startFreq)){
setBoundingRect(QwtDoubleRect(startFreq, 0,
stopFreq-startFreq,
boundingRect().height()));
_fftPoints = fftPoints;
delete[] _spectrumData;
_spectrumData = new double[_fftPoints * _historyLength];
}
#else
if((fftPoints != GetNumFFTPoints()) ||
(interval(Qt::XAxis).width() != (stopFreq - startFreq)) ||
(interval(Qt::XAxis).minValue() != startFreq)){
setInterval(Qt::XAxis, QwtInterval(startFreq, stopFreq));
_fftPoints = fftPoints;
delete[] _spectrumData;
_spectrumData = new double[_fftPoints * _historyLength];
}
#endif
Reset();
}
QwtRasterData *WaterfallData::copy() const
{
#if QWT_VERSION < 0x060000
WaterfallData* returnData = new WaterfallData(boundingRect().left(),
boundingRect().right(),
_fftPoints, _historyLength);
#else
WaterfallData* returnData = new WaterfallData(interval(Qt::XAxis).minValue(),
interval(Qt::XAxis).maxValue(),
_fftPoints, _historyLength);
#endif
returnData->Copy(this);
return returnData;
}
#if QWT_VERSION < 0x060000
QwtDoubleInterval WaterfallData::range() const
{
return _intensityRange;
}
void WaterfallData::setRange(const QwtDoubleInterval& newRange)
{
_intensityRange = newRange;
}
#endif
double WaterfallData::value(double x, double y) const
{
double returnValue = 0.0;
#if QWT_VERSION < 0x060000
const unsigned int intY = static_cast<unsigned int>((1.0 - (y/boundingRect().height())) *
static_cast<double>(_historyLength-1));
const unsigned int intX = static_cast<unsigned int>((((x - boundingRect().left()) / boundingRect().width()) *
static_cast<double>(_fftPoints-1)) + 0.5);
#else
double height = interval(Qt::YAxis).maxValue();
double left = interval(Qt::XAxis).minValue();
double right = interval(Qt::XAxis).maxValue();
double ylen = static_cast<double>(_historyLength-1);
double xlen = static_cast<double>(_fftPoints-1);
const unsigned int intY = static_cast<unsigned int>((1.0 - y/height) * ylen);
const unsigned int intX = static_cast<unsigned int>((((x - left) / (right-left)) * xlen) + 0.5);
#endif
const int location = (intY * _fftPoints) + intX;
if((location > -1) && (location < static_cast<int64_t>(_fftPoints * _historyLength))){
returnValue = _spectrumData[location];
}
return returnValue;
}
uint64_t WaterfallData::GetNumFFTPoints() const
{
return _fftPoints;
}
void WaterfallData::addFFTData(const double* fftData,
const uint64_t fftDataSize,
const int droppedFrames){
if(fftDataSize == _fftPoints){
int64_t heightOffset = _historyLength - 1 - droppedFrames;
uint64_t drawingDroppedFrames = droppedFrames;
// Any valid data rolled off the display so just fill in zeros and write new data
if(heightOffset < 0){
heightOffset = 0;
drawingDroppedFrames = static_cast<uint64_t>(_historyLength-1);
}
// Copy the old data over if any available
if(heightOffset > 0){
memmove( _spectrumData, &_spectrumData[(drawingDroppedFrames+1) * _fftPoints],
heightOffset * _fftPoints * sizeof(double)) ;
}
if(drawingDroppedFrames > 0){
// Fill in zeros data for dropped data
memset(&_spectrumData[heightOffset * _fftPoints], 0x00,
static_cast<int64_t>(drawingDroppedFrames) * _fftPoints * sizeof(double));
}
// add the new buffer
memcpy(&_spectrumData[(_historyLength - 1) * _fftPoints], fftData,
_fftPoints*sizeof(double));
}
}
double* WaterfallData::GetSpectrumDataBuffer() const
{
return _spectrumData;
}
void WaterfallData::SetSpectrumDataBuffer(const double* newData)
{
memcpy(_spectrumData, newData, _fftPoints * _historyLength * sizeof(double));
}
int WaterfallData::GetNumLinesToUpdate() const
{
return _numLinesToUpdate;
}
void WaterfallData::SetNumLinesToUpdate(const int newNum)
{
_numLinesToUpdate = newNum;
}
void WaterfallData::IncrementNumLinesToUpdate()
{
_numLinesToUpdate++;
}
#endif /* WATERFALL_GLOBAL_DATA_CPP */
|
