Wideband autonomous SDR analysis engine forked from sdr-visual-suite
您最多选择25个主题 主题必须以字母或数字开头,可以包含连字符 (-),并且长度不得超过35个字符

935 行
25KB

  1. package main
  2. import (
  3. "context"
  4. "encoding/json"
  5. "flag"
  6. "log"
  7. "math"
  8. "net/http"
  9. "os"
  10. "os/signal"
  11. "path/filepath"
  12. "runtime/debug"
  13. "sort"
  14. "sdr-visual-suite/internal/config"
  15. "strconv"
  16. "strings"
  17. "sync"
  18. "syscall"
  19. "time"
  20. "github.com/gorilla/websocket"
  21. "sdr-visual-suite/internal/classifier"
  22. "sdr-visual-suite/internal/detector"
  23. "sdr-visual-suite/internal/dsp"
  24. "sdr-visual-suite/internal/events"
  25. fftutil "sdr-visual-suite/internal/fft"
  26. "sdr-visual-suite/internal/fft/gpufft"
  27. "sdr-visual-suite/internal/mock"
  28. "sdr-visual-suite/internal/recorder"
  29. "sdr-visual-suite/internal/runtime"
  30. "sdr-visual-suite/internal/sdr"
  31. "sdr-visual-suite/internal/sdrplay"
  32. )
  33. type SpectrumFrame struct {
  34. Timestamp int64 `json:"ts"`
  35. CenterHz float64 `json:"center_hz"`
  36. SampleHz int `json:"sample_rate"`
  37. FFTSize int `json:"fft_size"`
  38. Spectrum []float64 `json:"spectrum_db"`
  39. Signals []detector.Signal `json:"signals"`
  40. }
  41. type client struct {
  42. conn *websocket.Conn
  43. send chan []byte
  44. done chan struct{}
  45. closeOnce sync.Once
  46. }
  47. type hub struct {
  48. mu sync.Mutex
  49. clients map[*client]struct{}
  50. frameCnt int64
  51. lastLogTs time.Time
  52. }
  53. type gpuStatus struct {
  54. mu sync.RWMutex
  55. Available bool `json:"available"`
  56. Active bool `json:"active"`
  57. Error string `json:"error"`
  58. }
  59. type signalSnapshot struct {
  60. mu sync.RWMutex
  61. signals []detector.Signal
  62. }
  63. func (s *signalSnapshot) set(sig []detector.Signal) {
  64. s.mu.Lock()
  65. defer s.mu.Unlock()
  66. s.signals = append([]detector.Signal(nil), sig...)
  67. }
  68. func (s *signalSnapshot) get() []detector.Signal {
  69. s.mu.RLock()
  70. defer s.mu.RUnlock()
  71. return append([]detector.Signal(nil), s.signals...)
  72. }
  73. func (g *gpuStatus) set(active bool, err error) {
  74. g.mu.Lock()
  75. defer g.mu.Unlock()
  76. g.Active = active
  77. if err != nil {
  78. g.Error = err.Error()
  79. } else {
  80. g.Error = ""
  81. }
  82. }
  83. func (g *gpuStatus) snapshot() gpuStatus {
  84. g.mu.RLock()
  85. defer g.mu.RUnlock()
  86. return gpuStatus{Available: g.Available, Active: g.Active, Error: g.Error}
  87. }
  88. func newHub() *hub {
  89. return &hub{clients: map[*client]struct{}{}, lastLogTs: time.Now()}
  90. }
  91. func (h *hub) add(c *client) {
  92. h.mu.Lock()
  93. defer h.mu.Unlock()
  94. h.clients[c] = struct{}{}
  95. log.Printf("ws connected (%d clients)", len(h.clients))
  96. }
  97. func (h *hub) remove(c *client) {
  98. c.closeOnce.Do(func() { close(c.done) })
  99. h.mu.Lock()
  100. defer h.mu.Unlock()
  101. delete(h.clients, c)
  102. log.Printf("ws disconnected (%d clients)", len(h.clients))
  103. }
  104. func (h *hub) broadcast(frame SpectrumFrame) {
  105. b, err := json.Marshal(frame)
  106. if err != nil {
  107. log.Printf("marshal frame: %v", err)
  108. return
  109. }
  110. h.mu.Lock()
  111. clients := make([]*client, 0, len(h.clients))
  112. for c := range h.clients {
  113. clients = append(clients, c)
  114. }
  115. h.mu.Unlock()
  116. for _, c := range clients {
  117. select {
  118. case c.send <- b:
  119. default:
  120. h.remove(c)
  121. }
  122. }
  123. h.frameCnt++
  124. if time.Since(h.lastLogTs) > 2*time.Second {
  125. h.lastLogTs = time.Now()
  126. log.Printf("broadcast frames=%d clients=%d", h.frameCnt, len(clients))
  127. }
  128. }
  129. type sourceManager struct {
  130. mu sync.RWMutex
  131. src sdr.Source
  132. newSource func(cfg config.Config) (sdr.Source, error)
  133. }
  134. func (m *sourceManager) Restart(cfg config.Config) error {
  135. m.mu.Lock()
  136. defer m.mu.Unlock()
  137. old := m.src
  138. _ = old.Stop()
  139. next, err := m.newSource(cfg)
  140. if err != nil {
  141. _ = old.Start()
  142. m.src = old
  143. return err
  144. }
  145. if err := next.Start(); err != nil {
  146. _ = next.Stop()
  147. _ = old.Start()
  148. m.src = old
  149. return err
  150. }
  151. m.src = next
  152. return nil
  153. }
  154. func (m *sourceManager) Stats() sdr.SourceStats {
  155. m.mu.RLock()
  156. defer m.mu.RUnlock()
  157. if sp, ok := m.src.(sdr.StatsProvider); ok {
  158. return sp.Stats()
  159. }
  160. return sdr.SourceStats{}
  161. }
  162. func (m *sourceManager) Flush() {
  163. m.mu.RLock()
  164. defer m.mu.RUnlock()
  165. if fl, ok := m.src.(sdr.Flushable); ok {
  166. fl.Flush()
  167. }
  168. }
  169. func newSourceManager(src sdr.Source, newSource func(cfg config.Config) (sdr.Source, error)) *sourceManager {
  170. return &sourceManager{src: src, newSource: newSource}
  171. }
  172. func (m *sourceManager) Start() error {
  173. m.mu.RLock()
  174. defer m.mu.RUnlock()
  175. return m.src.Start()
  176. }
  177. func (m *sourceManager) Stop() error {
  178. m.mu.RLock()
  179. defer m.mu.RUnlock()
  180. return m.src.Stop()
  181. }
  182. func (m *sourceManager) ReadIQ(n int) ([]complex64, error) {
  183. m.mu.RLock()
  184. defer m.mu.RUnlock()
  185. return m.src.ReadIQ(n)
  186. }
  187. func (m *sourceManager) ApplyConfig(cfg config.Config) error {
  188. m.mu.Lock()
  189. defer m.mu.Unlock()
  190. if updatable, ok := m.src.(sdr.ConfigurableSource); ok {
  191. if err := updatable.UpdateConfig(cfg.SampleRate, cfg.CenterHz, cfg.GainDb, cfg.AGC, cfg.TunerBwKHz); err == nil {
  192. return nil
  193. }
  194. }
  195. old := m.src
  196. _ = old.Stop()
  197. next, err := m.newSource(cfg)
  198. if err != nil {
  199. _ = old.Start()
  200. return err
  201. }
  202. if err := next.Start(); err != nil {
  203. _ = next.Stop()
  204. _ = old.Start()
  205. return err
  206. }
  207. m.src = next
  208. return nil
  209. }
  210. type dspUpdate struct {
  211. cfg config.Config
  212. det *detector.Detector
  213. window []float64
  214. dcBlock bool
  215. iqBalance bool
  216. useGPUFFT bool
  217. }
  218. func pushDSPUpdate(ch chan dspUpdate, update dspUpdate) {
  219. select {
  220. case ch <- update:
  221. default:
  222. select {
  223. case <-ch:
  224. default:
  225. }
  226. ch <- update
  227. }
  228. }
  229. func main() {
  230. var cfgPath string
  231. var mockFlag bool
  232. flag.StringVar(&cfgPath, "config", "config.yaml", "path to config YAML")
  233. flag.BoolVar(&mockFlag, "mock", false, "use synthetic IQ source")
  234. flag.Parse()
  235. cfg, err := config.Load(cfgPath)
  236. if err != nil {
  237. log.Fatalf("load config: %v", err)
  238. }
  239. cfgManager := runtime.New(cfg)
  240. gpuState := &gpuStatus{Available: gpufft.Available()}
  241. newSource := func(cfg config.Config) (sdr.Source, error) {
  242. if mockFlag {
  243. src := mock.New(cfg.SampleRate)
  244. if updatable, ok := interface{}(src).(sdr.ConfigurableSource); ok {
  245. _ = updatable.UpdateConfig(cfg.SampleRate, cfg.CenterHz, cfg.GainDb, cfg.AGC, cfg.TunerBwKHz)
  246. }
  247. return src, nil
  248. }
  249. src, err := sdrplay.New(cfg.SampleRate, cfg.CenterHz, cfg.GainDb, cfg.TunerBwKHz)
  250. if err != nil {
  251. return nil, err
  252. }
  253. if updatable, ok := src.(sdr.ConfigurableSource); ok {
  254. _ = updatable.UpdateConfig(cfg.SampleRate, cfg.CenterHz, cfg.GainDb, cfg.AGC, cfg.TunerBwKHz)
  255. }
  256. return src, nil
  257. }
  258. src, err := newSource(cfg)
  259. if err != nil {
  260. log.Fatalf("sdrplay init failed: %v (try --mock or build with -tags sdrplay)", err)
  261. }
  262. srcMgr := newSourceManager(src, newSource)
  263. if err := srcMgr.Start(); err != nil {
  264. log.Fatalf("source start: %v", err)
  265. }
  266. defer srcMgr.Stop()
  267. if err := os.MkdirAll(filepath.Dir(cfg.EventPath), 0o755); err != nil {
  268. log.Fatalf("event path: %v", err)
  269. }
  270. eventFile, err := os.OpenFile(cfg.EventPath, os.O_CREATE|os.O_APPEND|os.O_WRONLY, 0o644)
  271. if err != nil {
  272. log.Fatalf("open events: %v", err)
  273. }
  274. defer eventFile.Close()
  275. eventMu := &sync.RWMutex{}
  276. det := detector.New(cfg.Detector.ThresholdDb, cfg.SampleRate, cfg.FFTSize,
  277. time.Duration(cfg.Detector.MinDurationMs)*time.Millisecond,
  278. time.Duration(cfg.Detector.HoldMs)*time.Millisecond,
  279. cfg.Detector.EmaAlpha,
  280. cfg.Detector.HysteresisDb,
  281. cfg.Detector.MinStableFrames,
  282. time.Duration(cfg.Detector.GapToleranceMs)*time.Millisecond)
  283. window := fftutil.Hann(cfg.FFTSize)
  284. h := newHub()
  285. dspUpdates := make(chan dspUpdate, 1)
  286. ctx, cancel := context.WithCancel(context.Background())
  287. defer cancel()
  288. decodeMap := buildDecoderMap(cfg)
  289. recMgr := recorder.New(cfg.SampleRate, cfg.FFTSize, recorder.Policy{
  290. Enabled: cfg.Recorder.Enabled,
  291. MinSNRDb: cfg.Recorder.MinSNRDb,
  292. MinDuration: mustParseDuration(cfg.Recorder.MinDuration, 1*time.Second),
  293. MaxDuration: mustParseDuration(cfg.Recorder.MaxDuration, 300*time.Second),
  294. PrerollMs: cfg.Recorder.PrerollMs,
  295. RecordIQ: cfg.Recorder.RecordIQ,
  296. RecordAudio: cfg.Recorder.RecordAudio,
  297. AutoDemod: cfg.Recorder.AutoDemod,
  298. AutoDecode: cfg.Recorder.AutoDecode,
  299. MaxDiskMB: cfg.Recorder.MaxDiskMB,
  300. OutputDir: cfg.Recorder.OutputDir,
  301. ClassFilter: cfg.Recorder.ClassFilter,
  302. RingSeconds: cfg.Recorder.RingSeconds,
  303. }, cfg.CenterHz, decodeMap)
  304. sigSnap := &signalSnapshot{}
  305. go runDSP(ctx, srcMgr, cfg, det, window, h, eventFile, eventMu, dspUpdates, gpuState, recMgr, sigSnap)
  306. upgrader := websocket.Upgrader{CheckOrigin: func(r *http.Request) bool {
  307. origin := r.Header.Get("Origin")
  308. if origin == "" || origin == "null" {
  309. return true
  310. }
  311. // allow same-host or any local IP
  312. return true
  313. }}
  314. http.HandleFunc("/ws", func(w http.ResponseWriter, r *http.Request) {
  315. conn, err := upgrader.Upgrade(w, r, nil)
  316. if err != nil {
  317. log.Printf("ws upgrade failed: %v (origin: %s)", err, r.Header.Get("Origin"))
  318. return
  319. }
  320. c := &client{conn: conn, send: make(chan []byte, 32), done: make(chan struct{})}
  321. h.add(c)
  322. defer func() {
  323. h.remove(c)
  324. _ = conn.Close()
  325. }()
  326. conn.SetReadDeadline(time.Now().Add(60 * time.Second))
  327. conn.SetPongHandler(func(string) error {
  328. conn.SetReadDeadline(time.Now().Add(60 * time.Second))
  329. return nil
  330. })
  331. go func() {
  332. ping := time.NewTicker(30 * time.Second)
  333. defer ping.Stop()
  334. for {
  335. select {
  336. case msg, ok := <-c.send:
  337. if !ok {
  338. return
  339. }
  340. _ = conn.SetWriteDeadline(time.Now().Add(200 * time.Millisecond))
  341. if err := conn.WriteMessage(websocket.TextMessage, msg); err != nil {
  342. return
  343. }
  344. case <-ping.C:
  345. _ = conn.SetWriteDeadline(time.Now().Add(5 * time.Second))
  346. if err := conn.WriteMessage(websocket.PingMessage, nil); err != nil {
  347. log.Printf("ws ping error: %v", err)
  348. return
  349. }
  350. }
  351. }
  352. }()
  353. for {
  354. _, _, err := conn.ReadMessage()
  355. if err != nil {
  356. return
  357. }
  358. }
  359. })
  360. http.HandleFunc("/api/config", func(w http.ResponseWriter, r *http.Request) {
  361. w.Header().Set("Content-Type", "application/json")
  362. switch r.Method {
  363. case http.MethodGet:
  364. _ = json.NewEncoder(w).Encode(cfgManager.Snapshot())
  365. case http.MethodPost:
  366. var update runtime.ConfigUpdate
  367. if err := json.NewDecoder(r.Body).Decode(&update); err != nil {
  368. http.Error(w, "invalid json", http.StatusBadRequest)
  369. return
  370. }
  371. prev := cfgManager.Snapshot()
  372. next, err := cfgManager.ApplyConfig(update)
  373. if err != nil {
  374. http.Error(w, err.Error(), http.StatusBadRequest)
  375. return
  376. }
  377. sourceChanged := prev.CenterHz != next.CenterHz || prev.SampleRate != next.SampleRate || prev.GainDb != next.GainDb || prev.AGC != next.AGC || prev.TunerBwKHz != next.TunerBwKHz
  378. if sourceChanged {
  379. if err := srcMgr.ApplyConfig(next); err != nil {
  380. cfgManager.Replace(prev)
  381. http.Error(w, "failed to apply source config", http.StatusInternalServerError)
  382. return
  383. }
  384. }
  385. if err := config.Save(cfgPath, next); err != nil {
  386. log.Printf("config save failed: %v", err)
  387. }
  388. detChanged := prev.Detector.ThresholdDb != next.Detector.ThresholdDb ||
  389. prev.Detector.MinDurationMs != next.Detector.MinDurationMs ||
  390. prev.Detector.HoldMs != next.Detector.HoldMs ||
  391. prev.Detector.EmaAlpha != next.Detector.EmaAlpha ||
  392. prev.Detector.HysteresisDb != next.Detector.HysteresisDb ||
  393. prev.Detector.MinStableFrames != next.Detector.MinStableFrames ||
  394. prev.Detector.GapToleranceMs != next.Detector.GapToleranceMs ||
  395. prev.SampleRate != next.SampleRate ||
  396. prev.FFTSize != next.FFTSize
  397. windowChanged := prev.FFTSize != next.FFTSize
  398. var newDet *detector.Detector
  399. var newWindow []float64
  400. if detChanged {
  401. newDet = detector.New(next.Detector.ThresholdDb, next.SampleRate, next.FFTSize,
  402. time.Duration(next.Detector.MinDurationMs)*time.Millisecond,
  403. time.Duration(next.Detector.HoldMs)*time.Millisecond,
  404. next.Detector.EmaAlpha,
  405. next.Detector.HysteresisDb,
  406. next.Detector.MinStableFrames,
  407. time.Duration(next.Detector.GapToleranceMs)*time.Millisecond)
  408. }
  409. if windowChanged {
  410. newWindow = fftutil.Hann(next.FFTSize)
  411. }
  412. pushDSPUpdate(dspUpdates, dspUpdate{
  413. cfg: next,
  414. det: newDet,
  415. window: newWindow,
  416. dcBlock: next.DCBlock,
  417. iqBalance: next.IQBalance,
  418. useGPUFFT: next.UseGPUFFT,
  419. })
  420. _ = json.NewEncoder(w).Encode(next)
  421. default:
  422. http.Error(w, "method not allowed", http.StatusMethodNotAllowed)
  423. }
  424. })
  425. http.HandleFunc("/api/sdr/settings", func(w http.ResponseWriter, r *http.Request) {
  426. w.Header().Set("Content-Type", "application/json")
  427. if r.Method != http.MethodPost {
  428. http.Error(w, "method not allowed", http.StatusMethodNotAllowed)
  429. return
  430. }
  431. var update runtime.SettingsUpdate
  432. if err := json.NewDecoder(r.Body).Decode(&update); err != nil {
  433. http.Error(w, "invalid json", http.StatusBadRequest)
  434. return
  435. }
  436. prev := cfgManager.Snapshot()
  437. next, err := cfgManager.ApplySettings(update)
  438. if err != nil {
  439. http.Error(w, err.Error(), http.StatusBadRequest)
  440. return
  441. }
  442. if prev.AGC != next.AGC || prev.TunerBwKHz != next.TunerBwKHz {
  443. if err := srcMgr.ApplyConfig(next); err != nil {
  444. cfgManager.Replace(prev)
  445. http.Error(w, "failed to apply sdr settings", http.StatusInternalServerError)
  446. return
  447. }
  448. }
  449. if prev.DCBlock != next.DCBlock || prev.IQBalance != next.IQBalance {
  450. pushDSPUpdate(dspUpdates, dspUpdate{
  451. cfg: next,
  452. dcBlock: next.DCBlock,
  453. iqBalance: next.IQBalance,
  454. })
  455. }
  456. if err := config.Save(cfgPath, next); err != nil {
  457. log.Printf("config save failed: %v", err)
  458. }
  459. _ = json.NewEncoder(w).Encode(next)
  460. })
  461. http.HandleFunc("/api/stats", func(w http.ResponseWriter, r *http.Request) {
  462. w.Header().Set("Content-Type", "application/json")
  463. _ = json.NewEncoder(w).Encode(srcMgr.Stats())
  464. })
  465. http.HandleFunc("/api/gpu", func(w http.ResponseWriter, r *http.Request) {
  466. w.Header().Set("Content-Type", "application/json")
  467. _ = json.NewEncoder(w).Encode(gpuState.snapshot())
  468. })
  469. http.HandleFunc("/api/events", func(w http.ResponseWriter, r *http.Request) {
  470. w.Header().Set("Content-Type", "application/json")
  471. limit := 200
  472. if v := r.URL.Query().Get("limit"); v != "" {
  473. if parsed, err := strconv.Atoi(v); err == nil {
  474. limit = parsed
  475. }
  476. }
  477. var since time.Time
  478. if v := r.URL.Query().Get("since"); v != "" {
  479. if parsed, err := parseSince(v); err == nil {
  480. since = parsed
  481. } else {
  482. http.Error(w, "invalid since", http.StatusBadRequest)
  483. return
  484. }
  485. }
  486. snap := cfgManager.Snapshot()
  487. eventMu.RLock()
  488. evs, err := events.ReadRecent(snap.EventPath, limit, since)
  489. eventMu.RUnlock()
  490. if err != nil {
  491. http.Error(w, "failed to read events", http.StatusInternalServerError)
  492. return
  493. }
  494. _ = json.NewEncoder(w).Encode(evs)
  495. })
  496. http.HandleFunc("/api/signals", func(w http.ResponseWriter, r *http.Request) {
  497. w.Header().Set("Content-Type", "application/json")
  498. if sigSnap == nil {
  499. _ = json.NewEncoder(w).Encode([]detector.Signal{})
  500. return
  501. }
  502. _ = json.NewEncoder(w).Encode(sigSnap.get())
  503. })
  504. http.HandleFunc("/api/decoders", func(w http.ResponseWriter, r *http.Request) {
  505. w.Header().Set("Content-Type", "application/json")
  506. _ = json.NewEncoder(w).Encode(decoderKeys(cfgManager.Snapshot()))
  507. })
  508. http.HandleFunc("/api/recordings", func(w http.ResponseWriter, r *http.Request) {
  509. if r.Method != http.MethodGet {
  510. http.Error(w, "method not allowed", http.StatusMethodNotAllowed)
  511. return
  512. }
  513. w.Header().Set("Content-Type", "application/json")
  514. snap := cfgManager.Snapshot()
  515. list, err := recorder.ListRecordings(snap.Recorder.OutputDir)
  516. if err != nil {
  517. http.Error(w, "failed to list recordings", http.StatusInternalServerError)
  518. return
  519. }
  520. _ = json.NewEncoder(w).Encode(list)
  521. })
  522. http.HandleFunc("/api/recordings/", func(w http.ResponseWriter, r *http.Request) {
  523. w.Header().Set("Content-Type", "application/json")
  524. id := strings.TrimPrefix(r.URL.Path, "/api/recordings/")
  525. if id == "" {
  526. http.Error(w, "missing id", http.StatusBadRequest)
  527. return
  528. }
  529. snap := cfgManager.Snapshot()
  530. base := filepath.Clean(filepath.Join(snap.Recorder.OutputDir, id))
  531. if !strings.HasPrefix(base, filepath.Clean(snap.Recorder.OutputDir)) {
  532. http.Error(w, "invalid path", http.StatusBadRequest)
  533. return
  534. }
  535. if r.URL.Path == "/api/recordings/"+id+"/audio" {
  536. http.ServeFile(w, r, filepath.Join(base, "audio.wav"))
  537. return
  538. }
  539. if r.URL.Path == "/api/recordings/"+id+"/iq" {
  540. http.ServeFile(w, r, filepath.Join(base, "signal.cf32"))
  541. return
  542. }
  543. if r.URL.Path == "/api/recordings/"+id+"/decode" {
  544. mode := r.URL.Query().Get("mode")
  545. cmd := buildDecoderMap(cfgManager.Snapshot())[mode]
  546. if cmd == "" {
  547. http.Error(w, "decoder not configured", http.StatusBadRequest)
  548. return
  549. }
  550. meta, err := recorder.ReadMeta(filepath.Join(base, "meta.json"))
  551. if err != nil {
  552. http.Error(w, "meta read failed", http.StatusInternalServerError)
  553. return
  554. }
  555. audioPath := filepath.Join(base, "audio.wav")
  556. if _, errStat := os.Stat(audioPath); errStat != nil {
  557. audioPath = ""
  558. }
  559. res, err := recorder.DecodeOnDemand(cmd, filepath.Join(base, "signal.cf32"), meta.SampleRate, audioPath)
  560. if err != nil {
  561. http.Error(w, res.Stderr, http.StatusInternalServerError)
  562. return
  563. }
  564. _ = json.NewEncoder(w).Encode(res)
  565. return
  566. }
  567. // default: meta.json
  568. http.ServeFile(w, r, filepath.Join(base, "meta.json"))
  569. })
  570. http.HandleFunc("/api/demod", func(w http.ResponseWriter, r *http.Request) {
  571. if r.Method != http.MethodGet {
  572. http.Error(w, "method not allowed", http.StatusMethodNotAllowed)
  573. return
  574. }
  575. q := r.URL.Query()
  576. freq, _ := strconv.ParseFloat(q.Get("freq"), 64)
  577. bw, _ := strconv.ParseFloat(q.Get("bw"), 64)
  578. sec, _ := strconv.Atoi(q.Get("sec"))
  579. if sec < 1 {
  580. sec = 1
  581. }
  582. if sec > 10 {
  583. sec = 10
  584. }
  585. mode := q.Get("mode")
  586. data, _, err := recMgr.DemodLive(freq, bw, mode, sec)
  587. if err != nil {
  588. http.Error(w, err.Error(), http.StatusBadRequest)
  589. return
  590. }
  591. w.Header().Set("Content-Type", "audio/wav")
  592. _, _ = w.Write(data)
  593. })
  594. http.Handle("/", http.FileServer(http.Dir(cfg.WebRoot)))
  595. server := &http.Server{Addr: cfg.WebAddr}
  596. go func() {
  597. log.Printf("web listening on %s", cfg.WebAddr)
  598. if err := server.ListenAndServe(); err != nil && err != http.ErrServerClosed {
  599. log.Fatalf("server: %v", err)
  600. }
  601. }()
  602. stop := make(chan os.Signal, 1)
  603. signal.Notify(stop, os.Interrupt, syscall.SIGTERM)
  604. <-stop
  605. ctxTimeout, cancelTimeout := context.WithTimeout(context.Background(), 5*time.Second)
  606. defer cancelTimeout()
  607. _ = server.Shutdown(ctxTimeout)
  608. }
  609. 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) {
  610. defer func() {
  611. if r := recover(); r != nil {
  612. log.Printf("FATAL: runDSP goroutine panic: %v\n%s", r, debug.Stack())
  613. }
  614. }()
  615. ticker := time.NewTicker(cfg.FrameInterval())
  616. defer ticker.Stop()
  617. logTicker := time.NewTicker(5 * time.Second)
  618. defer logTicker.Stop()
  619. enc := json.NewEncoder(eventFile)
  620. dcBlocker := dsp.NewDCBlocker(0.995)
  621. dcEnabled := cfg.DCBlock
  622. iqEnabled := cfg.IQBalance
  623. plan := fftutil.NewCmplxPlan(cfg.FFTSize)
  624. useGPU := cfg.UseGPUFFT
  625. var gpuEngine *gpufft.Engine
  626. if useGPU && gpuState != nil {
  627. snap := gpuState.snapshot()
  628. if snap.Available {
  629. if eng, err := gpufft.New(cfg.FFTSize); err == nil {
  630. gpuEngine = eng
  631. gpuState.set(true, nil)
  632. } else {
  633. gpuState.set(false, err)
  634. useGPU = false
  635. }
  636. } else {
  637. gpuState.set(false, nil)
  638. useGPU = false
  639. }
  640. } else if gpuState != nil {
  641. gpuState.set(false, nil)
  642. }
  643. gotSamples := false
  644. for {
  645. select {
  646. case <-ctx.Done():
  647. return
  648. case <-logTicker.C:
  649. st := srcMgr.Stats()
  650. log.Printf("stats: buf=%d drop=%d reset=%d last=%dms", st.BufferSamples, st.Dropped, st.Resets, st.LastSampleAgoMs)
  651. case upd := <-updates:
  652. prevFFT := cfg.FFTSize
  653. prevUseGPU := useGPU
  654. cfg = upd.cfg
  655. if rec != nil {
  656. rec.Update(cfg.SampleRate, cfg.FFTSize, recorder.Policy{
  657. Enabled: cfg.Recorder.Enabled,
  658. MinSNRDb: cfg.Recorder.MinSNRDb,
  659. MinDuration: mustParseDuration(cfg.Recorder.MinDuration, 1*time.Second),
  660. MaxDuration: mustParseDuration(cfg.Recorder.MaxDuration, 300*time.Second),
  661. PrerollMs: cfg.Recorder.PrerollMs,
  662. RecordIQ: cfg.Recorder.RecordIQ,
  663. RecordAudio: cfg.Recorder.RecordAudio,
  664. AutoDemod: cfg.Recorder.AutoDemod,
  665. AutoDecode: cfg.Recorder.AutoDecode,
  666. MaxDiskMB: cfg.Recorder.MaxDiskMB,
  667. OutputDir: cfg.Recorder.OutputDir,
  668. ClassFilter: cfg.Recorder.ClassFilter,
  669. RingSeconds: cfg.Recorder.RingSeconds,
  670. }, cfg.CenterHz, buildDecoderMap(cfg))
  671. }
  672. if upd.det != nil {
  673. det = upd.det
  674. }
  675. if upd.window != nil {
  676. window = upd.window
  677. plan = fftutil.NewCmplxPlan(cfg.FFTSize)
  678. }
  679. dcEnabled = upd.dcBlock
  680. iqEnabled = upd.iqBalance
  681. if cfg.FFTSize != prevFFT || cfg.UseGPUFFT != prevUseGPU {
  682. srcMgr.Flush()
  683. gotSamples = false
  684. if gpuEngine != nil {
  685. gpuEngine.Close()
  686. gpuEngine = nil
  687. }
  688. useGPU = cfg.UseGPUFFT
  689. if useGPU && gpuState != nil {
  690. snap := gpuState.snapshot()
  691. if snap.Available {
  692. if eng, err := gpufft.New(cfg.FFTSize); err == nil {
  693. gpuEngine = eng
  694. gpuState.set(true, nil)
  695. } else {
  696. gpuState.set(false, err)
  697. useGPU = false
  698. }
  699. } else {
  700. gpuState.set(false, nil)
  701. useGPU = false
  702. }
  703. } else if gpuState != nil {
  704. gpuState.set(false, nil)
  705. }
  706. }
  707. dcBlocker.Reset()
  708. ticker.Reset(cfg.FrameInterval())
  709. case <-ticker.C:
  710. iq, err := srcMgr.ReadIQ(cfg.FFTSize)
  711. if err != nil {
  712. log.Printf("read IQ: %v", err)
  713. if strings.Contains(err.Error(), "timeout") {
  714. if err := srcMgr.Restart(cfg); err != nil {
  715. log.Printf("restart failed: %v", err)
  716. }
  717. }
  718. continue
  719. }
  720. if rec != nil {
  721. rec.Ingest(time.Now(), iq)
  722. }
  723. if !gotSamples {
  724. log.Printf("received IQ samples")
  725. gotSamples = true
  726. }
  727. if dcEnabled {
  728. dcBlocker.Apply(iq)
  729. }
  730. if iqEnabled {
  731. dsp.IQBalance(iq)
  732. }
  733. var spectrum []float64
  734. if useGPU && gpuEngine != nil {
  735. if len(window) == len(iq) {
  736. for i := 0; i < len(iq); i++ {
  737. v := iq[i]
  738. w := float32(window[i])
  739. iq[i] = complex(real(v)*w, imag(v)*w)
  740. }
  741. }
  742. out, err := gpuEngine.Exec(iq)
  743. if err != nil {
  744. if gpuState != nil {
  745. gpuState.set(false, err)
  746. }
  747. useGPU = false
  748. spectrum = fftutil.SpectrumWithPlan(iq, nil, plan)
  749. } else {
  750. spectrum = fftutil.SpectrumFromFFT(out)
  751. }
  752. } else {
  753. spectrum = fftutil.SpectrumWithPlan(iq, window, plan)
  754. }
  755. for i := range spectrum {
  756. if math.IsNaN(spectrum[i]) || math.IsInf(spectrum[i], 0) {
  757. spectrum[i] = -200
  758. }
  759. }
  760. now := time.Now()
  761. finished, signals := det.Process(now, spectrum, cfg.CenterHz)
  762. // enrich classification with temporal IQ features on per-signal snippet
  763. if len(iq) > 0 {
  764. for i := range signals {
  765. snip := extractSignalIQ(iq, cfg.SampleRate, cfg.CenterHz, signals[i].CenterHz, signals[i].BWHz)
  766. cls := classifier.Classify(classifier.SignalInput{FirstBin: signals[i].FirstBin, LastBin: signals[i].LastBin, SNRDb: signals[i].SNRDb}, spectrum, cfg.SampleRate, cfg.FFTSize, snip)
  767. signals[i].Class = cls
  768. }
  769. det.UpdateClasses(signals)
  770. }
  771. if sigSnap != nil {
  772. sigSnap.set(signals)
  773. }
  774. eventMu.Lock()
  775. for _, ev := range finished {
  776. _ = enc.Encode(ev)
  777. }
  778. eventMu.Unlock()
  779. if rec != nil {
  780. rec.OnEvents(finished)
  781. }
  782. h.broadcast(SpectrumFrame{
  783. Timestamp: now.UnixMilli(),
  784. CenterHz: cfg.CenterHz,
  785. SampleHz: cfg.SampleRate,
  786. FFTSize: cfg.FFTSize,
  787. Spectrum: spectrum,
  788. Signals: signals,
  789. })
  790. }
  791. }
  792. }
  793. func mustParseDuration(raw string, fallback time.Duration) time.Duration {
  794. if raw == "" {
  795. return fallback
  796. }
  797. if d, err := time.ParseDuration(raw); err == nil {
  798. return d
  799. }
  800. return fallback
  801. }
  802. func buildDecoderMap(cfg config.Config) map[string]string {
  803. out := map[string]string{}
  804. if cfg.Decoder.FT8Cmd != "" {
  805. out["FT8"] = cfg.Decoder.FT8Cmd
  806. }
  807. if cfg.Decoder.WSPRCmd != "" {
  808. out["WSPR"] = cfg.Decoder.WSPRCmd
  809. }
  810. if cfg.Decoder.DMRCmd != "" {
  811. out["DMR"] = cfg.Decoder.DMRCmd
  812. }
  813. if cfg.Decoder.DStarCmd != "" {
  814. out["D-STAR"] = cfg.Decoder.DStarCmd
  815. }
  816. if cfg.Decoder.FSKCmd != "" {
  817. out["FSK"] = cfg.Decoder.FSKCmd
  818. }
  819. if cfg.Decoder.PSKCmd != "" {
  820. out["PSK"] = cfg.Decoder.PSKCmd
  821. }
  822. return out
  823. }
  824. func decoderKeys(cfg config.Config) []string {
  825. m := buildDecoderMap(cfg)
  826. keys := make([]string, 0, len(m))
  827. for k := range m {
  828. keys = append(keys, k)
  829. }
  830. sort.Strings(keys)
  831. return keys
  832. }
  833. func extractSignalIQ(iq []complex64, sampleRate int, centerHz float64, sigHz float64, bwHz float64) []complex64 {
  834. if len(iq) == 0 || sampleRate <= 0 {
  835. return nil
  836. }
  837. offset := sigHz - centerHz
  838. shifted := dsp.FreqShift(iq, sampleRate, offset)
  839. cutoff := bwHz / 2
  840. if cutoff < 200 {
  841. cutoff = 200
  842. }
  843. if cutoff > float64(sampleRate)/2-1 {
  844. cutoff = float64(sampleRate)/2 - 1
  845. }
  846. taps := dsp.LowpassFIR(cutoff, sampleRate, 101)
  847. filtered := dsp.ApplyFIR(shifted, taps)
  848. decim := sampleRate / 200000
  849. if decim < 1 {
  850. decim = 1
  851. }
  852. return dsp.Decimate(filtered, decim)
  853. }
  854. func parseSince(raw string) (time.Time, error) {
  855. if raw == "" {
  856. return time.Time{}, nil
  857. }
  858. if ms, err := strconv.ParseInt(raw, 10, 64); err == nil {
  859. if ms > 1e12 {
  860. return time.UnixMilli(ms), nil
  861. }
  862. return time.Unix(ms, 0), nil
  863. }
  864. if t, err := time.Parse(time.RFC3339Nano, raw); err == nil {
  865. return t, nil
  866. }
  867. return time.Parse(time.RFC3339, raw)
  868. }