I'm trying to find the area inside the intersection of the region bounded by the ellipses $\frac{x^2}{a^2}+\frac{y^2}{b^2}=1$ and $\frac{x^2}{b^2}+\frac{y^2}{a^2}=1$ using vector calculus methods.
I understand how to calculate area of regions by taking the line integral around the boundary and using Green's Theorem. However what I'm stuck with is finding a paramterization of the boundary! Can anyone point me in the right direction?

I suppose $a<b$. The intersection points of the two ellipses are (see the figure where $a=2$ and $b=4$) $A,B,C,D =(\pm k,\pm k)$ with $$k=\frac{ab}{\sqrt{a^2+b^2}}$$
and, by symmetry, the searched area is $$ 4\left(k^2+\frac{2b}{a}\int_k^a \sqrt{a^2-x^2}dx\right) $$