In this paper we propose a method for the approximation of high-dimensional functions over finite intervals with respect to complete orthonormal systems polynomials. An important tool is multivariate classical analysis variance (ANOVA) decomposition. For low-dimensional structure, i.e., low superposition dimension, are able achieve reconstruction from scattered data and simultaneously understan...