Ultrasound imaging at high spatial resolution that makes use of both magnitude and phase of echoes in an image reconstruction process that applies unique constraints to a fitting of echoes from the object of interest to echoes from an array of known scatterers.