我一直在谷歌搜索一段时间来找到效率相当高的8x8(或nxn)DCT算法的伪代码,我找不到任何东西!
我实施了天真的方法,执行时间太长了。
如果您可以发布一些伪代码或引用一本好书/文档/网站,那将会有所帮助。
C或C ++示例会更好!
答案 0 :(得分:4)
根据评论中的要求,来源(稍微警告,它在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 x 512像素(单通道),需要半秒钟。当然这很慢,但远不是“永远”。
答案 1 :(得分:1)
FFTW具有开源高效实施