The treatment of hypercubic lattice artifacts is essential for the calculation of non-perturbative renormalization constants of RI-MOM schemes. It has been shown that for the RI-MOM scheme a large part of these artifacts can be calculated and subtracted with the help of diagrammatic Lattice Perturbation Theory (LPT). Such calculations are typically restricted to 1-loop order, but one may overcome this limitation and calculate hypercubic corrections for any operator and action beyond the 1-loop order using Numerical Stochastic Perturbation Theory (NSPT). In this study, we explore the practicability of such an approach and consider, as a first test, the case of Wilson fermion bilinear operators in a quenched theory. Our results allow us to compare boosted and unboosted perturbative corrections up to the 3-loop order.