The original C implementation of DTW for Stratigraphy
4. The original C implementation of DTW for Stratigraphy#
The application of DTW to the problem of stratigraphic depth and depth matching as implemented by Pälike [G2] (software) operates with the following schema:
Fig. 4.1 Generalised and extended constraints for DTW stratigraphic matching. The allowed region for the search path is shaded in grey. Each segment, which is defined by bio- or magnetostratigraphic constraints at its end points and denoted by the subscript \(s\), is computed separately by the algorithm. The local constraint is that of Itakura [2] , resulting in a parallelogram with slopes 0.5 and 2 for each segment. Uncertainty intervals for the tie points can be specified by using appropriate values for \(\gamma_s\) and \(\delta_s\). The last path point of segment 1 then forms the first path point for segment 2, and so on. (from [G2])#
This plain C
software tool is configured via text parameter files and has the following features:
Tie-points can be specified (with potential
"slack"
), which allows for realistic constraints.Multiple simultaneous fits (multiple proxies) can be applied at the same time, with potential weigting.
(Multiple) band-pass filterings can be applied to series during fitting (mostly applicable for cyclostratigraphy against astronomical age models.)
The software and an example dataset are available here: https://paloz.marum.de/git/mydtw/dtw-spm