In this paper, a numerical procedure of probabilistic slope stability analysis that considers the spatial variability of soil properties is presented. The procedure extends the deterministic analysis based on the limit equilibrium method of slices to a probabilistic approach that accounts for the uncertainties and spatial variation of the soil parameters. Making no a priori assumptions about the critical failure surface like the Random Finite Element Method (RFEM), the approach saves the amount of solution time required to perform the analysis. Two-dimensional random fields are generated based on a Karhunen-Lo$grave{e}$ve expansion in a fashion consistent with a specified marginal distribution function and an autocorrelation function. A Monte Carlo simulation is then used to determine the statistical response based on the random fields. A series of analyses were performed to verify the application potential of the proposed method and to study the effects of uncertainty caused by the spatial heterogeneity on the stability of slope. The results show that the proposed method can efficiently consider the various failure mechanisms caused by the spatial variability of soil property in the probabilistic slope stability assessment.