Raise of Slopes — User Guide

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.

Quick-start workflow

  1. Load a DEM raster into QGIS.
  2. Open the plugin and go to : select two points on the map to extract the cross-section.
  3. Go to and enter the geotechnical properties.
  4. Optionally configure a second layer and / or water table in .
  5. Run an analysis from or using the Run button at the bottom of the window.
  6. Repeat for multiple cross-sections and accumulate results into a .
  7. Export results (chart, CSV, DXF, text report) from the Export menu.

Acknowledgements

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.

How to cite

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

Tab: Elevation Profile

All analyses are anchored to the elevation profile. Compute this first.

Profile extraction

Chart controls

Export menu

OptionOutput
Export CSVDistance–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 DXFProfile polyline and slip surfaces on separate DXF layers, with a legend. Uses ezdxf if installed, otherwise writes a minimal ASCII DXF.

Project menu

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.

Tab: Soil Parameters

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.

Tab: Stratigraphy & Water Table

Second layer

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.

Water table

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.

💡 Preview the stratigraphy and water table directly on the profile chart as you adjust the parameters — the chart updates in real time.

Tab: Grid Analysis

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.

ParameterDescription
Calculation methodLEM solver: Bishop (simplified), Morgenstern–Price, or Spencer.
Entry pointsNumber of grid nodes along the entry (upslope) search interval.
Exit pointsNumber of grid nodes along the exit (downslope) search interval.
Min η incrementMinimum 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 displayNumber of surfaces with the lowest FoS to show, coloured on a red–green scale.
💡 For a typical downslope analysis set In min=0.0, In max=0.5 and Out min=0.5, Out max=1.0 so that the entry and exit search intervals do not overlap. Start with 10×10 grid points and refine once you know the critical zone.

Tab: Simplex Analysis

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.

ParameterDescription
Calculation methodBishop, 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 / maxAngular bounds for the arc curvature parameter (degrees).
Max iterationsIteration cap for the Nelder–Mead optimiser (default 300).
Surfaces to displayNumber of distinct optimised surfaces to show (deduplicated by convergence point).

Tab: Hazard Map

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.

Workflow for a multi-profile hazard map

  1. Analyse the first cross-section and switch to Hazard Map.
  2. Select Create new rasters, define the output extent, cell size, and output paths, then click Add current result to map.
  3. Enable Update rasters automatically after each analysis for subsequent profiles.
  4. Analyse each additional cross-section; the rasters update automatically.
  5. To restart from scratch, switch to Update existing rasters and enable Overwrite.
⚠ The Hazard Map requires GDAL (osgeo) to be available in the QGIS Python environment. This is included in standard QGIS installations.

Stability methods

MethodInterslice forcesEquilibrium satisfiedNotes
Bishop (simplified)Horizontal onlyMomentFast and accurate for circular surfaces. Recommended default.
Morgenstern–PriceInclined (constant f(x))Force + MomentRigorous. Suitable for irregular or non-circular surfaces.
SpencerInclined (constant angle)Force + MomentRigorous. Similar accuracy to Morgenstern–Price.

Results interpretation

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.

File formats

ExtensionContents
.rslopeFull project state (JSON): profile data, parameters, computed surfaces, hazard-map configuration.
.csvElevation profile: two columns, distance (m) and elevation (m).
.txtPlain-text analysis report including parameters, results, and profile CSV.
.png / .svgProfile chart image (raster or vector).
.dxfProfile and slip surfaces as DXF polylines on separate named layers, with a legend.