Background and ObjectiveSleep was highly imperative in human daily life. However, an increasing number of people were undergoing sleep deprivation and sleep disorders. Sleep stage classification became a highly essential process in sleep scoring. Nevertheless, the visual scoring of arousals during sleep routinely performed by sleep specialists was a challenging exercise. MethodsThis paper introduced a novel approach for sleep stage classification using covariance feature matrix architecture with multivariate phase space reconstruction (MPSR). The goal was to capture the geometric properties and the hidden dynamic characteristics of multiple physiological signals. The covariance matrices constructed through the MPSR approach were considered as symmetric and positive definite (SPD) matrices, forming a Riemannian manifold space. These SPD matrices in the Riemannian manifold were mapped to the matrices in the tangent space, allowing them to be vectorized and treated as feature vectors in Euclidean space. Finally, an ensemble learning classifier was applied to perform various sleep stage tasks. ResultsOur proposed method was evaluated on three benchmark datasets to assess its effectiveness and robustness. For the tasks of both Five class and Six class sleep stages, the proposed approach in ten-fold cross validation achieved high accuracy of 93.57% and 92.56% for the Sleep EDF dataset, 86.36% and 84.18% for the DREAMS Subjects dataset, as well as 88.93% and 88.42% for the Sleep EDF Expanded dataset (95 subjects), respectively. In leave-one-subject-out cross validation, our proposed approach for the tasks of both Five class and Six class sleep stages yielded an accuracy of 84.46% and 80.73% for the Sleep EDF dataset, 82.50% and 79.51% for the DREAMS Subjects dataset, as well as 93.25% and 92.06% for the Sleep EDF Expanded-20 dataset (20 subjects), respectively. ConclusionCompared to the traditional sample covariance matrix (SCM), the covariance feature matrix using the MPSR method successfully captured the distinction of spatial information among various sleep stages. Moreover, our proposed method obtained good performance without requiring computationally large artifact suppression or a long signal decomposition process.