1

我正在用 C 语言制作一个图像大小调整程序，目前我在使用双线性插值函数时遇到了问题（这是我正在使用的众多函数之一）。这个问题只出现在 16 位位图上，如果我使用 24 位版本，它会完美地调整它们的大小。

这是我的双线性插值代码。n_w 和 n_h 是图像的新宽度和高度：

#define getelm(x)   (((pix+index)->x)*(1-xdiff)*(1-ydiff))+(((pix+index+1)->x)*(xdiff)*(1-ydiff))+(((pix+index+o_w)->x)*(1-xdiff)*(ydiff))+((pix+index+o_w+1)->x)*xdiff*ydiff

int pad = (2*n_w) & 3;
if (pad)
pad = 4-pad;
uint16_t *buffer;
if (buf)
buffer = malloc(2*n_w);
for (i = 0; i < n_h; i++) {
    for (j = 0; j < n_w; j++) {
        x = (int)(j*xrat);
        y = (int)(i*yrat);
        xdiff = (xrat*j)-x;
        ydiff = (yrat*i)-y;
        index = y*o_w+x;
        uint16_t container = 0;
        container |= (int)(round(getelm(b))) << 11;
        container |= (int)(round(getelm(g))) << 5;
        container |= (int)round(getelm(r));
        if (buf)
            *(buffer+j) = container;
        else
            fwrite(&container, 1, 2, dest);
    }
    if (buf)
        fwrite(buffer, 1, 2*n_w, dest);
        fwrite(&pad, 1, pad, dest);
}

我的这段代码的 24 位版本（唯一的区别是不使用容器，而是 3 个 8 位整数保存 RGB 值）运行良好。

然而，这段代码给出了奇怪的结果。看下图：原始图像

当我调整它的大小时，它给了我这个：调整大小的图像

我不明白为什么会发生这种情况，特别是当它适用于 24 位位图时，而且其他一些调整大小的算法（例如最近邻）以与这应该相同的方式使用 16 位。

编辑：我不认为这是一个溢出问题，因为添加以下代码在运行时没有输出：

if (MAX((int)(getelm(b)), 31) > 31)
    printf("blue overflow: %.10f\n", (getelm(b)));
if (MAX((int)(getelm(g)), 63) > 63)
    printf("green overflow: %.10f\n", (getelm(g)));
if (MAX((int)(getelm(r)), 31) > 31)
    printf("red overflow: %.10f\n", (getelm(r)));

编辑2：我也不认为这是一个下溢问题，这什么都不做：

if ((getelm(b)) < 0 || (getelm(g)) < 0 || (getelm(r)) < 0)
    printf("Underflow\n");
4

1 回答 1

2

假设数据pix具有类型

struct
{
    uint16_t r : 5;
    uint16_t g : 6;
    uint16_t b : 5;
};

的计算中存在错误container。使用round不会总是防止溢出。下一个代码将：

uint16_t container = 0;
        container |= ((int)(round(getelm(b))) & 31) << 11; // optionally
        container |= ((int)(round(getelm(g))) & 63) << 5;
        container |= ((int)round(getelm(r)) & 31);

或保留最大的丢失信息：

uint16_t container = 0;
        container |= min((int)(round(getelm(b))) , 31) << 11;
        container |= min(((int)(round(getelm(g))) , 63) << 5;
        container |= min(((int)round(getelm(r)) , 31);

编辑

由于pix->r和来自 8 位值，pix->g因此pix->b同样的推理适用于它们并且需要检查它们的范围。

由于白色区域变为紫色，这意味着绿色由于溢出而被抑制，或者首先被读取为零。在这种情况下，检查读取颜色会有所帮助。类似地，黑色变成绿色意味着代表一个小值的位被移动并且颜色以某种方式反转。

为了找到错误，我建议将代码拆分为小函数并声明每个函数的输入。

于 2013-08-17T00:42:06.010 回答