This is a doubt that I carry since my PDE classes.
Some background (skippable):
In the multivariable calculus course at my university we made all sorts of standard calculations involving surface and volume integrals in $R^3$, jacobians and the generalizations of the fundamental theorem of calculus. In order to make those calculations we had to parametrize domains and calculate differentials.
A couple of years later I took a PDE course. We worked with Evans' Partial differential equations book. This was my first experience with calculus in $\mathbb R^n$ and manipulations like $$\text{average}\int_{B(x,r)}f(y)\,dy= \text{average}\int_{B(0,1)}f(x+rz)\,dz.$$ This was an ordinary change of variables. $y=x+rz,\,\,dy=r^n\,dz$ and the mystery was solved. Like in that case, I was able to justify most of these formal manipulations after disentangling definitions.
That aside, I found these quick formal calculations to be very powerful.
However,
I realized that I wasn't able to justify this: $$\text{average} \int_{\partial B(x,r)}f(y)dS(y)= \text{average}\int_{\partial B(0,1)}f(x+rz)\,dS(z).$$ I have some vague idea of what's happening: the same substitution as before, but this time the jacobian is $r^{n-1}$ because the transformation is actually happening between regions which "lack one dimension". Also, I see some kind of pattern: a piece of arc-length in the plane is $r\,d\theta$, a piece of sphere-area is $r^2 \sin\theta \, d\phi \,d\theta$, "and so on". Maybe some measure-theoretic argument can help me: I know, roughly speaking, that for any measure $\mu$, $$\int_\Omega f\circ \phi \,d\mu=\int_{\phi(\Omega)} f \, d(\mu\circ\phi^{-1}).$$ I'd say $\phi(z)=(z-x)/r$ and $\phi^{-1}(y)=ry+x$, but I actually don't know how $dS(y)$ looks like "as a measure" (It's not a product measure or a restriction of one, but it somehow relates to Lebesgue's in $\mathbb R^n$...). Why would I conclude that $dS(y)\circ \phi^{-1}=r^{n-1}dS(z)$? I have an intuition, but either I lack the mathematical concepts and definitions to express it or I'm just too confused. Is there some theory that I could learn in order to understand? Maybe something about the measure $dS$. Is it expressible in terms of the Lebesgue measure in some way? Or set-theoretically, maybe, without having to resort to $n-1$ parameters and complicated relations?
Maybe all of this would not have been a problem if I had ever mastered n-dimensional spherical coordinates. But even so, more generally, is there a way of changing variables when I'm integrating over a subregion of "dimension$<n$" without necessarily parametrizing?
Sorry for the vagueness, but I don't really know what to ask for exactly.
Note: I saw some of the answers to this post, but none of them were deep enough in the direction I intend.
Note II: If there are no general methods or theories, maybe restricting to linear transformations, to Lebesgue measure exclusively, or to subregions defined by simple expressions like $g(x)=C$ or $g(|x|)=C$ could get me somewhere.
Edit: I have not yet studied differential geometry, which has been mentioned in a comment. I added it to the tags.