Here's a general solution to the problem of defining the representations of gamma matrices.
A representation of gamma matrices is completely defined by defining the following five state vectors:
$$
D_0 = D_{++} = \left[\begin{matrix} 1 \\ 0 \\ 0 \\ 0 \end{matrix}\right],\;\;\;
D_1 = D_{+-} = \left[\begin{matrix} 0 \\ 1 \\ 0 \\ 0 \end{matrix}\right],\;\;\;
D_2 = D_{-+} = \left[\begin{matrix} 0 \\ 0 \\ 1 \\ 0 \end{matrix}\right],\\
D_3 = D_{--} = \left[\begin{matrix} 0 \\ 0 \\ 0 \\ 1 \end{matrix}\right],\;\;\;\;\;
V = \left[\begin{matrix} 1/2 \\ 1/2 \\ 1/2 \\ 1/2 \end{matrix}\right].\\
$$
The reason for calling these the "D" and "V" states is that the density matrices of the first four are diagonal while the last one is an example of what Schwinger called the "vacuum" in one of his papers on the Schwinger Measurement Algebra. I've given two notations for the D_xx in order to simplify the derivation.
The four diagonal states are defined by a "Complete Set of Commuting Observables". This is just enough information to determine a quantum state. For the case of Dirac gamma matrices, a complete set of commuting observables will have two (distinct) observables that happen to be "roots of unity". That is, the two observables $P$ and $Q$ have to satisfy:
$$PP = QQ = 1,\;\;P\;\textrm{not equal to}\;Q ,\;\;\textrm{and}\;\; PQ = QP,$$
with neither equal to 1 or -1. An example of two gamma algebra elements that satisfy these requirements are
$$P = \gamma^3\gamma^0,\;\;Q = i\gamma^1\gamma^2$$
Ignoring signs, there are sixteen possible products of gamma matrices:
$$\{1,\gamma^1,\gamma^2,\gamma^1\gamma^2, ... \gamma^1\gamma^2\gamma^3\gamma^0\}$$
To use one of these as a root of unity we have to make it square to +1. So if it squares instead to -1, simply multiply it by i. We can't use 1 as $P$, so after converting the above to roots of unity, we have 15 choices for $P$.
Once we choose $P$, there are 14 roots of unity left. Of these, just six will commute with $Q$ so we can choose one of those six. And having chosen $P$ and $Q$, their product $PQ=QP$ will also be a root of unity that commutes with $P$ and $Q$. So in addition to $\{P,Q\}$, we could as well have chosen $\{P,PQ\}$ or $\{Q,PQ\}$ and ended up with the same representation.
Next we have to choose $V$. This is a single state that is a "complementary state" to our Complete Set of Commuting Observables. This is a basic concept in quantum mechanics; the usual example is position and momentum. The requirement is that the transition probabilities between the state $V$ and the $D_n$ states all need to be the same.
In fact, the $V$ state is easy to choose. Of the 16 products of gamma matrices, one was trivial (i.e. 1), and three are the $\{P,Q,PQ\}$. The remaining 12 products are complementary observables. All we need to do is to choose two of them that are distinct and that commute. Call them $R$ and $S$. We can choose $R$ completely freely from those 12 products. After choosing $R$, there will be four choices for $S$ and, as before, they will come in pairs so there are really only two choices. Then the $V$ state will take +1 or -1 as its eigenvalues with respect to $R$ and $S$ so with $R$ and $S$ in hand, there will be four choices for $V$.
Things are easier if we write them in pure density matrix form. If we use the +1 eigenvalues for the state $V$, its pure density matrix state will be:
$$\rho_V = (1+R)(1+S)/4
=\frac{1}{4}\left[\begin{matrix}
1&1&1&1 \\ 1&1&1&1 \\ 1&1&1&1 \\ 1&1&1&1 \end{matrix}\right]$$
Now we can write the four (pure) density matrices corresponding to $D_n$ as:
$$
\rho_{D++} = (1+P)(1+Q)/4\\
\rho_{D+-} = (1+P)(1-Q)/4\\
\rho_{D-+} = (1-P)(1+Q)/4\\
\rho_{D--} = (1-P)(1-Q)/4\\$$
When written in the representation, the above are to be diagonal matrices, each with a single 1 on the diagonal and the rest zero.
Similarly, we can pick out the nm matrix element of the representation by multiplications of density matrices of the form:
$$M_{nm} = \rho_{Dn}\;\rho_V\;\rho_{Dm}$$
For example:
$$M_{02} = \rho_{D0}\;\rho_V\;\rho_{D2}$$
which in the representation is explicitly:
$$\frac{1}{4}\left[\begin{matrix}0&0&1&0\\0&0&0&0\\0&0&0&0\\0&0&0&0\end{matrix}\right]=
\frac{1}{4}\left[\begin{matrix}1&0&0&0\\0&0&0&0\\0&0&0&0\\0&0&0&0\end{matrix}\right]
\left[\begin{matrix}1&1&1&1\\1&1&1&1\\1&1&1&1\\1&1&1&1\end{matrix}\right]
\left[\begin{matrix}0&0&0&0\\0&0&0&0\\0&0&1&0\\0&0&0&0\end{matrix}\right]$$
Since we now know which algebra elements represent a particular spot in the matrix, we can read off the representation by seeing which algebra elements contribute to which matrix elements. For example, if $\gamma_0$ shows up in $M_{02}$, then the representation of $\gamma_0$ will have a 1 in its (0,2) position.
As it turns out, I've been doing research on gamma matrices recently and have written an application using Java that allows the user to work through the above steps by clicking buttons in a graphical user interface. It shows resulting representation and you can also show the results for other products of gamma matrices such as $i\gamma^1\gamma^2$. I put in some presets so you can get standard gamma matrix representations in a single key click and found this question while using a search engine to find more cases. I'll link it in when it's done.
Also, it should be noted that the above explicitly describes only the representations that happen to represent the products of gamma matrices as matrices whose 16 elements include four complex phases and 12 zeroes. By rotations and boosts, you can get more general representations.