Crests and throughs along tracks


The algorithm determines local minimum and maximum heights of cyclical features such as megadunes. It uses a moving filter along the track. This filter selects the elevations representing a local minimum/maximum or a flat within a four-elevation range.


Language: C++
Compiler: CodeBlocks - Windows
Version: 2009-07-20

Program input

The input data format is:
track_id x y z
The program asks for a trend and an angular tolerance value used to analyze only the tracks with orientations within trend plus/minus angular tolerance. These parameters allow to emphasize structures sub-perpendicular to the chosen trend.

Program output

The output is an ASCII file with resulting local minima/maxima. Its format is:
   track_id, x, y, z, minmax_type
The last field represents the type of found point, i.e., it represents a minimum or a maximum.

Example of program output from multiple parallel tracks of ICESat data. Red dots represent the crests and blue dots represent the troughs of the megadunes. Base image: MOA satellite mosaic.


My name is Mauro Alberti. I'm a Gis and scientific programmer. My primary specializations are in scientific and GIS programming (OS tools, mainly Python, R, Fortran and C++), Earth Sciences (structural geology, Polar Sciences), geostatistical analysis. About me


Mauro Alberti

mobile: +39 338 8022247
skype: mauro.alberti.siena
Address: Via Canonica 32,
Milano (MI) 20100, Italy

P.I. 01292740527