Гауссовский фильтр C++

Написал программу, которая применяет фильтр Гаусса на зашумленное одноканальное изображение. Проблема в том, что на выходе оно получается несколько темнее, чем ожидалось. Не подскажете, в чем может быть проблема?

double** getGaussian(int height, int width, double sigma) {
    double** filter = new double * [width];
    for (size_t i = 0; i < width; ++i) {
        filter[i] = new double[height];
    }
    double sum = 0;
    for (int i = 0; i < height; i++) {
        for (int j = 0; j < width; j++) {
            filter[i][j] = exp(-(i * i + j * j) / (2 * sigma * sigma)) / (2 * 3.14 * sigma * sigma);
            sum += filter[i][j];
        }
    }
    for (int i = 0; i < height; i++) {
        for (int j = 0; j < width; j++) {
            filter[i][j] /= sum;
        }
    }
    return filter;
}
Image GaussianBlur(Image img, double** filter, int c, int r) {
    Image newimg;
    newimg.cols = img.cols - c + 1;;
    newimg.rows = img.rows - r + 1;;
    newimg.data = new unsigned char * [newimg.rows];
    for (size_t i = 0; i < newimg.rows; ++i) {
        newimg.data[i] = new uchar[newimg.cols];
    }
    for (size_t y = 0; y < newimg.rows; y++) {
        for (size_t x = 0; x < newimg.cols; x++) {
            int sum = 0;
            for (int h = y; h < y + r; h++) {
                for (int w = x; w < x + c; w++) {
                    sum += filter[h - y][w - x] * img.data[h][w];
                }
            }
            newimg.data[y][x] = sum;
        }
    }
    for (size_t y = 0; y < c; ++y) {
        delete[] filter[y];
    }
    delete[] filter;
    return newimg;
}
double FindSigma(const Image& img) {
    double mean = 0;
    double sum = 0;
    double count = 0;
    for (size_t y = 0; y < img.rows; y++) {
        for (size_t x = 0; x < img.cols; x++) {
            count += 1;
            sum += img.data[y][x];
        }
    }
    mean = sum / count;
    double sigma = 0;
    for (size_t y = 0; y < img.rows; y++) {
        for (size_t x = 0; x < img.cols; x++) {
            sigma += (img.data[y][x] - mean) * (img.data[y][x] - mean);
        }
    }
    sigma /= count;
    return sigma;
}

Функция getGaussian создает массив заданного размера (height, width), в котором хранятся значения фильтра Гаусса. Формула для расчета значений такая -g(x, y) = 1 / (2 * pi * sigma ^ 2) * exp(-(x^2 + y ^2)/(2*sigma^2)) (из Википедии), где сигма - стандартное отклонение, а g(x,y) - значения в фильтре. Тип double выбран для того, чтобы при умножении на unsigned char (тип данных во входном изображении) не получались нули. Image - это такая структура:

struct Image {
    unsigned char** data = nullptr;
    size_t rows = 0;
    size_t cols = 0;
};

Вызов функции происходит так

double sigma = FindSigma(image);
double** filter = getGaussian(5, 5, sigma);
Image gauss = GaussianBlur(image, filter, 5, 5);

Помогите найти ошибку, пожалуйста.