The inverse problem we consider in this paper is to determine the shape of an obstacle from the knowledge of the far field pattern for scattering of time-harmonic plane waves. In the case of scattering from a sound-soft obstacle, we will interpret Huygens' principle as a system of two integral equations, named data and field equation, for the unknown boundary of the scatterer and the induced surface flux, i.e., the unknown normal derivative of the total field on the boundary. Reflecting the ill-posedness of the inverse obstacle scattering problem these integral equations are ill-posed. They are linear with respect to the unknown flux and nonlinear with respect to the unknown boundary and offer, in principle, three immediate possibilities for their iterative solution via linearization and regularization. In addition to presenting new results on injectivity and dense range for the linearized operators, the main purpose of this paper is to establish and illuminate relations between these three solution methods based on Huygens' principle in inverse obstacle scattering. Furthermore, we will exhibit connections and differences to the Communicated by guest editors: