Added function to interpolate direction vector fields (e.g. fibre directions): interpDirVec.m