Classification techniques play an important role in the analysis of polarimetric synthetic aperture radar (PolSAR) images. PolSAR image classification is widely used in the fields of information extraction and scene interpretation or is performed as a preprocessing step for further applications. However, inherent speckle noise of PolSAR images hinders its application on further classification. A novel supervised superpixel-based classification method is proposed in this study to suppress the influence of speckle noise on PolSAR images for the purpose of obtaining accurate and consistent classification results. This method combines statistical information with spatial context information based on the stochastic expectation maximization (SEM) algorithm. First, a modified simple linear iterative clustering (SLIC) algorithm is utilized to generate superpixels as classification elements. Second, class posterior probabilities of superpixels are calculated by a K distribution in iterations of SEM. Then, a neighborhood function is defined to express the spatial relationship among adjacent superpixels quantitatively, and the class posterior probabilities are updated by this predefined neighborhood function in a probabilistic label relaxation (PLR) procedure. The final classification result is obtained by the maximum a posteriori decision rule. A simulated image, a spaceborne RADARSAT-2 image, and an airborne AIRSAR image are used to evaluate the proposed method, and the classification accuracy of our proposed method is 99.28%, 93.16% and 89.70%, respectively. The experimental results indicate that the proposed method obtains more accurate and consistent results than other methods.