Why does the ideal $(a+bi)$ have index $a^2+b^2$ in $\mathbb{Z}[i]$? [duplicate]

In the comments on the question Why does this module structure have $352512$ elements?, it is mentioned that the index of the ideal generated by $a+bi$ in $\mathbb{Z}[i]$ has order $a^2+b^2$.

Is there a nice rigorous explanation of why this is so?

In Quotient ring of Gaussian integers it was shown that for $a,b$ coprime the quotient is actually isomorphic to $\mathbb Z / (a^2+b^2)$, but in general, it is still true that they have the same size.

If you look at the picture

taken from the answer Quotient ring of Gaussian integers by quanta in the above-mentioned thread, you see that we want to count the lattice points in the square spanned by $a+bi$ and $-b+ai$.

As the points on the border have to be partially identified, it turns out we want to count interior points plus half the border points -1. (Since opposite sides of the square are identified, we want to count only half the border points, but we want to count only 1 of the 4 corners, so we have to subtract one.)

This gives exactly the area $\left(\sqrt{a^2+b^2}\right)^2=a^2+b^2$ of the square by Pick's theorem.

One have $(a+ib)\mathbb Z[i]= (a+ib)\mathbb Z \oplus (-b+ia)\mathbb Z$, so the ideal $(a+ib)$ is the free $\mathbb Z$-submodule of $\mathbb Z[i] =\mathbb Z \oplus i\mathbb Z$ generated by $a+ib$ and $-b+ia$. Some basic algebra theory tells you that the index of this submodule is $\det \begin{pmatrix} a&-b\\b&a \end{pmatrix} = a^2+b^2$.

Edit : See here for the proof : Why is the determinant equal to the index?