One approach for approximating a discrete derivative in $2$D is based

on computing differences of the form:

f$($x$+1,$y$)-$ f$($x$,$y$)$ and f$($x$,$y$+1)-$ f$($x$,$y$)$