Crests and throughs along tracks

Aim

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.

Code

Language: C++
Compiler: CodeBlocks - Windows
Version: 2009-07-20
Code: minmax.zip

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.

About

My name is Mauro Alberti. I'm a freelance 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

Contact

Mauro Alberti

alberti.m65@gmail.com

mobile: +39 338 8022247
skype: mauro.alberti.siena
Address: Strada per Casale 56,
Verolengo (TO) 10038, Italy

P.I. 01292740527