rotsmoother - Get mean rotations and covariance matrices from set of finite rotations
gmt rotsmoother [ rottable ] [ -A ] [ -C ] [ -N ] [ -S ] [ -Tages ] [ -V[level] ] [ -W ] [ -Z ] [ -bbinary ] [ -dnodata ] [ -eregexp ] [ -hheaders ] [ -iflags ] [ -oflags ] [ -sflags ] [ -:[i|o] ] [ --PAR=value ]
Note: No space is allowed between the option flag and the associated arguments.
rotsmoother reads a table of total reconstructions and computes mean rotations (and optionally covariance matrices) for sub-groups of rotations based on rotation age.
- Name of a rotation table containing (lon, lat, time, angle, [weight]) values.
- Use opening angles as a proxy for time. Suitable when no time can be assigned to the rotations. In this case, input is expected to contain lon lat angle [weight] records and -T settings refer to angles instead of time. [Default expects lon lat time angle [weight] and -T refers to time].
- Compute covariance matrix for each mean rotation. This is done by converting each finite rotation to a quaternion, determining the mean quaternion (rotation) and the consider all rotations as perturbation to the mean rotation. From these perturbations we determine the covariance matrix.
- Ensure all poles are in northern hemisphere [Default ensures positive opening angles].
- Ensure all poles are in southern hemisphere [Default ensures positive opening angles].
- Sets the desired groups of ages. For a single time append the desired time. For an equidistant range of reconstruction ages give -Tstart/stop/inc or -Tstart/stop/npoints+n. For an non-equidistant set of reconstruction ages please pass them via the first column in a file, e.g., -Tagefile. The ages indicate read or generated becomes the bin-boundaries and we output the average time of all rotations inside each bin.
- -V[level] (more …)
- Select verbosity level [c].
- Expect weights in last column for a weighted mean rotation [no weights].
- Report negative opening angles [positive].
- -bi[ncols][t] (more …)
- Select native binary format for primary input. [Default is 2 input columns].
- -bo[ncols][type] (more …)
- Select native binary output. [Default is same as input].
- -d[i|o]nodata (more …)
- Replace input columns that equal nodata with NaN and do the reverse on output.
- -e[~]”pattern” | -e[~]/regexp/[i] (more …)
- Only accept data records that match the given pattern.
- -h[i|o][n][+c][+d][+rremark][+rtitle] (more …)
- Skip or produce header record(s).
- -icols[+l][+sscale][+ooffset][,…][,t[word]] (more …)
- Select input columns and transformations (0 is first column, t is trailing text, append word to read one word only).
- -ocols[,…][t[word]] (more …)
- Select output columns (0 is first column; t is trailing text, append word to write one word only).
- -s[cols][+a|+r] (more …)
- Set handling of NaN records.
- -:[i|o] (more …)
- Swap 1st and 2nd column on input and/or output.
- -n[b|c|l|n][+a][+bBC][+c][+tthreshold] (more …)
- Select interpolation mode for grids.
- -^ or just -
- Print a short message about the syntax of the command, then exits (NOTE: on Windows just use -).
- -+ or just +
- Print an extensive usage (help) message, including the explanation of any module-specific option (but not the GMT common options), then exits.
- -? or no arguments
- Print a complete usage (help) message, including the explanation of all options, then exits.
- Temporarily override a GMT default setting; repeatable. See gmt.conf for parameters.
Geodetic versus Geocentric Coordinates¶
All spherical rotations are applied to geocentric coordinates. This means that incoming data points and grids are considered to represent geodetic coordinates and must first be converted to geocentric coordinates. Rotations are then applied, and the final reconstructed points are converted back to geodetic coordinates. This default behavior can be bypassed if the ellipsoid setting PROJ_ELLIPSOID is changed to Sphere.
To smooth rotation groups in increments of 3 Myr and ensure northern hemisphere poles, try
gmt rotsmoother rotations.txt -N -T3/3/30 -V > rot_means.txt