Every elliptic surface on $\mathbb P^1$ with four semistable fibers and no other singular fibers has a $j$ invariant with ramification of order $3$ around $j=0$ and ramification of order $2$ around $j=1728$. By Riemann-Roch, one can check that this is the only ramification, so this means the $j$ invariant is a Belyi map coming from a clean dessin with all vertices of valence $3$, with four points in the fiber over $\infty$.

Conversely, from such a Belyi map one can construct an elliptic surface, uniquely defined up to quadratic twist. One can choose the quadratic twist so that all fibers with $j\neq \infty$ are smooth and all fibers with $j=\infty$ are smooth.

Beauville (http://math1.unice.fr/~beauvill/pubs/ellss.pdf) classified such surfaces by classifying their Belyi map. (The diagram on page 659 is a combinatorial way of specifying a Belyi map - though I think a less visually attractive way than Grothendick's.) He found explicit equations for the surfaces, from which you can extract the Belyi map by taking the j invariant.

In particular your map, as a function of $t$, is the $j$ invariant of the curve with equation $$(X+Y)(Y+Z)(Z+X)+tXYZ=0$$ because according to his table this is the unique one with ramification at $\infty$ of order $1,2,3,6$. Computing the $j$ invariant from this is routine, at least in theory.

He also related them to modular curves. For instance, your map can be described as the covering map from the modular curve $X_1(6)$ to $X(6)$. One can check directly from the group theory of $SL_2(\mathbb Z)$ that this modular map has the desired ramification properties, and use the uniqueness of the dessins to conclude it equals yours.