/// <summary>
/// Perform an inverse decompressed haar wavelet transform. E.g. perform an ihaar2d and return logSpectrogram
/// </summary>
/// <param name="wavelet">wavelet matrix</param>
/// <returns>matrix inverse wavelet'ed (e.g. logSpectrogram)</returns>
public Matrix InverseWaveletCompression(ref Matrix wavelet, int firstHeight, int firstWidth, int rows, int columns)
{
using (new DebugTimer("InverseWaveletCompression(wavelet)")) {
// Resize, e.g. Uncompress
wavelet = wavelet.Resize(rows, columns);
// 6. Perform the Inverse Wavelet Transform
Matrix m = wavelet.Copy();
WaveletDecompress.Decompress2D(m.MatrixData, numberWaveletTransforms, firstHeight, firstWidth);
return(m);
}
}
public static void TestHaarInputOutput(string imageInPath)
{
var image = ReadImageGrayscale(imageInPath);
#region Wavelet Compress Methods
// Test HaarWaveletTransform
var inputMatrix = new Matrix(image);
inputMatrix.WriteCSV("haar-before.csv");
HaarWaveletTransform.HaarTransform2D(image, inputMatrix.Rows, inputMatrix.Columns);
Matrix haarMatrixInverse = inputMatrix.Copy();
HaarWaveletTransform.InverseHaarTransform2D(haarMatrixInverse.MatrixData, haarMatrixInverse.Rows, haarMatrixInverse.Columns);
haarMatrixInverse.WriteCSV("haar-after.csv");
haarMatrixInverse.DrawMatrixImage("haar-transform-forward-and-backward.png", -1, -1, false);
// Test Wavelet Compress and Decompress in one step
const int compLevels = 8;
const int compTreshold = 150;
Matrix haarMatrixCompDecomp = haarMatrixInverse.Copy();
WaveletComDec.CompressDecompress2D(haarMatrixCompDecomp.MatrixData, compLevels, compTreshold);
haarMatrixCompDecomp.DrawMatrixImage("haar-compress-and-decompress-combined.png", -1, -1, false);
// Test Compress and Decompress in two steps
int lastHeight = 0;
int lastWidth = 0;
Matrix haarMatrixComp = haarMatrixInverse.Copy();
WaveletCompress.Compress2D(haarMatrixComp.MatrixData, compLevels, compTreshold, out lastHeight, out lastWidth);
WaveletDecompress.Decompress2D(haarMatrixComp.MatrixData, compLevels, lastHeight, lastWidth);
haarMatrixComp.DrawMatrixImage("haar-compress-and-decompress.png", -1, -1, false);
#endregion
#region Test using HaarCSharp
// Test HaarCSharp using iterations
const int haarCSharpIterations = 3;
Matrix haarMatrixCSharp = haarMatrixInverse.Copy();
ForwardWaveletTransform.Transform2D(haarMatrixCSharp.MatrixData, false, haarCSharpIterations);
haarMatrixCSharp.DrawMatrixImage("haar-forward.png", -1, -1, false);
InverseWaveletTransform.Transform2D(haarMatrixCSharp.MatrixData, false, haarCSharpIterations);
haarMatrixCSharp.DrawMatrixImage("haar-inverse.png", -1, -1, false);
// Test HaarCSharp using all levels and only 1 iteration
Matrix haarMatrixCSharpAll = haarMatrixInverse.Copy();
ForwardWaveletTransform.Transform2D(haarMatrixCSharpAll.MatrixData, true, 1);
haarMatrixCSharpAll.DrawMatrixImage("haar-forward-all.png", -1, -1, false);
InverseWaveletTransform.Transform2D(haarMatrixCSharpAll.MatrixData, true, 1);
haarMatrixCSharpAll.DrawMatrixImage("haar-inverse-all.png", -1, -1, false);
#endregion
}
