By using the polarization identity, we propose a family of quasi-interpolants based on bivariate cubic super splines defined on triangulations with a Powell-Sabin refinement. Their spline coefficients only depend on a set of local function values. The quasi-interpolants reproduce cubic polynomials and have an optimal approximation order.