img50_nd = scipy.ndimage.interpolation.zoom(img, (0.5, 0.5, 1))
plot(img50_order0[(50), :, (1)])
plot(img50_order1[(50), :, (1)])
plot(img50_lm[(50), :, (1)])
plot(img50_nd[(50), :, (1)])
Image gradients become inaccurate when downscaling using a variety of different methods