Abstract: A method for estimating a real function that describes a phenomenon occurring in a space of any dimensionality is disclosed. The function is estimated by taking a series of measurements of the phenomenon being described and using those measurements to construct an expansion that has a manageable number of terms. A reduction in the number of terms is achieved by using an approximation that is defined as an expansion on kernel functions, the kernel functions forming an inner product in Hilbert space. By finding the support vectors for the measurements one specifies the expansion functions. The number of terms in an estimation according to the present invention is generally much less than the number of observations of the real world phenomenon that is being estimated. In one embodiment, the function estimation method may be used to reconstruct a radiation density image using Positron Emission Tomography (PET) scan measurements.