Edge finding in magnetic resonance imaging applications: the calculation of the first order derivative of two dimensional images