public void Analyze()
{
// F0 estimate
var harvest = new Harvest();
harvest.FramePeriod = FramePeriod;
harvest.F0Floor = 40.0;
var f0Length = harvest.GetSamplesForHarvest(SampleRate, Wave.Length, FramePeriod);
F0 = new double[f0Length];
TimeAxis = new double[f0Length];
harvest.Estimate(Wave, SampleRate, TimeAxis, F0);
// spectral envelope estimate
var cheapTrick = new CheapTrick(SampleRate);
cheapTrick.F0Floor = 71.0;
cheapTrick.FFTSize = cheapTrick.GetFFTSizeForCheapTrick(SampleRate);
FFTSize = cheapTrick.FFTSize;
Spectrogram = Enumerable.Range(0, f0Length).Select((i) => new double[FFTSize / 2 + 1]).ToArray();
cheapTrick.Estimate(Wave, SampleRate, TimeAxis, F0, Spectrogram);
// aperiodicity estimate
var d4c = new D4C();
d4c.Threshold = 0.85;
Aperiodicity = Enumerable.Range(0, f0Length).Select((i) => new double[FFTSize / 2 + 1]).ToArray();
d4c.Estimate(Wave, SampleRate, TimeAxis, F0, f0Length, FFTSize, Aperiodicity);
}
public static void Main(string[] args)
{
var audioPath = @"vaiueo2d.wav";
var audioOutPath = @"vaiueo2d-out.wav";
var(x, fs, nbit) = AudioIO.WavRead(audioPath);
var(f0, t) = Dio.Compute(x, fs);
f0 = StoneMask.Compute(x, fs, t, f0);
var fftSize = CheapTrick.GetFftSize(fs);
var sp = CheapTrick.Compute(x, fs, t, f0, fftSize: fftSize);
var ap = D4C.Compute(x, fs, t, f0, fftSize: fftSize);
var ndim = 60;
var mgc = Codec.CodeSpectralEnvelope(sp, fs, ndim);
var bap = Codec.CodeAperiodicity(ap, fs);
Console.WriteLine($"{audioPath}:");
Console.WriteLine($" input samples count: {x.Length}");
Console.WriteLine($" sampling rate: {fs}");
Console.WriteLine($" bit rate: {nbit}");
Console.WriteLine();
Console.WriteLine($" frame count: {f0.Length}");
Console.WriteLine($" fft size: {fftSize}");
Console.WriteLine($" sp width: {sp.GetLength(1)}");
Console.WriteLine();
Console.WriteLine($" mgc width: {ndim}");
Console.WriteLine($" bap width: {bap.GetLength(1)}");
Console.WriteLine();
for (int i = 0; i < f0.Length; i++)
{
f0[i] *= 1.6789;
}
sp = Codec.DecodeSpectralEnvelope(mgc, fs, fftSize);
ap = Codec.DecodeAperiodicity(bap, fs, fftSize);
var y = Synthesis.Compute(f0, sp, ap, fs);
Console.WriteLine($"--> {audioOutPath}");
AudioIO.WavWrite(y, fs, nbit, audioOutPath);
}