The commutative Gelfand-Naimark theorem tells us that every unital commutative C* algebra is isometrically isomorphic to the space of continuous functions on its maximal ideal space. The non- commutative GNS construction, on the other hand, tells us that every C* algebra is isometrically embedded into a * subalgebra of $\mathfrak{B}(\mathbb{H})$ for some Hilbert space $\mathbb{H}$.
My question is the following: how does the commutative case fit into the non commutative one?
I could see that given a commutative C* algebra $C(X)$, where $X$ is a locally compact Hausdorff space, one can find a non zero positive linear functional $F$ on $C(X)$ such that $F(e)=1$, where $e$ is the identity function on $C(X)$. Then using Riesz Representation theorem, one can find a corresponding Borel $\sigma$ algebra on $X$ and a unique probability measure $\mu$ such that $F(f)= \int f\,d\mu \quad \forall f\in C(X)$. From here it is easy to deduce that $C(X)$ is isometrically embedded into a * subalgebra of $\mathfrak{B}(\mathcal{L^2}(X,\mu))$ via multiplication operators.
My question is: Is there a canonical way of finding such a probability measure? Clearly if we choose a different positive linear functional, we get a different probability measure. So is there some unique/ natural way of doing this? Also, is this where the Hahn-Hellinger theorem comes in?
I would be grateful for some clarity on this, and also references if possible.