I'll treat polar coordinates, just so that I can deal with fewer variables.
The important thing is that we have scalar fields $x, y, r, \theta$ satisfying the equations $x = r \cos \theta$ and $y = r \sin \theta$. At this point, it doesn't matter so much on what they are fields, just that they are, and have this relationship.
From this, we can compute equations that their differentials satisfy:
$$ dx = \cos \theta \, dr - r \sin \theta \, d\theta $$
$$ dy = \sin \theta \, dr + r \cos \theta \, d\theta $$
These are also called "1-forms". While you can multiply a scalar field (which can also be called 0-forms) by a 1-form in the "normal" way, you can't do the same with two 1-forms. Instead, there is something called the wedge product, and you can take the wedge product of an $m$-form with an $n$-form to get an $(m+n)$-form.
In particular, $dx \wedge dy$... or just $dx \, dy$ for short... is a $2$-form.
The wedge product behaves similarly to the ordinary product, with one key difference: if $\omega$ and $\eta$ are both $1$-forms, then $\omega \eta = -\eta \omega$. So you have to pay attention to the order in which things are multiplied. More generally, if $\omega$ is an $m$-form and $\eta$ is an $n$-form, then $\omega \eta = (-1)^{mn} \eta \omega$.
Note that this means $(dx)^2 = 0$, and similarly for any $1$-form.
In particular, you can expand products in the normal way. So we can write
$$ \begin{align}dx \, dy &= (\cos \theta \, dr - r \sin \theta \, d\theta)
(\sin \theta \, dr + r \cos \theta \, d\theta )
\\&= (\cos \theta \, dr) (\sin \theta \, dr) - (r \sin \theta \, d\theta)
(\sin \theta \, dr) + \ldots
\\&= \cos \theta \sin \theta \, (dr)^2 - r \sin^2 \theta \, d\theta \, dr + \ldots
\\&= (r \sin^2 \theta + r \cos^2 \theta) \, dr \, d\theta
\\ &= r \, dr \, d\theta
\end{align}$$
Now, sometimes we might consider two different spaces. We might consider $E$, the ordinary Euclidean plane, with coordinate functions (i.e. scalar fields) $x$ and $y$, with the property that any point satisfies
$$ P = (x_P, y_P) $$
where $x_P$ is notation for the value of $x$ at the point $P$: i.e. for $x(P)$ if we view $x$ as a function. So we might just say that a generic point is $(x,y)$.
We might also have the space $X$ which is again a Euclidean plane, but this one has the coordinate functions $r$ and $\theta$.
And, as you wanted to do, we have a function $X \to E$ defined by
$$ f(r, \theta) = (r \cos \theta, r \sin \theta) $$
Or more precisely,
$$ f(P) = f(r_P, \theta_P) = (r_P \cos \theta_P, r_P \sin \theta_P)$$
so that if $Q = f(P)$, then
$$ (x_Q, y_Q) = (r_P \cos \theta_P, r_P \sin \theta_P)$$
We can actually use this to define new scalar fields $\bar{x}$ and $\bar{y}$ on $X$, by the equations
$$\bar{x}_P = x_{f(P)} \qquad \qquad
\bar{y}_P = y_{f(P)} $$
Now, $\bar{x}, \bar{y}, r, \theta$ are all scalar fields on the same space $X$. Using $f$ in this way to take scalars on $E$ and turning them into scalars on $X$ is called a "pull back". This also extends to differential forms: in particular, $dx$ pulls back to $d\bar{x}$, and similarly.
So the area form $dx \, dy$ on the Euclidean plane $E$ pulls back to the $2$-form $d\bar{x} \, d\bar{y}$ on $X$. If $U,V$ are subsets of $X$ and $E$ respectively with the property that $f$ is a bijection from $U$ to $V$, then we can actually compute the integrals of a $2$-form (e.g. $dx \, dy$) over the region $V$ by instead computing the integral of the pullback (e.g. $d\bar{x} \, d\bar{y}$) over $U$.
This amounts to one means of doing a rigorous treatment of "change of variables".
By using the calculation above, this tells us
$$ \begin{align}
\iint_V f(x,y) \, dx \, dy &= \iint_U f(\bar{x}, \bar{y}) \, d\bar{x} \, d\bar{y}
\\ &= \iint_U f(r \cos \theta, r \sin \theta) r\, dr \, d\theta
\end{align}$$
This shows how we can use the framework of differential forms to rigorously treat change of variables.
A dual thing we can do is a push forward. If, rather than differentials, we work with "tangent vectors" -- e.g. a vector field $e_r$ on $X$ that, at every point, is a "unit" vector that points in the direction in which $r$ grows if you hold $\theta$ constant, and similarly for $e_\theta$.
then, you can push these forward, and get new vector fields $\bar{e}_r$ and $\bar{e}_\theta$ on $E$. These are more often written as $f_*(e_r)$, or maybe $df \cdot e_r$ (we can give $df$ a precise meaning as a vector of differential forms, if we wanted to).
This is a very useful tool for many things... but if you're trying to understand differentials so that you can integrate with them, I think it makes things much more complicated to try and express things in terms of this.