PPM 格式卷积的特殊问题

Peculiar problem regarding convolution in PPM format

调试后,尝试了不同的图像软件(xdg,gimp),我仍然有一个错误让我完全失望。

问题是关于 PPM 格式的卷积,对于纵横比不同的图像,我使用的是 1500x1000 像素的图像,其中 {0,0,0, 0,1,0, 0,0,0} 的掩码有效很好(它只是复制图像),但是对于第一行或最后一行不同于 0 的掩码,例如。 {0,1,0, 0,0,0, 0,0,0} 图像向右移动其大小的 1/3。我觉得它很奇怪,因为据我所知,我没有溢出或任何可能导致此问题的指针算法。

我已经将其缩小到卷积核。 Afaik 我在保存、读取图像时没有任何问题,在 运行 i_convolution 之后它只是按预定义的值移动图像?

void i_convolution(unsigned int **in, unsigned int ***out,
                int y_max, int x_max, int kernel_size)
{
    int kernel_sum = 0;
    for(int i = 0; i < kernel_size; i++)
    {
        for(int j = 0; j < kernel_size; j++)
        {
            kernel_sum += kernel[i * kernel_size + j];
        }
    }

    printf("kernel sum = %d\n", kernel_sum);

    for (int i = 1; i < y_max - 1; i++)
    {
        for (int j = 1; j < x_max - 1; j++)
        {
            int r = 0;
            int g = 0;
            int b = 0;

            for (int y_conv = -1; y_conv <= 1; y_conv++)
            {
                for (int x_conv = -1; x_conv <= 1; x_conv++)
                {
                    int y_index = i + y_conv;
                    int x_index = j + x_conv;

                    unsigned char rval = (unsigned char)(in[y_index][x_index] & 0xff);
                    unsigned char gval = (unsigned char)((in[y_index][x_index] & 0xff00) >> 8);
                    unsigned char bval = (unsigned char)((in[y_index][x_index] & 0xff0000) >> 16);

                    int kernel_val = kernel[(y_conv + 1)*kernel_size + (x_conv + 1)];

                    r += (int)(rval * kernel_val);
                    g += (int)(gval * kernel_val);
                    b += (int)(bval * kernel_val);
                }
            }

            r /= kernel_sum;//median filtration
            g /= kernel_sum;//median filtration
            b /= kernel_sum;//median filtration
            // b = abs(b);

            if (r > 255) r = 255;
            else if(r < 0) r = 0;
            if (g > 255) g = 255;
            else if(g < 0) g = 0;
            if (b > 255) b = 255;
            else if(b < 0) b = 0;

            unsigned int val;
            val = 0;

            val |= b & 0xff;
            val <<= 8;

            val |= g & 0xff;
            val <<= 8;

            val |= r & 0xff;
            (*out)[i][j] = val;

        }
    }
}

让我们使用内核 {0, 1, 0, 0, 0, 0,
结果是这样的,左边是原始的,右边是卷积后的 https://i.imgur.com/rzXKjUY.png

如有任何帮助,我将不胜感激。 最好的问候。

我将其标记为已解决,因为我误解了 PPM 格式的高度和宽度,导致了这种行为,将 y 与 x 交换(并这样分配内存)解决了这个问题!