In the monograph Kounchev, O. I., Multivariate Polysplines. Applications to Numerical and Wavelet Analysis, Academic Press, San Diego-London, 2001, and in the paper Kounchev O., Render, H., Cardinal interpolation with polysplines on annuli, Journal of Approximation Theory 137 (2005) 89--107, we have introduced and studied a new paradigm for cardinal interpolation which is related to the theory of multivariate polysplines. In the present paper we show that this is related to a new sampling paradigm in the multivariate case, whereas we obtain a Shannon type function $S(x) $ and the following Shannon type formula: $f(r��) =\sum_{j=-\infty}^{\infty}\int_{\QTR{Bbb}{S}^{n-1}}S(e^{-j}r��) f(e^{j}��) d��.$ This formula relies upon infinitely many Shannon type formulas for the exponential splines arising from the radial part of the polyharmonic operator $��^{p}$ for fixed $p\geq 1$. Acknowledgement. The first and the second author have been partially supported by the Institutes partnership project with the Alexander von Humboldt Foundation. The first has been partially sponsored by the Greek-Bulgarian bilateral project BGr-17, and the second author by Grant MTM2006-13000-C03-03 of the D.G.I. of Spain.

Submitted to the conference proceedings of SAMPTA07 held in Thessaloniki, Greece, 2007