Derivative of 2D dose (intensity) image wrt distance

I have a 2D image of size M×N representing a radiation dose map with a local dose , where is the spatial coordinate. I want to estimate both derivatives and with respect to r - how do I do that? If neccessary, the spatial resolution is 21 /pixel and isotropic.

Thanks in advance for any help or guidance!