The following holds for any nonzero Gaussian integer $a+bi$.
If $d=gcd(a,b)$, then $\mathbb{Z}[i]/(a+bi)$ is isomorphic (as a group) to $$\mathbb{Z}/d\mathbb{Z} \oplus \mathbb{Z}/ f\mathbb{Z},$$
where $f=\frac{a^2+b^2}{d}=aa'+bb',$ and $a=a'd, b=b'd.$
There might be a more explicit link, but here one sees how the quotient group is classified through the use of matrices in $GL_n(\mathbb{Z})$.
http://www.cs.uleth.ca/~holzmann/notes/abelian.pdf
Take the vectors $u_1=(a,b), u_2=(-b,a)$ that form a basis of the lattice $(a+bi)$. The Bézout identity reads
$$ma'+nb'=1,$$
and suggests a $2\times 2$ invertible matrix (over $\mathbb{Z}$); a new basis of our lattice is now given by
$$v_1=b' u_1+a'nu_2=(0,aa'+bb')=(0,f),$$ and the vector $v_2= m u_1-n u_2=(d, d(b'm-a'n)).$ Changing now the basis in $\mathbb{Z}\oplus\mathbb{Z}$ to $a_1=(1,0), a_2=(1,s)$ where $s=b'm-a'n$ one has:
$$(a+bi)= \mathbb{Z} fi \oplus \mathbb{Z} d.(1+si) \subset \mathbb{Z}i \oplus \mathbb{Z} (1+si)=\mathbb{Z}[i],$$ which yields the desired result.
Note that the cardinal of the group $=df$ equals the determinant of the initial matrix, which is $a^2+b^2$.
\paragraph{DECOMPOSITION AS A RING}
In order to study the decomposition as a ring, one need factor the number $$a+bi=\prod \pi_i^{e_i},$$ then apply the Chinese number theorem.
Primes in this ring are $1+i, p$ (where $p$ is an odd prime which is $3$ modulo $4$), and prime factors of prime integers $q$ which are $\equiv 1 (\mbox{ mod } 4)$; this last type are of the form $m+ni$ where $m^2+n^2=q.$ This is done as follows, for any prime $p$:
$$\mathbb{Z}[i]/(p)=\mathbb{Z}[x]/(p,x^2+1)=\mathbb{F}_p[x]/(x^2+1),$$
and this ring is nonreduced when $p=2.$ If $p$ is odd, the ring is a domain (hence a field) exactly when $x^2+1$ has no zeros in $\mathbb{F}_p$, which is precisely when $p\equiv 1 (\mbox{ mod }4)$ (look up the Legendre symbol).
Finally, if $4|p-1$ then the above ring has two different primes, which translates into saying that $(p)=(\pi_1)\cap (\pi_2)$. After working a bit more, one derives $$p=\pi_1 \pi_2$ and $N(\pi_i)=p$, which yields $\pi_1=\overline{\pi_2}.$ The result follows.
A more detailed discussion would considerably extend this paragraph.