Blog

Landslide probabilities

Short title: landprob

Inputs: DEM in metric projection, Source Landslide areas, H/L value as stopping parameter.
Outputs: pq_limi and all other parameters as described in Gruber et al., 2008 article in the Geomorphometry book.

Purpose and use:

Computes Landslide probabilities — generate landslide probabilities areas using a DEM according to Gruber et al 2008. &r landslide1 <sourcearea\_dem> {h/l stopping value}

Programming environment: Arc AML
Status of work: Public Domain
Reference: Geomorphometry: Concepts, Software, Applications Data set name: Baranja hill

Attachment:

landslide_probabilities.zip

Flow magnitude and drainage basins

Short title:flowMag.lsc

Inputs:A DEM (the Baranja Hill 25m DEM is used as an example).
Outputs:A raster containing flow magnitude values for the DEM and a raster containing drainage basins.

Purpose and use:

Script to calculate flow magnitude and drainage basins from a DEM. It demonstrates how recursive function calls in LandScript can be used to calculate zonal geomorphometric parameters using map algebra. To use this script, start LandSerf 2.3 or above and then open the LandScript editor (menu: Edit->LandSctipt Editor).
Part of this script appears in Chapter 14 of Hengl and Reuter (2008) [the Geomorphometry book]. This script runs slowly so is better as a demonstration of how recursive function calls can be made in LandScript rather than a production quality flow magnitude calculator.

Programming environment:Landserf
Status of work:Public Domain
Reference: Geomorphometry: Concepts, Software, Applications
Data set name: Baranja hill

Attachment:

flowMag.lsc_.zip

Convert to Landserf format

Short title:importData.lsc

Inputs:Baranja Hill data (DEM25m.asc, DEM25srtm.asc, orthophoto.asc).
Outputs:LandSerf files of the input files plus a difference map of the two DEMs.

Purpose and use:

LandScript to import Baranja Hill data and convert to LandSerf format. This script is also included in Chapter 14 of Hengl and Reuter (2008) [the Geomorphometry book]

Programming environment:Landserf
Status of work:Public Domain
Reference: Geomorphometry: Concepts, Software, Applications
Data set name: Baranja hill


Attachment:

importData.lsc_.zip

Characteristic scale

Short title: characteristicScale.lsc

Inputs: A DEM (the Baranja Hill 25m DEM is used as an example). The name of the parameter to measure and the minimum and maximum window sizes over which to measure it.
Outputs: Two rasters, one containing the measured parameter, the other the window size at which the parameter is most extreme.

Purpose and use:

Finds the scale at which a geomorphometric parameter is most extreme for each cell in a DEM. Part of this script appears in Hengl and Reuter (2008) [the Geomorphometry book].
Script to measure surface parameter at characteristic scales. It is designed to incorporate scale-based analysis into surface parameterisation. It measures the given parameter (e.g. slope, profile curvature etc.) at a range of scales, and finds the scale at which that parameter is most extreme. Can be used to explore scale sensitivity of a surface.

Programming environment: Landserf
Status of work: Public Domain
Reference: Geomorphometry: Concepts, Software, Applications
Data set name: Baranja hill


Attachment:

characteristicScale.lsc_.zip

Volcano Maungawhau

Maunga Whau (Mt Eden) is one of about 50 volcanos in the Auckland volcanic field. This data set gives topographic information for Maunga Whau on a 10 m by 10 m grid. A matrix with 87 rows and 61 columns, rows corresponding to grid lines running east to west and columns to grid lines running south to north.

Perspective view - volcano

Available layers:

- volcano_maungawhau.asc — the 10 m DEM digitized from the topo map;

Grid definition:

ncols: 61
nrows: 87
xllcorner: 2667400
yllcorner: 6478700
cellsize: 10 m

proj4:+init=epsg:27200

Lineage:

Digitized from a topographic map by Ross Ihaka. These data should not be regarded as accurate.

\> data(volcano)
>library(spatstat)
>LLC <- data.frame(E=174.761345, N=-36.879784)
>coordinates(LLC) <- ~E+N
>proj4string(LLC) <- CRS("+proj=longlat +datum=WGS84")
>LLC.NZGD49 <- spTransform(LLC, CRS("+init=epsg:27200"))
>volcano.r <- as.im(list(x=seq(from=2667405, length.out=61, by=10),
+     y=seq(from=6478705, length.out=87, by=10), 
+     z=t(volcano)\[61:1,\]))
>volcano.sp <- as(volcano.r, "SpatialGridDataFrame")
>proj4string(volcano.sp) <- CRS("+init=epsg:27200")
# str(volcano.sp)
# spplot(volcano.sp, at=seq(min(volcano.sp$v), max(volcano.sp$v),5),
+    col.regions=topo.colors(45))
>write.asciigrid(volcano.sp, "volcano\_maungawhau.asc", na.value=-1)

