Systems and methods of frequency-domain photoacoustic imaging are provided utilizing an ultrasonic phased array probe and intensity modulated optical excitation with coding to improve signal-to-noise ratio. Embodiments employ frequency-domain photoacoustic imaging methodologies such as the photoacoustic radar, coupled with a multi-element ultrasonic sensor array to deliver spatially-resolved correlation images of photoacoustic sources, which may be employed to image optical heterogeneities within tissue-like scattering media.