Since Einstein’s Lagrangian is only weakly diffeomorphism invariant, the situation would be even worse if we used
Einstein’s Lagrangian. The corresponding canonical quantities would still be coordinate dependent, though in certain ‘natural’
coordinate systems they yield reasonable results (see, e.g., [7] and references therein).