Method 1 The area of a region $R$ in $2$D is given by
$$
A=\iint_R1\,dA.\tag{1}
$$
Noting the symmetry in the $x$ and $y$ axes, the area in the first quadrant can be multiplied by 4. Rearranging the expression for positive y gives
$$
y = \frac{b}{a}\sqrt{a^2-x^2}.
$$
The region can also be reduced to a single integral, so that (1) is equivalent to
$$
A = \frac{4b}{a}\int_0^a\sqrt{a^2-x^2}\,dx.
$$
This can be evaluated using a trigonometric substitution, such as $x=a\sin t$. Calculating the corresponding limits, $x=0\Rightarrow t=0$ and $x=a\Rightarrow t=\pi/2$, and noting $dx=a\cos t\,dt$ gives the expression
$$
A = \frac{4b}{a}\int_0^{\pi/2}a^2\cos t\cdot\cos t\,dt=2ab\int_0^{\pi/2}\cos(2t)+1\,dt=2ab\cdot\frac{\pi}{2}=\pi ab.
$$
Method 2 The area can also be calculated using a double integral, but this is much more difficult to evaluate. Parameterising the ellipse as $x(t)=a\cos t$ and $y(t)=b\sin t$, (1) can be written as
$$
A=\int_0^{2\pi}\int_0^{R(t)}r\,dr\,dt,
$$
where $R(t)$ is the boundary of the ellipse, given by $R(t)=\sqrt{x^2+y^2}=\sqrt{a^2\cos^2t+b^2\sin^2t} = \frac{1}{2}\pi(a^2+b^2)$. I'll leave the details here for you to confirm that the result is the same as the other techniques.