如何用libjpeg-turbo实现YCCK/CMYK到BGRA的颜色空间转换？

一、libjpeg-turbo直接转换YCCK到BGRA #

新版本libjpeg-turbo（2.0及以上版本）支持直接将YCCK格式JPEG转换为BGRA，无需手动中转CMYK。只需在解压缩初始化阶段设置输出色彩空间为JCS_EXT_BGRA，库会自动完成色彩空间转换：

// 在jpeg_start_decompress之前设置输出格式
jpegDecompressWrapper.cinfo().out_color_space = JCS_EXT_BGRA;

// 后续正常执行解压缩流程
jpeg_start_decompress(&cinfo);
// 读取扫描线时，buffer中的数据已经是BGRA格式，可直接复制到输出缓冲区

如果你的libjpeg-turbo版本较低，建议升级到2.0+版本以简化流程、避免手动转换的误差。

二、手动CMYK转RGB算法优化（解决过暗问题） #

若无法升级库，当前算法过暗的核心原因是整数除法精度丢失，以及可能存在的分量对应错误，可通过以下方式优化：

1. 修正分量对应（关键前提） #

首先确认libjpeg-turbo输出的CMYK分量顺序是否为C(青)、M(品红)、Y(黄)、K(黑)。若YCCK转CMYK时分量映射错误，会直接导致色彩失真。可通过测试单一颜色（比如纯黑、纯白）验证分量对应关系是否正确。

2. 浮点运算+Gamma校正优化转换公式 #

改用浮点运算保留精度，并添加sRGB标准的Gamma校正（匹配多数显示设备的色彩空间）：

for (std::size_t i{ 0U }; i < rowStride / 4; ++i)
{
    auto const c = static_cast<std::uint8_t>(buffer[0][i * 4U]);
    auto const m = static_cast<std::uint8_t>(buffer[0][i * 4U + 1U]);
    auto const y = static_cast<std::uint8_t>(buffer[0][i * 4U + 2U]);
    auto const k = static_cast<std::uint8_t>(buffer[0][i * 4U + 3U]);

    // 归一化到0-1浮点范围
    float c_norm = c / 255.0f;
    float m_norm = m / 255.0f;
    float y_norm = y / 255.0f;
    float k_norm = k / 255.0f;

    // 标准CMYK转线性RGB公式
    float r_linear = (1.0f - c_norm) * (1.0f - k_norm);
    float g_linear = (1.0f - m_norm) * (1.0f - k_norm);
    float b_linear = (1.0f - y_norm) * (1.0f - k_norm);

    // sRGB Gamma校正（gamma=2.2）
    const float gamma = 2.2f;
    float r_gamma = pow(r_linear, 1.0f / gamma);
    float g_gamma = pow(g_linear, 1.0f / gamma);
    float b_gamma = pow(b_linear, 1.0f / gamma);

    // 转换回uint8_t并限制范围在0-255
    auto r = static_cast<std::uint8_t>(std::clamp(r_gamma * 255.0f, 0.0f, 255.0f));
    auto g = static_cast<std::uint8_t>(std::clamp(g_gamma * 255.0f, 0.0f, 255.0f));
    auto b = static_cast<std::uint8_t>(std::clamp(b_gamma * 255.0f, 0.0f, 255.0f));

    // 写入BGRA缓冲区
    pOutput[i * 4U] = b;
    pOutput[i * 4U + 1U] = g;
    pOutput[i * 4U + 2U] = r;
    pOutput[i * 4U + 3U] = 255;
}

3. 可选：整数运算精度优化 #

若不想使用浮点运算，可通过整数乘法加偏移量实现四舍五入，避免直接除法的截断误差：

// 整数运算版本，+254实现四舍五入
auto const r = static_cast<std::uint8_t>(((255 - c) * (255 - k) + 254) / 255);
auto const g = static_cast<std::uint8_t>(((255 - m) * (255 - k) + 254) / 255);
auto const b = static_cast<std::uint8_t>(((255 - y) * (255 - k) + 254) / 255);

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