knitr::opts_chunk$set( collapse = TRUE, comment = "#>" )
First, we look at the BirdLife Range Maps of some invasive bird species (European starling, House sparrow and Sacred ibis).
library(rasterSp) sturnus_vulgaris <- sf::st_read(system.file("extdata", "sturnus_vulgaris.shp", package = "rasterSp")) passer_domesticus <- sf::st_read(system.file("extdata", "passer_domesticus.shp", package = "rasterSp")) threskiornis_aethiopicus <- sf::st_read(system.file("extdata", "threskiornis_aethiopicus.shp", package = "rasterSp"))
We want to compare the native origin (origin = 1, blue) with the introduced origin (origin = 3, red).
library(sf); library(dplyr) sturnus_vulgaris %>% st_geometry() %>% plot(main="European starling") sturnus_vulgaris %>% filter(ORIGIN==1) %>% st_geometry() %>% plot(add=T, col="blue") sturnus_vulgaris %>% filter(ORIGIN==3) %>% st_geometry() %>% plot(add=T, col="blue") passer_domesticus %>% st_geometry() %>% plot(main="House sparrow") passer_domesticus %>% filter(ORIGIN==1) %>% st_geometry() %>% plot(add=T, col="blue") passer_domesticus %>% filter(ORIGIN==3) %>% st_geometry() %>% plot(add=T, col="red")
threskiornis_aethiopicus %>% filter(ORIGIN==1) %>% st_geometry() %>% plot(col="blue", main="Sacred ibis (native)") threskiornis_aethiopicus %>% filter(ORIGIN==3) %>% st_geometry() %>% plot(col="red", main="Sacred ibis (introduced)")
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.