Method feedHistogramMode

Panoramic/Scanner.cpp:224–270

}

void
SpectrumView::feedHistogramMode(
    const SUFLOAT *psdData,
    SUFREQ freqMin,
    SUFREQ freqMax)
{
  SUFREQ relBw = (freqMax - freqMin) / this->freqRange;
  SUFREQ fStart = (freqMin - this->freqMin) / this->freqRange;
  SUFREQ fEnd   = (freqMax - this->freqMin) / this->freqRange;
  SUFLOAT t;
  SUFLOAT inv = 1. / SIGDIGGER_SCANNER_SPECTRUM_SIZE;
  unsigned int i;
  SUFLOAT accum = 0;

  fStart *= SIGDIGGER_SCANNER_SPECTRUM_SIZE;
  fEnd   *= SIGDIGGER_SCANNER_SPECTRUM_SIZE;
  relBw  *= SIGDIGGER_SCANNER_SPECTRUM_SIZE;

  unsigned int j = static_cast<unsigned int>(fStart);

  // Now, relBw represents the relative size of the range
  // with respecto to the spectrum bin.

  for (i = 0; i < SIGDIGGER_SCANNER_SPECTRUM_SIZE; ++i)
    accum += psdData[i];

  accum *= inv;

  assert(!std::isnan(inv));
  assert(!std::isnan(accum));

  if (floor(fStart) != floor(fEnd)) {
    // Between two bins.
    t = static_cast<SUFLOAT>((fStart - floor(fStart)) / relBw);

    this->psdCount[j] += 1 - t;
    this->psdAccum[j] += (1 - t) * accum;

    if (j + 1 < SIGDIGGER_SCANNER_SPECTRUM_SIZE) {
      this->psdCount[j + 1] += t;
      this->psdAccum[j + 1] += t * accum;
    }
  } else {
    this->psdCount[j] += 1;
    this->psdAccum[j] += accum;
  }
}

void
SpectrumView::feed(