A signal-subspace method is derived for the localization and imaging of unknown scatterers using intensity-only wave field data (lacking field phase information). The method is an extension of the time-reversal multiple-signal-classification imaging approach to intensity-only data. Of importance, the derived methodology works within exact scattering theory including multiple scattering.