Note that it will be enough to check the identity when $q$ is a prime power, in which case we can choose a field $F$ with $|F|=q$, and vector spaces $A$ and $B$ of dimension $a$ and $b$ over $F$.  In this context it is more natural to consider the function $\pi_q(k)=\prod_{i=1}^k(q^i-1)=(-1)^k(q;q)_k$ rather than $(q;q)_k$, because $\pi_q(k)$ is a positive integer and is more closely related to counting problems.  We then have 
$$ {\binom{n}{k}_q}=\frac{\pi_q(n)}{\pi_q(k)\pi_q(n-k)} $$ 
and the stated identity becomes 
$$ q^{ab} = \sum_{k\geq 0} q^{\binom{k}{2}}\binom{a}{k}_q\binom{b}{k}_q\pi_q(k). $$
The left hand side is $|\text{Hom}(A,B)|$.  To give a homomorphism $\alpha\colon A\to B$ we have to choose a rank $k$, a kernel $U\leq A$ of dimension $a-k$, an image $V\leq B$ of dimension $k$, and an isomorphism $\alpha_1\colon A/U\to V$.  The number of choices for $U$ is $\binom{a}{a-k}_q=\binom{a}{k}_q$.  The number of choices for $V$ is $\binom{b}{k}_q$.  The number of choices for $\alpha_1$ is $|GL_k(F)|=q^{\binom{k}{2}}\pi_q(k)$.