In this article, we provide a numerical algorithm to reconstruct a convex sound-soft scatterer from phaseless backscattering data, assuming sufficiently high frequency. Certain uniqueness and existence results for the case of circular scatterers are given as well, based on the asymptotic expansion for the normal derivative of the total field.