Due to providing greater degrees-of-freedom (DOFs) and securer communication guarantees than traditional scalar array, polarized massive multiple-input multiple-output (MIMO) technique offers a prospective insight into millimeter-wave (mmWave) communication. In this paper, a compressed two-dimensional (2D) direction-of-arrival (DOA) and polarization estimation algorithm combining the compressive sampling propagator method with reduced-dimension multiple signal classification (CSPM-RDMUSIC) is developed for mmWave polarized massive MIMO systems. First, a synchronous compressive network is constructed to compress the high-dimensional output into a low-dimensional one. Thereafter, a propagator is structured to acquire the signal subspace. Subsequently, accurate 2D-DOA and polarization estimation are achieved using a coarse-refined strategy, where the vector cross-product technology is investigated for coarse 2D-DOA estimation and then RDMUSIC is exploited to perform 2D local search for accurate parameter estimation. Owing to integrating the CSPM framework with the RDMUSIC strategy, the developed approach provides high-precision parameter estimation with relative computational economy, supporting mmWave communication demands. What's more, it's suitable for arbitrary array geometry with high flexibility. Several experimental results validate the superiorities of the developed framework.