Computing pullbacks is always a bit of a hairy task. The key point to notice here is that if we write $z\mapsto w = i\frac{1-z}{1+z}$ express the non-constant constant portion of the metric $g_N$ at $w$ in terms of $z$, we find
$$ \frac{1}{\text{Im}(w)^2} = \frac{(1+z)^2(1+\bar{z})^2}{(1-|z|^2)^2}$$
because
\begin{align*}
\text{Im}(w) &= \text{Im}\left(i\frac{1-z}{1+z}\right)\\
&= \frac{1}{2i}\left[i\frac{1-z}{1+z} - \overline{\left(i\frac{1-z}{1+z}\right)}\right]\\
&= \frac{1}{2i}\left[i\frac{1-z}{1+z} + i\frac{1-\bar{z}}{1+\bar{z}}\right]\\
&= \frac{1}{2}\left[\frac{(1-z)(1+\bar{z}) + (1-\bar{z})(1+z)}{(1+z)(1+\bar{z})}\right]\\
&= \frac{1}{2}\left[\frac{1 -z + \bar{z} - z\bar{z} + 1- \bar{z} + z - z\bar{z} }{(1+z)(1+\bar{z})}\right]\\
&= \frac{1- |z|^2}{(1+z)(1+\bar{z})}
\end{align*}
This factor of $(1-|z|^2)$ is the first indication that we are on the right track! So, when formally computing the pullback of your metric you'll begin to see
\begin{align*}
f^*(g_N)_z &= f^*\left(\frac{1}{Im(z)^2}dz\otimes d\bar{z}\right)\\
&:= (g_N)_{f(z)}\\
&= \frac{1}{\text{Im}(w)}dw \otimes d\bar{w} \\
&= \frac{(1+z)^2(1+\bar{z})^2}{(1-|z|^2)^2}dw \otimes d\bar{w}\\
\end{align*}
It remains to compute $|dw|^2 = dw\otimes d\bar{w}$ which we can only hope will cancel with those nasty copies of $(1+z)(1+\bar{z})$ in the numerator. This will look something like
$$dw = d\left(i\frac{1-z}{1+z}\right) = i\frac{-(1+z)-(1-z)}{(1+z)^2}dz = -\frac{2i}{(1+z)^2}dz$$
and similarly $$d\bar{w} = d\left(\overline{i\frac{1-z}{1+z}}\right) = -i d\left(\frac{1-\bar{z}}{1+\bar{z}}\right) = \frac{2i}{(1+\bar{z})^2}d\bar{z}.$$
Combining these yields the pullback as
$$f^*(g_N)_z = \frac{4}{(1-|z|^2)^2}dz\otimes d\bar{z} = 4g_M.$$
Geometrically this is classical and well recognized conformal isometry is between the Poincaré disk and the upper-half plane!