Data owner: LINZ
Location:Volcano Maungawhau, Auckland, New Zealand
36° 50’ 50.586” S,174° 45’ 56.646” E
See map: Google Maps

Attachment:
volcano_maungawhau.zip

Latest Posts

Whitebox Workflows Next Gen released

Whitebox Workflows Next Gen, a complete rewrite of the Whitebox geospatial analysis platform, is now publicly available for Python, R, and QGIS.

Built from the ground up in pure Rust, Whitebox Next Gen provides more than 700 geospatial analysis tools for geomorphometry, terrain analysis, spatial hydrology, LiDAR processing, remote sensing, vector GIS, and spatial statistics. All three interfaces run on the same high-performance backend, allowing users to work in scripts, notebooks, statistical workflows, or a familiar desktop GIS environment.

Whitebox Next Gen is not an incremental update to the previous Whitebox architecture. It replaces the earlier monolithic system with a modular, full-stack geospatial platform. Core capabilities for raster and vector I/O, coordinate systems and reprojection, vector topology, spatial indexing, LiDAR processing, and other foundational operations are implemented directly within the Whitebox codebase rather than delegated to external C or C++ GIS libraries. This approach provides consistent cross-platform behaviour, fewer system-level dependencies, tighter control over performance and correctness, and greater flexibility for continued research and development

Highlights

  • More than 700 tools for geomorphometry, hydrology, LiDAR, remote sensing, vector analysis, spatial statistics, and general geospatial processing
  • Publicly available interfaces for Python, R, and QGIS
  • A modular, high-performance backend written entirely in Rust, including wbprojection, wbraster, wbvector, wblidar, wbspatialstats, and wbtopology open-source (MIT/Apache licensed) backend libraries
  • Most of the 80+ tools that existed within the previous extension product have been migrated to the new Whitebox open core, including the tools for advanced surface curvature analysis and DEM processing
  • Native coordinate-reference-system handling and reprojection workflows via wbprojection
  • Expanded raster support (19 formats), including GeoTIFF, Cloud-Optimized GeoTIFF (COGs), GeoPackage Raster, and JPEG2000
  • Expanded vector support (12 formats), including Shapefile, GeoPackage, FlatGeobuf, GeoParquet, GeoJSON, TopoJSON, and GML
  • Modern point-cloud support, including LAS, LAZ, COPC, E57, and PLY
  • A dedicated topology engine supporting robust vector analysis, network analysis, route-event workflows, and linear referencing
  • Local-first processing on Windows, macOS, and Linux
  • A consistent analysis platform across Python scripts and R workflows
  • A newly updated and more advanced QGIS Processing toolbox

Whitebox has particular strengths in terrain analysis, geomorphometry, spatial hydrology, and LiDAR processing. Whitebox Next Gen carries these areas forward while substantially expanding the project’s capabilities for vector analysis, remote sensing, modern spatial data formats, and reproducible geospatial workflows.

Get started

The geomorphometry community is invited to explore Whitebox Next Gen, test its tools and workflows, report issues, and contribute to the continued development of the project.

DEMIX reveals which DEMs perform best

First round of the DEMIX exercise published final report

Global digital elevation models (DEMs) have become routine operational inputs across mapping, environmental monitoring, modelling and Earth-observation workflows. Their broad coverage and global availability have made them data that many practitioners simply take for granted. But knowledge of the terrain and topography of Earth’s surface is fundamental for monitoring and understanding terrestrial ecosystems and the planet’s habitability. For that reason, the CEOS Working Group on Calibration and Validation (WGCV) maintains a subgroup dedicated to digital topography and the quality of DEMs, the ‘Terrain Mapping Sub-Group (TMSG).

DEMs are representations of elevation in the form of a georectified grid, at global scale commonly derived from space-based interferometric synthetic aperture radar (InSAR) or stereoscopic optical observations, while regionally and locally airborne Laser induced detection and ranging (LIDAR) has become the primary source. DEMs can be classified as digital surface models (DSMs) when depicting the lower surface of the atmosphere or digital terrain models (DTMs) when depicting the top of the lithosphere (the Earth’s crust).

A wide variety of DEM products exist – each having different characteristics which suit different applications. It can be challenging to understand which of the available DEM products are fit for purpose for certain applications or regions. That is the practical gap addressed by the Digital Elevation Model Intercomparison Exercise (DEMIX) undertaken by TMSG. By exposing where widely used global DEMs agree, where they differ and how those differences affect rankings, the exercise delivered evidence to users for selecting a DEM rather than defaulting to one by habit.

