Measuring the polarization of electroweak bosons at the LHC allows for important tests of the electroweak-symmetry-breaking mechanism that is realized in nature. Therefore, precise Standard Model predictions are needed for the production of polarized bosons in the presence of realistic kinematic selections. We formulate a method for the calculation of polarized cross-sections at NLO that relies on the pole approximation and the separation of polarized matrix elements at the amplitude level. In this framework, we compute NLO-accurate cross-sections for the production of two polarized Z bosons at the LHC, including for the first time NLO EW corrections and combining them with NLO QCD corrections and contributions from the gluon-induced process.