Get geographic data
Get geographic data for anywhere in the world. Data are read from files that are first downloaded if necessary. Function
ccodes returns country names and the ISO codes
getData(name, download=TRUE, path='', ...) ccodes()
Data set name, currently supported are 'GADM', 'countries', 'SRTM', 'alt', and 'worldclim'. See Details for more info
TRUEdata 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. '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 package .
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)
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)
name='worldclim' you must also provide arguments
var, and a resolution
res. Valid variables names are 'tmin', 'tmax', 'prec' and 'bio'. Valid resolutions are 0.5, 2.5, 5, and 10 (minutes of a degree). In the case of
res=0.5, you must also provide a
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)
To get (projected) future climate data (CMIP5), you must provide arguments
res as above. Only resolutions 2.5, 5, and 10 are currently available. In addition, you need to provide
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 www.worldclim.org for details.
A spatial object (Raster* or Spatial*)