/// <summary>
/// Use this method when want to match defined shape in target using cross-correlation.
/// This was the method used by Stewart Gage.
/// First set target and source to same dynamic range.
/// Then NormaliseMatrixValues target and source to unit-length.
/// </summary>
public static Tuple <double[]> Execute_StewartGage(double[,] target, double dynamicRange, SpectrogramStandard sonogram,
List <AcousticEvent> segments, int minHz, int maxHz, double minDuration)
{
Log.WriteLine("SEARCHING FOR EVENTS LIKE TARGET.");
if (segments == null)
{
return(null);
}
int minBin = (int)(minHz / sonogram.FBinWidth);
int maxBin = (int)(maxHz / sonogram.FBinWidth);
int targetLength = target.GetLength(0);
//adjust target's dynamic range to that set by user
target = SNR.SetDynamicRange(target, 0.0, dynamicRange); //set event's dynamic range
double[] v1 = DataTools.Matrix2Array(target);
v1 = DataTools.normalise2UnitLength(v1);
//var image = BaseSonogram.Data2ImageData(target);
//ImageTools.DrawMatrix(image, 1, 1, @"C:\SensorNetworks\Output\FELT_Currawong\target.png");
double[] scores = new double[sonogram.FrameCount];
foreach (AcousticEvent av in segments)
{
Log.WriteLine("SEARCHING SEGMENT.");
int startRow = (int)Math.Round(av.TimeStart * sonogram.FramesPerSecond);
int endRow = (int)Math.Round(av.TimeEnd * sonogram.FramesPerSecond);
if (endRow >= sonogram.FrameCount)
{
endRow = sonogram.FrameCount;
}
int stopRow = endRow - targetLength;
if (stopRow <= startRow)
{
stopRow = startRow + 1; //want minimum of one row
}
int offset = targetLength / 2;
for (int r = startRow; r < stopRow; r++)
{
double[,] matrix = DataTools.Submatrix(sonogram.Data, r, minBin, r + targetLength - 1, maxBin);
matrix = SNR.SetDynamicRange(matrix, 0.0, dynamicRange); //set event's dynamic range
//var image = BaseSonogram.Data2ImageData(matrix);
//ImageTools.DrawMatrix(image, 1, 1, @"C:\SensorNetworks\Output\FELT_CURLEW\compare.png");
double[] v2 = DataTools.Matrix2Array(matrix);
v2 = DataTools.normalise2UnitLength(v2);
scores[r] = DataTools.DotProduct(v1, v2); //the Cross Correlation
} // end of rows in segment
} // foreach (AcousticEvent av in segments)
var tuple = Tuple.Create(scores);
return(tuple);
}
public static Tuple <double[]> Execute_MFCC_XCOR(double[,] target, double dynamicRange, SpectrogramStandard sonogram,
List <AcousticEvent> segments, int minHz, int maxHz, double minDuration)
{
Log.WriteLine("SEARCHING FOR EVENTS LIKE TARGET.");
if (segments == null)
{
return(null);
}
int minBin = (int)(minHz / sonogram.FBinWidth);
int maxBin = (int)(maxHz / sonogram.FBinWidth);
int targetLength = target.GetLength(0);
//set up the matrix of cosine coefficients
int coeffCount = 12; //only use first 12 coefficients.
int binCount = target.GetLength(1); //number of filters in filter bank
double[,] cosines = MFCCStuff.Cosines(binCount, coeffCount + 1); //set up the cosine coefficients
//adjust target's dynamic range to that set by user
target = SNR.SetDynamicRange(target, 3.0, dynamicRange); //set event's dynamic range
target = MFCCStuff.Cepstra(target, coeffCount, cosines);
double[] v1 = DataTools.Matrix2Array(target);
v1 = DataTools.normalise2UnitLength(v1);
string imagePath2 = @"C:\SensorNetworks\Output\FELT_Currawong\target.png";
var result1 = BaseSonogram.Data2ImageData(target);
var image = result1.Item1;
ImageTools.DrawMatrix(image, 1, 1, imagePath2);
double[] scores = new double[sonogram.FrameCount];
foreach (AcousticEvent av in segments)
{
Log.WriteLine("SEARCHING SEGMENT.");
int startRow = (int)Math.Round(av.TimeStart * sonogram.FramesPerSecond);
int endRow = (int)Math.Round(av.TimeEnd * sonogram.FramesPerSecond);
if (endRow >= sonogram.FrameCount)
{
endRow = sonogram.FrameCount - 1;
}
endRow -= targetLength;
if (endRow <= startRow)
{
endRow = startRow + 1; //want minimum of one row
}
for (int r = startRow; r < endRow; r++)
{
double[,] matrix = DataTools.Submatrix(sonogram.Data, r, minBin, r + targetLength - 1, maxBin);
matrix = SNR.SetDynamicRange(matrix, 3.0, dynamicRange); //set event's dynamic range
//string imagePath2 = @"C:\SensorNetworks\Output\FELT_Gecko\compare.png";
//var image = BaseSonogram.Data2ImageData(matrix);
//ImageTools.DrawMatrix(image, 1, 1, imagePath2);
matrix = MFCCStuff.Cepstra(matrix, coeffCount, cosines);
double[] v2 = DataTools.Matrix2Array(matrix);
v2 = DataTools.normalise2UnitLength(v2);
double crossCor = DataTools.DotProduct(v1, v2);
scores[r] = crossCor;
} //end of rows in segment
} //foreach (AcousticEvent av in segments)
var tuple = Tuple.Create(scores);
return(tuple);
}