use savgol2d() for gradient
use savgol2d() to compute derivatives
use savgol2d() to compute gradient