inst/doc/spatial-data.R

In robotoolbox: Client for the 'KoboToolbox' API

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

## ----setup, echo = FALSE------------------------------------------------------
library(robotoolbox)
library(dplyr)
library(sf)
library(mapview)
mapview::mapviewOptions(
    fgb = FALSE,
    basemaps = c(
      "Esri.WorldImagery",
      "Esri.WorldShadedRelief",
      "OpenTopoMap",
      "OpenStreetMap"),
    layers.control.pos = "topright")

## ----asset_list, echo = FALSE-------------------------------------------------
l <- asset_list

## ----eval = FALSE-------------------------------------------------------------
# library(robotoolbox)
# library(dplyr)
# 
# # Retrieve a list of all assets (projects) from your KoboToolbox server
# asset_list <- kobo_asset_list()
# 
# # Filter the asset list to find the specific project and get its unique identifier (uid)
# uid <- filter(asset_list, name == "Spatial data") |>
#   pull(uid)
# 
# # Load the specific asset (project) using its uid
# asset <- kobo_asset(uid)
# asset

## ----echo = FALSE-------------------------------------------------------------
asset <- asset_spatial
asset

## ----eval = FALSE-------------------------------------------------------------
# df <- kobo_data(asset)
# glimpse(df)

## ----echo = FALSE-------------------------------------------------------------
df <- data_spatial
glimpse(df)

## -----------------------------------------------------------------------------
pull(df, point)
pull(df, point_wkt)

## -----------------------------------------------------------------------------
df |>
  select(starts_with("point_"))

## -----------------------------------------------------------------------------
pull(df, line)
pull(df, line_wkt)

## -----------------------------------------------------------------------------
pull(df, polygon)
pull(df, polygon_wkt)

## -----------------------------------------------------------------------------
point_sf <- st_as_sf(data_spatial,
                     wkt = "point_wkt", crs = 4326)
mapview(point_sf)

## -----------------------------------------------------------------------------
line_sf <- st_as_sf(data_spatial,
                     wkt = "line_wkt", crs = 4326)
mapview(line_sf)

## -----------------------------------------------------------------------------
poly_sf <- st_as_sf(data_spatial,
                    wkt = "polygon_wkt", crs = 4326)
mapview(poly_sf)