In this paper, we report a significant recovery of the linear baryonic acoustic oscillation (BAO) signature by applying the isobaric reconstruction algorithm to the non-linearly evolved matter density field. Assuming that only the longitudinal component of the displacement is cosmologically relevant, this algorithm iteratively solves the non-linear coordinate transform between the Lagrangian and Eulerian frames without requiring any specific knowledge of the dynamics. For dark matter field, it produces the non-linear displacement potential with very high fidelity. The reconstruction error at the pixel level is within a few percent caused only by the emergence of the transverse component after the shell-crossing. As this method circumvents one of the most strongest non-linearity in density field, the reconstructed field is well-described by linear theory and is immune from the bulk-flow smearing of the BAO signature, and therefore could be used to significantly improve the precision of measuring the sound horizon scale. For a perfect large-scale structure survey at redshift zero without Poisson or instrumental noise, the fractional error is reduced by a factor of 2.7, and is very close to the ideal limit one would ever achieve with linear power spectrum and Gaussian covariance matrix.