In general if $M$ is a smooth $n$-dimensional manifold and $g_1,g_2$ are two Riemannian metrics which are conformal, say $g_2=fg_1$ for some smooth positive function $f:M\to\Bbb{R}$, then the induced volume measures are related as
\begin{align}
dV_{g_2}&=f^{n/2}\cdot dV_{g_1},\tag{$1$}
\end{align}
meaning that for all (Lebesgue-measurable) sets $A\subset M$, we have
\begin{align}
\text{vol}_{g_2}(A)&=\int_Af^{n/2}\,dV_{g_1}.\tag{$2$}
\end{align}
Now, unless $f$ is constant, you can’t simplify this any further to get something like $\text{vol}_{g_2}(A)=f^{n/2}\cdot\text{vol}_{g_1}(A)$.
In order to prove this result, you simply need to recall that the volume measure $dV_g$ is essentially defined such that in local coordinates $(U,x)$, we have for all (Lebesgue-measurable) $A\subset U$, $\text{vol}_g(A):=\int_{x[A]}\sqrt{|\det g_{(x), ij}|}\,dx$. Note that if you don’t know the language of measure theory, that’s perfectly fine; you can simply assume $M$ is oriented, and focus on the Riemannian volume form, and you’ll still get the same result.
Now if you have two conformally related metrics, then you simply note that the determinant is simplified as follows:
\begin{align}
\sqrt{|\det (g_2)_{(x),ij}|}=\sqrt{|\det f\cdot (g_1)_{(x), ij}|}=\sqrt{f^n\cdot |\det (g_1)_{(x),ij}|}=f^{n/2}\sqrt{|\det (g_1)_{(x),ij}|},
\end{align}
where I simply used that if a matrix is multiplied by $f$, then by multi-linearity of the determinant in each of the $n$ columns, the determinant is multiplied by $f^n$. So that’s where the $f^{n/2}$ comes from in $(1)$.
Finally, in your special case, we have $n=2$ so $f^{n/2}=f$, hence
\begin{align}
S_g(A)&=\int_Af(x,y)\,dx\,dy.
\end{align}
Besides this, there isn’t a simpler way of expressing it. Even if you assumed the conformal factor is radial, or harmonic or whatever. You still need to evaluate the integral, and you need to know how the set $A$ looks.
But, one special case I can give you is if $f(x,y)=\lambda(x^2+y^2)$ is radial, and $A$ is a Euclidean annulus centered at the origin with radii $0\leq r_1<r_2\leq\infty$, then
\begin{align}
S_g(A)&=\int_Af(x,y)\,dx\,dy=2\pi\int_{r_1}^{r_2}\lambda(r^2)\cdot r\,dr,
\end{align}
which thus leaves you with a single variable integral to deal with (which could still be very complicated). Really, I would prefer you normalize the radial function as $f(x,y)=\phi(\sqrt{x^2+y^2})$, so the formula reads $S_g(A)=2\pi\int_{r_1}^{r_2}\phi(r)\cdot r\,dr$.