自研BMP转JPEG编码器时DCT与量化环节是否存在实现问题？

你的代码存在3个明确错误，按影响权重从高到低排序如下：

• 缺失DCT前必需的电平偏移步骤（核心错误）
  JPEG标准明确要求，8x8像素块送入DCT变换前，必须将每个像素值从[0,255]的无符号范围偏移到[-128, 127]的有符号范围，即每个像素值减128。你当前代码直接传入原始像素值计算DCT，直流分量结果会完全偏离正确值，这是输出图像和原图差异极大的首要原因。
• 输入块存在隐式类型转换错误
  你定义的块入参为char*类型：如果编译环境下char默认为有符号类型，像素值大于127时会被解析为负数，直接参与浮点运算会得到完全错误的DCT结果。建议将入参改为unsigned char*，或转double时显式和0xFF做与运算，保证取值正确。
• 量化环节直接截断未做四舍五入
  你当前量化逻辑直接将浮点除法结果强转为short，属于向零截断，会引入最大±1的量化误差。这个误差本身不会导致图像整体错乱，但会提升块效应明显程度，正确实现需要对量化结果做四舍五入处理。

关于你问的DCT实现精度问题：

行列分别执行一维DCT的实现方案是完全符合标准的，只要dct_tbl（正向DCT系数矩阵）和dct_tblt（DCT矩阵的转置）系数计算正确，用double类型存储中间计算结果的精度完全足够，不会造成肉眼可见的图像失真。
注意校验DCT系数的正确性，标准8x8 DCT系数计算公式为：

// c[u]为归一化系数
c[0] = sqrt(1.0/8);
for (int u = 1; u < 8; u++) c[u] = sqrt(2.0/8);
// 矩阵元素计算
dct_tbl[u*8 + v] = c[u] * cos( (2*v + 1) * u * M_PI / 16.0 );

如果系数计算错误、或转置矩阵和原矩阵不匹配，才会导致DCT结果异常。

量化环节的修正代码参考：

// 替换原有直接强转的逻辑
double q_res = dblock_2[i] / (double)std_quan_tabl[i];
data[zigzag[i]] = (short)(q_res >= 0 ? q_res + 0.5 : q_res - 0.5);

参考素材 #

• 编码生成JPEG文件
• 原始BMP文件

问题关联代码 #

void FdctAndQuant(char *block,u8 const* std_quan_tabl,short *data)
{
   double dblock_1[64];
   double dblock_2[64];

    memset(dblock_1, 0, sizeof(dblock_1));
    memset(dblock_2, 0, sizeof(dblock_2));


    for (int i = 0; i < 8; i++) {
        for (int j = 0; j < 8; j++) {
            for (int k = 0; k < 8; k++)
                dblock_1[i * 8 + j] += (dct_tbl[i * 8 + k] * (double)block[8 * k + j]);
        }
    }

    for (int i = 0; i < 8; i++) {
        for (int j = 0; j < 8; j++) {
            for (int k = 0; k < 8; k++)
                dblock_2[i * 8 + j] +=  (dblock_1[i * 8 + k] * dct_tblt[8 * k + j]);
        }
    }

//Zig-zig
    for (int i = 0; i < 64; i++) {
        data[zigzag[i]] = (short)(dblock_2[i]/((double)(std_quan_tabl[i])));
    }

}

内容的提问来源于stack exchange，提问作者TIger_zh