Quasi-interpolation operators based on the trivariate seven-direction $C^2$ quartic box spline