suppressPackageStartupMessages({
library(terra)
library(ggplot2)
library(spatialsample)
library(rsample)
library(sf)
library(dplyr)
})
This file contains hidden or bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
#' Create a composite raster with a VRT pixel function | |
#' @param src_files a character vector with the paths to the source rasters. | |
#' @param outfile a character path to the output raster. | |
#' @param fun a character with the pixel function to apply. | |
#' @param t_srs a character with the target SRS. If empty string "", the spatial | |
#' reference of src_files[1] will be used. | |
#' @param vrt_options a character vector with additional options to pass to the | |
#' gdalbuildVRT command. | |
#' @param warp_options a character vector with additional options to pass to the | |
#' gdalwarp command. |
This file contains hidden or bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
#' Download geospatial embeddings for Major-TOM datsets | |
#' @param bbox numeric vector of length 4 representing the bounding box in | |
#' longitude and latitude coordinates (WGS84). Ordered as c(xmin, xmax, ymin, | |
#' ymax). | |
#' @param embed character vector of length 1 representing the embedding to | |
#' download. | |
#' @param limit integer vector of length 1 representing the maximum number of | |
#' rows to return. Useful for testing, if NULL then all rows within the bbox | |
#' are returned. |
This file contains hidden or bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
library(reticulate) | |
py_require(c("dask", "distributed", "rioxarray")) | |
dask.distributed <- reticulate::import("dask.distributed") | |
rioxarray <- reticulate::import("rioxarray") | |
Sys.setenv(GDAL_VRT_ENABLE_PYTHON = "YES") | |
This file contains hidden or bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
# how to get at Digital Eath Africa GeoMAD datasets | |
library(vrtility) | |
mirai::daemons(10) | |
bbox <- gdalraster::bbox_from_wkt( | |
wkt = "POINT (46.33 -15.9)", | |
extend_x = 0.15, | |
extend_y = 0.2 | |
) |
This file contains hidden or bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
library(vrtility) | |
# only use 4 daemons - this is limit for concurrent requests: https://documentation.dataspace.copernicus.eu/Quotas.html | |
mirai::daemons(4) | |
Sys.setenv( | |
AWS_VIRTUAL_HOSTING = "FALSE", | |
AWS_ACCESS_KEY_ID = Sys.getenv("CDSE_ACCESS_KEY"), | |
AWS_SECRET_ACCESS_KEY = Sys.getenv("CDSE_SECRET_KEY"), | |
AWS_S3_ENDPOINT = "eodata.dataspace.copernicus.eu" |
OlderNewer