Wideband autonomous SDR analysis engine forked from sdr-visual-suite
You can not select more than 25 topics Topics must start with a letter or number, can include dashes ('-') and can be up to 35 characters long.

219 line
8.0KB

  1. package main
  2. import (
  3. "context"
  4. "encoding/json"
  5. "log"
  6. "os"
  7. "runtime/debug"
  8. "strings"
  9. "sync"
  10. "time"
  11. "sdr-wideband-suite/internal/config"
  12. "sdr-wideband-suite/internal/detector"
  13. "sdr-wideband-suite/internal/dsp"
  14. "sdr-wideband-suite/internal/pipeline"
  15. "sdr-wideband-suite/internal/recorder"
  16. )
  17. func runDSP(ctx context.Context, srcMgr *sourceManager, cfg config.Config, det *detector.Detector, window []float64, h *hub, eventFile *os.File, eventMu *sync.RWMutex, updates <-chan dspUpdate, gpuState *gpuStatus, rec *recorder.Manager, sigSnap *signalSnapshot, extractMgr *extractionManager, phaseSnap *phaseSnapshot) {
  18. defer func() {
  19. if r := recover(); r != nil {
  20. log.Printf("FATAL: runDSP goroutine panic: %v\n%s", r, debug.Stack())
  21. }
  22. }()
  23. rt := newDSPRuntime(cfg, det, window, gpuState)
  24. ticker := time.NewTicker(cfg.FrameInterval())
  25. defer ticker.Stop()
  26. logTicker := time.NewTicker(5 * time.Second)
  27. defer logTicker.Stop()
  28. enc := json.NewEncoder(eventFile)
  29. dcBlocker := dsp.NewDCBlocker(0.995)
  30. state := &phaseState{}
  31. for {
  32. select {
  33. case <-ctx.Done():
  34. return
  35. case <-logTicker.C:
  36. st := srcMgr.Stats()
  37. log.Printf("stats: buf=%d drop=%d reset=%d last=%dms", st.BufferSamples, st.Dropped, st.Resets, st.LastSampleAgoMs)
  38. case upd := <-updates:
  39. rt.applyUpdate(upd, srcMgr, rec, gpuState)
  40. dcBlocker.Reset()
  41. ticker.Reset(rt.cfg.FrameInterval())
  42. case <-ticker.C:
  43. art, err := rt.captureSpectrum(srcMgr, rec, dcBlocker, gpuState)
  44. if err != nil {
  45. log.Printf("read IQ: %v", err)
  46. if strings.Contains(err.Error(), "timeout") {
  47. if err := srcMgr.Restart(rt.cfg); err != nil {
  48. log.Printf("restart failed: %v", err)
  49. }
  50. }
  51. continue
  52. }
  53. if !rt.gotSamples {
  54. log.Printf("received IQ samples")
  55. rt.gotSamples = true
  56. }
  57. state.surveillance = rt.buildSurveillanceResult(art)
  58. state.refinement = rt.runRefinement(art, state.surveillance, extractMgr, rec)
  59. finished := state.surveillance.Finished
  60. thresholds := state.surveillance.Thresholds
  61. noiseFloor := state.surveillance.NoiseFloor
  62. var displaySignals []detector.Signal
  63. if len(art.detailIQ) > 0 {
  64. displaySignals = state.refinement.Result.Signals
  65. stableSignals := rt.det.StableSignals()
  66. streamSignals := displaySignals
  67. if len(stableSignals) > 0 {
  68. streamSignals = stableSignals
  69. }
  70. if rec != nil && len(art.allIQ) > 0 {
  71. log.Printf("LIVEAUDIO DSP: detailIQ=%d displaySignals=%d streamSignals=%d stableSignals=%d allIQ=%d", len(art.detailIQ), len(displaySignals), len(streamSignals), len(stableSignals), len(art.allIQ))
  72. aqCfg := extractionConfig{firTaps: rt.cfg.Recorder.ExtractionTaps, bwMult: rt.cfg.Recorder.ExtractionBwMult}
  73. streamSnips, streamRates := extractForStreaming(extractMgr, art.allIQ, rt.cfg.SampleRate, rt.cfg.CenterHz, streamSignals, rt.streamPhaseState, rt.streamOverlap, aqCfg)
  74. items := make([]recorder.StreamFeedItem, 0, len(streamSignals))
  75. for j, ds := range streamSignals {
  76. className := "<nil>"
  77. if ds.Class != nil {
  78. className = string(ds.Class.ModType)
  79. }
  80. snipLen := 0
  81. if j < len(streamSnips) {
  82. snipLen = len(streamSnips[j])
  83. }
  84. log.Printf("LIVEAUDIO DSP: streamSignal idx=%d id=%d center=%.3fMHz bw=%.0f class=%s snip=%d", j, ds.ID, ds.CenterHz/1e6, ds.BWHz, className, snipLen)
  85. if ds.ID == 0 || ds.Class == nil {
  86. continue
  87. }
  88. if j >= len(streamSnips) || len(streamSnips[j]) == 0 {
  89. continue
  90. }
  91. snipRate := rt.cfg.SampleRate
  92. if j < len(streamRates) && streamRates[j] > 0 {
  93. snipRate = streamRates[j]
  94. }
  95. items = append(items, recorder.StreamFeedItem{Signal: ds, Snippet: streamSnips[j], SnipRate: snipRate})
  96. }
  97. log.Printf("LIVEAUDIO DSP: feedItems=%d", len(items))
  98. if len(items) > 0 {
  99. rec.FeedSnippets(items)
  100. }
  101. }
  102. rt.maintenance(displaySignals, rec)
  103. } else {
  104. displaySignals = rt.det.StableSignals()
  105. }
  106. if rec != nil && len(displaySignals) > 0 {
  107. runtimeInfo := rec.RuntimeInfoBySignalID()
  108. for i := range displaySignals {
  109. if info, ok := runtimeInfo[displaySignals[i].ID]; ok {
  110. displaySignals[i].DemodName = info.DemodName
  111. displaySignals[i].PlaybackMode = info.PlaybackMode
  112. displaySignals[i].StereoState = info.StereoState
  113. }
  114. }
  115. }
  116. state.arbitration = rt.arbitration
  117. state.presentation = state.surveillance.DisplayLevel
  118. if phaseSnap != nil {
  119. phaseSnap.Set(*state)
  120. }
  121. if sigSnap != nil {
  122. sigSnap.set(displaySignals)
  123. }
  124. eventMu.Lock()
  125. for _, ev := range finished {
  126. _ = enc.Encode(ev)
  127. }
  128. eventMu.Unlock()
  129. if rec != nil && len(finished) > 0 {
  130. evCopy := make([]detector.Event, len(finished))
  131. copy(evCopy, finished)
  132. rec.OnEvents(evCopy)
  133. }
  134. var debugInfo *SpectrumDebug
  135. plan := state.refinement.Input.Plan
  136. windowSummary := buildWindowSummary(plan, state.refinement.Input.Windows, state.surveillance.Candidates, state.refinement.Input.WorkItems, state.refinement.Result.Decisions)
  137. var windowStats *RefinementWindowStats
  138. var monitorSummary []pipeline.MonitorWindowStats
  139. if windowSummary != nil {
  140. windowStats = windowSummary.Refinement
  141. monitorSummary = windowSummary.MonitorWindows
  142. }
  143. hasPlan := plan.TotalCandidates > 0 || plan.Budget > 0 || plan.DroppedBySNR > 0 || plan.DroppedByBudget > 0
  144. hasWindows := windowStats != nil && windowStats.Count > 0
  145. if len(thresholds) > 0 || len(displaySignals) > 0 || noiseFloor != 0 || hasPlan || hasWindows {
  146. scoreDebug := make([]map[string]any, 0, len(displaySignals))
  147. for _, s := range displaySignals {
  148. if s.Class == nil || len(s.Class.Scores) == 0 {
  149. scoreDebug = append(scoreDebug, map[string]any{"center_hz": s.CenterHz, "class": nil})
  150. continue
  151. }
  152. scores := make(map[string]float64, len(s.Class.Scores))
  153. for k, v := range s.Class.Scores {
  154. scores[string(k)] = v
  155. }
  156. scoreDebug = append(scoreDebug, map[string]any{
  157. "center_hz": s.CenterHz,
  158. "mod_type": s.Class.ModType,
  159. "confidence": s.Class.Confidence,
  160. "second_best": s.Class.SecondBest,
  161. "scores": scores,
  162. })
  163. }
  164. debugInfo = &SpectrumDebug{Thresholds: thresholds, NoiseFloor: noiseFloor, Scores: scoreDebug}
  165. candidateSources := buildCandidateSourceSummary(state.surveillance.Candidates)
  166. candidateEvidence := buildCandidateEvidenceSummary(state.surveillance.Candidates)
  167. candidateEvidenceStates := buildCandidateEvidenceStateSummary(state.surveillance.Candidates)
  168. candidateWindows := buildCandidateWindowSummary(state.surveillance.Candidates, plan.MonitorWindows)
  169. if len(candidateSources) > 0 {
  170. debugInfo.CandidateSources = candidateSources
  171. }
  172. if len(candidateEvidence) > 0 {
  173. debugInfo.CandidateEvidence = candidateEvidence
  174. }
  175. if candidateEvidenceStates != nil {
  176. debugInfo.CandidateEvidenceStates = candidateEvidenceStates
  177. }
  178. if len(candidateWindows) > 0 {
  179. debugInfo.CandidateWindows = candidateWindows
  180. }
  181. if len(monitorSummary) > 0 {
  182. debugInfo.MonitorWindowStats = monitorSummary
  183. }
  184. if windowSummary != nil {
  185. debugInfo.WindowSummary = windowSummary
  186. }
  187. if hasPlan {
  188. debugInfo.RefinementPlan = &plan
  189. }
  190. if hasWindows {
  191. debugInfo.Windows = windowStats
  192. }
  193. refinementDebug := &RefinementDebug{}
  194. if hasPlan {
  195. refinementDebug.Plan = &plan
  196. refinementDebug.Request = &state.refinement.Input.Request
  197. refinementDebug.WorkItems = state.refinement.Input.WorkItems
  198. }
  199. if hasWindows {
  200. refinementDebug.Windows = windowStats
  201. }
  202. if len(monitorSummary) > 0 {
  203. refinementDebug.MonitorWindowStats = monitorSummary
  204. }
  205. if windowSummary != nil {
  206. refinementDebug.WindowSummary = windowSummary
  207. }
  208. refinementDebug.Arbitration = buildArbitrationSnapshot(state.refinement, state.arbitration)
  209. debugInfo.Refinement = refinementDebug
  210. }
  211. h.broadcast(SpectrumFrame{Timestamp: art.now.UnixMilli(), CenterHz: rt.cfg.CenterHz, SampleHz: rt.cfg.SampleRate, FFTSize: rt.cfg.FFTSize, Spectrum: art.surveillanceSpectrum, Signals: displaySignals, Debug: debugInfo})
  212. }
  213. }
  214. }