Compressive direct imaging of a billion-dimensional optical phase space

Optical phase-spaces represent fields of any spatial coherence, and are typically measured through phase-retrieval methods involving a computational inversion, interference, or a resolution-limiting lenslet array. Recently, a weak-values technique demonstrated that a beamu0027s Dirac phase-space is proportional to the measurable complex weak-value, regardless of coherence. These direct measurements require scanning through all possible position-polarization couplings, limiting their dimensionality to less than 100,000. We circumvent these limitations using compressive sensing, a numerical protocol that allows us to undersample, yet efficiently measure high-dimensional phase-spaces. We also propose an improved technique that allows us to directly measure phase-spaces with high spatial resolution and scalable frequency resolution. With this method, we are able to easily measure a 1.07-billion-dimensional phase-space. The distributions are numerically propagated to an object placed in the beam path, with excellent agreement. This protocol has broad implications in signal processing and imaging, including recovery of Fourier amplitudes in any dimension with linear algorithmic solutions and ultra-high dimensional phase-space imaging.