1. A method for spatial reconstruction of a plenoptic image of an object, the method comprising: accessing a plenoptic image captured by a plenoptic imaging system, the plenoptic imaging system performing a plenoptic imaging process that is modeled by a pupil image function (PIF) that relates an object to a plenoptic image formed by that object; and applying a PIF inversion process to the captured plenoptic image to produce an estimate of the object; wherein: the object comprises K object components; the PIF comprises K component PIFs, each component PIF modeling the plenoptic imaging process for the corresponding object component; and the step of applying a PIF inversion process to the captured plenoptic image produces estimates of each object component; the plenoptic imaging process is modeled as a vector-matrix operation [Iv f ]=[[PIFv 1 ][PIFv 2 ]. .. [PIFv K ]][[Iv o1 ] T [Iv o2 ] T. .. [Iv oK ] T ] T , where the [Iv ok ] are vectors representing the object components, the [PIFv k ] are matrices representing the component PIFs and [Iv f ] is a vector representing the captured plenoptic image; and the step of applying a PIF inversion process comprises solving an optimization problem to minimize a metric with respect to the [Iv ok ] where the metric is a function of ∥[[PIFv 1 ][PIFv 2 ]. .. [PIFv K ]][[Iv o1 ] T [IV o2 ] T. .. [Iv oK ] T ] T −[Iv f ]∥ and ∥.∥ defines a norm or combination of norms.