It is well-known that we can construct an orthonormal basis at any point $p\in M$, where (M,g) is a Riemannian manifold using exponential mapping. This coordinate is called normal coordinates.
My question is: Apart from this construction, can we get a coordinate $\{x^i\}$ such that $$ g_{ij}(p)=\delta_{ij} ? $$ For example using some linear transformations. Here we don’t ask the Christoffel symbol vanishes. In other words we want to construct an orthonormal basis at a fixed point.
It’s worthy pointing out that even using Gram-Schmidt Algorithm, we can get a orthonormal basis but we cannot get the corresponding coordinates.