This is a follow up question to (Distribution of distance from origin to any point of Poisson point process)
Assume a homogeneous Poisson point process in a plane (2D) with density $\lambda$. Let $n$, the number of points, be random according to the homogeneous Poisson point process.
Let $\{r_1, r_2, \ldots, r_n\}$ be the set of radial distances of the points from the origin ordered in increasing distance. Let $\{(x_1,y_1), (x_2,y_2) \ldots, (x_n,y_n)\}$ be the coordinates of the points ordered by the same increasing radial distance.
In (Distribution of distance from origin to any point of Poisson point process) the distribution of $r_i$, $i \in \{1,2,...,n\}$ is found. I am interested in the distribution of $\theta_i=\arctan\frac{y_i}{x_i}$.
Intuitively I believe it should be uniform over $(-\pi/2,\pi/2)$ as I can subdivide the expanding disc used to find contact distances into equal spaced non-overlapping segments, and by the properties of a PPP the probability of points being inside each segment must be equal.
The uniform distribution also is confirmed via numerical simulation. I am hoping there is a formal way to show that it is a uniform random variable.