The reason for the dominance of the path-integral approach presenting QED in most modern textbooks is that it is much simpler than the operator formalism, concerning how to deal with gauge invariance. Particularly when it comes to non-Abelian gauge theories, the Faddeev-Popov formalism is way more convenient than the covariant operator formalism. Of course, finally the result of both methods is the same (at least in the sense of perturbative QFT; neither the one or the other formalism is strictly formulated in a mathematically strict sense).
It's also a bit of the applications you have in mind, when presenting the material. Quantum-optics texts usually deal a lot with free photons, and there the most simple way is to fix the gauge completely using the radiation gauge for free photons, i.e., A^0=0 and \vec{\nabla} \cdot \vec{A}=0 and then quantizing the theory with only physical (spacelike transverse) field modes canonically.
For HEP applications that's a bit inconvenient, because you get Feynman rules that are not manifestly covariant. That's why one would use the Gupta-Bleuler formalism in Lorenz (covariant) gauge, but that's a bit complicated. With the path-integral formalism it's easy to derive everything in any gauge you like. Of course, at the end you get out the same physics.
In quantum optics, there's usually no reason to develop a manifestly Lorentz-covariant formalism, because it doesn't make the calculations more simple in any way. Of course, at the end there's only one successful quantum theory of the electromagnetic interaction, and that's QED. It's only used in many different fields (HEP, quantum optics, nano physics, condensed-matter physics), and in each field different choices of gauge and manifestly Lorentz invariant or non-invariant formalisms are more or less convenient. That's all.