Estimate a ridge from a time-frequency representation, using the modified snake method (modified cost function).
snakoid(modulus, guessA, guessB, snakesize=length(guessB),
tfspec=numeric(dim(modulus)[2]), subrate=1, temprate=3, muA=1,
muB=muA, lambdaB=2 * muB, lambdaA=2 * muA, iteration=1000000,
seed=-7, costsub=1, stagnant=20000, plot=TRUE)
Returns a structure containing:
1D array (of same length as the signal) containing the ridge.
1D array containing the cost function.
when set(default), some results will be displayed.
Time-Frequency representation (real valued).
Initial guess for the algorithm (frequency variable).
Initial guess for the algorithm (time variable).
The length of the first guess of time variable.
Estimate for the contribution of srthe noise to modulus.
Subsampling rate for ridge estimation.
Initial value of temperature parameter.
Coefficient of the ridge's derivative in cost function (frequency component).
Coefficient of the ridge's derivative in cost function (time component).
Coefficient of the ridge's second derivative in cost function (time component).
Coefficient of the ridge's second derivative in cost function (frequency component).
Maximal number of moves.
Initialization of random number generator.
Subsampling of cost function in output.
Maximum number of stationary iterations before stopping.
when set(default), some results will be displayed
See discussions in the text of “Practical Time-Frequency Analysis”.
corona
, coronoid
, icm
,
snake
.