knitr::opts_chunk$set( collapse = TRUE, fig.width=7, fig.height=5, comment = "#>" )
We want to rasterize shapefile data from one bird species using the rasterizeRange() function with different argument options:
library(rasterSp); library(terra) r_pass_dom <- rasterizeRange(dsn=system.file("extdata", "passer_domesticus.shp", package = "rasterSp"), id="SCINAME", resolution=0.5, save=F, df=F, touches=T, seasonal=NA, origin=NA, presence=c(1,2)) plot(r_pass_dom) global(r_pass_dom, fun="sum", na.rm=T) r_pass_dom <- rasterizeRange(dsn=system.file("extdata", "passer_domesticus.shp", package = "rasterSp"), id="SCINAME", resolution=0.5, save=F, df=F, touches=T, seasonal=c(1,2), origin=1, presence=c(1,2)) plot(r_pass_dom) global(r_pass_dom, fun="sum", na.rm=T) r_pass_dom <- rasterizeRange(dsn=system.file("extdata", "passer_domesticus.shp", package = "rasterSp"), id="SCINAME", resolution=0.5, save=F, df=F, touches=F, seasonal=c(1,2), origin=1, presence=c(1,2)) global(r_pass_dom, fun="sum", na.rm=T) r_pass_dom <- rasterizeRange(dsn=system.file("extdata", "passer_domesticus.shp", package = "rasterSp"), id="SCINAME", resolution=0.5, save=F, df=T, touches=T, seasonal=c(1,2), origin=1, presence=c(1,2)) head(r_pass_dom) r_pass_dom <- rasterizeRange(dsn=system.file("extdata", "passer_domesticus.shp", package = "rasterSp"), id="SCINAME", resolution=0.5, save=F, df=F, touches=T, getCover=T, seasonal=c(1,2), origin=1, presence=c(1,2)) plot(r_pass_dom) r_pass_dom <- rasterizeRange(dsn=system.file("extdata", "passer_domesticus.shp", package = "rasterSp"), id="SCINAME", resolution=0.5, save=F, df=F, touches=F, getCover=T, seasonal=c(1,2), origin=1, presence=c(1,2)) plot(r_pass_dom)
Next we test if the function also works with multiple species:
library(dplyr) r_birds <- rasterizeRange(dsn=list.files(paste0(system.file(package = "rasterSp"), "/extdata"), pattern=".shp", full.names=T), id="SCINAME", resolution=0.5, save=F, df=F, touches=F, getCover=F, seasonal=c(1,2), origin=1, presence=c(1,2)) plot(r_birds) global(r_birds, fun="sum", na.rm=T) r_birds <- rasterizeRange(dsn=list.files(paste0(system.file(package = "rasterSp"), "/extdata"), pattern=".shp", full.names=T), id="SCINAME", resolution=0.5, save=F, df=F, touches=T, getCover=F, seasonal=c(1,2), origin=1, presence=c(1,2)) global(r_birds, fun="sum", na.rm=T) r_birds <- rasterizeRange(dsn=list.files(paste0(system.file(package = "rasterSp"), "/extdata"), pattern=".shp", full.names=T), id="SCINAME", resolution=0.5, save=F, df=F, touches=T, getCover=T, seasonal=c(1,2), origin=1, presence=c(1,2)) plot(r_birds) r_birds <- rasterizeRange(dsn=list.files(paste0(system.file(package = "rasterSp"), "/extdata"), pattern=".shp", full.names=T), id="SCINAME", resolution=0.5, save=F, df=T, touches=T, seasonal=c(1,2), origin=1, presence=c(1,2)) r_birds %>% group_by(species) %>% summarise(n=n())
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.