为什么在应用 Sobel 边缘检测器时我的图像会被裁剪成两半？

如何解决为什么在应用 Sobel 边缘检测器时我的图像会被裁剪成两半？

首先让我说我是 C++ 的新手，所以我真的不知道最佳实践或很好地处理语法。 我正在尝试读取黑白图像，并使用 sobel 算法检测其边缘并输出结果，但在执行过程中出现错误：

munmap_chunk(): invalid pointer
Aborted (core dumped)

虽然输出了图像，但只有一半，我似乎无法弄清楚是什么原因造成的。

我写了以下代码：

#include <iostream>
#include <omp.h>
#include "CImg.h"
using namespace std;
using namespace cimg_library;

int main() {
    const int x_mask [9] = {
        -1,1,-2,2,-1,1
    };

    const int y_mask [9] = {
        -1,1
    };

    const char* fileName = "test.png";
    CImg<float> img = CImg<float>(fileName);
    
    int cols = img.width();
    int lines = img.height();
    CImg<float> output = CImg<float>(cols,lines,0.0);
    printf("Loading %d x %d image...\n",cols,lines);
    const int mask_size = 3;

    int gradient_x;
    int gradient_y;
    // Loop through image ignoring borders
    for(int i = 1; i<cols-1; i++) {
        for(int j = 1; j<lines-1; j++){
            output(j,i) = 0;
            gradient_x = 0;
            gradient_y = 0;

            // Find the x_gradient and y_gradient
            for(int m = 0; m < mask_size; m++) {
                for(int n = 0; n < mask_size; n++) {
                    // Neighbourgh pixels
                    int np_x = j + (m - 1);
                    int np_y = i + (n - 1);

                    float v = img(np_x,np_y);

                    int mask_index = (m*3) + n;

                    gradient_x = gradient_x + (x_mask[mask_index] * v);
                    gradient_y = gradient_y + (y_mask[mask_index] * v);
                }
            }
            float gradient_sum = sqrt((gradient_x * gradient_x) + (gradient_y * gradient_y));
            if(gradient_sum >= 255) {
                gradient_sum = 255;
            } else if(gradient_sum <= 0) {
                gradient_sum = 0;
            }
            output(j,i) = gradient_sum;
        }
    }
    printf("Outputed image of size %d x %d\n",output.width(),output.height());
    output.save("test_edges.png");
 return 0;
}

应用于此图像：

我得到这个输出：