c# - 将小波滤波器添加到离散 Haar 小波变换 (DWT)

Question

我需要使用 c# 获得 1 级 5/3 离散 Haar 小波变换 (DWT) 源代码的帮助。

我使用这个项目，前向小波变换的方法在这里：

FWT(double[] data)
    {
        int h = data.Length >> 1;
        for (int i = 0; i < h; i++)
        {
            int k = (i << 1);
            temp[i] = data[k] * s0 + data[k + 1] * s1;
            temp[i + h] = data[k] * w0 + data[k + 1] * w1;
        }
    }
FWT(double[,] data)
{
for (int k = 0; k < 1; k++)
{
    for (int i = 0; i < rows / (k+1); i++)
    {
        for (int j = 0; j < row.Length / (k+1); j++)
            row[j] = data[i, j];

        FWT(row);

        for (int j = 0; j < row.Length / (k+1); j++)
            data[i, j] = row[j];
    }

    for (int j = 0; j < cols / (k+1); j++)
    {
        for (int i = 0; i < col.Length / (k+1); i++)
            col[i] = data[i, j];

        FWT(col);

        for (int i = 0; i < col.Length / (k+1); i++)
            data[i, j] = col[i];
    }
}

}

w0 = 0.5；w1 = -0.5;s0 = 0.5;s1 = 0.5;

我在论文中搜索了这个主题，但我不了解 5/3 或 9/7 小波滤波器的算法，如何更改此代码？任何帮助将非常感激

Answer 1

您可能会发现它很有用：在 pdf.js 中实现 jpeg2000 解码器。

5-3代码的核心实现：

function reversibleTransformFilter(x, offset, length) {
    var len = length >> 1;
    offset = offset | 0;
    var j, n;

    for (j = offset, n = len + 1; n--; j += 2) {
        x[j] -= (x[j - 1] + x[j + 1] + 2) >> 2;
    }

    for (j = offset + 1, n = len; n--; j += 2) {
        x[j] += (x[j - 1] + x[j + 1]) >> 1;
    }
};