/// <summary>
/// Haar Wavelet Transform and Compress
/// </summary>
/// <param name="m">matrix (logSpectrogram)</param>
/// <returns>matrix wavelet'ed</returns>
public Matrix ApplyWaveletCompression(ref Matrix m, out int lastHeight, out int lastWidth)
{
using (new DebugTimer("ApplyWaveletCompression(matrix)")) {
// Wavelet Transform
Matrix wavelet = m.Copy();
WaveletCompress.HaarTransform2D(wavelet.MatrixData, numberWaveletTransforms, out lastHeight, out lastWidth);
// Compress
Matrix waveletCompressed = wavelet.Resize(numberCoefficients, wavelet.Columns);
return(waveletCompressed);
}
}
public static Matrix GetWaveletTransformedMatrix(double[][] image, WaveletMethod waveletMethod)
{
int width = image[0].Length;
int height = image.Length;
Matrix dwtMatrix = null;
Stopwatch stopWatch = Stopwatch.StartNew();
long startS = stopWatch.ElapsedTicks;
switch (waveletMethod)
{
case WaveletMethod.Haar:
Haar.Haar2D(image, height, width);
dwtMatrix = new Matrix(image);
break;
case WaveletMethod.HaarTransformTensor: // This is using the tensor product layout
dwtMatrix = WaveletUtils.HaarWaveletTransform2D(image);
break;
case WaveletMethod.HaarWaveletDecompositionTensor: // This is using the tensor product layout
var haar = new StandardHaarWaveletDecomposition();
haar.DecomposeImageInPlace(image);
dwtMatrix = new Matrix(image);
break;
case WaveletMethod.NonStandardHaarWaveletDecomposition: // JPEG 2000
var haarNonStandard = new NonStandardHaarWaveletDecomposition();
haarNonStandard.DecomposeImageInPlace(image);
dwtMatrix = new Matrix(image);
break;
case WaveletMethod.HaarCSharp:
ForwardWaveletTransform.Transform2D(image, false, 2);
dwtMatrix = new Matrix(image);
break;
case WaveletMethod.HaarWaveletCompress:
int lastHeight = 0;
int lastWidth = 0;
WaveletCompress.HaarTransform2D(image, 10000, out lastHeight, out lastWidth);
dwtMatrix = new Matrix(image);
break;
default:
break;
}
long endS = stopWatch.ElapsedTicks;
Console.WriteLine("WaveletMethod: {0} Time in ticks: {1}", Enum.GetName(typeof(WaveletMethod), waveletMethod), (endS - startS));
// increase all values
const int mul = 50;
double[][] haarImageNormalized5k = dwtMatrix.MatrixData.Select(i => i.Select(j => j * mul).ToArray()).ToArray();
// convert to byte values (0 - 255)
// duplicate the octave/ matlab method uint8
var uint8 = new double[haarImageNormalized5k.Length][];
for (int i = 0; i < haarImageNormalized5k.Length; i++)
{
uint8[i] = new double[haarImageNormalized5k.Length];
for (int j = 0; j < haarImageNormalized5k[i].Length; j++)
{
double v = haarImageNormalized5k[i][j];
if (v > 255)
{
uint8[i][j] = 255;
}
else if (v < 0)
{
uint8[i][j] = 0;
}
else
{
uint8[i][j] = (byte)haarImageNormalized5k[i][j];
}
}
}
var uint8Matrix = new Matrix(uint8);
return(uint8Matrix);
}
