inst/doc/hwsdr-v2-vignette.R

In hwsdr: Interface to the 'HWSD' Web Services

## ----setup, include = FALSE---------------------------------------------------
knitr::opts_chunk$set(
  collapse = TRUE,
  comment = "#>"
)

# load the library
library(ncdf4)
library(terra)
library(sf)
library(hwsdr)

sand <- terra::rast(system.file(
  package = "hwsdr",
  "extdata/d1_sand.tif"))

# extract values
values <- readRDS(system.file(
  package = "hwsdr",
  "extdata/d1_values.rds")
)


## ----eval = FALSE-------------------------------------------------------------
#  # set the ws_path variable using a FULL path name
#  path <- ws_download(
#    ws_path = "/your/full/path",
#    verbose = TRUE
#  )

## ----eval = FALSE-------------------------------------------------------------
#  usethis::edit_r_environ()

## ----eval = FALSE-------------------------------------------------------------
#  values <- ws_subset(
#    site = "HWSD_V2",
#    location = c(-81, 34),
#    param = c("SAND","SILT"),
#    layer = "D1",
#    version = "2.0",
#    ws_path = "/your/full/path"
#  )

## -----------------------------------------------------------------------------
print(values)

## ----eval = FALSE-------------------------------------------------------------
#  sand <- ws_subset(
#      location = c(-81, 32, -80, 34),
#      param = "SAND",
#      layer = "D1",
#      version = "2.0",
#      ws_path = Sys.getenv("WS_PATH"),
#      # ws_path = "/your/full/path",
#      internal = TRUE
#    )

## -----------------------------------------------------------------------------
terra::plot(sand)

## ----eval = FALSE-------------------------------------------------------------
#    a <- sf::st_sf(a = 1:2,
#                  geom = sf::st_sfc(
#                    sf::st_point(c(34, -81)),
#                    sf::st_point(c(32, -80))),
#                  crs = 4326)
#  
#    t_sand <- a %>%
#      sf::st_bbox() %>%
#      ws_subset(
#        version = "2.0",
#        param = "SAND",
#        layer = "D1",
#        ws_path = Sys.getenv("WS_PATH")
#        )

## -----------------------------------------------------------------------------
terra::plot(sand)