I think you can get the result blowing-down the solution. I'll just sketch the argument since I don't want to be confusing with the details.
Set $u:=1-f$, then $u$ solves $-\Delta u = \frac{h}{2}(1-u)$ and satisfies $\lim_{|x|\to \infty} u(x)=0$. For $r>0$ consider $u_r(x):= ru(rx)$, you can check easily that $u_r$ solves
$$-\Delta u_r =\frac{r^3}{2}h(rx)(1-u_r) =: g_r(x)\,. $$
Since $h$ has compact support you have $-\Delta u_r= g_r \to 0$ in $C^2_{\rm loc}(\mathbb{R}^3 \setminus \{0\})$ and moreover
$$ \int_{\mathbb{R}^3} g_r \, dx = \frac{r^3}{2} \int_{\mathbb{R}^3} h(rx)(1-u_r(x)) \, dx = \frac{1}{2}\int_{\mathbb{R}^3} h(1-u) = \frac{1}{2}\int_{\mathbb{R}^3} hf =: c_0 \,.$$
Hence $g_r \to c_0\, \delta_{\{x=0\}}$ weakly in the sense of distributions, and from here by elliptic estimates you get that $u_r \to c_0 \, G_3$ where $G_3=C_n/|x|$ is the fundamental solution in $\mathbb{R}^3$. This holds since $G_3$ is the unique solution to $-\Delta G_3=\delta_{\{x=0\}}$ that goes to zero as $|x|\to \infty$. This implies that, for example, in the annulus $A_{1,2}:=B_2(0)\setminus B_1(0)$ there holds
$$ \sup_{x\in A_{1,2}} \left|ru(rx) -\frac{c_0 C_n}{|x|} \right| \to 0$$
as $r\to \infty$, or equivalently
$$ r \sup_{x \in A_{r,2r}} \left|u(x) -\frac{c_0 C_n}{|x|} \right| \to 0 \,.$$
This implies what you ask.