As you have said we have to perform a translation and a rotation.
For the translation note that the plane $ax+by+cz+d=0$ intersects the $z$ axis at $(0,0,-d/c)$ so the translation $\vec{t}:(x,y,z)\rightarrow (x,y,z-d/c)$ give you a plane $ax+by+cz=0$ that passes thorough the origin and is orthogonal to the vector $\vec{v}=(a,b,c)^T$
For the rotation note that the angle between $\vec{v}$ and $\vec{k}=(0,0,1)^T$ is given by:
$$
\cos \theta=\dfrac{(\vec{v},\vec{k})}{|\vec{v}|}=\dfrac{c}{\sqrt{a^2+b^2+c^2}}
$$
and the axis of rotation have to be orthogonal to $\vec{v}$ and $\vec{k}$ so its versor is:
$$
\vec{u}=\dfrac{\vec{v}\times\vec{k}}{|\vec{v}\times\vec{k}|}=\dfrac{1}{\sqrt{a^2+b^2}}\left(b,-a,0\right)^T=(u_1,u_2,0)^T
$$
the rotation (if I've not made some typo) is represented by the matrix:
$$
\left (
\begin{array}{cccc}
\cos \theta +u_1^2 (1-\cos \theta) &u_1u_2 (1-\cos \theta) & +u_2\sin \theta \\
u_1u_2 (1-\cos \theta)& \cos \theta+ u_2^2 (1-\cos \theta)& -u_1 \sin \theta \\
-u_2 \sin \theta & u_1 \sin \theta& \cos \theta
\end {array}
\right)
$$
see here.
EDIT:
To help you compute the rotation matrix more quickly, I've isolated all the quantities that appear:
$ cos \theta = \frac{c}{\sqrt{a^2 + b^2 + c^2}} $
$ sin \theta = \sqrt{\frac{a^2+b^2}{a^2 + b^2 + c^2}} $
$ u_1 = \frac{b}{\sqrt{a^2 + b^2 }} $
$ u_2 = -\frac{a}{\sqrt{a^2 + b^2}} $