As for the Heisenberg picture... it's still there! You can switch to a Heisenberg picture through a unitary transformation. In this picture the operators evolve through time as:
A_H(t) = e^{iHt} A_H(0) e^{-iHt}
The operator e^{-iHt) is the time-evolution operator. In the Schroedinger picture the states transform in the usual way:
|\psi_S(t)\rangle = e^{-iHt}|\psi_S(t)\rangle
And here is something you do not hear every day (but is certainly still true!). The equation of motion of the time evolution operator is nothing but Schroedingers equation!
\frac{d}{dt} e^{-iHt} = -i H e^{-iHt}
Yes, the Schroedinger equation still plays a role in relativistic quantum field theory!
In fact, for fields operators we can apply the time evolution for operators. Take for instance the conjugate field \pi(x,t)
\pi(x,t) = e^{iHt} \pi(x,0) e^{-iHt}
which -- when taking the time derivative -- reduces to Heisenbergs equations of motion:
\partial_t \phi(x,t) = i[H,\pi(x,t)]
But for a free scalar field theory, i.e. KG theory, this commutator can be computed explicitly. Namely, the Hamiltonian is given by:
H = \frac{1}{2} \int d^3 x' \left(\pi(x,0)^2 + (\nabla\phi(x,t))^2 +m^2\phi(x,t)^2\right)
Now, the equal-time commutation relations should be familiar
[\phi(x,t),\pi(x',t)] = i \delta(x-x')
with on the right hand side the three dimensional spatial delta function. All other commutators are zero. The commutator we are after is then...
[H,\pi(x,t)] = i\int d^3x' \left(\nabla^2 \phi(x',t) + m^2\phi(x',t)\right)\delta(x-x') = \nabla^2 \phi(x',t) + m^2\phi(x',t)
You need to apply some partial integration to obtain this, but it's quite straightforward. Anyways, now, this commutator determines the time-evolution of the conjugate field. So we have:
\partial_t \pi(x,t) = \nabla^2 \phi(x,t) - m^2\phi(x,t)
This is simply the Heisenberg equation of motion for the conjugate field. But remember, the conjugate field \pi(x,t) is nothing but the time derivative of the original KG field \phi(x,t). And so we end up with...
\partial_t^2 \phi(x,t) - \nabla^2 \phi(x,t) + m^2\phi(x,t) = 0
Ha, that should look familiar!