We propose a new algorithm for the estimation of functional activations in diffuse optical imaging. Our approach considers the activations to be support limited. We simultaneously estimate the function values as well as the support from the available measurements. Since this scheme exploit the structure inherent to functional imaging, it provide reconstructions with better spatial resolution and is more robust to noise.