A system and method for quantitative functional neuroimaging through thick brain tissue in live animals. A computational imaging method is disclosed that uses plenoptic image acquisition including a first initialization step that identifies individual neurons by their optical signature and provides a reliable estimate of their position in space and a second stimulation-based image processing step that used acquired calibration data to quickly quantify activity in each identified neuron at video frame-rate.