There are many rigorous works showing this semiclassical limit. However I'll show you an intuitive and qualitative method to see this limit.
Consider the matrix:
$$\frac{S_z}{\sqrt{s (s+1)}}$$
Where $S_z$ is the spin projection onto the $z$ axis. The eigenvalues of this matrix give the normalized heights of the small circles of the allowed spin states. Of course in the representation $s$, $S_z$ has the form:
$$S_z = \mathrm{diag}(-s, -s+1, ., ., ., ., s-1, s)$$
For example for $s=2$, the allowed normalized heights are: $\frac{1}{\sqrt{6}}(-2, -1, 0, 1, 2)$.
When the spin representation gets larger and larger, the allowed normalized spins will get denser and denser on the $z$ axis, until they reach a continuum as $s \rightarrow \infty$. We can identify this matrix in the limit with the $z$ coordinate on the sphere.
We can do the same thing with $S_x$ and $S_y$, and obtain the $x$ and $y$ coordinates. One can easily check that due to the normalization, the identity:
$$S_x^2 + S_y^2 + S_z^2 = 1$$
is satisfied throughout the limiting process.
Any function of the coordinates can be built from these infinite dimensional matrices (say $x^2+y+z$).
What needs a little more work to show, but is absolutely correct, is that the commutator of any two functions expressed as infinite dimensional matrices is exactly the Poisson bracket between the corresponding classical functions on the sphere, where the symplectic structure is the sphere area form.
Thus, we can see that this construction resulted in the Poisson algebra of the sphere expressed as infinite dimensional matrices.
One special example to see the correspondence between the infinite dimensional matrices and functions on the sphere is to compute the variance of the height function using both representations:
In the infinite dimensional representation:
$$\langle z^2 \rangle = \lim_{s \to \infty} \frac{1}{s(s+1)}\frac{2 \sum_0^s i^2}{2s+1}= \lim_{s \to \infty}\frac{2s(s+1)(2s+1)}{ 6s(s+1)(2s+1)} = \frac{1}{3}$$
In the continuum representation
$$\langle z^2 \rangle = \frac{\int z^2 dV_{S^2}}{ V_{S^2}} = \frac{2 \pi \int_{-1}^{1} \cos^2\theta \sin\theta d\theta }{ 4 \pi} = \frac{1}{3}$$