We derive the high frequency limit of the Helmholtz equations in terms of quadratic observables. We prove that it can be written as a stationary Liouville equation with source terms. Our method is based on the Wigner Transform, which is a classical tool for evolution dispersive equations. We extend its use to the stationary case after an appropriate scaling of the Helmholtz equation. Several specific difficulties arise here; first, the identification of the source term (which does not share the quadratic aspect) in the limit, then, the lack of L^{2} bounds which can be handled with homogeneous Morrey-Campanato estimates, and finally the problem of uniqueness which, at several stage of the proof, is related to outgoing conditions at infinity.