A continuum model for a population of self-propelled particles
interacting through nematic alignment is derived from an individual-based model. The methodology consists of introducing a hydrodynamic scaling of the corresponding mean field kinetic equation. The resulting perturbation problem is solved thanks to the concept of generalized collision invariants. It yields a hyperbolic but non-conservative system of equations for the nematic
mean direction of the ow and the densities of particles owing parallel or anti-parallel to this mean direction. Diffusive terms are introduced under a weakly non-local interaction assumption and the diffusion coefficient is proven to be positive. An application to the modeling of myxobacteria is outlined.