public double[] test2000A(string filenavn)
{
filename = filenavn;
// FileSystemAccess _fileSystem = new FileSystemAccess();
//_filePathFig = filename.Substring(filename.Length - 4,filename.Length-1);
int længde = filename.Length - 4;
_filePathFig = filename.Substring(0, længde) + ".png";
_filePathtxt = filename.Substring(0, længde) + ".txt";
_filePathwav = filename.Substring(0, længde) + ".wav";
(double[] audio, int sampleRate) = ReadWAV(_filePathwav);
// int fftSize1 = 16384;
int targetWidthPx = 300;
int fftSize1 = (int)MathF.Pow(2, 5);
int stepSize1 = 80;
//int stepSize1 = audio.Length / targetWidthPx;
var sg = new SpectrogramGenerator(sampleRate, fftSize: fftSize1, stepSize: stepSize1, maxFreq: 1000);
sg.Add(audio);
double[] testData = TakeSum(sg.GetFFTs());
sg.SaveImage(_filePathFig, intensity: 8, dB: true, roll: true);
line = line.Replace(',', '.');
return(testData);
}
public void Test_SaveEmpty_Throws()
{
(double[] audio, int sampleRate) = AudioFile.ReadWAV("../../../../../data/cant-do-that-44100.wav");
int fftSize = 4096;
var spec = new SpectrogramGenerator(sampleRate, fftSize, stepSize: 500);
//spec.Add(audio);
Assert.Throws <InvalidOperationException>(() => { spec.SaveImage("empty.png"); });
}
public void Test_Quickstart_Hal()
{
(double[] audio, int sampleRate) = AudioFile.ReadWAV("../../../../../data/cant-do-that-44100.wav");
var sg = new SpectrogramGenerator(sampleRate, fftSize: 4096, stepSize: 500, maxFreq: 3000);
sg.Add(audio);
sg.SaveImage("../../../../../dev/graphics/hal.png");
Console.WriteLine(sg);
}
public void Test_Readme_HeaderImage()
{
(double[] audio, int sampleRate) = AudioFile.ReadWAV("../../../../../data/cant-do-that-44100.wav");
int fftSize = 2048;
var spec = new SpectrogramGenerator(sampleRate, fftSize, stepSize: 400, maxFreq: 6000);
spec.Add(audio);
spec.SaveImage("../../../../../dev/graphics/hal-spectrogram.png", intensity: 10, dB: true, dBScale: .05);
Console.WriteLine(spec);
}
public void Test_AGC_off()
{
string wavFilePath = "../../../../../data/qrss-10min.wav";
(double[] audio, int sampleRate) = AudioFile.ReadWAV(wavFilePath);
int fftSize = 8192;
var spec = new SpectrogramGenerator(sampleRate, fftSize, stepSize: 2000, maxFreq: 3000);
spec.Add(audio);
spec.SaveImage("qrss-agc-off.png", intensity: 3);
}
void test()
{
FileSystemAccess _fileSystem = new FileSystemAccess();
_filePath = _fileSystem.GetCombinePath("TestH1.wav");
_filePathFig = _fileSystem.GetCombinePath("TestH1Fig.png");
(double[] audio, int sampleRate) = ReadWAV(_filePath);
var sg = new SpectrogramGenerator(sampleRate, fftSize: 4096, stepSize: 500, maxFreq: 3000);
sg.Add(audio);
sg.SaveImage(_filePathFig);
}
public void Test_Make_CommonColormaps()
{
(double[] audio, int sampleRate) = AudioFile.ReadWAV("../../../../../data/cant-do-that-44100.wav");
int fftSize = 1 << 12;
var spec = new SpectrogramGenerator(sampleRate, fftSize, stepSize: 700, maxFreq: 2000);
spec.SetWindow(FftSharp.Window.Hanning(fftSize / 3)); // sharper window than typical
spec.Add(audio);
// delete old colormap files
foreach (var filePath in System.IO.Directory.GetFiles("../../../../../dev/graphics/", "hal-*.png"))
{
System.IO.File.Delete(filePath);
}
foreach (var cmap in Colormap.GetColormaps())
{
spec.SetColormap(cmap);
spec.SaveImage($"../../../../../dev/graphics/hal-{cmap.Name}.png");
Debug.WriteLine($"");
}
}
public void Test_AGC_normWindow()
{
// strategy here is to create a weighted moving window mean and normalize to that
string wavFilePath = "../../../../../data/qrss-10min.wav";
(double[] audio, int sampleRate) = AudioFile.ReadWAV(wavFilePath);
int fftSize = 8192;
var spec = new SpectrogramGenerator(sampleRate, fftSize, stepSize: 2000, maxFreq: 3000);
spec.Add(audio);
var ffts = spec.GetFFTs();
for (int i = 0; i < ffts.Count; i++)
{
ffts[i] = SubtractMovingWindowFloor(ffts[i]);
}
spec.SaveImage("qrss-agc-norm-window.png", intensity: 3);
}
public void Test_Quickstart_Handel()
{
double[] audio = AudioFile.ReadMP3("../../../../../data/Handel - Air and Variations.mp3");
int sampleRate = 44100;
int fftSize = 16384;
int targetWidthPx = 3000;
int stepSize = audio.Length / targetWidthPx;
var sg = new SpectrogramGenerator(sampleRate, fftSize, stepSize, maxFreq: 2200);
sg.Add(audio);
sg.SaveImage("../../../../../dev/graphics/spectrogram-song.png", intensity: 5, dB: true);
Console.WriteLine(sg);
/*
* Spectrogram (2993, 817)
* Vertical (817 px): 0 - 2,199 Hz, FFT size: 16,384 samples, 2.69 Hz/px
* Horizontal (2993 px): 2.96 min, window: 0.37 sec, step: 0.06 sec, overlap: 84%
*/
}
public void Test_AGC_normToNoiseFloor()
{
// strategy here is to normalize to the magnitude of the quietest 20% of frequencies
string wavFilePath = "../../../../../data/qrss-10min.wav";
(double[] audio, int sampleRate) = AudioFile.ReadWAV(wavFilePath);
int fftSize = 8192;
var spec = new SpectrogramGenerator(sampleRate, fftSize, stepSize: 2000, maxFreq: 3000);
spec.Add(audio);
var ffts = spec.GetFFTs();
double normalIntensity = 2;
for (int i = 0; i < ffts.Count; i++)
{
double[] sorted = new double[ffts[i].Length];
ffts[i].CopyTo(sorted, 0);
Array.Sort(sorted);
double percentile = 0.25;
int percentileIndex = (int)(percentile * ffts[0].Length);
double floorValue = sorted[percentileIndex];
for (int y = 0; y < ffts[i].Length; y++)
{
ffts[i][y] = ffts[i][y] / floorValue * normalIntensity;
}
Console.WriteLine(floorValue);
}
spec.SaveImage("qrss-agc-norm-floor.png", intensity: 3);
}
public static void SaveSgAsPNG(SpectrogramGenerator sg, string filePathFig)
{
sg.SaveImage(filePathFig, intensity: 8, dB: false, roll: true);
}
