We obtain a model independent expression for the muonic hydrogen Lamb shift up to O(m µ α 6 , m µ α 5 m 2 µ m 2 ρ ). The hadronic effects are controlled by the chiral theory, which allows for their model independent determination. We give their complete expression including the pion and Delta particles. Out of this analysis and the experimental measurement of the muonic hydrogen Lamb shift we determine the electromagnetic proton radius: r p =0.8412(15) fm. This number is at 6.8σ variance with respect to the CODATA value. The parametric control of the uncertainties allows us to obtain a model independent determination of the error, which is dominated by hadronic effects.