Laura Toma and I wrote a chapter on algorithms for elevation models (“Terrain modeling for the geosciences”) for the computing handbook set published by CRC Press in 2014. If you would like to read it, please mail me.
My remaining work on algorithms for elevation models concerns algorithms to compute optimal triangulations of the sea surface; algorithms to analyse the flow of water across a terrain, algorithms to compute visibility maps: which part of the terrain is visible from a given viewpoint? The odd one out (of which the applications are not restricted to geographical elevation data) is our paper:
To be able to estimate the sea surface height in the past, we need to work with measurements from tide gauge stations located on the coast. We can triangulate the sea surface between these stations, and use the triangulation for interpolation. But what triangulation gives the best interpolations? For an answer, we look at recent times for which sea surface measurements from tide gauge stations and satellites are available, and we look for a triangulation that best matches those measurements:
In hydrological studies we analyse the flow of water on a terrain, for example to study the effects of possible human intervention or to estimate risks of erosion. Assuming that flowing water always takes the way down that descends most steeply, we can predict the flow of water on a terrain from a digital elevation model. Thus we may try to compute, for every point (or river) q of a terrain, how much water flows to q, and where it comes from (the watershed of q). In principle, watersheds have a hierarchical structure: if water flows from p to q, then the watershed of p is contained in the watershed of q. Simple as it may sound, it is far from trivial to actually compute flow paths, watersheds, and their hierarchical structure. Some of the challenges, which we try to deal with in our work, are:
When analysing the flow of water across a terrain, represented by a digital elevation model, two main approaches can be distinguished. One approach treats the terrain as a continuous surface, of which the exact shape is estimated by interpolating between sample points for which an elevation is given. The other approach treats the terrain as a network of cells, in which each cell is treated as an atomic unit, and water that arrives in any cell continues its course down-hill to one or more neighboring cells.