This paper proposes an algorithm to solve the phase retrieval problem in optical imaging in two stages. First, we introduce an initialization to estimate the support of the sparse representation θ in some sparse basis Ψ of an image x. Second, we solve a smoothing optimization problem to reconstruct x. Simulations show that the proposed algorithm requires less number of measurements compared with existing methods.