# planar Poisson line process & angles of inclination

Let $L_j$ denote the lines of a planar Poisson line process. By Example 20 of

inclination angles $\omega_{j}$ of the lines $L_j$ relative to the $x$-axis are independent and identically distributed with density $\sin(\omega)/2$ on $[0,\pi]$. In particular, the vertical lines are weighted more than horizontal lines.

The R package spatstat has planar random process simulation capabilities. I can generate Poisson lines in a window via rpoisline and determine their inclination angles $\omega_{j}$ via angles.psp (with option directed=FALSE). To avoid edge effects, I restrict attention to only those lines fairly close to the window center. Surprisingly, the angles seem to follow not the sine density, but instead the uniform density on $[0,\pi]$.

What am I missing please? I would be grateful for help, e.g., R, Matlab or Mathematica code giving experimental results that match theoretical prediction.

Thank you very much.

The following code shows the paradox mentioned by @adrian-baddeley experimentally. In a circular sampling window we see a nearly uniform distribution of angles:

library(spatstat)
X  <- rpoisline(400, win = disc())
aa <- angles.psp(X)
plot(density(aa))


In a long thin sampling window we see a nearly sinusodial distribution of angles:

X  <- rpoisline(1000, win = owin(c(0,10), c(0,.01)))
aa <- angles.psp(X)
plot(density(aa))