The ergodic theory of the open KPZ equation has seen significant progress in recent years, with explicit invariant measures described in a series of works by Corwin--Knizel, Barraquand--Le Doussal, and Bryc--Kuznetsov--Wang--Wesołowski. In this paper, we provide a stochastic analytic proof of the formula for the invariant measures. Our approach starts from the Gaussian invariant measure for the case of homogeneous boundary conditions. We approximate the inhomogeneous problem by a homogeneous one with a singular boundary potential. Using tools including change of measure, time reversal for Markov processes, and Itô's formula, we then reduce the problem to analyzing the KPZ nonlinearity in a thin boundary layer. Finally, using the theory of regularity structures, we establish a central limit theorem for the time-integrated nonlinearity near the boundary, which completes the proof of the invariance. Although it is known that different boundary parameters give rise to distinct physical regimes for the invariant measures, our method is robust and does not rely on any particular choice of boundary parameters.