In the context of dynamical systems, nonlinearity measures quantify the strength of their nonlinearity by means of the distance of their input-output behaviour to a set of linear input-output mappings. In this paper, we establish a framework to determine nonlinearity measures and other optimal input-output properties for nonlinear polynomial systems without explicitly identifying a model but directly from a finite number of input-state measurements which are subject to noise. To this end, we deduce from noisy data for the unidentified ground-truth system three set-membership representations for which we prove asymptotic consistency with respect to the number of samples. Then, we leverage these representations to compute guaranteed upper bounds on nonlinearity measures and the corresponding optimal linear approximation model via semi-definite programming. Furthermore, we apply the framework to determine optimal input-output properties described by certain classes of time domain hard integral quadratic inequalities.