How do I numerically calculate a function from its noisy gradient?
You can filter either the gradient itself or the result, $s$.  You would need to know the characteristics of the true gradients fairly well to know what the frequency bandwidth is.  At that point you could design a low-pass filter that would preserve the signal but get rid of higher frequency noise.