В поисках дискретного косинусного преобразования 8x8 (или nxn) псевдокода IDCT [закрыто

StackOverflow https://stackoverflow.com/questions/8363753

  •  27-10-2019
  •  | 
  •  

Вопрос

Я некоторое время искал Google, чтобы найти псевдо-код прилично эффективного алгоритма DCT 8x8 (или NXN), и я ничего не могу найти!

Я реализовал наивный подход, и это заняло слишком много времени, чтобы выполнить.

Если бы вы могли опубликовать немного псевдокода или ссылаться на хорошую книгу/документ/веб-сайт, это было бы полезно.

Примеры C или C ++ еще лучше!

Это было полезно?

Решение

Как запрошено в комментариях, источник (будьте слегка осторожен, это в C#, но разница с C ++ должна быть минимальной, и да, я знаю, что код хромает):

Главный цикл (a = результат, b = вход):

for (int y = 0; y < 8; y++)
{
    for (int x = 0; x < 8; x++)
    {
        A[y * 8 + x] = 0;
        for (int u = 0; u < 8; u++)
            for (int v = 0; v < 8; v++)
                A[y * 8 + x] += alpha(u) * alpha(v) * B[u, v] *
                    cosine[u, x] * cosine[v, y];
    }
}

Поддержка:

static double alpha(int i)
{
    if (i == 0)
        return SQRT2o2 * 0.5;
    return 0.5;
}
const double SQRT2o2 = 1.414213562373095048801688724209 * 0.5;
cosine = new double[8, 8];
const double inv16 = 1.0 / 16.0;
for (int i = 0; i < 8; i++)
{
     for (int j = 0; j < 8; j++)
     {
         cosine[j, i] = Math.Cos(Math.PI * j * (2.0 * i + 1) * inv16);
     }
}

РЕДАКТИРОВАТЬ: Я рассчитывал это - на 512 на 512 пикселей (один канал). Это займет полсекунды. Конечно, это медленно, но нигде не «навсегда».

Другие советы

FFTW имеет эффективную реализацию с открытым исходным кодом

Лицензировано под: CC-BY-SA с атрибуция
Не связан с StackOverflow
scroll top