I am self-learning elliptic functions and modular form but struggling with the very basics. I have programmed the Eisenstein series $$E_6 = \sum_{m,n \in \mathbb{Z}} \frac{1}{(m +n \tau)^6} $$ and produced the following image (color denotes the argument)
which is similar to that on wiki.
Eisenstein series and automorphic representations by P. Fleig introduced the fundamental domain, as the figure below.

A similar figure finds in L. Vepstak’s The Minkowski question Mark, GL(2,Z) and the Modular Group.

My question is how to transform from the fundamental domain $$ |Re(z)|<\frac{1}{2} \quad and \quad |z|>1| \quad and \quad Im(z)\geqslant 0 $$ to other domains (appears as distorted triangles), or how to map a point in the fundamental domain to a corresponding point in another domain?
I guess that SL(2, $\mathbb{Z}$)-invariance should associate with the symmetries between the triangular domain $$ G_{2k}\left( \frac{az+b}{cz+d} \right) = (cz+d)^{2k}G_{2k}(z) $$ but no clues how this equation divides into triangles that correspond to the fundamental domain.
To answer your question, and to understand the picture of the "fundamental domains in the upper half plane" (your third picture), it's probably best if we forget about $E_6$ and $G_{2k}$ for the time being. They are not directly relevant to your question, although we will return to them.
Instead, I would recommend that we concentrate on what the picture is actually showing.
The picture is showing the "fundamental domains" of the modular group. The modular group is a group of transformations on the (extended) upper half plane: $$ \left\{ z \mapsto\frac{az + b}{cz+d} \ : \ a, b, c, d \in \mathbb Z, \ \ ad-bc=1\right\}.$$ In the picture, the upper half plane is divided into many regions (also known as "fundamental domains"). These regions obey a nice property: Given any two regions $R_1$, $R_2$, there exists a unique transformation $g$ in the group of modular transformations that maps $R_1$ onto $R_2$.
Shall we work through an explicit example? Let's consider the transformation $$ g : z \mapsto - \frac 1 {z}$$ (i.e. $a = 0, b = 1, c = - 1, d = 0$). And let's consider the region $$ R_1 = \left\{ |{\rm Re}(z) | \leq \tfrac 1 2, \ \ |z| \geq 1 \right\}.$$ Can we work out the image of $R_1$ under the action of $g$?
Well, $g$ is a Mobius transformation, and we know that Mobius transformations map (generalised) circles to (generalised) circles. The boundary of $R_1$ is made up of segments of generalised circles. So the boundary of $g(R_1)$ should be made up of segments of generalised circles too. To be more specific:
The part of the boundary of $R_1$ that is a line segment containing the points $\frac 1 2 + \frac{i \sqrt 3}{2}$, $\frac 1 2 + i$ and $\infty$ is mapped to a circle segment containing the points $-\frac 1 2 + \frac{i\sqrt 3}{2}$, $- \frac 2 5 + \frac{4i}{5}$ and $0$.
The part of the boundary of $R_1$ that is a line segment containing the points $-\frac 1 2 + \frac{i \sqrt 3}{2}$, $-\frac 1 2 + i$ and $\infty$ is mapped to a circle segment containing the points $\frac 1 2 + \frac{i\sqrt 3}{2}$, $ \frac 2 5 + \frac{4i}{5}$ and $0$.
The part of the boundary of $R_1$ that is a circle segment containing the points $- \frac 1 2 + \frac{i\sqrt 3}{2}$ and $i$, $-\frac 1 2 + \frac{i\sqrt 3}{2}$ is mapped to itself.
Thus $g$ maps our region $R_1$ to the region in your diagram with vertices at $0$, $- \frac 1 2 + \frac{i \sqrt 3}{2}$ and $\frac 1 2 + \frac{i \sqrt 3}{2}$. And for the remainder of this answer, I'll call this new region $R_2$.
So what does all this have to do with your Eisenstein series? Well, the Eisenstein series obey the relation $$ G_{2k} \left( \frac{az+b}{cz+d}\right) = (cz+d)^{2k}G_{2k}(z).$$ You may like to think of this equation as a simple relationship between the values of the function $G_{2k}$ in one region of your picture and the values in another region.
For example, we've seen that $g : z \mapsto - 1 / z$ maps the region $R_1$ to the region $R_2$. The equation $$ G_{2k} \left( - \frac 1 z \right) = z^{2k} G_{2k}(z)$$ is then a nice relationship between the values that $G_{2k}(z)$ takes in the region $R_1$ and the values it takes in the region $R_2$: if you calculate the values of $G_{2k}$ in one region, you get the values of $G_{2k}$ in the other region for free.