I think there is scope to add higher level functions, i.e. to take in ndcube / map objects and return velocity cubes.

One function to take two 2D nddata objects and specify a δt or a function to take a 3D (space, space, time) nddata and output a (space, space, time, velocity) nddata