We consider diffusion fields induced by a finite number of spatially localized sources and address the problem of estimating these sources using spatiotemporal samples of the field obtained with a sensor network. Within this framework, we consider two different time evolutions: the case where the sources are instantaneous, as well as, the case where the sources decay exponentially in time after activation. We first derive novel exact inversion formulas, for both source distributions, through the use of Green's second theorem and a family of sensing functions to compute generalized field samples. These generalized samples can then be inverted using variations of existing algebraic methods such as Prony's method. Next, we develop a novel and robust reconstruction method for diffusion fields by properly extending these formulas to operate on the spatiotemporal samples of the field. Finally, we present numerical results using both synthetic and real data to verify the algorithms proposed herein.