Estimate a (single) ridge from a time-frequency representation, using the ICM minimization method.
icm(modulus, guess, tfspec=numeric(dim(modulus)[2]), subrate=1,
mu=1, lambda=2 * mu, iteration=100)
Returns the estimated ridge and the cost function.
1D array (of same length as the signal) containing the ridge.
1D array containing the cost function.
Time-Frequency representation (real valued).
Initial guess for the algorithm.
Estimate for the contribution of the noise to modulus.
Subsampling rate for ridge estimation.
Coefficient of the ridge's second derivative in cost function.
Coefficient of the ridge's derivative in cost function.
Maximal number of moves.
To accelerate convergence, it is useful to preprocess modulus before
running annealing method. Such a preprocessing (smoothing and
subsampling of modulus) is implemented in icm
. The
parameter subrate specifies the subsampling rate.
See discussions in the text of ``Practical Time-Frequency Analysis''.
corona
, coronoid
, and snake
,
snakoid
.