NF: function giving time difference for images

New function that accepts an image as input and returns time series
difference analysis on the image.
WIP: refactor tsdiffana script to select slice axis