A well known change of variable theorem of integral calculus says the following.
Let $X$ and $Y$ be two open sets in $\mathbb{R}^2$ and $\phi: X \rightarrow Y$ be a diffeomorphism. Then \begin{eqnarray}\tag{1} \int\limits_{Y} f(y_1,y_2) dy_1 dy_2=\int\limits_{X} f(\phi(x_1,x_2)) \left|J_{\phi}(x_1,x_2)\right| dx_1 dx_2, \end{eqnarray} where $J_{\phi}$ denotes the Jacobian matrix of $\phi.$
What is the corresponding change of variable formula if the integrand involves derivatives. More precisely, for $f \in C^1$ what is the change of variable formula if the integrand is \begin{eqnarray} \int\limits_{Y} \left(\partial_1f(y_1,y_2)+\partial_2f(y_1,y_2)\right) dy_1 dy_2=\int\limits_{X} \cdot\cdot\cdot\cdot dx_1 dx_2. \end{eqnarray} How to prove the change of variable formula for such integrand using the change of formula (1) ?