With $a \in \mathbb{C}$, let $h: D_1 \to D_2$ and $f: D_2\setminus{\{h(a)\}} \to \mathbb{C}$ be analytic functions. Also require $h'(a) \neq 0$.
Claim: $$ \operatorname{Res}(f; h(a)) = \operatorname{Res}((f \circ h)h'; a)$$
Attempt
For $\rho$ small enough, we have: $$ \operatorname{Res}(f; h(a)) = \frac{1}{2\pi i}\int_{|u-h(a)|=\rho} f(u)\ \mathrm{d}u$$ and \begin{align}\operatorname{Res}((f \circ h)h', a)& = \frac{1}{2\pi i} \int_{|u-a|=\rho} (f \circ h) h'\ \mathrm{d}u = \frac{1}{2\pi i} \int_0^{2\pi} f(h(a + \rho e^{it})) h'(a + \rho e^{it})\rho e^{it} i\ \mathrm{d}t\\&= \frac{1}{2\pi i} \int_\alpha f(u)\ \mathrm{d}u \end{align}
where $\alpha: [0, 2\pi] \to \mathbb{C}, \alpha(t) = h(a+\rho e^{it})$.
We see that the curve of the LHS defined by $|u-h(a)|=\rho$ (i.e. $h(a) + \rho e^{it}$) is "approximated" by the curve defined by $\alpha$ for $\rho \to 0$ because analytic functions approximately retain circles (see here).
(Note that we already used $h'(a) \neq 0$: Writing the residue for the RHS was only possible because $(f \circ h)$ is analytic on $U_\rho(a)\setminus{\{a\}}$. Since $h'(a) \neq 0$, $h$ is injective on a small enough neighborhood around $a$.)
Capturing the approximation
I tried to generalize the claim above by introducing the following
Lemma: If $f: D \to \mathbb{C}$ is a continuous function and $(\alpha_n)_n, (\beta_n)_n \subseteq D$ series of piecewise smooth curves $\alpha_n, \beta_n: [a, b] \to D$ such that $(f \circ \alpha_n) \alpha_n' - (f \circ \beta_n) \beta_n' \to 0$ uniformly, then:
If for a fixed $z \in \mathbb{C}$ the equation holds: $z = \int_{\alpha_n} f(u)\ \mathrm{d}u$ for all $n \ge 0$
then $z = \int_{\beta_n} f(u)\ \mathrm{d}u$.
Proof: For every $\varepsilon > 0$ we have $|z - \int_{\beta_n} f(u)\ \mathrm{d}u| = \left|\int_a^b ((f \circ \alpha_n) \alpha_n' - (f \circ \beta_n) \beta_n')(t)\ \mathrm{d}t\right| \le (b - a) \sup_{t \in [a, b]} \left|...\right| \le \varepsilon$ for a chosen $n$ big enough.
Using the lemma
Choose $\alpha_n(t) := h(a + \rho_n e^{it})$ (from the RHS), $\beta_n(t) := h(a) + \rho_n e^{it}$ (from the LHS). We have $$|\alpha_n(t) - \beta_n(t)| = \left|f(h(a + \rho_n e^{it})) h'(a + \rho_n e^{it})\rho_n e^{it} i - f(h(a) + \rho_n e^{it}) \rho_n e^{it} i\right| = |\rho_n| \cdot |*|$$ where $$|*| = \left|f(h(a + \rho_n e^{it})) h'(a + \rho_n e^{it}) - f(h(a) + \rho_n e^{it})\right|.$$
If $|*|$ was bounded, we would be able to apply the theorem. This is the place where I am currently stuck. I tried expressing $f$ with its Laurent series around $h(a)$, but that did not lead to easier-looking equations.