Ideally, the full characterization of coherent terahertz (THz) pulses would provide information on the amplitude and direction of its THz electric field, in space and in time, with unlimited dynamic range. Here, we propose and demonstrate a new approach based on the Stokes-Mueller formalism. Our approach can measure the full temporal and spatial variation of coherent THz fields, as well as its polarization state with a high dynamic range. This method employs a simple configuration, using a polarization state analyzer after the electro-optic sampling crystal. This technique could allow high sensitivity due to its ability to use thick detection crystals, which also would lead to improved spectral resolution by allowing longer scans in the time domain.