/// <summary>
/// This method compares the acoustic indices derived from two different long duration recordings of the same length.
/// It takes as input six csv files of acoustic indices in spectrogram columns, three csv files for each of the original recordings to be compared.
/// The method produces one spectrogram image files:
/// 1) A false-colour difference spectrogram, where the difference is shown as a plus/minus departure from grey.
/// </summary>
public static void DrawDifferenceSpectrogram(DirectoryInfo ipdir, FileInfo ipFileName1, FileInfo ipFileName2, DirectoryInfo opdir)
{
var cs1 = new LDSpectrogramRGB(minuteOffset, xScale, sampleRate, frameWidth, colorMap)
{
FileName = ipFileName1.Name,
ColorMode = colorMap,
BackgroundFilter = backgroundFilterCoeff,
};
string[] keys = colorMap.Split('-');
cs1.ReadCsvFiles(ipdir, ipFileName1.Name, keys);
if (cs1.GetCountOfSpectrogramMatrices() == 0)
{
LoggedConsole.WriteLine("There are no spectrogram matrices in cs1.dictionary.");
return;
}
var cs2 = new LDSpectrogramRGB(minuteOffset, xScale, sampleRate, frameWidth, colorMap)
{
FileName = ipFileName2.Name,
ColorMode = colorMap,
BackgroundFilter = backgroundFilterCoeff,
};
cs2.ReadCsvFiles(ipdir, ipFileName2.Name, keys);
if (cs2.GetCountOfSpectrogramMatrices() == 0)
{
LoggedConsole.WriteLine("There are no spectrogram matrices in cs2.dictionary.");
return;
}
//string title1 = String.Format("DIFFERENCE SPECTROGRAM ({0} - {1}) (scale:hours x kHz) (colour: R-G-B={2})", ipFileName1, ipFileName2, cs1.ColorMODE);
//Image deltaSp1 = LDSpectrogramDifference.DrawDifferenceSpectrogram(cs1, cs2, colourGain);
//Image titleBar1 = LDSpectrogramRGB.DrawTitleBarOfFalseColourSpectrogram(title1, deltaSp1.Width, SpectrogramConstants.HEIGHT_OF_TITLE_BAR);
//deltaSp1 = LDSpectrogramRGB.FrameSpectrogram(deltaSp1, titleBar1, minuteOffset, cs1.X_interval, cs1.Y_interval);
//string opFileName1 = ipFileName1 + ".Difference.COLNEG.png";
//deltaSp1.Save(Path.Combine(opdir.FullName, opFileName1));
//Draw positive difference spectrograms in one image.
double colourGain = 2.0;
Image <Rgb24>[] images = DrawPositiveDifferenceSpectrograms(cs1, cs2, colourGain);
int nyquist = cs1.SampleRate / 2;
int herzInterval = 1000;
string title =
$"DIFFERENCE SPECTROGRAM where {ipFileName1} > {ipFileName2}. (scale:hours x kHz) (colour: R-G-B={cs1.ColorMode})";
var titleBar = LDSpectrogramRGB.DrawTitleBarOfFalseColourSpectrogram(title, images[0].Width);
images[0] = LDSpectrogramRGB.FrameLDSpectrogram(images[0], titleBar, cs1, nyquist, herzInterval);
title = string.Format("DIFFERENCE SPECTROGRAM where {1} > {0} (scale:hours x kHz) (colour: R-G-B={2})", ipFileName1, ipFileName2, cs1.ColorMode);
titleBar = LDSpectrogramRGB.DrawTitleBarOfFalseColourSpectrogram(title, images[1].Width);
images[1] = LDSpectrogramRGB.FrameLDSpectrogram(images[1], titleBar, cs1, nyquist, herzInterval);
Image combinedImage = ImageTools.CombineImagesVertically(images);
string opFileName = ipFileName1 + "-" + ipFileName2 + ".Difference.png";
combinedImage.Save(Path.Combine(opdir.FullName, opFileName));
}
public static Image DrawTStatisticSpectrogramsOfSingleIndex(string key, LDSpectrogramRGB cs1, LDSpectrogramRGB cs2, double tStatThreshold)
{
var image1 = cs1.DrawGreyscaleSpectrogramOfIndex(key);
var image2 = cs2.DrawGreyscaleSpectrogramOfIndex(key);
if (image1 == null || image2 == null)
{
Console.WriteLine("WARNING: From method ColourSpectrogram.DrawTStatisticGreyscaleSpectrogramOfIndex()");
Console.WriteLine(" Null image returned with key: {0}", key);
return(null);
}
//frame image 1
int nyquist = cs1.SampleRate / 2;
int herzInterval = 1000;
string title = $"{key} SPECTROGRAM for: {cs1.FileName}. (scale:hours x kHz)";
var titleBar = LDSpectrogramRGB.DrawTitleBarOfGrayScaleSpectrogram(title, image1.Width);
image1 = LDSpectrogramRGB.FrameLDSpectrogram(image1, titleBar, cs1, nyquist, herzInterval);
//frame image 2
title = $"{key} SPECTROGRAM for: {cs2.FileName}. (scale:hours x kHz)";
titleBar = LDSpectrogramRGB.DrawTitleBarOfGrayScaleSpectrogram(title, image2.Width);
image2 = LDSpectrogramRGB.FrameLDSpectrogram(image2, titleBar, cs1, nyquist, herzInterval);
//get matrices required to calculate matrix of t-statistics
double[,] avg1 = cs1.GetSpectrogramMatrix(key);
if (key.Equals("ENT"))
{
avg1 = MatrixTools.SubtractValuesFromOne(avg1);
}
double[,] std1 = cs1.GetStandarDeviationMatrix(key);
double[,] avg2 = cs2.GetSpectrogramMatrix(key);
if (key.Equals("ENT"))
{
avg2 = MatrixTools.SubtractValuesFromOne(avg2);
}
double[,] std2 = cs2.GetStandarDeviationMatrix(key);
//draw a spectrogram of t-statistic values
//double[,] tStatMatrix = SpectrogramDifference.GetTStatisticMatrix(avg1, std1, cs1.SampleCount, avg2, std2, cs2.SampleCount);
//Image image3 = SpectrogramDifference.DrawTStatisticSpectrogram(tStatMatrix);
//titleBar = SpectrogramDifference.DrawTitleBarOfTStatisticSpectrogram(cs1.BaseName, cs2.BaseName, image1.Width, titleHt);
//image3 = ColourSpectrogram.FrameSpectrogram(image3, titleBar, minOffset, cs2.X_interval, cs2.Y_interval);
//draw a difference spectrogram derived from by thresholding a t-statistic matrix
var image4 = DrawDifferenceSpectrogramDerivedFromSingleTStatistic(key, cs1, cs2, tStatThreshold, ColourGain);
title = string.Format("{0} DIFFERENCE SPECTROGRAM (thresholded by t-statistic={3}) for: {1} - {2}. (scale:hours x kHz)", key, cs1.FileName, cs2.FileName, tStatThreshold);
titleBar = LDSpectrogramRGB.DrawTitleBarOfGrayScaleSpectrogram(title, image2.Width);
image4 = LDSpectrogramRGB.FrameLDSpectrogram(image4, titleBar, cs2, nyquist, herzInterval);
var combinedImage = ImageTools.CombineImagesVertically(image1, image2, image4);
return(combinedImage);
}
public static Image DrawTStatisticSpectrogramsOfMultipleIndices(LDSpectrogramRGB cs1, LDSpectrogramRGB cs2, double tStatThreshold, double colourGain)
{
string[] keys = cs1.ColorMap.Split('-'); //assume both spectorgrams have the same acoustic indices in same order
double[,] m1 = GetDifferenceSpectrogramDerivedFromSingleTStatistic(keys[0], cs1, cs2, tStatThreshold);
double[,] m2 = GetDifferenceSpectrogramDerivedFromSingleTStatistic(keys[1], cs1, cs2, tStatThreshold);
double[,] m3 = GetDifferenceSpectrogramDerivedFromSingleTStatistic(keys[2], cs1, cs2, tStatThreshold);
int rows = m1.GetLength(0); //number of rows
int cols = m1.GetLength(1); //number
var spg1Image = new Image <Rgb24>(cols, rows);
var spg2Image = new Image <Rgb24>(cols, rows);
int maxRgbValue = 255;
for (int row = 0; row < rows; row++)
{
for (int column = 0; column < cols; column++)
{
var dR = m1[row, column] * colourGain;
var dG = m2[row, column] * colourGain;
var dB = m3[row, column] * colourGain;
byte iR1 = 0;
byte iR2 = 0;
byte iG1 = 0;
byte iG2 = 0;
byte iB1 = 0;
byte iB2 = 0;
var value = Convert.ToByte(Math.Abs(dR) * maxRgbValue);
value = (byte)Math.Min(maxRgbValue, value);
if (dR > 0.0)
{
iR1 = value;
}
else
{
iR2 = value;
}
value = Convert.ToByte(Math.Abs(dG) * maxRgbValue);
value = (byte)Math.Min(maxRgbValue, value);
if (dG > 0.0)
{
iG1 = value;
}
else
{
iG2 = value;
}
value = Convert.ToByte(Math.Abs(dB) * maxRgbValue);
value = (byte)Math.Min(maxRgbValue, value);
if (dB > 0.0)
{
iB1 = value;
}
else
{
iB2 = value;
}
var colour1 = Color.FromRgb(iR1, iG1, iB1);
var colour2 = Color.FromRgb(iR2, iG2, iB2);
spg1Image[column, row] = colour1;
spg2Image[column, row] = colour2;
}
}
var images = new Image <Rgb24> [2];
int nyquist = cs1.SampleRate / 2;
int herzInterval = 1000;
string title = string.Format("DIFFERENCE SPECTROGRAM (thresholded by t-Statistic={2}) where {0} > {1} (scale:hours x kHz) (colour: R-G-B={2})", cs1.FileName, cs2.FileName, tStatThreshold);
var titleBar = LDSpectrogramRGB.DrawTitleBarOfFalseColourSpectrogram(title, spg1Image.Width);
images[0] = LDSpectrogramRGB.FrameLDSpectrogram(spg1Image, titleBar, cs1, nyquist, herzInterval);
title = string.Format("DIFFERENCE SPECTROGRAM (thresholded by t-Statistic={2}) where {1} > {0} (scale:hours x kHz) (colour: R-G-B={2})", cs1.FileName, cs2.FileName, tStatThreshold);
titleBar = LDSpectrogramRGB.DrawTitleBarOfFalseColourSpectrogram(title, spg2Image.Width);
images[1] = LDSpectrogramRGB.FrameLDSpectrogram(spg2Image, titleBar, cs1, nyquist, herzInterval);
var compositeImage = ImageTools.CombineImagesVertically(images);
return(compositeImage);
}
/// <summary>
/// This method compares the acoustic indices derived from two different long duration recordings of the same length.
/// It takes as input any number of csv files of acoustic indices in spectrogram columns.
/// Typically there will be at least three indices csv files for each of the original recordings to be compared.
/// The method produces four spectrogram image files:
/// 1) A negative false-color spectrogram derived from the indices of recording 1.
/// 2) A negative false-color spectrogram derived from the indices of recording 2.
/// 3) A spectrogram of euclidean distances between the two input files.
/// 4) The above three spectrograms combined in one image.
/// </summary>
public static void DrawDistanceSpectrogram(
DirectoryInfo inputDirectory,
FileInfo inputFileName1,
FileInfo inputFileName2,
DirectoryInfo outputDirectory)
{
// PARAMETERS
string outputFileName1 = inputFileName1.Name;
var cs1 = new LDSpectrogramRGB(minuteOffset, xScale, sampleRate, frameWidth, colorMap);
cs1.ColorMode = colorMap;
cs1.BackgroundFilter = backgroundFilterCoeff;
string[] keys = colorMap.Split('-');
cs1.ReadCsvFiles(inputDirectory, inputFileName1.Name, keys);
double blueEnhanceParameter = 0.0;
cs1.DrawNegativeFalseColorSpectrogram(outputDirectory, outputFileName1, blueEnhanceParameter);
string imagePath = Path.Combine(outputDirectory.FullName, outputFileName1 + ".COLNEG.png");
var spg1Image = Image.Load <Rgb24>(imagePath);
if (spg1Image == null)
{
LoggedConsole.WriteLine("SPECTROGRAM IMAGE DOES NOT EXIST: {0}", imagePath);
return;
}
int nyquist = cs1.SampleRate / 2;
int hertzInterval = 1000;
string title =
$"FALSE COLOUR SPECTROGRAM: {inputFileName1}. (scale:hours x kHz) (colour: R-G-B={cs1.ColorMode})";
var titleBar = LDSpectrogramRGB.DrawTitleBarOfFalseColourSpectrogram(title, spg1Image.Width);
spg1Image = LDSpectrogramRGB.FrameLDSpectrogram(
spg1Image,
titleBar,
cs1,
nyquist,
hertzInterval);
string outputFileName2 = inputFileName2.Name;
var cs2 = new LDSpectrogramRGB(minuteOffset, xScale, sampleRate, frameWidth, colorMap)
{
ColorMode = colorMap,
BackgroundFilter = backgroundFilterCoeff,
};
cs2.ReadCsvFiles(inputDirectory, inputFileName2.Name, keys);
// cs2.DrawGreyScaleSpectrograms(opdir, opFileName2);
cs2.DrawNegativeFalseColorSpectrogram(outputDirectory, outputFileName2, blueEnhanceParameter);
imagePath = Path.Combine(outputDirectory.FullName, outputFileName2 + ".COLNEG.png");
var spg2Image = Image.Load <Rgb24>(imagePath);
if (spg2Image == null)
{
LoggedConsole.WriteLine("SPECTROGRAM IMAGE DOES NOT EXIST: {0}", imagePath);
return;
}
title =
$"FALSE COLOUR SPECTROGRAM: {inputFileName2}. (scale:hours x kHz) (colour: R-G-B={cs2.ColorMode})";
titleBar = LDSpectrogramRGB.DrawTitleBarOfFalseColourSpectrogram(title, spg2Image.Width);
spg2Image = LDSpectrogramRGB.FrameLDSpectrogram(
spg2Image,
titleBar,
cs1,
nyquist,
hertzInterval);
string outputFileName4 = inputFileName1 + ".EuclideanDistance.png";
var deltaSp = DrawDistanceSpectrogram(cs1, cs2);
Color[] colorArray = LDSpectrogramRGB.ColourChart2Array(GetDifferenceColourChart());
titleBar = DrawTitleBarOfEuclidianDistanceSpectrogram(
inputFileName1.Name,
inputFileName2.Name,
colorArray,
deltaSp.Width,
SpectrogramConstants.HEIGHT_OF_TITLE_BAR);
deltaSp = LDSpectrogramRGB.FrameLDSpectrogram(deltaSp, titleBar, cs2, nyquist, hertzInterval);
deltaSp.Save(Path.Combine(outputDirectory.FullName, outputFileName4));
string outputFileName5 = inputFileName1 + ".2SpectrogramsAndDistance.png";
var combinedImage = ImageTools.CombineImagesVertically(spg1Image, spg2Image, deltaSp);
combinedImage.Save(Path.Combine(outputDirectory.FullName, outputFileName5));
}
