Smoothing is implemented as simple convolution in which the convolution matrix is a matrix of constant, i.e., every pixel is replaced by the average of its nearest nine neighbors