psxy: Add Lines/Symbols to GMT Map

In gmt: Interface Between GMT Map-Making Software and R

Add Lines/Symbols to GMT Map

Description

Call GMT to add lines/symbols to a map and save in postscript format.

Usage

psxy(x, cmd="-J -R -Scp -W2p -O -K", file=getOption("gmt.file"))

Arguments

x	data frame, matrix, or filename containing the data to be plotted.
cmd	string of arguments passed to GMT shell command psxy.
file	filename where the map is saved.

Details

The data are arranged in two (Lon, Lat) or more columns, depending on the -S argument.

If x is a filename, the data should be tabular with or without a header, separated by commas or whitespace. The first line is interpreted as header if the first non-whitespace character is not minus, point, or number.

See the GMT documentation for details on psxy and other GMT commands.

Value

NULL, but the map is annotated and saved in postscript format.

The temporary GMT input file ‘xy.gmt’ is saved in directory dirname(tempdir()), for the user to view or edit. It is later removed by psclose().

See Also

Similar to points, lines, and postscript in native R graphics.

gmt, pscoast, psxy, pstext, psbar, and psclose work together to draw maps.

gmt-package gives an overview of the package.

Examples

## Not run: 
# Draw map and save as "map.ps" in current working directory
gmt(demo.par)
pscoast(demo.coast)
psxy(demo.xy)
pstext(demo.text, "-J -R -F+f+a+j -O -K")
psbar(demo.bar, ref=66)
psclose()
# See directory gmt/example for details

# Simple map and circles
pscoast("-JM12c -R158/192/-42/-8 -Di -Gdarkgreen -B10f5 -A2000 -K", "quake.ps")
psxy(data.frame(lon=quakes$lon, lat=quakes$lat, mag=10^quakes$mag/1e6),
     "-J -R -Scp -W0.5p,red -O", "quake.ps")

## End(Not run)

gmt documentation built on June 7, 2022, 1:11 a.m.