In this study, we enhance the numerical approach of Lee et al.$^{(1)}$ to spherical indentation technique for property evaluation of hyper-elastic rubber. We first determine the friction coefficient between rubber and indenter in a practical viewpoint. We perform finite element numerical simulations for deeper indentation depth. An optimal data acquisition spot is selected, which features sufficiently large strain energy density and negligible frictional effect. We then improve two normalized functions mapping an indentation load vs. deflection curve into a strain energy density vs. first invariant curve, the latter of which in turn gives the Yeoh-model constants. The enhanced spherical indentation approach produces the rubber material properties with an average error of less than 3%.