Applications of soft materials are customarily linked to complex deformation scenarios and material nonlinearities. In the bioengineering field, soft materials typically mimic the low stiffness of biological matter subjected to extreme deformations. Computational frameworks surge as a versatile tool to assist the design of functional applications. The constitutive model lies at the core of such frameworks. In this regard, the customary extreme non-linear behavior of elastomers poses an additional challenge to thoroughly capture the material behavior. Here, data-driven methodologies hold considerable promise for enhancing constitutive modeling when contrasted with phenomenological approaches. In this investigation, we introduce a versatile data-adaptive method tailored to the modeling of hyperelastic soft materials at finite strains. Specifically, our method substitutes an a priori chosen model for the strain energy function by a flexible interpolant defined on a discretized invariant space. Within this framework, the interpolation values assume the role of material parameters and are determined through finite element model updating to conform to measured experimental data — comprising full-field displacements coming from Digital-Image-Correlation and global reaction forces. We validate the method on uniaxial experimental tests of soft elastomers, encompassing ELASTOSILTM, DOWSILTM, and V HBTM. Overall, we aim to establish a new route for the construction of hyperelastic energy functions, untethered from any predefined existing models or assumptions regarding the shape of the energy.