/// <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 }