How to deduce the area of sphere in polar coordinates?

$A = r \int_{0}^{\pi}\int_{0}^{2\pi} e^{i (\alpha+\theta)} d\alpha d\theta = r \int_{0}^{\pi} [\frac{-i}{\alpha+\theta} e^{i(\alpha + \theta)}]_{0}^{2\pi} d\theta = ... = r \int_{0}^{\pi} e^{2i\alpha} ( \frac{1}{\alpha} - \frac{1}{\alpha+2\pi} ) d\theta$

and according to wolfram it is wrong, here. So how do you deduce the area of a sphere in polar coordinates?

[Edit] Trevor gave us formula $r^2 \int_{2π}^{0}\int^{π}_{0}sin(\phi) d\phi d\theta$ without deduction, by deduction, I mean things on which it is based on.

1. Why is it odd function $sin(\phi)$ and not instead an even function $cos(\phi)$? (I know that it is and I can calculate it but some proof here missing)
2. Why no term with $\theta$ but with $\phi$?
3. Why $r^{2}$? -Okay it is area and by symmetry it sounds logical but why $r^{2}$ and not for example $rd$ where $d$ is some lenght or even a path?
4. What about an arbitrary case?
5. What determines the order in integration?

Many things open to really deduce it.

[Edit] Zhen Lin has an excellent answer where s/he noted that 2-sphere, i.e. the sphere in three dimensional euclidean space, is two-dimensional. Next deductions should be on 0-sphere and then incrementally on others. More about n-sphere here.

1. How do you deduce area in 0-sphere, 1-sphere, 2-sphere, and n-sphere?
2. Why does the Jacobian matrix in 2-sphere-in-3-dim-euclidean-case have only two two variables $(\phi, \theta)$?
3. Are $\theta$ and $\phi$ non-free-variables because they are bounded somehow by the manifold (maybe misusing terminology)?
4. What is $r$ in the sphere? For example in the 2-sphere?
5. I can sense that the var $r$ is somehow different to $\phi$ and $\theta$ because what is it if you derivate it, zero so perhaps it is a trivial case (again probably misusing terminology, sorry)?
6. What about if $r$ depends on $t$, time, what is it then? $\frac{\partial r}{\partial t}$ is now something. But would the extra dimension be now $t$ instead of $r$?

Please, use the term dimension in your replies, I feel it is very important to even understand what a sphere is, let alone a dimension.

[Conclusion] The general solution is still open. The $r$ is a homothetic transformation. The terms $(\theta, \phi)$ are apparently tied to the manifold, (ALERT! perhaps misusing terminology). The acceptance of a question does not mean that all things appeared here are solved. Zhen`s solution is really the solution to look at, big thanks for it!

Here is a diagram that might help you to visually understand that the element of area $\text{d}A = r^2 \sin \phi \, \text{d} \phi \, \text{d} \theta.$

Note that in the diagram $0 \le \phi \le \pi $ and $ 0 \le \theta \le 2 \pi.$

The easiest proof is to use the formula $$A = \int_S \sqrt{\det g} \, du \, dv $$ from Riemannian geometry, where $g = \begin{pmatrix} E & F \\ F & G \end{pmatrix}$ is the metric tensor / first fundamental form and $(u, v)$ are coordinates for the manifold $S$. Unfortunately this doesn't really explain anything, because now the question is how to justify this formula!

So I'll attempt to justify it now. First, we need to restrict to the case where we already know how to calculate areas, namely subsets of the Euclidean plane. We know that the area of a subset $S \subset \mathbb{R}^2$ is just $$A = \int_S dx \, dy$$ where $(x, y)$ are the standard Cartesian coordinates. So, for example, the area of a square with corners $(0, 0), (a, 0), (a, b), (0, b)$ is just $$A = \int_0^b \int_0^a dx \, dy = ab$$ But what about shapes which we can't easily parametrise using Cartesian coordinates? The change of variables formula from multivariable calculus tells us we need to introduce a factor — the Jacobian determinant — when we do this. For example, if we use plane polar coordinates $(r, \theta)$, then $$\begin{aligned} x & = r \cos \theta \\ y & = r \sin \theta \end{aligned}$$ so the Jacobian matrix is $$J = \begin{pmatrix} \cos \theta & -r \sin \theta \\ \sin \theta & r \cos \theta \end{pmatrix}$$ and the Jacobian determinant is $\det J = r$. Hence, the polar-coordinate form of the general formula is $$A = \int_S r \, dr \, d\theta$$ and we can use this to calculate, say, the area of a circle of radius $a$: $$A = \int_0^{2\pi} \int_{0}^{a} r \, dr \, d\theta = \pi a^2$$

But what does any of this have to do with $g$? Well, it turns out that $g = J^T J$, so for example $$g = \begin{pmatrix} 1 & 0 \\ 0 & r^2 \end{pmatrix}$$ for polar coordinates. Now, by the various properties of the determinant, it follows that in general $\det g = (\det J)^2$, i.e. $\det J = \sqrt{\det g}$. So this justifies the formula at the very beginning of my post.

Now let's apply this formula to the sphere. We have spherical polar coordinates $(\phi, \theta)$ such that $$\begin{aligned} x & = r \cos \theta \sin \phi \\ y & = r \sin \theta \sin \phi \\ z & = r \cos \phi \end{aligned}$$ and this gives the Jacobian $$J = \begin{pmatrix} r \cos \theta \cos \phi & -r \sin \theta \sin \phi \\ r \sin \theta \cos \phi & r \cos \theta \sin \phi \\ -r \sin \phi & 0 \end{pmatrix}$$ so the metric is $$g = J^T J = \begin{pmatrix} r^2 & 0 \\ 0 & (r \sin \phi)^2 \end{pmatrix}$$ hence the area of the sphere is just $$A = \int_0^{2\pi} \int_0^\pi r^2 \sin \phi \, d\phi \, d\theta = 4 \pi r^2$$ as we all know.

I'll try and answer 1, 2 and 3.. It's probably best to make a drawing, but you can also understand it 'algebraically'. An area element "$d\phi \; d\theta$" close to one of the poles is really small, tending to zero as you approach the North or South pole of the sphere. Close to the equator, the area tends to resemble a flat surface. This corresponds (qualitavely) to the factor $\sin(\phi)$, which goes to zero as $\phi$ goes to 0 or $\pi$, and equals 1 for $\pi / 2$.

There is no $\theta$-dependence, since the area element does not depend on your $\theta$-coordinate (i.e. walking parallel to the equator should not change your area element). Your coordinates are axisymmetrical with respect to the axe that goes through both poles.

Finally, why should the result depend on some other length $d$? It doesn't appear in the problem description, so finding it in your final result should be very suspect (all you're doing is some calculus, there's no exotic math involved which could introduce such a parameter).

The surface area of a sphere would be $r^2 \int_{0}^{2 \pi} \int_{0}^{\pi} \sin(\phi) \ d \phi \ d \theta$. This is in spherical coordinates.