mapGriddedData: Produce maps of global gridded data at half degree resolution
In rworldmap: Mapping Global Data

mapGriddedData

R Documentation

Produce maps of global gridded data at half degree resolution

Description

Produce maps of global gridded data at half degree resolution

Usage

mapGriddedData(
  dataset = "",
  nameColumnToPlot = "",
  numCats = 5,
  catMethod = "quantiles",
  colourPalette = "heat",
  xlim = c(-180, 180),
  ylim = c(-80, 90),
  mapRegion = "world",
  addLegend = TRUE,
  addBorders = "coarse",
  borderCol = "grey",
  oceanCol = NA,
  landCol = NA,
  plotData = TRUE,
  aspect = 1,
  lwd = 1
)

Arguments

`dataset`	gridded data either as a : SpatialGridDataFrame (R object defined in package sp) file name of a GridAscii file - this is an Esri format 2D R matrix or array (rows by columns)
`nameColumnToPlot`	name of column containing the data to plot
`numCats`	number of categories to put the data in, may be overidden if catMethod ='pretty'
`catMethod`	method for categorisation of data "pretty", "fixedWidth","diverging", "logFixedWidth","quantiles","categorical", or a numeric vector defining breaks
`colourPalette`	a string describing the colour palette to use, choice of : "palette" for the current palette a vector of valid colours, e.g. =c('red','white','blue') or output from RColourBrewer one of "heat", "diverging", "white2Black", "black2White", "topo", "rainbow", "terrain", "negpos8", "negpos9"
`xlim`	map extents c(west,east), can be overidden by mapRegion
`ylim`	map extents c(south,north), can be overidden by mapRegion
`mapRegion`	a country name from getMap()[['NAME']] or 'world','africa','oceania','eurasia','uk' sets map extents, overrides xlim,ylim
`addLegend`	whether to add a legend or not
`addBorders`	options for country borders, 'low','coarse' = low or coarse resolution, 'coasts' = coasts only, 'none' or NA for none
`borderCol`	the colour for country borders
`oceanCol`	a colour for the ocean if the grid values are NA
`landCol`	a colour to fill countries if the grid values are NA over land
`plotData`	whether to plotData, if FALSE a legend can be added on its own
`aspect`	aspect for the map, defaults to 1, if set to 'variable' uses same method as plot.Spatial in sp
`lwd`	line width for country borders

Details

Plots a map of global half degree gridded data, allowing classification, colours and regions to be set.

Certain catMethod and colourPalette options go well together. e.g. "diverging" and "diverging", "categorical" and "rainbow"

Value

invisibly returns a list containing the data and main options used for the map, the list can be passed to addMapLegend along with additional options to allow greater flexibility in legend creation.

Author(s)

andy south #@importFrom maptools readAsciiGrid

Examples


## mapping continuous data
data(gridExData,envir=environment(),package="rworldmap")
gridExData <- get("gridExData")
mapGriddedData(gridExData)

## reclassing continuous data to categorical & mapping
data(gridExData,envir=environment(),package="rworldmap")
#find quartile breaks
cutVector <- quantile(gridExData@data[,1],na.rm=TRUE)
#classify the data to a factor
gridExData@data$categories <- cut( gridExData@data[,1]
                                      , cutVector, include.lowest=TRUE)
#rename the categories
levels(gridExData@data$categories) <- c('low', 'med', 'high', 'vhigh')
#mapping
mapGriddedData( gridExData, nameColumnToPlot= 'categories'
              , catMethod='categorical')

rworldmap documentation built on Oct. 16, 2023, 5:07 p.m.