Let
$$\Delta(\theta) = a^2 \cos^2(\theta) + b^2\sin^2(\theta)$$
In terms of the parametrization $(x,y) = (ar\cos\theta,br\sin\theta)$,
the integral $\mathscr{I}$ we want can be rewritten as
$$\begin{align}
\mathscr{I} = & ab \int_0^{\frac{\pi}{2}} d\theta \int_0^1 \frac{rdr}{\sqrt{1-r^2}\sqrt{\Delta(\theta) r^2 + 1 - r^2}^3}\\
= & \frac{ab}{2}\int_0^{\frac{\pi}{2}} d\theta \int_0^1 \frac{dr^2}{(1-r^2)^2}\frac{1}{\sqrt{\Delta(\theta) \frac{r^2}{1-r^2} + 1}^3}
\end{align}$$
Let $u = \frac{r^2}{1-r^2}$, $v = \Delta(\theta) u + 1$ and $t = \tan\theta$.
Notice $\displaystyle du = \frac{dr^2}{(1-r^2)^2} $, we have
$$\mathscr{I} = \frac{ab}{2}\int_0^{\frac{\pi}{2}} d\theta \int_0^{\infty} \frac{du}{\sqrt{\Delta(\theta) u + 1}^3}
= \frac{ab}{2}\int_0^{\frac{\pi}{2}} \frac{d\theta}{\Delta(\theta)} \int_1^{\infty} \frac{dv}{\sqrt{v}^3}\\
= ab \int_0^{\frac{\pi}{2}} \frac{d\theta}{a^2\cos^2(\theta) + b^2\sin^2(\theta)}
= ab \int_0^{\infty} \frac{dt}{a^2 + b^2 t^2}
= \frac{\pi}{2}
$$