hexPlot: Generate hexmaps from position data

Description Usage Arguments Value

Description

Generates various intermediary files (detailed in other functions) for use in plotting If proportionByAnimal writes weights for each animal to 'proportionWeights.csv' Generates postscript hexmap plot in outFile

Usage

1
2
3
4
5
6
7
8
hexPlot(daybrs, outFile, hexPerDegree = 1, limits = NULL, extraCmd = NULL,
  landMaskCmd = NULL, hardLimit = NULL, hexMax = NULL,
  logCounts = FALSE, debug = FALSE, maxInterp = 7, gmtDir = "",
  contourFile = NULL, contourDepth = -200, dayScale = 1, interp = 1,
  addNumberToName = TRUE, proportion = FALSE, showMax = TRUE,
  proportionByAnimal = FALSE, seasonLimit = 0, uniqueAnimals = FALSE,
  outlineCount = NULL, propCount = NULL, addPlus = TRUE,
  animalWeights = NULL, ...)

Arguments

daybrs

A data frame containing 'animal','rdate','lat','lon' (e.g. the output from readDouglas)

outFile

File location to write postscript map generated by GMT to

hexPerDegree

How many hexs should fit horizontally in 1 degree of longitude at the equator

limits

A vector of (southernMostLatitude,northernMostLatitude,westernMostLongitude,easternMostLongitude) for final hexMap or NULL to select automatically

extraCmd

If not NULL, vector of system commands to run after other GMT commands and before scale (e.g. extra labels with pstext)

landMaskCmd

If not NULL, vector of system commands to run after hex output and before pscoast (e.g. polygons over hexes)

hardLimit

A vector of (southernMostLatitude,northernMostLatitude,westernMostLongitude,easternMostLongitude) for cropping hexs or NULL to not crop

hexMax

Sets hexs with counts > hexMax to hexMax (for creating maps with the highest hex counts being something like "100+")

logCounts

If true, base color scale on logged hex counts

debug

If true, browser() before returning from makeHexs()

maxInterp

Maximum number of missing days to fill between two data points (missing days<=maxInterp)

gmtDir

System directory for calling GMT functions (e.g. if pscoast is in /home/bin/ and /home/bin/ isn't in path then gmtDir="/home/bin/")

contourFile

If not NULL, run pscontour at contourDepth meters based on this GMT elevation file

contourDepth

Depth to draw contour in meters

dayScale

Divide day totals by this number (useful for intervals less than 1 day)

interp

number of positions per day to interpolate

addNumberToName

Add number of days and turtles to file name?

proportion

Divide counts by sum(counts)

showMax

show maximum on scale ticks

proportionByAnimal

weight points to add up to one for each animal

seasonLimit

When using proportionByAnimal, if less than seasonLimit points then total animal weight = # points/seasonLimit otherwise 1

uniqueAnimals

Count unique occurrences of animal in a hex

outlineCount

Draw an outline around all hexes with count >= outlineCount

propCount

Draw an outline around enough hexs to cover propCount proportion of the data (overrides outlineCount)

addPlus

Add a "+" to maximum label in scale (if using hexMax)

animalWeights

Vector (with names of animals) of custom weighting for animals if proportionByAnimal is TRUE (e.g. monthly plots of proportion)

...

Any other arguments to be passed to runGMT

Value

Vector of (lowXLim,highXlim,lowYLim,highYLim,maxHexCount) can be stored and passed as limits=return, hardLimit=return to further hexPlot calls to use the same range of latitude and longitude


sherrillmix/GMTHexMapR documentation built on May 29, 2019, 9:23 p.m.