FormantFilter Sound_and_Pitch_to_FormantFilter (Sound me, Pitch thee, double analysisWidth, double dt,
double f1_hz, double fmax_hz, double df_hz, double relative_bw) {
try {
double t1, windowDuration = 2 * analysisWidth; /* gaussian window */
double nyquist = 0.5 / my dx, samplingFrequency = 2 * nyquist, fmin_hz = 0;
long nt, f0_undefined = 0;
if (my xmin > thy xmin || my xmax > thy xmax) Melder_throw
("The domain of the Sound is not included in the domain of the Pitch.");
double f0_median = Pitch_getQuantile (thee, thy xmin, thy xmax, 0.5, kPitch_unit_HERTZ);
if (f0_median == NUMundefined || f0_median == 0) {
f0_median = 100;
Melder_warning (L"Pitch values undefined. Bandwith fixed to 100 Hz. ");
}
if (f1_hz <= 0) {
f1_hz = 100;
}
if (fmax_hz <= 0) {
fmax_hz = nyquist;
}
if (df_hz <= 0) {
df_hz = f0_median / 2;
}
if (relative_bw <= 0) {
relative_bw = 1.1;
}
fmax_hz = MIN (fmax_hz, nyquist);
long nf = floor ( (fmax_hz - f1_hz) / df_hz + 0.5);
Sampled_shortTermAnalysis (me, windowDuration, dt, &nt, &t1);
autoFormantFilter him = FormantFilter_create (my xmin, my xmax, nt, dt, t1,
fmin_hz, fmax_hz, nf, df_hz, f1_hz);
// Temporary objects
autoSound sframe = Sound_createSimple (1, windowDuration, samplingFrequency);
autoSound window = Sound_createGaussian (windowDuration, samplingFrequency);
autoMelderProgress progress (L"Sound & Pitch: To FormantFilter");
for (long i = 1; i <= nt; i++) {
double t = Sampled_indexToX (him.peek(), i);
double b, f0 = Pitch_getValueAtTime (thee, t, kPitch_unit_HERTZ, 0);
if (f0 == NUMundefined || f0 == 0) {
f0_undefined++; f0 = f0_median;
}
b = relative_bw * f0;
Sound_into_Sound (me, sframe.peek(), t - windowDuration / 2);
Sounds_multiply (sframe.peek(), window.peek());
Sound_into_FormantFilter_frame (sframe.peek(), him.peek(), i, b);
if ( (i % 10) == 1) {
Melder_progress ( (double) i / nt, L"Frame ", Melder_integer (i), L" out of ",
Melder_integer (nt), L".");
}
}
double ref = FilterBank_DBREF * gaussian_window_squared_correction (window -> nx);
NUMdmatrix_to_dBs (his z, 1, his ny, 1, his nx, ref, FilterBank_DBFAC, FilterBank_DBFLOOR);
return him.transfer();
} catch (MelderError) {
Melder_throw ("FormantFilter not created from Pitch & FormantFilter.");
}
}
autoSpectrogram Sound_and_Pitch_to_Spectrogram (Sound me, Pitch thee, double analysisWidth, double dt, double f1_hz, double fmax_hz, double df_hz, double relative_bw) {
try {
double t1, windowDuration = 2.0 * analysisWidth; /* gaussian window */
double nyquist = 0.5 / my dx, samplingFrequency = 2.0 * nyquist, fmin_hz = 0.0;
long numberOfFrames, f0_undefined = 0.0;
if (my xmin > thy xmin || my xmax > thy xmax) Melder_throw
(U"The domain of the Sound is not included in the domain of the Pitch.");
double f0_median = Pitch_getQuantile (thee, thy xmin, thy xmax, 0.5, kPitch_unit_HERTZ);
if (f0_median == NUMundefined || f0_median == 0.0) {
f0_median = 100.0;
Melder_warning (U"Pitch values undefined. Bandwith fixed to 100 Hz. ");
}
if (f1_hz <= 0.0) {
f1_hz = 100.0;
}
if (fmax_hz <= 0.0) {
fmax_hz = nyquist;
}
if (df_hz <= 0.0) {
df_hz = f0_median / 2.0;
}
if (relative_bw <= 0.0) {
relative_bw = 1.1;
}
fmax_hz = MIN (fmax_hz, nyquist);
long numberOfFilters = lround ( (fmax_hz - f1_hz) / df_hz);
Sampled_shortTermAnalysis (me, windowDuration, dt, &numberOfFrames, &t1);
autoSpectrogram him = Spectrogram_create (my xmin, my xmax, numberOfFrames, dt, t1, fmin_hz, fmax_hz, numberOfFilters, df_hz, f1_hz);
// Temporary objects
autoSound sframe = Sound_createSimple (1, windowDuration, samplingFrequency);
autoSound window = Sound_createGaussian (windowDuration, samplingFrequency);
autoMelderProgress progress (U"Sound & Pitch: To FormantFilter");
for (long iframe = 1; iframe <= numberOfFrames; iframe++) {
double t = Sampled_indexToX (him.get(), iframe);
double b, f0 = Pitch_getValueAtTime (thee, t, kPitch_unit_HERTZ, 0);
if (f0 == NUMundefined || f0 == 0.0) {
f0_undefined ++;
f0 = f0_median;
}
b = relative_bw * f0;
Sound_into_Sound (me, sframe.get(), t - windowDuration / 2.0);
Sounds_multiply (sframe.get(), window.get());
Sound_into_Spectrogram_frame (sframe.get(), him.get(), iframe, b);
if (iframe % 10 == 1) {
Melder_progress ( (double) iframe / numberOfFrames, U"Frame ", iframe, U" out of ",
numberOfFrames, U".");
}
}
_Spectrogram_windowCorrection (him.get(), window -> nx);
return him;
} catch (MelderError) {
Melder_throw (U"FormantFilter not created from Pitch & FormantFilter.");
}
}
