Reformulemos la pregunta de la OP de la siguiente manera:
Demuestre que una transformación local de coordenadas $x^{\mu} \to y^{\rho}=y^{\rho}(x)$ entre dos sistemas de coordenadas locales (en una variedad lorentziana de 3+1 dimensiones) debe ser afín si la métrica $g_{\mu\nu}$ en ambos sistemas de coordenadas resultan estar en la forma plana constante de Minkowski $\eta_{\mu\nu}$ .
Aquí presentaremos una prueba que funciona tanto con la firma de Minkowski como con la de Euclides; de hecho, para cualquier firma y para cualquier número finito no nulo de dimensiones, siempre que la métrica $g_{\mu\nu}$ es invertible.
1) Recordemos primero la propiedad de transformación de la métrica inversa $g^{\mu\nu}$ que es una contravariante $(2,0)$ tensor simétrico,
$$ \frac{\partial y^{\rho}}{\partial x^{\mu}} g^{\mu\nu}_{(x)}\frac{\partial y^{\sigma}}{\partial x^{\nu}}~=~g^{\rho\sigma}_{(y)}, $$
donde $x^{\mu} \to y^{\rho}=y^{\rho}(x)$ es una transformación de coordenadas local. Recordemos que la métrica $g_{\mu\nu}=\eta_{\mu\nu}$ es la métrica constante plana en ambos sistemas de coordenadas. Así que podemos escribir
$$ \frac{\partial y^{\rho}}{\partial x^{\mu}} \eta^{\mu\nu}\frac{\partial y^{\sigma}}{\partial x^{\nu}}~=~\eta^{\rho\sigma}. \qquad (1) $$
2) Supongamos que la transformación local de coordenadas es analítica real
$$y^{\rho} ~=~ a^{(0)\rho} + a^{(1)\rho}_{\mu} x^{\mu} + \frac{1}{2} a^{(2)\rho}_{\mu\nu}x^{\mu}x^{\nu} + \frac{1}{3!} a^{(3)\rho}_{\mu\nu\lambda}x^{\mu} x^{\nu} x^{\lambda} + \ldots. $$
Realizando eventualmente una traslación adecuada, supondremos a partir de ahora, sin pérdida de generalidad, que el desplazamiento constante $ a^{(0)\rho} =0 $ es cero.
3) A la orden zerótica en $x$ la ecuación $(1)$ lee
$$ a^{(1)\rho}_{\mu} \eta^{\mu\nu}a^{(1)\sigma}_{\nu}~=~\eta^{\rho\sigma}, $$
que no es sorprendente que la matriz $a^{(1)\rho}_{\mu}$ es una matriz de Lorentz (o una ortogonal), respectivamente. Realizando eventualmente una "rotación" adecuada, supondremos a partir de ahora, sin pérdida de generalidad, que la matriz constante
$$ a^{(1)\rho}_{\mu}~=~\delta^{\rho}_{\mu} $$
es la matriz unitaria.
4) En lo sucesivo, será conveniente bajar el índice del $y^{\sigma}$ coordinar como
$$y_{\rho}~:=~\eta_{\rho\sigma}y^{\sigma}.$$
Entonces la transformación local de coordenadas se convierte en
$$y_{\rho} ~=~ \eta_{\rho\mu} x^{\mu} + \frac{1}{2} a^{(2)}_{\rho,\mu\nu}x^{\mu}x^{\nu} + \frac{1}{3!} a^{(3)}_{\rho,\mu\nu\lambda}x^{\mu} x^{\nu} x^{\lambda}+ \ldots$$ $$+\frac{1}{n!} a^{(n)}_{\rho,\mu_1\ldots\mu_n}x^{\mu_1} \cdots x^{\mu_n}+ \ldots. $$
5) Al primer orden en $x$ la ecuación $(1)$ lee
$$ a^{(2)}_{\rho,\sigma\mu}+a^{(2)}_{\sigma,\rho\mu}~=~0.$$
Eso es, $a^{(2)}_{\rho,\mu\nu}$ es simétrico en $\mu\leftrightarrow \nu$ pero antisimétrico en $\rho\leftrightarrow \mu$ . No es difícil ver (aplicando la propiedad de simetría y antisimetría en orden alterno tres veces cada una), que los coeficientes de segundo orden $a^{(2)}_{\rho,\mu\nu}=0$ debe desaparecer.
6) Al segundo orden en $x$ la ecuación $(1)$ lee
$$ a^{(3)}_{\rho,\sigma\mu\nu}+a^{(3)}_{\sigma,\rho\mu\nu}~=~0.$$
Eso es, $a^{(3)}_{\rho,\mu\nu\lambda}$ es simétrica en $\mu\leftrightarrow \nu\leftrightarrow \lambda $ pero antisimétrico en $\rho\leftrightarrow \mu$ . Para las instalaciones fijas $\lambda$ podemos llegar de nuevo a la conclusión $a^{(3)}_{\rho,\mu\nu\lambda}=0$ .
7) Del mismo modo, concluimos inductivamente que los coeficientes de orden superior $a^{(n)}_{\rho,\mu_1\ldots\mu_n}=0$ debe desaparecer también. Así que $y^{\mu}= x^{\mu}$ . Q.E.D.