[step:Verify $\Psi$ is a $*$-homomorphism]
We check linearity, the involution-preservation $\Psi(\bar{f}) = \Psi(f)^*$, and multiplicativity $\Psi(fg) = \Psi(f) \Psi(g)$ for $f, g \in L^\infty(K, \mathcal{B}(K))$.
*Linearity.* The form $\theta_f$ depends linearly on $f$ (linearity of the integral in the integrand), and the operator-from-form map is itself linear, so $\Psi(\alpha f + g) = \alpha \Psi(f) + \Psi(g)$.
*Involution.* For $x, y \in H$,
\begin{align*}
(\Psi(f)^* x, y)_H = \overline{(\Psi(f) y, x)_H} = \overline{\int_K f \, d\mu_{y, x}} = \int_K \bar{f} \, d\overline{\mu_{y, x}} = \int_K \bar{f} \, d\mu_{x,y} = (\Psi(\bar{f}) x, y)_H,
\end{align*}
using $\overline{\mu_{y,x}} = \mu_{x,y}$ (Step 2). Holding for all $x, y$, this gives $\Psi(f)^* = \Psi(\bar{f})$.
*Multiplicativity.* This is the deepest step; we prove it by showing two pushforward identities of measures.
[claim:For all $T \in A$ and $x, y \in H$, $\hat{T} \, d\mu_{x,y} = d\mu_{Tx, y}$ as regular complex Borel measures on $K$.]
[proof]
Both measures are regular complex Borel measures on $K$. By RRT uniqueness it suffices to show
\begin{align*}
\int_K f \, \hat{T} \, d\mu_{x,y} = \int_K f \, d\mu_{Tx, y} \quad \text{for all } f \in C(K).
\end{align*}
Both sides involve continuous integrands, and both are linear continuous functionals of $f \in C(K)$. They agree provided they agree on the dense subspace $\{f = \hat{T}'\} = C(K)$ itself — i.e.\ on all of $C(K)$. We compute both sides using the defining identity of $\mu$:
*RHS.* $\int_K f \, d\mu_{Tx, y} = (\gamma(f) Tx, y)_H$.
*LHS.* For $f = \hat{T}'$ with $T' \in A$, $f \cdot \hat{T} = \hat{T}' \cdot \hat{T} = \widehat{T' T}$ (Gelfand transform is multiplicative on $A$). Hence
\begin{align*}
\int_K f \cdot \hat{T} \, d\mu_{x,y} = \int_K \widehat{T'T} \, d\mu_{x,y} = ((T'T)x, y)_H = (T'(Tx), y)_H = (\gamma(\hat{T}')(Tx), y)_H = (\gamma(f)(Tx), y)_H.
\end{align*}
Since $\gamma$ inverts the Gelfand transform on $A$, every $f \in C(K)$ has the form $\hat{T}'$ for a unique $T' \in A$, so the identity holds for all $f \in C(K)$. RRT uniqueness gives $\hat{T} \, d\mu_{x,y} = d\mu_{Tx, y}$.
[/proof]
[/claim]
[claim:For all $f \in L^\infty(K, \mathcal{B}(K))$ and $x, y \in H$, $f \, d\mu_{x,y} = d\mu_{x, \Psi(f)^* y}$.]
[proof]
Both measures are regular complex Borel measures on $K$. By RRT uniqueness it suffices to test against $g \in C(K)$:
\begin{align*}
\int_K g \, d\mu_{x, \Psi(f)^* y} = (\gamma(g)x, \Psi(f)^* y)_H = (\Psi(f) \gamma(g) x, y)_H.
\end{align*}
We need this to equal $\int_K g f \, d\mu_{x,y}$, which by the defining identity for $\Psi$ equals $(\Psi(gf) x, y)_H$. So we must show $\Psi(f) \gamma(g) = \Psi(gf)$ as operators in $\mathcal{L}(H)$, applied to $x$.
Equivalently — using $\gamma(g) = \Psi(g)$ from Step 3 — we must show $\Psi(f) \Psi(g) = \Psi(gf)$ when $g \in C(K)$. We prove this by another application of RRT.
For $T \in A$ (so $g = \hat{T}$ for some $T$), by the previous claim, for any $w, z \in H$,
\begin{align*}
(\Psi(f) T w, z)_H = \int_K f \, d\mu_{Tw, z} = \int_K f \cdot \hat{T} \, d\mu_{w, z} = (\Psi(f \hat{T}) w, z)_H,
\end{align*}
holding for all $w, z$, so $\Psi(f) T = \Psi(f \hat{T})$, i.e.\ $\Psi(f) \gamma(\hat{T}) = \Psi(f \hat{T})$. Since every $g \in C(K)$ has the form $\hat{T}$, this gives $\Psi(f) \gamma(g) = \Psi(fg)$ for all $g \in C(K)$.
Since pointwise multiplication on $L^\infty(K, \mathcal{B}(K))$ is commutative, $fg = gf$, so $\Psi(fg) = \Psi(gf)$. Hence $\Psi(f) \gamma(g) = \Psi(fg) = \Psi(gf)$, which gives
\begin{align*}
(\Psi(f) \gamma(g) x, y)_H = (\Psi(gf) x, y)_H = \int_K gf \, d\mu_{x,y},
\end{align*}
proving the claim.
[/proof]
[/claim]
*Multiplicativity, finally.* Combine the two claims. For $f, g \in L^\infty(K, \mathcal{B}(K))$ and $x, y \in H$, applying the second claim with the integrand $g$ gives $g \, d\mu_{x,y} = d\mu_{x, \Psi(g)^* y}$, hence
\begin{align*}
(\Psi(fg) x, y)_H = \int_K f g \, d\mu_{x,y} = \int_K f \, d\mu_{x, \Psi(g)^*y} = (\Psi(f) x, \Psi(g)^* y)_H = (\Psi(g) \Psi(f) x, y)_H.
\end{align*}
Holding for all $x, y \in H$, we obtain $\Psi(fg) = \Psi(g) \Psi(f)$. Since $L^\infty(K, \mathcal{B}(K))$ is a commutative algebra under pointwise multiplication, $fg = gf$, so $\Psi(g) \Psi(f) = \Psi(gf) = \Psi(fg)$. Reversing the roles of $f$ and $g$ in the displayed identity gives $\Psi(gf) = \Psi(f) \Psi(g)$ as well, so $\Psi(fg) = \Psi(f) \Psi(g) = \Psi(g) \Psi(f)$. Multiplicativity holds and the image $\Psi(L^\infty(K, \mathcal{B}(K)))$ is commutative.
Hence $\Psi$ is a unital $*$-homomorphism.
[/step]