sdr-visual-suite/internal/cfar/caso.go

package cfar

import "math"

type caso struct {
	guard      int
	train      int
	scaleDb    float64
	wrapAround bool
}

func newCASO(cfg Config) CFAR {
	return &caso{guard: cfg.GuardCells, train: cfg.TrainCells, scaleDb: cfg.ScaleDb, wrapAround: cfg.WrapAround}
}

func (c *caso) Thresholds(spectrum []float64) []float64 {
	n := len(spectrum)
	if n == 0 {
		return nil
	}
	out := make([]float64, n)
	train := c.train
	guard := c.guard
	if train == 0 {
		return out
	}
	inv := 1.0 / float64(train)

	at := func(i int) float64 {
		if c.wrapAround {
			return spectrum[((i%n)+n)%n]
		}
		if i < 0 || i >= n {
			return spectrum[clampInt(i, 0, n-1)]
		}
		return spectrum[i]
	}

	toLinear := func(db float64) float64 {
		return math.Pow(10, db/10.0)
	}

	var leftSum, rightSum float64
	for k := 1; k <= train; k++ {
		leftSum += toLinear(at(0 - guard - k))
		rightSum += toLinear(at(0 + guard + k))
	}
	lm := leftSum * inv
	rm := rightSum * inv
	noise := lm
	if rm < noise {
		noise = rm
	}
	out[0] = 10*math.Log10(noise) + c.scaleDb

	for i := 1; i < n; i++ {
		leftSum -= toLinear(at(i - 1 - guard - train))
		leftSum += toLinear(at(i - guard - 1))
		rightSum -= toLinear(at(i - 1 + guard + 1))
		rightSum += toLinear(at(i + guard + train))
		lm = leftSum * inv
		rm = rightSum * inv
		noise = lm
		if rm < noise {
			noise = rm
		}
		out[i] = 10*math.Log10(noise) + c.scaleDb
	}
	return out
}