What is the notion of a spatial angle in general relativity? Is there a notion of spatial angles in general relativity?
Example:
The world line of a photon is given by $x^{\mu}(\lambda)$. Suppose it flies into my lab where I have a mirror. I align the mirror in such a way, that I measure a right angle between the incoming and outgoing photon. How can I now calculate the wordlline of the outgoing photon?
Of course I cheated a bit since I just used the old euclidean notion of an angle when I say "I measure a right angle between the incoming and outgoing photon". But I hope that this is allowed locally.
 A: If it is a photon, then you know that $x^{a}x_{a}$ is always zero.  You have measured the spatial direction of the photon with your mirror apparatus, so you know the values of the $x^{i}(0)$ at some time${}^{1}$, which we will call zero.  Furthermore, we can infer the value of $x^{0}$ at this time from the fact that $x^{a}x_{a} = 0$.
Then, all we need is to know the metric tensor, and we can simply solve the set of ODEs given by,
$$\frac{\mathrm{d}^2 x^a (\lambda)}{\mathrm{d}\lambda^2} + \Gamma^a_{bc} \, \frac{\mathrm{d} x^b(\lambda)}{\mathrm{d}\lambda} \frac{\mathrm{d}x^c (\lambda)}{\mathrm{d}\lambda} = 0$$
for all subsequent times${}^{2}$.  Typically, this will have to be done numerically, but there are some known exact solutions for simplified cases.  $\lambda$ is the null affine parameter, and the $\Gamma$'s are the Christoffel symbols, which are defined from the metric tensor.  For the familiar example of a spherically symmeteric gravitating body AND for weak gravity in nonrelativistic motion,  the only relevant component of the metric tensor is $g_{tt} = -\left(1-\frac{2M}{r}\right)$, and the only relevant Christoffel symbol is $\Gamma_{tt}{}^{r} = \frac{M}{r^{2}}$, and our theory becomes Newtonian gravity.
 
 
 
${}^{1}$Here, I take letters from the beginning of the alphabet to be spacetime indices, while letters from the middle of the alphabet are sums over only the spatial components.
${}^{2}$ note that "time" (given by $\lambda$) here doesn't really mean time in the proper sense.  Null geodesics don't experience any clock motion.  Null geodesics, however, DO trace out "null affine parameter", which you can interpret as how many perpendicular null geodesics they cross as the photon moves.
A: As a geometric quantity, the value of an "angle" can be determined and expressed in a coordinate-free way:
Given three pairwise space-like events, "$A$", "$B$" and "$C$", and given the positive real numbers
$\frac{s^2[ A C ]}{s^2[ A B ]}$, $\frac{s^2[ A C ]}{s^2[ B C ]}$ and $\frac{s^2[ B C ]}{s^2[ A B ]} = \frac{s^2[ A C ]}{s^2[ A B ]} / \frac{s^2[ A C ]}{s^2[ B C ]}$
as ratios between squares of interval magnitudes (or minimum arc lengths) between event pairs,
then the value of "angle at $B$, between $A$ and $C$" may either be expressed (directly) as     
$\angle [ A B C ] := $
$\text{ArcSin} \left[ \frac{1}{2} \sqrt{ 2 + 2 \frac{s^2[ A C ]}{s^2[ A B ]} + 2 \frac{s^2[ A C ]}{s^2[ B C ]} - \frac{s^2[ B C ]}{s^2[ A B ]} - \frac{s^2[ A B ]}{s^2[ B C ]} - \frac{s^2[ A C ]}{s^2[ A B ]} \frac{s^2[ A C ]}{s^2[ B C ]}  } \right];$    
or, more generally, in terms of squares of interval magnitudes (or minimum arc lengths) involving additional events "$F$", "$G$" (named in the sense of variables) which are space-like to each other as well as to "$A$", "$B$" and "$C$", as     
$\angle [ A B C ] := \text{Limit}_{\{ F, G \}} {\huge[} $
${\Large \{ } \frac{s^2[ B F ]}{s^2[ A B ]} \rightarrow 0, \frac{s^2[ B G ]}{s^2[ B C ]} \rightarrow 0, $
$2 + 2 \frac{s^2[ B F ]}{s^2[ A B ]} + 2 \frac{s^2[ B F ]}{s^2[ A F ]} - \frac{s^2[ A F ]}{s^2[ A B ]} - \frac{s^2[ A B ]}{s^2[ A F ]} - \frac{s^2[ B F ]}{s^2[ A B ]} \frac{s^2[ B F ]}{s^2[ A F ]} \rightarrow 0, $
$2 + 2 \frac{s^2[ B G ]}{s^2[ B C ]} + 2 \frac{s^2[ B G ]}{s^2[ C G ]} - \frac{s^2[ C G ]}{s^2[ B C ]} - \frac{s^2[ B C ]}{s^2[ C G ]} - \frac{s^2[ B G ]}{s^2[ B C ]} \frac{s^2[ B G ]}{s^2[ C G ]} \rightarrow 0 {\Large \} };  $
$\text{ArcSin} \left[ \frac{1}{2} \sqrt{ 2 + 2 \frac{s^2[ F G ]}{s^2[ B F ]} + 2 \frac{s^2[ F G ]}{s^2[ B G ]} - \frac{s^2[ B G ]}{s^2[ B F ]} - \frac{s^2[ B F ]}{s^2[ B G ]} - \frac{s^2[ F G ]}{s^2[ B F ]} \frac{s^2[ F G ]}{s^2[ B G ]}  } \right] {\huge]} . $ 
If the region containing the events $A$, $B$ and $C$ is flat then these two "angle" values are equal.
