Consider the problem of sampling signals that are nonbandlimited but have finite number of degrees of freedom per unit of time and call this number the rate of innovation. Streams of Diracs and piecewise polynomials are the examples of such signals, and thus are known as signals with finite rate of innovation (FRI) [3]. We know that the classical ('bandlimited-sinc') sampling theory does not enable perfect reconstruction of such signals from their samples since they are not bandlimited. However, the recent results on FRI sampling [3], [4] suggest that it is possible to sample and perfectly reconstruct such nonbandlimited signals using a rich class of kernels.In this paper, we extend the results of [4] in higher dimensions using compactly supported kernels that reproduce polynomials (satisfy Strang-Fix conditions). In fact, the polynomial reproduction property of the kernel makes it possible to obtain the continuous-moments of the signal from its samples. Using these moments and the annihilating filter method (Prony's method), the innovative part of the signal, and therefore, the signal itself is perfectly reconstructed. In particular, we present local (directional derivatives based) and global (complex-moments, Radon transform based) sampling schemes for classes of FRI signals such as sets of Diracs, bilevel and planar polygons, quadrature domains (e.g. circles, ellipses, cardioids), 2-D polynomials with polygonal boundaries, and n-dimensional Diracs and convex polytopes. This research has been promisingly explored in super-resolution algorithms [5] and distributed compression [6], and might find its applications in photogrammetry, computer graphics, and machine vision.