PSWIGGLE
Section: Misc. Reference Manual Pages (l)
Updated: DATE
Index
Return to Main Contents
#ifdef SI
#define UNIT cm
#define LEN 15
#else
#define UNIT inch
#define LEN 5
#endif
NAME
pswiggle - Plot anomaly along track on a map
SYNOPSIS
pswiggle xyz_files -Jparameters -Rwest/east/south/north[r]
-Zscale [ -Adirection ] [ -Btickinfo ] [ -Ccenter ] [ -Dgap ] [ -Eazimuth/elevation ]
[ -Fred/green/blue ] [ -Gfillrgb ] [ -H ]
[ -Ifix_az ] [ -K ] [ -M[flag] ] [ -N ] [ -O ] [ -P ]
[ -S[x]lon0/lat0/length[units] ] [ -Ttrack_pen ]
[ -U[/dx/dy/][label] ] [ -V ] [ -Wwiggle_pen ] [ -Xx-shift ]
[ -Yy-shift ] [ -ccopies ] [ -: ]
DESCRIPTION
pswiggle reads (x,y,z) triplets from files [or standard input] and plots z as a function of
distance along track. This means that two consecutive (x,y) points define the local distance axis,
and the local z axis is then perpendicular to the distance axis. The user may set a preferred
positive anomaly plot direction, and if the positive normal is outside
the plus/minus 90 degree window around the preferred direction, then 180 degrees are added to the
direction. Either the positive or the negative wiggle may be shaded. The resulting PostScript
code is written to standard output.
- files
-
List one or more file-names. If no files are given, psxy will read standard input.
#include "explain_-j.txt"
#include "explain_-R.txt"
- -Z
-
Gives anomaly scale in data-units/UNIT.
OPTIONS
No space between the option flag and the associated arguments.
- -A
-
Sets the preferred positive direction. Positive wiggles will "gravitate" towards that
direction.
#include "explain_-b.txt"
- -C
-
Subtract center from the data set before plotting [0].
- -D
-
Means there is a data gap if 2 consecutive points are more than gap distance
units apart. For longitude/latitude data gap is in km, else it is in the
user's units.
- -E
-
Sets the viewpoint's azimuth and elevation [180/90]
- -F
-
Sets the color used for Frame and annotation. [Default is 0/0/0 (black)]
- -G
-
Gray shade or color for filling of positive anomalies. Grayshade range from 0 (black) to 255 (white).
Color range from 0-255 for each of the rgb components. [Default is black]
#include "explain_-H.txt"
- -I
-
Set a fixed azimuth projection for wiggles [Default uses track azimuth, but see -A].
#include "explain_-K.txt"
- -M
-
Multiple segment file. Segments are separated by a record whose first character is <flag>.
[Default is '>'].
- -N
-
Paint negative wiggles instead of positive [Default].
#include "explain_-O.txt"
#include "explain_-P.txt"
- -S
-
Draws a simple vertical scale centered on lon0/lat0. Use -Sx to specify cartesian coordinates instead.
length is in z units, append unit name for labeling
- -T
-
Draw track [Default is no track]. Append pen attributes to use [Defaults: width = 1, color = 0/0/0, texture = solid].
#include "explain_-U.txt"
#include "explain_-V.txt"
- -W
-
Draw wiggle outline [Default is no outline]. Append pen attributes to use
[Defaults: width = 1, color = 0/0/0, texture = solid].
#include "explain_-XY.txt"
#include "explain_-t.txt"
#include "explain_-c.txt"
EXAMPLES
To plot the magnetic anomaly stored in the file track.xym along track @ 1000 nTesla/UNIT (after removing a mean value of 32000 Tesla),
using a LEN-UNIT-wide Polar Stereographic map ticked
every 5 degrees in Portrait mode, with positive anomalies in red on a blue track of width 1, try
pswiggle track.xym -R-20/10/-80/-60 -JP0/90/LEN -Z1000 -B5 -P -G255/0/0 -T1/0/0/255 -S1000 -V > track_xym.ps
BUGS
Sometimes the (x,y) coordinates are not printed with enough significant
digits, so the local perpendicular to the track swings around a lot. To
see if this is the problem, you should do this:
awk '{ if (NR > 1) print atan2(y-$1, x-$2); y=$1; x=$2; }' yourdata.xyz | more
(note that output is in radians; on some machines you need "nawk" to do this).
Then if these numbers jump around a lot, you may do this:
gmtset D_FORMAT %.12lg
awk '{ print NR, $0 }' yourdata.xyz | filter1d -Fb5 -N4/0 > smoothed.xyz
and plot this data set instead.
SEE ALSO
gmt, filter1d, psbasemap, splitxyz
#include "refs.i"
Index
- NAME
-
- SYNOPSIS
-
- DESCRIPTION
-
- OPTIONS
-
- EXAMPLES
-
- BUGS
-
- SEE ALSO
-
This document was created by
man2html,
using the manual pages.
Time: 07:11:16 GMT, January 07, 2025