  Let $A_1$ and $A_2$ be expansive dilations, respectively, on ${\mathbb R}^n$
and ${\mathbb R}^m$. Let $\vec A\equiv(A_1, A_2)$ and $\mathcal A_p(\vec A)$ be
the class of product Muckenhoupt weights on ${\mathbb R}^n\times{\mathbb R}^m$
for $p\in(1, \infty]$. When $p\in(1, \infty)$ and $w\in{\mathcal A}_p(\vec A)$,
the authors characterize the weighted Lebesgue space $L^p_w({\mathbb
R}^n\times{\mathbb R}^m)$ via the anisotropic Lusin-area function associated
with $\vec A$. When $p\in(0, 1]$, $w\in {\mathcal A}_\infty(\vec A)$, the
authors introduce the weighted anisotropic product Hardy space $H^p_w({\mathbb
R}^n\times{\mathbb R}^m; \vec A)$ via the anisotropic Lusin-area function and
establish its atomic decomposition. Moreover, the authors prove that finite
atomic norm on a dense subspace of $H^p_w({\mathbb R}^n\times{\mathbb R}^m;\vec
A)$ is equivalent with the standard infinite atomic decomposition norm. As an
application, the authors prove that if $T$ is a sublinear operator and maps all
atoms into uniformly bounded elements of a quasi-Banach space $\mathcal B $,
then $T$ uniquely extends to a bounded sublinear operator from $H^p_w({\mathbb
R}^n\times{\mathbb R}^m;\vec A)$ to $\mathcal B$. The results of this paper
improve the existing results for weighted product Hardy spaces and are new even
in the unweighted anisotropic setting.