DEM products were compared through a ‘wine contest,’ which identified requirements of openness, reproducibility, adaptability, and statistical rigour. The characteristics, evaluations, and ranking for each DEM were collected in a GIS database with over 50,000 entries. A quantitative assessment was carried out for each global DEM based on pixel-by-pixel differences of geomorphometric parameters against finer spatial resolution (1-5 m resolution) reference DEMs.

The study incorporates an unprecedented amount of high quality reference data covering a broad range of landforms and surface types, comparing their characteristics without resampling or interpolation. Users are invited to consult the results of DEMIX in order to make informed choices when needing to use a DEM, or in some cases a combination of the most relevant DEMs for their applications.

The final rankings of the wine contest, presented below, compared the six DEMs based on land cover (forest, urban, or barren), slope (cliff, steep, gentle, or flat), differences in elevation, slope and roughness, and statistical metrics.

Led by the European Commission’s Joint Research Centre (JRC), the CEOS Working Group on Calibration and Validation (WGCV) Terrain Mapping Subgroup (TMSG) established DEMIX in partnership with the International Society for Geomorphometry. Several members of the International Society of the Geomorphometry contributed significantly and we hosted various secssions at our latest conference. The exercise was performed over a three-year period starting in 2020, with community-wide calls for participation and an assembly of experts producing a number of peer-reviewed publications, culminating with a JRC reference report published in July 2026.

The objective of DEMIX was to propose a procedure to rank the available free and open global DEMs, taking into consideration user needs while promoting the implementation of FAIR (Findable, Accessible, Interoperable and Reusable) principles. The exercise compared DEMs based on criteria for land cover and terrain slope categories, representative testing, and a statistically sound ranking approach. The report outputs tailored recommendations regarding available DEM products that are not limited to one domain, geographic area, or landscape type, with flexibility for different user needs and applications.

A major challenge in the comparison of the global DEMs were the varying formats, data, and metadata contents. The adoption of a common set of standards would enhance the quality and interoperability of global DEMs and streamline the exchange and utilisation of DEM data for both providers and users. DEMIX makes the following recommendations for data providers:

● Adopt a standardised grid layout, with complete grid definitions and encodings following ISO/OGC rules.

● Include vertical datum information as part of raster DEM files.

● Clearly indicate the product version in the file names.

● Participate in future DEMIX rounds, which serve as a neutral platform for independently evaluating products before their release.

● Provide high quality, finer resolution, accurate, and multi-temporal elevation data to be used as a reference alongside DEM products.

Global DEMs are indispensable operational geospatial datasets, but they should not be treated as error-free or interchangeable. DEMIX delivered to the global community a systematic and transparent method to help choose elevation products more intelligently and help define what better elevation measurements should look like in the future.

Further Reading:

Read more in the DEMIX Final Report, or see other DEMIX publications below - mainly by our ISG members:

  • Benchmarking Elevation Plus Land Surface Parameters Finds FathomDEM and Copernicus DEM Win as Best Global DEMs (2025) and subsequent discussion (2026a, 2026b)
  • Ranking of 10 Global One-Arc-Second DEMs Reveals Limitations in Terrain Morphology Representation 2025
  • Novel Approach for Ranking DEMs: Copernicus DEM Improves One Arc Second Open Global Topography (2024)
  • Digital Elevation Models: Terminology and Definitions (2021)
  • The Digital Elevation Model Intercomparison eXperiment DEMIX, a community-based approach at global DEM benchmarking (2021)

Coffee Talk - GeoNadir

Building the world’s largest repository of FAIR drone mapping data

Paul Mead GeoNadir

September 2nd, 2026
11:30 (UTC)

Bio: Paul Mead is Co-founder and Head of Business & Strategy at GeoNadir, a platform building the world’s largest repository of FAIR drone mapping data to support environmental decision-making. A non-technical co-founder, Paul focuses on strategy, commercial development, and building high-performing teams — a discipline shaped by 11 years in the New Zealand Army, including leading a bomb disposal team in Lebanon and service in Afghanistan and the Asia-Pacific region. Since 2013 he has built and led several small businesses and start-ups, with drones as a central technology since 2016. Alongside GeoNadir, Paul co-founded She Maps, an internationally recognised drone and geospatial STEM education organisation, and has helped raise over $7M in equity and non-equity funding across his ventures. He holds a Master of Philosophy from the University of Queensland and is a Chartered Manager.

Abstract: GeoNadir set out to build the world’s largest repository of FAIR drone mapping data, and ended up building something more powerful. Drones have democratised data capture, but the same democratised ability to generate decision-grade insights from the data is lacking. I’ll walk through how we’re closing that gap: the data infrastructure behind hosting large volumes of UAV imagery and derived products, and how we’re moving from manual analysis toward automated, repeatable geospatial analytics.