The group of orthogonal transformations carries an invariant probability measure. This means that we can average a function over the group in a natural way. In particular, if $f$ is a function on the sphere and $θ$ is some point on the sphere, the average over orthogonal $U$ of the value $f(Uθ)$ is just the average of $f$ on the sphere: averaging over $U$ mimics averaging over the sphere: $$\text{avg}_U f(Uθ) = \int_{S^{n-1}}f (\phi) dσ(\phi)$$
See Pg. $22-23$ of these notes for context.
Is there an easy way to understand/prove this? I'm not able to figure out what it means!
My intuition:
Consider $\theta$ on the sphere $S^{n-1}$. Orthogonal transformations rotate this without changing the length. All orthogonal transformations send this point to someplace on the sphere, and for any two points on the sphere, we can always find an orthogonal transformation that relates them. Hence, averaging over all orthogonal transformations is the same as averaging over the sphere.
Is that correct?