[guided]We first reduce the statement to a calculation in one complex variable. The assumption $df_q \neq 0$ for every $q \in U$ means that, at each point of $U$, the complex-[linear map](/page/Linear%20Map)
$df_q:T_q^{1,0}\Sigma \to T_{f(q)}^{1,0}Y$
is nonzero. Since $T_q^{1,0}\Sigma$ is one-dimensional, this implies that $f^*\omega_Y$ is positive on nonzero tangent vectors at $q$. Therefore $f^*\omega_Y$ is a smooth Kähler metric on $U$.
Fix $p \in U$, choose a holomorphic coordinate $z$ on $\Sigma$ with $z(p)=0$, and choose holomorphic coordinates $w_1,\dots,w_m$ on $Y$ near $f(p)$. After shrinking the coordinate neighbourhood of $p$, the map $f$ is represented by a holomorphic map
\begin{align*}
F: z(V) &\to \mathbb{C}^m, \\
\zeta &\mapsto (F_1(\zeta),\dots,F_m(\zeta)).
\end{align*}
In these coordinates the target metric is written as
\begin{align*}
\omega_Y
=
i\sum_{\alpha,\beta=1}^m g_{\alpha\bar\beta}(w)\,dw_\alpha\wedge d\bar w_\beta,
\end{align*}
where each $g_{\alpha\bar\beta}:W\to\mathbb{C}$ is smooth and the Hermitian matrix
$(g_{\alpha\bar\beta})$ is positive definite.
Pulling back $dw_\alpha$ gives
\begin{align*}
f^*(dw_\alpha)=\frac{\partial F_\alpha}{\partial z}\,dz,
\end{align*}
because $F_\alpha$ is holomorphic. Hence
\begin{align*}
f^*\omega_Y
=
i\sum_{\alpha,\beta=1}^m
g_{\alpha\bar\beta}(F(z))
\frac{\partial F_\alpha}{\partial z}
\overline{\frac{\partial F_\beta}{\partial z}}
\,dz\wedge d\bar z.
\end{align*}
Thus the pullback metric has the conformal coefficient
\begin{align*}
\lambda(\zeta)
=
\sum_{\alpha,\beta=1}^m
g_{\alpha\bar\beta}(F(\zeta))\,
\frac{\partial F_\alpha}{\partial z}(\zeta)\,
\overline{\frac{\partial F_\beta}{\partial z}(\zeta)}.
\end{align*}
The positivity of $\lambda$ follows from positive definiteness of $(g_{\alpha\bar\beta})$ and from $df \neq 0$ on $U$.[/guided]