knitr::opts_chunk$set(echo = TRUE)
library(sf) library(terra) carteQc <- mapselector::regions_simplified_Ouranos carteQc_r <- vect(carteQc) plot(carteQc_r) night_light <- rast("./data/VNL_v2_npp_2021_global_vcmslcfg_c202203152300.average_masked.tif") night_light_Qc <- crop(night_light, carteQc_r) values(night_light_Qc)[values(night_light_Qc) == 0] <- NA plot(night_light_Qc, add = TRUE)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.