psxy function

Add Lines/Symbols to GMT Map

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

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.

Returns

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 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)