In this paper on page 7 it says that the difference between the whole length of the hyperbola and that of its asymptotes, $\Delta$, is given by $$\Delta=4a\int_0^{\alpha}\sqrt{1-e^2\sin^2\theta}d\theta,$$ where $\sin\alpha=1/e$, $e$ is the eccentricity, and this is for a hyperbola with equation $x^2/a^2-y^2/b^2=1$. No proof was given and I couldn't find one online.
I have tried to find expressions for the length of the asymptote and the length of the hyperbola in the 1st quadrant using polar coordinates with $L=\int_{\theta_0}^{\theta_1}\sqrt{r^2+(dr/d\theta)^2}d\theta$, but the difference was always an ugly expression instead of what was given(I'm not sure if converges if I do it this way). I also tried taking $\theta$ to be the acute angle between the tangent to the hyperbola and its asymptote(in the 1st quadrant) but I do not know how to interpret the integral geometrically or otherwise. I also tried substituting $u=\sqrt{1-e^2\sin^2\theta}$, but I cannot make sense of the resulting integral.
I have no idea how to proceed now, can someone please help?