Yale Center for Global Change

Overview

  • Spatial Objects
  • Points
  • Rasters
  • Workflow

Spatial objects

sp = SpatialPoints(coords)
str(sp)
## Formal class 'SpatialPoints' [package "sp"] with 3 slots
##   ..@ coords     : num [1:100, 1:2] -1.163 0.3609 -0.4343 -0.6074 -0.0481 ...
##   .. ..- attr(*, "dimnames")=List of 2
##   .. .. ..$ : NULL
##   .. .. ..$ : chr [1:2] "x" "y"
##   ..@ bbox       : num [1:2, 1:2] -2.51 -2.42 2.29 2.87
##   .. ..- attr(*, "dimnames")=List of 2
##   .. .. ..$ : chr [1:2] "x" "y"
##   .. .. ..$ : chr [1:2] "min" "max"
##   ..@ proj4string:Formal class 'CRS' [package "sp"] with 1 slot
##   .. .. ..@ projargs: chr NA

Rasters

  • Elevation (SRTM 90m resolution raster)
  • World Climate (Tmin, Tmax, Precip, BioClim rasters)
  • Countries from CIA factsheet (vector!)
  • Global Administrative boundaries (vector!)

alt text

Workflow

  • Download a global Maximum Temperature dataset
  • Crop it to a country you download
  • Calculate the overall range for each variable
  • Calculate the focal median with an 11x11 window
  • Create a transect across the region and extract the temperature data.

Using this tutorial

alt text