getData: Get geographic data

Description Usage Arguments Details Value References


Get global geographic data. Data are read from files that are first downloaded if not present in the indicated path. Function ccodes returns country names and the ISO codes


getData(name, download=TRUE, path='', ...)



Data set name, currently supported are 'GADM', 'countries', 'SRTM', 'alt', 'worldclim', and 'worldclim_past'. See Details for more info


Logical. If TRUE data will be downloaded if not locally available


Character. Path name indicating where to store the data. Default is the current working directory


Additional required (!) parameters. These are data set specific. See Details


'alt' stands for altitude (elevation); the data were aggregated from SRTM 90 m resolution data between -60 and 60 latitude. 'GADM' is a database of global administrative boundaries. 'worldclim' is a database of global interpolated climate data from 1960-1990 averages, and 'worldclim_past' contains inferred and interpolated climate data from mid-holocene (6,000 years ago) and last glacial maximum (22,000 years ago). 'SRTM' refers to the hole-filled CGIAR-SRTM (90 m resolution). 'countries' has polygons for all countries at a higher resolution than the 'wrld_simpl' data in the maptools pacakge .

If name is 'alt' or 'GADM' you must provide a 'country=' argument. Countries are specified by their 3 letter ISO codes. Use getData('ISO3') to see these codes. In the case of GADM you must also provide the level of administrative subdivision (0=country, 1=first level subdivision). In the case of alt you can set 'mask' to FALSE. If it is TRUE values for neighbouring countries are set to NA. For example:

getData('GADM', country='FRA', level=1)

getData('alt', country='FRA', mask=TRUE)

If name is 'SRTM' you must provide 'lon' and 'lat' arguments (longitude and latitude). These should be single numbers somewhere within the SRTM tile that you want.

getData('SRTM', lon=5, lat=45)

If name='worldclim' or name='worldclim_past' you must also provide arguments var, and a resolution res. Valid variables names are 'tmin', 'tmax', 'prec' and 'bio'. For 'worldclim_past' you also need to specify the past time, 'lgm' or 'mid', as past="lgm". Valid resolutions are 0.5 (not for past climates), 2.5, 5, and 10 (minutes of a degree). In the case of res=0.5, you must also provide a lon and lat argument for a tile; for the lower resolutions global data will be downloaded. In all cases there are 12 (monthly) files for each variable except for 'bio' which contains 19 files.

getData('worldclim', var='tmin', res=0.5, lon=5, lat=45)

getData('worldclim', var='bio', res=10)

getData('worldclim_past', var='bio', res=10, past="lgm")

To get (projected) future climate data (CMIP5), you must provide arguments var and res as above. Only resolutions 2.5, 5, and 10 are currently available. In addition, you need to provide model, rcp and year. For example,

getData('CMIP5', var='tmin', res=10, rcp=85, model='AC', year=70)

function (var, model, rcp, year, res, lon, lat, path, download = TRUE)

'model' should be one of "AC", "BC", "CC", "CE", "CN", "GF", "GD", "GS", "HD", "HG", "HE", "IN", "IP", "MI", "MR", "MC", "MP", "MG", or "NO".

'rcp' should be one of 26, 45, 60, or 85.

'year' should be 50 or 70

Not all combinations are available. See for details.


A spatial object (Raster* or Spatial*)


MoisesExpositoAlonso/rbioclim documentation built on May 7, 2019, 4:58 p.m.