We assess the ability to image the photosphere of red supergiants and, in particular Betelgeuse, through the modelling of the observed linear polarization in atomic spectral lines. We also aim to analyse the resulting images over time, to measure the size and dynamics of the convective structures in these stars. Rayleigh scattering polarizes the continuum and spectral lines depolarize it.This depolarization is seen as a linear polarization signal parallel to the radial direction on the stellar disk. Integrated over the disk, it would result in a null signal, except if brightness asymmetries/inhomogeneities are present. This is the basic concept behind our imaging technique. The several tests and comparisons performed prove that our technique reliably retrieves the salient brightness structures in the photosphere of Betelgeuse, and should be relevant to other red supergiants. We demonstrate that these structures are convective cells, with a characteristic size of more than 60% of the stellar radius. We also derive the characteristic upflow and downflow speeds, 22 and 10 km/s respectively. We find weak magnetic fields concentrated in the downflow lanes in between granules. We follow those convective structures in time. Changes happen on timescales of one week, but individual structures can be tracked over four years of observations.