#include #include #include #include #include #include #include #include #include #include #include #include #include #include #include std::thread worker; std::mutex fft_mtx; ImGui::WaterFall wtf; FrequencySelect fSel; fftwf_complex *fft_in, *fft_out; fftwf_plan p; float* tempData; int fftSize = 8192 * 8; bool dcbias = true; io::SoapyWrapper soapy; SignalPath sigPath; std::vector _data; std::vector fftTaps; void fftHandler(dsp::complex_t* samples) { fftwf_execute(p); int half = fftSize / 2; for (int i = 0; i < half; i++) { _data.push_back(log10(std::abs(std::complex(fft_out[half + i][0], fft_out[half + i][1])) / (float)fftSize) * 10.0f); } for (int i = 0; i < half; i++) { _data.push_back(log10(std::abs(std::complex(fft_out[i][0], fft_out[i][1])) / (float)fftSize) * 10.0f); } for (int i = 5; i < fftSize; i++) { _data[i] = (_data[i - 4] + _data[i - 3] + _data[i - 2] + _data[i - 1] + _data[i]) / 5.0f; } wtf.pushFFT(_data, fftSize); _data.clear(); } void windowInit() { int sampleRate = 8000000; wtf.setBandwidth(sampleRate); wtf.setCenterFrequency(90500000); wtf.setVFOBandwidth(200000); wtf.setVFOOffset(0); fSel.init(); fSel.setFrequency(90500000); fft_in = (fftwf_complex*) fftw_malloc(sizeof(fftwf_complex) * fftSize); fft_out = (fftwf_complex*) fftw_malloc(sizeof(fftwf_complex) * fftSize); p = fftwf_plan_dft_1d(fftSize, fft_in, fft_out, FFTW_FORWARD, FFTW_ESTIMATE); printf("Starting DSP Thread!\n"); sigPath.init(sampleRate, 20, fftSize, &soapy.output, (dsp::complex_t*)fft_in, fftHandler); sigPath.start(); } int devId = 0; int _devId = -1; int srId = 0; int _srId = -1; bool showExample = false; long freq = 90500000; long _freq = 90500000; int demod = 0; bool state = false; bool mulstate = true; float vfoFreq = 92000000.0f; float lastVfoFreq = 92000000.0f; float volume = 1.0f; float lastVolume = 1.0f; float fftMin = -70.0f; float fftMax = 0.0f; float offset = 0.0f; float lastOffset = -1.0f; float bw = 8000000.0f; float lastBW = -1.0f; int sampleRate = 1000000; bool playing = false; void setVFO(float freq) { float currentOff = wtf.getVFOOfset(); float currentTune = wtf.getCenterFrequency() + currentOff; float delta = freq - currentTune; float newVFO = currentOff + delta; float vfoBW = wtf.getVFOBandwidth(); float vfoBottom = newVFO - (vfoBW / 2.0f); float vfoTop = newVFO + (vfoBW / 2.0f); float view = wtf.getViewOffset(); float viewBW = wtf.getViewBandwidth(); float viewBottom = view - (viewBW / 2.0f); float viewTop = view + (viewBW / 2.0f); float wholeFreq = wtf.getCenterFrequency(); float BW = wtf.getBandwidth(); float bottom = -(BW / 2.0f); float top = (BW / 2.0f); // VFO still fints in the view if (vfoBottom > viewBottom && vfoTop < viewTop) { sigPath.setVFOFrequency(newVFO); wtf.setVFOOffset(newVFO); return; } // VFO too low for current SDR tuning if (vfoBottom < bottom) { wtf.setViewOffset((BW / 2.0f) - (viewBW / 2.0f)); float newVFOOffset = (BW / 2.0f) - (vfoBW / 2.0f) - (viewBW / 10.0f); sigPath.setVFOFrequency(newVFOOffset); wtf.setVFOOffset(newVFOOffset); wtf.setCenterFrequency(freq - newVFOOffset); soapy.setFrequency(freq - newVFOOffset); return; } // VFO too high for current SDR tuning if (vfoTop > top) { wtf.setViewOffset((viewBW / 2.0f) - (BW / 2.0f)); float newVFOOffset = (vfoBW / 2.0f) - (BW / 2.0f) + (viewBW / 10.0f); sigPath.setVFOFrequency(newVFOOffset); wtf.setVFOOffset(newVFOOffset); wtf.setCenterFrequency(freq - newVFOOffset); soapy.setFrequency(freq - newVFOOffset); return; } // VFO is still without the SDR's bandwidth if (delta < 0) { float newViewOff = vfoTop - (viewBW / 2.0f) + (viewBW / 10.0f); float newViewBottom = newViewOff - (viewBW / 2.0f); float newViewTop = newViewOff + (viewBW / 2.0f); if (newViewBottom > bottom) { wtf.setVFOOffset(newVFO); wtf.setViewOffset(newViewOff); sigPath.setVFOFrequency(newVFO); return; } wtf.setViewOffset((BW / 2.0f) - (viewBW / 2.0f)); float newVFOOffset = (BW / 2.0f) - (vfoBW / 2.0f) - (viewBW / 10.0f); sigPath.setVFOFrequency(newVFOOffset); wtf.setVFOOffset(newVFOOffset); wtf.setCenterFrequency(freq - newVFOOffset); soapy.setFrequency(freq - newVFOOffset); } else { float newViewOff = vfoBottom + (viewBW / 2.0f) - (viewBW / 10.0f); float newViewBottom = newViewOff - (viewBW / 2.0f); float newViewTop = newViewOff + (viewBW / 2.0f); if (newViewTop < top) { wtf.setVFOOffset(newVFO); wtf.setViewOffset(newViewOff); sigPath.setVFOFrequency(newVFO); return; } wtf.setViewOffset((viewBW / 2.0f) - (BW / 2.0f)); float newVFOOffset = (vfoBW / 2.0f) - (BW / 2.0f) + (viewBW / 10.0f); sigPath.setVFOFrequency(newVFOOffset); wtf.setVFOOffset(newVFOOffset); wtf.setCenterFrequency(freq - newVFOOffset); soapy.setFrequency(freq - newVFOOffset); } } void drawWindow() { if (fSel.frequencyChanged) { fSel.frequencyChanged = false; setVFO(fSel.frequency); } if (wtf.centerFreqMoved) { wtf.centerFreqMoved = false; soapy.setFrequency(wtf.getCenterFrequency()); fSel.setFrequency(wtf.getCenterFrequency() + wtf.getVFOOfset()); } if (wtf.vfoFreqChanged) { wtf.vfoFreqChanged = false; sigPath.setVFOFrequency(wtf.getVFOOfset()); fSel.setFrequency(wtf.getCenterFrequency() + wtf.getVFOOfset()); } if (volume != lastVolume) { lastVolume = volume; sigPath.setVolume(volume); } if (devId != _devId) { _devId = devId; soapy.setDevice(soapy.devList[devId]); } if (srId != _srId) { _srId = srId; sampleRate = soapy.sampleRates[srId]; printf("Setting sample rate to %f\n", (float)soapy.sampleRates[srId]); soapy.setSampleRate(sampleRate); wtf.setBandwidth(sampleRate); wtf.setViewBandwidth(sampleRate); sigPath.setSampleRate(sampleRate); bw = sampleRate; } ImVec2 vMin = ImGui::GetWindowContentRegionMin(); ImVec2 vMax = ImGui::GetWindowContentRegionMax(); int width = vMax.x - vMin.x; int height = vMax.y - vMin.y; // To Bar if (playing) { if (ImGui::ImageButton(icons::STOP_RAW, ImVec2(30, 30))) { soapy.stop(); playing = false; } } else { if (ImGui::ImageButton(icons::PLAY_RAW, ImVec2(30, 30))) { soapy.start(); soapy.setFrequency(wtf.getCenterFrequency()); playing = true; } } ImGui::SameLine(); ImGui::SetCursorPosY(ImGui::GetCursorPosY() + 8); ImGui::SetNextItemWidth(200); ImGui::SliderFloat("##_2_", &volume, 0.0f, 1.0f, ""); ImGui::SameLine(); fSel.draw(); ImGui::SetCursorPosY(ImGui::GetCursorPosY() - 8); ImGui::SameLine(ImGui::GetWindowWidth() - 40); ImGui::Image(icons::LOGO, ImVec2(30, 30)); ImGui::Columns(3, "WindowColumns", false); ImVec2 winSize = ImGui::GetWindowSize(); ImGui::SetColumnWidth(0, 300); ImGui::SetColumnWidth(1, winSize.x - 300 - 60); ImGui::SetColumnWidth(2, 60); // Left Column ImGui::BeginChild("Left Column"); if (ImGui::CollapsingHeader("Source")) { ImGui::PushItemWidth(ImGui::GetWindowSize().x); ImGui::Combo("##_0_", &devId, soapy.txtDevList.c_str()); if (!playing) { ImGui::Combo("##_1_", &srId, soapy.txtSampleRateList.c_str()); } else { ImGui::Text("%s Samples/s", soapy.txtSampleRateList.c_str()); } if (ImGui::Button("Refresh")) { soapy.refresh(); } } if (ImGui::CollapsingHeader("Radio")) { ImGui::BeginGroup(); ImGui::Columns(4, "RadioModeColumns", false); if (ImGui::RadioButton("NFM", demod == 0) && demod != 0) { sigPath.setDemodulator(SignalPath::DEMOD_NFM); demod = 0; wtf.setVFOBandwidth(12500); wtf.setVFOReference(ImGui::WaterFall::REF_CENTER); } if (ImGui::RadioButton("WFM", demod == 1) && demod != 1) { sigPath.setDemodulator(SignalPath::DEMOD_FM); demod = 1; wtf.setVFOBandwidth(200000); wtf.setVFOReference(ImGui::WaterFall::REF_CENTER); } ImGui::NextColumn(); if (ImGui::RadioButton("AM", demod == 2) && demod != 2) { sigPath.setDemodulator(SignalPath::DEMOD_AM); demod = 2; wtf.setVFOBandwidth(12500); wtf.setVFOReference(ImGui::WaterFall::REF_CENTER); } if (ImGui::RadioButton("DSB", demod == 3) && demod != 3) { demod = 3; }; ImGui::NextColumn(); if (ImGui::RadioButton("USB", demod == 4) && demod != 4) { sigPath.setDemodulator(SignalPath::DEMOD_USB); demod = 4; wtf.setVFOBandwidth(3000); wtf.setVFOReference(ImGui::WaterFall::REF_LOWER); } if (ImGui::RadioButton("CW", demod == 5) && demod != 5) { demod = 5; }; ImGui::NextColumn(); if (ImGui::RadioButton("LSB", demod == 6) && demod != 6) { sigPath.setDemodulator(SignalPath::DEMOD_LSB); demod = 6; wtf.setVFOBandwidth(3000); wtf.setVFOReference(ImGui::WaterFall::REF_UPPER); } if (ImGui::RadioButton("RAW", demod == 7) && demod != 7) { demod = 7; }; ImGui::Columns(1, "EndRadioModeColumns", false); //ImGui::InputInt("Frequency (kHz)", &freq); ImGui::Checkbox("DC Bias Removal", &dcbias); ImGui::EndGroup(); } ImGui::CollapsingHeader("Audio"); ImGui::CollapsingHeader("Display"); ImGui::CollapsingHeader("Recording"); if(ImGui::CollapsingHeader("Debug")) { ImGui::Text("Frame time: %.3f ms/frame", 1000.0f / ImGui::GetIO().Framerate); ImGui::Text("Framerate: %.1f FPS", ImGui::GetIO().Framerate); } ImGui::EndChild(); // Right Column ImGui::NextColumn(); ImGui::BeginChild("Waterfall"); wtf.draw(); ImGui::EndChild(); ImGui::NextColumn(); ImGui::Text("Zoom"); ImGui::VSliderFloat("##_7_", ImVec2(20.0f, 150.0f), &bw, 1000.0f, sampleRate, ""); ImGui::NewLine(); ImGui::Text("Max"); ImGui::VSliderFloat("##_8_", ImVec2(20.0f, 150.0f), &fftMax, -100.0f, 0.0f, ""); ImGui::NewLine(); ImGui::Text("Min"); ImGui::VSliderFloat("##_9_", ImVec2(20.0f, 150.0f), &fftMin, -100.0f, 0.0f, ""); if (bw != lastBW) { lastBW = bw; wtf.setViewBandwidth(bw); wtf.setViewOffset(wtf.getVFOOfset()); } wtf.setFFTMin(fftMin); wtf.setFFTMax(fftMax); wtf.setWaterfallMin(fftMin); wtf.setWaterfallMax(fftMax); }