From the Riemannian metric, you can calculate the curvature, which reflects, roughly speaking, how non-flat the space is. For a Riemannian 2-manifold, this is just a scalar function, namely the Gauss curvature. This tells you how curved the space is itself intrinsically.
The Riemannian metric tells you how to measure the length of a tangent vector. This in turn allows you to compute the length of a curve. By minimizing the lengths of curves joining two given points, you can define the distance between any two points. This gives the Riemannian manifold a metric space structure. If it is complete, then there will be a curve whose length is equal to the distance between the two points. That's called a length-minimizing geodesic.
The metric also allows you to define the angle between two curves that meet at a common point. You can now define the concept of a right geodesic triangle and ask whether it satisfies the Pythagorean theorem. The fundamental fact is that all geodesic right triangles satisfy the Pythagorean theorem if and only if the Gauss curvature is zero everywhere.
This is all done intrinsically, which means it's all done using only the Riemannian metric and otherwise not on how the surface looks in $\mathbb{R}^3$.
Knowing Gauss curvature does give information about what the surface looks like in $\mathbb{R}^3$ but the shape is not uniquely determined. You can see this already when $K$ is zero everywhere. Any cylindrical surface is flat because you can change the shape of an inflexible piece of paper by bending it. Right triangles drawn on the paper still satisfy the Pythagorean theorem even if you bend the paper.
If you're learning differential geometry, you have or will learn how to calculate the Gauss curvature of a sphere and the spherical version of the Pythagorean theorem. More generally, there is a Pythagorean theorem for any curved surface, where the "error term" is the integral of the Gauss curvature over the interior of the triangle. This is a version of the Gauss-Bonnet theorem for a geodesic right triangle.
Whenever you have a tensor field $T$ on a manifold $M$ (let's say twice-covariant, for simplicity), and coordinates $(x^1,\ldots, x^n)$ on some open subset of $M$, the components of $T$ relative to this coordinate system are defined by $$T_{ij} = T\left(\frac{\partial}{\partial x^i} ,\frac{\partial}{\partial x^j}\right),\qquad \mbox{ for }i,j=1,\ldots,n,$$and we have $$T = \sum_{i,j=1}^n T_{ij}\,{\rm d}x^i\otimes {\rm d}x^j.$$
Here, $M = \Bbb S^2$, $T = \iota^*g_E$, and the coordinates $(\theta,\varphi)$ are given by the inverse of the parametrization $X$. By taking derivatives, we see that $$\iota_*\left(\frac{\partial}{\partial \theta}\bigg|_{X(\theta,\varphi)}\right) = \frac{\partial X}{\partial\theta}(\theta,\varphi) = -\sin\theta\sin\varphi \frac{\partial}{\partial x}\bigg|_{X(\theta,\varphi)} + \cos\theta\sin\varphi\frac{\partial}{\partial y}\bigg|_{X(\theta,\varphi)}, $$and so $$(\iota^*g_E)\left(\frac{\partial}{\partial\theta},\frac{\partial}{\partial\theta}\right) = g_E\left(\frac{\partial X}{\partial\theta},\frac{\partial X}{\partial \theta}\right) = \sin^2\varphi.$$This means that $$\iota^*g_E = \sin^2\varphi\,{\rm d}\theta \otimes {\rm d}\theta + \cdots.$$
You can compute $\iota_* (\partial/\partial \varphi) = \partial X/\partial\varphi$ and the remaining components and inner products like above.
Best Answer
Essentially, the metric will have the form $${\rm d}s^2 = E(\theta,\varphi){\rm d}\theta^2 +2F(\theta,\varphi){\rm d}\theta\,{\rm d}\varphi + G(\theta,\varphi){\rm d}\varphi^2,$$where $$E(\theta,\phi)=\left\langle\frac{\partial\phi}{\partial \theta}(\theta,\varphi),\frac{\partial\phi}{\partial \theta}(\theta,\varphi)\right\rangle,\quad F(\theta,\varphi)=\left\langle\frac{\partial\phi}{\partial \theta}(\theta,\varphi),\frac{\partial\phi}{\partial \varphi}(\theta,\varphi)\right\rangle,\quad\mbox{and}\quad G(\theta,\varphi)=\left\langle\frac{\partial\phi}{\partial \varphi}(\theta,\varphi),\frac{\partial\phi}{\partial \varphi}(\theta,\varphi)\right\rangle,$$where these inner products are computed with the ambient metric. This is a general mechanism: if $(M^n,g)$ is a Riemannian manifold, $\iota:S \to M$ is a submanifold, and $(u^1,...,u^k)$ are coordinates for $S$, then $\iota^*g$ is described in these coordinates as $$\sum_{i,j=1}^k a_{ij}\,{\rm d}u^i\,{\rm d}u^j,$$where $$a_{ij} = g\left({\rm d}\iota\left(\frac{\partial}{\partial u^i}\right),{\rm d}\iota\left(\frac{\partial}{\partial u^j}\right)\right).$$