/// <summary>
/// Compress a 2D matrix
/// </summary>
/// <param name="data_input">data matrix</param>
/// <param name="level">number of wavelet levels</param>
/// <param name="threshold">threshold where all absolute values less than this is set to zero (set to 0 to avoid compression)</param>
/// <param name="lastHeight">return last height processed</param>
/// <param name="lastWidth">return last width processed</param>
public static void Compress2D(double[][] data_input, int level, int threshold, out int lastHeight, out int lastWidth)
{
int temp_level = level;
int ex_height = data_input.Length;
int ex_width = data_input[0].Length;
int temp_ex_height = ex_height;
int temp_ex_width = ex_width;
while (temp_level > 0 && ex_height > 1 && ex_width > 1)
{
HaarWaveletTransform.HaarTransform2D(data_input, ex_height, ex_width);
if (ex_width > 1)
{
ex_width = ex_width / 2;
}
if (ex_height > 1)
{
ex_height = ex_height / 2;
}
temp_level--;
}
if (threshold > 0)
{
Quantize.DataQuantize2D(data_input, temp_ex_height, temp_ex_width, threshold);
}
lastHeight = ex_height;
lastWidth = ex_width;
}
/// <summary>
/// Compress and decompress 2D matrix
/// </summary>
/// <param name="data_input">data matrix</param>
/// <param name="level">number of wavelet levels</param>
/// <param name="threshold">threshold where all absolute values less than this is set to zero (set to 0 to avoid compression)</param>
public static void CompressDecompress2D(double[][] data_input, int level, int threshold)
{
int temp_level = level;
int ex_height = data_input.Length;
int ex_width = data_input[0].Length;
int temp_ex_height = ex_height;
int temp_ex_width = ex_width;
while (temp_level > 0 && ex_height > 1 && ex_width > 1)
{
HaarWaveletTransform.HaarTransform2D(data_input, ex_height, ex_width);
if (ex_width > 1)
{
ex_width = ex_width / 2;
}
if (ex_height > 1)
{
ex_height = ex_height / 2;
}
temp_level--;
}
if (threshold > 0)
{
Quantize.DataQuantize2D(data_input, temp_ex_height, temp_ex_width, threshold);
}
while (temp_level < level && ex_height > 1 && ex_width > 1)
{
if (ex_width > 1)
{
ex_width = ex_width * 2;
}
if (ex_height > 1)
{
ex_height = ex_height * 2;
}
HaarWaveletTransform.InverseHaarTransform2D(data_input, ex_height, ex_width);
temp_level++;
}
}