Bivariate $C^2$ cubic spline quasi-interpolants on uniform Powell-Sabin triangulations of a rectangular domain