added functions to compute spatial center-surround filter and working on simulation