The Galilean-invariant field theories are quantized by using the canonical method and the five-dimensional Lorentz-like covariant expressions of non-relativistic field equations. This method is motivated by the fact that the extended Galilei group in 3 + 1 dimensions is a subgroup of the inhomogeneous Lorentz group in 4 + 1 dimensions. First, we consider complex scalar fields, where the Schrödinger field follows from a reduction of the Klein–Gordon equation in the extended space. The underlying discrete symmetries are discussed, and we calculate the scattering cross-sections for the Coulomb interaction and for the self-interacting term λΦ 4. Then, we turn to the Dirac equation, which, upon dimensional reduction, leads to the Lévy-Leblond equations. Like its relativistic analogue, the model allows for the existence of antiparticles. Scattering amplitudes and cross-sections are calculated for the Coulomb interaction, the electron–electron and the electron–positron scattering. These examples show that the so-called ‘non-relativistic’ approximations, obtained in low-velocity limits, must be treated with great care to be Galilei-invariant. The non-relativistic Proca field is discussed briefly.