(possibly later)
scipy.spatial has similar functions for distance matrices, some in cython and faster.

general comment: multidimensional `time`: if we mix real time and one or two spatial indices, then we have an asymmetry in the distance metric across indices. Autocorrelation across time might/will be different then spatial autocorrelation.