double Ltas_getLocalPeakHeight (Ltas me, double environmentMin, double environmentMax, double peakMin, double peakMax, int averagingUnits) {
double environmentLow = Sampled_getMean (me, environmentMin, peakMin, 0, averagingUnits, FALSE);
double environmentHigh = Sampled_getMean (me, peakMax, environmentMax, 0, averagingUnits, FALSE);
double peak = Sampled_getMean (me, peakMin, peakMax, 0, averagingUnits, FALSE);
if (environmentLow == NUMundefined || environmentHigh == NUMundefined || peak == NUMundefined) return NUMundefined;
return averagingUnits == 3 ? peak - 0.5 * (environmentLow + environmentHigh) :
Function_convertSpecialToStandardUnit (me, peak / (0.5 * (environmentLow + environmentHigh)), 0, averagingUnits);
}
autoPitch Pitch_subtractLinearFit (Pitch me, int unit) {
try {
autoPitch thee = Pitch_interpolate (me);
/*
* Find the first and last voiced frame.
*/
long imin = thy nx + 1, imax = 0;
for (long i = 1; i <= my nx; i ++)
if (Pitch_isVoiced_i (thee.peek(), i)) { imin = i; break; }
for (long i = imin + 1; i <= my nx; i ++)
if (! Pitch_isVoiced_i (thee.peek(), i)) { imax = i - 1; break; }
long n = imax - imin + 1;
if (n < 3) return thee;
/*
* Compute average pitch and time.
*/
double sum = 0.0;
for (long i = imin; i <= imax; i ++) {
sum += Sampled_getValueAtSample (thee.peek(), i, Pitch_LEVEL_FREQUENCY, unit);
}
double fmean = sum / n;
double tmean = thy x1 + (0.5 * (imin + imax) - 1) * thy dx;
/*
* Compute slope.
*/
double numerator = 0.0, denominator = 0.0;
for (long i = imin; i <= imax; i ++) {
double t = thy x1 + (i - 1) * thy dx - tmean;
double f = Sampled_getValueAtSample (thee.peek(), i, Pitch_LEVEL_FREQUENCY, unit) - fmean;
numerator += f * t;
denominator += t * t;
}
double slope = numerator / denominator;
/*
* Modify frequencies.
*/
for (long i = imin; i <= imax; i ++) {
Pitch_Frame myFrame = & my frame [i], thyFrame = & thy frame [i];
double t = thy x1 + (i - 1) * thy dx - tmean, myFreq = FREQUENCY (myFrame);
double f = Sampled_getValueAtSample (thee.peek(), i, Pitch_LEVEL_FREQUENCY, unit);
f -= slope * t;
if (NOT_VOICED (myFreq))
FREQUENCY (thyFrame) = 0.0;
else
FREQUENCY (thyFrame) = Function_convertSpecialToStandardUnit (me, f, Pitch_LEVEL_FREQUENCY, unit);
}
return thee;
} catch (MelderError) {
Melder_throw (me, U": linear fit not subtracted.");
}
}
double Sampled_getStandardDeviation_standardUnit (Sampled me, double xmin, double xmax, long ilevel, int averagingUnit, int interpolate) {
return Function_convertSpecialToStandardUnit (me, Sampled_getStandardDeviation (me, xmin, xmax, ilevel, averagingUnit, interpolate), ilevel, averagingUnit);
}
double Ltas_getSlope (Ltas me, double f1min, double f1max, double f2min, double f2max, int averagingUnits) {
double low = Sampled_getMean (me, f1min, f1max, 0, averagingUnits, FALSE);
double high = Sampled_getMean (me, f2min, f2max, 0, averagingUnits, FALSE);
if (low == NUMundefined || high == NUMundefined) return NUMundefined;
return averagingUnits == 3 ? high - low : Function_convertSpecialToStandardUnit (me, high / low, 0, averagingUnits);
}