Let $D$ be a closed, bounded domain in $\mathbb R^2$ and let $\vec r(u,v) = \langle x(u,v), y(u,v), z(u,v) \rangle$ for $(u, v) \in D$ be a parametrization of a smooth surface $S \subseteq \mathbb R^3$. Then the area of $S$ is
$$\iint_D \| \vec r_u \times \vec r_v \| \, dA$$
I am interested in what might be considered a kind of converse of this: Is it always possible to interpret an integral $\iint_D f \, dA$ as the area of some surface?
More precisely, given a (non-negative) function $f:\mathbb R^2 \to \mathbb R$, defined on some domain $D$, under what conditions can we find a surface $S \subseteq \mathbb R^3$ and a parametrization $\vec r(u,v): D \to S$ such that
$$f(u,v) = \| \vec r_u \times \vec r_v \|$$ so that consequently $$\iint_D f \, dA = \iint_D \| \vec r_u \times \vec r_v \| \, dA \, ?$$