Question

我需要c＃的1级5/3离散Haar小波变换（DWT）源代码的帮助。

我使用this project，前向小波变换的方法在这里：

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

您可能会发现它很有用：implementation of jpeg2000 decoder in pdf.js。

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

将小波滤波器添加到离散Haar小波变换（DWT）

1 个答案: