The exact analytical gradient of the perturbed field can be calculated by simply replacing the convolution kernel, , with its gradient and neglecting the first term, since it is proportional to and has zero gradient at the centre of the voxel. This can be also be seen from equation 24 which gives
(26) |