Function DifferenceOfGaussian

examples/image_processing/filters.cpp:121–131

}

array DifferenceOfGaussian(const array &in, int window_radius1,
                           int window_radius2) {
    array ret_val;
    int w1   = 2 * window_radius1 + 1;
    int w2   = 2 * window_radius2 + 1;
    array g1 = gaussianKernel(w1, w1);
    array g2 = gaussianKernel(w2, w2);
    ret_val  = (convolve(in, g1) - convolve(in, g2));
    normalizeImage(ret_val);
    return ret_val;
}

array medianfilter(const array &in, int window_width, int window_height) {
    array ret_val(in.dims());