pub fn centered(img: &DMatrix<u8>) -> (DMatrix<i16>, DMatrix<i16>)
Compute a centered gradient.
1/2 * ( img(i+1,j) - img(i-1,j), img(i,j+1) - img(i,j-1) )
Gradients of pixels at the border of the image are set to 0.