## ---- echo=FALSE, include=FALSE----------------------------------------------- library(knitr) opts_chunk$set(fig.width = 5, fig.height = 5, fig.cap = '', collapse = TRUE ) library(raster) library(rgdal) ## ---- echo=FALSE-------------------------------------------------------------- knitr::include_graphics("img/sphere.png") ## ----------------------------------------------------------------------------- library(rgdal) epsg <- make_EPSG() i <- grep("France", epsg$note, ignore.case=TRUE) # first three epsg[i[1:3], ] ## ---- crs1-------------------------------------------------------------------- library(raster) library(rgdal) f <- system.file("external/lux.shp", package="raster") p <- shapefile(f) p ## ----------------------------------------------------------------------------- crs(p) ## ----------------------------------------------------------------------------- pp <- p crs(pp) <- NA crs(pp) crs(pp) <- CRS("+proj=longlat +datum=WGS84") crs(pp) ## ----------------------------------------------------------------------------- newcrs <- CRS("+proj=robin +datum=WGS84") ## ----------------------------------------------------------------------------- rob <- spTransform(p, newcrs) rob ## ----------------------------------------------------------------------------- p2 <- spTransform(rob, CRS("+proj=longlat +datum=WGS84")) ## ---- prj1-------------------------------------------------------------------- r <- raster(xmn=-110, xmx=-90, ymn=40, ymx=60, ncols=40, nrows=40) r <- setValues(r, 1:ncell(r)) r plot(r) ## ----------------------------------------------------------------------------- newproj <- "+proj=lcc +lat_1=48 +lat_2=33 +lon_0=-100 +ellps=WGS84" ## ----------------------------------------------------------------------------- pr1 <- projectRaster(r, crs=newproj) crs(pr1) ## ----------------------------------------------------------------------------- pr2 <- projectRaster(r, crs=newproj, res=20000) pr2 ## ----------------------------------------------------------------------------- pr3 <- projectExtent(r, newproj) # Set the cell size res(pr3) <- 200000 ## ---- prj3-------------------------------------------------------------------- pr3 <- projectRaster(r, pr3) pr3 plot(pr3)