This QGIS plugin performs slope stability analysis based on the Limit Equilibrium Method (Bishop, Spencer, Morgenstern-Price) with circular slip surfaces. Starting from a DEM, it extracts an elevation cross-section, models the subsurface stratigraphy and phreatic surface, and finds the critical slip surface that minimises the Factor of Safety (FoS) using either an exhaustive grid search or a Nelder–Mead simplex optimiser.
This study was funded by the European Union — NextGenerationEU and by the Ministry of University and Research (MUR), National Recovery and Resilience Plan (NRRP), Mission 4, Component 2, Investment 1.5, project “RAISE — Robotics and AI for Socio-economic Empowerment” (ECS00000035). Andrea Bressan, Simone Pittaluga, Lorenzo Tamellini and Domenico Gallipoli are part of the RAISE Innovation Ecosystem.
If you use this plugin in your research or professional work, please cite the underlying methodology:
Lalicata, L.M., Bressan, A., Pittaluga, S. et al. An Efficient Slope Stability Algorithm with Physically Consistent Parametrisation of Slip Surfaces. Int J Civ Eng 23, 671–682 (2025).
https://doi.org/10.1007/s40999-024-01053-1
All analyses are anchored to the elevation profile. Compute this first.
| Option | Output |
|---|---|
| Export CSV | Distance–elevation table of the sampled profile. |
| Export results (TXT) | Plain-text report with all analysis results, parameters, and surface summary. |
| Export image (PNG/SVG) | The profile chart as a raster or vector image. Optionally includes the stratigraphy legend. |
| Export DXF | Profile polyline and slip surfaces on separate DXF layers, with a legend. Uses ezdxf if installed, otherwise writes a minimal ASCII DXF. |
Save or load the complete plugin state to a .rslope file (JSON). The file stores the sampled profile data, all parameters, computed surfaces, surface visibility, and hazard-map configuration. On load, any raster layers referenced by the project are looked up in the current QGIS project by name; missing layers are reported with an option to continue.
Geotechnical properties for the primary (surface) soil layer. These values are shared between the Grid and Simplex analyses.
The Second Layer parameter block (γ₂, c₂, n₂, φ₂) is used when a second stratum is enabled in the Stratigraphy & Water Table tab.
Enables a two-layer model. The interface between layer 1 (surface soil) and layer 2 (deeper formation) can be defined in three ways:
The interface is drawn as a dashed brown line on the profile chart when the second layer is active.
Enables pore-pressure calculations. The water table can be defined using the same three modes as the second layer. When active, the phreatic surface is drawn as a blue dashed line on the chart, and pore pressures modify the effective stresses in the LEM solver.
Evaluates the Factor of Safety for a large discrete set of circular slip surfaces sampled on a regular grid of (entry, exit) pairs. Provides a global view of the stability landscape at the cost of computation time.
| Parameter | Description |
|---|---|
| Calculation method | LEM solver: Bishop (simplified), Morgenstern–Price, or Spencer. |
| Entry points | Number of grid nodes along the entry (upslope) search interval. |
| Exit points | Number of grid nodes along the exit (downslope) search interval. |
| Min η increment | Minimum angular step for the slip-circle curvature parameter η. Smaller values explore shallower arcs. |
| In – min / max (fraction) | Fractional position along the profile [0 = P1, 1 = P2] where slip surfaces may enter the slope. |
| Out – min / max (fraction) | Fractional position where slip surfaces may exit the slope. |
| Surfaces to display | Number of surfaces with the lowest FoS to show, coloured on a red–green scale. |
Uses a Nelder–Mead optimiser seeded from a coarse grid to find the minimum-FoS surface. Much faster than a full grid for the same resolution, and well-suited for fine-tuning around a known critical zone.
| Parameter | Description |
|---|---|
| Calculation method | Bishop, Morgenstern–Price, or Spencer. |
| x_in min / max (fraction) | Search bounds for the entry point, as a fraction of the profile length. |
| x_out min / max (fraction) | Search bounds for the exit point. |
| η min / max | Angular bounds for the arc curvature parameter (degrees). |
| Max iterations | Iteration cap for the Nelder–Mead optimiser (default 300). |
| Surfaces to display | Number of distinct optimised surfaces to show (deduplicated by convergence point). |
Accumulates stability results from many profile analyses into two georeferenced GeoTIFF rasters covering the study area:
The rasters are styled automatically when loaded into QGIS: FoS uses the red–yellow–green colour ramp; depth uses a white–blue ramp.
osgeo) to be available in the QGIS Python environment. This is included in standard QGIS installations.| Method | Interslice forces | Equilibrium satisfied | Notes |
|---|---|---|---|
| Bishop (simplified) | Horizontal only | Moment | Fast and accurate for circular surfaces. Recommended default. |
| Morgenstern–Price | Inclined (constant f(x)) | Force + Moment | Rigorous. Suitable for irregular or non-circular surfaces. |
| Spencer | Inclined (constant angle) | Force + Moment | Rigorous. Similar accuracy to Morgenstern–Price. |
On the chart, slip surfaces are coloured on a continuous red–yellow–green scale spanning the FoS range of all currently displayed surfaces. The critical surface (lowest FoS) is drawn with a thicker line and labelled #1.
| Extension | Contents |
|---|---|
.rslope | Full project state (JSON): profile data, parameters, computed surfaces, hazard-map configuration. |
.csv | Elevation profile: two columns, distance (m) and elevation (m). |
.txt | Plain-text analysis report including parameters, results, and profile CSV. |
.png / .svg | Profile chart image (raster or vector). |
.dxf | Profile and slip surfaces as DXF polylines on separate named layers, with a legend. |