For two measurable sets $A,B$, let $p=\mu(A)$ and $q=\nu(B)$. Consider any coupling of a Bernoulli$(p)$ and a Bernoulli$(q)$, say, $C:\{0,1\}^2 \to [0,1]$. Then we can find a coupling $(X,Y)$ of $\mu$ and $\nu$ with
\begin{align}
(X\in A, Y\in B) = q(1,1),\qquad
(X\in A, Y\in B^c)= q(1,0),
\\(X\in A^c, Y\in B^)= q(0,1),\qquad
(X\in A^c, Y\in B^c)= q(0,0).
\end{align}
For convenience, denote $A^1=A$ and $A^0=A^c$, similarly $B^1=B$ and $B^0=B^c$.
Sample an iid sequence with the product measure $\mu\times \nu$, say $(X_t,Y_t)_{t\ge 1}$.
Sample $(U,V)\sim q$, return $(X_T,Y_T)$ where $T=\min\{t\ge 1: X_t\in A^U, Y_t\in B^V\}$. This is a way to construct conditional distributions given that a an independent pair belong to one of $(A,B)$, $(A^c,B)$, $(A, B^c)$, $(A^c,B^c)$.
The conditional distribution is given by $P(X_T\in D|U=u, V=v)=\mu(D\cap A^u)/\mu(A^u)$. The law of $X_T$ is the desired marginal:
$$
P(X_T\in D)=\sum_{u=0}^1\sum_{v=0}^1
P(X_T\in D|U=u, V=v)
q(u,v)=\sum_{u=0}^1\sum_{v=0}^1
\frac{\mu(D\cap A^u)}{\mu(A^u)}
q(u,v)
$$
which is $\mu(D)$ as the denominator cancels with $\sum_v q(u,v)$.
Furthermore the pushforward $(\mathbf 1\{X_T\in A\}, \mathbf 1\{Y_T\in B\})$ is distributed according to $q$ by construction.
As as long as we can find two sets $A,B$ with $\mu(A)\in (0,1)$ and $\nu(B)\in(0,1)$, we can construct infinitely many distinct couplings of $\mu$ and $\nu$ corresponding to distinct couplings of Bernoulli$(\mu(A))$ and Bernoulli$(\mu(B))$. One of these couplings is the same as that in the accepted answer, so the above construction is likely just different words for the same argument.