Deterministic single factor of safety method cannot explicitly account for the influences of various sources of uncertainties (e.g., spatial variability of geomaterials, measurement and transformation uncertainties) in stability design of slopes. Many probabilistic methods have been applied to the reliability‐based design (RBD) of slopes, but they typically require performing numerous deterministic slope stability analyses. In this paper, an efficient reliability‐based design method for spatially varying slopes based on field data is proposed. Here, the RBD of a slope angle is concerned. Reliability‐based design is implemented using an inverse First Order Reliability Method (inverse FORM ‐ IFORM). A sandy slope and a cohesive slope are investigated as examples, respectively, to illustrate the proposed method. The results indicate that the proposed method can quickly obtain rational design schemes of slope angles accounting for the spatial variability of soil properties, measurement and transformation uncertainties based on the field data. It can act as a practical and effective tool for the RBD of slopes in two‐dimensional spatially variable soils. Additionally, it is found that the random field mesh size affects the RBD results significantly, while the probability distribution and horizontal autocorrelation of soil parameters have slight influences on the RBD results.