I wonder if just searching for local maxima first would help? Quick
search yielded something similar with a threshold added:

https://www.dsprelated.com/showcode/217.php

Not my area of expertise at all, so take with a grain of salt. ;)

On Mon, Oct 19, 2020 at 9:27 AM Victor Lazzarini <Victor.Lazzarini@mu.ie> wrote:
>
> Not really a solution, more of a diagnostic. The lumping is due to the fact that windowing makes a sine component appear
> at three or four bins (in the case of a Hanning window) around the peak. So I suppose, depending on the intensity of the peak
> relative to everything else, this would explain it.
> ========================
> Prof. Victor Lazzarini
> Maynooth University
> Ireland
>
> > On 19 Oct 2020, at 14:00, Oeyvind Brandtsegg <obrandts@gmail.com> wrote:
> >
> > WARNINGThis email originated from outside of Maynooth University's Mail System. Do not reply, click links or open attachments unless you recognise the sender and know the content is safe.
> > Hi
> > I seem to recall that someone (Joachim?) did some work with pvstrace, extracting the strongest unique frequencies from a spectrum. With a signal that contains only a few pure sines, pvsanal seems to "lump" the bins around those frequencies. I would like to set a minimum threshold of separation (in Hz), so that I could read the strongest partials of a spectrum with pvstrace. Before delving into doing this myself, I wanted to ask if someone had already come up with a solution?
> > all best
> > Oeyvind
> > Csound mailing list Csound@listserv.heanet.ie https://listserv.heanet.ie/cgi-bin/wa?A0=CSOUND Send bugs reports to https://github.com/csound/csound/issues Discussions of bugs and features can be posted here
>
> Csound mailing list
> Csound@listserv.heanet.ie
> https://listserv.heanet.ie/cgi-bin/wa?A0=CSOUND
> Send bugs reports to
>         https://github.com/csound/csound/issues
> Discussions of bugs and features can be posted here

Csound mailing list
Csound@listserv.heanet.ie
https://listserv.heanet.ie/cgi-bin/wa?A0=CSOUND
Send bugs reports to
        https://github.com/csound/csound/issues
Discussions of bugs and features can be posted here

Thanks for the tip. It is much more elegant than what I had in mind.

Here's my implementation (skipped the low thresh in the original example):

<CsoundSynthesizer>
<CsOptions>
-n --displays
</CsOptions>
<CsInstruments>

ksmps = 64
nchnls = 2
0dbfs = 1

giShape ftgen 0, 0, 256, 10, .1, 0.5, 0.2, 0.1,0.7

;***************************************************
; find local maxima
; from matlab Author: sparafucile17 10/02/2003
; https://www.dsprelated.com/showcode/217.php
;***************************************************
opcode LocalMax, k[], i
giData xin
kndx init 1
kpeakindx init 0
kPeaks[] init ftlen(giData)/2
kprev_amp table 0, giData
kprev_slope init 1; allow a maxima point at the second sample
while kndx < ftlen(giData) do
kamp table kndx, giData
kslope = kamp - kprev_amp
if (kamp < kprev_amp) && (kprev_slope > 0) then
kPeaks[kpeakindx] = kndx-1
kpeakindx += 1
endif
kprev_amp = kamp
kprev_slope = kslope
kndx += 1
od
xout kPeaks
endop

instr 1
 kTest[] LocalMax giShape
 kndx init 0
 while kTest[kndx] > 0 do
   printf "amp %f at index %i \n", kndx+1, table(kTest[kndx],giShape), kTest[kndx]
   kndx += 1
 od
 endin

</CsInstruments>
<CsScore>
; start dur
i1 0 1
e
</CsScore>
</CsoundSynthesizer>

man. 19. okt. 2020 kl. 16:40 skrev Steven Yi <stevenyi@gmail.com>:

I wonder if just searching for local maxima first would help? Quick
search yielded something similar with a threshold added:

https://www.dsprelated.com/showcode/217.php

Not my area of expertise at all, so take with a grain of salt. ;)

On Mon, Oct 19, 2020 at 9:27 AM Victor Lazzarini <Victor.Lazzarini@mu.ie> wrote:
>
> Not really a solution, more of a diagnostic. The lumping is due to the fact that windowing makes a sine component appear
> at three or four bins (in the case of a Hanning window) around the peak. So I suppose, depending on the intensity of the peak
> relative to everything else, this would explain it.
> ========================
> Prof. Victor Lazzarini
> Maynooth University
> Ireland
>
> > On 19 Oct 2020, at 14:00, Oeyvind Brandtsegg <obrandts@gmail.com> wrote:
> >
> > WARNINGThis email originated from outside of Maynooth University's Mail System. Do not reply, click links or open attachments unless you recognise the sender and know the content is safe.
> > Hi
> > I seem to recall that someone (Joachim?) did some work with pvstrace, extracting the strongest unique frequencies from a spectrum. With a signal that contains only a few pure sines, pvsanal seems to "lump" the bins around those frequencies. I would like to set a minimum threshold of separation (in Hz), so that I could read the strongest partials of a spectrum with pvstrace. Before delving into doing this myself, I wanted to ask if someone had already come up with a solution?
> > all best
> > Oeyvind
> > Csound mailing list Csound@listserv.heanet.ie https://listserv.heanet.ie/cgi-bin/wa?A0=CSOUND Send bugs reports to https://github.com/csound/csound/issues Discussions of bugs and features can be posted here
>
> Csound mailing list
> Csound@listserv.heanet.ie
> https://listserv.heanet.ie/cgi-bin/wa?A0=CSOUND
> Send bugs reports to
>         https://github.com/csound/csound/issues
> Discussions of bugs and features can be posted here

Csound mailing list
Csound@listserv.heanet.ie
https://listserv.heanet.ie/cgi-bin/wa?A0=CSOUND
Send bugs reports to
        https://github.com/csound/csound/issues
Discussions of bugs and features can be posted here